Comprehensive identification of pathogenic variants in retinoblastoma by long- and short-read sequencing.

Retinoblastoma (RB) is the most common intraocular malignancy in childhood. The causal variants in RB are mostly characterized by previously used short-read sequencing (SRS) analysis, which has technical limitations in identifying structural variants (SVs) and phasing information. Long-read sequencing (LRS) technology has advantages over SRS in detecting SVs, phased genetic variants, and methylation. In this study, we comprehensively characterized the genetic landscape of RB using combinatorial LRS and SRS of 16 RB tumors and 16 matched blood samples. We detected a total of 232 somatic SVs, with an average of 14.5 SVs per sample across the whole genome in our cohort. We identified 20 distinct pathogenic variants disrupting RB1 gene, including three novel small variants and five somatic SVs. We found more somatic SVs were detected from LRS than SRS (140 vs. 122) in RB samples with WGS data, particularly the insertions (18 vs. 1). Furthermore, our analysis shows that, with the exception of one sample who lacked the methylation data, all samples presented biallelic inactivation of RB1 in various forms, including two cases with the biallelic hypermethylated promoter and four cases with compound heterozygous mutations which were missing in SRS analysis. By inferring relative timing of somatic events, we reveal the genetic progression that RB1 disruption early and followed by copy number changes, including amplifications of Chr2p and deletions of Chr16q, during RB tumorigenesis. Altogether, we characterize the comprehensive genetic landscape of RB, providing novel insights into the genetic alterations and mechanisms contributing to RB initiation and development. Our work also establishes a framework to analyze genomic landscape of cancers based on LRS data.

Human pangenome analysis of sequences missing from the reference genome reveals their widespread evolutionary, phenotypic, and functional roles.

Nonreference sequences (NRSs) are DNA sequences present in global populations but absent in the current human reference genome. However, the extent and functional significance of NRSs in the human genomes and populations remains unclear. Here, we de novo assembled 539 genomes from five genetically divergent human populations using long-read sequencing technology, resulting in the identification of 5.1 million NRSs. These were merged into 45284 unique NRSs, with 29.7% being novel discoveries. Among these NRSs, 38.7% were common across the five populations, and 35.6% were population specific. The use of a graph-based pangenome approach allowed for the detection of 565 transcript expression quantitative trait loci on NRSs, with 426 of these being novel findings. Moreover, 26 NRS candidates displayed evidence of adaptive selection within human populations. Genes situated in close proximity to or intersecting with these candidates may be associated with metabolism and type 2 diabetes. Genome-wide association studies revealed 14 NRSs to be significantly associated with eight phenotypes. Additionally, 154 NRSs were found to be in strong linkage disequilibrium with 258 phenotype-associated SNPs in the GWAS catalogue. Our work expands the understanding of human NRSs and provides novel insights into their functions, facilitating evolutionary and biomedical researches.

Remodeling of the immune and stromal cell compartment by PD-1 blockade in mismatch repair-deficient colorectal cancer.

Immune checkpoint inhibitor (ICI) therapy can induce complete responses in mismatch repair-deficient and microsatellite instability-high (d-MMR/MSI-H) colorectal cancers (CRCs). However, the underlying mechanism for pathological complete response (pCR) to immunotherapy has not been completely understood. We utilize single-cell RNA sequencing (scRNA-seq) to investigate the dynamics of immune and stromal cells in 19 patients with d-MMR/MSI-H CRC who received neoadjuvant PD-1 blockade. We found that in tumors with pCR, there is a concerted decrease in CD8+ Trm-mitotic, CD4+ Tregs, proinflammatory IL1B+ Mono and CCL2+ Fibroblast following treatment, while the proportions of CD8+ Tem, CD4+ Th, CD20+ B, and HLA-DRA+ Endothelial cells increase. Proinflammatory features in the tumor microenvironment mediate the persistence of residual tumors by modulating CD8+ T cells and other response-associated immune cell populations. Our study provides valuable resources and biological insights into the mechanism of successful ICI therapy and potential targets for improving treatment efficacy.

Single-cell characterization of malignant phenotypes and microenvironment alteration in retinoblastoma.

