Mitochondrial genome analysis in penile carcinoma

Penile cancer is a rare neoplasm that seems to be linked to socio-economic differences. Mitochondrial genome alterations are common in many tumors types and are reported as regulating oxidative metabolism and impacting tumorigenesis. In this study, we evaluate for the first time the mitochondrial genome in penile carcinoma (PeCa), aiming to evaluate heteroplasmy, mitochondrial DNA (mtDNA) mutational load and mtDNA content in Penile tumors. Using next generation sequencing (NGS), we sequenced the mitochondrial genome of 13 penile tumors and 12 non-neoplastic tissue samples, which allowed us to identify mtDNA variants and heteroplasmy. We further evaluated variant’s pathogenicity using Mutpred predictive software and calculated mtDNA content using quantitative PCR. Mitochondrial genome sequencing revealed an increase number of non-synonymous variants in the tumor tissue, along with higher frequency of heteroplasmy and mtDNA depletion in penile tumors, suggesting an increased mitochondrial instability in penile tumors. We also described a list of mitochondrial variants found in penile tumor and normal tissue, including five novel variants found in the tumoral tissue. Our results showed an increased mitochondrial genome instability in penile tumors. We also suggest that mitochondrial DNA copy number (mtDNAcn) and mtDNA variants may act together to imbalance mitochondrial function in PeCa. The better understanding of mitochondrial biology can bring new insights on mechanisms and open a new field for therapy in PeCa.     

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.     

Mitochondrial DNA Variant Discovery and Evaluation in Human Cardiomyopathies through Next-Generation Sequencing

Mutations in mitochondrial DNA (mtDNA) may cause maternally-inherited cardiomyopathy and heart failure. In homoplasmy all mtDNA copies contain the mutation. In heteroplasmy there is a mixture of normal and mutant copies of mtDNA. The clinical phenotype of an affected individual depends on the type of genetic defect and the ratios of mutant and normal mtDNA in affected tissues. We aimed at determining the sensitivity of next-generation sequencing compared to Sanger sequencing for mutation detection in patients with mitochondrial cardiomyopathy. We studied 18 patients with mitochondrial cardiomyopathy and two with suspected mitochondrial disease. We “shotgun” sequenced PCR-amplified mtDNA and multiplexed using a single run on Roche''s 454 Genome Sequencer. By mapping to the reference sequence, we obtained 1,300× average coverage per case and identified high-confidence variants. By comparing these to >400 mtDNA substitution variants detected by Sanger, we found 98% concordance in variant detection. Simulation studies showed that >95% of the homoplasmic variants were detected at a minimum sequence coverage of 20× while heteroplasmic variants required >200× coverage. Several Sanger “misses” were detected by 454 sequencing. These included the novel heteroplasmic 7501T>C in tRNA serine 1 in a patient with sudden cardiac death. These results support a potential role of next-generation sequencing in the discovery of novel mtDNA variants with heteroplasmy below the level reliably detected with Sanger sequencing. We hope that this will assist in the identification of mtDNA mutations and key genetic determinants for cardiomyopathy and mitochondrial disease.     

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.     

Simultaneous A8344G heteroplasmy and mitochondrial DNA copy number quantification in myoclonus epilepsy and ragged-red fibers (MERRF) syndrome by a multiplex molecular beacon based real-time fluorescence PCR

The association of a particular mitochondrial DNA (mtDNA) mutation with different clinical phenotypes is a well-known feature of mitochondrial diseases. A simple genotype–phenotype correlation has not been found between mutation load and disease expression. Tissue and intercellular mosaicism as well as mtDNA copy number are thought to be responsible for the different clinical phenotypes. As disease expression of mitochondrial tRNA mutations is mostly in postmitotic tissues, studies to elucidate disease mechanisms need to be performed on patient material. Heteroplasmy quantitation and copy number estimation using small patient biopsy samples has not been reported before, mainly due to technical restrictions. In order to resolve this problem, we have developed a robust assay that utilizes Molecular Beacons to accurately quantify heteroplasmy levels and determine mtDNA copy number in small samples carrying the A8344G tRNA^Lys mutation. It provides the methodological basis to investigate the role of heteroplasmy and mtDNA copy number in determining the clinical phenotypes.     

