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.

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.

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.

HIrisPlex-S system for eye, hair, and skin color prediction from DNA: Massively parallel sequencing solutions for two common forensically used platforms.

Forensic DNA Phenotyping (FDP) provides the ability to predict externally visible characteristics from minute amounts of crime scene DNA, which can help find unknown perpetrators who are typically unidentifiable via conventional forensic DNA profiling. Fundamental human genetics research has led to a better understanding of the specific DNA variants responsible for physical appearance characteristics, particularly eye, hair, and skin color. Recently, we introduced the HIrisPlex-S system for the simultaneous prediction of eye, hair, and skin color based on 41 DNA variants generated from two forensically validated SNaPshot multiplex assays using capillary electrophoresis (CE). Here we introduce massively parallel sequencing (MPS) solutions for the HIrisPlex-S (HPS) system on two MPS platforms commonly used in forensics, Ion Torrent and MiSeq, that cover all 41 DNA variants in a single assay, respectively. Additionally, we present the forensic developmental validation of the two HPS-MPS assays. The Ion Torrent MPS assay, based on Ion AmpliSeq technology, illustrated the successful generation of full HIrisPlex-S genotypic profiles from 100 pg of input control DNA, while the MiSeq MPS assay based on an in-house design yielded complete profiles from 250 pg of input DNA. Assessing simulated forensic casework samples such as saliva, hair (bulb), blood, semen, and low quantity touch DNA, as well as artificially damaged DNA samples, concordance testing, and samples from numerous species, all illustrated the ability of both versions of the HIrisPlex-S MPS assay to produce results that motivate forensic applications. By also providing an integrated bioinformatics analysis pipeline, MPS data can now be analyzed and a file generated for upload to the publically accessible HIrisPlex online webtool (https://hirisplex.erasmusmc.nl). In addition, we updated the website to accept VCF input data for those with genome sequence data. We thus provide a user-friendly and semi-automated MPS workflow from DNA sample to individual eye, hair, and skin color prediction probabilities. Furthermore, we present a 2-person mixture separation tool that not only assesses genotype reliability with regards genotyping confidence but also provides the most fitting mixture scenario for both minor and major contributors, including profile separation. We envision this MPS implementation of the HIrisPlex-S system for eye, hair, and skin color prediction from DNA as a starting point for further expanding MPS-based forensic DNA phenotyping. This may include the future addition of SNPs predictive for more externally visible characteristics, as well as SNPs for bio-geographic ancestry inference, provided the statistical framework for DNA prediction of these traits is in place.

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.

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.

Dog breed affiliation with a forensically validated canine STR set.

We tested a panel of 13 highly polymorphic canine short tandem repeat (STR) markers for dog breed assignment using 392 dog samples from the 23 most popular breeds in Austria, Germany, and Switzerland. This STR panel had originally been selected for canine identification. The dog breeds sampled in this study featured a population frequency ≥1% and accounted for nearly 57% of the entire pedigree dog population in these three countries. Breed selection was based on a survey comprising records for nearly 1.9 million purebred dogs belonging to more than 500 different breeds. To derive breed membership from STR genotypes, a range of algorithms were used. These methods included discriminant analysis of principal components (DAPC), STRUCTURE, GeneClass2, and the adegenet package for R. STRUCTURE analyses suggested 21 distinct genetic clusters. Differentiation between most breeds was clearly discernable. Fourteen of 23 breeds (61%) exhibited maximum mean cluster membership proportions of more than 0.70 with a highest value of 0.90 found for Cavalier King Charles Spaniels. Dogs of only 6 breeds (26%) failed to consistently show only one major cluster. The DAPC method yielded the best assignment results in all 23 declared breeds with 97.5% assignment success. The frequency-based assignment test also provided a high success rate of 87%. These results indicate the potential viability of dog breed prediction using a well-established and sensitive set of 13 canine STR markers intended for forensic routine use.

Inferring biogeographic ancestry with compound markers of slow and fast evolving polymorphisms.

Bio-geographic ancestry is an area of considerable interest in the medical genetics, anthropology and forensics. Although genome-wide panels are ideal as they provide dense genotyping data, small sets of ancestry informative marker provide a cost-effective way to investigate genetic ancestry and population structure. Here, we investigate the performance of a reduced marker set that combine different types of autosomal markers through haplotype analysis. In particular, recently described DIP-STR markers should offer the advantage of comprising both, low mutation rate Indels (DIPs), to study human history over longer time scale; and high mutation rate STRs, to trace relatively recent demographic events. In this study, we assessed the ability of an initial set of 23 DIP-STRs to distinguish major population groups using the HGDP-CEPH reference samples. The results obtained applying the STRUCTURE algorithm show that the discrimination capacity of the DIP-STRs is comparable to currently used small-scale ancestry informative markers by approaching seven major demographic groups. Yet, the DIP-STRs show an improved success rate in assigning individuals to populations of Europe and Middle East. These data show a remarkable ability of a preliminary set of 23 DIP-STR markers to infer major biogeographic origins. A novel set of DIP-STRs preselected to contain ancestry information should lead to further improvements.

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.

Evaluation of the predictive capacity of DNA variants associated with straight hair in Europeans.

DNA-based prediction of hair morphology, defined as straight, curly or wavy hair, could contribute to an improved description of an unknown offender and allow more accurate forensic reconstructions of physical appearance in the field of forensic DNA phenotyping. Differences in scalp hair morphology are significant at the worldwide scale and within Europe. The only genome-wide association study made to date revealed the Trichohyalin gene (TCHH) to be significantly associated with hair morphology in Europeans and reported weaker associations for WNT10A and FRAS1 genes. We conducted a study that centered on six SNPs located in these three genes with a sample of 528 individuals from Poland. The predictive capacity of the candidate DNA variants was evaluated using logistic regression; classification and regression trees; and neural networks, by applying a 10-fold cross validation procedure. Additionally, an independent test set of 142 males from six European populations was used to verify performance of the developed prediction models. Our study confirmed association of rs11803731 (TCHH), rs7349332 (WNT10A) and rs1268789 (FRAS1) SNPs with hair morphology. The combined genotype risk score for straight hair had an odds ratio of 2.7 and these predictors explained ∼ 8.2% of the total variance. The selected three SNPs were found to predict straight hair with a high sensitivity but low specificity when a 10-fold cross validation procedure was applied and the best results were obtained using the neural networks approach (AUC=0.688, sensitivity=91.2%, specificity=23.0%). Application of the neural networks model with 65% probability threshold on an additional test set gave high sensitivity (81.4%) and improved specificity (50.0%) with a total of 78.7% correct calls, but a high non-classification rate (66.9%). The combined TTGGGG SNP genotype for rs11803731, rs7349332, rs1268789 (European frequency=4.5%) of all six straight hair-associated alleles was identified as the best predictor, giving >80% probability of straight hair. Finally, association testing of 44 SNPs previously identified to be associated with male pattern baldness revealed a suggestive association with hair morphology for rs4679955 on 3q25.1. The study results reported provide the starting point for the development of a predictive test for hair morphology in Europeans. More studies are now needed to discover additional determinants of hair morphology to improve the predictive accuracy of this trait in forensic analysis.

Comparison of the genetic background of different Colombian populations using the SNPforID 52plex identification panel.

Various strategies for analysing SNP markers and genotyping have been published with the goal of obtaining informative profiles from biological samples that contain only small amounts of template and/or degraded DNA. In this study, a multiplex assay of 52 autosomal single-nucleotide polymorphisms (SNPs) was used to analyse 438 individuals from urban populations from different regions of Colombia, as well as a sample of 50 Native American individuals of the Pastos ethnic group from Nariño. To determine if significant differences in these 52 SNPs exist between the distinct regions of Colombia, genetic distance and admixture analyses were performed based on the available data for 17 different Colombian population groups and for population groups from Africa, Europe and America. The results demonstrate significant differences between the populations from the Southwest Andean, Central-West Andean, Central-East Andean, Orinoquian and northern Colombian Pacific Coast regions. Most of the regions exhibited a European and Native American admixture. One exception is the population from the region of Chocó (on the northern Pacific Coast), which exhibits a high proportion of African admixture (54 %). From the observed genetic backgrounds, it is possible to conclude that a single reference database for the entire country would not be suitable for forensic purposes. The allele frequencies and the forensically relevant parameters were calculated for all of the markers in each Colombian region with significant values for the combined matching probability (power of discrimination ≥0.99999999999999990) and the combined probability of exclusion (≥0.9990) in trios that were obtained from all of the population groups.