mtDNA Next-generation Sequencing Research Articles

A Customized Human Mitochondrial DNA Database (hMITO DB v1.0) for Rapid Sequence Analysis, Haplotyping and Geo-Mapping.

The field of mitochondrial genomics has advanced rapidly and has revolutionized disciplines such as molecular anthropology, population genetics, and medical genetics/oncogenetics. However, mtDNA next-generation sequencing (NGS) analysis for matrilineal haplotyping and phylogeographic inference remains hindered by the lack of a consolidated mitogenome database and an efficient bioinformatics pipeline. To address this, we developed a customized human mitogenome database (hMITO DB) embedded in a CLC Genomics workflow for read mapping, variant analysis, haplotyping, and geo-mapping. The database was constructed from 4286 mitogenomes. The macro-haplogroup (A to Z) distribution and representative phylogenetic tree were found to be consistent with published literature. The hMITO DB automated workflow was tested using mtDNA-NGS sequences derived from Pap smears and cervical cancer cell lines. The auto-generated read mapping, variants track, and table of haplotypes and geo-origins were completed in 15 min for 47 samples. The mtDNA workflow proved to be a rapid, efficient, and accurate means of sequence analysis for translational mitogenomics.

MitoSomatic: a tool for accurate identification of mitochondrial DNA somatic mutations without paired controls.

Mitochondrial DNA (mtDNA) somatic mutations play important roles in the initiation and progression of cancer. Although next-generation sequencing (NGS) of paired tumor and control samples has become a common practice to identify tumor-specific mtDNA mutations, the unique nature of mtDNA and NGS-associated sequencing bias could cause false-positive/-negative somatic mutation calling. Additionally, there are clinical scenarios where matched control tissues are unavailable for comparison. Therefore, a novel approach for accurately identifying somatic mtDNA variants is greatly needed, particularly in the absence of matched controls. In this study, the ground truth mtDNA variants orthogonally validated by triple-paired tumor, adjacent nontumor, and blood samples were used to develop mitoSomatic, a random forest-based machine learning tool. We demonstrated that mitoSomatic achieved area under the curve (AUC) values over 0.99 for identifying somatic mtDNA variants without paired control in three tumor types. In addition, mitoSomatic was also applicable in nontumor tissues such as adjacent nontumor and blood samples, suggesting the flexibility of mitoSomatic's classification capability. Furthermore, analysis of triple-paired samples identified a small group of variants with uncertain somatic/germline origin, whereas application of mitoSomatic significantly facilitated the prediction of their possible source. Finally, a control-free evaluation of the public pan-cancer NGS dataset with mitoSomatic revealed a substantial number of variants that were probably misclassified by conventional tumor-control comparison, further emphasizing the usefulness of mitoSomatic in application. Taken together, our study demonstrates that mitoSomatic is valuable for accurately identifying somatic mtDNA variants in mtDNA NGS data without paired controls, applicable for both tumor and nontumor tissues.

Assessment of the Role of NuclearENDOGGene and mtDNA Variations on Paternal Mitochondrial Elimination (PME) in Infertile Men: An Experimental Study.

In humans and most animals, maternal inheritance of mitochondria and mitochondrial DNA (mtDNA) is considered as an universal assumption. Recently, several lines of evidence suggest that different species seem to employ distinct mechanisms to prevent the inheritance of paternal mtDNA. There are few studies in the literature on the molecular basis of sperm mtDNA elimination in mammals and paternal mtDNA transmission in humans. Endonuclease G (ENDOG) is a mitochondrial nuclease encoded by nuclear ENDOG gene. The critical importance of ENDOG gene on paternal mitochondrial elimination (PME) has been previously demonstrated in model organisms such as C. elegans and D. melanogaster. However, its mechanism in human is still unclear. Therefore, we aimed to evaluate whether nuclear ENDOG gene copy number could be a potential marker of paternal mtDNA transmission or not.Male factor infertility patients diagnosed with different infertility subgroups such as azoospermia, oligoteratozoospermia, astheno-teratozoospermia were included in this study: 13 infertile men and 25 healthy men as control group. Quantitative real-time polymerase chain reaction (qPCR) analysis and dual-color Fluorescence in situ hybridization (FISH) method were used to compare the groups. FISH method was applied to verify qPCR results and two signals were observed in nearly all patients. ENDOG gene copy number data were evaluated by comparing them with entire human mtDNA next-generation sequencing (NGS) analysis results obtained through bioinformatics and proteomics tools. Mitochondrial whole genome sequencing (WGS) data allowed determination of novel and reported variations such as single nucleotide polymorphisms (SNPs), multiple nucleotide polymorphism (MNP), insertion/deletion (INDEL). Missense variants causing amino acid substitution were filtered out from patients' mtDNA WGS data.Relative copy number of target ENDOG gene in male infertility patients [0.49 (0.31 - 0.77)] was lower than healthy controls [1.00 (0.66 - 1.51)], and statistical results showed significant differences between the groups (p < 0.01). A total of 38 missense variants were detected in the genes encoding the proteins involved in the respiratory chain complex. Moreover, we detected paternal mtDNA transmissions in the children of these patients who applied to assisted reproductive techniques.In conclusion, this study reveals that ENDOG gene may be an important factor for the PME mechanism in humans. To the best of our knowledge, this is the first study in humans about this topic and assessment of ENDOG gene sequencing and gene expression studies in a larger sample size including patients with male factor infertility would be our future project.

