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.     

Mutation analysis of the entire mitochondrial genome using denaturing high performance liquid chromatography

In patients with mitochondrial disease a continuously increasing number of mitochondrial DNA (mtDNA) mutations and polymorphisms have been identified. Most pathogenic mtDNA mutations are heteroplasmic, resulting in heteroduplexes after PCR amplification of mtDNA. To detect these heteroduplexes, we used the technique of denaturing high performance liquid chromatography (DHPLC). The complete mitochondrial genome was amplified in 13 fragments of 1–2 kb, digested in fragments of 90–600 bp and resolved at their optimal melting temperature. The sensitivity of the DHPLC system was high with a lowest detection of 0.5% for the A8344G mutation. The muscle mtDNA from six patients with mitochondrial disease was screened and three mutations were identified. The first patient with a limb-girdle-type myopathy carried an A3302G substitution in the tRNA^Leu(UUR) gene (70% heteroplasmy), the second patient with mitochondrial myopathy and cardiomyopathy carried a T3271C mutation in the tRNA^Leu(UUR) gene (80% heteroplasmy) and the third patient with Leigh syndrome carried a T9176C mutation in the ATPase6 gene (93% heteroplasmy). We conclude that DHPLC analysis is a sensitive and specific method to detect heteroplasmic mtDNA mutations. The entire automatic procedure can be completed within 2 days and can also be applied to exclude mtDNA involvement, providing a basis for subsequent investigation of nuclear genes.     

A twin study of mitochondrial DNA polymorphisms shows that heteroplasmy at multiple sites is associated with mtDNA variant 16093 but not with zygosity

The mitochondrial theory of ageing proposes that damage to mitochondria and diminished mitochondrial DNA (mtDNA) repair are major contributors to cellular dysfunction and age-related diseases. We investigate the prevalence of heteroplasmy in the mtDNA control region in buccal swab and blood derived samples for 178 women from the TwinsUK cohort (41 DZ pair 39 MZ pairs, 18 singletons, mean age 57.5 range 28–82) and its relationship to age, BMI and fasting insulin and glucose serum levels. The overall estimated prevalence of heteroplasmy for both tissues in the control region measured for 37 sites was 17%. The prevalence of heteroplasmy was higher among the older half of the study subjects than in the younger half (23% vs 10% p<0.03), primarily reflecting the increase in the prevalence of a heteroplasmic dinucleotide CA repeat in variable region II (VRII) with age. The VRII 523–524 heteroplasmic site (heteroplasmic in 25 subjects) was also associated with a decrease in BMI. In addition, concordance rates for common heteroplasmy were observed to be near complete for both dizygotic (DZ = 94%) and monozygotic twin pairs (MZ = 100%), consistent with previous reports that suggest variation in heteroplasmy rates between generations are determined by bottlenecks in maternal transmission of mitochondria. Differences in the prevalence of heteroplasmy were observed overall between samples derived from buccal swabs (19%) and blood (15%, p<0.04). These were particularly marked at position 16093 of hypervariable region I (HVI, 7% vs 0%, respectively, p<4×10⁻¹¹). The presence of the C allele at position 16093 in blood was associated with the presence of heteroplasmy in buccal swabs at this position (p = 3.5×10⁻¹⁴) and also at VRII (p = 2×10⁻⁴) suggesting a possible predisposing role for this site in the accumulation of heteroplasmy. Our data indicate that BMI is potentially associated with control region heteroplasmy.     

Next-Generation Sequencing of Human Mitochondrial Reference Genomes Uncovers High Heteroplasmy Frequency

We describe methods for rapid sequencing of the entire human mitochondrial genome (mtgenome), which involve long-range PCR for specific amplification of the mtgenome, pyrosequencing, quantitative mapping of sequence reads to identify sequence variants and heteroplasmy, as well as de novo sequence assembly. These methods have been used to study 40 publicly available HapMap samples of European (CEU) and African (YRI) ancestry to demonstrate a sequencing error rate <5.63×10⁻⁴, nucleotide diversity of 1.6×10⁻³ for CEU and 3.7×10⁻³ for YRI, patterns of sequence variation consistent with earlier studies, but a higher rate of heteroplasmy varying between 10% and 50%. These results demonstrate that next-generation sequencing technologies allow interrogation of the mitochondrial genome in greater depth than previously possible which may be of value in biology and medicine.     

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.     

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.     

