Forensic identity SNPs: Characterisation of flanking region variation using massively parallel sequencing.

Single nucleotide polymorphisms (SNPs) can be analysed for identity or kinship applications in forensic genetics to either provide an adjunct to traditional STR typing or as a stand-alone approach. The advent of massively parallel sequencing technology (MPS) has provided a useful opportunity to more easily deploy SNP typing in a forensic context, given the ability to simultaneously amplify a large number of markers. Furthermore, MPS also provides valuable sequence data for the targeted regions, which enables the detection of any additional variation seen in the flanking regions of amplicons. In this study we genotyped 977 samples across five UK-relevant population groups (White British, East Asian, South Asian, North-East African and West African) for 94 identity-informative SNP markers using the ForenSeq DNA Signature Prep Kit. Examination of flanking region variation allowed for the identification of 158 additional alleles across all populations studied. Here we present allele frequencies for all 94 identity-informative SNPs, both including and excluding the flanking region sequence of these markers. We also present information on the configuration of these SNPs in the ForenSeq DNA Signature Prep Kit, including performance metrics for the markers and investigation of bioinformatic and chemistry-based discordances. Overall, the inclusion of flanking region variation in the analysing workflow for these markers reduced the average combined match probability 2175 times across all populations, with a maximum reduction of 675,000-fold in the West African population. The gain due to flanking region-based discrimination increased the heterozygosity of some loci above that of some of the least useful forensic STR loci; thus demonstrating the benefit of enhanced analysis of currently targeted SNP markers for forensic applications.

The Y chromosome and its use in forensic DNA analysis.

Originally relatively ignored in forensic investigations because its genetic analysis lacks inference of individual identification, the value of Y chromosome analysis has been proven in cases of sexual assault, particularly where the amount of material left by a male assailant is limited in comparison with female DNA. All routine analysis of autosomal DNA, however, targets a gene (AMELY) on the Y chromosome in order to identify the sex of the DNA source and this is discussed in the context of the genetic structure of this male-specific chromosome. Short-tandem repeat markers on the chromosome are tested in dedicated multiplexes that have developed over time and these are described alongside international guidance as to their use in a forensic setting. As a marker of lineage, the Y chromosome provides additional tools to assist in the inference of ancestry, both geographical and familial and the value of Y chromosome testing is illustrated through descriptions of cases of criminal and historical interest. A decision to analyse the Y chromosome has to be considered in the context, not only of the circumstances of the case, but also with regard to the ethical questions it might raise, and these are discussed in relation to the cases that have been described in more detail in the accompanying online supplementary material.

Mitochondrial DNA in forensic use.

Genetic analysis of mitochondrial DNA (mtDNA) has always been a useful tool for forensic geneticists, mainly because of its ubiquitous presence in biological material, even in the absence of nuclear DNA. Sequencing, however, is not a skill that is part of the routine forensic analysis because of the relative rarity of requests, and the need for retention of necessary skill sets and associated accreditation issues. While standard Sanger sequencing may be relatively simple, many requests are made in the face of compromised biological samples. Newer technologies, provided through massively parallel sequencing (MPS), will increase the opportunity for scientists to include this tool in their routine, particularly for missing person investigations. MPS has also enabled a different approach to sequencing that can increase sensitivity in a more targeted approach. In these circumstances it is likely that only a laboratory that specialises in undertaking forensic mtDNA analysis will be able to take these difficult cases forward, more so because reviews of the literature have revealed significantly high levels of typing errors in publications reporting mtDNA sequences. The forensic community has set out important guidelines, not only in the practical aspects of analysis, but also in the interpretation of that sequence to ensure that accurate comparisons can be made. Analysis of low-level, compromised and ancient DNA is not easy, however, as contamination is extremely difficult to eliminate and circumstances leading to sequencing errors are all too easily introduced. These problems, and solutions, are discussed in the article in relation to several historic cases.

Reply to Response to Vacuous standards - Subversion of the OSAC standards-development process.

This Letter to the Editor is a reply to Mohammed et al. (2021) https://doi.org/10.1016/j.fsisyn.2021.100145, which in turn is a response to Morrison et al. (2020) "Vacuous standards - subversion of the OSAC standards-development process" https://doi.org/10.1016/j.fsisyn.2020.06.005.