Retinoblastoma (RB) is the most common primary intraocular malignancy of childhood. It is known that the tumor microenvironment (TME) regulates tumorigenesis and metastasis. However, how the malignant progression in RB is determined by the heterogeneity of tumor cells and TME remains uncharacterized. Here, we conducted integrative single-cell transcriptome and whole-exome sequencing analysis of RB patients with detailed pathological and clinical measurements. By single-cell transcriptomic sequencing, we profiled around 70,000 cells from tumor samples of seven RB patients. We identified that the major cell types in RB were cone precursor-like (CP-like) and MKI67+ cone precursor (MKI67+ CP) cells. By integrating copy number variation (CNV) analysis, we found that RB samples had large clonal heterogeneity, where the malignant MKI67+ CP cells had significantly larger copy number changes. Enrichment analysis revealed that the conversion of CP-like to MKI67+ CP resulted in the loss of photoreceptor function and increased cell proliferation ability. The TME in RB was composed of tumor-associated macrophages (TAMs), astrocyte-like, and cancer-associated fibroblasts (CAFs). Particularly, during the invasion process, TAMs created an immunosuppressive environment, in which the proportion of TAMs decreased, M1-type macrophage was lost, and the TAMs-related immune functions were depressed. Finally, we identified that TAMs regulated tumor cells through GRN and MIF signaling pathways, while TAMs self-regulated through inhibition of CCL and GALECTIN signaling pathways during the invasion process. Altogether, our study creates a detailed transcriptomic map of RB with single-cell characterization of malignant phenotypes and provides novel molecular insights into the occurrence and progression of RB.

Structural variants in the Chinese population and their impact on phenotypes, diseases and population adaptation.

A complete characterization of genetic variation is a fundamental goal of human genome research. Long-read sequencing has improved the sensitivity of structural variant discovery. Here, we conduct the long-read sequencing-based structural variant analysis for 405 unrelated Chinese individuals, with 68 phenotypic and clinical measurements. We discover a landscape of 132,312 nonredundant structural variants, of which 45.2% are novel. The identified structural variants are of high-quality, with an estimated false discovery rate of 3.2%. The concatenated length of all the structural variants is approximately 13.2% of the human reference genome. We annotate 1,929 loss-of-function structural variants affecting the coding sequence of 1,681 genes. We discover rare deletions in HBA1/HBA2/HBB associated with anemia. Furthermore, we identify structural variants related to immunity which differentiate the northern and southern Chinese populations. Our study describes the landscape of structural variants in the Chinese population and their contribution to phenotypes and disease.

Classification of gastric cancers based on immunogenomic profiling.

BACKGROUND: Extensive evidence showed that gastric cancer (GC) is heterogeneous, and many studies have been focused on identifying GC subtypes based on genomic profiles. However, few studies have specifically explored the GC classification and predicted the classification accuracy that may help facilitate the optimal stratification of GC patients responsive to immunotherapy. METHODS: Using two publicly available GC genomics datasets, we classified GC on the basis of 797 immune related genes. Unsupervised and supervised machine learning methods were used to predict the classification. RESULTS: We identified two GC subtypes that we named as Immunity-High (IM-H) and Immunity- Low (IM-L), and demonstrated that this classification was duplicable and predictable by analyzing other datasets. IM-H subtype was characterized by greater immune cell infiltration, stronger immune activities, lower tumor purity, as well as worse survival prognosis compared to IM-L subtype. Besides the immune signatures, some cancer-associated pathways were hyperactivated in IM-H, including TGF-beta signaling pathway, Focal adhesion, Cell adhesion molecules (CAMs), Calcium signaling pathway, mTOR signaling pathway, MAPK signaling pathway and Wnt signaling pathway. In contrast, IM-L presented depressed immune signatures and increased activation of base excision repair, DNA replication, homologous recombination, non-homologous end-joining and nucleotide excision repair pathways. Furthermore, we identified subtype-specific genomic or clinical features, and subtype-specific gene ontology and networks in IM-H and IM-L subtype. CONCLUSIONS: We proposed and validated two reproducible immune molecular subtypes of GC, which has potential clinical implications for GC patient selection of immunotherapy.

Correlate the TP53 Mutation and the HRAS Mutation with Immune Signatures in Head and Neck Squamous Cell Cancer.