Cloning of neuronal mtDNA variants in cultured cells by synaptosome fusion with mtDNA-less cells

Synaptosome cybrids were used to confirm the presence of heteroplasmic mtDNA sequence variants in the human brain. Synaptosomes contain one to several mitochondria, and when fused to mtDNA-deficient (ρ°) mouse or human cell lines result in viable cybrid cell lines. The brain origin of mouse synaptosome cybrid mtDNAs was confirmed using sequence polymorphisms in the mtDNA COIII, ND3 and tRNA^Arg genes. The brain origin of the human synaptosome cybrids was confirmed using a rare mtDNA MboI polymorphism. Fusion of synaptosomes from the brain of a 35-year-old woman resulted in 71 synaptosome cybrids. Sequencing the mtDNA control region of these cybrid clones revealed differences in the number of Cs in a poly C track between nucleotide pairs (nps) 301 and 309. Three percent of the cybrid clones had mtDNAs with 10 Cs, 76% had nine, 18% had eight and 3% had seven Cs. Comparable results were obtained by PCR amplification, cloning and sequencing of mtDNA control regions directly from the patient’s brain tissue, but not when the control region was amplified and cloned from a synaptosome cybrid homoplasmic for a mtDNA with nine Cs. Thus, we have clonally recovered mtDNA control region length variants from an adult human brain without recourse to PCR, and established the variant mtDNAs within living cultured cells. This confirms that some mtDNA heteroplasmy can exist in human neurons, and provides the opportunity to study its functional significance.     

Somatic mtDNA Mutation Spectra in the Aging Human Putamen

The accumulation of heteroplasmic mitochondrial DNA (mtDNA) deletions and single nucleotide variants (SNVs) is a well-accepted facet of the biology of aging, yet comprehensive mutation spectra have not been described. To address this, we have used next generation sequencing of mtDNA-enriched libraries (Mito-Seq) to investigate mtDNA mutation spectra of putamen from young and aged donors. Frequencies of the “common” deletion and other “major arc” deletions were significantly increased in the aged cohort with the fold increase in the frequency of the common deletion exceeding that of major arc deletions. SNVs also increased with age with the highest rate of accumulation in the non-coding control region which contains elements necessary for translation and replication. Examination of predicted amino acid changes revealed a skew towards pathogenic SNVs in the coding region driven by mutation bias. Levels of the pathogenic m.3243A>G tRNA mutation were also found to increase with age. Novel multimeric tandem duplications that resemble murine control region multimers and yeast ρ⁻ mtDNAs, were identified in both young and aged specimens. Clonal ∼50 bp deletions in the control region were found at high frequencies in aged specimens. Our results reveal the complex manner in which the mitochondrial genome alters with age and provides a foundation for studies of other tissues and disease states.     

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.     

Comparison of the Illumina Genome Analyzer and Roche 454 GS FLX for Resequencing of Hypertrophic Cardiomyopathy-Associated Genes

Next-generation sequencing (NGS) is widely used in biomedical research, but its adoption has been limited in molecular diagnostics. One application of NGS is the targeted resequencing of genes whose mutations lead to an overlapping clinical phenotype. This study evaluated the comparative performance of the Illumina Genome Analyzer and Roche 454 GS FLX for the resequencing of 16 genes associated with hypertrophic cardiomyopathy (HCM). Using a single human genomic DNA sample enriched by long-range PCR (LR-PCR), 40 GS FLX and 31 Genome Analyzer exon variants were identified using ≥30-fold read-coverage and ≥20% read-percentage selection criteria. Twenty-seven platform concordant variants were Sanger-confirmed. The discordant variants segregated into two categories: variants with read coverages ≥30 on one platform but <30-fold on the alternate platform and variants with read percentages ≥20% on one platform but <20% on the alternate platform. All variants with <30-fold coverage were Sanger-confirmed, suggesting that the coverage criterion of ≥30-fold is too stringent for variant discovery. The variants with <20% read percentage were identified as reference sequence based on Sanger sequencing. These variants were found in homopolymer tracts and short-read misalignments, specifically in genes with high identity. The results of the current study demonstrate the feasibility of combining LR-PCR with the Genome Analyzer or GS FLX for targeted resequencing of HCM-associated genes.     

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.     

