A collaborative exercise on DNA methylation-based age prediction and body fluid typing.

DNA methylation has become one of the most useful biomarkers for age prediction and body fluid identification in the forensic field. Therefore, several assays have been developed to detect age-associated and body fluid-specific DNA methylation changes. Among the many methods developed, SNaPshot-based assays should be particularly useful in forensic laboratories, as they permit multiplex analysis and use the same capillary electrophoresis instrumentation as STR analysis. However, technical validation of any developed assays is crucial for their proper integration into routine forensic workflow. In the present collaborative exercise, two SNaPshot multiplex assays for age prediction and a SNaPshot multiplex for body fluid identification were tested in twelve laboratories. The experimental set-up of the exercise was designed to reflect the entire workflow of SNaPshot-based methylation analysis and involved four increasingly complex tasks designed to detect potential factors influencing methylation measurements. The results of body fluid identification from each laboratory provided sufficient information to determine appropriate age prediction methods in subsequent analysis. In age prediction, systematic measurement differences resulting from the type of genetic analyzer used were identified as the biggest cause of DNA methylation variation between laboratories. Also, the use of a buffer that ensures a high ratio of specific to non-specific primer binding resulted in changes in DNA methylation measurement, especially when using degenerate primers in the PCR reaction. In addition, high input volumes of bisulfite-converted DNA often caused PCR failure, presumably due to carry-over of PCR inhibitors from the bisulfite conversion reaction. The proficiency of the analysts and experimental conditions for efficient SNaPshot reactions were also important for consistent DNA methylation measurement. Several bisulfite conversion kits were used for this study, but differences resulting from the use of any specific kit were not clearly discerned. Even when different experimental settings were used in each laboratory, a positive outcome of the study was a mean absolute age prediction error amongst participant's data of only 2.7 years for semen, 5.0 years for blood and 3.8 years for saliva.

Development and Evaluation of the Ancestry Informative Marker Panel of the VISAGE Basic Tool.

We detail the development of the ancestry informative single nucleotide polymorphisms (SNPs) panel forming part of the VISAGE Basic Tool (BT), which combines 41 appearance predictive SNPs and 112 ancestry predictive SNPs (three SNPs shared between sets) in one massively parallel sequencing (MPS) multiplex, whereas blood-based age analysis using methylation markers is run in a parallel MPS analysis pipeline. The selection of SNPs for the BT ancestry panel focused on established forensic markers that already have a proven track record of good sequencing performance in MPS, and the overall SNP multiplex scale closely matched that of existing forensic MPS assays. SNPs were chosen to differentiate individuals from the five main continental population groups of Africa, Europe, East Asia, America, and Oceania, extended to include differentiation of individuals from South Asia. From analysis of 1000 Genomes and HGDP-CEPH samples from these six population groups, the BT ancestry panel was shown to have no classification error using the Bayes likelihood calculators of the Snipper online analysis portal. The differentiation power of the component ancestry SNPs of BT was balanced as far as possible to avoid bias in the estimation of co-ancestry proportions in individuals with admixed backgrounds. The balancing process led to very similar cumulative population-specific divergence values for Africa, Europe, America, and Oceania, with East Asia being slightly below average, and South Asia an outlier from the other groups. Comparisons were made of the African, European, and Native American estimated co-ancestry proportions in the six admixed 1000 Genomes populations, using the BT ancestry panel SNPs and 572,000 Affymetrix Human Origins array SNPs. Very similar co-ancestry proportions were observed down to a minimum value of 10%, below which, low-level co-ancestry was not always reliably detected by BT SNPs. The Snipper analysis portal provides a comprehensive population dataset for the BT ancestry panel SNPs, comprising a 520-sample standardised reference dataset; 3445 additional samples from 1000 Genomes, HGDP-CEPH, Simons Foundation and Estonian Biocentre genome diversity projects; and 167 samples of six populations from in-house genotyping of individuals from Middle East, North and East African regions complementing those of the sampling regimes of the other diversity projects.

