New hopes for the breast cancer treatment: perspectives on the oncolytic virus therapy.

Oncolytic virus (OV) therapy has emerged as a promising frontier in cancer treatment, especially for solid tumours. While immunotherapies like immune checkpoint inhibitors and CAR-T cells have demonstrated impressive results, their limitations in inducing complete tumour regression have spurred researchers to explore new approaches targeting tumours resistant to current immunotherapies. OVs, both natural and genetically engineered, selectively replicate within cancer cells, inducing their lysis while sparing normal tissues. Recent advancements in clinical research and genetic engineering have enabled the development of targeted viruses that modify the tumour microenvironment, triggering anti-tumour immune responses and exhibiting synergistic effects with other cancer therapies. Several OVs have been studied for breast cancer treatment, including adenovirus, protoparvovirus, vaccinia virus, reovirus, and herpes simplex virus type I (HSV-1). These viruses have been modified or engineered to enhance their tumour-selective replication, reduce toxicity, and improve oncolytic properties.Newer generations of OVs, such as Oncoviron and Delta-24-RGD adenovirus, exhibit heightened replication selectivity and enhanced anticancer effects, particularly in breast cancer models. Clinical trials have explored the efficacy and safety of various OVs in treating different cancers, including melanoma, nasopharyngeal carcinoma, head and neck cancer, and gynecologic malignancies. Notably, Talimogene laherparepvec (T-VEC) and Oncorine have. been approved for advanced melanoma and nasopharyngeal carcinoma, respectively. However, adverse effects have been reported in some cases, including flu-like symptoms and rare instances of severe complications such as fistula formation. Although no OV has been approved specifically for breast cancer treatment, ongoing preclinical clinical trials focus on four groups of viruses. While mild adverse effects like low-grade fever and nausea have been observed, the effectiveness of OV monotherapy in breast cancer remains insufficient. Combination strategies integrating OVs with chemotherapy, radiotherapy, or immunotherapy, show promise in improving therapeutic outcomes. Oncolytic virus therapy holds substantial potential in breast cancer treatment, demonstrating safety in trials. Multi-approach strategies combining OVs with conventional therapies exhibit more promising therapeutic effects than monotherapy, signalling a hopeful future for OV-based breast cancer treatments.

Predicting Mortality for COVID-19 Patients Admitted to an Emergency Department Using Early Warning Scores in Poland.

COVID-19 disease is characterised by a wide range of symptoms that in most cases resemble flu or cold. Early detection of infections, monitoring of patients' conditions, and identification of patients with worsening symptoms became crucial during the peak of pandemic. The aim of this study was to assess and compare the performance of common early warning scores at the time of admission to an emergency department in predicting in-hospital mortality in patients with COVID-19. The study was based on a retrospective analysis of patients with SARS-CoV-2 infection admitted to an emergency department between March 2020 and April 2022. The prognostic value of early warning scores in predicting in-hospital mortality was assessed using the receiver operating characteristic (ROC) curve. Patients' median age was 59 years, and 52.33% were male. Among all the EWS we assessed, REMS had the highest overall accuracy (AUC 0.84 (0.83-0.85)) and the highest NPV (97.4%). REMS was the most accurate scoring system, characterised by the highest discriminative power and negative predictive value compared to the other analysed scoring systems. Incorporating these tools into clinical practice in a hospital emergency department could provide more effective assessment of mortality and, consequently, avoid delayed medical assistance.

Genetic predisposition to differentiated thyroid cancer in the Polish population.

