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

Background

Oral 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.

Methods

We 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).

Results

We 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.

Conclusions

We 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.     

Age-Related and Heteroplasmy-Related Variation in Human mtDNA Copy Number

The mitochondrial (mt) genome is present in many copies in human cells, and intra-individual variation in mtDNA sequences is known as heteroplasmy. Recent studies found that heteroplasmies are highly tissue-specific, site-specific, and allele-specific, however the functional implications have not been explored. This study investigates variation in mtDNA copy numbers (mtCN) in 12 different tissues obtained at autopsy from 152 individuals (ranging in age from 3 days to 96 years). Three different methods to estimate mtCN were compared: shotgun sequencing (in 4 tissues), capture-enriched sequencing (in 12 tissues) and droplet digital PCR (ddPCR, in 2 tissues). The highest precision in mtCN estimation was achieved using shotgun sequencing data. However, capture-enrichment data provide reliable estimates of relative (albeit not absolute) mtCNs. Comparisons of mtCN from different tissues of the same individual revealed that mtCNs in different tissues are, with few exceptions, uncorrelated. Hence, each tissue of an individual seems to regulate mtCN in a tissue-related rather than an individual-dependent manner. Skeletal muscle (SM) samples showed an age-related decrease in mtCN that was especially pronounced in males, while there was an age-related increase in mtCN for liver (LIV) samples. MtCN in SM samples was significantly negatively correlated with both the total number of heteroplasmic sites and with minor allele frequency (MAF) at two heteroplasmic sites, 408 and 16327. Heteroplasmies at both sites are highly specific for SM, accumulate with aging and are part of functional elements that regulate mtDNA replication. These data support the hypothesis that selection acting on these heteroplasmic sites is reducing mtCN in SM of older individuals.     

Mitochondrial Mutations in Subjects with Psychiatric Disorders

A considerable body of evidence supports the role of mitochondrial dysfunction in psychiatric disorders and mitochondrial DNA (mtDNA) mutations are known to alter brain energy metabolism, neurotransmission, and cause neurodegenerative disorders. Genetic studies focusing on common nuclear genome variants associated with these disorders have produced genome wide significant results but those studies have not directly studied mtDNA variants. The purpose of this study is to investigate, using next generation sequencing, the involvement of mtDNA variation in bipolar disorder, schizophrenia, major depressive disorder, and methamphetamine use. MtDNA extracted from multiple brain regions and blood were sequenced (121 mtDNA samples with an average of 8,800x coverage) and compared to an electronic database containing 26,850 mtDNA genomes. We confirmed novel and rare variants, and confirmed next generation sequencing error hotspots by traditional sequencing and genotyping methods. We observed a significant increase of non-synonymous mutations found in individuals with schizophrenia. Novel and rare non-synonymous mutations were found in psychiatric cases in mtDNA genes: ND6, ATP6, CYTB, and ND2. We also observed mtDNA heteroplasmy in brain at a locus previously associated with schizophrenia (T16519C). Large differences in heteroplasmy levels across brain regions within subjects suggest that somatic mutations accumulate differentially in brain regions. Finally, multiplasmy, a heteroplasmic measure of repeat length, was observed in brain from selective cases at a higher frequency than controls. These results offer support for increased rates of mtDNA substitutions in schizophrenia shown in our prior results. The variable levels of heteroplasmic/multiplasmic somatic mutations that occur in brain may be indicators of genetic instability in mtDNA.     

Assessing Mitochondrial DNA Variation and Copy Number in Lymphocytes of ~2,000 Sardinians Using Tailored Sequencing Analysis Tools

