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.

Combining current knowledge on DNA methylation-based age estimation towards the development of a superior forensic DNA intelligence tool.

The estimation of chronological age from biological fluids has been an important quest for forensic scientists worldwide, with recent approaches exploiting the variability of DNA methylation patterns with age in order to develop the next generation of forensic 'DNA intelligence' tools for this application. Drawing from the conclusions of previous work utilising massively parallel sequencing (MPS) for this analysis, this work introduces a DNA methylation-based age estimation method for blood that exhibits the best combination of prediction accuracy and sensitivity reported to date. Statistical evaluation of markers from 51 studies using microarray data from over 4000 individuals, followed by validation using in-house generated MPS data, revealed a final set of 11 markers with the greatest potential for accurate age estimation from minimal DNA material. Utilising an algorithm based on support vector machines, the proposed model achieved an average error (MAE) of 3.3 years, with this level of accuracy retained down to 5 ng of starting DNA input (~ 1 ng PCR input). The accuracy of the model was retained (MAE = 3.8 years) in a separate test set of 88 samples of Spanish origin, while predictions for donors of greater forensic interest (< 55 years of age) displayed even higher accuracy (MAE = 2.6 years). Finally, no sex-related bias was observed for this model, while there were also no signs of variation observed between control and disease-associated populations for schizophrenia, rheumatoid arthritis, frontal temporal dementia and progressive supranuclear palsy in microarray data relating to the 11 markers.

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.

Broadening the Applicability of a Custom Multi-Platform Panel of Microhaplotypes: Bio-Geographical Ancestry Inference and Expanded Reference Data.

The development of microhaplotype (MH) panels for massively parallel sequencing (MPS) platforms is gaining increasing relevance for forensic analysis. Here, we expand the applicability of a 102 autosomal and 11 X-chromosome panel of MHs, previously validated with both MiSeq and Ion S5 MPS platforms and designed for identification purposes. We have broadened reference population data for identification purposes, including data from 240 HGDP-CEPH individuals of native populations from North Africa, the Middle East, Oceania and America. Using the enhanced population data, the panel was evaluated as a marker set for bio-geographical ancestry (BGA) inference, providing a clear differentiation of the five main continental groups of Africa, Europe, East Asia, Native America, and Oceania. An informative degree of differentiation was also achieved for the population variation encompassing North Africa, Middle East, Europe, South Asia, and East Asia. In addition, we explored the potential for individual BGA inference from simple mixed DNA, by simulation of mixed profiles followed by deconvolution of mixture components.

Performance of ancestry-informative SNP and microhaplotype markers.

The use of microhaplotypes (MHs) for ancestry inference has added to an increasing number of ancestry-informative markers (AIMs) for forensic application that includes autosomal single nucleotide polymorphisms (SNPs) and insertions/deletions (indels). This study compares bi-allelic and tri-allelic SNPs as well as MH markers for their ability to differentiate African, European, South Asian, East Asian, and American population groups from the 1000 Genomes Phase 3 database. A range of well-established metrics were applied to rank each marker according to the population differentiation potential they measured. These comprised: absolute allele frequency differences (δ); Rosenberg's informativeness for (ancestry) assignment (In); the fixation index (FST); and the effective number of alleles (Ae). A panel consisting of all three marker types resulted in the lowest mean divergence per population per individual (MDPI = 2.16%) when selected by In. However, when marker types were not mixed, MHs were the highest performing markers by most metrics (MDPI < 4%) for differentiation between the five continental populations.

Global patterns of STR sequence variation: Sequencing the CEPH human genome diversity panel for 58 forensic STRs using the Illumina ForenSeq DNA Signature Prep Kit.

