Evolution analysis of heterogeneous non-small cell lung carcinoma by ultra-deep sequencing of the mitochondrial genome

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.

Mito-SiPE is a sequence-independent and PCR-free mtDNA enrichment method for accurate ultra-deep mitochondrial sequencing

The analysis of somatic variation in the mitochondrial genome requires deep sequencing of mitochondrial DNA. This is ordinarily achieved by selective enrichment methods, such as PCR amplification or probe hybridization. These methods can introduce bias and are prone to contamination by nuclear-mitochondrial sequences (NUMTs), elements that can introduce artefacts into heteroplasmy analysis. We isolated intact mitochondria using differential centrifugation and alkaline lysis and subjected purified mitochondrial DNA to a sequence-independent and PCR-free method to obtain ultra-deep (>80,000X) sequencing coverage of the mitochondrial genome. This methodology avoids false-heteroplasmy calls that occur when long-range PCR amplification is used for mitochondrial DNA enrichment. Previously published methods employing mitochondrial DNA purification did not measure mitochondrial DNA enrichment or utilise high coverage short-read sequencing. Here, we describe a protocol that yields mitochondrial DNA and have quantified the increased level of mitochondrial DNA post-enrichment in 7 different mouse tissues. This method will enable researchers to identify changes in low frequency heteroplasmy without introducing PCR biases or NUMT contamination that are incorrectly identified as heteroplasmy when long-range PCR is used.

Epigenetic modifications precede molecular alterations and drive human hepatocarcinogenesis

Development of primary liver cancer is a multi-stage process. Detailed understanding of sequential epigenetic alterations is largely missing. Here, we performed Infinium Human Methylation 450k BeadChips and RNA sequencing analyses for genome-wide methylome and transcriptome profiling of cirrhotic liver (n=7), low- (n=4) and high-grade (n=9) dysplastic lesions, early (n=5) and progressed (n=3) hepatocellular carcinomas (HCC) synchronously detected in eight HCC patients with chronic hepatitis B infection. Integrative analyses of epigenetically driven molecular changes were identified and validated in two independent cohorts comprising 887 HCC. Mitochondrial DNA sequencing was further employed for clonality analyses and indicates multi-clonal origin in the majority of investigated HCC. Alterations in DNA methylation progressively increased from CL to dysplastic lesions and reached a maximum in early HCC. Associated early alterations identified by IPA pathway analyses involved apoptosis, immune regulation and stemness pathways, while late changes centered on cell survival, proliferation and invasion. We further validated putative 23 epi-drivers with concomitant expression changes and associated with overall survival. Functionally, Striatin 4 (STRN4) was demonstrated to be epigenetically regulated and inhibition of STRN4 significantly suppressed tumorigenicity of HCC cell lines.Overall, application of integrative genomic analyses defines epigenetic driver alterations and provides promising targets for novel therapeutic approaches.

Epigenetic Regulation of Mitochondrial Quality Control Genes in Multiple Myeloma: A Sequenom MassARRAY Pilot Investigation on HMCLs

The mitochondrial quality control network includes several epigenetically-regulated genes involved in mitochondrial dynamics, mitophagy, and mitochondrial biogenesis under physiologic conditions. Dysregulated expression of such genes has been reported in various disease contexts, including cancer. However, their expression pattern and the possible underlying epigenetic modifications remain to be defined within plasma cell (PC) dyscrasias. Herein, we compared the mRNA expression of mitochondrial quality control genes from multiple myeloma, plasma cell leukemia patients and human myeloma cell lines (HMCLs) with healthy plasma cells; moreover, by applying the Sequenom MassARRAY EpiTYPER technology, we performed a pilot investigation of their CpG methylation status in HMCLs. Overall, the results provided indicate dysregulated expression of several mitochondrial network’s genes, and alteration of the CpG methylation profile, underscoring novel potential myeloma biomarkers deserving in-depth functional investigation in the future.

Expression profiling on subclasses of primary parotid gland carcinomas

INTRODUCTION

The underlying molecular mechanisms of parotid gland carcinomas (PGC) are still unknown. Knowledge about the tumor-driving signaling pathways is necessary either for diagnostics or developing new therapeutic options in this heterogeneous and rare entity.

MATERIAL AND METHODS

94 matching RNA formalin-fixed and paraffin-embedded tissue samples from PGC and the corresponding non-tumor area, RNA quality and quantity were sufficient for gene expression profiling of 770 genes using the NanoString's nCounter technology. Oncogenic and tumor suppressor genes were examined in the three common PGC tumor entities: adenoid cystic carcinoma (ACC), adenocarcinoma NOS (AC-NOS), and mucoepidermoid carcinoma (MEC).

RESULTS

Expression profiling and subsequent hierarchical cluster analysis clearly differentiated between non-tumor gland tissue samples and PGC. In addition expression pattern of all three entities differed. The extensive pathway analysis proved a prominent dysregulation of the Wnt signaling pathway in the three PGC entities. Moreover, transcript upstream analysis demonstrated a pronounced activation of the PI3K pathway in ACC and MEC.

DISCUSSION

