Investigation into length heteroplasmy in the mitochondrial DNA control region after treatment with bisulfite

Forensic evaluation of mitochondrial DNA heteroplasmy in Gujarat population, India

BACKGROUND

Owing to its high copy number and its small size, mtDNA analysis is the most reliable choice when biological materials from crime scenes are degraded or have mixed STR profiles.

AIM

To examine the occurrence of heteroplasmy along with its frequency and pattern in both HV1 and HV2 regions of the mtDNA among unrelated individuals from India.

SUBJECTS AND METHODS

Mitochondrial DNA control region [hypervariable region one (HV1) and hypervariable region two (HV2)] were analysed in blood and buccal tissues of 104 unrelated individuals from the Indian state of Gujarat.

RESULTS

A high frequency of point heteroplasmy (PH) and length heteroplasmy (LH) was revealed. PH was detected in 7.69% of the population, with a higher frequency observed in blood than in buccal samples. However, there were no statistically significant differences in PH between the two tissues (Chi-square = 0.552, p ≥ 0.05). A total of six PH positions were detected: three at HV1, and another three at HV2. The studied population showed 46.15% LH in the HV1 and HV2 regions of both tissues. The LH positions observed in the Gujarat population were the same as those previously reported at HV1 np16184-16193 and HV2 np303-315.

CONCLUSIONS

Our findings suggest that differences in the pattern of heteroplasmy found in different tissues can complicate the forensic analysis, on the other hand, the probability of a match between the questioned and reference samples increases when the heteroplasmy is identical in both tissues. Variability of PH among persons and even within tissues recommends analysing multiple tissue samples before drawing a conclusion in forensic mtDNA analyses.

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.