INTRODUCTION: Genome sequencing technologies reveal molecular mechanisms of differentiated thyroid cancer (DTC). Unlike somatic mutation analysis from thyroidectomy samples, germline mutations showing genetic susceptibility to DTC are less understood. OBJECTIVES: The study aimed to assess the prevalence of germline mutations predisposing to DTC in a cohort of Polish individuals based on their whole genome sequencing data. PATIENTS AND METHODS: We analyzed sequencing data from 1076 unrelated individuals totaling over 1018 billion read pairs and yielding an average 35.26 × read depth per genome, released openly for academic and clinical research as the Thousand Polish Genomes database (https://1000polishgenomes.com). The list of genes chosen for further analysis was based on the review of previous studies. RESULTS: The cohort contained 104 variants located within the coding and noncoding DNA sequences of 90 genes selected by ClinVar classification as pathogenic and potentially pathogenic. The frequency of variants in the Polish cohort was compared with the frequency estimated for the non­Finnish European population obtained from the gnomAD database (gnomad.broadinstitute.org). Significant differences in variant frequency were found for the APC, ARSB, ATM, BRCA1, CHEK2, DICER1, GPD1L, INSR, KCNJ10, MYH9, PALB2, PLCB1, PLEKHG5, PTEN, RET, SEC23B, SERPINA1, SLC26A4, SMAD3, STK11, TERT, TOE1, and WRN genes. CONCLUSIONS: Even though the Polish population is genetically similar to the other European populations, there are significant differences in variant frequencies contributing to the disease development and progression, such as those in the RET, CHEK2, BRCA1, SLC26A4, or TERT genes. Further studies are needed to identify genomic variants associated directly with DTC.

SARS-CoV-2 infection activating a novel variant of the NOTCH3 gene and subsequently causing development of CADASIL.

Introduction: In the following study we describe the diagnostic process and further case analysis of a 30-year-old woman admitted with typical COVID-19 symptoms, who subsequently developed additional symptoms suggesting cerebral autosomal dominant arteriopathy with sub-cortical infarcts and leukoencephalopathy (CADASIL). Material and methods: Other than the standard diagnostic procedures, whole genome sequencing (WGS) was used, which led to following findings. A new variant of the NOTCH3 gene, which led to CADASIL-like symptoms, was found, and it had been most likely activated by the SARS-CoV-2 infection. This novel variant in NOTCH3 has not been found in existing databases and has never been mentioned in research concerning CADASIL before. Results: Furthermore, after subjecting the patient's close relatives to WGS it was found that no other examined person demonstrated the same genetic mutation. Conclusions: It seems therefore that the new variant of NOTCH3 is of de novo origin in the patient's genome. Additionally, the relatively early onset of CADASIL and the unexpectedly severe COVID-19 infection suggest that the two occurred simultaneously: the infection with SARS-CoV-2 accelerated development of CADASIL symptoms and the unusual variant of the NOTCH3 gene contributed to the more severe course of COVID-19.

The Spectrum of the Heterozygous Effect in Biallelic Mendelian Diseases-The Symptomatic Heterozygote Issue.

Heterozygous carriers of pathogenic/likely pathogenic variants in autosomal recessive disorders seem to be asymptomatic. However, in recent years, an increasing number of case reports have suggested that mild and unspecific symptoms can occur in some heterozygotes, as symptomatic heterozygotes have been identified across different disease types, including neurological, neuromuscular, hematological, and pulmonary diseases. The symptoms are usually milder in heterozygotes than in biallelic variants and occur "later in life". The status of symptomatic heterozygotes as separate entities is often disputed, and alternative diagnoses are considered. Indeed, often only a thin line exists between dual, dominant, and recessive modes of inheritance and symptomatic heterozygosity. Interestingly, recent population studies have found global disease effects in heterozygous carriers of some genetic variants. What makes the few heterozygotes symptomatic, while the majority show no symptoms? The molecular basis of this phenomenon is still unknown. Possible explanations include undiscovered deep-splicing variants, genetic and environmental modifiers, digenic/oligogenic inheritance, skewed methylation patterns, and mutational burden. Symptomatic heterozygotes are rarely reported in the literature, mainly because most did not undergo the complete diagnostic procedure, so alternative diagnoses could not be conclusively excluded. However, despite the increasing accessibility to high-throughput technologies, there still seems to be a small group of patients with mild symptoms and just one variant of autosomes in biallelic diseases. Here, we present some examples, the current state of knowledge, and possible explanations for this phenomenon, and thus argue against the existing dominant/recessive classification.

The cancer-risk variant frequency among Polish population reported by the first national whole-genome sequencing study.

Introduction: Population-based cancer screening has raised many controversies in recent years, not only regarding the costs but also regarding the ethical nature and issues related to variant interpretation. Nowadays, genetic cancer screening standards are different in every country and usually encompass only individuals with a personal or family history of relevant cancer. Methods: Here we performed a broad genetic screening for cancer-related rare germline variants on population data from the Thousand Polish Genomes database based on 1076 Polish unrelated individuals that underwent whole genome sequencing (WGS). Results: We identified 19 551 rare variants in 806 genes related to oncological diseases, among them 89% have been located in non-coding regions. The combined BRCA1/BRCA2 pathogenic/likely pathogenic according to ClinVar allele frequency in the unselected population of 1076 Poles was 0.42%, corresponding to nine carriers. Discussion: Altogether, on the population level, we found especially problematic the assessment of the pathogenicity of variants and the relation of ACMG guidelines to the population frequency. Some of the variants may be overinterpreted as disease-causing due to their rarity or lack of annotation in the databases. On the other hand, some relevant variants may have been overseen given that there is little pooled population whole genome data on oncology. Before population WGS screening will become a standard, further studies are needed to assess the frequency of the variants suspected to be pathogenic on the population level and with reporting of likely benign variants.

WGS Data Collections: How Do Genomic Databases Transform Medicine?

As a scientific community we assumed that exome sequencing will elucidate the basis of most heritable diseases. However, it turned out it was not the case; therefore, attention has been increasingly focused on the non-coding sequences that encompass 98% of the genome and may play an important regulatory function. The first WGS-based datasets have already been released including underrepresented populations. Although many databases contain pooled data from several cohorts, recently the importance of local databases has been highlighted. Genomic databases are not only collecting data but may also contribute to better diagnostics and therapies. They may find applications in population studies, rare diseases, oncology, pharmacogenetics, and infectious and inflammatory diseases. Further data may be analysed with Al technologies and in the context of other omics data. To exemplify their utility, we put a highlight on the Polish genome database and its practical application.

Follicular Thyroid Adenoma and Follicular Thyroid Carcinoma-A Common or Distinct Background? Loss of Heterozygosity in Comprehensive Microarray Study.

Pre- and postsurgical differentiation between follicular thyroid adenoma (FTA) and follicular thyroid cancer (FTC) represents a significant diagnostic challenge. Furthermore, it remains unclear whether they share a common or distinct background and what the mechanisms underlying follicular thyroid lesions malignancy are. The study aimed to compare FTA and FTC by the comprehensive microarray and to identify recurrent regions of loss of heterozygosity (LOH). We analyzed formalin-fixed paraffin-embedded (FFPE) samples acquired from 32 Caucasian patients diagnosed with FTA (16) and FTC (16). We used the OncoScan™ microarray assay (Affymetrix, USA), using highly multiplexed molecular inversion probes for single nucleotide polymorphism (SNP). The total number of LOH was higher in FTC compared with FTA (18 vs. 15). The most common LOH present in 21 cases, in both FTA (10 cases) and FTC (11 cases), was 16p12.1, which encompasses many cancer-related genes, such as TP53, and was followed by 3p21.31. The only LOH present exclusively in FTA patients (56% vs. 0%) was 11p11.2-p11.12. The alteration which tended to be detected more often in FTC (6 vs. 1 in FTA) was 12q24.11-q24.13 overlapping FOXN4, MYL2, PTPN11 genes. FTA and FTC may share a common genetic background, even though differentiating rearrangements may also be detected.

Genetics, Genomics and Emerging Molecular Therapies of Pancreatic Cancer.

The number of cases of pancreatic cancers in 2019 in Poland was 3852 (approx. 2% of all cancers). The course of the disease is very fast, and the average survival time from the diagnosis is 6 months. Only <2% of patients live for 5 years from the diagnosis, 8% live for 2 years, and almost half live for only about 3 months. A family predisposition to pancreatic cancer occurs in about 10% of cases. Several oncogenes in which somatic changes lead to the development of tumours, including genes BRCA1/2 and PALB2, TP53, CDKN2A, SMAD4, MLL3, TGFBR2, ARID1A and SF3B1, are involved in pancreatic cancer. Between 4% and 10% of individuals with pancreatic cancer will have a mutation in one of these genes. Six percent of patients with pancreatic cancer have NTRK pathogenic fusion. The pathogenesis of pancreatic cancer can in many cases be characterised by homologous recombination deficiency (HRD)-cell inability to effectively repair DNA. It is estimated that from 24% to as many as 44% of pancreatic cancers show HRD. The most common cause of HRD are inactivating mutations in the genes regulating this DNA repair system, mainly BRCA1 and BRCA2, but also PALB2, RAD51C and several dozen others.

Better safe than sorry-Whole-genome sequencing indicates that missense variants are significant in susceptibility to COVID-19.

Undoubtedly, genetic factors play an important role in susceptibility and resistance to COVID-19. In this study, we conducted the GWAS analysis. Out of 15,489,173 SNPs, we identified 18,191 significant SNPs for severe and 11,799 SNPs for resistant phenotype, showing that a great number of loci were significant in different COVID-19 representations. The majority of variants were synonymous (60.56% for severe, 58.46% for resistant phenotype) or located in introns (55.77% for severe, 59.83% for resistant phenotype). We identified the most significant SNPs for a severe outcome (in AJAP1 intron) and for COVID resistance (in FIG4 intron). We found no missense variants with a potential causal function on resistance to COVID-19; however, two missense variants were determined as significant a severe phenotype (in PM20D1 and LRP4 exons). None of the aforementioned SNPs and missense variants found in this study have been previously associated with COVID-19.

Population WGS-based spinal muscular atrophy carrier screening in a cohort of 1076 healthy Polish individuals.

Spinal muscular atrophy is a severe neuromuscular disorder with an autosomal recessive inheritance pattern. The disease-causing gene is SMN1, and its paralogue, SMN2, is a disease course modifier. Both genes SMN1 and SMN2 show over 99.9% sequence identity and a high rate of crossing over in the genomic region. Due to this reason, SMN1/SMN2 is usually excluded from the whole-genome sequencing (WGS) analysis and investigated with traditional methods, such as MLPA and qPCR. Recently, novel bioinformatic algorithms dedicated to analyzing this particular genomic region have been developed. Here, we analyze the SMN1/SMN2 genomic region with a dedicated program, SMNCopyNumberCaller. We report a similar prevalence of SMN1 gene deletion carrier status (1 per 41 people) to published data from the Polish population (1 per 35 people). Additionally, SMNCopyNumberCaller can identify SMN2 CNVs and SMN2Δ7-8 present in 153 healthy Polish individuals. Two other programs for the CNV analysis in standard genomic regions were not able to provide reliable results. Using WGS-based tools for SMN1/2 genomic region analysis is not only an efficient method in terms of time but will also enable more complex analysis such screening for markers related with a silent carrier status and identification of further genetic modifiers. Although still an experimental method, soon WGS-based SMN1/SMN2 carrier identification may become a standard method for patients screened with WGS for other purposes.

Gene Variants Related to Cardiovascular and Pulmonary Diseases May Correlate with Severe Outcome of COVID-19.

Background: Severe outcomes of COVID-19 account for up to 15% of all cases. The study aims to check if any gene variants related to cardiovascular (CVD) and pulmonary diseases (PD) are correlated with a severe outcome of COVID-19 in a Polish cohort of COVID-19 patients. Methods: In this study, a subset of 747 samples from unrelated individuals collected across Poland in 2020 and 2021 was used and whole-genome sequencing was performed. Results: The GWAS analysis of SNPs and short indels located in genes related to CVD identified one variant significant in COVID-19 severe outcome in the HADHA gene, while for the PD gene panel, we found two significant variants in the DRC1 gene. In this study, both potentially protective and risk variants were identified, of which variants in the HADHA gene deserve the most attention. Conclusions: This is the first study reporting the association between the HADHA and DRC1 genetic variants and COVID-19 severe outcome based on the cohort WGS analysis. Although all the identified variants are localised in introns, they may be correlated and therefore inherited along with other risk variants, potentially causative to severe outcome of COVID-19 but not discovered yet.

Beyond GWAS-Could Genetic Differentiation within the Allograft Rejection Pathway Shape Natural Immunity to COVID-19?

COVID-19 infections pose a serious global health concern so it is crucial to identify the biomarkers for the susceptibility to and resistance against this disease that could help in a rapid risk assessment and reliable decisions being made on patients' treatment and their potential hospitalisation. Several studies investigated the factors associated with severe COVID-19 outcomes that can be either environmental, population based, or genetic. It was demonstrated that the genetics of the host plays an important role in the various immune responses and, therefore, there are different clinical presentations of COVID-19 infection. In this study, we aimed to use variant descriptive statistics from GWAS (Genome-Wide Association Study) and variant genomic annotations to identify metabolic pathways that are associated with a severe COVID-19 infection as well as pathways related to resistance to COVID-19. For this purpose, we applied a custom-designed mixed linear model implemented into custom-written software. Our analysis of more than 12.5 million SNPs did not indicate any pathway that was significant for a severe COVID-19 infection. However, the Allograft rejection pathway (hsa05330) was significant (p = 0.01087) for resistance to the infection. The majority of the 27 SNP marking genes constituting the Allograft rejection pathway were located on chromosome 6 (19 SNPs) and the remainder were mapped to chromosomes 2, 3, 10, 12, 20, and X. This pathway comprises several immune system components crucial for the self versus non-self recognition, but also the components of antiviral immunity. Our study demonstrated that not only single variants are important for resistance to COVID-19, but also the cumulative impact of several SNPs within the same pathway matters.

CAPN3 c.1746-20C>G variant is hypomorphic for LGMD R1 calpain 3-related.

The investigated intronic CAPN3 variant NM_000070.3:c.1746-20C>G occurs in the Central and Eastern Europe with a frequency of >1% and there are conflicting interpretations on its pathogenicity. We collected data on 14 patients carrying the CAPN3 c.1746-20C>G variant in trans position with another CAPN3 pathogenic/likely pathogenic variant. The patients compound heterozygous for the CAPN3 c.1746-20C>G variant presented a phenotype consistent with calpainopathy of mild/medium severity. This variant is most frequent in the North/West regions of Russia and may originate from that area. Molecular studies revealed that different splicing isoforms are produced in the muscle. We hypothesize that c.1746-20C>G is a hypomorphic variant with a reduction of RNA and protein expression and only individuals having a higher ratio of abnormal isoforms are affected. Reclassification of the CAPN3 variant c.1746-20C>G from variant with a conflicting interpretation of pathogenicity to hypomorphic variant explains many unidentified cases of limb girdle muscular dystrophy R1 calpain 3-related in Eastern and Central Europe.

The Thousand Polish Genomes-A Database of Polish Variant Allele Frequencies.

Although Slavic populations account for over 4.5% of world inhabitants, no centralised, open-source reference database of genetic variation of any Slavic population exists to date. Such data are crucial for clinical genetics, biomedical research, as well as archeological and historical studies. The Polish population, which is homogenous and sedentary in its nature but influenced by many migrations of the past, is unique and could serve as a genetic reference for the Slavic nations. In this study, we analysed whole genomes of 1222 Poles to identify and genotype a wide spectrum of genomic variation, such as small and structural variants, runs of homozygosity, mitochondrial haplogroups, and de novo variants. Common variant analyses showed that the Polish cohort is highly homogenous and shares ancestry with other European populations. In rare variant analyses, we identified 32 autosomal-recessive genes with significantly different frequencies of pathogenic alleles in the Polish population as compared to the non-Finish Europeans, including C2, TGM5, NUP93, C19orf12, and PROP1. The allele frequencies for small and structural variants, calculated for 1076 unrelated individuals, are released publicly as The Thousand Polish Genomes database, and will contribute to the worldwide genomic resources available to researchers and clinicians.

Transcriptional and post-transcriptional regulation of checkpoint genes on the tumour side of the immunological synapse.

Cancer is a disease of the genome, therefore, its development has a clear Mendelian component, demonstrated by well-studied genes such as BRCA1 and BRCA2 in breast cancer risk. However, it is known that a single genetic variant is not enough for cancer to develop leading to the theory of multistage carcinogenesis. In many cases, it is a sequence of events, acquired somatic mutations, or simply polygenic components with strong epigenetic effects, such as in the case of brain tumours. The expression of many genes is the product of the complex interplay between several factors, including the organism's genotype (in most cases Mendelian-inherited), genetic instability, epigenetic factors (non-Mendelian-inherited) as well as the immune response of the host, to name just a few. In recent years the importance of the immune system has been elevated, especially in the light of the immune checkpoint genes discovery and the subsequent development of their inhibitors. As the expression of these genes normally suppresses self-immunoreactivity, their expression by tumour cells prevents the elimination of the tumour by the immune system. These discoveries led to the rapid growth of the field of immuno-oncology that offers new possibilities of long-lasting and effective treatment options. Here we discuss the recent advances in the understanding of the key mechanisms controlling the expression of immune checkpoint genes in tumour cells.

Challenges of the Immunotherapy: Perspectives and Limitations of the Immune Checkpoint Inhibitor Treatment.

Immunotherapy is a quickly developing type of treatment and the future of therapy in oncology. This paper is a review of recent findings in the field of immunotherapy with an emphasis on immune checkpoint inhibitors. The challenges that immunotherapy might face in near future, such as primary and acquired resistance and the irAEs, are described in this article, as well as the perspectives such as identification of environmental modifiers of immunity and development of anti-cancer vaccines and combined therapies. There are multiple factors that may be responsible for immunoresistance, such as genomic factors, factors related to the immune system cells or to the cancer microenvironment, factors emerging from the host cells, as well as other factors such as advanced age, biological sex, diet, many hormones, existing comorbidities, and the gut microbiome.

Comprehensive Analysis of Correlations in the Expression of miRNA Genes and Immune Checkpoint Genes in Bladder Cancer Cells.

Personalised medicine is the future and hope for many patients, including those with cancers. Early detection, as well as rapid, well-selected treatment, are key factors leading to a good prognosis. MicroRNA mediated gene regulation is a promising area of development for new diagnostic and therapeutic methods, crucial for better prospects for patients. Bladder cancer is a frequent neoplasm, with high lethality and lacking modern, advanced therapeutic modalities, such as immunotherapy. MicroRNAs are involved in bladder cancer pathogenesis, proliferation, control and response to treatment, which we summarise in this perspective in response to lack of recent review publications in this field. We further performed a correlation-based analysis of microRNA and gene expression data in bladder cancer (BLCA) TCGA dataset. We identified 27 microRNAs hits with opposite expression profiles to genes involved in immune response in bladder cancer, and 24 microRNAs hits with similar expression profiles. We discuss previous studies linking the functions of these microRNAs to bladder cancer and assess if they are good candidates for personalised medicine therapeutics and diagnostics. The discussed functions include regulation of gene expression, interplay with transcription factors, response to treatment, apoptosis, cell proliferation and angiogenesis, initiation and development of cancer, genome instability and tumour-associated inflammatory reaction.

Associations between T cell infiltration, T cell receptor clonality, histology and recurrence in renal cell carcinoma.

Renal cell carcinoma (RCC) is comprised of clear-cell (ccRCC) and non-clear-cell (nccRCC) tumors. Despite definitive surgical resection in localized disease, recurrence often occurs. A commercial method based on a multiplex polymerase chain reaction (PCR) assay exclusively targets rearranged T cell receptor (TCR) genes to generate high-throughput sequencing-based data, allowing characterization of the immune repertoire within tumors. In this study we performed a retrospective analysis on archived tumor samples from patients with recurring versus non-recurring T3 ccRCC and on samples from early nccRCC versus ccRCC. Following genomic DNA extraction and multiplex PCR, the fraction of T cells within tumors, the number of unique receptors ('richness') and their relative abundances ('clonality') were calculated. Statistical significance and correlations were calculated using Student's t-test and Spearman's rho, respectively. Average fraction and clonality of T cells in tumors from non-recurring patients was 2.5- and 4.3-fold higher than in recurring patients (P = 0.025 and P = 0.043, respectively). A significant positive correlation was found between T cell fraction and clonality (Spearman's rho = 0.78, P = 0.008). The average fraction of T cells in ccRCC tumors was 2.8-fold higher than in nccRCC tumors (P = 0.015). Clonality and estimated richness were similar between ccRCC and nccRCC tumors. In summary, recurrence of ccRCC is associated with a lower fraction and clonality of T cells within tumors; nccRCC tumors are more 'deserted' than ccRCC, but similar in their ability to generate a clonal T cell repertoire. Our work suggests associations between the characteristics of T cell infiltrate, histology and tumor recurrence.

Checkpoint Genes at the Cancer Side of the Immunological Synapse in Bladder Cancer.

Immune checkpoint inhibitors have revolutionized cancer therapy, but not all cancers respond to the currently available drugs, and even within cancers considered responsive to such modality, response rates range between 15 and 40%, depending on the cancer type, the line of treatment, and yet unknown clinical/molecular factors. Coordinated expression of checkpoint proteins was shown to occur on T cells, probably allowing fine-tuning of the signal transmitted to the cell. We performed a bioinformatic analysis of the expression of putative checkpoint mRNAs at the cancer side of the immunological synapse from the bladder cancer tumorgenome atlas (TCGA) database. Fifteen mRNAs, corresponding to both coinhibitory and costimulatory checkpoints, were shown to be expressed above a designated threshold. Of these, seven mRNAs were found to be coexpressed: CD277, PD-1L, CD48, CD86, galectin-9, TNFRSF14 (HVEM), and CD40. The expression of 2 of these mRNAs-BTN3A1 (CD277) and TNFRSF14 (HVEM)-was positively correlated with overall survival in the TCGA database. All these seven mRNA share putative binding sites of a few transcription factors (TFs). Of these, the expression of the TF BACH-2 was positively correlated with the expression of checkpoint mRNAs from the network. This suggests a joint transcriptional regulation on the expression of checkpoint mRNAs at the bladder tumor side of the immunological synapse.