Entire Mitochondrial DNA Sequencing on Massively Parallel Sequencing for the Korean Population

High prevalence of m.1555A > G in patients with hearing loss in the Baikal Lake region of Russia as a result of founder effect

Mitochondrial forms account approximately 1-2% of all nonsyndromic cases of hearing loss (HL). One of the most common causative variants of mtDNA is the m.1555A > G variant of the MT-RNR1 gene (OMIM 561000). Currently the detection of the m.1555A > G variant of the MT-RNR1 gene is not included in all research protocols. In this study this variant was screened among 165 patients with HL from the Republic of Buryatia, located in the Baikal Lake region of Russia. In our study, the total contribution of the m.1555A > G variant to the etiology of HL was 12.7% (21/165), while the update global prevalence of this variant is 1.8% (863/47,328). The m.1555A > G variant was notably more prevalent in Buryat (20.2%) than in Russian patients (1.3%). Mitogenome analysis in 14 unrelated Buryat families carrying the m.1555A > G variant revealed a predominant lineage: in 13 families, a cluster affiliated with sub-haplogroup A5b (92.9%) was identified, while one family had the D5a2a1 lineage (7.1%). In a Russian family with the m.1555A > G variant the lineage affiliated with sub-haplogroup F1a1d was found. Considering that more than 90% of Buryat families with the m.1555A > G variant belong to the single maternal lineage cluster we conclude that high prevalence of this variant in patients with HL in the Baikal Lake region can be attributed to a founder effect.

Genetic Diversity Analysis of the Chinese Daur Ethnic Group in Heilongjiang Province by Complete Mitochondrial Genome Sequencing

Mitochondrial DNA (mtDNA) has the characteristics of maternal inheritance, high mutation rate, high copy number, and no recombination. As the most powerful tool for studying the origin and evolution of modern humans, mtDNA has great significance in the research of population genetics and evolutionary genetics. Here, we provide new insights into the maternal genetic history of the Daur ethnic group by generating complete mitochondrial genomes from a total of 146 Daur individuals in China. We also collected the published complete mitochondrial genome sequences of 5,094 individuals from 56 worldwide populations as reference data to further explore the matrilineal genetic landscape of the Daur ethnic group. First, the haplotype diversity was 0.9943 ± 0.0019 and nucleotide diversity was 0.0428 ± 0.0210. The neutrality tests of the Daur group showed significant negative values and the mismatch distribution curve was obviously distributed in a unimodal pattern. The results showed that the Daur ethnic group has high genetic diversity and may have experienced recent population expansion. In addition, the main haplogroups of the Daur population were haplogroup D (31.51%), M* (20.55%), C (10.28%), F (7.53%), and B (6.85%), all of which were prevalent in northern China. It probably implies the northern Chinese origin of the Daur population. The PCA, F_ST, and phylogenetic analysis results indicated that the Daur group formed a cluster with East Asian populations, and had few genetic differences with the populations in northern China. More importantly, we found that disease-related mutation sites of the mitochondrial genome may be related to ethnic groups, which may have important implications for the prevention and occurrence of specific diseases. Overall, this study revealed the complexity and diversity of the matrilineal genetic background of the Daur ethnic group. Meanwhile, it provided meaningful data for the research on the diversity of the human genome.

The frequency of the known mitochondrial variants associated with drug-induced toxicity in a Korean population

BACKGROUND

Few studies have annotated the whole mitochondrial DNA (mtDNA) genome associated with drug responses in Asian populations. This study aimed to characterize mtDNA genetic profiles, especially the distribution and frequency of well-known genetic biomarkers associated with diseases and drug-induced toxicity in a Korean population.

METHOD

Whole mitochondrial genome was sequenced for 118 Korean subjects by using a next-generation sequencing approach. The bioinformatic pipeline was constructed for variant calling, haplogroup classification and annotation of mitochondrial mutation.

RESULTS