Non-Random mtDNA Segregation Patterns Indicate a Metastable Heteroplasmic Segregation Unit in m.3243A>G Cybrid Cells

Many pathogenic mitochondrial DNA mutations are heteroplasmic, with a mixture of mutated and wild-type mtDNA present within individual cells. The severity and extent of the clinical phenotype is largely due to the distribution of mutated molecules between cells in different tissues, but mechanisms underpinning segregation are not fully understood. To facilitate mtDNA segregation studies we developed assays that measure m.3243A>G point mutation loads directly in hundreds of individual cells to determine the mechanisms of segregation over time. In the first study of this size, we observed a number of discrete shifts in cellular heteroplasmy between periods of stable heteroplasmy. The observed patterns could not be parsimoniously explained by random mitotic drift of individual mtDNAs. Instead, a genetically metastable, heteroplasmic mtDNA segregation unit provides the likely explanation, where stable heteroplasmy is maintained through the faithful replication of segregating units with a fixed wild-type/m.3243A>G mutant ratio, and shifts occur through the temporary disruption and re-organization of the segregation units. While the nature of the physical equivalent of the segregation unit remains uncertain, the factors regulating its organization are of major importance for the pathogenesis of mtDNA diseases.     

Exome Sequencing of a Multigenerational Human Pedigree

Over the next few years, the efficient use of next-generation sequencing (NGS) in human genetics research will depend heavily upon the effective mechanisms for the selective enrichment of genomic regions of interest. Recently, comprehensive exome capture arrays have become available for targeting approximately 33 Mb or ∼180,000 coding exons across the human genome. Selective genomic enrichment of the human exome offers an attractive option for new experimental designs aiming to quickly identify potential disease-associated genetic variants, especially in family-based studies. We have evaluated a 2.1 M feature human exome capture array on eight individuals from a three-generation family pedigree. We were able to cover up to 98% of the targeted bases at a long-read sequence read depth of ≥3, 86% at a read depth of ≥10, and over 50% of all targets were covered with ≥20 reads. We identified up to 14,284 SNPs and small indels per individual exome, with up to 1,679 of these representing putative novel polymorphisms. Applying the conservative genotype calling approach HCDiff, the average rate of detection of a variant allele based on Illumina 1 M BeadChips genotypes was 95.2% at ≥10x sequence. Further, we propose an advantageous genotype calling strategy for low covered targets that empirically determines cut-off thresholds at a given coverage depth based on existing genotype data. Application of this method was able to detect >99% of SNPs covered ≥8x. Our results offer guidance for “real-world” applications in human genetics and provide further evidence that microarray-based exome capture is an efficient and reliable method to enrich for chromosomal regions of interest in next-generation sequencing experiments.     

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.     

Private Mitochondrial DNA Variants in Danish Patients with Hypertrophic Cardiomyopathy

Hypertrophic cardiomyopathy (HCM) is a genetic cardiac disease primarily caused by mutations in genes coding for sarcomeric proteins. A molecular-genetic etiology can be established in ~60% of cases. Evolutionarily conserved mitochondrial DNA (mtDNA) haplogroups are susceptibility factors for HCM. Several polymorphic mtDNA variants are associated with a variety of late-onset degenerative diseases and affect mitochondrial function. We examined the role of private, non-haplogroup associated, mitochondrial variants in the etiology of HCM. In 87 Danish HCM patients, full mtDNA sequencing revealed 446 variants. After elimination of 312 (69.9%) non-coding and synonymous variants, a further 109 (24.4%) with a global prevalence > 0.1%, three (0.7%) haplogroup associated and 19 (2.0%) variants with a low predicted in silico likelihood of pathogenicity, three variants: MT-TC: m.5772G>A, MT-TF: m.644A>G, and MT-CYB: m.15024G>A, p.C93Y remained. A detailed analysis of these variants indicated that none of them are likely to cause HCM. In conclusion, private mtDNA mutations are frequent, but they are rarely, if ever, associated with HCM.     

