Search

© 2017 Elsevier. This manuscript version is made available under the CC-BY-NC-ND 4.0 license http://creativecommons.org/licenses/by-nc-nd/4.0/ This author accepted manuscript is made available following 12 month embargo from date of publication (July 2017) in accordance with the publisher’s archiving policy, 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.

We report on a method to analyze length heteroplasmy within the human mitochondrial genome in which there are polycytosine [poly(C)] stretches. These poly(C) tracts induce heteroplasmy with the resultant inherent problems of accurate sequence designations. In this study, 20 samples that exhibited length heteroplasmy due to variation in the C-tracts within hypervariable region I (HVI) were treated with bisulfite, and one or more cytosine bases in these C-tracts were converted randomly to uracil. This resulted in an accurate sequence designation for nearly all samples. The only exceptions in which the DNA sequence could still not be determined occurred when there was total conversion, or a lack of conversion, of the cytosine bases. Replicate tests on the same samples showed that individual cytosine bases were randomly converted to uracil. This simple method was useful for investigating length heteroplasmy due to 16189C and 310C transitions in the mitochondrial-DNA control region. It is valuable for medical and forensic investigations.