A total of 681 variants was identified among all subjects. The MT-TRNP gene and displacement loop showed the highest numbers of variants (113 and 74 variants, respectively). The m.16189T > C allele, which is known to reduce the mtDNA copy number in human cells was detected in 25.4% of subjects. The variants (m.2706A > G, m.3010A > G, and m.1095T > C), which are associated with drug-induced toxicity, were observed with the frequency of 99.15%, 30.51%, and 0.08%, respectively. The m.2150T > A, a genotype associated with highly disruptive effects on mitochondrial ribosomes, was identified in five subjects. The D and M groups were the most dominant groups with the frequency of 34.74% and 16.1%, respectively.

CONCLUSIONS

Our finding was consistent with Korean Genome Project and well reflected the unique profile of mitochondrial haplogroup distribution. It was the first study to annotate the whole mitochondrial genome with drug-induced toxicity to predict the ADRs event in clinical implementation for Korean subjects. This approach could be extended for further study for validation of the potential ethnic-specific mitochondrial genetic biomarkers in the Korean population.

Mitochondrial genome sequencing with short overlapping amplicons on MiSeq FGx system

With the development and maturation of massively parallel sequencing (MPS) technology, the mitochondrial genome (mitogenome) sequencing is increasingly applied in the forensic field. In this study, we employed the strategy of short overlapping amplicons for the whole mitogenome, library preparation with tagmentation using the Nextera® XT DNA Library Preparation Kit, sequencing on the MiSeq FGxTM Forensic Genomics System and analyzing data using the mitochondrial(mtDNA) MSR Plug-in and the mtDNA Variant Analyzer. A total of 27 libraries and 56 libraries were sequenced in a run using MiSeq Reagent Kit v2 and v3, respectively. Results showed more than 1800 × of averaged depth of coverage (DoC) at each position. Concordant haplotypes of 9947 A and 2800 M were obtained at 32 variants. Cross-reactivity was observed with 1 ng primate DNA and 10 ng non-primate DNA but could be easily distinguished. Full and accurate variants were obtained from at least 50 pg input DNA and from minor contributors between 19:1 and 1:19 mixed ratios with known reference profiles. More than 86% variants were detected from ≥200-bp degraded samples but its haplotype was assigned to more ancestral haplogroup. Further, a total of 3 962 variants were observed at 613 nucleotide positions from 103 Xibe mitogenomes with 25:1 ratio of transitions to transversions. Two new transversions (C13735A and A14755C) and two tri-alleles at nps 9824 and 16092 were identified. There were 103 unique mitogenome haplotypes from 103 Chinese Xibe that were assigned to 79 haplogroups. Haplogroup D was the preponderant top-level haplogroup in Xibe followed by F, B, M, A, N, G, C, Z, Y, HV and J. Random match probability (RMP) and haplotype diversity (HD) of the whole mitogenome was calculated as 0.0097 and 1.0000, respectively. Compared with HVS-I only, RMP decreased 33.56%, while the number of haplotypes and HD increased 15.73% and 0.49%, respectively. Principal component analysis (PCA) showed that Xibe was clustered to East and Southeast Asian. As a whole, this MPS strategy is suitable for the whole mitogenome sequencing especially for degraded samples and can facilitate generating mitogenome data to support the routine application in forensic sciences. EMP00726 is the first whole mitogenome dataset from Xibe contributed to the EMPOP. Supplemental data for this article are available online at.

Massively parallel sequencing of human skeletal remains in Vietnam using the precision ID mtDNA control region panel on the Ion S5™ system

Mitochondrial DNA (mtDNA) analysis using Sanger sequencing has been a routine practice for the identification of human skeletal remains. However, this process is usually challenging since DNA from the remains is highly degraded and at low concentration. Recently, the advent and implementation of massively parallel sequencing (MPS) have been offered the ability to improve mtDNA sequence data for forensic analysis. To assess the utility of the Ion S5™ system - an MPS platform for mtDNA analysis in challenging samples, we sequenced the mitochondrial control region of 52 age-old skeletal remains. Using the Precision ID mtDNA Control Region Panel, 50 full and two partial control region haplotypes at relatively high mean coverage of 2494 × were achieved for variant calling. Further variant analysis at 10% threshold for point heteroplasmy showed high degradation degree in terms of DNA damage in our bone samples. A higher point heteroplasmy threshold of 20% was required to diminish most of background noise caused by the damage. The results from this study indicated the potential application of the Ion S5™ system in sequencing degraded samples in Vietnam and provided valuable data sources for forensic analyses in the future.

