A single multiplex PCR and SNaPshot minisequencing reaction of 42 SNPs to classify admixture populations into mitochondrial DNA haplogroups

Forensic efficiency evaluation of a mtDNA whole genome sequencing system constructed with long fragment amplification strategy on DNA nanoball sequencing platform

Mitochondrial DNA (mtDNA) is an important genetic marker for degraded biological sample identification, maternal pedigree tracing, and population genetic structure study owing to its characteristics of high copy number, anti-degradable ring structure, and maternal inheritance. Whole mtDNA genome sequencing is an optimal method for the analysis of mtDNA polymorphism and heterogeneity because it allows for the comprehensive use of maternal genetic information. However, because of lacking quantitative evaluations for sequencing data, the scientific interpretation standards for mtDNA sequencing results of the previously used sequencing systems are often different, and false positive or false negative results are prone to occur when faced with the interference of nuclear genomic DNA, or the heterogeneities of mtDNA sequence and structure. In this study, we evaluated a novel mtDNA whole genome sequencing system using long fragment amplification strategy on the DNA nanoball (DNB) sequencing platform. This system demonstrated high sequencing quality and specific mtDNA sequencing efficiencies on positive control DNA and FTA bloodstain samples, as the average Q20 and Q30 values of the corresponding samples were 97.17 % and 91.93 %; 97.37 % and 92.48 %, respectively. The mean mapping percentages for the reference sequences of whole genome DNA (wgDNA), mtDNA, and nuclear genomic DNA (ngDNA) in the corresponding samples were 99.98 %, 99.97 %, 0.03 %, and 99.91 %, 99.40 %, 0.60 %; respectively. The average error calling rates for the bases A, C, G, and T of the whole mtDNA genome were 0.2519 %, 0.2550 %, 0.2906 %; and 0.2392 %, respectively. The efficacy of heteroplasmy identification was assessed using a set of theoretical sites with predetermined rates. These sites were created by combining the samples with known mtDNA haplotypes in certain proportions. The absolute errors between observed and theoretical heteroplasmy values were 89.59 %, 74.68 %, 50.20 %, 12.65 %, 8.31 %, and 4.85 %, while the theoretical heteroplasmy values were 5 %, 10 %, 20 %, 80 %, 90 %, and 95 %, respectively. The absolute error exhibited relative stability when the mtDNA sequencing depth exceeded 500×. Furthermore, the system sequencing efficiency was also confirmed among different kinds of samples, and these samples included natural samples (e.g., peripheral blood samples preserved on FTA cards for 2 and 11 years, and on filter paper for 6 and 9 years), degraded samples, sensitivity samples, samples derived from various bodily fluids, and maternal pedigree samples. In summary, the whole mtDNA genome sequencing system used for forensic identification demonstrated high performance in analyzing mtDNA sequence information, and showed significant prospects for forensic application and maternal genetic research.

Hi-MC: a novel method for high-throughput mitochondrial haplogroup classification

Effective approaches for assessing mitochondrial DNA (mtDNA) variation are important to multiple scientific disciplines. Mitochondrial haplogroups characterize branch points in the phylogeny of mtDNA. Several tools exist for mitochondrial haplogroup classification. However, most require full or partial mtDNA sequence which is often cost prohibitive for studies with large sample sizes. The purpose of this study was to develop Hi-MC, a high-throughput method for mitochondrial haplogroup classification that is cost effective and applicable to large sample sizes making mitochondrial analysis more accessible in genetic studies. Using rigorous selection criteria, we defined and validated a custom panel of mtDNA single nucleotide polymorphisms that allows for accurate classification of European, African, and Native American mitochondrial haplogroups at broad resolution with minimal genotyping and cost. We demonstrate that Hi-MC performs well in samples of European, African, and Native American ancestries, and that Hi-MC performs comparably to a commonly used classifier. Implementation as a software package in R enables users to download and run the program locally, grants greater flexibility in the number of samples that can be run, and allows for easy expansion in future revisions. Hi-MC is available in the CRAN repository and the source code is freely available at https://github.com/vserch/himc.

Development of a new screening method to determine the main 52 mitochondrial haplogroups through a single minisequencing reaction

This work presents the design, development and optimization of a screening method based on single-base extension sequencing to simultaneously analyze a panel of 52 mitochondrial SNPs. This enables to recognize the main mitochondrial haplogroups and to discriminate even between lineages from the same phylogenetic branch that diverged in different continents. The unavailability of individuals harboring infrequent variants was a limitation to optimize the panel. To overcome this, we have modified DNA by site-directed mutagenesis to create the unavailable allelic variants. This allowed us to verify the reliability of this panel and its usefulness to be applied in biomedicine, forensic and population genetic studies.