DNA sequencing identifies common and rare genetic variants for association studies, but studies typically focus on variants in nuclear DNA and ignore the mitochondrial genome. In fact, analyzing variants in mitochondrial DNA (mtDNA) sequences presents special problems, which we resolve here with a general solution for the analysis of mtDNA in next-generation sequencing studies. The new program package comprises 1) an algorithm designed to identify mtDNA variants (i.e., homoplasmies and heteroplasmies), incorporating sequencing error rates at each base in a likelihood calculation and allowing allele fractions at a variant site to differ across individuals; and 2) an estimation of mtDNA copy number in a cell directly from whole-genome sequencing data. We also apply the methods to DNA sequence from lymphocytes of ~2,000 SardiNIA Project participants. As expected, mothers and offspring share all homoplasmies but a lesser proportion of heteroplasmies. Both homoplasmies and heteroplasmies show 5-fold higher transition/transversion ratios than variants in nuclear DNA. Also, heteroplasmy increases with age, though on average only ~1 heteroplasmy reaches the 4% level between ages 20 and 90. In addition, we find that mtDNA copy number averages ~110 copies/lymphocyte and is ~54% heritable, implying substantial genetic regulation of the level of mtDNA. Copy numbers also decrease modestly but significantly with age, and females on average have significantly more copies than males. The mtDNA copy numbers are significantly associated with waist circumference (p-value = 0.0031) and waist-hip ratio (p-value = 2.4×10^-5), but not with body mass index, indicating an association with central fat distribution. To our knowledge, this is the largest population analysis to date of mtDNA dynamics, revealing the age-imposed increase in heteroplasmy, the relatively high heritability of copy number, and the association of copy number with metabolic traits.     

Fidelity of capture-enrichment for mtDNA genome sequencing: influence of NUMTs

Enriching target sequences in sequencing libraries via capture hybridization to bait/probes is an efficient means of leveraging the capabilities of next-generation sequencing for obtaining sequence data from target regions of interest. However, homologous sequences from non-target regions may also be enriched by such methods. Here we investigate the fidelity of capture enrichment for complete mitochondrial DNA (mtDNA) genome sequencing by analyzing sequence data for nuclear copies of mtDNA (NUMTs). Using capture-enriched sequencing data from a mitochondria-free cell line and the parental cell line, and from samples previously sequenced from long-range PCR products, we demonstrate that NUMT alleles are indeed present in capture-enriched sequence data, but at low enough levels to not influence calling the authentic mtDNA genome sequence. However, distinguishing NUMT alleles from true low-level mutations (e.g. heteroplasmy) is more challenging. We develop here a computational method to distinguish NUMT alleles from heteroplasmies, using sequence data from artificial mixtures to optimize the method.     

A sensitive denaturing gradient-Gel electrophoresis assay reveals a high frequency of heteroplasmy in hypervariable region 1 of the human mtDNA control region

A population study of heteroplasmy in the hypervariable region 1 (HV1) portion of the human mtDNA control region was performed. Blood samples from 253 randomly chosen individuals were examined using a sensitive denaturing gradient-gel electrophoresis (DGGE) system. This method is capable of detecting heteroplasmic proportions as low as 1% and virtually all heteroplasmy where the minor component is > or = 5%. Heteroplasmy was observed in 35 individuals (13.8%; 95% confidence interval [CI] 9.6-18.0). Of these individuals, 33 were heteroplasmic at one nucleotide position, whereas 2 were heteroplasmic at two different positions (a condition known as "triplasmy"). Although heteroplasmy occurred at a total of 16 different positions throughout HV1, it was most frequently observed at positions 16093 (n=13) and 16129 (n=6). In addition, the majority of heteroplasmic variants occurred at low proportions and could not be detected by direct sequencing of PCR products. This study indicates that low-level heteroplasmy in HV1 is relatively common and that it occurs at a broad spectrum of sites. Our results corroborate those of other recent reports indicating that heteroplasmy in the control region is more common than was previously believed-a finding that is of potential importance to evolutionary studies and forensic applications that are based on mtDNA variation.     

Analysis of heteroplasmy in the major noncoding region of mitochondrial DNA in the blood and atherosclerotic plaques of carotid arteries

For identification of somatic mitochondrial DNA (mtDNA) mutations, the mtDNA major noncoding region (D-loop) sequence in blood samples and carotid atherosclerosis plaques from patients with atherosclerosis was analyzed. Five point heteroplasmic positions were observed in 4 of 23 individuals (17%). Only in two cases could heteroplasmy have resulted from somatic mutation, whereas three heteroplasmic positions were found in both vascular tissue and blood. In addition, length heteroplasmy in a polycytosine stretches was registered at nucleotide positions 303–315 in 16 individuals, and also in the 16184–16193 region in four patients. The results suggest that somatic mtDNA mutations can occur during atherosclerosis, but some heteroplasmic mutations may appear in all tissues, possibly being inherited.     