Regional variation in mitochondrial DNA copy number in mouse brain

Mitochondria have their own DNA (mitochondrial DNA [mtDNA]). Although mtDNA copy number is dependent on tissues and its decrease is associated with various neuromuscular diseases, detailed distribution of mtDNA copies in the brain remains uncertain. Using real-time quantitative PCR assay, we examined regional variation in mtDNA copy number in 39 brain regions of male mice. A significant regional difference in mtDNA copy number was observed (P<4.8×10(-35)). High levels of mtDNA copies were found in the ventral tegmental area and substantia nigra, two major nuclei containing dopaminergic neurons. In contrast, cerebellar vermis and lobes had significantly lower copy numbers than other regions. Hippocampal dentate gyrus also had a relatively low mtDNA copy number. This study is the first quantitative analysis of regional variation in mtDNA copy number in mouse brain. Our findings are important for the physiological and pathophysiological studies of mtDNA in the brain.     

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.     

MitoScape: A big-data,machine-learning platform for obtaining mitochondrial DNA from next-generation sequencing data

The growing number of next-generation sequencing (NGS) data presents a unique opportunity to study the combined impact of mitochondrial and nuclear-encoded genetic variation in complex disease. Mitochondrial DNA variants and in particular, heteroplasmic variants, are critical for determining human disease severity. While there are approaches for obtaining mitochondrial DNA variants from NGS data, these software do not account for the unique characteristics of mitochondrial genetics and can be inaccurate even for homoplasmic variants. We introduce MitoScape, a novel, big-data, software for extracting mitochondrial DNA sequences from NGS. MitoScape adopts a novel departure from other algorithms by using machine learning to model the unique characteristics of mitochondrial genetics. We also employ a novel approach of using rho-zero (mitochondrial DNA-depleted) data to model nuclear-encoded mitochondrial sequences. We showed that MitoScape produces accurate heteroplasmy estimates using gold-standard mitochondrial DNA data. We provide a comprehensive comparison of the most common tools for obtaining mtDNA variants from NGS and showed that MitoScape had superior performance to compared tools in every statistically category we compared, including false positives and false negatives. By applying MitoScape to common disease examples, we illustrate how MitoScape facilitates important heteroplasmy-disease association discoveries by expanding upon a reported association between hypertrophic cardiomyopathy and mitochondrial haplogroup T in men (adjusted p-value = 0.003). The improved accuracy of mitochondrial DNA variants produced by MitoScape will be instrumental in diagnosing disease in the context of personalized medicine and clinical diagnostics.     

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.     

Segregation of Naturally Occurring Mitochondrial DNA Variants in a Mini-Pig Model

The maternally inherited mitochondrial genome (mtDNA) is present in multimeric form within cells and harbors sequence variants (heteroplasmy). While a single mtDNA variant at high load can cause disease, naturally occurring variants likely persist at low levels across generations of healthy populations. To determine how naturally occurring variants are segregated and transmitted, we generated a mini-pig model, which originates from the same maternal ancestor. Following next-generation sequencing, we identified a series of low-level mtDNA variants in blood samples from the female founder and her daughters. Four variants, ranging from 3% to 20%, were selected for validation by high-resolution melting analysis in 12 tissues from 31 animals across three generations. All four variants were maintained in the offspring, but variant load fluctuated significantly across the generations in several tissues, with sex-specific differences in heart and liver. Moreover, variant load was persistently reduced in high-respiratory organs (heart, brain, diaphragm, and muscle), which correlated significantly with higher mtDNA copy number. However, oocytes showed increased heterogeneity in variant load, which correlated with increased mtDNA copy number during in vitro maturation. Altogether, these outcomes show that naturally occurring mtDNA variants segregate and are maintained in a tissue-specific manner across generations. This segregation likely involves the maintenance of selective mtDNA variants during organogenesis, which can be differentially regulated in oocytes and preimplantation embryos during maturation.     