Polymorphisms of mitochondrial DNA control region are associated to endometriosis

PURPOSE

Polymorphisms in the control region of mitochondrial DNA (mtDNA) can affect generation of reactive oxygen species and impact in the pathogenesis of endometriosis. This study investigated the association of mtDNA polymorphisms with endometriosis.

METHODS

Patients were divided in two groups: endometriosis (n = 90) and control (n = 92). Inclusion criteria were as follows: women between 18 and 50 years, with histological diagnosis and surgical staging of endometriosis (endometriosis group) or undergoing gynecological surgery for tubal ligation, leiomyoma, or ovarian cysts, with no evidence of endometriosis (control group). DNA extraction was performed from peripheral blood. Sanger sequencing of mtDNA control region was performed, and polymorphisms were determined comparing the sequences obtained with the Cambridge Reference Sequence.

RESULTS

The frequency of polymorphisms T16217C (14.4 and 5.4% of endometriosis and control group, respectively; p = 0.049) and G499A (13.3 vs. 4.3%; p = 0.038) was higher in the endometriosis group, while T146C (32.6 vs. 18.9%; p = 0.042) and 573.2C (5.6 vs. 29.3%; p < 0.001) were lower. No difference was observed in haplogroups between groups.

CONCLUSION

mtDNA polymorphisms T16217C and G499A were associated with endometriosis, while T416C and 573.2C were shown to be associated with an absence of disease.

Mitochondrial Haplogroups and Risk of Pulmonary Arterial Hypertension

Pulmonary arterial hypertension (PAH) is a serious and often fatal disease. It is a panvasculopathy of the pulmonary microcirculation characterized by vasoconstriction and arterial obstruction due to vascular proliferation and remodeling and ultimately right ventricular failure. Mitochondrial dysfunction is a universal finding in pulmonary vascular cells of patients with PAH, and is mechanistically linked to disease origins in animal models of pulmonary hypertension. Mitochondria have their own circular DNA (mtDNA), which can be subgrouped into polymorphic haplogroup variants, some of which have been identified as at-risk or protective from cardiovascular and/or neurodegenerative diseases. Here, we hypothesized that mitochondrial haplogroups may be associated with PAH. To test this, mitochondrial haplogroups were determined in a cohort of PAH patients and controls [N = 204 Caucasians (125 PAH and 79 controls) and N = 46 African Americans (13 PAH and 33 controls)]. Haplogroup L was associated with a lower rate of PAH as compared to macrohaplogroups N and M. When haplogroups were nested based on ancestral inheritance and controlled for age, gender and race, haplogroups M and HV, JT and UK of the N macro-haplogroup had significantly higher rates of PAH compared to the ancestral L (L0/1/2 and L3) (all p ≤ 0.05). Overall, the findings suggest that mitochondrial haplogroups influence risk of PAH and that a vulnerability to PAH may have emerged under the selective enrichment of specific haplogroups that occurred with the migration of populations out of Africa.

A Simple Algorithm for Population Classification

A single-nucleotide polymorphism (SNP) is a variation in the DNA sequence that occurs when a single nucleotide in the genome differs across members of the same species. Variations in the DNA sequences of humans are associated with human diseases. This makes SNPs as a key to open up the door of personalized medicine. SNP(s) can also be used for human identification and forensic applications. Compared to short tandem repeat (STR) loci, SNPs have much lower statistical testing power for individual recognition due to the fact that there are only 3 possible genotypes for each SNP marker, but it may provide sufficient information to identify the population to which a certain samples may belong. In this report, using eight SNP markers for 641 samples, we performed a standard statistical classification procedure and found that 86% of the samples could be classified accurately under a two-population model. This study suggests the potential use of SNP(s) in population classification with a small number (n ≤ 8) of genetic markers for forensic screening, biodiversity and disaster victim controlling.

Development of Multiple Assays using 46 SNPs for Comprehensive mtDNA Haplogrouping and Application to Highly Degraded DNA

Six multiplex PCR systems using single-base extension reactions to analyze 46 mitochondrial DNA (mtDNA)-coding region single nucleotide polymorphisms (SNPs) that define 42 haplogroups, that is, 24 major mtDNA haplogroups and 18 subclades, were devised. To improve the usefulness of the established systems for the analysis of degraded DNA samples, novel primers to render amplicons with sizes <150 bp were designed. By applying these systems to 214 Japanese individuals, 24 different haplogroups (power of discrimination = 93.4%) were found. To assess the effectiveness of our systems in grouping degraded DNA, an ancient bone sample of a Jomon skeleton was analyzed and then classified as haplogroup N9b. We conclude that the present systems are powerful screening tools for major haplogroups of mtDNA in addition to the prevalent subhaplogroups in the Japanese population and that these systems are capable of analyzing highly degraded DNA samples in forensic studies.

