Generating population data for the EMPOP database - an overview of the mtDNA sequencing and data evaluation processes considering 273 Austrian control region sequences as example

Refining the Global Phylogeny of Mitochondrial N1a, X, and HV2 Haplogroups Based on Rare Mitogenomes from Croatian Isolates

Mitochondrial DNA (mtDNA) has been used for decades as a predominant tool in population genetics and as a valuable addition to forensic genetic research, owing to its unique maternal inheritance pattern that enables the tracing of individuals along the maternal lineage across numerous generations. The dynamic interplay between evolutionary forces, primarily genetic drift, bottlenecks, and the founder effect, can exert significant influence on genetic profiles. Consequently, the Adriatic islands have accumulated a subset of lineages that exhibits remarkable absence or rarity within other European populations. This distinctive genetic composition underscores the islands' potential as a significant resource in phylogenetic research, with implications reaching beyond regional boundaries to contribute to a global understanding. In the initial attempt to expand the mitochondrial forensic database of the Croatian population with haplotypes from small isolated communities, we sequenced mitogenomes of rare haplogroups from different Croatian island and mainland populations using next-generation sequencing (NGS). In the next step and based on the obtained results, we refined the global phylogeny of haplogroup N1a, HV2, and X by analyzing rare haplotypes, which are absent from the current phylogenetic tree. The trees were based on 16 novel and 52 previously published samples, revealing completely novel branches in the X and HV2 haplogroups and a new European cluster in the ancestral N1a variant, previously believed to be an exclusively African-Asian haplogroup. The research emphasizes the importance of investigating geographically isolated populations and their unique characteristics within a global context.

Helena's Many Daughters: More Mitogenome Diversity behind the Most Common West Eurasian mtDNA Control Region Haplotype in an Extended Italian Population Sample

The high number of matching haplotypes of the most common mitochondrial (mt)DNA lineages are considered to be the greatest limitation for forensic applications. This study investigates the potential to solve this constraint by massively parallel sequencing a large number of mitogenomes that share the most common West Eurasian mtDNA control region (CR) haplotype motif (263G 315.1C 16519C). We augmented a pilot study on 29 to a total of 216 Italian mitogenomes that represents the largest set of the most common CR haplotype compiled from a single country. The extended population sample confirmed and extended the huge coding region diversity behind the most common CR motif. Complete mitogenome sequencing allowed for the detection of 163 distinct haplotypes, raising the power of discrimination from 0 (CR) to 99.6% (mitogenome). The mtDNAs were clustered into 61 named clades of haplogroup H and did not reveal phylogeographic trends within Italy. Rapid individualization approaches for investigative purposes are limited to the most frequent H clades of the dataset, viz. H1, H3, and H7.

Exploring statistical weight estimates for mitochondrial DNA matches involving heteroplasmy

Massively parallel sequencing (MPS) of mitochondrial (mt) DNA allows forensic laboratories to report heteroplasmy on a routine basis. Statistical approaches will be needed to determine the relative frequency of observing an mtDNA haplotype when including the presence of a heteroplasmic site. Here, we examined 1301 control region (CR) sequences, collected from individuals in four major population groups (European, African, Asian, and Latino), and covering 24 geographically distributed haplogroups, to assess the rates of point heteroplasmy (PHP) on an individual and nucleotide position (np) basis. With a minor allele frequency (MAF) threshold of 2%, the data was similar across population groups, with an overall PHP rate of 37.7%, and the majority of heteroplasmic individuals (77.3%) having only one site of heteroplasmy. The majority (75.2%) of identified PHPs had an MAF of 2-10%, and were observed at 12.6% of the nps across the CR. Both the broad and phylogenetic testing suggested that in many cases the low number of observations of heteroplasmy at any one np results in a lack of statistical association. The posterior frequency estimates, which skew conservative to a degree depending on the sample size in a given haplogroup, had a mean of 0.152 (SD 0.134) and ranged from 0.031 to 0.83. As expected, posterior frequency estimates decreased in accordance with 1/n as the sample size (n) increased. This provides a proposed conservative statistical framework for assessing haplotype/heteroplasmy matches when applying an MPS technique in forensic cases and will allow for continual refinement as more data is generated, both within the CR and across the mitochondrial genome.

Development of the VISAGE enhanced tool and statistical models for epigenetic age estimation in blood, buccal cells and bones

DNA methylation is known as a biomarker for age with applications in forensics. Here we describe the VISAGE (VISible Attributes through GEnomics) Consortium's enhanced tool for epigenetic age estimation in somatic tissues. The tool is based on eight DNA methylation markers (44 CpGs), bisulfite multiplex PCR followed by sequencing on the MiSeq FGx platform, and three statistical prediction models for blood, buccal cells and bones. The model for blood is based on six CpGs from ELOVL2, MIR29B2CHG, KLF14, FHL2, TRIM59 and PDE4C, and predicts age with a mean absolute error (MAE) of 3.2 years, while the model for buccal cells includes five CpGs from PDE4C, MIR29B2CHG, ELOVL2, KLF14 and EDARADD and predicts age with MAE of 3.7 years, and the model for bones has six CpGs from ELOVL2, KLF14, PDE4C and ASPA and predicts age with MAE of 3.4 years. The VISAGE enhanced tool for age estimation in somatic tissues enables reliable collection of DNA methylation data from small amounts of DNA using a sensitive multiplex MPS assay that provides accurate estimation of age in blood, buccal swabs, and bones using the statistical model tailored to each tissue.

Evidence for multi-copy Mega-NUMTs in the human genome

The maternal mode of mitochondrial DNA (mtDNA) inheritance is central to human genetics. Recently, evidence for bi-parental inheritance of mtDNA was claimed for individuals of three pedigrees that suffered mitochondrial disorders. We sequenced mtDNA using both direct Sanger and Massively Parallel Sequencing in several tissues of eleven maternally related and other affiliated healthy individuals of a family pedigree and observed mixed mitotypes in eight individuals. Cells without nuclear DNA, i.e. thrombocytes and hair shafts, only showed the mitotype of haplogroup (hg) V. Skin biopsies were prepared to generate ρ° cells void of mtDNA, sequencing of which resulted in a hg U4c1 mitotype. The position of the Mega-NUMT sequence was determined by fluorescence in situ hybridization and two different quantitative PCR assays were used to determine the number of contributing mtDNA copies. Thus, evidence for the presence of repetitive, full mitogenome Mega-NUMTs matching haplogroup U4c1 in various tissues of eight maternally related individuals was provided. Multi-copy Mega-NUMTs mimic mixtures of mtDNA that cannot be experimentally avoided and thus may appear in diverse fields of mtDNA research and diagnostics. We demonstrate that hair shaft mtDNA sequencing provides a simple but reliable approach to exclude NUMTs as source of misleading results.