Heteroplasmic mitochondrial DNA variants in cardiovascular diseases

Mitochondria are implicated in the pathogenesis of cardiovascular diseases (CVDs) but the reasons for this are not well understood. Maternally-inherited population variants of mitochondrial DNA (mtDNA) which affect all mtDNA molecules (homoplasmic) are associated with cardiometabolic traits and the risk of developing cardiovascular disease. However, it is not known whether mtDNA mutations only affecting a proportion of mtDNA molecules (heteroplasmic) also play a role. To address this question, we performed a high-depth (~1000-fold) mtDNA sequencing of blood DNA in 1,399 individuals with hypertension (HTN), 1,946 with ischemic heart disease (IHD), 2,146 with ischemic stroke (IS), and 723 healthy controls. We show that the per individual burden of heteroplasmic single nucleotide variants (mtSNVs) increases with age. The age-effect was stronger for low-level heteroplasmies (heteroplasmic fraction, HF, 5–10%), likely reflecting acquired somatic events based on trinucleotide mutational signatures. After correcting for age and other confounders, intermediate heteroplasmies (HF 10–95%) were more common in hypertension, particularly involving non-synonymous variants altering the amino acid sequence of essential respiratory chain proteins. These findings raise the possibility that heteroplasmic mtSNVs play a role in the pathophysiology of hypertension.     

Restriction fragment analysis as a source of error in detection of heteroplasmic mtDNA mutations.

The transition from A to G at nt 5656 (5656A-->G) in mitochondrial DNA has been suggested to be a pathogenic mutation and, furthermore, a heteroplasmic one. We found that the mutation was present in 14 out of 83 healthy controls from northern Finland and that 5656A-->G was exclusively associated with mtDNA haplogroup U. Interestingly, 5656A-->G appeared to be heteroplasmic in NheI digestion of PCR fragments that were amplified by using a mismatched oligonucleotide primer creating a digestion site in the presence of the mutant variant. However, we did not detect the wild type genome in clones from such a sample and subsequent experiments revealed that the apparent heteroplasmy was due to inhibition of NheI by NaCl. Our results suggest that 5656A-->G is a polymorphism and it may be highly characteristic for Finns. Furthermore, new heteroplasmic mutations identified by restriction fragment analysis should be adequately controlled for any false positive results that may be due to incomplete digestion.     

Ultra-deep next generation mitochondrial genome sequencing reveals widespread heteroplasmy in Chinese hamster ovary cells

Recent sequencing of the Chinese hamster ovary (CHO) cell and Chinese hamster genomes has dramatically advanced our ability to understand the biology of these mammalian cell factories. In this study, we focus on the powerhouse of the CHO cell, the mitochondrion. Utilizing a high-resolution next generation sequencing approach we sequenced the Chinese hamster mitochondrial genome for the first time and surveyed the mutational landscape of CHO cell mitochondrial DNA (mtDNA). Depths of coverage ranging from ~3,319X to 8,056X enabled accurate identification of low frequency mutations (>1%), revealing that mtDNA heteroplasmy is widespread in CHO cells. A total of 197 variants at 130 individual nucleotide positions were identified across a panel of 22 cell lines with 81% of variants occurring at an allele frequency of between 1% and 99%. 89% of the heteroplasmic mutations identified were cell line specific with the majority of shared heteroplasmic SNPs and INDELs detected in clones from 2 cell line development projects originating from the same host cell line. The frequency of common predicted loss of function mutations varied significantly amongst the clones indicating that heteroplasmic mtDNA variation could lead to a continuous range of phenotypes and play a role in cell to cell, production run to production run and indeed clone to clone variation in CHO cell metabolism. Experiments that integrate mtDNA sequencing with metabolic flux analysis and metabolomics have the potential to improve cell line selection and enhance CHO cell metabolic phenotypes for biopharmaceutical manufacturing through rational mitochondrial genome engineering.     

Restriction fragment analysis as a source of error in detection of heteroplasmic mtDNA mutations