Although immunotherapy has emerged as an effective therapeutic strategy for various cancers including head and neck squamous cell carcinomas (HNSCCs), only a subset of patients can benefit from such therapy. Hence, it is pressing to discover predictive biomarkers for cancer immunotherapy response. TP53 and HRAS mutations frequently occur in HNSCC and correlate with a worse prognosis in HNSCC. We extensively characterized the associations of TP53 mutations and HRAS mutations with HNSCC immunity based on multiple cancer genomics datasets. We compared the enrichment levels of 20 immune signatures between TP53-mutated and TP53-wildtype HNSCCs, and between HRAS-mutated and HRAS-wildtype HNSCCs, and found that TP53 mutations were associated with depressed immune signatures while HRAS mutations were associated with enhanced immune signatures in HNSCC. Moreover, we found multiple p53- and RAS-mediated pathways showing significant correlations with HNSCC immunity. Furthermore, we demonstrated that the association between TP53 mutation and tumor immunity was independent of the human papillomavirus (HPV) infection and smoking status in HNSCC. These data suggest that p53 and RAS may play important roles in regulating HNSCC immunity and that the TP53 and HRAS mutation status could be useful biomarkers for stratifying HNSCC patients responsive to immunotherapy.

TP53 Mutations Promote Immunogenic Activity in Breast Cancer.

BACKGROUND: Although immunotherapy has recently achieved clinical successes in a variety of cancers, thus far there is no immunotherapeutic strategy for breast cancer (BC). Thus, it is important to discover biomarkers for identifying BC patients responsive to immunotherapy. TP53 mutations were often associated with worse clinical outcomes in BC whose triple-negative subtype has a high TP53 mutation rate (approximately 80%). To explore a potentially promising therapeutic option for the TP53-mutated BC subtype, we studied the association between TP53 mutations and immunogenic activity in BC. METHODS: We compared the enrichment levels of 26 immune signatures that indicated activities of diverse immune cells, functions, and pathways between TP53-mutated and TP53-wildtype BCs based on two large-scale BC multiomics datasets. Moreover, we explored the molecular cues associated with the differences in immunogenic activity between TP53-mutated and TP53-wildtype BCs. Furthermore, we performed experimental validation of the findings from bioinformatics analysis. RESULTS: Bioinformatics analysis showed that almost all analyzed immune signatures showed significantly higher enrichment levels in TP53-mutated BCs than in TP53-wildtype BCs. Moreover, in vitro experiments confirmed that mutant p53 could increase BC immunogenicity. Both computational and experimental results demonstrated that TP53 mutations could promote BC immunogenicity via regulation of the p53-mediated pathways including cell cycle, apoptosis, Wnt, Jak-STAT, NOD-like receptor, and glycolysis. Furthermore, we found that elevated immune activity was likely associated with a better survival prognosis in TP53-mutated BCs, but not necessarily in TP53-wildtype BCs. CONCLUSIONS: TP53 mutations may promote immunogenic activity in BC, suggesting that the TP53 mutation status could be a useful biomarker for stratifying BC patients responsive to immunotherapy.

ARID1A Mutations Are Associated with Increased Immune Activity in Gastrointestinal Cancer.

Because traditional treatment strategies for advanced gastrointestinal (GI) cancers often have a limited therapeutic effect, immunotherapy could be a viable approach for the therapy of advanced GI cancers, considering the recent success of immunotherapy in treating various refractory malignancies, including the DNA mismatch repair-deficient GI cancers. However, only a subset of cancer patients currently respond to immunotherapy. Thus, it is important to identify useful biomarkers for predicting cancer immunotherapy response. The tumor suppressor gene ARID1A has a high mutation rate in GI cancers and its deficiency is correlated with the microsatellite instability (MSI) genomic feature of cancer. We investigated the correlation between ARID1A mutations and tumor immunity using three GI cancer genomics datasets by the bioinformatic approach, and found that diverse antitumor immune signatures were more highly enriched in ARID1A-mutated GI cancers than in ARID1A-wildtype GI cancers. The elevated immune activity in ARID1A-mutated GI cancers was associated with the higher tumor mutation burden and lower tumor aneuploidy level, as well as a higher proportion of MSI cancers in this GI cancer subtype. Moreover, we found that ARID1A-mutated GI cancers more highly expressed PD-L1 than ARID1A-wildtype GI cancers. The elevated antitumor immune signatures and PD-L1 expression could contribute to the more active immunotherapeutic responsiveness and better survival prognosis in ARID1A-mutated GI cancers than in ARID1A-wildtype GI cancers in the immunotherapy setting, as evidenced in three cancer cohorts receiving immunotherapy. Thus, the ARID1A mutation could be a useful biomarker for identifying GI cancer patients responsive to immunotherapy.