MitoDataclean: A machine learning approach for the accurate identification of cross-contamination-derived tumor mitochondrial DNA mutations.

Next-generation sequencing (NGS) of mitochondrial DNA (mtDNA) has widespread applications in aging and cancer studies. However, cross-contamination of mtDNA constitutes a major concern. Previous methods for the detection of mtDNA contamination mainly focus on haplogroup-level phylogeny, but neglect haplotype-level differences, leading to limited sensitivity and accuracy. In our study, we present mitoDataclean, a random-forest-based machine learning package for accurate identification of cross-contamination, evaluation of contamination levels and detection of contamination-derived variants in mtDNA NGS data. Comprehensive optimization of mitoDataclean revealed that training simulation with mixtures of small haplogroup distance and low polymorphic difference was critical for optimal modeling. Compared to existing methods, mitoDataclean exhibited significantly improved sensitivity and accuracy for the detection of sample contamination in simulated data. In addition, mitoDataclean achieved area under the curve values of 0.91 and 0.97 for discerning genuine and contamination-derived mtDNA variants in a simulated Western dataset and private sequencing contamination data, respectively, suggesting that this tool may be applicable for different populations and samples with different sources of contamination. Finally, mitoDataclean was further evaluated in several private and public datasets and showed a robust ability for contamination detection. Altogether, our study demonstrates that mitoDataclean may be used for accurate detection of contaminated samples and contamination-derived variants in mtDNA NGS data.

Identification of Variants in Mitochondrial D-Loop and OriL Region and Analysis of Mitochondrial DNA Copy Number in Women with Polycystic Ovary Syndrome.

Polycystic ovary syndrome (PCOS) is a multifactorial disorder characterized by irregular menstrual problems, hyperandrogenism, and presence of polycystic ovaries. Till date, molecular mechanism underlying PCOS remains elusive. Recently mitochondrial displacement loop (D-loop) variants have been identified to be novel players in the pathogenesis of PCOS. At present, rare variants, besides common variants, are also the focus of research as it is believed to make essential contribution to the risk of complex diseases. However, rare and low hetroplasmic variants in mitochondrial D-loop are still not investigated in PCOS women. Furthermore, variants in light-strand origin of DNA replication (OriL) of mitochondrial DNA (mtDNA) have not been explored in PCOS. Hence, in this study, we investigated rare to common mitochondrial D-loop and OriL region variants obtained using mtDNA next-generation sequencing in women with PCOS. Furthermore, we also assessed mtDNA copy number, a biomarker of mitochondrial dysfunction (MD) in women with PCOS, as the variants in mtDNA are known to be associated with low mtDNA copy number in PCOS women. A total of 67 D-loop variants including 6 novel variants were identified in 30 PCOS women. Among 67 variants, 29 variants were reported in PCOS women. A single variant, 5746A was found in OriL region in two PCOS women. Both transition and transversion variants were found but transition variants occur at very high frequency compared with transversions (82.35% vs. 17.64%, respectively). As transition variants in mtDNA are known to arise because of polymerase γ errors, occurrence of high transition rates indicates that most mutation arises because of defect in replication errors that causes mtDNA damage leading to MD. Furthermore, mtDNA copy number was found to be low in women with PCOS compared with healthy control women suggesting that MD may be the contributing factor in the pathogenesis of PCOS.

Novel homozygous mutation in the FBXL4 gene is associated with mitochondria DNA depletion syndrome-13