The transition from A to G at nt 5656 (5656A→G) in mitochondrial DNA has been suggested to be a pathogenic mutation and, furthermore, a heteroplasmic one. We found that the mutation was present in 14 out of 83 healthy controls from northern Finland and that 5656A→G was exclusively associated with mtDNA haplogroup U. Interestingly, 5656A→G appeared to be heteroplasmic in NheI digestion of PCR fragments that were amplified by using a mismatched oligonucleotide primer creating a digestion site in the presence of the mutant variant. However, we did not detect the wild type genome in clones from such a sample and subsequent experiments revealed that the apparent heteroplasmy was due to inhibition of NheI by NaCl. Our results suggest that 5656A→G is a polymorphism and it may be highly characteristic for Finns. Furthermore, new heteroplasmic mutations identified by restriction fragment analysis should be adequately controlled for any false positive results that may be due to incomplete digestion.     

Frequency and Pattern of Heteroplasmy in the Complete Human Mitochondrial Genome

Determining the levels of human mitochondrial heteroplasmy is of utmost importance in several fields. In spite of this, there are currently few published works that have focused on this issue. In order to increase the knowledge of mitochondrial DNA (mtDNA) heteroplasmy, the main goal of this work is to investigate the frequency and the mutational spectrum of heteroplasmy in the human mtDNA genome. To address this, a set of nine primer pairs designed to avoid co-amplification of nuclear DNA (nDNA) sequences of mitochondrial origin (NUMTs) was used to amplify the mitochondrial genome in 101 individuals. The analysed individuals represent a collection with a balanced representation of genders and mtDNA haplogroup distribution, similar to that of a Western European population. The results show that the frequency of heteroplasmic individuals exceeds 61%. The frequency of point heteroplasmy is 28.7%, with a widespread distribution across the entire mtDNA. In addition, an excess of transitions in heteroplasmy were detected, suggesting that genetic drift and/or selection may be acting to reduce its frequency at population level. In fact, heteroplasmy at highly stable positions might have a greater impact on the viability of mitochondria, suggesting that purifying selection must be operating to prevent their fixation within individuals. This study analyses the frequency of heteroplasmy in a healthy population, carrying out an evolutionary analysis of the detected changes and providing a new perspective with important consequences in medical, evolutionary and forensic fields.     

Human aging and somatic point mutations in mtDNA: A comparative study of generational differences (grandparents and grandchildren)

The accumulation of somatic mutations in mtDNA is correlated with aging. In this work, we sought to identify somatic mutations in the HVS-1 region (D-loop) of mtDNA that might be associated with aging. For this, we compared 31 grandmothers (mean age: 63 ± 2.3 years) and their 62 grandchildren (mean age: 15 ± 4.1 years), the offspring of their daughters. Direct DNA sequencing showed that mutations absent in the grandchildren were detected in a presumably homoplasmic state in three grandmothers and in a heteroplasmic state in an additional 13 grandmothers; no mutations were detected in the remaining 15 grandmothers. However, cloning followed by DNA sequencing in 12 grandmothers confirmed homoplasia in only one of the three mutations previously considered to be homoplasmic and did not confirm heteroplasmy in three out of nine grandmothers found to be heteroplasmic by direct sequencing. Thus, of 12 grandmothers in whom mtDNA was analyzed by cloning, eight were heteroplasmic for mutations not detected in their grandchildren. In this study, the use of genetically related subjects allowed us to demonstrate the occurrence of age-related (> 60 years old) mutations (homoplasia and heteroplasmy). It is possible that both of these situations (homoplasia and heteroplasmy) were a long-term consequence of mitochondrial oxidative phosphorylation that can lead to the accumulation of mtDNA mutations throughout life.     

Detection of mutations in mitochondrial DNA by droplet digital PCR

Mutations in mitochondrial DNA (mtDNA) may result in various pathological processes. Detection of mutant mtDNAs is a problem for diagnostic practice that is complicated by heteroplasmy – a phenomenon of the inferring presence of at least two allelic variants of the mitochondrial genome. Also, the level of heteroplasmy largely determines the profile and severity of clinical manifestations. Here we discuss detection of mutations in heteroplasmic mtDNA using up-todate methods that have not yet been introduced as routine clinical assays. These methods can be used for detecting mutations in mtDNA to verify diagnosis of “mitochondrial disease”, studying dynamics of mutant mtDNA in body tissues of patients, as well as investigating structural features of mtDNAs. Original data on allele-specific discrimination of m.11778G>A mutation by droplet digital PCR are presented, which demonstrate an opportunity for simultaneous detection and quantitative assessment of mutations in mtDNAs.     