Identification of genomic features associated with immunotherapy response in gastrointestinal cancers.

Gastrointestinal (GI) cancers prevail and account for an extremely high number of cancer deaths worldwide. The traditional treatment strategies, including surgery, chemotherapy, radiotherapy, and targeted therapy, have a limited therapeutic effect for advanced GI cancers. Recently, immunotherapy has shown promise in treating various refractory malignancies, including the GI cancers with mismatch repair deficiency (dMMR) or microsatellite instability (MSI). Thus, immunotherapy could be a promising treatment approach for GI cancers. Unfortunately, only a small proportion of GI cancer patients currently respond to immunotherapy. Therefore, it is important to discover predictive biomarkers for stratifying GI cancer patients response to immunotherapy. Certain genomic features, such as dMMR/MSI, tumor mutation burden (TMB), and tumor aneuploidy have been associated with tumor immunity and im-munotherapy response and may serve as predictive biomarkers for cancer immunotherapy. In this review, we examined the correlations between tumor immunity and three genomic features: dMMR/MSI, TMB, and tumor aneuploidy. We also explored their correlations using The Cancer Genome Atlas data and confirmed that the dMMR/MSI status, high TMB, and low tumor aneuploidy are associated with elevated tumor immunity in GI cancers. To improve the immunotherapeutic potential in GI cancers, more genetic or genomic features associated with tumor immune response need to be identified. Furthermore, it is worth exploring the combination of different immunotherapeutic methods and the combination of immunotherapy with other therapeutic approaches for cancer therapy.

Increased glycolysis correlates with elevated immune activity in tumor immune microenvironment.

BACKGROUND: Prior studies showed that tumor glycolysis and tumor immune evasion are interdependent. However, a systematic investigation of the association between tumor glycolysis and tumor immunity in various cancers remains lacking. METHODS: Using the Cancer Genome Atlas (TCGA) datasets, we explored the association between glycolytic activity and immune signatures in 14 cancer types. We also explored the associations between glycolytic activity and tumor immunity associated genetic features, including PD-L1 expression, tumor mutation burden (TMB), and tumor aneuploidy. Moreover, we performed in vitro experiments to verify some findings from bioinformatics analysis. Furthermore, we explored the association between tumor glycolytic activity and immunotherapy response. FINDINGS: Glycolytic activity was likely correlated with active immune signatures in various cancers and highly glycolytic tumors presented an immune-stimulatory tumor microenvironment. Compared to TMB and aneuploidy, glycolytic activity was a stronger and more consistent predictor for immune signatures in diverse cancers. Both computational and experimental analyses showed that glycolysis could increase PD-L1 expression in tumor. Glycolytic activity had a strong correlation with apoptosis which was a strong positive predictor for immune signatures, suggesting that apoptosis could be an important medium connecting glycolytic activity with immune activity in cancer. Finally, highly glycolytic tumors exhibited a better immunotherapy response and a favorable survival in the immunotherapy setting. INTERPRETATION: Tumor glycolysis may increase tumor immunity in diverse cancers. Glycolytic activity enhances PD-L1 expression on tumor cells and thus promotes anti-PD-1/PD-L1 immunotherapy response. Thus, the tumor glycolytic activity could be a predictive biomarker for immunotherapy response in diverse cancers. FUND: This work was supported by the China Pharmaceutical University (grant numbers 3150120001, 2632018YX01 to XW).

Classification of triple-negative breast cancers based on Immunogenomic profiling.