Complete mitogenome data for the Serbian population: the contribution to high-quality forensic databases

Mitochondrial genome (mtDNA) is a valuable resource in resolving various human forensic casework. The usage of variability of complete mtDNA genomes increases their discriminatory power to the maximum and enables ultimate resolution of distinct maternal lineages. However, their wider employment in forensic casework is nowadays limited by the lack of appropriate reference database. In order to fill in the gap in the reference data, which, considering Slavic-speaking populations, currently comprises only mitogenomes of East and West Slavs, we present mitogenome data for 226 Serbians, representatives of South Slavs from the Balkan Peninsula. We found 143 (sub)haplogroups among which West Eurasian ones were dominant. The percentage of unique haplotypes was 85%, and the random match probability was as low as 0.53%. We support previous findings on both high levels of genetic diversity in the Serbian population and patterns of genetic differentiation among this and ten studied European populations. However, our high-resolution data supported more pronounced genetic differentiation among Serbians and two Slavic populations (Russians and Poles) as well as expansion of the Serbian population after the Last Glacial Maximum and during the Migration period (fourth to ninth century A.D.), as inferred from the Bayesian skyline analysis. Phylogenetic analysis of haplotypes found in Serbians contributed towards the improvement of the worldwide mtDNA phylogeny, which is essential for the interpretation of the mtDNA casework.

Development and optimization of the VISAGE basic prototype tool for forensic age estimation

The VISAGE (VISible Attributes through GEnomics) consortium aims to develop, optimize and validate prototype tools to broaden the use of DNA intelligence methods in forensic routine laboratories. This includes age estimation based on the quantification of DNA methylation at specific CpG sites. Here, we present the VISAGE basic prototype tool for age estimation targeting 32 CpGs from five genes ELOVL2, MIR29B2CHG (herein, MIR29B2C), FHL2, TRIM59 and KLF14. The assay interrogates these well described age markers by multiplex PCR for bisulfite converted DNA and massively parallel sequencing on a MiSeq FGx instrument. We describe protocol optimizations including tests on five bisulfite conversion kits and an evaluation of the assay's reproducibility and sensitivity with artificially methylated DNA standards. We observed robust quantification of methylation levels with a mean standard deviation of 1.4 % across ratios. Sensitivity tests showed no increase of variability down to 20 ng DNA input into bisulfite conversion with a median difference below 1.6 % between technical replicates.

Genetic identification of a war-evacuated child in search of his own identity for more than seventy years

V. M. E. was evacuated when he was a young boy in 1939. He left an aunt and cousins in Spain (G. E. family). He was adopted in Belgium by the D. family and thus his new name became V. D. He has been unable to remember his childhood before his adoption, a symptomatology compatible with amnesia for personal identity, presumably because he may have suffered a head contusion before or during his exodus. Identification tests were performed on blood samples from V. D. and V. G. E., a mitochondrial cousin of the missing boy. V. G. E. and the missing boy have a common mitochondrial ancestor, their maternal grandmother. The mitochondrial profile of both samples turned out to be highly specific, which allowed the genetic identification of V. D. as V. M. E. As a result, V. D. has reclaimed his past and reunited with his former family in Spain after more than seven decades. As far as we know, this is the first report describing the application of mitochondrial DNA in the identification of a person evacuated during the Spanish Civil War suffering from amnesia for personal identity.

The genus Hebeloma in the Rocky Mountain Alpine Zone

Numerous taxa of Hebeloma have been reported in association with Salix, Dryas, and Betula in arctic-alpine habitats. However, species are notoriously difficult to delineate because morphological features overlap, and previously there was little reliable molecular data available. Recent progress in ITS-sequencing within the genus, coupled with an extensive database of parametrically described collections, now allows comparisons between species and their distributions. Here we report 16 species of Hebeloma from the Rocky Mountain alpine zone from some of the lowest latitudes (latitude 36°–45°N) and highest elevations (3000–4000 m) for arctic-alpine fungi in the northern hemisphere. Twelve of these species have been reported from arctic-alpine habitats in Europe and Greenland and are now molecularly confirmed from the Middle and Southern Rockies, greatly expanding their distribution. These are: Hebelomaalpinum, H.aurantioumbrinum, H.dunense, H.hiemale, H.marginatulum, H.mesophaeum, H.nigellum, H.oreophilum, H.subconcolor, H.spetsbergense, H.vaccinum, and H.velutipes. Hebelomahygrophilum is known from subalpine habitats in Europe, but was never recorded in arctic-alpine ecology. Three species recorded from the Rockies, but as yet not reported from Europe, are H.alpinicola, H.avellaneum, and H.excedens. The last two have never previously been reported from an arctic-alpine habitat. For all three of these species, the holotypes have been studied morphologically and molecularly, and have been incorporated into the analysis.

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

Haplotype-specific extraction (HSE) is a new field of application for the separation of mitochondrial DNA (mtDNA) mixtures and is developed to identify the mtDNA haplotypes of the contributors subsequently by sequencing. Here we show the validation of HSE with an exemplary mitochondrial DNA mixture into its individual haplotypes according to our laboratory standards. These specify several critical areas of assay performance to be tested, such as sensitivity, robustness and mixture studies comprising varying proportions of their components,degraded samples and samples of different qualities and material. Wereport the successful and unambiguous analysis of the exemplary separated mitochondrial DNA mixture under various conditions as well as simulated casework samples, which manifest as mixed nucleotide calls at single base positions previously. Here we demonstrate that the HSE assay is high sensitive, stable against degradation and applicable in a wide range of sample qualities. Based on our findings from the validation study, we believe that this assay has great potential power and may be useful for distinguishing among the mtDNA of individuals and their geographical origin in mixed DNA samples.

Defining mtDNA origins and population stratification in Rio de Janeiro