The 944 individuals of the CEPH human genome diversity panel (HGDP-CEPH), a standard sample set of 51 globally distributed populations, were sequenced using the Illumina ForenSeq™ DNA Signature Prep Kit. The ForenSeq™ system is a single multiplex for the MiSeq/FGx™ massively parallel sequencing instrument, comprising: amelogenin, 27 autosomal STRs, 24 Y-STRs, 7 X-STRs, and 94 SNPforID+Kiddlab autosomal ID-SNPs (plus optionally detected ancestry and phenotyping SNP sets). We report in detail the patterns of sequence variation observed in the repeat regions of the 58 forensic STR loci typed by the ForenSeq™ system. Sequence alleles were characterized and repeat region structures annotated by aligning the ForenSeq™ sequence output to the latest GRCh38 human reference sequence, necessitating the reversal and re-alignment of STR allele sequences reported by the Forenseq™ system in 20 of 58 STRs (plus the reverse alleles in two Y-STRs with duplicated-inverted repeat regions). Individual population sample sizes of the HGDP-CEPH panel do not allow reliable inferences to be made about levels of genetic variability in low frequency STR alleles-where particular sequence variants are found in only a few individuals; but we assessed the occurrence of both population-specific sequence variants and singleton observations; finding each of these in a sizeable proportion of HGDP-CEPH samples, with consequences for planning the co-ordinated compilation of sequence variation on a much larger scale than was required before by forensic laboratories now adopting massively parallel sequencing.

Ancestry analysis in rural Brazilian populations of African descent.

Rural communities comprise around 20% of Caribbean and South American populations, but are under-represented in autosomal marker databases. That deficiency is problematic for forensic genetics, as it relies on accurate descriptions of genetic variation and population structure. Brazilian populations were shaped by an intense, complex and heterogeneous process of admixture encompassing mainly Amerindians, Sub-Saharan Africans and Europeans. Quilombos are Brazilian populations with significant African descent that have remained genetically isolated to some extent from surrounding populations. In the reported study, we analyzed three rural Quilombo populations: Kalunga; Riacho de Sacutiaba e Sacutiaba; and Mocambo, along with a dataset from the HGDP-CEPH panel. Aiming to contribute to representative genetic databases of forensic interest, we analyse the three rural Quilombos populations and investigate how their genetic makeup relates to their history by analyzing an established forensic test, comprising 46 ancestry-informative (AIM) Indels. The panel was chosen for its high power in differentiating the main contributing populations of Brazil. Parental populations were selected from HGDP-CEPH data available at the forInDel allele frequency browser based on historic patterns applicable to the study populations and the amount of variability observed within and between continents. Our results show the main admixture components in the Quilombos are African and European. Those estimates are in accordance with previous analyses for both uniparental and autosomal markers. PCA, structure analysis and ancestry estimates indicate a correlation between the extent of isolation and the degree of admixture in the Quilombos: Kalunga is the most isolated population and accordingly has a higher African admixture component (67.3%). Sacutiaba is the smallest and most impacted by migration, with the highest European component (46.8%). Mocambo neighbors a Native American population and therefore has the highest Amerindian contribution (12.2%). Our results are consistent with the history and demography of Quilombos. The heterogeneity observed in these populations stresses the genetic diversity that Latin American and Caribbean rural populations can have and reiterates the need to describe them in greater detail.

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.

Modified aging of elite athletes revealed by analysis of epigenetic age markers.

Recent progress in epigenomics has led to the development of prediction systems that enable accurate age estimation from DNA methylation data. Our objective was to track responses to intense physical exercise of individual age-correlated DNA methylation markers and to infer their potential impact on the aging processes. The study showed accelerated DNA hypermethylation for two CpG sites in TRIM59 and KLF14. Both markers predicted the investigated elite athletes to be several years older than controls and this effect was more substantial in subjects involved in power sports. Accordingly, the complete 5-CpG model revealed age acceleration of elite athletes (P=1.503x10-7) and the result was more significant amongst power athletes (P=1.051x10-9). The modified methylation of TRIM59 and KLF14 in top athletes may be accounted for by the biological roles played by these genes. Their known anti-tumour and anti-inflammatory activities suggests that intense physical training has a complex influence on aging and potentially launches signalling networks that contribute to the observed lower risk of elite athletes to develop cardiovascular disease and cancer.

Inference of Ancestry in Forensic Analysis I: Autosomal Ancestry-Informative Marker Sets.