BACKGROUND: Abundant evidence shows that triple-negative breast cancer (TNBC) is heterogeneous, and many efforts have been devoted to identifying TNBC subtypes on the basis of genomic profiling. However, few studies have explored the classification of TNBC specifically based on immune signatures that may facilitate the optimal stratification of TNBC patients responsive to immunotherapy. METHODS: Using four publicly available TNBC genomics datasets, we classified TNBC on the basis of the immunogenomic profiling of 29 immune signatures. Unsupervised and supervised machine learning methods were used to perform the classification. RESULTS: We identified three TNBC subtypes that we named Immunity High (Immunity_H), Immunity Medium (Immunity_M), and Immunity Low (Immunity_L) and demonstrated that this classification was reliable and predictable by analyzing multiple different datasets. Immunity_H was characterized by greater immune cell infiltration and anti-tumor immune activities, as well as better survival prognosis compared to the other subtypes. Besides the immune signatures, some cancer-associated pathways were hyperactivated in Immunity_H, including apoptosis, calcium signaling, MAPK signaling, PI3K-Akt signaling, and RAS signaling. In contrast, Immunity_L presented depressed immune signatures and increased activation of cell cycle, Hippo signaling, DNA replication, mismatch repair, cell adhesion molecule binding, spliceosome, adherens junction function, pyrimidine metabolism, glycosylphosphatidylinositol (GPI)-anchor biosynthesis, and RNA polymerase pathways. Furthermore, we identified a gene co-expression subnetwork centered around five transcription factor (TF) genes (CORO1A, STAT4, BCL11B, ZNF831, and EOMES) specifically significant in the Immunity_H subtype and a subnetwork centered around two TF genes (IRF8 and SPI1) characteristic of the Immunity_L subtype. CONCLUSIONS: The identification of TNBC subtypes based on immune signatures has potential clinical implications for TNBC treatment.

Immunogenomics Analysis Reveals that TP53 Mutations Inhibit Tumor Immunity in Gastric Cancer.

Although immunotherapy continues to demonstrate efficacy in a variety of refractory cancers, currently, no any immunotherapeutic strategy is clinically used for gastric cancer (GC) except its microsatellite instable subtype. Thus, it is important to identify molecular biomarkers for predicting the responders to GC immunotherapy. TP53 mutations frequently occur in GC and are associated with unfavorable clinical outcomes in GC. We performed a comprehensive characterization of the associations between TP53 mutations and immune activities in GC based on two large-scale GC cancer genomics data. We compared expression and enrichment levels of 787 immune-related genes and 23 immune gene-sets among TP53-mutated GCs, TP53-wildtype GCs, and normal tissue, and explored the correlations between p53-mediated pathways and immune activities in GC. Strikingly, almost all analyzed immune gene-sets were significantly downregulated in enrichment levels in TP53-mutated GCs compared to TP53-wildtype GCs. These less active immune pathways and cell types in TP53-mutated GCs included 15 immune cell types and function, tumor-infiltrating lymphocytes, regulatory T cells, immune checkpoint, cytokine and cytokine receptor, human leukocyte antigen, pro-inflammatory, and parainflammation. Moreover, we identified a number of p53-mediated pathways and proteins that were significantly associated with immune activities in GC. Furthermore, we demonstrated that the TP53 mutation itself could result in the depressed immune activities in GC and other cancer types. We revealed that chromosomal instability was an important mechanism for the depressed tumor immunity in TP53-mutated cancers. Finally, we showed that immune cell infiltration and immune activities were likely positively associated with survival prognosis in GC. Our findings suggest that p53 may play an important role in activating tumor immunity in GC and other cancer types and that the TP53 mutation status could be useful in stratifying cancer patients responsive to a certain immunotherapy.

A Comprehensive Immunologic Portrait of Triple-Negative Breast Cancer.

Triple-negative breast cancer (TNBC) is a high-risk malignancy due to its high capacity for invasion and lack of targeted therapy. Immunotherapy continues to demonstrate efficacy in a variety of cancers, and thus may be a promising strategy for TNBC given the limited therapeutic options currently available for TNBC. In this study, we performed an exhaustive analysis of immunogenic signatures in TNBC based on 2 large-scale breast cancer (BC) genomic data. We compared enrichment levels of 26 immune cell activities and pathways among TNBC, non-TNBC, and normal tissue, and within TNBCs of different genotypic or phenotypic features. We found that almost all analyzed immune activities and pathways had significantly higher enrichment levels in TNBC than non-TNBC. Elevated enrichment of these immune activities and pathways was likely to be associated with better survival prognosis in TNBC. This study demonstrated that TNBC likely exhibits the strongest immunogenicity among BC subtypes, and thus warrants the immunotherapeutic option for TNBC.