The mitochondrial DNA makeup of Romanians: A forensic mtDNA control region database and phylogenetic characterization

Analysis of maternal genetic structure of mitochondrial DNA control region from Tai-Kadai-speaking Buyei population in southwestern China

BACKGROUND

Even though the Buyei are a recognised ethnic group in southwestern China, there hasn't been much work done on forensic population genetics, notably using mitochondrial DNA. The sequences and haplogroups of mitochondrial DNA control regions of the Buyei peoples were studied to provide support for the establishment of a reference database for forensic DNA analysis in East Asia.

METHODS AND RESULTS

The mitochondrial DNA control region sequences of 200 Buyei individuals in Guizhou were investigated. The haplotype frequencies and haplogroup distribution of the Buyei nationality in Guizhou were calculated. At the same time, the paired Fst values of the study population and other populations around the world were computed, to explore their genetic polymorphism and population relationship. A total of 179 haplotypes were detected in the Buyei population, with frequencies of 0.005-0.015. All haplotypes were assigned to 89 different haplogroups. The haplotype diversity and random matching probability were 0.999283 and 0.0063, respectively. The paired Fst genetic distances and correlation p-values among the 54 populations revealed that the Guizhou Buyei was most closely related to the Henan Han and the Guizhou Miao, and closer to the Hazara population in Pakistan and the Chiang Mai population.

CONCLUSIONS

The study of mitochondrial DNA based on the maternal genetic structure of the Buyei nationality in Guizhou will benefit the establishment of an East Asian forensic DNA reference database and provide a reference for anthropological research in the future.

Genetic analysis of the mitochondrial DNA control region in Tai-Kadai-speaking Dong population in southwest China

BACKGROUND

Dong people in Southwest China are officially recognised as an ethnic group, but there has been a lack of population genetic research on this group, especially based on mitochondrial DNA data.

AIM

To study the sequences and haplogroups of the mitochondrial DNA control region in a typical Dong population, and to provide help for the construction of a forensic mitochondrial DNA analysis reference database in East Asia.

SUBJECTS AND METHODS

The sequences of the mitochondrial DNA control region were analysed in 200 individuals of Dong in Guizhou. The haplotype frequencies, haplogroup distribution and paired Fst values of Guizhou Dong and 51 other populations in the world were calculated and explained to explore the genetic polymorphism and population relationships.

RESULTS

A total of 180 haplotypes were detected, with frequencies of 0.005-0.02. All haplotypes were assigned to 97 different haplogroups. The haplotype diversity and random matching probability were 0.998643 and 0.00635, respectively. The paired Fst values and correlation p values of 52 populations showed that the Guizhou Dong had the closest genetic relationship with the Henan Han and the Guizhou Miao in China, and were closest to the Punjab population in Pakistan and the Kashmiri population when compared with the world populations.

CONCLUSIONS

Our study was based on the matrilineal genetic structure of Guizhou Dong to study mitochondrial DNA, which was helpful to promote the establishment of the forensic DNA reference database in East Asia and provide reference for anthropological research.

Evaluating the suitability of current mitochondrial DNA interpretation guidelines for multigenerational whole mitochondrial genome comparisons

Sanger sequencing of the mitochondrial DNA (mtDNA) control region was previously the only method available for forensic casework involving degraded samples from skeletal remains. The introduction of Next Generation Sequencing (NGS) has transformed genetic data generation and human identification using mtDNA. Whole mitochondrial genome (mtGenome) analysis is now being introduced into forensic laboratories around the world to analyze historical remains. Research into large pedigrees using the mtGenome is critical to evaluate currently available interpretation guidelines for mtDNA analysis, which were developed for comparisons using the control region. This study included mtGenomes from 225 individuals from the last four generations of the Norfolk Island (NI) genetic isolate pedigree consisting of 49 distinct maternal lineages. The data from these individuals were arranged into 2339 maternally related pairs separated by up to 18 meioses. Our results show that 97.3% of maternally related pairs were concordant at all nucleotide positions, resulting in the correct interpretation of “Cannot Exclude”; 2.7% of pairs produced an “Inconclusive” result, and there were no instances of false exclusion. While these results indicate that existing guidelines are suitable for multigenerational whole mtGenome analysis, we recommend caution be taken when classifying heteroplasmic changes as differences for human identification. Our data showed the classification of heteroplasmic changes as differences increases the prevalence of inconclusive identification by 6%, with false exclusions observed in 0.34% of pairs examined. Further studies of multigenerational pedigrees, however, are needed to validate mtGenome interpretation guidelines for historical case work to more fully utilize emerging advancements.