Investigation of control region sequences of mtDNA in Naqu Tibetan population from Northwestern China

BACKGROUND

The sequence polymorphisms of mitochondrial DNA (mtDNA) are valuable in forensic medicine and anthropological genetics.

AIM

We analysed the sequences of the mtDNA control region in 207 unrelated Tibetan individuals from the Naqu region, Tibet Autonomous Region in the People's Republic of China, and investigated the population structure of the region by population comparison with other groups.

SUBJECTS AND METHODS

Genomic DNA was extracted and hypervariable regions (HVS-I and HVS-II) were amplified and sequenced. Subsequently, sequences were aligned and compared with the revised Cambridge sequence. Moreover, population comparison was performed between the Naqu Tibetan group and the other groups.

CONCLUSION

Our study provided available data for exploring the mtDNA haplotype of the Tibetan population in the Naqu region, and population comparisons found that the Naqu Tibetan population has its own unique structure.

Delineation of Mitochondrial DNA Variants From Exome Sequencing Data and Association of Haplogroups With Obesity in Kuwait

BACKGROUND/OBJECTIVES

Whole-exome sequencing is a valuable tool to determine genetic variations that are associated with rare and common health conditions. A limited number of studies demonstrated that mitochondrial DNA can be captured using whole-exome sequencing. Previous studies have suggested that mitochondrial DNA variants and haplogroup lineages are associated with obesity. Therefore, we investigated the role of mitochondrial variants and haplogroups contributing to the risk of obesity in Arabs in Kuwait using exome sequencing data.

SUBJECTS/METHODS

Indirect mitochondrial genomes were extracted from exome sequencing data from 288 unrelated native Arab individuals from Kuwait. The cohort was divided into obese [body mass index (BMI) ≥ 30 kg/m2] and non-obese (BMI < 30 kg/m2) groups. Mitochondrial variants were identified, and haplogroups were classified and compared with other sequencing technologies. Statistical analysis was performed to determine associations and identify mitochondrial variants and haplogroups affecting obesity.

RESULTS

Haplogroup R showed a protective effect on obesity [odds ratio (OR) = 0.311; P = 0.006], whereas haplogroup L individuals were at high risk of obesity (OR = 2.285; P = 0.046). Significant differences in mitochondrial variants between the obese and non-obese groups were mainly haplogroup-defining mutations and were involved in processes in energy generation. The majority of mitochondrial variants and haplogroups extracted from exome were in agreement with technical replica from Sanger and whole-genome sequencing.

CONCLUSIONS

This is the first to utilize whole-exome data to extract entire mitochondrial haplogroups to study its association with obesity in an Arab population.

Mitochondrial haplogroup J associated with higher risk of obesity in the Qatari population

Obesity, a major risk factor for metabolic disorders, is highly prevalent in Qatari population. Maternal transmission of obesity traits can be significant; for example, X haplogroup is known to be associated with lower BMI and body fat mass in Northern Europeans and T haplogroup which is a sister haplogroup of J is known to be associated with obesity in Caucasian subjects from Austria and Southern Italy. We aimed to delineate the mitochondrial haplogroups and variants associated with obesity in Qatari population. Mitochondrial genomes of 864 Qatari individuals were extracted from whole exome sequencing data with an average coverage of 77X. We distributed the participants into 2 sub-cohorts: obese (BMI ≥ 30) and non-obese (BMI < 30); the mean value of BMI from these two groups were 36.5 ± 5.7 and 26.5 ± 2.6, respectively. Mitochondrial haplogroup profiling followed by uni- and multivariant association tests adjusted for covariates were performed. Qatari individuals with mitochondrial haplogroup J had an increased (twofold) risk of obesity (odds ratio [OR] 1.925; 95% CI 1.234–3.002; P = 0.0038; the Bonferroni adjusted P value threshold is 0.0041), whereas the individuals with haplogroup X were at low risk of obesity (OR 0.387; 95% CI 0.175–0.857; P = 0.019). Further, a set of 38 mitochondrial variants were found to be associated (at P ≤ 0.05) with obesity in models adjusted for age, sex and haplogroup.

