Validation of computational determination of microsatellite status using whole exome sequencing data from colorectal cancer patients

Clonal origin and genomic diversity in Lynch syndrome-associated endometrial cancer with multiple synchronous tumors: Identification of the pathogenicity of MLH1 p.L582H

Lynch syndrome-associated endometrial cancer patients often present multiple synchronous tumors and this assessment can affect treatment strategies. We present a case of a 27-year-old woman with tumors in the uterine corpus, cervix, and ovaries who was diagnosed with endometrial cancer and exhibited cervical invasion and ovarian metastasis. Her family history suggested Lynch syndrome, and genetic testing identified a variant of uncertain significance, MLH1 p.L582H. We conducted immunohistochemical staining, microsatellite instability analysis, and Sanger sequencing for Lynch syndrome-associated cancers in three generations of the family and identified consistent MLH1 loss. Whole-exome sequencing for the corpus, cervical, and ovarian tumors of the proband identified a copy-neutral loss of heterozygosity (LOH) occurring at the MLH1 position in all tumors. This indicated that the germline variant and the copy-neutral LOH led to biallelic loss of MLH1 and was the cause of cancer initiation. All tumors shared a portion of somatic mutations with high mutant allele frequencies, suggesting a common clonal origin. There were no mutations shared only between the cervix and ovary samples. The profiles of mutant allele frequencies shared between the corpus and cervix or ovary indicated that two different subclones originating from the corpus independently metastasized to the cervix or ovary. Additionally, all tumors presented unique mutations in endometrial cancer-associated genes such as ARID1A and PIK3CA. In conclusion, we demonstrated clonal origin and genomic diversity in a Lynch syndrome-associated endometrial cancer, suggesting the importance of evaluating multiple sites in Lynch syndrome patients with synchronous tumors.

DNA polymerase ε and δ variants drive mutagenesis in polypurine tracts in human tumors

Alterations in the exonuclease domain of DNA polymerase ε cause ultramutated cancers. These cancers accumulate AGA>ATA transversions; however, their genomic features beyond the trinucleotide motifs are obscure. We analyze the extended DNA context of ultramutation using whole-exome sequencing data from 524 endometrial and 395 colorectal tumors. We find that G>T transversions in POLE-mutant tumors predominantly affect sequences containing at least six consecutive purines, with a striking preference for certain positions within polypurine tracts. Using this signature, we develop a machine-learning classifier to identify tumors with hitherto unknown POLE drivers and validate two drivers, POLE-E978G and POLE-S461L, by functional assays in yeast. Unlike other pathogenic variants, the E978G substitution affects the polymerase domain of Pol ε. We further show that tumors with POLD1 drivers share the extended signature of POLE ultramutation. These findings expand the understanding of ultramutation mechanisms and highlight peculiar mutagenic properties of polypurine tracts in the human genome.

Ultra-Deep Sequencing Reveals the Mutational Landscape of Classical Hodgkin Lymphoma

The malignant Hodgkin and Reed Sternberg (HRS) cells of classical Hodgkin lymphoma (cHL) are scarce in affected lymph nodes, creating a challenge to detect driver somatic mutations. As an alternative to cell purification techniques, we hypothesized that ultra-deep exome sequencing would allow genomic study of HRS cells, thereby streamlining analysis and avoiding technical pitfalls. To test this, 31 cHL tumor/normal pairs were exome sequenced to approximately 1,000× median depth of coverage. An orthogonal error-corrected sequencing approach verified >95% of the discovered mutations. We identified mutations in genes novel to cHL including: CDH5 and PCDH7, novel stop gain mutations in IL4R, and a novel pattern of recurrent mutations in pathways regulating Hippo signaling. As a further application of our exome sequencing, we attempted to identify expressed somatic single-nucleotide variants (SNV) in single-nuclei RNA sequencing (snRNA-seq) data generated from a patient in our cohort. Our snRNA analysis identified a clear cluster of cells containing a somatic SNV identified in our deep exome data. This cluster has differentially expressed genes that are consistent with genes known to be dysregulated in HRS cells (e.g., PIM1 and PIM3). The cluster also contains cells with an expanded B-cell clonotype further supporting a malignant phenotype. This study provides proof-of-principle that ultra-deep exome sequencing can be utilized to identify recurrent mutations in HRS cells and demonstrates the feasibility of snRNA-seq in the context of cHL. These studies provide the foundation for the further analysis of genomic variants in large cohorts of patients with cHL.