Classification of STR allelic variation using massively parallel sequencing and assessment of flanking region power.

The application of massively parallel sequencing (MPS) to forensic genetics has led to improvements in multiple aspects of DNA analysis, however, additional complexities are concurrently associated with these advances. In relation to short tandem repeat (STR) typing, the move to sequence rather than length-based methodologies has highlighted the extent to which previous allelic variation was masked - both within and outside of the repeat regions (the flanking regions). In order to fully implement MPS for autosomal STR analysis, sequence-based allelic frequencies must be available, and concordance with previous typing techniques needs to be assessed. In this work, a series of samples (n = 1007) from five different population groups were genotyped using the MiSeq FGx™ Forensic Genomics System. Results were compared to those obtained using capillary electrophoresis (CE), and sequence variation has been characterised both within and outside STR repeat regions, with allelic frequencies provided for all variants observed within this database. Analysing and characterising flanking region sequence is currently less straightforward than studying repeat region variation alone, and the added value of doing so remains largely unexplored - this paper provides data to show that the gain in polymorphism achieved when analysing flanking regions is less than might be expected. In the White British population for example, including the sequence variation within repeat regions of 26 autosomal STRs made the average combined random match probability (RMP) over 700 times lower than with length-based alleles alone. Including the sequence variation within the flanking regions only resulted in a combined RMP that was a further 4 times lower.

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.

DNA methylation-based age prediction using massively parallel sequencing data and multiple machine learning models.

The field of DNA intelligence focuses on retrieving information from DNA evidence that can help narrow down large groups of suspects or define target groups of interest. With recent breakthroughs on the estimation of geographical ancestry and physical appearance, the estimation of chronological age comes to complete this circle of information. Recent studies have identified methylation sites in the human genome that correlate strongly with age and can be used for the development of age-estimation algorithms. In this study, 110 whole blood samples from individuals aged 11-93 years were analysed using a DNA methylation quantification assay based on bisulphite conversion and massively parallel sequencing (Illumina MiSeq) of 12 CpG sites. Using this data, 17 different statistical modelling approaches were compared based on root mean square error (RMSE) and a Support Vector Machine with polynomial function (SVMp) model was selected for further testing. For the selected model (RMSE = 4.9 years) the mean average error (MAE) of the blind test (n = 33) was calculated at 4.1 years, with 52% of the samples predicting with less than 4 years of error and 86% with less than 7 years. Furthermore, the sensitivity of the method was assessed both in terms of methylation quantification accuracy and prediction accuracy in the first validation of this kind. The described method retained its accuracy down to 10 ng of initial DNA input or ∼2 ng bisulphite PCR input. Finally, 34 saliva samples were analysed and following basic normalisation, the chronological age of the donors was predicted with less than 4 years of error for 50% of the samples and with less than 7 years of error for 70%.

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.

Forensic genealogy: Some serious concerns.

Genealogical databases have provided links to possible perpetrators of crimes in several cold cases in the US. This commentary discusses some of the ethical issues associated with this approach while recognising the underlying value of the identifications.

UK and Irish Y-STR population data-A catalogue of variant alleles.

A total of 3128 Y-STR profiles from three UK and one Irish population have been analysed with the PowerPlex Y23 system and are reported here. Instances of haplotype sharing between apparently unrelated individuals were identified and further investigated with the use of the 5 additional markers within the Yfiler Plus kit, resulting in a reduction by 76% in the number of shared haplotypes. Furthermore, Yfiler Plus was also employed to verify locus deletions and duplications observed in Y23 genotypes while inconsistencies between the two kits were sequenced, revealing underlying Y23 primer binding site mutations in loci DYS392 and DYS576. Finally, the mechanism behind a previously reported population specific peak shift observed in DYS481 in South Asian samples has been evaluated and further investigated in a novel case of this phenomenon seen in a Black British individual featuring a different flanking region mutation.

Concordance of the ForenSeq™ system and characterisation of sequence-specific autosomal STR alleles across two major population groups.