Development and validation of the VISAGE AmpliSeq basic tool to predict appearance and ancestry from DNA.

Forensic DNA phenotyping is gaining interest as the number of applications increases within the forensic genetics community. The possibility of providing investigative leads in addition to conventional DNA profiling for human identification provides new insights into otherwise "cold" police investigations. The ability of reporting on the bio-geographical ancestry (BGA), appearance characteristics and age based on DNA obtained from a crime scene sample of an unknown donor makes the exploration of such markers and the development of new methods meaningful for criminal investigations. The VISible Attributes through GEnomics (VISAGE) Consortium aims to disseminate and broaden the use of predictive markers and develop fully optimized and validated prototypes for forensic casework implementation. Here, the first VISAGE appearance and ancestry tool development, performance and validation is reported. A total of 153 SNPs (96.84 % assay conversion rate) were successfully incorporated into a single multiplex reaction using the AmpliSeq™ design pipeline, and applied for massively parallel sequencing with the Ion S5 platform. A collaborative effort involving six VISAGE laboratory partners was devised to perform all validation tests. An extensive validation plan was carefully organized to explore the assay's overall performance with optimum and low-input samples, as well as with challenging and casework mock samples. In addition, forensic validation studies such as concordance and mixture tests recurring to the Coriell sample set with known genotypes were performed. Finally, inhibitor tolerance and specificity were also evaluated. Results showed a robust, highly sensitive assay with good overall concordance between laboratories.

Towards broadening Forensic DNA Phenotyping beyond pigmentation: Improving the prediction of head hair shape from DNA.

Human head hair shape, commonly classified as straight, wavy, curly or frizzy, is an attractive target for Forensic DNA Phenotyping and other applications of human appearance prediction from DNA such as in paleogenetics. The genetic knowledge underlying head hair shape variation was recently improved by the outcome of a series of genome-wide association and replication studies in a total of 26,964 subjects, highlighting 12 loci of which 8 were novel and introducing a prediction model for Europeans based on 14 SNPs. In the present study, we evaluated the capacity of DNA-based head hair shape prediction by investigating an extended set of candidate SNP predictors and by using an independent set of samples for model validation. Prediction model building was carried out in 9674 subjects (6068 from Europe, 2899 from Asia and 707 of admixed European and Asian ancestries), used previously, by considering a novel list of 90 candidate SNPs. For model validation, genotype and phenotype data were newly collected in 2415 independent subjects (2138 Europeans and 277 non-Europeans) by applying two targeted massively parallel sequencing platforms, Ion Torrent PGM and MiSeq, or the MassARRAY platform. A binomial model was developed to predict straight vs. non-straight hair based on 32 SNPs from 26 genetic loci we identified as significantly contributing to the model. This model achieved prediction accuracies, expressed as AUC, of 0.664 in Europeans and 0.789 in non-Europeans; the statistically significant difference was explained mostly by the effect of one EDAR SNP in non-Europeans. Considering sex and age, in addition to the SNPs, slightly and insignificantly increased the prediction accuracies (AUC of 0.680 and 0.800, respectively). Based on the sample size and candidate DNA markers investigated, this study provides the most robust, validated, and accurate statistical prediction models and SNP predictor marker sets currently available for predicting head hair shape from DNA, providing the next step towards broadening Forensic DNA Phenotyping beyond pigmentation traits.

Tracking age-correlated DNA methylation markers in the young.

DNA methylation is the most extensively studied epigenetic signature, with a large number of studies reporting age-correlated CpG sites in overlapping genes. However, most of these studies lack sample coverage of individuals under 18 years old and therefore little is known about the progression of DNA methylation patterns in children and adolescents. In the present study we aimed to select candidate age-correlated DNA methylation markers based on public datasets from Illumina BeadChip arrays and previous publications, then to explore the resulting markers in 209 blood samples from donors aged between 2 to 18 years old using the EpiTYPER® DNA methylation analysis system. Results from our analyses identified six genes highly correlated with age in the young, in particular the gene KCNAB3, which indicates its potential as a highly informative and specific age biomarker for childhood and adolescence. We outline a preliminary age prediction model based on quantile regression that uses data from the six CpG sites most strongly correlated with age ranges extended to include children and adolescents.