BackgroundMitochondrial DNA depletion syndrome-13 (MTDPS13) is caused by mutations in FBXL4 (F-box and leucine-rich repeat protein 4), a nuclear gene encoding an F-box protein that plays a role in maintaining mtDNA integrity and stability. MethodsWe identified a novel homozygous FBXL4 gene mutation, c.993dupA (p.L332Tfs*3), in a 1-year-old girl of Han Chinese descent. We performed three-dimensional protein structural analysis and targeted mtDNA next-generation sequencing. We analysed FBXL4 expression and mitochondrial DNA level, and reviewed mutations reported in FBXL4-related literature. ResultsThis mutation resulted in premature termination of translation and loss of 288 amino acids from C-terminus. A three-dimensional structural analysis revealed that conserved LRR domains were lost in mutant FBXL4 protein, which likely affected its ability to form protein-protein interactions. There were no differences in FBXL4 mRNA expression levels between the patient and her parents. There were no mtDNA mutations in either the patient or her parents. However, ND1/GAPDH ratio in lymphocytes and urine, which represents mtDNA/nuclear DNA ratio, showed that the number of mitochondrial genomes was significantly lower in the patient than in her parents or wild-type subjects. ConclusionHomozygous FBXL4 gene mutation, c.993dupA, can cause mitochondrial dysfunction, and LRR region is especially important for FBXL4 protein function.

Mitochondrial DNA mutations in Malaysian female breast cancer patients.

Cancer development has been ascribed with diverse genetic variations which are identified in both mitochondrial and nuclear genomes. Mitochondrial DNA (mtDNA) alterations have been detected in several tumours which include lung, colorectal, renal, pancreatic and breast cancer. Several studies have explored the breast tumour-specific mtDNA alteration mainly in Western population. This study aims to identify mtDNA alterations of 20 breast cancer patients in Malaysia by next generation sequencing analysis. Twenty matched tumours with corresponding normal breast tissues were obtained from female breast cancer patients who underwent mastectomy. Total DNA was extracted from all samples and the entire mtDNA (16.6kb) was amplified using long range PCR amplification. The amplified PCR products were sequenced using mtDNA next-generation sequencing (NGS) on an Illumina Miseq platform. Sequencing involves the entire mtDNA (16.6kb) from all pairs of samples with high-coverage (~ 9,544 reads per base). MtDNA variants were called and annotated using mtDNA-Server, a web server. A total of 18 of 20 patients had at least one somatic mtDNA mutation in their tumour samples. Overall, 65 somatic mutations were identified, with 30 novel mutations. The majority (59%) of the somatic mutations were in the coding region, whereas only 11% of the mutations occurred in the D-loop. Notably, somatic mutations in protein-coding regions were non-synonymous (49%) in which 15.4% of them are potentially deleterious. A total of 753 germline mutations were identified and four of which were novel mutations. Compared to somatic alterations, less than 1% of germline missense mutations are harmful. The findings of this study may enhance the current knowledge of mtDNA alterations in breast cancer. To date, the catalogue of mutations identified in this study is the first evidence of mtDNA alterations in Malaysian female breast cancer patients.

An Effective Strategy to Eliminate Inherent Cross-Contamination in mtDNA Next-Generation Sequencing of Multiple Samples

Heteroplasmic mutations in mitochondrial DNA (mtDNA) play critical roles in mitochondrial disease, aging, and cancer. Recently, next-generation sequencing (NGS) has been widely used to detect mtDNA mutations for diagnosis and monitoring of the above-mentioned diseases. However, little attention is paid on inherent cross-contamination generated during mtDNA capture and sequencing of mixed samples, which may seriously reduce the detection accuracy of mtDNA heteroplasmic mutations. In this study, a novel sequencing strategy based on a unique double-barcode design was established. The results showed that when single barcode-based analysis strategy was used, cross-contamination level of 20 DNA samples ranged from 0.27% to 11.90% on HiSeq 2500 and from 0.93% to 17.70% on HiSeq X ten, whereas double barcode-based strategy could effectively eliminate cross-contamination. Moreover, the data indicated that cross-contamination was mainly derived from capture process and was significantly affected by different NGS platforms. In addition, contamination level was negatively related to sequencing depth. Moreover, cross-contamination significantly increased the false-positive calling of mtDNA heteroplasmic mutations and remarkably affected the heteroplasmy level of mtDNA mutations. In contrast, cross-contamination had no notable effect on classification of mtDNA haplogroup. Taken together, our novel double barcode-based sequencing strategy is effective in eliminating cross-contamination, enhancing the detection accuracy of mtDNA NGS, and improving its application in diagnosis or monitoring of diseases associated with mtDNA mutations.

Mitochondrial DNA Mismatches in Thoracic Transplantation: Potential Triggers of Allograft Immunogenicity