An expanding choice of ancestry-informative marker single nucleotide polymorphisms (AIM-SNPs) is becoming available for the forensic user in the form of sensitive SNaPshot-based tests or in alternative single-base extension genotyping systems (e.g., Sequenom iPLEX) that can be adapted for analysis with SNaPshot. In addition, alternative ancestry-informative variation: Indels and STRs can be analyzed using direct PCR-to-CE techniques that offer the possibility to detect mixed profiles. We review the current forensically viable AIM panels, their optimized PCR multiplexes, and the population differentiation power they offer. We also describe how improved population divergence balance can be achieved with the enlarged multiplex scales of next-generation sequencing approaches to enable analysis of admixed individuals without biased estimation of co-ancestry proportions.

Inference of Ancestry in Forensic Analysis II: Analysis of Genetic Data.

Three approaches applicable to the analysis of forensic ancestry-informative marker data-STRUCTURE, principal component analysis, and the Snipper Bayesian classification system-are reviewed. Detailed step-by-step guidance is provided for adjusting parameter settings in STRUCTURE with particular regard to their effect when differentiating populations. Several enhancements to the Snipper online forensic classification portal are described, highlighting the added functionality they bring to particular aspects of ancestry-informative SNP analysis in a forensic context.

Completion of a worldwide reference panel of samples for an ancestry informative Indel assay.

The use of ancestry informative markers (AIMs) in forensic analysis is of considerable utility since ancestry inference can progress an investigation when no identification has been made of DNA from the crime-scene. Short-amplicon markers, including insertion deletion polymorphisms, are particularly useful in forensic analysis due to their mutational stability, capacity to amplify degraded samples and straightforward amplification technique. In this study we report the completion of H952 HGDP-CEPH panel genotyping with a set of 46 AIM-Indels. The study adds Central South Asian and Middle Eastern population data, allowing a comparison of patterns of variation in Eurasia for these markers, in order to enhance their use in forensic analyses, particularly when combined with sets of ancestry informative SNPs. Ancestry analysis using principal component analysis and Bayesian methods indicates that a proportion of classification error occurs with European-Middle East population comparisons, but the 46 AIM-Indels have the capability to differentiate six major population groups when European-Central South Asian comparisons are made. These findings have relevance for forensic ancestry analyses in countries where South Asians form much of the demographic profile, including the UK, USA and South Africa. A novel third allele detected in MID-548 was characterized - despite a low frequency in the HGDP-CEPH panel samples, it appears confined to Central South Asian populations, increasing the ability to differentiate this population group. The H952 data set was implemented in a new open access SPSmart frequency browser - forInDel: Forensic Indel browser.

SNP Markers as Additional Information to Resolve Complex Kinship Cases.

BACKGROUND: DNA profiling with sets of highly polymorphic autosomal short tandem repeat (STR) markers has been applied in various aspects of human identification in forensic casework for nearly 20 years. However, in some cases of complex kinship investigation, the information provided by the conventionally used STR markers is not enough, often resulting in low likelihood ratio (LR) calculations. In these cases, it becomes necessary to increment the number of loci under analysis to reach adequate LRs. Recently, it has been proposed that single nucleotide polymorphisms (SNPs) could be used as a supportive tool to STR typing, eventually even replacing the methods/markers now employed. METHODS: In this work, we describe the results obtained in 7 revised complex paternity cases when applying a battery of STRs, as well as 52 human identification SNPs (SNPforID 52plex identification panel) using a SNaPshot methodology followed by capillary electrophoresis. RESULTS: Our results show that the analysis of SNPs, as complement to STR typing in forensic casework applications, would at least increase by a factor of 4 total PI values and correspondent Essen-Möller's W value. CONCLUSIONS: We demonstrated that SNP genotyping could be a key complement to STR information in challenging casework of disputed paternity, such as close relative individualization or complex pedigrees subject to endogamous relations.

"New turns from old STaRs": enhancing the capabilities of forensic short tandem repeat analysis.

The field of research and development of forensic STR genotyping remains active, innovative, and focused on continuous improvements. A series of recent developments including the introduction of a sixth dye have brought expanded STR multiplex sizes while maintaining sensitivity to typical forensic DNA. New supplementary kits complimenting the core STRs have also helped improve analysis of challenging identification cases such as distant pairwise relationships in deficient pedigrees. This article gives an overview of several recent key developments in forensic STR analysis: availability of expanded core STR kits and supplementary STRs, short-amplicon mini-STRs offering practical options for highly degraded DNA, Y-STR enhancements made from the identification of rapidly mutating loci, and enhanced analysis of genetic ancestry by analyzing 32-STR profiles with a Bayesian forensic classifier originally developed for SNP population data. As well as providing scope for genotyping larger numbers of STRs optimized for forensic applications, the launch of compact next-generation sequencing systems provides considerable potential for genotyping the sizeable proportion of nucleotide variation existing in forensic STRs, which currently escapes detection with CE.