The genetic composition of the Brazilian population was shaped by interethnic admixture between autochthonous Native Americans, Europeans settlers and African slaves. This structure, characteristic of most American populations, implies the need for large population forensic databases to capture the high diversity that is usually associated with admixed populations. In the present work, we sequenced the control region of mitochondrial DNA from 205 non-related individuals living in the Rio de Janeiro metropolitan region. Overall high haplotype diversity (0.9994 ± 0.0006) was observed, and pairwise comparisons showed a high proportion of haplotype pairs with more than one-point differences. When ignoring homopolymeric tracts, pairwise comparisons showed no differences 0.18% of the time, and differences in a single position were found with a frequency of 0.32%. A high percentage of African mtDNA was found (42%), with lineages showing a major South West origin. For the West Eurasian and Native American haplogroups (representing 32% and 26%, respectively) it was not possible to evaluate a clear geographic or linguistic affiliation. When grouping the mtDNA lineages according to their continental origin (Native American, European and African), differences were observed for the ancestry proportions estimated with autosomal ancestry-informative markers, suggesting some level of genetic substructure. The results from this study are in accordance with historical data where admixture processes are confirmed with a strong maternal contribution of African maternal ancestry and a relevant contribution of Native American maternal ancestry. Moreover, the evidence for some degree of association between mtDNA and autosomal information should be considered when combining these types of markers in forensic analysis.

Detection of mitochondrial haplogroups in a small avar-slavic population from the eigth-ninth century AD

OBJECTIVES

In the sixth century AD, Avars came to Central Europe from middle Eurasian steppes and founded a strong Empire called the Avar Khagante (568-799/803 AD) in the Pannonian basin. During the existence of this empire, they undertook many military and pugnacious campaigns. In the seventh century, they conquered the northern territory inhabited by Slavs, who were further recruited in Avar military and were commissioned with obtaining food supplies. During almost 200 years of Avar domination, a significant influence by the Avar culture (especially on the burial rite) and assimilation with indigenous population (occurrence of "East Asian"cranial features) could be noticed in this mixed area, which is supported by achaeological and anthropologcal research. Therefore we expected higher incidence of east Eurasian haplogroups (introduced by Avars) than the frequencies detected in present-day central European populations.

MATERIALS AND METHODS

Mitochondrial DNA from 62 human skeletal remains excavated from the Avar-Slavic burial site Cífer-Pác (Slovakia) dated to the eighth and ninth century was analyzed by the sequencing of hypervariable region I and selected parts of coding region. Obtained haplotypes were compared with other present-day and historical populations and genetic distances were calculated using standard statistical method.

RESULTS AND DISCUSSION

In total, the detection of mitochondrial haplogroups was possible in 46 individuals. Our results prooved a higher frequency of east Eurasian haplogroups in our analyzed population (6.52%) than in present-day central European populations. However, it is almost three times lower than the frequency of east Eurasian haplogroups detected in other medieval Avar populations. The statistical analysis showed a greater similarity and the lowest genetic distances between the Avar-Slavic burial site Cifer-Pac and medieval European populations than the South Siberian, East and Central Asian populations.

CONCLUSION

Our results indicate that the transfer of Avar genetic variation through their mtDNA was rather weak in the analyzed mixed population.

Evaluation of GeneMarker® HTS for improved alignment of mtDNA MPS data, haplotype determination, and heteroplasmy assessment

Existing software has not allowed for effective alignment of mitochondrial (mt) DNA sequence data generated using a massively parallel sequencing (MPS) approach, combined with the ability to perform a detailed assessment of the data. The regions of sequence that are typically difficult to align are homopolymeric stretches, isolated patterns of SNPs (single nucleotide polymorphisms), and INDELs (insertions/deletions). A custom software solution, GeneMarker^® HTS, was developed and evaluated to address these limitations, and to provide a user-friendly interface for forensic practitioners and others interested in mtDNA analysis of MPS data. GeneMarker^® HTS generates an exportable consensus mtDNA sequence that produces phylogenetically correct SNP and INDEL calls using a customizable motif-based alignment algorithm. Sequence data from 500 individuals, with various alignment asymmetries and levels of heteroplasmy, were used to assess the software. Accuracy in producing mtDNA haplotypes, the ability to correctly identify low-level heteroplasmic sequence variants, and the user-based features of the software were evaluated. Analyzed sequences yielded correct mtDNA haplotypes, and heteroplasmic variants were properly identified with minimal manual interpretation. The software offers numerous user-defined parameters for filtering the data that address the interests of researchers and practitioners, and provides multiple options for viewing and navigating through the data.

High-quality mtDNA control region sequences from 680 individuals sampled across the Netherlands to establish a national forensic mtDNA reference database

The use of mitochondrial DNA (mtDNA) for maternal lineage identification often marks the last resort when investigating forensic and missing-person cases involving highly degraded biological materials. As with all comparative DNA testing, a match between evidence and reference sample requires a statistical interpretation, for which high-quality mtDNA population frequency data are crucial. Here, we determined, under high quality standards, the complete mtDNA control-region sequences of 680 individuals from across the Netherlands sampled at 54 sites, covering the entire country with 10 geographic sub-regions. The complete mtDNA control region (nucleotide positions 16,024-16,569 and 1-576) was amplified with two PCR primers and sequenced with ten different sequencing primers using the EMPOP protocol. Haplotype diversity of the entire sample set was very high at 99.63% and, accordingly, the random-match probability was 0.37%. No population substructure within the Netherlands was detected with our dataset. Phylogenetic analyses were performed to determine mtDNA haplogroups. Inclusion of these high-quality data in the EMPOP database (accession number: EMP00666) will improve its overall data content and geographic coverage in the interest of all EMPOP users worldwide. Moreover, this dataset will serve as (the start of) a national reference database for mtDNA applications in forensic and missing person casework in the Netherlands.

Simultaneous Whole Mitochondrial Genome Sequencing with Short Overlapping Amplicons Suitable for Degraded DNA Using the Ion Torrent Personal Genome Machine