SIGNIFICANCE

Our data demonstrate the utility of ultra-deep exome sequencing in uncovering somatic variants in Hodgkin lymphoma, creating new opportunities to define the genes that are recurrently mutated in this disease. We also show for the first time the successful application of snRNA-seq in Hodgkin lymphoma and describe the expression profile of a putative cluster of HRS cells in a single patient.

Tumor monocyte content predicts immunochemotherapy outcomes in esophageal adenocarcinoma

For inoperable esophageal adenocarcinoma (EAC), identifying patients likely to benefit from recently approved immunochemotherapy (ICI+CTX) treatments remains a key challenge. We address this using a uniquely designed window-of-opportunity trial (LUD2015-005), in which 35 inoperable EAC patients received first-line immune checkpoint inhibitors for four weeks (ICI-4W), followed by ICI+CTX. Comprehensive biomarker profiling, including generation of a 65,000-cell single-cell RNA-sequencing atlas of esophageal cancer, as well as multi-timepoint transcriptomic profiling of EAC during ICI-4W, reveals a novel T cell inflammation signature (INCITE) whose upregulation correlates with ICI-induced tumor shrinkage. Deconvolution of pre-treatment gastro-esophageal cancer transcriptomes using our single-cell atlas identifies high tumor monocyte content (TMC) as an unexpected ICI+CTX-specific predictor of greater overall survival (OS) in LUD2015-005 patients and of ICI response in prevalent gastric cancer subtypes from independent cohorts. Tumor mutational burden is an additional independent and additive predictor of LUD2015-005 OS. TMC can improve patient selection for emerging ICI+CTX therapies in gastro-esophageal cancer.

Exome-Based Genomic Markers Could Improve Prediction of Checkpoint Inhibitor Efficacy Independently of Tumor Type

Immune checkpoint inhibitors (ICIs) have improved the care of patients in multiple cancer types. However, PD-L1 status, high Tumor Mutational Burden (TMB), and mismatch repair deficiency are the only validated biomarkers of efficacy for ICIs. These markers remain imperfect, and new predictive markers represent an unmet medical need. Whole-exome sequencing was carried out on 154 metastatic or locally advanced cancers from different tumor types treated by immunotherapy. Clinical and genomic features were investigated using Cox regression models to explore their capacity to predict progression-free survival (PFS). The cohort was split into training and validation sets to assess validity of observations. Two predictive models were estimated using clinical and exome-derived variables, respectively. Stage at diagnosis, surgery before immunotherapy, number of lines before immunotherapy, pleuroperitoneal, bone or lung metastasis, and immune-related toxicity were selected to generate a clinical score. KRAS mutations, TMB, TCR clonality, and Shannon entropy were retained to generate an exome-derived score. The addition of the exome-derived score improved the prediction of prognosis compared with the clinical score alone. Exome-derived variables could be used to predict responses to ICI independently of tumor type and might be of value in improving patient selection for ICI therapy.

S. B, Chatterjee M, Prabhash K, Sreekanthreddy P, Rishi KD, Goswami HM, Veldore VH. An absolute approach to using whole exome DNA and RNA workflow for cancer biomarker testing

Introduction

The concept of personalized medicine in cancer has emerged rapidly with the advancement of genome sequencing and the identification of clinically relevant variants that contribute to disease prognosis and facilitates targeted therapy options. In this study, we propose to validate a whole exome-based tumor molecular profiling for DNA and RNA from formalin-fixed paraffin-embedded (FFPE) tumor tissue.

Methods

The study included 166 patients across 17 different cancer types. The scope of this study includes the identification of single-nucleotide variants (SNVs), insertions/deletions (INDELS), copy number alterations (CNAs), gene fusions, tumor mutational burden (TMB), and microsatellite instability (MSI). The assay yielded a mean read depth of 200×, with >80% of on-target reads and a mean uniformity of >90%. Clinical maturation of whole exome sequencing (WES) (DNA and RNA)- based assay was achieved by analytical and clinical validations for all the types of genomic alterations in multiple cancers. We here demonstrate a limit of detection (LOD) of 5% for SNVs and 10% for INDELS with 97.5% specificity, 100% sensitivity, and 100% reproducibility.

