Comprehensive genomic analysis of primary malignant melanoma of the esophagus reveals similar genetic patterns compared with epithelium-associated melanomas.

Primary malignant melanoma of the esophagus (PMME) is an exceedingly rare disease with a poor prognosis. The etiology of PMME remains largely unknown and genetic characteristics are yet to be clarified, essential for identifying potential therapeutic targets and defining treatment guidelines. Here, we performed whole-exome sequencing on 47 formalin-fixed paraffin-embedded specimens from 18 patients with PMME, including 23 tumor samples, 6 metastatic lymph nodes, and 18 tumor-adjacent normal tissues. The genomic features of PMME were comprehensively characterized, and comparative genomic analysis was further performed between these specimens and 398 skin cutaneous melanomas (SKCM), 67 non-esophagus mucosal melanomas (NEMM), and 79 uveal melanomas (UVM). In the PMME cohort, recurrently mutated driver genes, such as MUC16, RANBP2, NRAS, TP53, PTPRT, NF1, MUC4, KMT2C, and BRAF, were identified. All RANBP2 mutations were putatively deleterious, and most affected samples had multipoint mutations. Furthermore, RANBP2 showed parallel evolution by multiregional analysis. Whole-genome doubling was an early truncal event that occurred before most driver mutations, except for in TP53. An ultraviolet radiation-related mutational signature, SBS38, was identified as specific to epithelial melanomas and could predict inferior survival outcomes in both PMME and SKCM patients. Comparing the mutational and copy number landscapes between PMME and other subtypes of melanoma revealed that PMME has a similar genomic pattern and biological characteristics to SKCM. In summary, we comprehensively defined the key genomic aberrations and mutational processes driving PMME and suggested for the first time that PMME may share similar genomic patterns with SKCM; therefore, patients with rare melanomas, such as PMME, may benefit from the current treatment used for common cutaneous melanoma.

Integrated Transcriptomic Landscape and Deep Learning Based Survival Prediction in Uterine Sarcomas.

Purpose: The genomic characteristics of uterine sarcomas have not been fully elucidated. This study aimed to explore the genomic landscape of the USs. Materials and Methods: Comprehensive genomic analysis through RNA-sequencing was conducted. Gene fusion, differentially expressed genes (DEGs), signaling pathway enrichment, immune cell infiltration, and prognosis were analyzed. A deep learning model was constructed to predict the survival of US patients. Results: A total of 71 US samples were examined, including 47 endometrial stromal sarcomas (ESS), 18 uterine leiomyosarcomas (uLMS), 3 adenosarcomas, 2 carcinosarcomas, and 1 uterine tumor resembling an ovarian sex-cord tumor (UTROSCT). ESS (including high-grade ESS and low-grade ESS) and uLMS showed distinct gene fusion signatures; a novel gene fusion site, MRPS18A - PDC-AS1 could be a potential diagnostic marker for the pathology differential diagnosis of uLMS and ESS; 797 and 477 uDEGs were identified in the ESS vs. uLMS and HGESS vs. LGESS groups, respectively. The uDEGs were enriched in multiple pathways. Fifteen genes including LAMB4 were confirmed with prognostic value in USs; immune infiltration analysis revealed the prognositic value of myeloid dendritic cells, plasmacytoid dendritic cells, natural killer cells, macrophage M1, monocytes and hematopoietic stem cells in USs; the deep learning model named MMN-MIL showed satisfactory performance in predicting the survival of US patients, with the area under the receiver operating curve curve reached 0.909 and accuracy achieved 0.804. Conclusion: USs harbored distinct gene fusion characteristics and gene expression features between HGESS, LGESS, and uLMS. The MMN-MIL model could effectively predict the survival of US patients.

Next-generation sequencing shows the genomic features of ovarian clear cell cancer and compares the genetic architectures of high-grade serous ovarian cancer and clear cell carcinoma in ovarian and endometrial tissues.