Length and restriction site heteroplasmy in the mitochondrial DNA of american shad (alosa sapidissima)

Restriction endonuclease analysis was used to assess mitochondrial DNA (mtDNA) variation in American shad (Alosa sapidissima) collected from 14 rivers ranging from Florida to Quebec. Two types of heteroplasmy were observed, one involving a major length polymorphism and the other a single restriction site. Shad mtDNA occurred in two principal size classes, 18.3 and 19.8 kb. Of 244 shad examined, 30 were heteroplasmic and carried both size classes of mtDNA in varying proportions; the remainder were homoplasmic for the smaller size class of mtDNA. The large mtDNA variant occurred most frequently at the southern end of the range, and except for two individuals from Nova Scotia, was not detected among shad from rivers north of the Delaware. In contrast, ten shad heteroplasmic for a SalI restriction site originated from rivers ranging from South Carolina to Nova Scotia. DNA mapping and hybridization experiments indicated that the length polymorphism is in the D-loop-containing region and consists of a tandemly repeated 1.5-kb DNA sequence occurring in two and three copies, respectively, in the two major size classes of shad mtDNA. Continuous length variation up to approximately 40 bp occurs among copies of the repeat both within and among individuals. Restriction site data support the conclusion that both forms of heteroplasmy in shad mtDNA have originated more than once.     

Mitochondrial DNA somatic mutation burden and heteroplasmy are associated with chronological age,smoking, and HIV infection

The gradual accumulation of mitochondrial DNA (mtDNA) mutations is implicated in aging and may contribute to the accelerated aging phenotype seen with tobacco smoking and HIV infection. mtDNA mutations are thought to arise from oxidative damage; however, recent reports implicate polymerase γ errors during mtDNA replication. Investigations of somatic mtDNA mutations have been hampered by technical challenges in measuring low‐frequency mutations. We use primer ID‐based next‐generation sequencing to quantify both somatic and heteroplasmic blood mtDNA point mutations within the D‐loop, in 164 women and girls aged 2–72 years, of whom 35% were smokers and 56% were HIV‐positive. Somatic mutations and the occurrence of heteroplasmic mutations increased with age. While transitions are theorized to result from polymerase γ errors, transversions are believed to arise from DNA oxidative damage. In our study, both transition and transversion mutations were associated with age. However, transition somatic mutations were more prevalent than transversions, and no heteroplasmic transversions were observed. We also measured elevated somatic mutations, but not heteroplasmy, in association with high peak HIV viremia. Conversely, heteroplasmy was higher among smokers, but somatic mutations were not, suggesting that smoking promotes the expansion of preexisting mutations rather than de novo mutations. Taken together, our results are consistent with blood mtDNA mutations increasing with age, inferring a greater contribution of polymerase γ errors in mtDNA mutagenesis. We further suggest that smoking and HIV infection both contribute to the accumulation of mtDNA mutations, though in different ways.     

Real-Time PCR Quantification of Heteroplasmy in a Mouse Model with Mitochondrial DNA of C57BL/6 and NZB/BINJ Strains

Mouse models are widely employed to study mitochondrial inheritance, which have implications to several human diseases caused by mutations in the mitochondrial genome (mtDNA). These mouse models take advantage of polymorphisms between the mtDNA of the NZB/BINJ and the mtDNA of common inbred laboratory (i.e., C57BL/6) strains to generate mice with two mtDNA haplotypes (heteroplasmy). Based on PCR followed by restriction fragment length polymorphism (PCR-RFLP), these studies determine the level of heteroplasmy across generations and in different cell types aiming to understand the mechanisms underlying mitochondrial inheritance. However, PCR-RFLP is a time-consuming method of low sensitivity and accuracy that dependents on the use of restriction enzyme digestions. A more robust method to measure heteroplasmy has been provided by the use of real-time quantitative PCR (qPCR) based on allelic refractory mutation detection system (ARMS-qPCR). Herein, we report an ARMS-qPCR assay for quantification of heteroplasmy using heteroplasmic mice with mtDNA of NZB/BINJ and C57BL/6 origin. Heteroplasmy and mtDNA copy number were estimated in germline and somatic tissues, providing evidence of the reliability of the approach. Furthermore, it enabled single-step quantification of heteroplasmy, with sensitivity to detect as low as 0.1% of either NZB/BINJ or C57BL/6 mtDNA. These findings are relevant as the ARMS-qPCR assay reported here is fully compatible with similar heteroplasmic mouse models used to study mitochondrial inheritance in mammals.     