Results

The results were >98% concordant with other orthogonal techniques and appeared to be more robust and comprehensive in detecting all the clinically relevant alterations. Our study demonstrates the clinical utility of the exome-based approach of comprehensive genomic profiling (CGP) for cancer patients at diagnosis and disease progression.

Discussion

The assay provides a consolidated picture of tumor heterogeneity and prognostic and predictive biomarkers, thus helping in precision oncology practice. The primary intended use of WES (DNA+RNA) assay would be for patients with rare cancers as well as for patients with unknown primary tumors, and this category constitutes nearly 20–30% of all cancers. The WES approach may also help us understand the clonal evolution during disease progression to precisely plan the treatment in advanced stage disease.

Pelvic Carcinosarcoma Showing a Diverse Histology and Hierarchical Gene Mutation with a Common POLE Mutation to Endometrial Endometroid Carcinoma: A Case Report

POLE mutation-type endometrial cancer is characterized by an extremely high tumor mutation burden. Most POLE mutation-type endometrial cancers are histologically endometrioid carcinomas, and POLE mutation-type carcinosarcomas are rare among endometrial cancers. We report a case of endometrial and pelvic cancer in a 53-year-old woman who was analyzed using next-generating sequencing. The endometrial lesion harbored a p.T457del POLE mutation with an elevated tumor mutation burden and low microsatellite instability. The pelvic lesion showed divergent histological features, consisting of high-grade endometrioid carcinoma, neuroendocrine carcinoma, and chondrosarcoma. In addition to the common POLE mutation detected in the endometrial lesion, the pelvic lesion in each element showed additional gene mutations in a hierarchical manner. Therefore, it is indicated that the p.T457del POLE mutation is a pathogenic mutation and may be related to POLE mutation-induced carcinogenesis and divergent morphogenesis in endometrial cancer.

Molecular Evaluation of Endometrial Dedifferentiated Carcinoma, Endometrioid Carcinoma, Carcinosarcoma, and Serous Carcinoma Using a Custom-Made Small Cancer Panel

It is often difficult to histologically differentiate among endometrial dedifferentiated carcinoma (DC), endometrioid carcinoma (EC), serous carcinoma (SC), and carcinosarcoma (CS) due to the presence of solid components. In this study, we aimed to categorize these carcinomas according to The Cancer Genome Atlas (TCGA) classification using a small custom-made cancer genome panel (56 genes and 17 microsatellite regions) for integrated molecular diagnosis. A total of 36 endometrial cancer cases with solid components were assessed using IHC, next-generation sequencing (NGS), and the custom-made panel. Among 19 EC cases, six were categorized as MMR-deficient (MMR-d) and eight were classified as having a nonspecific molecular profile. Three EC cases were classified as POLE mutation (POLEmut)-type, which had a very high tumor mutation burden (TMB) and low microsatellite instability (MSI). Increased TMB and MSI were observed in all three DC cases, classified as MMR-d with mutations in MLH1 and POLD1. Except for one case classified as MMR-d, all SC cases exhibited TP53 mutations and were classified as p53 mutation-type. SC cases also exhibited amplification of CCND1, CCNE1, and MYC. CS cases were classified as three TCGA types other than the POLEmut-type. The IHC results for p53 and ARID1A were almost consistent with their mutation status. NGS analysis using a small panel enables categorization of endometrial cancers with solid proliferation according to TCGA classification. As TCGA molecular classification does not consider histological findings, an integrated analytical procedure including IHC and NGS may be a practical diagnostic tool for endometrial cancers.

Experimental and Meta-Analytic Validation of RNA Sequencing Signatures for Predicting Status of Microsatellite Instability