Ovarian clear cell carcinoma (OCCC) is a special histological type of epithelial ovarian cancer (EOC) that is not derived from epithelial cells of the ovarian or fallopian tube as the most common type of ovarian cancer, high-grade serous ovarian carcinoma (HGSOC), but is closely related to endometriosis and similar to endometrial clear cell carcinoma (ECCC) at morphologic and phenotypic features. However, limited data was shown in OCCC genomic features and compared with that in OCCC, HGSOC and ECCC. Herein, we utilized next-generation sequencing analysis of a panel of 1,021 genes to profile the mutational alterations in 34 OCCC and compared them to those from HGSOC (402 cases) and ECCC (30 cases). In result, the ARID1A and PIK3CA are high-frequency mutations of OCCC. Clonal architectures showed that all the mutations of genes occur in the later stage in the OCCC progress, whereas KRAS mutation is the earlier event compared with mutation of ARID1A or PIK3CA, which usually occurs in a group of ARID1A or PIK3CA mutations. The mutation frequency of main driver genes is similar between OCCC and ECCC, while TP53 is the main mutation in HGSOC and ECCC. Shared mutational signatures between OCCC and ECCC tissues with commonly observed a C>T change indicated a common carcinogens-exposed between these two carcinomas, but HGSOC and ECCC have common and distinct mutational signatures across cohorts respectively. In addition, we identified some novel CNV gains in NF1, ASXL1, TCF7L2, CREBBP and LRP1B and loss in ATM, FANCM, RB1 and FLT in OCCC. Our study offered a new perspective for OCCC tumorigenesis from two organs, the ovary and uterus, at genomic architectures and revealed novel CNV events for helping to provide theoretical support for OCCC treatment.

A potential treatment option for transformed small-cell lung cancer on PD-L1 inhibitor-based combination therapy improved survival.

OBJECTIVES: Transformed small-cell lung cancer (T-SCLC) has an extremely poor prognosis, and no remedies based on immunotherapy have been evaluated among T-SCLC patients. We retrospectively analysed the efficacy and safety of combining atezolizumab with chemotherapy for T-SCLC. METHODS: Forty-seven patients harbouring EGFR mutations who developed T-SCLC were enrolled. Eleven patients who used immunotherapy were defined as the I/O group, and the remaining 36 were defined as the Non-I/O group. Clinical characteristics, pathological data, and survival outcomes were collected. RNA sequencing and whole-exome sequencing (WES) were performed for in-depth analysis. RESULTS: All patients received at least one line of EGFR-TKI before rebiopsy to confirm T-SCLC. Nine patients received atezolizumab-bevacizumab-carboplatin-paclitaxel (albumin-bound) (ABCP), and the remaining 2 received atezolizumab-etoposide-carboplatin (ECT) in the I/O group. The objective response rate was 73 % (8/11). The median progression-free survival (mPFS) of T-SCLC on post-transformation therapy with I/O group and Non-I/O group was 5.1 m and 4.1 m, respectively. The median post-T-SCLC overall survival of the I/O group was significantly longer than that Non-I/O group (20.2 m vs 7.9 m, P ＜ 0.01). T-SCLC harbouring EGFR L858R tended to be longer than EGFR 19del (mPFS: not reached vs 3.7 m, P = 0.11). Positive PD-L1 status was also associated with PFS benefits (mPFS: 6.0 m vs 3.7 m, P = 0.20). Furthermore, RNA sequencing revealed that expression of SFTPA1 is significantly higher in the durable clinical benefit group. WES showed that STC2 mutation is more frequently observed at the time-point immunotherapy acquired resistance. Combination therapy based on a PD-L1 inhibitor was well tolerated, and the safety profile was consistent with previously reported studies. CONCLUSION: Our study first demonstrated that a PD-L1 inhibitor combined with chemotherapy ± bevacizumab could be a potential safe option for specific SCLC-transformed patients. Subsequent studies with more patients are essential to verify the efficacy and potential biomarkers.

Identification of heritable rare variants associated with early-stage lung adenocarcinoma risk.