Detecting Heteroplasmy from High-Throughput Sequencing of Complete Human Mitochondrial DNA Genomes

Heteroplasmy, the existence of multiple mtDNA types within an individual, has been previously detected by using mostly indirect methods and focusing largely on just the hypervariable segments of the control region. Next-generation sequencing technologies should enable studies of heteroplasmy across the entire mtDNA genome at much higher resolution, because many independent reads are generated for each position. However, the higher error rate associated with these technologies must be taken into consideration to avoid false detection of heteroplasmy. We used simulations and phiX174 sequence data to design criteria for accurate detection of heteroplasmy with the Illumina Genome Analyzer platform, and we used artificial mixtures and replicate data to test and refine the criteria. We then applied these criteria to mtDNA sequence reads for 131 individuals from five Eurasian populations that had been generated via a parallel tagged approach. We identified 37 heteroplasmies at 10% frequency or higher at 34 sites in 32 individuals. The mutational spectrum does not differ between heteroplasmic mutations and polymorphisms in the same individuals, but the relative mutation rate at heteroplasmic mutations is significantly higher than that estimated for all mutable sites in the human mtDNA genome. Moreover, there is also a significant excess of nonsynonymous mutations observed among heteroplasmies, compared to polymorphism data from the same individuals. Both mutation-drift and negative selection influence the fate of heteroplasmies to determine the polymorphism spectrum in humans. With appropriate criteria for avoiding false positives due to sequencing errors, next-generation technologies can provide novel insights into genome-wide aspects of mtDNA heteroplasmy.     

The strength and timing of the mitochondrial bottleneck in salmon suggests a conserved mechanism in vertebrates

In most species mitochondrial DNA (mtDNA) is inherited maternally in an apparently clonal fashion, although how this is achieved remains uncertain. Population genetic studies show not only that individuals can harbor more than one type of mtDNA (heteroplasmy) but that heteroplasmy is common and widespread across a diversity of taxa. Females harboring a mixture of mtDNAs may transmit varying proportions of each mtDNA type (haplotype) to their offspring. However, mtDNA variants are also observed to segregate rapidly between generations despite the high mtDNA copy number in the oocyte, which suggests a genetic bottleneck acts during mtDNA transmission. Understanding the size and timing of this bottleneck is important for interpreting population genetic relationships and for predicting the inheritance of mtDNA based disease, but despite its importance the underlying mechanisms remain unclear. Empirical studies, restricted to mice, have shown that the mtDNA bottleneck could act either at embryogenesis, oogenesis or both. To investigate whether the size and timing of the mitochondrial bottleneck is conserved between distant vertebrates, we measured the genetic variance in mtDNA heteroplasmy at three developmental stages (female, ova and fry) in chinook salmon and applied a new mathematical model to estimate the number of segregating units (N(e)) of the mitochondrial bottleneck between each stage. Using these data we estimate values for mtDNA Ne of 88.3 for oogenesis, and 80.3 for embryogenesis. Our results confirm the presence of a mitochondrial bottleneck in fish, and show that segregation of mtDNA variation is effectively complete by the end of oogenesis. Considering the extensive differences in reproductive physiology between fish and mammals, our results suggest the mechanism underlying the mtDNA bottleneck is conserved in these distant vertebrates both in terms of it magnitude and timing. This finding may lead to improvements in our understanding of mitochondrial disorders and population interpretations using mtDNA data.     

Absolute quantitation of a heteroplasmic mitochondrial DNA deletion using a multiplex three-primer real-time PCR assay