Meta-Analysis of Mitochondrial DNA Variation in the Iberian Peninsula.

The Iberian Peninsula has been the focus of attention of numerous studies dealing with mitochondrial DNA (mtDNA) variation, most of them targeting the control region segment. In the present study we sequenced the control region of 3,024 Spanish individuals from areas where available data were still limited. We also compiled mtDNA haplotypes from the literature involving 4,588 sequences and 28 population groups or small regions. We meta-analyzed all these data in order to shed further light on patterns of geographic variation, taking advantage of the large sample size and geographic coverage, in contrast with the atomized sampling strategy of previous work. The results indicate that the main mtDNA haplogroups show primarily clinal geographic patterns across the Iberian geography, roughly along a North-South axis. Haplogroup HV0 (where haplogroup U is nested) is more prevalent in the Franco Cantabrian region, in good agreement with previous findings that identified this area as a climate refuge during the Last Glacial Maximum (LGM), prior to a subsequent demographic re-expansion towards Central Europe and the Mediterranean. Typical sub-Saharan and North African lineages are slightly more prevalent in South Iberia, although at low frequencies; this pattern has been shaped mainly by the transatlantic slave trade and the Arab invasion of the Iberian Peninsula. The results also indicate that summary statistics that aim to measure molecular variation, or AMOVA, have limited sensitivity to detect population substructure, in contrast to patterns revealed by phylogeographic analysis. Overall, the results suggest that mtDNA variation in Iberia is substantially stratified. These patterns might be relevant in biomedical studies given that stratification is a common cause of false positives in case-control mtDNA association studies, and should be also considered when weighting the DNA evidence in forensic casework, which is strongly dependent on haplotype frequencies.

Genomic insights on the ethno-history of the Maya and the 'Ladinos' from Guatemala.

BACKGROUND: Guatemala is a multiethnic and multilingual country located in Central America. The main population groups separate 'Ladinos' (mixed Native American-African-Spanish), and Native indigenous people of Maya descent. Among the present-day Guatemalan Maya, there are more than 20 different ethnic groups separated by different languages and cultures. Genetic variation of these communities still remains largely unexplored. The principal aim of this study is to explore the genetic variability of the Maya and 'Ladinos' from Guatemala by means of uniparental and ancestry informative markers (AIMs). RESULTS: Analyses of uniparental genetic markers indicate that Maya have a dominant Native American ancestry (mitochondrial DNA [mtDNA]: 100%; Y-chromosome: 94%). 'Ladino', however, show a clear gender-bias as indicated by the large European ancestry observed in the Y-chromosome (75%) compared to the mtDNA (0%). Autosomal polymorphisms (AIMS) also mirror this marked gender-bias: (i) Native American ancestry: 92% for the Maya vs. 55% for the 'Ladino', and (ii) European ancestry: 8% for the Maya vs. 41% for the 'Ladino'. In addition, the impact of the Trans-Atlantic slave trade on the present-day Guatemalan population is very low (and only occurs in the 'Ladino'; mtDNA: 9%; AIMS: 4%), in part mirroring the fact that Guatemala has a predominant orientation to the Pacific Ocean instead of a Caribbean one. Sequencing of entire Guatemalan mitogenomes has led to improved Native American phylogeny via the addition of new haplogroups that are mainly observed in Mesoamerica and/or the North of South America. CONCLUSIONS: The data reveal the existence of a fluid gene flow in the Mesoamerican area and a predominant unidirectional flow towards South America, most likely occurring during the Pre-Classic (1800 BC-200 AD) and the Classic (200-1000 AD) Eras of the Mesoamerican chronology, coinciding with development of the most distinctive and advanced Mesoamerican civilization, the Maya. Phylogenetic features of mtDNA data also suggest a demographic scenario that is compatible with moderate local endogamy and isolation in the Maya combined with episodes of gene exchange between ethnic groups, suggesting an ethno-genesis in the Guatemalan Maya that is recent and supported on a cultural rather than a biological basis.