Whole mitochondrial (mt) genome analysis enables a considerable increase in analysis throughput, and improves the discriminatory power to the maximum possible phylogenetic resolution. Most established protocols on the different massively parallel sequencing (MPS) platforms, however, invariably involve the PCR amplification of large fragments, typically several kilobases in size, which may fail due to mtDNA fragmentation in the available degraded materials. We introduce a MPS tiling approach for simultaneous whole human mt genome sequencing using 161 short overlapping amplicons (average 200 bp) with the Ion Torrent Personal Genome Machine. We illustrate the performance of this new method by sequencing 20 DNA samples belonging to different worldwide mtDNA haplogroups. Additional quality control, particularly regarding the potential detection of nuclear insertions of mtDNA (NUMTs), was performed by comparative MPS analysis using the conventional long-range amplification method. Preliminary sensitivity testing revealed that detailed haplogroup inference was feasible with 100 pg genomic input DNA. Complete mt genome coverage was achieved from DNA samples experimentally degraded down to genomic fragment sizes of about 220 bp, and up to 90% coverage from naturally degraded samples. Overall, we introduce a new approach for whole mt genome MPS analysis from degraded and nondegraded materials relevant to resolve and infer maternal genetic ancestry at complete resolution in anthropological, evolutionary, medical, and forensic applications.

Demographic History of Indigenous Populations in Mesoamerica Based on mtDNA Sequence Data

The genetic characterization of Native American groups provides insights into their history and demographic events. We sequenced the mitochondrial D-loop region (control region) of 520 samples from eight Mexican indigenous groups. In addition to an analysis of the genetic diversity, structure and genetic relationship between 28 Native American populations, we applied Bayesian skyline methodology for a deeper insight into the history of Mesoamerica. AMOVA tests applying cultural, linguistic and geographic criteria were performed. MDS plots showed a central cluster of Oaxaca and Maya populations, whereas those from the North and West were located on the periphery. Demographic reconstruction indicates higher values of the effective number of breeding females (Nef) in Central Mesoamerica during the Preclassic period, whereas this pattern moves toward the Classic period for groups in the North and West. Conversely, Nef minimum values are distributed either in the Lithic period (i.e. founder effects) or in recent periods (i.e. population declines). The Mesomerican regions showed differences in population fluctuation as indicated by the maximum Inter-Generational Rate (IGRmax): i) Center-South from the lithic period until the Preclassic; ii) West from the beginning of the Preclassic period until early Classic; iii) North characterized by a wide range of temporal variation from the Lithic to the Preclassic. Our findings are consistent with the genetic variations observed between central, South and Southeast Mesoamerica and the North-West region that are related to differences in genetic drift, structure, and temporal survival strategies (agriculture versus hunter-gathering, respectively). Interestingly, although the European contact had a major negative demographic impact, we detect a previous decline in Mesoamerica that had begun a few hundred years before.

Traces of forgotten historical events in mountain communities in Central Italy: A genetic insight

OBJECTIVES

Analysis of human genetic variation in mountain communities can shed light on the peopling of mountainous regions, perhaps revealing whether the remote geographic location spared them from outside invasion and preserved their gene pool from admixture. In this study, we created a model to assess genetic traces of historical events by reconstructing the paternal and maternal genetic history of seven small mountain villages in inland valleys of Central Italy.

METHODS

The communities were selected for their geographic isolation, attested biodemographic stability, and documented history prior to the Roman conquest. We studied the genetic structure by analyzing two hypervariable segments (HVS-I and HVS-II) of the mtDNA D-loop and several informative single nucleotide polymorphisms (SNPs) of the mtDNA coding region in 346 individuals, in addition to 17 short tandem repeats (STRs) and Y-chromosome SNPs in 237 male individuals.

RESULTS

For both uniparental markers, most of the haplogroups originated in Western Europe while some Near Eastern haplogroups were identified at low frequencies. However, there was an evident genetic similarity between the Central Italian samples and Near Eastern populations mainly in the male genetic pool.

CONCLUSIONS

The samples highlight an overall European genetic pattern both for mtDNA and Y chromosome. Notwithstanding this scenario, Y chromosome haplogroup Q, a common paternal lineage in Central/Western Asia but almost Europe-wide absent, was found, suggesting that Central Italy could have hosted a settlement from Anatolia that might be supported by cultural, topographic and genetic evidence.

Evidence for frequent and tissue-specific sequence heteroplasmy in human mitochondrial DNA

Mitochondrial point heteroplasmy is a common event observed not only in patients with mitochondrial diseases but also in healthy individuals. We here report a comprehensive investigation of heteroplasmy occurrence in human including the whole mitochondrial control region from nine different tissue types of 100 individuals. Sanger sequencing was used as a standard method and results were supported by cloning, minisequencing, and massively parallel sequencing. Only 12% of all individuals showed no heteroplasmy, whereas 88% showed at least one heteroplasmic position within the investigated tissues. In 66% of individuals up to 8 positions were affected. The highest relative number of heteroplasmies was detected in muscle and liver (79%, 69%), followed by brain, hair, and heart (36.7%-30.2%). Lower percentages were observed in bone, blood, lung, and buccal cells (19.8%-16.2%). Accumulation of position-specific heteroplasmies was found in muscle (positions 64, 72, 73, 189, and 408), liver (position 72) and brain (partial deletion at position 71). Deeper analysis of these specific positions in muscle revealed a non-random appearance and position-specific dependency on age. MtDNA heteroplasmy frequency and its potential functional importance have been underestimated in the past and its occurrence is ubiquitous and dependent at least on age, tissue, and position-specific mutation rates.

Massively parallel pyrosequencing of the mitochondrial genome with the 454 methodology in forensic genetics

RESULTS

of sequencing of whole mitochondrial genome, HV1 and HV2 DNA with the second generation system (SGS) Roche 454 GS Junior were compared with results of Sanger sequencing and SNP typing with SNaPshot single base extension detected with MALDI-TOF and capillary electrophoresis. We investigated the performance of the software analysis of the data, reproducibility, ability to sequence homopolymeric regions, detection of mixtures and heteroplasmy as well as the implications of the depth of coverage. We found full reproducibility between samples sequenced twice with SGS. We found close to full concordance between the mtDNA sequences of 26 samples obtained with (1) the 454 SGS method using a depth of coverage above 100 and (2) Sanger sequencing and SNP typing. The discrepancies were primarily observed in homopolymeric regions. The 454 SGS method was able to sequence 95% of the reads correctly in homopolymers up to 4 bases, and up to 6 bases could be sequenced with similar success if the results were carefully, visually inspected. The 454 technology was able to detect mixtures or heteroplasmy of approximately 10%. We detected previously unreported heteroplasmy in the GM9947A component of the NIST human mitochondrial DNA SRM-2392 standard reference material.