Quantitation of wild-type and deleted mitochondrial DNA (mtDNA) coexisting within the same cell (a.k.a., heteroplasmy) is important in mitochondrial disease and aging. We report the development of a multiplex three-primer PCR assay that is capable of absolute quantitation of wild-type and deleted mtDNA simultaneously. Molecular beacons were designed to hybridize with either type of mtDNA molecule, allowing real-time detection during PCR amplification. The assay is specific and can detect down to six copies of mtDNA, making it suitable for single-cell analyses. The relative standard deviation in the threshold cycle number is approximately 0.6%. Heteroplasmy was quantitated in individual cytoplasmic hybrid cells (cybrids), containing a large mtDNA deletion, and bulk cell samples. Individual cybrid cells contained 100-2600 copies of wild-type mtDNA and 950-4700 copies of deleted mtDNA, and the percentage of heteroplasmy ranged from 43+/-16 to 95+/-16%. The average amount of total mtDNA was 3800+/-1600 copies/cybrid cell, and the average percentage of heteroplasmy correlated well with the bulk cell sample. The single-cell analysis also revealed that heteroplasmy in individual cells is highly heterogeneous. This assay will be useful for monitoring clonal expansions of mtDNA deletions and investigating the role of heteroplasmy in cell-to-cell heterogeneity in cellular models of mitochondrial disease and aging.     

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.     

Very Low-Level Heteroplasmy mtDNA Variations Are Inherited in Humans

Little is known about the inheritance of very low heteroplasmy mitochondria DNA （mtDNA） variations. Even with the development of new next-generation sequencing methods, the practical lower limit of measured heteroplasmy is still about 1% due to the inherent noise level of the sequencing. In this study, we sequenced the mitochondrial genome of 44 individuals using Illumina high-throughput sequencing technology and obtained high-coverage mitochondria sequencing data. Our study population contains many mother-offspring pairs. This unique study design allows us to bypass the usual heteroplasmy limitation by analyzing the correlation of mutation levels at each position in the mtDNA sequence between maternally related pairs and non-related pairs. The study showed that very low heteroplasmy variants, down to almost 0.1%, are inherited maternally and that this inheritance begins to decrease at about 0.5%, cor- resnondin to abottleneck of about 200 mtDNA.     

Presequence-Independent Mitochondrial Import of DNA Ligase Facilitates Establishment of Cell Lines with Reduced mtDNA Copy Number

Due to the essential role played by mitochondrial DNA (mtDNA) in cellular physiology and bioenergetics, methods for establishing cell lines with altered mtDNA content are of considerable interest. Here, we report evidence for the existence in mammalian cells of a novel, low- efficiency, presequence-independent pathway for mitochondrial protein import, which facilitates mitochondrial uptake of such proteins as Chlorella virus ligase (ChVlig) and Escherichia coli LigA. Mouse cells engineered to depend on this pathway for mitochondrial import of the LigA protein for mtDNA maintenance had severely (up to >90%) reduced mtDNA content. These observations were used to establish a method for the generation of mouse cell lines with reduced mtDNA copy number by, first, transducing them with a retrovirus encoding LigA, and then inactivating in these transductants endogenous Lig3 with CRISPR-Cas9. Interestingly, mtDNA depletion to an average level of one copy per cell proceeds faster in cells engineered to maintain mtDNA at low copy number. This makes a low-mtDNA copy number phenotype resulting from dependence on mitochondrial import of DNA ligase through presequence-independent pathway potentially useful for rapidly shifting mtDNA heteroplasmy through partial mtDNA depletion.     

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.     

Length Variation in Mitochondrial DNA of the Minnow Cyprinella Spiloptera