No evidence of association between common European mitochondrial DNA variants in Alzheimer, Parkinson, and migraine in the Spanish population.

Certain mitochondrial DNA (mtDNA) variants and haplogroups have been found to be associated with neurological disorders. Several studies have suggested that mtDNA variation could have an etiologic role in these disorders by affecting the ATP production on high-energy demanding organs, such as the brain. We have analyzed 15 mtDNA SNPs (mtSNPs) in five cohorts of cases presenting Alzheimer disease (AD), Parkinson disease (PD), and migraine, and in controls, to evaluate the role mtDNA variation in disease risk. Association tests were undertaken both for mtSNPs and mitochondrial haplogroups. No significant association was detected for any mtSNP or haplogroup in AD and PD cohorts. Two mtSNPs were associated with one migraine cohort after correcting for multiple tests, namely, T4216C and G13708A and haplogroup J (FDR q-value = 0.02; Santiago's cohort). However, this association was not confirmed in a second replication migraine series. A review of the literature reveals the existence of inconsistent findings and methodological shortcomings affecting a large proportion of mtDNA association studies on AD, PD, and migraine. A detailed inspection of the literature highlights the need for performing more rigorous methodological and statistical standards in mtDNA genetic association studies aimed to avoid false positive results of association between mtDNA variants and neurological diseases.

Cuba: exploring the history of admixture and the genetic basis of pigmentation using autosomal and uniparental markers.

We carried out an admixture analysis of a sample comprising 1,019 individuals from all the provinces of Cuba. We used a panel of 128 autosomal Ancestry Informative Markers (AIMs) to estimate the admixture proportions. We also characterized a number of haplogroup diagnostic markers in the mtDNA and Y-chromosome in order to evaluate admixture using uniparental markers. Finally, we analyzed the association of 16 single nucleotide polymorphisms (SNPs) with quantitative estimates of skin pigmentation. In the total sample, the average European, African and Native American contributions as estimated from autosomal AIMs were 72%, 20% and 8%, respectively. The Eastern provinces of Cuba showed relatively higher African and Native American contributions than the Western provinces. In particular, the highest proportion of African ancestry was observed in the provinces of Guantánamo (40%) and Santiago de Cuba (39%), and the highest proportion of Native American ancestry in Granma (15%), Holguín (12%) and Las Tunas (12%). We found evidence of substantial population stratification in the current Cuban population, emphasizing the need to control for the effects of population stratification in association studies including individuals from Cuba. The results of the analyses of uniparental markers were concordant with those observed in the autosomes. These geographic patterns in admixture proportions are fully consistent with historical and archaeological information. Additionally, we identified a sex-biased pattern in the process of gene flow, with a substantially higher European contribution from the paternal side, and higher Native American and African contributions from the maternal side. This sex-biased contribution was particularly evident for Native American ancestry. Finally, we observed that SNPs located in the genes SLC24A5 and SLC45A2 are strongly associated with melanin levels in the sample.

Evaluating the role of mitochondrial DNA variation to the genetic predisposition to radiation-induced toxicity.

BACKGROUND AND PURPOSE: Mitochondrial DNA common variants have been reported to be associated with the development of radiation-induced toxicity. Using a large cohort of patients, we aimed to validate these findings by investigating the potential role of common European mitochondrial DNA SNPs (mtSNPs) to the development of radio-toxicity. MATERIAL AND METHODS: Overall acute and late toxicity data were assessed in a cohort of 606 prostate cancer patients by means of Standardized Total Average Toxicity (STAT) score. We carried out association tests between radiation toxicity and a selection of 15 mtSNPs (and the haplogroups defined by them). RESULTS: Statistically significant association between mtSNPs and haplogroups with toxicity could not be validated in our Spanish cohort. CONCLUSIONS: The present study suggests that the mtDNA common variants analyzed are not associated with clinically relevant increases in risk of overall radiation-induced toxicity in prostate cancer patients.

No evidence that major mtDNA European haplogroups confer risk to schizophrenia.

Previous studies suggest that genetic factors could be involved in mitochondrial dysfunction observed in schizophrenia (SZ), some of them claiming a role of mtDNA common variants (mtSNPs) and/or haplogroups (hgs) in developing this disorder. These studies, however, have mainly been undertaken on relatively small cohorts of patients and control individuals and most have not yet been replicated. To further analyze the role of mtSNPs in SZ risk, we have carried out the largest genotyping effort to date using two Spanish case-control samples comprising a total of 942 schizophrenic patients and 1,231 unrelated controls: 454 patients and 616 controls from Santiago de Compostela (Galicia) and 488 patients and 615 controls from Reus (Catalonia). A set of 25 mtSNPs representing main branches of the European mtDNA phylogeny were genotyped in the Galician cohort and a subset of 16 out of these 25 mtSNPs was genotyped in the Catalan cohort. These 16 common variants characterize the most common European branches of the mtDNA phylogeny. We did not observe any positive association of mtSNPs and hgs with SZ. We discuss several deficiencies of previous studies that might explain the false positive nature of previous findings, including the confounding effect of population sub-structure and deficient statistical methodologies. It is unlikely that mtSNPs defining the most common European mtDNA haplogroups are related to SZ.

The impact of modern migrations on present-day multi-ethnic Argentina as recorded on the mitochondrial DNA genome.

BACKGROUND: The genetic background of Argentineans is a mosaic of different continental ancestries. From colonial to present times, the genetic contribution of Europeans and sub-Saharan Africans has superposed to or replaced the indigenous genetic 'stratum'. A sample of 384 individuals representing different Argentinean provinces was collected and genotyped for the first and the second mitochondrial DNA (mtDNA) hypervariable regions, and selectively genotyped for mtDNA SNPs. This data was analyzed together with additional 440 profiles from rural and urban populations plus 304 from Native American Argentineans, all available from the literature. A worldwide database was used for phylogeographic inferences, inter-population comparisons, and admixture analysis. Samples identified as belonging to hg (hg) H2a5 were sequenced for the entire mtDNA genome. RESULTS: Phylogenetic and admixture analyses indicate that only half of the Native American component in urban Argentineans might be attributed to the legacy of extinct ancestral Argentineans and that the Spanish genetic contribution is slightly higher than the Italian one. Entire H2a5 genomes linked these Argentinean mtDNAs to the Basque Country and improved the phylogeny of this Basque autochthonous clade. The fingerprint of African slaves in urban Argentinean mtDNAs was low and it can be phylogeographically attributed predominantly to western African. The European component is significantly more prevalent in the Buenos Aires province, the main gate of entrance for Atlantic immigration to Argentina, while the Native American component is larger in North and South Argentina. AMOVA, Principal Component Analysis and hgs/haplotype patterns in Argentina revealed an important level of genetic sub-structure in the country. CONCLUSIONS: Studies aimed to compare mtDNA frequency profiles from different Argentinean geographical regions (e.g., forensic and case-control studies) should take into account the important genetic heterogeneity of the country in order to prevent false positive claims of association in disease studies or inadequate evaluation of forensic evidence.

The mitochondrial genome is a "genetic sanctuary" during the oncogenic process.

Since Otto Warburg linked mitochondrial physiology and oncogenesis in the 1930s, a number of studies have focused on the analysis of the genetic basis for the presence of aerobic glycolysis in cancer cells. However, little or no evidence exists today to indicate that mtDNA mutations are directly responsible for the initiation of tumor onset. Based on a model of gliomagenesis in the mouse, we aimed to explore whether or not mtDNA mutations are associated with the initiation of tumor formation, maintenance and aggressiveness. We reproduced the different molecular events that lead from tumor initiation to progression in the mouse glioma. In human gliomas, most of the genetic alterations that have been previously identified result in the aberrant activation of different signaling pathways and deregulation of the cell cycle. Our data indicates that mitochondrial dysfunction is associated with reactive oxygen species (ROS) generation, leading to increased nuclear DNA (nDNA) mutagenesis, but maintaining the integrity of the mitochondrial genome. In addition, mutational stability has been observed in entire mtDNA of human gliomas; this is in full agreement with the results obtained in the cancer mouse model. We use this model as a paradigm of oncogenic transformation due to the fact that mutations commonly found in gliomas appear to be the most common molecular alterations leading to tumor development in most types of human cancer. Our results indicate that the mtDNA genome is kept by the cell as a "genetic sanctuary" during tumor development in the mouse and humans. This is compatible with the hypothesis that the mtDNA molecule plays an essential role in the control of the cellular adaptive survival response to tumor-induced oxidative stress. The integrity of mtDNA seems to be a necessary element for responding to the increased ROS production associated with the oncogenic process.

Reassessing the role of mitochondrial DNA mutations in autism spectrum disorder.

BACKGROUND: There is increasing evidence that impairment of mitochondrial energy metabolism plays an important role in the pathophysiology of autism spectrum disorders (ASD; OMIM number: 209850). A significant proportion of ASD cases display biochemical alterations suggestive of mitochondrial dysfunction and several studies have reported that mutations in the mitochondrial DNA (mtDNA) molecule could be involved in the disease phenotype. METHODS: We analysed a cohort of 148 patients with idiopathic ASD for a number of mutations proposed in the literature as pathogenic in ASD. We also carried out a case control association study for the most common European haplogroups (hgs) and their diagnostic single nucleotide polymorphisms (SNPs) by comparing cases with 753 healthy and ethnically matched controls. RESULTS: We did not find statistical support for an association between mtDNA mutations or polymorphisms and ASD. CONCLUSIONS: Our results are compatible with the idea that mtDNA mutations are not a relevant cause of ASD and the frequent observation of concomitant mitochondrial dysfunction and ASD could be due to nuclear factors influencing mitochondrion functions or to a more complex interplay between the nucleus and the mitochondrion/mtDNA.

New population and phylogenetic features of the internal variation within mitochondrial DNA macro-haplogroup R0.

BACKGROUND: R0 embraces the most common mitochondrial DNA (mtDNA) lineage in West Eurasia, namely, haplogroup H (approximately 40%). R0 sub-lineages are badly defined in the control region and therefore, the analysis of diagnostic coding region polymorphisms is needed in order to gain resolution in population and medical studies. METHODOLOGY/PRINCIPAL FINDINGS: We sequenced the first hypervariable segment (HVS-I) of 518 individuals from different North Iberian regions. The mtDNAs belonging to R0 (approximately 57%) were further genotyped for a set of 71 coding region SNPs characterizing major and minor branches of R0. We found that the North Iberian Peninsula shows moderate levels of population stratification; for instance, haplogroup V reaches the highest frequency in Cantabria (north-central Iberia), but lower in Galicia (northwest Iberia) and Catalonia (northeast Iberia). When compared to other European and Middle East populations, haplogroups H1, H3 and H5a show frequency peaks in the Franco-Cantabrian region, declining from West towards the East and South Europe. In addition, we have characterized, by way of complete genome sequencing, a new autochthonous clade of haplogroup H in the Basque country, named H2a5. Its coalescence age, 15.6+/-8 thousand years ago (kya), dates to the period immediately after the Last Glacial Maximum (LGM). CONCLUSIONS/SIGNIFICANCE: In contrast to other H lineages that experienced re-expansion outside the Franco-Cantabrian refuge after the LGM (e.g. H1 and H3), H2a5 most likely remained confined to this area till present days.