Searching the undetected mtDNA variants in forensic MPS data

The efficiency of MPS in forensic mtDNA analysis has been thoroughly proven, although a reliable and well established data evaluation still remains a critical point. Numerous bioinformatics tools have been developed, but most of them require specific operating systems and high costs, while free open-source programs with user-friendly interfaces are few. In this study, 43 full mtGenomes were sequenced using the Ion Personal Genome Machine™ (PGM™) System and analyzed utilizing the plug-in Variant Caller (TVC) of the Ion Torrent Software Suite and the mtDNA-Server (mDS), a free web-based mitochondrial analysis tool for MPS data. The outcomes of these two different analysis tools were compared to variants noted after manual inspection of the aligned reads performed using Integrative Genomics Viewer (IGV). The comparison highlighted the presence of thirty-nine discordant variant calls, which were resolved by Sanger sequencing that confirmed the presence of all variants, except for 7 deletions. The combined adoption of IGV and Sanger type sequencing confirmatory steps, in addition of TVC and mDS analysis, resulted in a more accurate variants assignment with the detection of 32 additional true polymorphisms, which were noted in the final dataset. Regarding the heteroplasmy issue, out of a total of thirty heteroplasmic variants, twenty-eight were detected by the TVC, while the mDS detected twenty-two. Overall, none of the used bioinformatics tools were the perfect choice and a secondary analysis with an expert's opinion in complete mtGenome MPS data evaluation is still required in forensic genetic analysis.

Genetic legacy of cultures indigenous to the Northeast Asian coast in mitochondrial genomes of nearly extinct maritime tribes

Background

We have described the diversity of complete mtDNA sequences from ‘relic’ groups of the Russian Far East, primarily the Nivkhi (who speak a language isolate with no clear relatedness to any others) and Oroki of Sakhalin, as well as the sedentary Koryak from Kamchatka and the Udegey of Primorye. Previous studies have shown that most of their traditional territory was dramatically reshaped by the expansion of Tungusic-speaking groups.

Results

Overall, 285 complete mitochondrial sequences were selected for phylogenetic analyses of published, revised and new mitogenomes. To highlight the likely role of Neolithic expansions in shaping the phylogeographical landscape of the Russian Far East, we focus on the major East Eurasian maternal lineages (Y1a, G1b, D4m2, D4e5, M7a2, and N9b) that are restricted to the coastal area. To obtain more insight into autochthonous populations, we removed from the phylogeographic analysis the G2a, G3a2, M8a1, M9a1, and C4b1 lineages, also found within our samples, likely resulting from admixture between the expanding proto-Tungus and the indigenous Paleoasiatic groups with whom they assimilated. Phylogenetic analysis reveals that unlike the relatively diverse lineage spectrum observed in the Amur estuary and northwestern Sakhalin, the present-day subpopulation on the northeastern coast of the island is relatively homogenous: a sole Y1a sublineage, conspicuous for its nodal mutation at m.16189 T > C!, includes different haplotypes. Sharing of the Y1a-m.16189 T > C! sublineages and haplotypes among the Nivkhi, Ulchi and sedentary Koryak is also evident. Aside from Y1a, the entire tree approach expands our understanding of the evolutionary history of haplogroups G1, D4m, N9b, and M7a2. Specifically, we identified the novel haplogroup N9b1 in Primorye, which implies a link between a component of the Udegey ancestry and the Hokkaido Jomon.

Conclusions