Length differences in animal mitochondrial DNA (mtDNA) are common, frequently due to variation in copy number of direct tandem duplications. While such duplications appear to form without great difficulty in some taxonomic groups, they appear to be relatively short-lived, as typical duplication products are geographically restricted within species and infrequently shared among species. To better understand such length variation, we have studied a tandem and direct duplication of approximately 260 bp in the control region of the cyprinid fish, Cyprinella spiloptera. Restriction site analysis of 38 individuals was used to characterize population structure and the distribution of variation in repeat copy number. This revealed two length variants, including individuals with two or three copies of the repeat, and little geographic structure among populations. No standard length (single copy) genomes were found and heteroplasmy, a common feature of length variation in other taxa, was absent. Nucleotide sequence of tandem duplications and flanking regions localized duplication junctions in the phenylalanine tRNA and near the origin of replication. The locations of these junctions and the stability of folded repeat copies support the hypothesized importance of secondary structures in models of duplication formation.     

Genotypic stability, segregation and selection in heteroplasmic human cell lines containing np 3243 mutant mtDNA

The mitochondrial genotype of heteroplasmic human cell lines containing the pathological np 3243 mtDNA mutation, plus or minus its suppressor at np 12300, has been followed over long periods in culture. Cell lines containing various different proportions of mutant mtDNA remained generally at a consistent, average heteroplasmy value over at least 30 wk of culture in nonselective media and exhibited minimal mitotic segregation, with a segregation number comparable with mtDNA copy number (>/=1000). Growth in selective medium of cells at 99% np 3243 mutant mtDNA did, however, allow the isolation of clones with lower levels of the mutation, against a background of massive cell death. As a rare event, cell lines exhibited a sudden and dramatic diversification of heteroplasmy levels, accompanied by a shift in the average heteroplasmy level over a short period (<8 wk), indicating selection. One such episode was associated with a gain of chromosome 9. Analysis of respiratory phenotype and mitochondrial genotype of cell clones from such cultures revealed that stable heteroplasmy values were generally reestablished within a few weeks, in a reproducible but clone-specific fashion. This occurred independently of any straightforward phenotypic selection at the individual cell-clone level. Our findings are consistent with several alternate views of mtDNA organization in mammalian cells. One model that is supported by our data is that mtDNA is found in nucleoids containing many copies of the genome, which can themselves be heteroplasmic, and which are faithfully replicated. We interpret diversification and shifts of heteroplasmy level as resulting from a reorganization of such nucleoids, under nuclear genetic control. Abrupt remodeling of nucleoids in vivo would have major implications for understanding the developmental consequences of heteroplasmy, including mitochondrial disease phenotype and progression.     

Stable heteroplasmy at the single-cell level is facilitated by intercellular exchange of mtDNA

Eukaryotic cells carry two genomes, nuclear (nDNA) and mitochondrial (mtDNA), which are ostensibly decoupled in their replication, segregation and inheritance. It is increasingly appreciated that heteroplasmy, the occurrence of multiple mtDNA haplotypes in a cell, plays an important biological role, but its features are not well understood. Accurately determining the diversity of mtDNA has been difficult, due to the relatively small amount of mtDNA in each cell (<1% of the total DNA), the intercellular variability of mtDNA content and mtDNA pseudogenes (Numts) in nDNA. To understand the nature of heteroplasmy, we developed Mseek, a novel technique to purify and sequence mtDNA. Mseek yields high purity (>90%) mtDNA and its ability to detect rare variants is limited only by sequencing depth, providing unprecedented sensitivity and specificity. Using Mseek, we confirmed the ubiquity of heteroplasmy by analyzing mtDNA from a diverse set of cell lines and human samples. Applying Mseek to colonies derived from single cells, we find heteroplasmy is stably maintained in individual daughter cells over multiple cell divisions. We hypothesized that the stability of heteroplasmy could be facilitated by intercellular exchange of mtDNA. We explicitly demonstrate this exchange by co-culturing cell lines with distinct mtDNA haplotypes. Our results shed new light on the maintenance of heteroplasmy and provide a novel platform to investigate features of heteroplasmy in normal and diseased states.