Haplotype diversity in mitochondrial DNA hypervariable region in a population of southeastern Brazil

Brazilian population derives from Native Amerindians, Europeans, and Africans. Southeastern Brazil is the most populous region of the country. The present study intended to characterize the maternal genetic ancestry of 290 individuals from southeastern (Brazil) population. Thus, we made the sequencing of the three hypervariable regions (HV1, HV2, and HV3) of the mitochondrial DNA (mtDNA). The statistical analyses were made using Arlequin software, and the median-joining haplotype networks were generated using Network software. The analysis of three hypervariable regios showed 230 (79.3 %) unique haplotypes and the most common haplotype was "263G" carried by 12 (4.1 %) individuals. The strikingly high variability generated by intense gene flow is mirrored in a high sequence diversity (0.9966 ± 0.0010), and the probability of two random individuals showing identical mtDNA haplotypes were 0.0068. The analysis of haplogroup distribution revealed that 36.9 % (n = 107) presented Amerindian haplogroups, 35.2 % (n = 102) presented African haplogroups, 27.6 % (n = 80) presented European haplogroups, and one (0.3 %) individual presented East Asian haplogroup, evidencing that the southeastern population is extremely heterogeneous and the coexistence of matrilineal lineages with three different phylogeographic origins. The genetic diversity found in the mtDNA control region in the southeastern Brazilian population reinforces the importance of increased national database in order to be important and informative in forensic cases.

AmericaPlex26: a SNaPshot multiplex system for genotyping the main human mitochondrial founder lineages of the Americas

Phylogeographic studies have described a reduced genetic diversity in Native American populations, indicative of one or more bottleneck events during the peopling and prehistory of the Americas. Classical sequencing approaches targeting the mitochondrial diversity have reported the presence of five major haplogroups, namely A, B, C, D and X, whereas the advent of complete mitochondrial genome sequencing has recently refined the number of founder lineages within the given diversity to 15 sub-haplogroups. We developed and optimized a SNaPshot assay to study the mitochondrial diversity in pre-Columbian Native American populations by simultaneous typing of 26 single nucleotide polymorphisms (SNPs) characterising Native American sub-haplogroups. Our assay proved to be highly sensitive with respect to starting concentrations of target DNA and could be applied successfully to a range of ancient human skeletal material from South America from various time periods. The AmericaPlex26 is a powerful assay with enhanced phylogenetic resolution that allows time- and cost-efficient mitochondrial DNA sub-typing from valuable ancient specimens. It can be applied in addition or alternative to standard sequencing of the D-loop region in forensics, ancestry testing, and population studies, or where full-resolution mitochondrial genome sequencing is not feasible.

Genotyping of BCL11A and HBS1L-MYB SNPs associated with fetal haemoglobin levels: a SNaPshot minisequencing approach

BACKGROUND

B-thalassaemia and sickle cell disease (SCD) are two of the most common monogenic diseases that are found in many populations worldwide. In both disorders the clinical severity is highly variable, with the persistence of fetal haemoglobin (HbF) being one of the major ameliorating factors. HbF levels are affected by, amongst other factors, single nucleotide polymorphisms (SNPs) at the BCL11A gene and the HBS1L-MYB intergenic region, which are located outside the β-globin locus. For this reason, we developed two multiplex assays that allow the genotyping of SNPs at these two genomic regions which have been shown to be associated with variable HbF levels in different populations.

RESULTS

Two multiplex assays based on the SNaPshot minisequencing approach were developed. The two assays can be used to simultaneous genotype twelve SNPs at the BCL11A gene and sixteen SNPs at HBS1L-MYB intergenic region which were shown to modify HbF levels. The different genotypes can be determined based on the position and the fluorescent colour of the peaks in a single electropherogram. DNA sequencing and restriction fragment length polymorphism (PCR-RFLP) assays were used to verify genotyping results obtained by SNaPshot minisequencing.

CONCLUSIONS

In summary, we propose two multiplex assays based on the SNaPshot minisequencing approach for the simultaneous identification of SNPs located at the BCL11A gene and HBS1L-MYB intergenic region which have an effect on HbF levels. The assays can be easily applied for accurate, time and cost efficient genotyping of the selected SNPs in various populations.

Characterization of mitochondrial haplogroups in a large population-based sample from the United States

