Subclonal evolution of pulmonary adenocarcinomas delineated by spatially distributed somatic mitochondrial mutations

Comparison of Two Variant Analysis Programs for Next-Generation Sequencing Data of Whole Mitochondrial Genome

BACKGROUND

With advance of next-generation sequencing (NGS) techniques, the need for mitochondrial DNA analysis is increasing not only in the forensic area, but also in medical fields.

METHODS

Two commercial programs, Converge Software (CS) and Torrent Variant Caller for variant calling of NGS data, were compared with a considerable amount of sequence data of 50 samples with a homogeneous ethnicity.

RESULTS

About 2,300 variants were identified and the two programs showed about 90% of consistency. CS, a dedicated analysis program for mitochondrial DNA, showed some advantages for forensic use. By additional visual inspection, several causes of discrepancy in variant calling results were identified. Application of different notation rules for mitochondrial sequence and the minor allele frequency close to detection threshold were the two most significant reasons.

CONCLUSION

With prospective improvement of each program, researchers and practitioners should be aware of characteristics of the analysis program they use and prepare their own strategies to determine variants.

Metastatic pattern of ovarian cancer delineated by tracing the evolution of mitochondrial DNA mutations

Ovarian cancer (OC) is the most lethal gynecologic tumor and is characterized by a high rate of metastasis. Challenges in accurately delineating the metastatic pattern have greatly restricted the improvement of treatment in OC patients. An increasing number of studies have leveraged mitochondrial DNA (mtDNA) mutations as efficient lineage-tracing markers of tumor clonality. We applied multiregional sampling and high-depth mtDNA sequencing to determine the metastatic patterns in advanced-stage OC patients. Somatic mtDNA mutations were profiled from a total of 195 primary and 200 metastatic tumor tissue samples from 35 OC patients. Our results revealed remarkable sample-level and patient-level heterogeneity. In addition, distinct mtDNA mutational patterns were observed between primary and metastatic OC tissues. Further analysis identified the different mutational spectra between shared and private mutations among primary and metastatic OC tissues. Analysis of the clonality index calculated based on mtDNA mutations supported a monoclonal tumor origin in 14 of 16 patients with bilateral ovarian cancers. Notably, mtDNA-based spatial phylogenetic analysis revealed distinct patterns of OC metastasis, in which a linear metastatic pattern exhibited a low degree of mtDNA mutation heterogeneity and a short evolutionary distance, whereas a parallel metastatic pattern showed the opposite trend. Moreover, a mtDNA-based tumor evolutionary score (MTEs) related to different metastatic patterns was defined. Our data showed that patients with different MTESs responded differently to combined debulking surgery and chemotherapy. Finally, we observed that tumor-derived mtDNA mutations were more likely to be detected in ascitic fluid than in plasma samples. Our study presents an explicit view of the OC metastatic pattern, which sheds light on efficient treatment for OC patients.

Spatial profiling of the microenvironment reveals low intratumoral heterogeneity and STK11-associated immune evasion in therapy-naïve lung adenocarcinomas

OBJECTIVE

Intratumoral heterogeneity was found to be a significant factor causing resistance to lung cancer therapies, including immune checkpoint blockade. Lesser is known about spatial heterogeneity of the tumor microenvironment (TME) and its association with genetic properties of the tumor, which is of particular interest in the therapy-naïve setting.

MATERIALS AND METHODS

We performed multi-region sampling (2-4 samples per tumor; total of 55 samples) from a cohort of 19 untreated stage IA-IIIB lung adenocarcinomas (n = 11 KRAS mutant, n = 1 ERBB2 mutant, n = 7 KRAS wildtype). For each sample the expression of 770 immunooncology-related genes was analyzed using the nCounter platform, while the mutational status was determined by hybrid capture-based next-generation sequencing (NGS) using a large panel covering more than 500 genes.

RESULTS

