Validation of haplotype-specific extraction for separating a mitochondrial DNA model mixture and application to simulated casework

Whole mitochondrial genome analysis of Tai-Kadai-speaking populations in Southwest China

As a single matrilineal gene, human mitochondrial DNA plays a very important role in the study of population genetics. The whole mitogenome sequences of 287 individuals of the Tai-Kadai-speaking population in Guizhou were obtained. It was discovered that there were 82, 104, and 94 haplotypes in 83 Bouyei individuals, 107 Dong individuals, and 97 Sui individuals, respectively; and the haplotype diversity in Bouyei, Dong, and Sui groups was 1.000 ± 0.02, 0.9993 ± 0.0015, and 0.999 ± 0.002, respectively. The result of neutrality tests of the Tai-Kadai-speaking population in Guizhou showed significant negative values, and the analysis of mismatch distribution showed an obvious unimodal distribution. The results implied that Guizhou Tai-Kadai-speaking populations had high genetic diversities and may have experienced recent population expansion. In addition, the primary haplogroups of studied populations were M*, F, B, D, and R*, implying that they may origin from Southern China. The matrilineal genetic structure of the Tai-Kadai-speaking populations in Guizhou was analyzed by merging the mitogenome data of 79 worldwide populations as reference data. The results showed that there were close relationships between studied populations and other Tai-Kadai as well as some Austronesian populations in East and Southeast Asia. Overall, the mitogenome data generated in this study will provide important data for the study of genetic structure of Tai-Kadai speaking populations.

Analyzing Low-Level mtDNA Heteroplasmy-Pitfalls and Challenges from Bench to Benchmarking

Massive parallel sequencing technologies are promising a highly sensitive detection of low-level mutations, especially in mitochondrial DNA (mtDNA) studies. However, processes from DNA extraction and library construction to bioinformatic analysis include several varying tasks. Further, there is no validated recommendation for the comprehensive procedure. In this study, we examined potential pitfalls on the sequencing results based on two-person mtDNA mixtures. Therefore, we compared three DNA polymerases, six different variant callers in five mixtures between 50% and 0.5% variant allele frequencies generated with two different amplification protocols. In total, 48 samples were sequenced on Illumina MiSeq. Low-level variant calling at the 1% variant level and below was performed by comparing trimming and PCR duplicate removal as well as six different variant callers. The results indicate that sensitivity, specificity, and precision highly depend on the investigated polymerase but also vary based on the analysis tools. Our data highlight the advantage of prior standardization and validation of the individual laboratory setup with a DNA mixture model. Finally, we provide an artificial heteroplasmy benchmark dataset that can help improve somatic variant callers or pipelines, which may be of great interest for research related to cancer and aging.

Whole mitochondrial genome analysis of highland Tibetan ethnicity using massively parallel sequencing

Mitochondrial DNA (mtDNA) is a key player in numerous multifaceted and intricate biological processes and plays a pivotal role in dissecting the peopling of different populations, due to its maternally inherited property and comparatively high mutation rate. In this study, 119 Tibetan individuals from the Muli Tibetan Autonomous County of China (average altitude above 3,000 m) were employed in mitochondrial genome (mitogenome) sequencing by massively parallel sequencing (MPS) techniques using the Precision ID mtDNA Whole Genome Panel on an Ion S5XL system. The dataset presented 88 distinct haplotypes, resulting in the haplotype diversity of 0.9909. The majority of haplotypes were assigned to East Asian lineages and the distribution of haplogroups of Muli Tibetan significantly differed from reference Tibetan populations. The maximum parsimony phylogeny reconstructed by 119 newly generated mitogenomes revealed 12 major Muli Tibetan lineages. Intriguingly, a Sherpa-specific sub-haplogroup A15c1 with the lack of mutations at 4216 and 15,924 was discerned in our dataset, which suggested that the maternal gene pool of Sherpas may derive from Tibetan populations. The shared haplogroups between Muli Tibetan and lowland Han Chinese hinted that these lineages may derive from non-Tibetans and have already differentiated before their arrival on the Tibetan Plateau. Furthermore, extensive pairwise population comparisons displayed that Muli Tibetan had a closer genetic relationship with ethnically or linguistically close Nyingtri Tibetan, Nyingtri Lhoba and Chamdo Tibetan populations. Genetic affinity was also observed between the Muli Tibetan and North Han Chinese. Collectively, the results generated in this study enriched the existing forensic mtDNA database and raised additional interest in the application of whole mitogenome sequencing in forensic investigations.

Discrimination of haplotype in mitochondrial DNA mixtures using LNA-mediated PCR clamping

Locked nucleic acid (LNA) has been widely used for various genetic analyses, and has many benefits, in terms of the specificity or sensitivity of amplification, because LNA-containing primers/probes form more stable duplexes with template DNA than probes lacking LNA. Here, we developed a new method for discriminating HV1 haplotypes from mitochondrial DNA (mtDNA) mixtures by applying PCR clamping using LNA. PCR clamping is based on the selective inhibition of amplification using LNA-containing probes, which can discriminate single-nucleotide differences. Before designing probes, we selected 171 sequences with single-nucleotide variations from the HV1 region, and evaluated the specificity of LNA-containing probes for them by predicting Tm values. The differences of Tm between mismatched and exactly matched probe-template duplexes depended markedly on the type of LNA nucleotides for discriminating single-nucleotide differences, and the cytosine LNA nucleotide at the site of variations in the probes was most effective to discriminate these differences. For mixture analysis, each probe targeted one or two variations (16209C, 16217C, 16257A/16261T, 16297C/16298C, 16304C, 16362C, or 16362T) that are particularly common in the Japanese population, and seven designed probes completely inhibited the amplification of exactly matched templates. We prepared mixed samples by mixing DNA from two individuals at a ratio of 1:9, 1:4, 1:1, 4:1, or 9:1, and then performed Sanger sequencing analysis after PCR clamping with each probe. Our method distinguished each haplotype at lower ratios from two-person mixtures, and enabled sensitive detection at 12 pg of total DNA including 600 copies of mtDNA. Moreover, we analyzed three-person mixtures with representative sequences, and detected the minor haplotype of one individual present at a rate of 10% by adding two selected probes. The ability to discriminate haplotypes in mixed samples by using LNA-mediated PCR clamping indicates the potential value of mtDNA analysis in criminal investigations.