Mitochondrial DNA (mtDNA) haplogroups are valuable for investigations in forensic science, molecular anthropology, and human genetics. In this study, we developed a custom panel of 61 mtDNA markers for high-throughput classification of European, African, and Native American/Asian mitochondrial haplogroup lineages. Using these mtDNA markers, we constructed a mitochondrial haplogroup classification tree and classified 18,832 participants from the National Health and Nutrition Examination Surveys (NHANES). To our knowledge, this is the largest study to date characterizing mitochondrial haplogroups in a population-based sample from the United States, and the first study characterizing mitochondrial haplogroup distributions in self-identified Mexican Americans separately from Hispanic Americans of other descent. We observed clear differences in the distribution of maternal genetic ancestry consistent with proposed admixture models for these subpopulations, underscoring the genetic heterogeneity of the United States Hispanic population. The mitochondrial haplogroup distributions in the other self-identified racial/ethnic groups within NHANES were largely comparable to previous studies. Mitochondrial haplogroup classification was highly concordant with self-identified race/ethnicity (SIRE) in non-Hispanic whites (94.8 %), but was considerably lower in admixed populations including non-Hispanic blacks (88.3 %), Mexican Americans (81.8 %), and other Hispanics (61.6 %), suggesting SIRE does not accurately reflect maternal genetic ancestry, particularly in populations with greater proportions of admixture. Thus, it is important to consider inconsistencies between SIRE and genetic ancestry when performing genetic association studies. The mitochondrial haplogroup data that we have generated, coupled with the epidemiologic variables in NHANES, is a valuable resource for future studies investigating the contribution of mtDNA variation to human health and disease.

MtDNA ancestry of Rio de Janeiro population, Brazil

Polymorphism studies concerning HVI and HVII regions of mitochondrial DNA (mtDNA) have improved the understanding of the admixture genetic process related to the occupation of the continents by human population groups. We have analyzed the mtDNA lineages of 190 healthy and maternally unrelated individuals born in the metropolitan region of the Rio de Janeiro city, the capital of the State of Rio de Janeiro, southeastern Brazil. The data showing that 57.9, 25.3 and 16.8 % of the matrilineages found in Rio de Janeiro come from African, Amerindian and European population groups. They are, respectively, in close agreement with historical records which indicate that the admixture population of Brazil is the resulting of interethnic asymmetry crosses between individuals from those population groups. The high proportion of African mtDNA lineages in the population of Rio de Janeiro is in accordance with studies related to other Brazilian states.

Development and validation of a D-loop mtDNA SNP assay for the screening of specimens in forensic casework

Mitochondrial DNA (mtDNA) analysis is usually a last resort in routine forensic DNA casework. However, it has become a powerful tool for the analysis of highly degraded samples or samples containing too little or no nuclear DNA, such as old bones and hair shafts. The gold standard methodology still constitutes the direct sequencing of polymerase chain reaction (PCR) products or cloned amplicons from the HVS-1 and HVS-2 (hypervariable segment) control region segments. Identifications using mtDNA are time consuming, expensive and can be very complex, depending on the amount and nature of the material being tested. The main goal of this work is to develop a less labour-intensive and less expensive screening method for mtDNA analysis, in order to aid in the exclusion of non-matching samples and as a presumptive test prior to final confirmatory DNA sequencing. We have selected 14 highly discriminatory single nucleotide polymorphisms (SNPs) based on simulations performed by Salas and Amigo (2010) to be typed using SNaPShot(TM) (Applied Biosystems, Foster City, CA, USA). The assay was validated by typing more than 100 HVS-1/HVS-2 sequenced samples. No differences were observed between the SNP typing and DNA sequencing when results were compared, with the exception of allelic dropouts observed in a few haplotypes. Haplotype diversity simulations were performed using 172 mtDNA sequences representative of the Brazilian population and a score of 0.9794 was obtained when the 14 SNPs were used, showing that the theoretical prediction approach for the selection of highly discriminatory SNPs suggested by Salas and Amigo (2010) was confirmed in the population studied. As the main goal of the work is to develop a screening assay to skip the sequencing of all samples in a particular case, a pair-wise comparison of the sequences was done using the selected SNPs. When both HVS-1/HVS-2 SNPs were used for simulations, at least two differences were observed in 93.2% of the comparisons performed. The assay was validated with casework samples. Results show that the method is straightforward and can be used for exclusionary purposes, saving time and laboratory resources. The assay confirms the theoretic prediction suggested by Salas and Amigo (2010). All forensic advantages, such as high sensitivity and power of discrimination, as also the disadvantages, such as the occurrence of allele dropouts, are discussed throughout the article.