Purpose Organ transplant (Tx) recipients of non-European ancestry experience higher rejection rates. Stem cell research shows that single nucleotide variant mismatches (SNVMs) in mitochondria DNA (mtDNA) of the donor and host cells trigger alloimmune responses. However, no prior studies have addressed mtDNA SNVMs in organ transplantation. We hypothesized that SNVMs between donor-recipient (D-R) mtDNA sequences contribute to allograft immunogenicity. Our purpose was to conduct a proof-of-concept study examining the relationship between mtDNA SNVMs and immunogenicity in human thoracic Tx recipients. Methods We assayed 364 mtDNAs in pre-Tx peripheral blood mononuclear cells (PBMCs) samples from 182 Tx recipients (19 heart and 163 lung) and their respective donors, using customized standard techniques for mtDNA next-generation sequencing to a mean depth of 9750x. MtDNA sequence reads for D-R pairs were compared to a reference: rCRS database. SNVMs were identified and analyzed for relationships to donor and recipient self-reported sex and race. To determine if nonsynonymous SNVMs trigger allo-specific immunity, we constructed peptides with and without SNVMs and ran Elispot assays to quantify Interferon gamma (INFγ) release from recipients’ re-activated PBMCs collected 1 year post-transplant. Results We identified on average 46 (SD=20) mtDNA SNVMs per D-R pair, 13 (SD=7) of which were nonsynonymous and could form allo-specific immunogenic proteins. Nonsynonymous SNVMs did not differ significantly in D-R sex-matched (N=133, M=13, SD=6) and sex-mismatched pairs (N=46, M=14, SD=7), t(72)=0.87, p = 0.386. In contrast, nonsynonymous SNVMs were significantly higher in D-R race-discordant transplant pairs (N=50, M=16, SD=9) compared to race-concordant pairs (N=100, M=12, SD=5); t(64)=3.07, p=0.003. This difference persisted even after excluding 69 mtDNA haplogroups hot spots: race-discordant pairs (N=50, M=11, SD=7) and race-concordant pairs (N=100, M=8, SD=4); t(68)=2.14, p=0.029. Elispot data showed 30-200 fold increase in INFγ release from peptides with SNVM compared to control, p Conclusion The higher frequency of D-R mtDNA SNVMs in race-discordant transplant pairs, and the immunogenicity of peptides generated from patients with mtDNA SNVMs supports the hypothesis that SNVMs may contribute to allograft immunogenicity.

Natural underlying mtDNA heteroplasmy as a potential source of intra-person hiPSC variability.

Functional variability among human clones of induced pluripotent stem cells (hiPSCs) remains a limitation in assembling high-quality biorepositories. Beyond inter-person variability, the root cause of intra-person variability remains unknown. Mitochondria guide the required transition from oxidative to glycolytic metabolism in nuclear reprogramming. Moreover, mitochondria have their own genome (mitochondrial DNA [mtDNA]). Herein, we performed mtDNA next-generation sequencing (NGS) on 84 hiPSC clones derived from a cohort of 19 individuals, including mitochondrial and non-mitochondrial patients. The analysis of mtDNA variants showed that low levels of potentially pathogenic mutations in the original fibroblasts are revealed through nuclear reprogramming, generating mutant hiPSCs with a detrimental effect in their differentiated progeny. Specifically, hiPSC-derived cardiomyocytes with expanded mtDNA mutations non-related with any described human disease, showed impaired mitochondrial respiration, being a potential cause of intra-person hiPSC variability. We propose mtDNA NGS as a new selection criterion to ensure hiPSC quality for drug discovery and regenerative medicine.

MtDNA-Server: next-generation sequencing data analysis of human mitochondrial DNA in the cloud.

Next generation sequencing (NGS) allows investigating mitochondrial DNA (mtDNA) characteristics such as heteroplasmy (i.e. intra-individual sequence variation) to a higher level of detail. While several pipelines for analyzing heteroplasmies exist, issues in usability, accuracy of results and interpreting final data limit their usage. Here we present mtDNA-Server, a scalable web server for the analysis of mtDNA studies of any size with a special focus on usability as well as reliable identification and quantification of heteroplasmic variants. The mtDNA-Server workflow includes parallel read alignment, heteroplasmy detection, artefact or contamination identification, variant annotation as well as several quality control metrics, often neglected in current mtDNA NGS studies. All computational steps are parallelized with Hadoop MapReduce and executed graphically with Cloudgene. We validated the underlying heteroplasmy and contamination detection model by generating four artificial sample mix-ups on two different NGS devices. Our evaluation data shows that mtDNA-Server detects heteroplasmies and artificial recombinations down to the 1% level with perfect specificity and outperforms existing approaches regarding sensitivity. mtDNA-Server is currently able to analyze the 1000G Phase 3 data (n = 2,504) in less than 5 h and is freely accessible at https://mtdna-server.uibk.ac.at.