Maternal genetic heritage of Southeastern Europe reveals a new Croatian isolate and a novel, local sub-branching in the x2 haplogroup

High mtDNA variation in Southeastern Europe (SEE) is a reflection of the turbulent and complex demographic history of this area, influenced by gene flow from various parts of Eurasia and a long history of intermixing. Our results of 1035 samples (488 from Croatia, 239 from Bosnia and 130 from Herzegovina, reported earlier, and 97 Slovenians and 81 individuals from Žumberak, reported here for the first time) show that the SEE maternal genetic diversity fits within a broader European maternal genetic landscape. The study also shows that the population of Žumberak, located in the continental part of Croatia, developed some unique mtDNA haplotypes and elevated haplogroup frequencies due to distinctive demographic history and can be considered a moderate genetic isolate. We also report seven samples from the Bosnian population and one Herzegovinian sample designated as X2* individuals that could not be assigned to any of its sublineages (X2a'o) according to the existing X2 phylogeny. In an attempt to clarify the phylogeny of our X2 samples, their mitochondrial DNA has been completely sequenced. We suppose that these lineages are signs of local microdifferentiation processes that occurred in the recent demographic past in this area and could possibly be marked as SEE-specific X2 sublineages.

Large-scale mitochondrial DNA analysis in Southeast Asia reveals evolutionary effects of cultural isolation in the multi-ethnic population of Myanmar

Background

Myanmar is the largest country in mainland Southeast Asia with a population of 55 million people subdivided into more than 100 ethnic groups. Ruled by changing kingdoms and dynasties and lying on the trade route between India and China, Myanmar was influenced by numerous cultures. Since its independence from British occupation, tensions between the ruling Bamar and ethnic minorities increased.

Results

Our aim was to search for genetic footprints of Myanmar’s geographic, historic and sociocultural characteristics and to contribute to the picture of human colonization by describing and dating of new mitochondrial DNA (mtDNA) haplogroups. Therefore, we sequenced the mtDNA control region of 327 unrelated donors and the complete mitochondrial genome of 44 selected individuals according to highest quality standards.

Conclusion

Phylogenetic analyses of the entire mtDNA genomes uncovered eight new haplogroups and three unclassified basal M-lineages. The multi-ethnic population and the complex history of Myanmar were reflected in its mtDNA heterogeneity. Population genetic analyses of Burmese control region sequences combined with population data from neighboring countries revealed that the Myanmar haplogroup distribution showed a typical Southeast Asian pattern, but also Northeast Asian and Indian influences. The population structure of the extraordinarily diverse Bamar differed from that of the Karen people who displayed signs of genetic isolation. Migration analyses indicated a considerable genetic exchange with an overall positive migration balance from Myanmar to neighboring countries. Age estimates of the newly described haplogroups point to the existence of evolutionary windows where climatic and cultural changes gave rise to mitochondrial haplogroup diversification in Asia.

Reviewing population studies for forensic purposes: Dog mitochondrial DNA

The identification of dog hair through mtDNA analysis has become increasingly important in the last 15 years, as it can provide associative evidence connecting victims and suspects. The evidential value of an mtDNA match between dog hair and its potential donor is determined by the random match probability of the haplotype. This probability is based on the haplotype’s population frequency estimate. Consequently, implementing a population study representative of the population relevant to the forensic case is vital to the correct evaluation of the evidence. This paper reviews numerous published dog mtDNA studies and shows that many of these studies vary widely in sampling strategies and data quality. Therefore, several features influencing the representativeness of a population sample are discussed. Moreover, recommendations are provided on how to set up a dog mtDNA population study and how to decide whether or not to include published data. This review emphasizes the need for improved dog mtDNA population data for forensic purposes, including targeting the entire mitochondrial genome. In particular, the creation of a publicly available database of qualitative dog mtDNA population studies would improve the genetic analysis of dog traces in forensic casework.

Association of polymorphisms in pharmacogenetic candidate genes (OPRD1, GAL, ABCB1, OPRM1) with opioid dependence in European population: a case-control study

It is becoming increasingly evident that genetic variants contribute to the development of opioid addiction. An elucidation of these genetic factors is crucial for a better understanding of this chronic disease and may help to develop novel therapeutic strategies. In recent years, several candidate genes were implicated in opioid dependence. However, most study findings have not been replicated and additional studies are required before reported associations can be considered robust. Thus, the major objective of this study was to replicate earlier findings and to identify new genetic polymorphisms contributing to the individual susceptibility to opioid addiction, respectively. Therefore, a candidate gene association study was conducted including 142 well-phenotyped long-term opioid addicts undergoing opioid maintenance therapy and 142 well-matched healthy controls. In both study groups, 24 single nucleotide polymorphisms predominantly located in pharmacogenetic candidate genes have been genotyped using an accurate mass spectrometry based method. The most significant associations with opioid addiction (remaining significant after adjustment for multiple testing) were observed for the rs948854 SNP in the galanin gene (GAL, p = 0.001) and the rs2236861 SNP in the delta opioid receptor gene (OPRD1, p = 0.001). Moreover, an association of the ATP binding cassette transporter 1 (ABCB1) variant rs1045642 and the Mu Opioid receptor (OPRM1) variant rs9479757 with opioid addiction was observed. The present study provides further support for a contribution of GAL and OPRD1 variants to the development of opioid addiction. Furthermore, our results indicate a potential contribution of OPRM1 and ABCB1 SNPs to the development of this chronic relapsing disease. Therefore it seems important that these genes are addressed in further addiction related studies.

High resolution mapping of Y haplogroup G in Tyrol (Austria)

The distribution of Y-chromosomal haplogroup G2a (G-P15) in present-day paternal lineages in Tyrol (Austria) was analyzed by applying a high-density regional sampling scheme that also covered remote mountain areas. There is evidence from ancient genetic data for a high frequency of Y-chromosomal haplogroup G in prehistoric populations of Central Europe, whilst nowadays levels well below 10% are routinely observed. A population sample comprising ∼3700 specimens was analyzed for Y-chromosomal variation by genotyping Y-SNPs and Y-STRs. The set of binary markers included nine SNPs specific for sub-lineages of haplogroup G. The frequency of haplogroup G in 2379 unrelated men born in Tyrol amounted to 11.3%. Nearly all of these Y chromosomes belonged to haplogroup G2a. The main sub-haplogroup within G2a was defined by the SNP L497 (G2a3b1c) and reached a population frequency of 8.6%. Although this average level is higher than reported for other countries the geographical distribution of haplogroup G-L497 showed a differentiated pattern with a clustered distribution within some alpine valleys, where maxima above 40% were found. Both, the estimation of coalescent times and a principle coordinates analysis based on RST values derived from Y-STR haplotypes from different sub-regions of Tyrol revealed evidence for an old settlement history associated with Y chromosomes belonging to haplogroup G in the Tyrolean Alps.

A novel amplification strategy for genotyping with liquid chromatography-electrospray ionization mass spectrometry

Among numerous available genotyping techniques, mass spectrometry (MS) based methods play a major role in providing high quality genotype data at reasonable costs for research and diagnostics, e.g. for pharmacogenetic applications. Ion-pair reversed-phase liquid chromatography hyphenated to electrospray ionization time-of-flight MS (ICEMS) is, for example, a powerful instrument that allows a direct characterization of complex mixtures of polymerase chain reaction (PCR) amplified DNA fragments. Current limitations of PCR-ICEMS genotyping are mainly concerned with the multiplex PCR set-up. Assay development often requires time-consuming primer design and intensive optimization of PCR conditions. To overcome this restraint, a robust amplification strategy originally combined with arrayed primer extension genotyping was transferred and adapted to ICEMS genotyping. The modifications involved limitation of the primer length, application of two universal sequences and amplification with an appropriate DNA polymerase. To demonstrate the applicability of the novel amplification strategy for ICEMS, a 23-plex pharmacogenetic genotyping assay was developed. After slight optimization steps, an efficient and quantitatively balanced amplification of all targeted markers was achieved, resulting in a convenient characterization of the multiplexed PCR fragments with ICEMS. Expenditure of time, costs and hands-on work associated with assay design and optimization was dramatically lowered compared to previous multiplex PCR-ICEMS assays. The developed 23-plex assay was applied in a pharmacogenetic study including 284 individuals (genotype call rate 99.0%). A total of 399 SNPs were retyped by Sanger sequencing (concordance rate 99.8%). The PCR-ICEMS assay turned out to be an accurate, reliable, cost-effective and a ready-to-use tool for pharmacogenetic genotyping.

Polymorphisms in the gene regions of the adaptor complex LAMTOR2/LAMTOR3 and their association with breast cancer risk

Background

The late endosomal LAMTOR complex serves as a convergence point for both the RAF/MEK/ERK and the PI3K/AKT/mTOR pathways. Interestingly, both of these signalling cascades play a significant role in the aetiology of breast cancer. Our aim was to address the possible role of genetic polymorphisms in LAMTOR2 and LAMTOR3 as genetic risk factors for breast cancer.

Methodology/Results

We sequenced the exons and exon–intron boundaries of LAMTOR2 (p14) and LAMTOR3 (MP1) in 50 prospectively collected pairs of cancerous tissue and blood samples from breast cancer patients and compared their genetic variability. We found one single nucleotide polymorphism (SNP) in LAMTOR2 (rs7541) and two SNPs in LAMTOR3 (rs2298735 and rs148972953) in both tumour and blood samples, but no somatic mutations in cancerous tissues. In addition, we genotyped all three SNPs in 296 samples from the Risk Prediction of Breast Cancer Metastasis Study and found evidence of a genetic association between rs148972953 and oestrogen (ER) and progesterone receptor negative status (PR) (ER: OR = 3.60 (1.15–11.28); PR: OR = 4.27 (1.43–12.72)). However, when we additionally genotyped rs148972953 in the MARIE study including 2,715 breast cancer cases and 5,216 controls, we observed neither a difference in genotype frequencies between patients and controls nor was the SNP associated with ER or PR. Finally, all three SNPs were equally frequent in breast cancer samples and female participants (n = 640) of the population-based SAPHIR Study.

Conclusions

The identified polymorphisms in LAMTOR2 and LAMTOR3 do not seem to play a relevant role in breast cancer. Our work does not exclude a role of other not yet identified SNPs or that the here annotated polymorphism may in fact play a relevant role in other diseases. Our results underscore the importance of replication in association studies.

Assessing the genetic influence of ancient sociopolitical structure: micro-differentiation patterns in the population of Asturias (Northern Spain)

The human populations of the Iberian Peninsula are the varied result of a complex mixture of cultures throughout history, and are separated by clear social, cultural, linguistic or geographic barriers. The stronger genetic differences between closely related populations occur in the northern third of Spain, a phenomenon commonly known as "micro-differentiation". It has been argued and discussed how this form of genetic structuring can be related to both the rugged landscape and the ancient societies of Northern Iberia, but this is difficult to test in most regions due to the intense human mobility of previous centuries. Nevertheless, the Spanish autonomous community of Asturias shows a complex history which hints of a certain isolation of its population. This, joined together with a difficult terrain full of deep valleys and steep mountains, makes it suitable for performing a study of genetic structure, based on mitochondrial DNA and Y-Chromosome markers. Our analyses do not only show that there are micro-differentiation patterns inside the Asturian territory, but that these patterns are strikingly similar between both uniparental markers. The inference of barriers to gene flow also indicates that Asturian populations from the coastal north and the mountainous south seem to be relatively isolated from the rest of the territory. These findings are discussed in light of historic and geographic data and, coupled with previous evidence, show that the origin of the current genetic patterning might indeed lie in Roman and Pre-Roman sociopolitical divisions.

Mitochondrial diversity patterns and the Magdalenian resettlement of Europe: new insights from the edge of the Franco-Cantabrian refuge

Phylogeography of the mitochondrial lineages commonly found in Western Europe can be interpreted in the light of a postglacial resettlement of the continent. The center of this proposal lies in the Franco-Cantabrian glacial refuge, located in the northern Iberian Peninsula and Southwestern France. Recently, this interpretation has been confronted by the unexpected patterns of diversity found in some European haplogroups. To shed new lights on this issue, research on Iberian populations is crucial if events behind the actual genetics of the European continent are to be untangled. In this regard, the region of Asturias has not been extensively studied, despite its convoluted history with prolonged periods of isolation. As mitochondrial DNA is a kind of data that has been commonly used in human population genetics, we conducted a thorough regional study in which we collected buccal swabs from 429 individuals with confirmed Asturian ancestry. The joint analysis of these sequences with a large continent-wide database and previously published diversity patterns allowed us to discuss a new explanation for the population dynamics inside the Franco-Cantabrian area, based on range expansion theory. This approximation to previously contradictory findings has made them compatible with most proposals about the postglacial resettlement of Western Europe.

Pasture names with Romance and Slavic roots facilitate dissection of Y chromosome variation in an exclusively German-speaking alpine region

The small alpine district of East Tyrol (Austria) has an exceptional demographic history. It was contemporaneously inhabited by members of the Romance, the Slavic and the Germanic language groups for centuries. Since the Late Middle Ages, however, the population of the principally agrarian-oriented area is solely Germanic speaking. Historic facts about East Tyrol's colonization are rare, but spatial density-distribution analysis based on the etymology of place-names has facilitated accurate spatial mapping of the various language groups' former settlement regions. To test for present-day Y chromosome population substructure, molecular genetic data were compared to the information attained by the linguistic analysis of pasture names. The linguistic data were used for subdividing East Tyrol into two regions of former Romance (A) and Slavic (B) settlement. Samples from 270 East Tyrolean men were genotyped for 17 Y-chromosomal microsatellites (Y-STRs) and 27 single nucleotide polymorphisms (Y-SNPs). Analysis of the probands' surnames revealed no evidence for spatial genetic structuring. Also, spatial autocorrelation analysis did not indicate significant correlation between genetic (Y-STR haplotypes) and geographic distance. Haplogroup R-M17 chromosomes, however, were absent in region A, but constituted one of the most frequent haplogroups in region B. The R-M343 (R1b) clade showed a marked and complementary frequency distribution pattern in these two regions. To further test East Tyrol's modern Y-chromosomal landscape for geographic patterning attributable to the early history of settlement in this alpine area, principal coordinates analysis was performed. The Y-STR haplotypes from region A clearly clustered with those of Romance reference populations and the samples from region B matched best with Germanic speaking reference populations. The combined use of onomastic and molecular genetic data revealed and mapped the marked structuring of the distribution of Y chromosomes in an alpine region that has been culturally homogeneous for centuries.

Influence of evolutionary forces and demographic processes on the genetic structure of three Croatian populations: a maternal perspective

BACKGROUND

Many Croatian islands are examples of genetic isolates, with low level of heterozygosity and high level of inbreeding, due to practice of endogamy.

AIM

The aim was to study the genetic structure of two insular and one mainland population through high-resolution phylogenetic analysis of mitochondrial DNA (mtDNA).

SUBJECTS AND METHODS

MtDNA polymorphisms were explored in 300 unrelated individuals from Mljet, Lastovo and the coastal city of Dubrovnik, based on SNP polymorphisms.

RESULTS

All mtDNA haplogroups found in the sample were of typical European origin. However, the frequency distribution of their subclades differed significantly from other Croatian and European populations. MtDNA haplotype analysis revealed only two possible founder lineages on Mljet and six on Lastovo, accounting for almost half of the sample on both islands. The island of Mljet also has the lowest reported haplotype and nucleotide diversity among Croatian isolates and the island of Lastovo, a new sublineage of a usually quite rare U1b clade.

CONCLUSION

The results can be explained by the effect evolutionary forces have on genetic structure, which is in line with the specific demographic histories of the islands. An additional research value of these two island isolates is the appearance of certain Mendelian disorders, highlighting their importance in epidemiological studies.

Multiplex genotyping system for efficient inference of matrilineal genetic ancestry with continental resolution

Background

In recent years, phylogeographic studies have produced detailed knowledge on the worldwide distribution of mitochondrial DNA (mtDNA) variants, linking specific clades of the mtDNA phylogeny with certain geographic areas. However, a multiplex genotyping system for the detection of the mtDNA haplogroups of major continental distribution that would be desirable for efficient DNA-based bio-geographic ancestry testing in various applications is still missing.

Results

Three multiplex genotyping assays, based on single-base primer extension technology, were developed targeting a total of 36 coding-region mtDNA variants that together differentiate 43 matrilineal haplo-/paragroups. These include the major diagnostic haplogroups for Africa, Western Eurasia, Eastern Eurasia and Native America. The assays show high sensitivity with respect to the amount of template DNA: successful amplification could still be obtained when using as little as 4 pg of genomic DNA and the technology is suitable for medium-throughput analyses.

Conclusions

We introduce an efficient and sensitive multiplex genotyping system for bio-geographic ancestry inference from mtDNA that provides resolution on the continental level. The method can be applied in forensics, to aid tracing unknown suspects, as well as in population studies, genealogy and personal ancestry testing. For more complete inferences of overall bio-geographic ancestry from DNA, the mtDNA system provided here can be combined with multiplex systems for suitable autosomal and, in the case of males, Y-chromosomal ancestry-sensitive DNA markers.

Somatic mutations throughout the entire mitochondrial genome are associated with elevated PSA levels in prostate cancer patients

The genetic etiology of prostate cancer, the most common form of male cancer in western countries, is complex and the interplay of disease genes with environmental factors is far from being understood. Studies on somatic mitochondrial DNA (mtDNA) mutations have become an important aspect of cancer research because these mutations might have functional consequences and/or might serve as biosensors for tumor detection and progression. We sequenced the entire mitochondrial genome (16,569 bp) from 30 prospectively collected pairs of macrodissected cancerous and benign cells from prostate cancer patients and compared their genetic variability. Given recent concerns regarding the authenticity of newly discovered mtDNA mutations, we implemented a high-quality procedure for mtDNA whole-genome sequencing. In addition, the mitochondrial genes MT-CO2, MT-CO3, MT-ATP6, and MT-ND6 were sequenced in further 35 paired samples from prostate cancer patients. We identified a total of 41 somatic mutations in 22 out of 30 patients: the majority of these mutations have not previously been observed in the human phylogeny. The presence of somatic mutations in transfer RNAs (tRNAs) was found to be associated with elevated PSA levels (14.25 ± 5.44 versus 7.15 ± 4.32 ng/ml; p = 0.004). The level and degree of heteroplasmy increased with increasing tumor activity. In summary, somatic mutations in the mitochondrial genome are frequent events in prostate cancer. Mutations mapping to mitochondrial tRNAs, ribosomal RNAs, and protein coding genes might impair processes that occur within the mitochondrial compartment (e.g., transcription, RNA processing, and translation) and might finally affect oxidative phosphorylation.

HaploGrep: a fast and reliable algorithm for automatic classification of mitochondrial DNA haplogroups

An ongoing source of controversy in mitochondrial DNA (mtDNA) research is based on the detection of numerous errors in mtDNA profiles that led to erroneous conclusions and false disease associations. Most of these controversies could be avoided if the samples' haplogroup status would be taken into consideration. Knowing the mtDNA haplogroup affiliation is a critical prerequisite for studying mechanisms of human evolution and discovering genes involved in complex diseases, and validating phylogenetic consistency using haplogroup classification is an important step in quality control. However, despite the availability of Phylotree, a regularly updated classification tree of global mtDNA variation, the process of haplogroup classification is still time-consuming and error-prone, as researchers have to manually compare the polymorphisms found in a population sample to those summarized in Phylotree, polymorphism by polymorphism, sample by sample. We present HaploGrep, a fast, reliable and straight-forward algorithm implemented in a Web application to determine the haplogroup affiliation of thousands of mtDNA profiles genotyped for the entire mtDNA or any part of it. HaploGrep uses the latest version of Phylotree and offers an all-in-one solution for quality assessment of mtDNA profiles in clinical genetics, population genetics and forensics. HaploGrep can be accessed freely at http://haplogrep.uibk.ac.at.

Application of a west Eurasian-specific filter for quasi-median network analysis: Sharpening the blade for mtDNA error detection

The application of quasi-median networks provides an effective tool to check the quality of mtDNA data. Filtering of highly recurrent mutations prior to network analysis is required to simplify the data set and reduce the complexity of the network. The phylogenetic background determines those mutations that need to be filtered. While the traditional EMPOPspeedy filter was based on the worldwide mtDNA phylogeny, haplogroup-specific filters can more effectively highlight potential errors in data of the respective (sub)-continental region. In this study we demonstrate the performance of a new, west Eurasian filter EMPOPspeedyWE for the fine-tuned examination of data sets belonging to macrohaplogroup N that constitutes the main portion of mtDNA lineages in Europe. The effects on the resulting network of different database sizes, high-quality and flawed data, as well as the examination of a phylogenetically distant data set, are presented by examples. The analyses are based on a west Eurasian etalon data set that was carefully compiled from more than 3500 control region sequences for network purposes. Both, etalon data and the new filter file, are provided through the EMPOP database (www.empop.org).

Inspecting close maternal relatedness: Towards better mtDNA population samples in forensic databases

Reliable data are crucial for all research fields applying mitochondrial DNA (mtDNA) as a genetic marker. Quality control measures have been introduced to ensure the highest standards in sequence data generation, validation and a posteriori inspection. A phylogenetic alignment strategy has been widely accepted as a prerequisite for data comparability and database searches, for forensic applications, for reconstructions of human migrations and for correct interpretation of mtDNA mutations in medical genetics.

There is continuing effort to enhance the number of worldwide population samples in order to contribute to a better understanding of human mtDNA variation. This has often lead to the analysis of convenience samples collected for other purposes, which might not meet the quality requirement of random sampling for mtDNA data sets. Here, we introduce an additional quality control means that deals with one aspect of this limitation: by combining autosomal short tandem repeat (STR) marker with mtDNA information, it helps to avoid the bias introduced by related individuals included in the same (small) sample. By STR analysis of individuals sharing their mitochondrial haplotype, pedigree construction and subsequent software-assisted calculation of likelihood ratios based on the allele frequencies found in the population, closely maternally related individuals can be identified and excluded. We also discuss scenarios that allow related individuals in the same set.

An ideal population sample would be representative for its population: this new approach represents another contribution towards this goal.

SAM: String-based sequence search algorithm for mitochondrial DNA database queries

The analysis of the haploid mitochondrial (mt) genome has numerous applications in forensic and population genetics, as well as in disease studies. Although mtDNA haplotypes are usually determined by sequencing, they are rarely reported as a nucleotide string. Traditionally they are presented in a difference-coded position-based format relative to the corrected version of the first sequenced mtDNA. This convention requires recommendations for standardized sequence alignment that is known to vary between scientific disciplines, even between laboratories. As a consequence, database searches that are vital for the interpretation of mtDNA data can suffer from biased results when query and database haplotypes are annotated differently. In the forensic context that would usually lead to underestimation of the absolute and relative frequencies. To address this issue we introduce SAM, a string-based search algorithm that converts query and database sequences to position-free nucleotide strings and thus eliminates the possibility that identical sequences will be missed in a database query. The mere application of a BLAST algorithm would not be a sufficient remedy as it uses a heuristic approach and does not address properties specific to mtDNA, such as phylogenetically stable but also rapidly evolving insertion and deletion events. The software presented here provides additional flexibility to incorporate phylogenetic data, site-specific mutation rates, and other biologically relevant information that would refine the interpretation of mitochondrial DNA data. The manuscript is accompanied by freeware and example data sets that can be used to evaluate the new software (http://stringvalidation.org).

Mitochondrial haplogroup N1a phylogeography, with implication to the origin of European farmers

BACKGROUND

Tracing the genetic origin of central European farmer N1a lineages can provide a unique opportunity to assess the patterns of the farming technology spread into central Europe in the human prehistory. Here, we have chosen twelve N1a samples from modern populations which are most similar with the farmer N1a types and performed the complete mitochondrial DNA genome sequencing analysis. To assess the genetic and phylogeographic relationship, we performed a detailed survey of modern published N1a types from Eurasian and African populations.

RESULTS

The geographic origin and expansion of farmer lineages related N1a subclades have been deduced from combined analysis of 19 complete sequences with 166 N1a haplotypes. The phylogeographic analysis revealed that the central European farmer lineages have originated from different sources: from eastern Europe, local central Europe, and from the Near East via southern Europe.

CONCLUSIONS

The results obtained emphasize that the arrival of central European farmer lineages did not occur via a single demic diffusion event from the Near East at the onset of the Neolithic spread of agriculture into Europe. Indeed these results indicate that the Neolithic transition process was more complex in central Europe and possibly the farmer N1a lineages were a result of a 'leapfrog' colonization process.