Development of a novel forensic STR multiplex for ancestry analysis and extended identity testing.

There is growing interest in developing additional DNA typing techniques to provide better investigative leads in forensic analysis. These include inference of genetic ancestry and prediction of common physical characteristics of DNA donors. To date, forensic ancestry analysis has centered on population-divergent SNPs but these binary loci cannot reliably detect DNA mixtures, common in forensic samples. Furthermore, STR genotypes, forming the principal DNA profiling system, are not routinely combined with forensic SNPs to strengthen frequency data available for ancestry inference. We report development of a 12-STR multiplex composed of ancestry informative marker STRs (AIM-STRs) selected from 434 tetranucleotide repeat loci. We adapted our online Bayesian classifier for AIM-SNPs: Snipper, to handle multiallele STR data using frequency-based training sets. We assessed the ability of the 12-plex AIM-STRs to differentiate CEPH Human Genome Diversity Panel populations, plus their informativeness combined with established forensic STRs and AIM-SNPs. We found combining STRs and SNPs improves the success rate of ancestry assignments while providing a reliable mixture detection system lacking from SNP analysis alone. As the 12 STRs generally show a broad range of alleles in all populations, they provide highly informative supplementary STRs for extended relationship testing and identification of missing persons with incomplete reference pedigrees. Lastly, mixed marker approaches (combining STRs with binary loci) for simple ancestry inference tests beyond forensic analysis bring advantages and we discuss the genotyping options available.

SNPs as Supplements in Simple Kinship Analysis or as Core Markers in Distant Pairwise Relationship Tests: When Do SNPs Add Value or Replace Well-Established and Powerful STR Tests?

BACKGROUND: Genetic tests for kinship testing routinely reach likelihoods that provide virtual proof of the claimed relationship by typing microsatellites-commonly consisting of 12-15 standard forensic short tandem repeats (STRs). Single nucleotide polymorphisms (SNPs) have also been applied to kinship testing but these binary markers are required in greater numbers than multiple-allele STRs. However SNPs offer certain advantageous characteristics not found in STRs, including, much higher mutational stability, good performance typing highly degraded DNA, and the ability to be readily up-scaled to very high marker numbers reaching over a million loci. This article outlines kinship testing applications where SNPs markedly improve the genetic data obtained. In particular we explore the minimum number of SNPs that will be required to confirm pairwise relationship claims in deficient pedigrees that typify missing persons' identification or war grave investigations where commonly few surviving relatives are available for comparison and the DNA is highly degraded. METHODS: We describe the application of SNPs alongside STRs when incomplete profiles or allelic instability in STRs create ambiguous results, we review the use of high density SNP arrays when the relationship claim is very distant, and we outline simulations of kinship analyses with STRs supplemented with SNPs in order to estimate the practical limit of pairwise relationships that can be differentiated from random unrelated pairs from the same population. RESULTS: The minimum number of SNPs for robust statistical inference of parent-offspring relationships through to those of second cousins (S-3-3) is estimated for both simple, single multiplex SNP sets and for subsets of million-SNP arrays. CONCLUSIONS: There is considerable scope for resolving ambiguous STR results and for improving the statistical power of kinship analysis by adding small-scale SNP sets but where the pedigree is deficient the pairwise relationships must be relatively close. For more distant relationships it is possible to reduce chip-based SNP arrays from the million+ markers down to ∼7,000. However, such numbers indicate that current genotyping approaches will not be able to deliver sufficient data to resolve distant pairwise relationships from the limited DNA typical of the most challenging identification cases.

Differentiation of African components of ancestry to stratify groups in a case-control study of a Brazilian urban population.