Microsatellite instability (MSI) is an important diagnostic and prognostic cancer biomarker. In colorectal, cervical, ovarian, and gastric cancers, it can guide the prescription of chemotherapy and immunotherapy. In laboratory diagnostics of susceptible tumors, MSI is routinely detected by the size of marker polymerase chain reaction products encompassing frequent microsatellite expansion regions. Alternatively, MSI status is screened indirectly by immunohistochemical interrogation of microsatellite binding proteins. RNA sequencing (RNAseq) profiling is an emerging source of data for a wide spectrum of cancer biomarkers. Recently, three RNAseq-based gene signatures were deduced for establishing MSI status in tumor samples. They had 25, 15, and 14 gene products with only one common gene. However, they were developed and tested on the incomplete literature of The Cancer Genome Atlas (TCGA) sampling and never validated experimentally on independent RNAseq samples. In this study, we, for the first time, systematically validated these three RNAseq MSI signatures on the literature colorectal cancer (CRC) (n = 619), endometrial carcinoma (n = 533), gastric cancer (n = 380), uterine carcinosarcoma (n = 55), and esophageal cancer (n = 83) samples and on the set of experimental CRC RNAseq samples (n = 23) for tumors with known MSI status. We found that all three signatures performed well with area under the curve (AUC) ranges of 0.94-1 for the experimental CRCs and 0.94-1 for the TCGA CRC, esophageal cancer, and uterine carcinosarcoma samples. However, for the TCGA endometrial carcinoma and gastric cancer samples, only two signatures were effective with AUC 0.91-0.97, whereas the third signature showed a significantly lower AUC of 0.69-0.88. Software for calculating these MSI signatures using RNAseq data is included.

Performance of Next-Generation Sequencing for the Detection of Microsatellite Instability in Colorectal Cancer With Deficient DNA Mismatch Repair

BACKGROUND & AIMS

Next-generation sequencing (NGS) was recently approved by the United States Food and Drug Administration to detect microsatellite instability (MSI) arising from defective mismatch repair (dMMR) in patients with metastatic colorectal cancer (mCRC) before treatment with immune checkpoint inhibitors (ICI). In this study, we aimed to evaluate and improve the performance of NGS to identify MSI in CRC, especially dMMR mCRC treated with ICI.

METHODS

CRC samples used in this post hoc study were reassessed centrally for MSI and dMMR status using the reference methods of pentaplex polymerase chain reaction and immunohistochemistry. Whole-exome sequencing (WES) was used to evaluate MSISensor, the Food and Drug Administration-approved and NGS-based method for assessment of MSI. This was performed in (1) a prospective, multicenter cohort of 102 patients with mCRC (C1; 25 dMMR/MSI, 24 treated with ICI) from clinical trials NCT02840604 and NCT033501260, (2) an independent retrospective, multicenter cohort of 113 patients (C2; 25 mCRC, 88 non-mCRC, all dMMR/MSI untreated with ICI), and (3) a publicly available series of 118 patients with CRC from The Cancer Genome Atlas (C3; 51 dMMR/MSI). A new NGS-based algorithm, namely MSICare, was developed. Its performance for assessment of MSI was compared with MSISensor in C1, C2, and C3 at the exome level or after downsampling sequencing data to the MSK-IMPACT gene panel. MSICare was validated in an additional retrospective, multicenter cohort (C4) of 152 patients with new CRC (137 dMMR/MSI) enriched in tumors deficient in MSH6 (n = 35) and PMS2 (n = 9) after targeted sequencing of samples with an optimized set of microsatellite markers (MSIDIAG).

RESULTS

At the exome level, MSISensor was highly specific but failed to diagnose MSI in 16% of MSI/dMMR mCRC from C1 (4 of 25; sensitivity, 84%; 95% confidence interval [CI], 63.9%-95.5%), 32% of mCRC (8 of 25; sensitivity, 68%; 95% CI, 46.5%-85.1%), and 9.1% of non-mCRC from C2 (8 of 88; sensitivity, 90.9%; 95% CI, 82.9%-96%), and 9.8% of CRC from C3 (5 of 51; sensitivity, 90.2%; 95% CI, 78.6%-96.7%). Misdiagnosis included 4 mCRCs treated with ICI, of which 3 showed an overall response rate without progression at this date. At the exome level, reevaluation of the MSI genomic signal using MSICare detected 100% of cases with true MSI status among C1 and C2. Further validation of MSICare was obtained in CRC tumors from C3, with 96.1% concordance for MSI status. Whereas misdiagnosis with MSISensor even increased when analyzing downsampled WES data from C1 and C2 with microsatellite markers restricted to the MSK-IMPACT gene panel (sensitivity, 72.5%; 95% CI, 64.2%-79.7%), particularly in the MSH6-deficient setting, MSICare sensitivity and specificity remained optimal (100%). Similar results were obtained with MSICare after targeted NGS of tumors from C4 with the optimized microsatellite panel MSIDIAG (sensitivity, 99.3%; 95% CI, 96%-100%; specificity, 100%).