A Wide Range of 3243A>G/tRNALeu(UUR) (MELAS) Mutation Loads May Segregate in Offspring through the Female Germline Bottleneck

Segregation of mutant mtDNA in human tissues and through the germline is debated, with no consensus about the nature and size of the bottleneck hypothesized to explain rapid generational shifts in mutant loads. We investigated two maternal lineages with an apparently different inheritance pattern of the same pathogenic mtDNA 3243A>G/tRNA^Leu(UUR) (MELAS) mutation. We collected blood cells, muscle biopsies, urinary epithelium and hair follicles from 20 individuals, as well as oocytes and an ovarian biopsy from one female mutation carrier, all belonging to the two maternal lineages to assess mutant mtDNA load, and calculated the theoretical germline bottleneck size (number of segregating units). We also evaluated “mother-to-offspring” segregations from the literature, for which heteroplasmy assessment was available in at least three siblings besides the proband. Our results showed that mutation load was prevalent in skeletal muscle and urinary epithelium, whereas in blood cells there was an inverse correlation with age, as previously reported. The histoenzymatic staining of the ovarian biopsy failed to show any cytochrome-c-oxidase defective oocyte. Analysis of four oocytes and one offspring from the same unaffected mother of the first family showed intermediate heteroplasmic mutant loads (10% to 75%), whereas very skewed loads of mutant mtDNA (0% or 81%) were detected in five offspring of another unaffected mother from the second family. Bottleneck size was 89 segregating units for the first mother and 84 for the second. This was remarkably close to 88, the number of “segregating units” in the “mother-to-offspring” segregations retrieved from literature. In conclusion, a wide range of mutant loads may be found in offspring tissues and oocytes, resulting from a similar theoretical bottleneck size.     

Next‐generation DNA barcoding: using next‐generation sequencing to enhance and accelerate DNA barcode capture from single specimens

DNA barcoding is an efficient method to identify specimens and to detect undescribed/cryptic species. Sanger sequencing of individual specimens is the standard approach in generating large‐scale DNA barcode libraries and identifying unknowns. However, the Sanger sequencing technology is, in some respects, inferior to next‐generation sequencers, which are capable of producing millions of sequence reads simultaneously. Additionally, direct Sanger sequencing of DNA barcode amplicons, as practiced in most DNA barcoding procedures, is hampered by the need for relatively high‐target amplicon yield, coamplification of nuclear mitochondrial pseudogenes, confusion with sequences from intracellular endosymbiotic bacteria (e.g. Wolbachia) and instances of intraindividual variability (i.e. heteroplasmy). Any of these situations can lead to failed Sanger sequencing attempts or ambiguity of the generated DNA barcodes. Here, we demonstrate the potential application of next‐generation sequencing platforms for parallel acquisition of DNA barcode sequences from hundreds of specimens simultaneously. To facilitate retrieval of sequences obtained from individual specimens, we tag individual specimens during PCR amplification using unique 10‐mer oligonucleotides attached to DNA barcoding PCR primers. We employ 454 pyrosequencing to recover full‐length DNA barcodes of 190 specimens using 12.5% capacity of a 454 sequencing run (i.e. two lanes of a 16 lane run). We obtained an average of 143 sequence reads for each individual specimen. The sequences produced are full‐length DNA barcodes for all but one of the included specimens. In a subset of samples, we also detected Wolbachia, nontarget species, and heteroplasmic sequences. Next‐generation sequencing is of great value because of its protocol simplicity, greatly reduced cost per barcode read, faster throughout and added information content.     

Skewed Segregation of the mtDNA nt 8993 (TrG) Mutation in Human Oocytes

Rapid changes in mtDNA variants between generations have led to the bottleneck theory, which proposes a dramatic reduction in mtDNA numbers during early oogenesis. We studied oocytes from a woman with heteroplasmic expression of the mtDNA nt 8993 (T-->G) mutation. Of seven oocytes analyzed, one showed no evidence of the mutation, and the remaining six had a mutant load > 95%. This skewed expression of the mutation in oocytes is not compatible with the conventional bottleneck theory. A possible explanation is that, during amplification of mtDNA in the developing oocyte, mtDNA from one mitochondrion is preferentially amplified. Thus, subsequent mature oocytes may contain predominantly wild-type or mutant mitochondrial genomes.