Validation of Next-Generation Sequencing of Entire Mitochondrial Genomes and the Diversity of Mitochondrial DNA Mutations in Oral Squamous Cell Carcinoma.

BackgroundOral squamous cell carcinoma (OSCC) is mainly caused by smoking and alcohol abuse and shows a five-year survival rate of ~50%. We aimed to explore the variation of somatic mitochondrial DNA (mtDNA) mutations in primary oral tumors, recurrences and metastases.MethodsWe performed an in-depth validation of mtDNA next-generation sequencing (NGS) on an Illumina HiSeq 2500 platform for its application to cancer tissues, with the goal to detect low-level heteroplasmies and to avoid artifacts. Therefore we genotyped the mitochondrial genome (16.6 kb) from 85 tissue samples (tumors, recurrences, resection edges, metastases and blood) collected from 28 prospectively recruited OSCC patients applying both Sanger sequencing and high-coverage NGS (~35,000 reads per base).ResultsWe observed a strong correlation between Sanger sequencing and NGS in estimating the mixture ratio of heteroplasmies (r = 0.99; p<0.001). Non-synonymous heteroplasmic variants were enriched among cancerous tissues. The proportions of somatic and inherited variants in a given gene region were strongly correlated (r = 0.85; p<0.001). Half of the patients shared mutations between benign and cancerous tissue samples. Low level heteroplasmies (<10%) were more frequent in benign samples compared to tumor samples, where heteroplasmies >10% were predominant. Four out of six patients who developed a local tumor recurrence showed mutations in the recurrence that had also been observed in the primary tumor. Three out of five patients, who had tumor metastases in the lymph nodes of their necks, shared mtDNA mutations between primary tumors and lymph node metastases. The percentage of mutation heteroplasmy increased from the primary tumor to lymph node metastases.ConclusionsWe conclude that Sanger sequencing is valid for heteroplasmy quantification for heteroplasmies ≥10% and that NGS is capable of reliably detecting and quantifying heteroplasmies down to the 1%-level. The finding of shared mutations between primary tumors, recurrences and metastasis indicates a clonal origin of malignant cells in oral cancer.

An integrated pipeline for next generation sequencing and annotation of the complete mitochondrial genome of the giant intestinal fluke,Fasciolopsis buski(Lankester, 1857) Looss, 1899

Helminths include both parasitic nematodes (roundworms) and platyhelminths (trematode and cestode flatworms) that are abundant, and are of clinical importance. The genetic characterization of parasitic flatworms using advanced molecular tools is central to the diagnosis and control of infections. Although the nuclear genome houses suitable genetic markers (e.g., in ribosomal (r) DNA) for species identification and molecular characterization, the mitochondrial (mt) genome consistently provides a rich source of novel markers for informative systematics and epidemiological studies. In the last decade, there have been some important advances in mtDNA genomics of helminths, especially lung flukes, liver flukes and intestinal flukes. Fasciolopsis buski, often called the giant intestinal fluke, is one of the largest digenean trematodes infecting humans and found primarily in Asia, in particular the Indian subcontinent. Next-generation sequencing (NGS) technologies now provide opportunities for high throughput sequencing, assembly and annotation within a short span of time. Herein, we describe a high-throughput sequencing and bioinformatics pipeline for mt genomics for F. buski that emphasizes the utility of short read NGS platforms such as Ion Torrent and Illumina in successfully sequencing and assembling the mt genome using innovative approaches for PCR primer design as well as assembly. We took advantage of our NGS whole genome sequence data (unpublished so far) for F. buski and its comparison with available data for the Fasciola hepatica mtDNA as the reference genome for design of precise and specific primers for amplification of mt genome sequences from F. buski. A long-range PCR was carried out to create an NGS library enriched in mt DNA sequences. Two different NGS platforms were employed for complete sequencing, assembly and annotation of the F. buski mt genome. The complete mt genome sequences of the intestinal fluke comprise 14,118 bp and is thus the shortest trematode mitochondrial genome sequenced to date. The noncoding control regions are separated into two parts by the tRNA-Gly gene and don’t contain either tandem repeats or secondary structures, which are typical for trematode control regions. The gene content and arrangement are identical to that of F. hepatica. The F. buski mtDNA genome has a close resemblance with F. hepatica and has a similar gene order tallying with that of other trematodes. The mtDNA for the intestinal fluke is reported herein for the first time by our group that would help investigate Fasciolidae taxonomy and systematics with the aid of mtDNA NGS data. More so, it would serve as a resource for comparative mitochondrial genomics and systematic studies of trematode parasites.