CONCLUSIONS

In contrast to MSISensor, the new MSICare test we propose performs at least as efficiently as the reference method, MSI polymerase chain reaction, to detect MSI in CRC regardless of the defective MMR protein under both WES and targeted NGS conditions. We suggest MSICare may rapidly become a reference method for NGS-based testing of MSI in CRC, especially in mCRC, where accurate MSI status is required before the prescription of ICI.

Construction of a novel methylation-related prognostic model for colorectal cancer based on microsatellite status

The present study aimed to construct a novel methylation-related prognostic model based on microsatellite status that may enhance the prognosis of colorectal cancer (CRC) from methylation and microsatellite status perspective. DNA methylation and mRNA expression data with clinical information were downloaded from The Cancer Genome Atlas (TCGA) data set. The samples were divided into microsatellite stability and microsatellite instability group, and CIBERSORT was used to assess the immune cell infiltration characteristics. After identifying the differentially methylated genes and differentially expression genes using R packages, the methylation-driven genes were further identified. Prognostic genes that were used to establish the methylation-related risk score model were generated by the univariate and multivariate Cox regression model. Finally, we established and evaluated the methylation-related prognostic model for CRC patients. A total of 69 MDGs were obtained and three of these genes (MIOX, TH, DKFZP434K028) were selected to construct the prognostic model. Patients in the low-risk score group had a conspicuously better overall survival than those in the high-risk score group (p < .0001). The area under the receiver operating characteristic curve for this model was 0.689 at 3 years, 0.674 at 4 years, and 0.658 at 5 years. The Wilcoxon test showed that higher risk score was associated with higher T stage (p = .01), N stages (p = .0028), metastasis (p = .013), and advanced pathological stage (p = .0013). However, the more instability of microsatellite status, the lower risk score of CRC patients (p = .0048). Our constructed methylation-related prognostic model based on microsatellite status presents potential significance in assessing recurrence risk stratification, tumor staging, and immunotherapy for CRC patients.

Direct next-generation sequencing analysis using endometrial liquid-based cytology specimens for rapid cancer genomic profiling

BACKGROUND

Genomic examination of cytology specimens is often performed on cell blocks or conventional smears rather than on liquid-based cytology (LBC) specimens. Since LBC specimens preserve high-quality DNA, cancer genome profiling using next-generation sequencing (NGS) is also attainable from residual LBC specimens. One of the advantages of using LBC specimens for NGS is that it allows direct extraction of DNA from residual specimens, avoiding a sacrifice of smear slides and minimizing genomic profiling processing time.

METHODS

Endometrial LBC specimens were subjected to NGS analysis to validate the practicality of rapid cancer genomic profiling in a pathology laboratory. The extracted DNA was subjected to NGS using a customized cancer gene panel comprising 56 genes and 17 microsatellite regions. The workflow strategy was defined, and the processing time estimated for specimen sampling, cell counting, NGS run, and genome profiling.

RESULTS

NGS analysis of most LBC specimens revealed somatic mutations, tumor mutation burden, and microsatellite instability, which were almost identical to those obtained from formalin-fixed paraffin-embedded tissues. The processing time for direct NGS analysis and cancer genomic profiling of the residual LBC specimens was approximately 5 days.

CONCLUSION

The residual LBC specimens collected using endometrial cytology were verified to carry a high tumor fraction for NGS analysis and could serve as an alternate source for rapid molecular classification and diagnosis of endometrial cancers, as a routine process in a pathology laboratory.

FDA-Approved and Emerging Next Generation Predictive Biomarkers for Immune Checkpoint Inhibitors in Cancer Patients