The mitochondrial DNA control region sequences from the Chinese Sui population of southwestern China

BACKGROUND

Sui people are officially recognised people living in southwest China, but there has been a lack of genetic research, especially based on mitochondrial DNA data.

AIM

To study the sequences and haplogroups of the mitochondrial DNA control region in a typical Sui population, with the aim of helping to promote the establishment of a forensic DNA analysis reference database in East Asia.

SUBJECTS AND METHODS

We analysed 201 Sui individuals and observed the sequences of the mitochondrial DNA control region. We calculated and explained the haplotype frequencies, haplogroup distribution and pairwise Fst values between the Sui and 47 other populations in the world, in order to explore genetic polymorphisms and population relationships.

RESULTS

161 haplotypes were found in the Sui population, with frequencies of 0.0049-0.0199. All samples were assigned to 80 different haplogroups. The haplotype diversity and random matching probability were 0.999938 and 0.024729, respectively. The pairwise Fst values and correlation p-values of 48 populations showed that the Sui population was most closely related to the Miao population in Guizhou and the Han population in Henan, and closer to the Punjab population and Pukhtunkhwa population in Pakistan, and was significantly different from the other 43 groups. Compared with the other 43 groups, it is relatively isolated.

CONCLUSION

Our results show that the study of mitochondrial DNA based on the analysis of matrilineal genetic structure of the Sui population can help to promote the establishment of a forensic DNA reference database in East Asia and provide reference for future anthropological research.

Genetic characterization of a collection of Tsantsas from Ecuadorian museums

Tsantsas are shrunken human heads originally made for ceremonial purposes by Amazonian indigenous groups of the Shuar and Achuar family, previously called Jivaroan tribes. A significant demand of these objects during the first half of the 20th century led to the manufacture of counterfeit shrunken heads for commercial purposes. For museums where these collections are held, as well as for the indigenous groups who claim their ownership, it is important to identify the origin and authenticity of these tsantsas. We hypothesized that a collection of 14 tsantsas from 3 different museum collections in Ecuador are human and aimed to characterize their sex and potential origin. We amplified the amelogenin gene and performed a high resolution melting analysis to determine their human origin and characterize their sex. We also analyzed a fragment (16209-16402) from the HVR-1 region to identify the mtDNA haplogroups present in the tsantsa collection. Our exploratory results show that all the tsantsas are human and that the collection is comprised of 13 males and 1 female. A total of seven mtDNA haplogroups were found among the tsantsa collection using the mtDNA EMPOP database. These results show a predominance of the Amerindian mtDNA haplogroups B, C and D. Additional principal component analysis, genetic distance tree and haplotype network analyses suggest a relationship between the tsantsa specimens and Native American groups.

Mitochondrial DNA Profiles of Individuals from a 12th Century Necropolis in Feldioara (Transylvania)

The genetic signature of modern Europeans is the cumulated result of millennia of discrete small-scale exchanges between multiple distinct population groups that performed a repeated cycle of movement, settlement, and interactions with each other. In this study we aimed to highlight one such minute genetic cycle in a sea of genetic interactions by reconstructing part of the genetic story of the migration, settlement, interaction, and legacy of what is today the Transylvanian Saxon. The analysis of the mitochondrial DNA control region of 13 medieval individuals from Feldioara necropolis (Transylvania region, Romania) reveals a genetically heterogeneous group where all identified haplotypes are different. Most of the perceived maternal lineages are of Western Eurasian origin, except for the Central Asiatic haplogroup C seen in only one sample. Comparisons with historical and modern populations describe the contribution of the investigated Saxon settlers to the genetic history of this part of Europe.