BACKGROUND: Balancing the subject composition of case and control groups to create homogenous ancestries between each group is essential for medical association studies. METHODS: We explored the applicability of single-tube 34-plex ancestry informative markers (AIM) single nucleotide polymorphisms (SNPs) to estimate the African Component of Ancestry (ACA) to design a future case-control association study of a Brazilian urban sample. RESULTS: One hundred eighty individuals (107 case group; 73 control group) self-described as white, brown-intermediate or black were selected. The proportions of the relative contribution of a variable number of ancestral population components were similar between case and control groups. Moreover, the case and control groups demonstrated similar distributions for ACA <0.25 and >0.50 categories. Notably a high number of outlier values (23 samples) were observed among individuals with ACA <0.25. These individuals presented a high probability of Native American and East Asian ancestral components; however, no individuals originally giving these self-described ancestries were observed in this study. CONCLUSIONS: The strategy proposed for the assessment of ancestry and adjustment of case and control groups for an association study is an important step for the proper construction of the study, particularly when subjects are taken from a complex urban population. This can be achieved using a straight forward multiplexed AIM-SNPs assay of highly discriminatory ancestry markers.

Differentiation of african components of ancestry to stratify groups in a Case–control study of a brazilian urban population

Background: Balancing the subject composition of case and control groups to create homogenous ancestries between each group is essential for medical association studies. Methods: We explored the applicability of single-tube 34-plex ancestry informative markers (AIM) single nucleotide polymorphisms (SNPs) to estimate the African Component of Ancestry (ACA) to design a future case-control association study of a Brazilian urban sample. Results: One hundred eighty individuals (107 case group; 73 control group) self-described as white, brown-intermediate or black were selected. The proportions of the relative contribution of a variable number of ancestral population components were similar between case and control groups. Moreover, the case and control groups demonstrated similar distributions for ACA 0.50 categories. Notably a high number of outlier values (23 samples) were observed among individuals with ACA <0.25. These individuals presented a high probability of Native American and East Asian ancestral components; however, no individuals originally giving these self-described ancestries were observed in this study. Conclusions: The strategy proposed for the assessment of ancestry and adjustment of case and control groups for an association study is an important step for the proper construction of the study, particularly when subjects are taken from a complex urban population. This can be achieved using a straight forward multiplexed AIM-SNPs assay of highly discriminatory ancestry markers.

Mitochondrial echoes of first settlement and genetic continuity in El Salvador.

BACKGROUND: From Paleo-Indian times to recent historical episodes, the Mesoamerican isthmus played an important role in the distribution and patterns of variability all around the double American continent. However, the amount of genetic information currently available on Central American continental populations is very scarce. In order to shed light on the role of Mesoamerica in the peopling of the New World, the present study focuses on the analysis of the mtDNA variation in a population sample from El Salvador. METHODOLOGY/PRINCIPAL FINDINGS: We have carried out DNA sequencing of the entire control region of the mitochondrial DNA (mtDNA) genome in 90 individuals from El Salvador. We have also compiled more than 3,985 control region profiles from the public domain and the literature in order to carry out inter-population comparisons. The results reveal a predominant Native American component in this region: by far, the most prevalent mtDNA haplogroup in this country (at approximately 90%) is A2, in contrast with other North, Meso- and South American populations. Haplogroup A2 shows a star-like phylogeny and is very diverse with a substantial proportion of mtDNAs (45%; sequence range 16090-16365) still unobserved in other American populations. Two different Bayesian approaches used to estimate admixture proportions in El Salvador shows that the majority of the mtDNAs observed come from North America. A preliminary founder analysis indicates that the settlement of El Salvador occurred about 13,400+/-5,200 Y.B.P.. The founder age of A2 in El Salvador is close to the overall age of A2 in America, which suggests that the colonization of this region occurred within a few thousand years of the initial expansion into the Americas. CONCLUSIONS/SIGNIFICANCE: As a whole, the results are compatible with the hypothesis that today's A2 variability in El Salvador represents to a large extent the indigenous component of the region. Concordant with this hypothesis is also the observation of a very limited contribution from European and African women ( approximately 5%). This implies that the Atlantic slave trade had a very small demographic impact in El Salvador in contrast to its transformation of the gene pool in neighbouring populations from the Caribbean facade.