Recurrent Tissue-Specific mtDNA Mutations Are Common in Humans

Mitochondrial DNA (mtDNA) variation can affect phenotypic variation; therefore, knowing its distribution within and among individuals is of importance to understanding many human diseases. Intra-individual mtDNA variation (heteroplasmy) has been generally assumed to be random. We used massively parallel sequencing to assess heteroplasmy across ten tissues and demonstrate that in unrelated individuals there are tissue-specific, recurrent mutations. Certain tissues, notably kidney, liver and skeletal muscle, displayed the identical recurrent mutations that were undetectable in other tissues in the same individuals. Using RFLP analyses we validated one of the tissue-specific mutations in the two sequenced individuals and replicated the patterns in two additional individuals. These recurrent mutations all occur within or in very close proximity to sites that regulate mtDNA replication, strongly implying that these variations alter the replication dynamics of the mutated mtDNA genome. These recurrent variants are all independent of each other and do not occur in the mtDNA coding regions. The most parsimonious explanation of the data is that these frequently repeated mutations experience tissue-specific positive selection, probably through replication advantage.     

Mitochondrial DNA Length Variation and Heteroplasmy in Populations of White Sturgeon (Acipenser Transmontanus)

Southern blot analysis was used to quantify the extent of mtDNA D-loop length variation in two populations of white sturgeon, Acipenser transmontanus. Over 42% of individuals were heteroplasmic for up to six different mtDNA length variants attributable to varying copy numbers of an 82-bp repeat sequence. Chi-square analyses revealed that the frequencies of length genotypes and the incidence of heteroplasmy were significantly different between Fraser and Columbia River sturgeon populations but not between restriction site haplotypes. Heteroplasmic fish have, on average, higher copy number than homoplasmic fish. Forty-five of 101 homoplasmic individuals carry only a single copy of the repeat, while none of the 73 heteroplasmic fish has the single repeat as the predominant variant. On the basis of differences in frequency distributions of copy number within and between fish, we suggest that (1) heteroplasmy is maintained by high recurrent mutation of multiple copy genomes, favoring increased copy number and (2) the mutation pressure toward higher copy number heteroplasmy is partially offset by selection to reduced genome size and segregation to the homoplasmic condition.     

Dynamics of mitochondrial heteroplasmy in three families investigated via a repeatable re-sequencing study

Background

Originally believed to be a rare phenomenon, heteroplasmy - the presence of more than one mitochondrial DNA (mtDNA) variant within a cell, tissue, or individual - is emerging as an important component of eukaryotic genetic diversity. Heteroplasmies can be used as genetic markers in applications ranging from forensics to cancer diagnostics. Yet the frequency of heteroplasmic alleles may vary from generation to generation due to the bottleneck occurring during oogenesis. Therefore, to understand the alterations in allele frequencies at heteroplasmic sites, it is of critical importance to investigate the dynamics of maternal mtDNA transmission.

Results

Here we sequenced, at high coverage, mtDNA from blood and buccal tissues of nine individuals from three families with a total of six maternal transmission events. Using simulations and re-sequencing of clonal DNA, we devised a set of criteria for detecting polymorphic sites in heterogeneous genetic samples that is resistant to the noise originating from massively parallel sequencing technologies. Application of these criteria to nine human mtDNA samples revealed four heteroplasmic sites.

Conclusions

Our results suggest that the incidence of heteroplasmy may be lower than estimated in some other recent re-sequencing studies, and that mtDNA allelic frequencies differ significantly both between tissues of the same individual and between a mother and her offspring. We designed our study in such a way that the complete analysis described here can be repeated by anyone either at our site or directly on the Amazon Cloud. Our computational pipeline can be easily modified to accommodate other applications, such as viral re-sequencing.