A patient's response to immune checkpoint inhibitors (ICIs) is a complex quantitative trait, and determined by multiple intrinsic and extrinsic factors. Three currently FDA-approved predictive biomarkers (progra1mmed cell death ligand-1 (PD-L1); microsatellite instability (MSI); tumor mutational burden (TMB)) are routinely used for patient selection for ICI response in clinical practice. Although clinical utility of these biomarkers has been demonstrated in ample clinical trials, many variables involved in using these biomarkers have poised serious challenges in daily practice. Furthermore, the predicted responders by these three biomarkers only have a small percentage of overlap, suggesting that each biomarker captures different contributing factors to ICI response. Optimized use of currently FDA-approved biomarkers and development of a new generation of predictive biomarkers are urgently needed. In this review, we will first discuss three widely used FDA-approved predictive biomarkers and their optimal use. Secondly, we will review four novel gene signature biomarkers: T-cell inflamed gene expression profile (GEP), T-cell dysfunction and exclusion gene signature (TIDE), melanocytic plasticity signature (MPS) and B-cell focused gene signature. The GEP and TIDE have shown better predictive performance than PD-L1, and PD-L1 or TMB, respectively. The MPS is superior to PD-L1, TMB, and TIDE. The B-cell focused gene signature represents a previously unexplored predictive biomarker to ICI response. Thirdly, we will highlight two combined predictive biomarkers: TMB+GEP and MPS+TIDE. These integrated biomarkers showed improved predictive outcomes compared to a single predictor. Finally, we will present a potential nucleic acid biomarker signature, allowing DNA and RNA biomarkers to be analyzed in one assay. This comprehensive signature could represent a future direction of developing robust predictive biomarkers, particularly for the cold tumors, for ICI response.

Hypermutated phenotype in gliosarcoma of the spinal cord

Gliosarcoma is a variant of glioblastoma with equally poor prognosis and characterized by mixed glial and mesenchymal pathology. Metastasis is not uncommon but the involvement of the spinal cord is rare, and comprehensive genetic characterization of spinal gliosarcoma is lacking. We describe a patient initially diagnosed with a low-grade brain glioma via biopsy, followed by adjuvant radiation and temozolomide treatment. Nearly 2 years after diagnosis, she developed neurological deficits from an intradural, extramedullary tumor anterior to the spinal cord at T4, which was resected and diagnosed as gliosarcoma. Whole-exome sequencing (WES) of this tumor revealed a hypermutated phenotype, characterized by somatic mutations in key DNA mismatch repair (MMR) pathway genes, an abundance of C>T transitions within the identified somatic single nucleotide variations, and microsatellite stability, together consistent with temozolomide-mediated hypermutagenesis. This is the first report of a hypermutator phenotype in gliosarcoma, which may represent a novel genomic mechanism of progression from lower grade glioma.

Next-generation sequencing analysis of endometrial screening liquid-based cytology specimens: a comparative study to tissue specimens

Background

Liquid-based cytology (LBC) is now a widely used method for cytologic screening and cancer diagnosis. Since the cells are fixed with alcohol-based fixatives, and the specimens are stored in a liquid condition, LBC specimens are suitable for genetic analyses.

Methods

Here, we established a small cancer gene panel, including 60 genes and 17 microsatellite markers for next-generation sequencing, and applied to residual LBC specimens obtained by endometrial cancer screening to compare with corresponding formalin-fixed paraffin-embedded (FFPE) tissues.

Results

A total of 49 FFPE and LBC specimens (n = 24) were analyzed, revealing characteristic mutations for endometrial cancer, including PTEN, CTNNB1, PIK3CA, and PIK3R1 mutations. Eight cases had higher scores for both tumor mutation burden (TMB) and microsatellite instability (MSI), which agree with defective mismatch repair (MMR) protein expression. Paired endometrial LBC, and biopsied and/or resected FFPE tissues from 7 cases, presented almost identical mutations, TMB, and MSI profiles in all cases.

Conclusion

These findings demonstrate that our ad hoc cancer gene panel enabled the detection of therapeutically actionable gene mutations in endometrial LBC and FFPE specimens. Endometrial cancer LBC specimens offer an alternative and affordable source of molecular testing materials.