Through a comprehensive dataset of mitochondrial genomes retained in autochthonous populations along the coast between Primorye and the Bering Strait, we considerably extended the sequence diversity of these populations to provide new features based on the number and timing of founding lineages. We emphasize the value of integrating genealogical information with genetic data for reconstructing the population history of indigenous groups dramatically impacted by twentieth century resettlement and social upheavals.

Validation of NGS for mitochondrial DNA casework at the FBI Laboratory

As a first step towards integrating next generation sequencing (NGS) technology into the FBI Laboratory's operational casework, the PowerSeq™ CRM Nested System, an NGS-based mitochondrial DNA (mtDNA) control region assay, was developmentally and internally validated. The validation studies were conducted in accordance with the Scientific Working Group on DNA Analysis Methods (SWGDAM) Validation Guidelines for Forensic DNA Analysis Methods, and the FBI's Quality Assurance Standards (QAS) for Forensic DNA Testing Laboratories. The assay was shown to be highly reproducible, with variant frequencies across intra and inter-run replicates of the same sample differing, on average, by just 0.3% for substitutions and point heteroplasmies and 1.5% for insertions and deletions. The assay was also shown to be extremely sensitive, yielding complete control region sequence data from as few as 2000 copies of mtDNA. This is a more than 20-fold increase in sensitivity when compared to the FBI Laboratory's current Sanger sequencing-based protocols and, based on mtDNA quantitation values of samples routinely encountered in mtDNA casework, suggests that the percentage of questioned samples from which full control region data can be recovered will increase from our current 20% to approximately 90% success with NGS technology. In addition, the assay requires on average only 30% of the extract volume typically required to develop control region profiles from degraded samples via Sanger sequencing. Overall, these studies establish the reliability of the PowerSeq™ CRM Nested System for accurate mtDNA control region typing and can serve as a model for laboratories seeking to validate NGS protocols for forensic mtDNA analysis.

Mitochondrial DNA in human identification: a review

Mitochondrial DNA (mtDNA) presents several characteristics useful for forensic studies, especially related to the lack of recombination, to a high copy number, and to matrilineal inheritance. mtDNA typing based on sequences of the control region or full genomic sequences analysis is used to analyze a variety of forensic samples such as old bones, teeth and hair, as well as other biological samples where the DNA content is low. Evaluation and reporting of the results requires careful consideration of biological issues as well as other issues such as nomenclature and reference population databases. In this work we review mitochondrial DNA profiling methods used for human identification and present their use in the main cases of humanidentification focusing on the most relevant issues for forensics.

Full mtDNA genome sequencing of Brazilian admixed populations: A forensic-focused evaluation of a MPS application as an alternative to Sanger sequencing methods

The use of Massive Parallel Sequencing (MPS) techniques have been proposed by the forensic community as an alternative to Sanger sequencing methods in routine forensic casework analysis regarding mitochondrial DNA (mtDNA). Interesting features of MPS include high throughput, ability to simultaneously genotype a significant number of samples by barcoding techniques, processing automation, reduced time and costs, among others. Advantages include the capability of generating full mtDNA genome sequences versus usual techniques, usually limited to hypervariable or control regions exclusively. In this work, 96 reference single-source samples from three different Brazilian cities were subjected to full mtDNA genome sequencing by MPS techniques using an early-access version of Precision ID mtDNA Whole Genome Panel on an Ion Torrent PGM platform (Thermo Fisher Scientific, Waltham, MA, USA). Complete, high-quality sequences were obtained and sequencing performance was evaluated via four different metrics. As a subset of evaluated samples have been previously submitted for Sanger sequencing of the control region, a comparative analysis of both methods' results was conducted in order to compare technique adequacy within a forensic context. Even though this study is one of the first to report full mtDNA genome sequences for Brazilian admixed populations, the observed haplotypes exhibit a predominance of Native American and African maternal lineages in the studied sample set, reproducing results described in the literature for control regions only. Interpopulation analysis among Brazilian and 26 worldwide populations was also carried out. The results indicate that MPS-generated full mtDNA genome sequences may have great utility in forensic real casework applications, with a pronounced gain of genetic information and discrimination power provided by coding region evaluation and the enhanced capacity of heteroplasmies determination. Database construction and other relevant factors concerning implementation of such techniques in Brazilian forensic laboratories are also discussed.