Background: In East Asia, the number of patients with adenocarcinoma, especially those presenting with ground-glass nodules (GGNs), is gradually increasing. Family aggregation of pulmonary GGNs is not uncommon; however, genetic predisposition in these patients remains poorly understood and identification of genes involved in the cause of these early-stage lung cancers might contribute to understanding of the underlying mechanisms and potential prevention strategies. Methods: Fifty patients with early-stage lung adenocarcinoma (LUAD) presenting as GGNs and a first-degree family history of lung cancer (FHLC) from 34 independent families were enrolled into this study. Germline mutations of these patients were analyzed with whole exome sequencing (WES) and compared with age- and sex-matched 39 patients with sporadic lung cancer and 689 local healthy people. We used a stepwise variant filtering strategy, gene-based burden testing, and enrichment analysis to investigate rare but potentially pathogenic heritable mutations. Somatic tumor mutations were analyzed to consolidate germline findings. Results: In total, 1,571 single nucleotide variants (SNVs) and 238 frameshifts with a minor allele frequency (MAF) <0.01, which were rare, recurrent, and potentially damaging candidates, were finally identified through the filtering in the GGN cohort. Pathway analysis showed the extracellular matrix to be the top dysregulated pathway. Gene-based burden testing of these highly disruptive risk-conferring heritable variants showed that MSH5 [odds ratio (OR), 9.28, 95% confidence interval (CI): 2.49-35.87], MMP9 (OR, 8.11, 95% CI: 2.22-28.43), and CYP2D6 (OR, 8.09, 95% CI: 2.68-24.92) were significantly enriched in our cohort (P<0.05). The number of rare damaging germline variants in non-smoking patients was significantly higher than that of smoking-affected patients (Spearman's ρ=-0.39, P=0.02). Conclusions: Heritable, potentially deleterious, and rare candidate variants of MSH5, MMP9 and CYP2D6 were significantly associated with early-stage LUAD presenting with GGNs. Nonsmoking patients likely have a higher genetic predisposition to this type of cancer than smoking-affected patients. These results have extended our understanding of the underlying mechanisms of early-stage LUAD.

Homologous recombination deficiency (HRD) can predict the therapeutic outcomes of immuno-neoadjuvant therapy in NSCLC patients.

BACKGROUND: Neoadjuvant immunotherapy is emerging as novel effective intervention in lung cancer, but study to unearth effective surrogates indicating its therapeutic outcomes is limited. We investigated the genetic changes between non-small cell lung cancer (NSCLC) patients with varied response to neoadjuvant immunotherapy and discovered highly potential biomarkers with indicative capability in predicting outcomes. METHODS: In this study, 3 adenocarcinoma and 11 squamous cell carcinoma NSCLC patients were treated by neoadjuvant immunotherapy with variated regimens followed by surgical resection. Treatment-naive FFPE or fresh tissues and blood samples were subjected to whole-exome sequencing (WES). Genetic alternations were compared between differently-responded patients. Findings were further validated in multiple public cohorts. RESULTS: DNA damage repair (DDR)-related InDel signatures and DDR-related gene mutations were enriched in better-responded patients, i.e., major pathological response (MPR) group. Besides, MPR patients exhibited provoked genome instability and unique homologous recombination deficiency (HRD) events. By further inspecting alternation status of homology-dependent recombination (HR) pathway genes, the clonal alternations were exclusively enriched in MPR group. Additionally, associations between HR gene alternations, percentage of viable tumor cells and HRD event were identified, which orchestrated tumor mutational burden (TMB), mutational intratumor heterogeneity (ITH), somatic copy number alteration (SCNA) ITH and clonal neoantigen load in patients. Validations in public cohorts further supported the generality of our findings. CONCLUSIONS: We reported for the first time the association between HRD event and enhanced neoadjuvant immunotherapy response in lung cancer. The power of HRD event in patient therapeutic stratification persisted in multifaceted public cohorts. We propose that HR pathway gene status could serve as novel and additional indicators guiding immune-neoadjuvant and immunotherapy treatment decisions for NSCLC patients.

Whole exome and transcriptome sequencing reveal clonal evolution and exhibit immune-related features in metastatic colorectal tumors.

Liver is the most common site where metastatic lesions of colorectal cancer (CRC) arise. Although researches have shown mutations in driver genes, copy number variations (CNV) and alterations in relevant signaling pathways promoted the tumor evolution and immune escape during colorectal liver metastasis (CLM), the underlying mechanism remains largely elusive. Tumor and matched metastatic tissues were collected from 16 patients diagnosed with colorectal cancer and subjected to whole-exome sequencing (WES) and RNA sequencing (RNA-seq) for studying colorectal cancer clonal evolution and immune escape during CLM. Shared somatic mutations between primary and metastatic tissues with a commonly observed subclonal-clonal (S-C) changing pattern indicated a common clonal origin between two lesions. The recurrent mutations with S-C changing pattern included those in KRAS, SYNE1, CACNA1H, PCLO, FBXL2, and DNAH11. The main CNV events underwent clonal-clonal evolution (20q amplification (amp), 17p deletion (del), 18q del and 8p del), subclonal-clonal evolution (8q amp, 13q amp, 8p del) and metastasis-specific evolution (8q amp) during the process of CLM. In addition, we revealed a potential mechanism of tumor cell immune escape by analyzing human leukocytes antigens (HLA) related clonal neoantigens and immune cell components in CLM. Our study proposed a novel liver metastasis-related evolutionary process in colorectal cancer and emphasized the theory of neo-immune escape in colorectal liver metastasis.

Construction and Validation of a Platinum Sensitivity Predictive Model With Multiple Genomic Variations for Epithelial Ovarian Cancer.

Platinum-based chemotherapy is still the standard of care after cytoreductive surgery in the first-line treatment for epithelial ovarian cancer. This study aims to integrate novel biomarkers for predicting platinum sensitivity in EOC after initial cytoreductive surgery precisely. To this end, 60 patients were recruited from September 2014 to October 2019. Based on the duration of progress-free survival, 44 and 16 patients were assigned to platinum-sensitive and platinum-resistant group, respectively. Next generation sequencing was performed to dissect the genomic features of ovarian tumors obtained from surgery. Multiple genomic variations were compared between two groups, including single-nucleotide variant, single base or indel signature, loss of heterozygosity (LOH), whole-genome duplication (WGD), and others. The results demonstrated that patients with characteristics including positive SBS10a signature (p < 0.05), or FAM175A LOH (p < 0.01), or negative WGD (p < 0.01) were significantly enriched in platinum-sensitive group. Consistently, patients with positive SBS10a signature (15.8 vs. 10.1 months, p < 0.05), or FAM175A LOH (16.5 vs. 9.2 months, p < 0.05), or negative WGD (16.5 vs. 9.1 months, p < 0.05) have significantly longer PFS than those without these genetic features. By integrating these three biomarkers, a lasso regression model was employed to train and test for all patients, with the AUC value 0.864 in platinum sensitivity prediction. Notably, 388 ovarian cancer patients from TCGA dataset were leveraged as independent validation cohort with AUC value 0.808, suggesting the favorable performance and reliability of this model.

TCR-Seq Identifies Distinct Repertoires of Distant-Metastatic and Nondistant-Metastatic Thyroid Tumors.

CONTEXT: Malignant thyroid tumor with distant metastasis is associated with poor outcome. Early detection of distant metastasis is of great clinical importance. OBJECTIVE: Thyroid tumor infiltrated with T cells can serve as a biomarker for monitoring metastasis. DESIGN: A retrospective analysis was performed of patient clinical samples collected between 2012 to 2018, using T-cell receptor sequencing (TCR-seq) for clinical exploration. SETTING: This study took place at Zhejiang Cancer Hospital. PATIENTS: Sixty-eight patients with papillary thyroid cancer (PTC) (distinct metastatic status) and 21 patients with benign nodules were enrolled. All patients had not received any treatment before surgery. MAIN OUTCOME MEASURE: The characteristics of TCRß complementary-determining region 3 (CDR3) for each patient were determined by high-throughput sequencing. RESULTS: The TCRß diversity of malignant tumors is significantly higher than benign nodules both in blood and tumor samples (Shannon index, blood, P < .01; tumor, P < .001). The malignant tumors with distant metastasis or invasiveness showed lower TCRß diversity than nonmetastasis (Shannon index, P < .01) or noninvasive (Shannon index, P < .01) malignant tumors. Analysis of the Morisita-Horn similarity index indicated significant TCRß repertoire similarity between tumor and blood in distant-metastatic patients (comparison with nonmetastasis, P < .05). According to the discrepancy of the CDR3 among patients with different clinicopathological status, the classifier was constructed to discriminate distant-metastatic individuals. A promising area under the curve value of 83.8% was obtained with the number of overlapping CDR3 clonotypes. CONCLUSION: The availability and reliability of TCR-seq render it prospective to translate these intrinsic attributes into clinical practice for monitoring distant metastasis in PTC patients.