By using sequencing technology to genotype loci of forensic interest it is possible to simultaneously target autosomal, X and Y STRs as well as identity, ancestry and phenotypic informative SNPs, resulting in a breadth of data obtained from a single run that is considerable when compared to that generated with standard technologies. It is important however that this information aligns with the genotype data currently obtained using commercially available kits for CE-based investigations such that results are compatible with existing databases and hence can be of use to the forensic community. In this work, 400 samples were typed using commercially available STR kits and CE, as well as using the Ilumina ForenSeq™ DNA Signature Prep Kit and MiSeq® FGx to assess concordance of autosomal STRs and population variability. Results show a concordance rate between the two technologies exceeding 99.98% while numerous novel sequence based alleles are described. In order to make use of the sequence variation observed, sequence specific allele frequencies were generated for White British and British Chinese populations.

DNA methylation-based forensic age prediction using artificial neural networks and next generation sequencing.

The ability to estimate the age of the donor from recovered biological material at a crime scene can be of substantial value in forensic investigations. Aging can be complex and is associated with various molecular modifications in cells that accumulate over a person's lifetime including epigenetic patterns. The aim of this study was to use age-specific DNA methylation patterns to generate an accurate model for the prediction of chronological age using data from whole blood. In total, 45 age-associated CpG sites were selected based on their reported age coefficients in a previous extensive study and investigated using publicly available methylation data obtained from 1156 whole blood samples (aged 2-90 years) analysed with Illumina's genome-wide methylation platforms (27K/450K). Applying stepwise regression for variable selection, 23 of these CpG sites were identified that could significantly contribute to age prediction modelling and multiple regression analysis carried out with these markers provided an accurate prediction of age (R2=0.92, mean absolute error (MAE)=4.6 years). However, applying machine learning, and more specifically a generalised regression neural network model, the age prediction significantly improved (R2=0.96) with a MAE=3.3 years for the training set and 4.4 years for a blind test set of 231 cases. The machine learning approach used 16 CpG sites, located in 16 different genomic regions, with the top 3 predictors of age belonged to the genes NHLRC1, SCGN and CSNK1D. The proposed model was further tested using independent cohorts of 53 monozygotic twins (MAE=7.1 years) and a cohort of 1011 disease state individuals (MAE=7.2 years). Furthermore, we highlighted the age markers' potential applicability in samples other than blood by predicting age with similar accuracy in 265 saliva samples (R2=0.96) with a MAE=3.2 years (training set) and 4.0 years (blind test). In an attempt to create a sensitive and accurate age prediction test, a next generation sequencing (NGS)-based method able to quantify the methylation status of the selected 16 CpG sites was developed using the Illumina MiSeq® platform. The method was validated using DNA standards of known methylation levels and the age prediction accuracy has been initially assessed in a set of 46 whole blood samples. Although the resulted prediction accuracy using the NGS data was lower compared to the original model (MAE=7.5years), it is expected that future optimization of our strategy to account for technical variation as well as increasing the sample size will improve both the prediction accuracy and reproducibility.

Discovery of potential DNA methylation markers for forensic tissue identification using bisulphite pyrosequencing.

The presence of specific body fluids at crime scenes could be linked with particular types of crime, therefore attributing a DNA profile to a specific tissue could increase the evidential significance of a match with a suspect. Current methodologies such as tissue-specific mRNA profiling are useful but drawbacks include low tissue specificity and applicability to degraded samples. In this study, the potential of 11 tissue-specific differentially methylated regions, initially identified following large-scale methylation analysis of whole blood, buccal cells and sperm, was explored in order to identify markers for blood, saliva and semen. Bisulphite pyrosequencing analysis supported previous findings, but tissue-specific differentially methylated regions for blood and buccal cells did not show enough specificity to be proposed as markers for blood and saliva, respectively. For some CpGs, a large inter-individual variation in methylation levels was also observed. Two of the semen markers (cg04382920 and cg11768416) were used for further validation on a large set of stains. These two semen-specific assays showed high sensitivity (as low as 50 pg) and stability. Future experiments will shed light on the usefulness of these markers in forensic casework.

Comprehensive Analysis of Pan-African Mitochondrial DNA Variation Provides New Insights into Continental Variation and Demography.