Detection of Innate and Artificial Mitochondrial DNA Heteroplasmy by Massively Parallel Sequencing: Considerations for Analysis

BACKGROUND

Mitochondrial heteroplasmy, the co-existence of different mitochondrial polymorphisms within an individual, has various forensic and clinical implications. But there is still no guideline on the application of massively parallel sequencing (MPS) in heteroplasmy detection. We present here some critical issues that should be considered in heteroplasmy studies using MPS.

METHODS

Among five samples with known innate heteroplasmies, two pairs of mixture were generated for artificial heteroplasmies with target minor allele frequencies (MAFs) ranging from 50% to 1%. Each sample was amplified by two-amplicon method and sequenced by Ion Torrent system. The outcomes of two different analysis tools, Torrent Suite Variant Caller (TVC) and mtDNA-Server (mDS), were compared.

RESULTS

All the innate heteroplasmies were detected correctly by both analysis tools. Average MAFs of artificial heteroplasmies correlated well to the target values. The detection rates were almost 90% for high-level heteroplasmies, but decreased for low-level heteroplasmies. TVC generally showed lower detection rates than mDS, which seems to be due to their own computation algorithms which drop out some reference-dominant heteroplasmies. Meanwhile, mDS reported several unintended low-level heteroplasmies which were suggested as nuclear mitochondrial DNA sequences. The average coverage depth of each sample placed on the same chip showed considerable variation. The increase of coverage depth had no effect on the detection rates.

CONCLUSION

In addition to the general accuracy of the MPS application on detecting heteroplasmy, our study indicates that the understanding of the nature of mitochondrial DNA and analysis algorithm would be crucial for appropriate interpretation of MPS results.

Evaluation of the precision ID whole MtDNA genome panel for forensic analyses

Mitochondrial DNA (mtDNA) amplification and Massively Parallel Sequencing (MPS) using an early access version of the Precision ID Whole MtDNA Genome Panel (Thermo Fisher Scientific) and the Ion Personal Genome Machine (PGM) were evaluated using 15 forensically relevant samples. Samples were selected to represent typical forensic specimens for mtDNA analysis including hairs, hair shafts, swabs and ancient solid tissue samples (bones and teeth) that were stored in the freezer for up to several years after having been typed with conventional Sanger-type Sequencing and Capillary Electrophoresis. The MPS haplotypes confirmed the earlier results in all samples and provided additional sequence information that improved discrimination power and haplogroup estimation. The results raised the appetite for further experiments to validate and apply the new technology in forensic practice.

Massively parallel sequencing-enabled mixture analysis of mitochondrial DNA samples

The mitochondrial genome has a number of characteristics that provide useful information to forensic investigations. Massively parallel sequencing (MPS) technologies offer improvements to the quantitative analysis of the mitochondrial genome, specifically the interpretation of mixed mitochondrial samples. Two-person mixtures with nuclear DNA ratios of 1:1, 5:1, 10:1, and 20:1 of individuals from different and similar phylogenetic backgrounds and three-person mixtures with nuclear DNA ratios of 1:1:1 and 5:1:1 were prepared using the Precision ID mtDNA Whole Genome Panel and Ion Chef, and sequenced on the Ion PGM or Ion S5 sequencer (Thermo Fisher Scientific, Waltham, MA, USA). These data were used to evaluate whether and to what degree MPS mixtures could be deconvolved. Analysis was effective in identifying the major contributor in each instance, while SNPs from the minor contributor's haplotype only were identified in the 1:1, 5:1, and 10:1 two-person mixtures. While the major contributor was identified from the 5:1:1 mixture, analysis of the three-person mixtures was more complex, and the mixed haplotypes could not be completely parsed. These results indicate that mixed mitochondrial DNA samples may be interpreted with the use of MPS technologies.

