Alterations of length heteroplasmy in mitochondrial DNA under various amplification conditions

Investigation of length heteroplasmy in mitochondrial DNA control region by massively parallel sequencing

Accurate sequencing of the control region of the mitochondrial genome is notoriously difficult due to the presence of polycytosine bases, termed C-tracts. The precise number of bases that constitute a C-tract and the bases beyond the poly cytosines may not be accurately defined when analyzing Sanger sequencing data separated by capillary electrophoresis. Massively parallel sequencing has the potential to resolve such poor definition and provides the opportunity to discover variants due to length heteroplasmy. In this study, the control region of mitochondrial genomes from 20 samples was sequenced using both standard Sanger methods with separation by capillary electrophoresis and also using massively parallel DNA sequencing technology. After comparison of the two sets of generated sequence, with the exception of the C-tracts where length heteroplasmy was observed, all sequences were concordant. Sequences of three segments 16184-16193, 303-315 and 568-573 with C-tracts in HVI, II and III can be clearly defined from the massively parallel sequencing data using the program SEQ Mapper. Multiple sequence variants were observed in the length of C-tracts longer than 7 bases. Our report illustrates the accurate designation of all the length variants leading to heteroplasmy in the control region of the mitochondrial genome that can be determined by SEQ Mapper based on data generated by massively parallel DNA sequencing.

Length heteroplasmy of the polyC-polyT-polyC stretch in the dog mtDNA control region

Previously, the mitochondrial control region of 214 Belgian dogs was sequenced. Analysis of this data indicated length heteroplasmy of the polyT stretch in the polyC-polyT-polyC stretch from positions 16661 to 16674. Nine polyC-polyT-polyC haplotype combinations were observed, consisting of seven major haplotypes (highest signal intensity) combined with minor haplotypes (lower signal intensity) one T shorter than the major haplotype in all but three dogs. The longer the polyT stretch, the smaller was the difference in signal intensity between the major and minor haplotype peaks. Additional sequencing, cloning, and PCR trap experiments were performed to further study the intra-individual variation of this mitochondrial DNA (mtDNA) region. Cloning experiments demonstrated that the proportion of clones displaying the minor haplotypes also increased with the length of the polyT stretch. Clone amplification showed that in vitro polymerase errors might contribute to the length heteroplasmy of polyT stretches with at least 10 Ts. Although major and minor polyC-polyT-polyC haplotypes did not differ intra-individually within and between tissues in this study, interpretation of polyT stretch variation should be handled with care in forensic casework.

The complete mitochondrial genome analysis of the tiger (Panthera tigris)

The complete mitochondrial genomes of five tiger samples from three subspecies (P. t. sumatrae, P. t. altica, and P. t. tigris) were successfully obtained by using 26 specifically designed Panthera-specific primer sets. The genome organization and gene arrangement of the five tiger samples were similar to each other; however polymorphic tandem repeat sequences were observed in the control region (CR). This led to a difference in the genome lengths obtained from these five samples with an average size of 16,994 bp for the five tiger mitochondrial genomes. The nucleotide base composition was on average as follows: A, 31.8%; T, 27.0%; C, 26.6%; G, 14.6% and exhibited compositional asymmetry. Most of tiger mitochondrial genome characteristics are similar to those of other common vertebrate species; however, some distinctive features were observed in the CR. First, the repetitive sequence 2 (RS 2) contained two repeat units of 80 bp and the first 15 bp of what would be the third repeat motif. The repetitive sequence 3 (RS 3) contained 47-50 repeat motifs of a shorter 8 bp (ACGTAYAC)(n). Second, length heteroplasmy polycystosine (poly-C) stretches was observed at the end of the HV I locus in all tiger samples.

High-throughput sequencing of core STR loci for forensic genetic investigations using the Roche Genome Sequencer FLX platform

The analysis and profiling of short tandem repeat (STR) loci is routinely used in forensic genetics. Current methods to investigate STR loci, including PCR-based standard fragment analyses and capillary electrophoresis, only provide amplicon lengths that are used to estimate the number of STR repeat units. These methods do not allow for the full resolution of STR base composition that sequencing approaches could provide. Here we present an STR profiling method based on the use of the Roche Genome Sequencer (GS) FLX to simultaneously sequence multiple core STR loci. Using this method in combination with a bioinformatic tool designed specifically to analyze sequence lengths and frequencies, we found that GS FLX STR sequence data are comparable to conventional capillary electrophoresis-based STR typing. Furthermore, we found DNA base substitutions and repeat sequence variations that would not have been identified using conventional STR typing.

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.