Africa is the cradle of all human beings, and although it has been the focus of a number of genetic studies, there are many questions that remain unresolved. We have performed one of the largest and most comprehensive meta-analyses of mitochondrial DNA (mtDNA) lineages carried out in the African continent to date. We generated high-throughput mtDNA single nucleotide polymorphism (SNP) data (230 SNPs) from 2024 Africans, where more than 500 of them were additionally genotyped for the control region. These data were analyzed together with over 12,700 control region profiles collected from the literature, representing more than 300 population samples from Africa. Insights into the African homeland of humans are discussed. Phylogeographic patterns for the African continent are shown at a high phylogeographic resolution as well as at the population and regional levels. The deepest branch of the mtDNA tree, haplogroup L0, shows the highest sub-haplogroup diversity in Southeast and East Africa, suggesting this region as the homeland for modern humans. Several demographic estimates point to the coast as a facilitator of human migration in Africa, but the data indicate complex patterns, perhaps mirroring the effect of recent continental-scaled demographic events in re-shaping African mtDNA variability.

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.

Evaluation of DNA variants associated with androgenetic alopecia and their potential to predict male pattern baldness.

Androgenetic alopecia, known in men as male pattern baldness (MPB), is a very conspicuous condition that is particularly frequent among European men and thus contributes markedly to variation in physical appearance traits amongst Europeans. Recent studies have revealed multiple genes and polymorphisms to be associated with susceptibility to MPB. In this study, 50 candidate SNPs for androgenetic alopecia were analyzed in order to verify their potential to predict MPB. Significant associations were confirmed for 29 SNPs from chromosomes X, 1, 5, 7, 18 and 20. A simple 5-SNP prediction model and an extended 20-SNP model were developed based on a discovery panel of 305 males from various European populations fitting one of two distinct phenotype categories. The first category consisted of men below 50 years of age with significant baldness and the second; men aged 50 years or older lacking baldness. The simple model comprised the five best predictors: rs5919324 near AR, rs1998076 in the 20p11 region, rs929626 in EBF1, rs12565727 in TARDBP and rs756853 in HDAC9. The extended prediction model added 15 SNPs from five genomic regions that improved overall prevalence-adjusted predictive accuracy measured by area under the receiver characteristic operating curve (AUC). Both models were evaluated for predictive accuracy using a test set of 300 males reflecting the general European population. Applying a 65% probability threshold, high prediction sensitivity of 87.1% but low specificity of 42.4% was obtained in men aged <50 years. In men aged ≥50, prediction sensitivity was slightly lower at 67.7% while specificity reached 90%. Overall, the AUC=0.761 calculated for men at or above 50 years of age indicates these SNPs offer considerable potential for the application of genetic tests to predict MPB patterns, adding a highly informative predictive system to the emerging field of forensic analysis of externally visible characteristics.

A novel application of real-time RT-LAMP for body fluid identification: using HBB detection as the model.

We report on a novel application of real-time reverse transcription-loop-mediated isothermal amplification (real-time RT-LAMP) to identify the presence of a specific body fluid using blood as a proof-of-concept model. By comparison with recently developed methods of body fluid identification, the RT-LAMP assay is rapid and requires only one simple heating-block maintained at a single temperature, circumventing the need for dedicated equipment. RNA was extracted from different body fluids (blood, semen, saliva, menstrual blood, sweat, and urine) for use in real-time RT-LAMP reaction. The 18S rRNA locus was used as the internal control and hemoglobin beta (HBB) as the blood-specific marker. Reverse transcription and LAMP reaction were performed in the same tube using a turbidimeter for real-time monitoring the reaction products within a threshold of 60 min. HBB LAMP products were only detected in blood and not in any of the other body fluid, but products from the 18S rRNA gene were detected in all the tested body fluids as expected. The limit of detection was a minimum of 10(-5) ng total RNA for detection of both 18S rRNA and HBB. Augmenting the detection of RT-LAMP products was performed by separation of the products using gel electrophoresis and collecting the fluorescence of calcein. The data collected indicated complete concordance with the body fluid tested regardless of the method of detection used. This is the first application of real-time RT-LAMP to detect body fluid specific RNA and indicates the use of this method in forensic biology.