The shaping of immunological responses through natural selection after the Roma Diaspora

The Roma people are the largest transnational ethnic minority in Europe and can be considered the last human migration of South Asian origin into the continent. They left Northwest India approximately 1,000 years ago, reaching the Balkan Peninsula around the twelfth century and Romania in the fourteenth century. Here, we analyze whole-genome sequencing data of 40 Roma and 40 non-Roma individuals from Romania. We performed a genome-wide scan of selection comparing Roma, their local host population, and a Northwestern Indian population, to identify the selective pressures faced by the Roma mainly after they settled in Europe. We identify under recent selection several pathways implicated in immune responses, among them cellular metabolism pathways known to be rewired after immune stimulation. We validated the interaction between PIK3-mTOR-HIF-1α and cytokine response influenced by bacterial and fungal infections. Our results point to a significant role of these pathways for host defense against the most prevalent pathogens in Europe during the last millennium.

Searching the undetected mtDNA variants in forensic MPS data

The efficiency of MPS in forensic mtDNA analysis has been thoroughly proven, although a reliable and well established data evaluation still remains a critical point. Numerous bioinformatics tools have been developed, but most of them require specific operating systems and high costs, while free open-source programs with user-friendly interfaces are few. In this study, 43 full mtGenomes were sequenced using the Ion Personal Genome Machine™ (PGM™) System and analyzed utilizing the plug-in Variant Caller (TVC) of the Ion Torrent Software Suite and the mtDNA-Server (mDS), a free web-based mitochondrial analysis tool for MPS data. The outcomes of these two different analysis tools were compared to variants noted after manual inspection of the aligned reads performed using Integrative Genomics Viewer (IGV). The comparison highlighted the presence of thirty-nine discordant variant calls, which were resolved by Sanger sequencing that confirmed the presence of all variants, except for 7 deletions. The combined adoption of IGV and Sanger type sequencing confirmatory steps, in addition of TVC and mDS analysis, resulted in a more accurate variants assignment with the detection of 32 additional true polymorphisms, which were noted in the final dataset. Regarding the heteroplasmy issue, out of a total of thirty heteroplasmic variants, twenty-eight were detected by the TVC, while the mDS detected twenty-two. Overall, none of the used bioinformatics tools were the perfect choice and a secondary analysis with an expert's opinion in complete mtGenome MPS data evaluation is still required in forensic genetic analysis.

Massively parallel sequencing of mitogenome sequences reveals the forensic features and maternal diversity of tai-kadai-speaking hlai islanders

As a single maternally inherited locus, human mitochondrial DNA (mtDNA) is geographically arranged and plays a key role in forensic applications. Hlai population has been evidenced as the most typical and unmixed representative of the Tai-Kadai-speaking populations via genome-wide analyses. However, forensic features and maternal diversity of the complete mitogenomes in this Tai-Kadai ancestrally related population are scarce. Thus, we sequenced the complete mitogenomes in 127 Hainan Hlais and found 109 distinct haplotypes belonging to 43 terminal haplogroups resulting in the haplotype diversity of 0.9970. Our results of comprehensive population comparisons showed that Hlai islanders had a close genetic affinity with Tai-Kadai-speaking populations from Southeast Asia, which is consistent with the back-migration of Chinese Neolithic farmers into this region via the inland route. Besides, maternally genetic evidence further revealed a close genetic relationship between Tai-Kadai-speaking and Austronesian-speaking populations when only East Asian dataset was considered, which is consistent with the common origin from Yangtze rice farmers and then spread southward along the inland and coastal routes, respectively. In the reconstructed phylogenetic tree and median-joining networks, the vast majority of Hlais were clustered in exclusive clades, which demonstrated that Hlai people probably had undergone founder effect or genetic bottleneck in their history, and remained genetically isolated for a long time. Collectively, Hainan Hlai did not exhibit detectable maternal gene flow from surrounding or incoming populations. Mitogenome information generated in this study is a contribution in mitigating the underrepresentation of Chinese data in forensic mitogenetics and will assist geography-, metapopulation-, as well as phylogeny-based queries.