Performance of the Early Access AmpliSeq™ Mitochondrial Panel with degraded DNA samples using the Ion Torrent™ platform

The Early Access AmpliSeq™ Mitochondrial Panel amplifies whole mitochondrial genomes for phylogenetic and kinship identifications, using Ion Torrent™ technology. There is currently limited information on its performance with degraded DNA, a common occurrence in forensic samples. This study evaluated the performance of the Panel with DNA samples degraded in vitro, to mimic conditions commonly found in forensic investigations. Purified DNA from five individuals was heat-treated at five time points each (125°C for 0, 30, 60, 120, and 240 min; total n = 25). The quality of DNA was assessed via a real-time DNA assay of genomic DNA and prepared for massively parallel sequencing on the Ion Torrent™ platform. Mitochondrial sequences were obtained for all samples and had an amplicon coverage averaging between 66X to 2803X. Most amplicons (157/162) displayed high coverages (452 ± 333X), while reads with less than 100X coverage were recorded in five amplicons only (90 ± 5X). Amplicon coverage was decreased with prolonged heating. At 72% strand balance, reads were well balanced between forward and reverse strands. Using a coverage threshold of ten reads per SNP, complete sequences were recovered in all samples and resolved kinship and, haplogroup relations. Additionally, the HV1 and HV2 regions of the reference and 240-min heat-treated samples (n = 10) were Sanger-sequenced for concordance. Overall, this study demonstrates the efficacy of a novel forensic Panel that recovers high quality mitochondrial sequences from degraded DNA samples.

Mitogenomic differences between the normal and tumor cells of colorectal cancer patients

So far, a reliable spectrum of mitochondrial DNA mutations in colorectal cancer cells is still unknown, and neither is their significance in carcinogenesis. Indeed, it remains debatable whether mtDNA mutations are "drivers" or "passengers" of colorectal carcinogenesis. Thus, we analyzed 200 mitogenomes from normal and cancer tissues of 100 colorectal cancer patients. Minority variant mutations were detected at the 1% level. We showed that somatic mutations frequently occur in colorectal cancer cells (75%) and are randomly distributed across the mitochondrial genome. Mutational signatures of somatic mitogenome mutations suggest that they might arise through nucleotide deamination due to oxidative stress. The majority of somatic mutations localized within the coding region (in positions not known from the human phylogeny) and was potentially pathogenic to cell metabolism. Further analysis suggested that the relaxation of negative selection in the mitogenomes of colorectal cancer cells may allow accumulation of somatic mutations. Thus, a shift in glucose metabolism from oxidative phosphorylation to glycolysis may create advantageous conditions for accumulation of mtDNA mutations. Considering the fact that the presence of somatic mtDNA mutations was not associated with any clinicopathological features, we suggested that mtDNA somatic mutations are "passengers" rather than the cause of colorectal carcinogenesis.

Assessment of mitochondrial DNA heteroplasmy detected on commercial panel using MPS system with artificial mixture samples

Mitochondrial DNA (mtDNA) heteroplasmy is a potential genetic marker for forensic mtDNA analysis as well as phylogenic studies. Frequency of mtDNA heteroplasmy has been investigated in different populations through massively parallel sequencing (MPS) analysis, revealing various levels of frequency based on different MPS systems. For accurate heteroplasmy identification, it is essential to explore reliable detection threshold on various MPS systems. In addition, software solutions and pipelines need to be evaluated to analyze data effectively. In this study, heteroplasmy analysis was conducted on a commercially available mtDNA analysis system developed for forensic caseworks with artificially mixed DNA samples known for ratios and variant positions for assessment. mtDNA heteroplasmy > 10% was detectable with Torrent Variant Caller (TVC) while lower levels were identified using GeneMarker® HTS specialized software for minor variant detection. This study implies that analytical parameters and tools need to be optimized and evaluated for low-level heteroplasmy identification. Automated system with simple and efficient workflow is needed for forensic caseworks.