Our findings revealed divergent molecular expression profiles in MEC, ACC and AC-NOS that are presently studied for their potential application in PGC diagnostics. Importantly, identification of Wnt and PI3K signaling in PGC revealed novel options of PGC therapy.

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.

Optimized PCR-Based Enrichment Improves Coverage Uniformity and Mutation Detection in Mitochondrial DNA Next-Generation Sequencing

Next-generation sequencing-based methods have been commonly used for detecting mutations of mitochondrial genome (mtDNA). PCR amplification is a highly effective method of mtDNA enrichment before sequencing. However, it has been observed that highly variable sequencing depth within PCR amplicons severely reduces the coverage uniformity and accuracy of mutation calling. Therefore, it is necessary to develop an optimized PCR-based strategy for mtDNA sequencing. Herein, the effect of DNA quality on the efficiency of PCR amplification was analyzed and the effects of different primer-design methods, including the number of primer pairs, overlap length of amplicons, and modification of primers, on coverage uniformity and mutation calling in mtDNA sequencing were assessed. Results showed that DNA quality significantly affected the efficiency of PCR amplification. Importantly, over- and under-representation of coverage depth at overlap regions of amplicons were observed when amplicons were not modified and overlap was shorter than two sequencing fragment sizes (800 bp). Then, under-representation was overcome by increasing the overlap length of the amplicons, and over-representation was effectively reduced by 5'-block modification of primers and sticky-end ligation of amplicons. Moreover, findings indicated that these two optimized PCR-based sequencing strategies effectively improved mutation calling in primer-binding regions. Optimized PCR-based mtDNA enrichment and sequencing approaches have been established, which laid a foundation for accurate mutation detection of mtDNA in diseases.

Mitochondrion: I am more than a fuel server

Apart from reliable management of the "powerhouse" of the cell, mitochondria faithfully orchestrate a diverse array of important and critical functions in governing cellular signaling, apoptosis, autophagy, mitophagy and innate and adaptive immune system. Introduction of instability and imbalance in the mitochondrial own genome or the nuclear encoded mitochondrial proteome would result in the manifestation of various diseases through alterations in the oxidative phosphorylation system (OXPHOS) and nuclear-mitochondria retrograde signaling. Understanding mitochondrial biology and dynamism are thus of paramount importance to develop strategies to prevent or treat various diseases caused due to mitochondrial alterations.

The Optimal Sequencing Depth of Tumor Biopsies for Identifying Clonal Cell Populations

The tumor content of a biopsy and the average depth of coverage are two essential aspects when performing DNA sequencing using next-generation sequencing technologies. The heterogeneous nature of cancer necessitates the identification of distinct clonal cell populations to better understand and treat cancer. Deep sequencing enables researchers to identify these populations, but no consensus on an optimal depth exists for identifying clonal populations. Data from eight deep-sequenced oral squamous cell carcinoma biopsies obtained from three stage IV patients, with various degrees of tumor content, were used to randomly down sample the depth before being subjected to cluster analysis. An increase in coverage resulted in an increase in resolution for clusters of mutations, enabling the identification of distinct clonal cell populations and clonal events. From a depth of 800×, limited gain in resolution can be achieved; and from a depth of 1200×, the resolution stabilizes. Overall, researchers should aim for an average depth of 1000× to 1200× when performing deep sequencing. The tumor content will, however, dictate the resolution and fidelity of the analysis, as an increase in tumor complexity increases the need for higher tumor content.

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

OBJECTIVES

The potential role of cancer associated somatic mutations of the mitochondrial genome (mtDNA) is controversial and still poorly understood. Our group and others recently challenged a direct tumorigenic impact and suggested a passenger-like character. In combination with the known increased mutation rate, somatic mtDNA mutations account for an interesting tool to delineate tumor evolution. Here, we comprehensively analyzed the spatial distribution of somatic mtDNA mutations throughout whole tumor sections of pulmonary adenocarcinoma (ADC).

MATERIALS AND METHODS

Central sections of 19 ADC were analyzed in a segmented manner (11-34 segments/tumor) together with non-neoplastic tissue samples and lymph node metastasis, if present. We performed whole mtDNA sequencing and real-time PCR based quantification of mtDNA copy numbers for all samples. Further, histological growth patterns were determined on H&E sections and the tumor cell content was quantified by digital pathology analyses.

RESULTS

Somatic mtDNA mutations were present in 96% (18/19) of the analyzed tumors, either ubiquitously or restricted to specific tumor regions. Spatial and histological mapping of the mutations enabled the identification of subclonal structures and phylogenetic relations within a tumor section indicating different progression levels. In this regard, lymph node metastases seem to be related to early events in ADC development. There was no concurrence between histological and mtDNA mutation based clusters. However, micropapillary patterns occurred only in tumors with ubiquitous mutations. ADC with more than two ubiquitous mutations were associated with shorter disease-free survival (p < 0.01).

CONCLUSION

Cancer related mtDNA mutations are interesting candidates for the understanding of subclonal ADC evolution and perspectively for monitoring tumor progression. Our data reveal a potential prognostic relevance of somatic mtDNA mutations.