The mitochondrial DNA control region sequences from the Chinese Miao population of southeastern China

Background: Miao people are an officially recognised ethnic group living in southwest China, but have seldom been studied genetically, especially with respect to mtDNA data.Aim: To investigate the sequences and haplogroups of the mtDNA control region in a typical Miao population, with the aim of providing a good start for the expansion of the East Asian mtDNA reference database for forensic DNA analysis.Subjects and methods: We analysed 203 Miao individuals, looking at mtDNA control region sequences. We calculated and illustrated the haplotype frequencies, haplogroup distribution and pairwise Fst values between the Miao and six other worldwide populations to explore genetic polymorphisms and population relationships.Results: We observed 121 haplotypes with corresponding frequencies ranging from 0.0049 to 0.0690 in the Miao population. All the samples were assigned to 71 different haplogroups. The haplotype diversity and the random match probability were estimated to be 0.9844 and 0.0204, respectively. The pairwise Fst values and associated p values among seven populations suggest that the Miao population has significant differences to the other six populations, and is relatively isolated compared with them.Conclusions: Our results suggest that frequency estimates for mtDNA haplotypes in Miao ethnic groups should be determined independently rather than being pooled with other populations.

A probable case of infantile cortical hyperostosis in 2nd-4th centuries AD Romania

OBJECTIVE

This study aims to discuss the differential diagnosis for the pathological alterations displayed on an infant skeleton from Romania.

MATERIALS

One infant skeleton retrieved form the bathhouse of an abandoned Roman fort and dated between the 2nd and the 4th centuries AD.

METHODS

All available skeletal elements were analyzed macroscopically. In addition, the isotopic signatures (δ13C and δ15N) and the control region of the human mitochondrial genome for this archaeological sample were analyzed.

RESULTS

Based on dental development and long bone length, the skeleton was aged between birth and 2 months of age. Pathological lesions were noted on the mandible and diaphyses of long bones, but spared the metaphyses.

CONCLUSIONS

The perinatal age of the individual, along with lesion morphology and location, suggests a diagnosis of infantile cortical hyperostosis.

LIMITATIONS

The analysis would benefit from further stable isotope and mitochondrial genome analyses, which was limited due to the absence of comparative human and faunal remains from the site.

SUGGESTIONS FOR FURTHER RESEARCH

Further multidisciplinary research on human archaeological remains from Romania would provide a clearer image of past disease and life histories in this geographic area.

Mitochondrial DNA in human identification: a review

Mitochondrial DNA (mtDNA) presents several characteristics useful for forensic studies, especially related to the lack of recombination, to a high copy number, and to matrilineal inheritance. mtDNA typing based on sequences of the control region or full genomic sequences analysis is used to analyze a variety of forensic samples such as old bones, teeth and hair, as well as other biological samples where the DNA content is low. Evaluation and reporting of the results requires careful consideration of biological issues as well as other issues such as nomenclature and reference population databases. In this work we review mitochondrial DNA profiling methods used for human identification and present their use in the main cases of humanidentification focusing on the most relevant issues for forensics.

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.

A genetic portrait of the South-Eastern Carpathians based on autosomal short tandem repeats loci used in forensics

OBJECTIVES

This work aimed to describe the genetic landscape of the Balkan Peninsula, as revealed by STR markers commonly used in forensics and spatial methods specifically developed for genetic data.

METHODS