Global unsupervised analyses revealed clustering of the samples into two groups corresponding to a 'hot' or 'cold' immunologic tumor contexture based on the abundance of immune cell infiltrates. All analyzed specific immune cell signatures (ICsig) showed a significantly higher intertumoral than intratumoral heterogeneity (p < 0.02), as most of the analyzed cases (14/19) showed a very homogenous spatial immune cell profile. PD-L1 exhibited a significantly higher intertumoral than intratumoral heterogeneity (p = 1.03e-13). We found a specific association with 'cold' TME for STK11 (11/14, p < 0.07), but not KRAS, TP53, LRP1B, MTOR, U2AF1 co-mutations, and validated this finding using The Cancer Genome Atlas (TCGA) data.

CONCLUSION

Early-stage lung adenocarcinomas show considerable intertumoral, but limited intratumoral heterogeneity, which is clinically highly relevant as assessment before neoadjuvant treatment is based on small biopsies. STK11 mutations are specifically associated with a 'cold' TME, which could affect the efficacy of perioperative immunotherapy.

Mitochondrial DNA Mutation Analysis in Breast Cancer: Shifting From Germline Heteroplasmy Toward Homoplasmy in Tumors

Studies have suggested a potential role of somatic mitochondrial mutations in cancer development. To analyze the landscape of somatic mitochondrial mutation in breast cancer and to determine whether mitochondrial DNA (mtDNA) mutational burden is correlated with overall survival (OS), we sequenced whole mtDNA from 92 matched-paired primary breast tumors and peripheral blood. A total of 324 germline variants and 173 somatic mutations were found in the tumors. The most common germline allele was 663G (12S), showing lower heteroplasmy levels in peripheral blood lymphocytes than in their matched tumors, even reaching homoplasmic status in several cases. The heteroplasmy load was higher in tumors than in their paired normal tissues. Somatic mtDNA mutations were found in 73.9% of breast tumors; 59% of these mutations were located in the coding region (66.7% non-synonymous and 33.3% synonymous). Although the CO1 gene presented the highest number of mutations, tRNA genes (T,C, and W), rRNA 12S, and CO1 and ATP6 exhibited the highest mutation rates. No specific mtDNA mutational profile was associated with molecular subtypes of breast cancer, and we found no correlation between mtDNA mutational burden and OS. Future investigations will provide insight into the molecular mechanisms through which mtDNA mutations and heteroplasmy shifting contribute to breast cancer development.

Approaching sites of action of drugs in clinical pharmacology: New analytical options and their challenges

Clinical pharmacology is an important discipline for drug development aiming to define pharmacokinetics (PK), pharmacodynamics (PD) and optimum exposure to drugs, i.e. the concentration-response relationship and its modulators. For this purpose, information on drug concentrations at the anatomical, cellular and molecular sites of action is particularly valuable. In pharmacological assays, the limited accessibility of target cells in readily available samples (i.e. blood) often hampers mass spectrometry-based monitoring of the absolute quantity of a compound and the determination of its molecular action at the cellular level. Recently, new sample collection methods have been developed for the specific capture of rare circulating cells, especially for the diagnosis of circulating tumour cells. In parallel, new advances and developments in mass spectrometric instrumentation now allow analyses to be scaled down to the cellular level. Together, these developments may permit the monitoring of minute drug quantities and show their effect at the cellular level. In turn, such PK/PD associations on a cellular level would not only enrich our pharmacological knowledge of a given compound but also expand the basis for PK/PD simulations. In this review, we describe novel concepts supporting clinical pharmacology at the anatomical, cellular and molecular sites of action, and highlight the new challenges in mass spectrometry-based monitoring. Moreover, we present methods to tackle these challenges and define future needs.

Somatic mitochondrial mutation discovery using ultra-deep sequencing of the mitochondrial genome reveals spatial tumor heterogeneity in head and neck squamous cell carcinoma