A rapid method for simultaneous screening of multi-gene mutations associated with hearing loss in the Korean population

Hearing loss (HL) is a congenital disease with a high prevalence, and patients with hearing loss need early diagnosis for treatment and prevention. The GJB2, MT-RNR1, and SLC26A4 genes have been reported as common causative genes of hearing loss in the Korean population and some mutations of these genes are the most common mutations associated with hearing loss. Accordingly, we developed a method for the simultaneous detection of seven mutations (c.235delC of GJB2, c.439A>G, c.919-2A>G, c.1149+3A>G, c.1229C>T, c.2168A>G of SLC26A4, and m.1555A>G of the MT-RNR1 gene) using multiplex SNaPshot minisequencing to enable rapid diagnosis of hereditary hearing loss. This method was confirmed in patients with hearing loss and used for genetic diagnosis of controls with normal hearing and neonates. We found that 4.06% of individuals with normal hearing and 4.32% of neonates were heterozygous carriers. In addition, we detected that an individual is heterozygous for two different mutations of GJB2 and SLC26A4 gene, respectively and one normal hearing showing the heteroplasmy of m.1555A>G. These genotypes corresponded to those determined by direct sequencing. Overall, we successfully developed a robust and cost-effective diagnosis method that detects common causative mutations of hearing loss in the Korean population. This method will be possible to detect up to 40% causative mutations associated with prelingual HL in the Korean population and serve as a useful genetic technique for diagnosis of hearing loss for patients, carriers, neonates, and fetuses.

MtDNA SNP multiplexes for efficient inference of matrilineal genetic ancestry within Oceania

Human mitochondrial DNA (mtDNA) is a convenient marker for tracing matrilineal bio-geographic ancestry and is widely applied in forensic, genealogical and anthropological studies. In forensic applications, DNA-based ancestry inference can be useful for finding unknown suspects by concentrating police investigations in cases where autosomal STR profiling was unable to provide a match, or can help provide clues in missing person identification. Although multiplexed mtDNA single nucleotide polymorphism (SNP) assays to infer matrilineal ancestry at a (near) continental level are already available, such tools are lacking for the Oceania region. Here, we have developed a hierarchical system of three SNaPshot multiplexes for genotyping 26 SNPs defining all major mtDNA haplogroups for Oceania (including Australia, Near Oceania and Remote Oceania). With this system, it was possible to conclusively assign 74% of Oceanian individuals to their Oceanian matrilineal ancestry in an established literature database (after correcting for obvious external admixture). Furthermore, in a set of 161 genotyped individuals collected in Australia, Papua New Guinea and Fiji, 87.6% were conclusively assigned an Oceanian matrilineal origin. For the remaining 12.4% of the genotyped samples either a Eurasian origin was detected indicating likely European admixture (1.9%), the identified haplogroups are shared between Oceania and S/SE-Asia (5%), or the SNPs applied did not allow a geographic inference to be assigned (5.6%). Sub-regional assignment within Oceania was possible for 32.9% of the individuals genotyped: 49.5% of Australians were assigned an Australian origin and 13.7% of the Papua New Guineans were assigned a Near Oceanian origin, although none of the Fijians could be assigned a specific Remote Oceanian origin. The low assignment rates of Near and Remote Oceania are explained by recent migrations from Asia via Near Oceania into Remote Oceania. Combining the mtDNA multiplexes for Oceania introduced here with those we developed earlier for all other continental regions, global matrilineal bio-geographic ancestry assignment from DNA is now achievable in a highly efficient way that is also suitable for applications with limited material such as forensic case work.

An allelic ladder based upon reference alleles for mtDNA SNP analysis using the SNaPshot technique

The analysis of mitochondrial DNA (mtDNA) single-nucleotide polymorphisms (SNPs) using the SNaPshot technique (Applied Biosystems) is a fast and sensitive method for the reliable identification of disease-associated mtDNA SNPs, genetic ancestry mtDNA SNPs and forensically important mtDNA SNPs. The detection of many SNPs in one multiplex PCR and one subsequent multiplex minisequencing reaction is challenging for laboratories who want to establish this technique, due to the problem that there is no allelic ladder available for mtDNA SNP analysis via SNaPshot technique. Normally, the laboratory has to invent long-term testing and studies. The interpretation of false and correct alleles is up to some specialists knowing the expected and the estimated size of each allele SNP. We here present a protocol to assemble up to 84 alleles of 42 different mtDNA SNPs in an allelic ladder that is based upon reference alleles. We recommend using allelic ladders/reference alleles for SNP analysis to maintain high-quality analysis standards.