We generated and analyzed 16 short tandem repeats (STRs) autosomal genotypes in 287 subjects from ten administrative/geographical regions of Eastern Europe (Romania and the Republic of Moldova). We report estimates of the allele frequencies in these sub-populations, their fixation indexes, and use these results to complement previous spatial analyses of Southern Europe.

RESULTS

In seven out of ten analyzed regional samples the heterozygosity, averaged across loci, was lower than expected. The average Fis was 0.011. Among the 16 loci, five returned a significant fixation index Fst. The composite Fst across the 16 loci, among the 10 regional samples, was 0.00417, a figure twice as large as that obtained with the same markers across the entire Northern Mediterranean. The first spatial principal component (sPC1) returned the picture of a Central-European pattern of frequencies for the Carpathians, which extended to the Southern boundary of the Balkan Peninsula. However, the 8 alleles extracted by sPC1 returned a picture of a strong reduction of the migration rate in the Carpathian region, mostly between the inner locations.

CONCLUSIONS

Our results revealed an unexpected heterogeneity in the area. We believe that populations from some regions will require treatment as distinct entities when considered in forensic applications.

Maternal DNA lineages at the gate of Europe in the 10th century AD

Given the paucity of archaeogenetic data available for medieval European populations in comparison to other historical periods, the genetic landscape of this age appears as a puzzle of dispersed, small, known pieces. In particular, Southeastern Europe has been scarcely investigated to date. In this paper, we report the study of mitochondrial DNA in 10^th century AD human samples from Capidava necropolis, located in Dobruja (Southeastern Romania, Southeastern Europe). This geographical region is particularly interesting because of the extensive population flux following diverse migration routes, and the complex interactions between distinct population groups during the medieval period. We successfully amplified and typed the mitochondrial control region of 10 individuals. For five of them, we also reconstructed the complete mitochondrial genomes using hybridization-based DNA capture combined with Next Generation Sequencing. We have portrayed the genetic structure of the Capidava medieval population, represented by 10 individuals displaying 8 haplotypes (U5a1c2a, V1a, R0a2’3, H1, U3a, N9a9, H5e1a1, and H13a1a3). Remarkable for this site is the presence of both Central Asiatic (N9a) and common European mtDNA haplotypes, establishing Capidava as a point of convergence between East and West. The distribution of mtDNA lineages in the necropolis highlighted the existence of two groups of two individuals with close maternal relationships as they share the same haplotypes. We also sketch, using comparative statistical and population genetic analyses, the genetic relationships between the investigated dataset and other medieval and modern Eurasian populations.

Genetic affinities among the historical provinces of Romania and Central Europe as revealed by an mtDNA analysis

Background

As a major crossroads between Asia and Europe, Romania has experienced continuous migration and invasion episodes. The precise routes may have been shaped by the topology of the territory and had diverse impacts on the genetic structure of mitochondrial DNA (mtDNA) in historical Romanian provinces. We studied 714 Romanians from all historical provinces, Wallachia, Dobrudja, Moldavia, and Transylvania, by analyzing the mtDNA control region and coding markers to encompass the complete landscape of mtDNA haplogroups.

Results

We observed a homogenous distribution of the majority of haplogroups among the Romanian provinces and a clear association with the European populations. A principal component analysis and multidimensional scaling analysis supported the genetic similarity of the Wallachia, Moldavia, and Dobrudja groups with the Balkans, while the Transylvania population was closely related to Central European groups. These findings could be explained by the topology of the Romanian territory, where the Carpathian Arch played an important role in migration patterns. Signals of Asian maternal lineages were observed in all Romanian historical provinces, indicating gene flow along the migration routes through East Asia and Europe.

Conclusions

Our current findings based on the mtDNA analysis of populations in historical provinces of Romania suggest similarity between populations in Transylvania and Central Europe, supported both by the observed clines in haplogroup frequencies for several European and Asian maternal lineages and MDS analyses.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0487-5) contains supplementary material, which is available to authorized users.