Mutations in mitochondrial DNA (mtDNA) have been linked to risk, progression, and treatment response of head and neck squamous cell carcinoma (HNSCC). Due to their clonal nature and high copy number, mitochondrial mutations could serve as powerful molecular markers for detection of cancer cells in bodily fluids, surgical margins, biopsies and lymph node (LN) metastasis, especially at sites where tumor involvement is not histologically apparent. Despite a pressing need for high-throughput, cost-effective mtDNA mutation profiling system, current methods for library preparation are still imperfect for detection of low prevalence heteroplasmic mutations. To this end, we have designed an ultra-deep amplicon-based sequencing library preparation approach that covers the entire mitochondrial genome. We sequenced mtDNA in 28 HNSCCs, matched LNs, surgical margins and bodily fluids, and applied multiregional sequencing approach on 14 primary tumors. Our results demonstrate that this quick, sensitive and cost-efficient method allows obtaining a snapshot on the mitochondrial heterogeneity, and can be used for detection of low frequency tumor-associated mtDNA mutations in LNs, sputum and serum specimens. These findings provide the foundation for using mitochondrial sequencing for risk assessment, early detection, and tumor surveillance.

Automated sample preparation with SP3 for low-input clinical proteomics

High-throughput and streamlined workflows are essential in clinical proteomics for standardized processing of samples from a variety of sources, including fresh-frozen tissue, FFPE tissue, or blood. To reach this goal, we have implemented single-pot solid-phase-enhanced sample preparation (SP3) on a liquid handling robot for automated processing (autoSP3) of tissue lysates in a 96-well format. AutoSP3 performs unbiased protein purification and digestion, and delivers peptides that can be directly analyzed by LCMS, thereby significantly reducing hands-on time, reducing variability in protein quantification, and improving longitudinal reproducibility. We demonstrate the distinguishing ability of autoSP3 to process low-input samples, reproducibly quantifying 500-1,000 proteins from 100 to 1,000 cells. Furthermore, we applied this approach to a cohort of clinical FFPE pulmonary adenocarcinoma (ADC) samples and recapitulated their separation into known histological growth patterns. Finally, we integrated autoSP3 with AFA ultrasonication for the automated end-to-end sample preparation and LCMS analysis of 96 intact tissue samples. Collectively, this constitutes a generic, scalable, and cost-effective workflow with minimal manual intervention, enabling reproducible tissue proteomics in a broad range of clinical and non-clinical applications.

Spatial and Temporal Heterogeneity of Panel-Based Tumor Mutational Burden in Pulmonary Adenocarcinoma: Separating Biology From Technical Artifacts

BACKGROUND

Tumor mutational burden (TMB) is an emerging biomarker used to identify patients who are more likely to benefit from immuno-oncology therapy. Aside from various unsettled technical aspects, biological variables such as tumor cell content and intratumor heterogeneity may play an important role in determining TMB.

METHODS

TMB estimates were determined applying the TruSight Oncology 500 targeted sequencing panel. Spatial and temporal heterogeneity was analyzed by multiregion sequencing (two to six samples) of 24 pulmonary adenocarcinomas and by sequencing a set of matched primary tumors, locoregional lymph node metastases, and distant metastases in five patients.

RESULTS

On average, a coding region of 1.28 Mbp was covered with a mean read depth of 609x. Manual validation of the mutation-calls confirmed a good performance, but revealed noticeable misclassification during germline filtering. Different regions within a tumor showed considerable spatial TMB variance in 30% (7 of 24) of the cases (maximum difference, 14.13 mut/Mbp). Lymph node-derived TMB was significantly lower (p = 0.016). In 13 cases, distinct mutational profiles were exclusive to different regions of a tumor, leading to higher values for simulated aggregated TMB. Combined, intratumor heterogeneity and the aggregated TMB could result in divergent TMB designation in 17% of the analyzed patients. TMB variation between primary tumor and distant metastases existed but was not profound.

CONCLUSIONS

Our data show that, in addition to technical aspects such as germline filtering, the tumor content and spatially divergent mutational profiles within a tumor are relevant factors influencing TMB estimation, revealing limitations of single-sample-based TMB estimations in a clinical context.