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.

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.

A global metagenomic map of urban microbiomes and antimicrobial resistance

We present a global atlas of 4,728 metagenomic samples from mass-transit systems in 60 cities over 3 years, representing the first systematic, worldwide catalog of the urban microbial ecosystem. This atlas provides an annotated, geospatial profile of microbial strains, functional characteristics, antimicrobial resistance (AMR) markers, and genetic elements, including 10,928 viruses, 1,302 bacteria, 2 archaea, and 838,532 CRISPR arrays not found in reference databases. We identified 4,246 known species of urban microorganisms and a consistent set of 31 species found in 97% of samples that were distinct from human commensal organisms. Profiles of AMR genes varied widely in type and density across cities. Cities showed distinct microbial taxonomic signatures that were driven by climate and geographic differences. These results constitute a high-resolution global metagenomic atlas that enables discovery of organisms and genes, highlights potential public health and forensic applications, and provides a culture-independent view of AMR burden in cities.

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.

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

•

Blood DNA methylation profiles of 1156 individuals were assessed for age correlation.

•

Stepwise regression identified 23 age-associated CpG sites in DNA from blood.

•

A machine learning model based on 16 markers predicted age with a mean error of 3.8 years.

•

The model predicted age successfully for twins and ‘diseased’ individuals.

•

A new NGS-based method was combined with machine learning for age prediction.

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 (27 K/450 K). 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 (R² = 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 (R² = 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 (R² = 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.5 years), 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.

Brief report: isogenic induced pluripotent stem cell lines from an adult with mosaic down syndrome model accelerated neuronal ageing and neurodegeneration

Trisomy 21 (T21), Down Syndrome (DS) is the most common genetic cause of dementia and intellectual disability. Modeling DS is beginning to yield pharmaceutical therapeutic interventions for amelioration of intellectual disability, which are currently being tested in clinical trials. DS is also a unique genetic system for investigation of pathological and protective mechanisms for accelerated ageing, neurodegeneration, dementia, cancer, and other important common diseases. New drugs could be identified and disease mechanisms better understood by establishment of well-controlled cell model systems. We have developed a first nonintegration-reprogrammed isogenic human induced pluripotent stem cell (iPSC) model of DS by reprogramming the skin fibroblasts from an adult individual with constitutional mosaicism for DS and separately cloning multiple isogenic T21 and euploid (D21) iPSC lines. Our model shows a very low number of reprogramming rearrangements as assessed by a high-resolution whole genome CGH-array hybridization, and it reproduces several cellular pathologies seen in primary human DS cells, as assessed by automated high-content microscopic analysis. Early differentiation shows an imbalance of the lineage-specific stem/progenitor cell compartments: T21 causes slower proliferation of neural and faster expansion of hematopoietic lineage. T21 iPSC-derived neurons show increased production of amyloid peptide-containing material, a decrease in mitochondrial membrane potential, and an increased number and abnormal appearance of mitochondria. Finally, T21-derived neurons show significantly higher number of DNA double-strand breaks than isogenic D21 controls. Our fully isogenic system therefore opens possibilities for modeling mechanisms of developmental, accelerated ageing, and neurodegenerative pathologies caused by T21.

Worldwide F(ST) estimates relative to five continental-scale populations

We estimate the population genetics parameter (also referred to as the fixation index) from short tandem repeat (STR) allele frequencies, comparing many worldwide human subpopulations at approximately the national level with continental-scale populations. is commonly used to measure population differentiation, and is important in forensic DNA analysis to account for remote shared ancestry between a suspect and an alternative source of the DNA. We estimate comparing subpopulations with a hypothetical ancestral population, which is the approach most widely used in population genetics, and also compare a subpopulation with a sampled reference population, which is more appropriate for forensic applications. Both estimation methods are likelihood-based, in which is related to the variance of the multinomial-Dirichlet distribution for allele counts. Overall, we find low values, with posterior 97.5 percentiles when comparing a subpopulation with the most appropriate population, and even for inter-population comparisons we find . These are much smaller than single nucleotide polymorphism-based inter-continental estimates, and are also about half the magnitude of STR-based estimates from population genetics surveys that focus on distinct ethnic groups rather than a general population. Our findings support the use of up to 3% in forensic calculations, which corresponds to some current practice.

Species detection using HyBeacon(®) probe technology: Working towards rapid onsite testing in non-human forensic and food authentication applications

Identifying individual species or determining species' composition in an unknown sample is important for a variety of forensic applications. Food authentication, monitoring illegal trade in endangered species, forensic entomology, sexual assault case work and counter terrorism are just some of the fields that can require the detection of the biological species present. Traditional laboratory based approaches employ a wide variety of tools and technologies and exploit a number of different species specific traits including morphology, molecular differences and immuno-chemical analyses. A large number of these approaches require laboratory based apparatus and results can take a number of days to be returned to investigating authorities. Having a presumptive test for rapid identification could lead to savings in terms of cost and time and allow sample prioritisation if confirmatory testing in a laboratory is required later. This model study describes the development of an assay using a single HyBeacon(®) probe and melt curve analyses allowing rapid screening and authentication of food products labelled as Atlantic cod (Gadus morhua). Exploiting melt curve detection of species specific SNP sites on the COI gene the test allows detection of a target species (Atlantic cod) and closely related species which may be used as substitutes. The assay has been designed for use with the Field Portable ParaDNA system, a molecular detection platform for non-expert users. The entire process from sampling to result takes approximately 75min. Validation studies were performed on both single source genomic DNA, mixed genomic DNA and commercial samples. Data suggests the assay has a lower limit of detection of 31 pg DNA. The specificity of the assay to Atlantic cod was measured by testing highly processed food samples including frozen, defrosted and cooked fish fillets as well as fish fingers, battered fish fillet and fish pie. Ninety-six (92.7%) of all Atlantic cod food products, tested, provided a correct single species result with the remaining samples erroneously identified as containing non-target species. The data shows that the assay was quick to design and characterise and is also capable of yielding results that would be beneficial in a variety of fields, not least the authentication of food.

Persistence of DNA from laundered semen stains: Implications for child sex trafficking cases

In sexual assault cases, particularly those involving internal child sex trafficking (ICST), victims often hide their semen-stained clothing. This can result in a lag time of several months before the items are laundered and subsequently seized during a criminal investigation. Although it has been demonstrated previously that DNA can be recovered from clothing washed immediately after semen deposition, laundered items of clothing are not routinely examined in ICST cases, due to the assumption that the time delay and washing would result in no detectable DNA. The aim of this study was to examine whether viable DNA profiles could be recovered from laundered semen stains where there has been a significant lag time between semen deposition from one or more individuals and one or more washes of the stained clothing. Items of UK school uniform (T-shirts, trousers, tights) were stained with fresh semen (either from a single donor or a 1:1 mixture from two donors) and stored in a wardrobe for eight months. Stained and unstained items (socks) were then washed at 30 °C or 60 °C and with non-biological or biological detergent. DNA samples extracted from the semen-stained sites and from the unstained socks were quantified and profiled. High quantities of DNA, (6-18 μg) matching the DNA profiles of the semen donors, were recovered from all semen-stained clothing that had been laundered once, irrespective of wash conditions. This quantity,and profile quality,did not decline significantly with multiple washes. The two donor semen samples yielded ∼ 10-fold more DNA from the T-shirts than from the trousers. This disparity resulted in the T-shirts yielding a ∼ 1:1 mixture of DNA from the two donors, whereas the trousers yielded a major DNA profile matching only that of the second donor. The quantities of DNA recovered from the unstained socks were an order of magnitude lower, with most of the DNA being attributable to the donor of the semen on the stained clothing within the same wash, demonstrating the transfer of semen-derived DNA among items of clothing in the washing machine. This study demonstrates that complete DNA profiles can be obtained from laundered semen stains on school uniform-type clothing, with an eight-month lag time between semen deposition and laundering, despite multiple washes and stains from two semen donors. These data emphasise the need to recover and examine the clothing of victims for semen and DNA evidence, even if the clothing has been stored for several months or washed multiple times since the sexual offence took place.

A global analysis of Y-chromosomal haplotype diversity for 23 STR loci

In a worldwide collaborative effort, 19,630 Y-chromosomes were sampled from 129 different populations in 51 countries. These chromosomes were typed for 23 short-tandem repeat (STR) loci (DYS19, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS385ab, DYS437, DYS438, DYS439, DYS448, DYS456, DYS458, DYS635, GATAH4, DYS481, DYS533, DYS549, DYS570, DYS576, and DYS643) and using the PowerPlex Y23 System (PPY23, Promega Corporation, Madison, WI). Locus-specific allelic spectra of these markers were determined and a consistently high level of allelic diversity was observed. A considerable number of null, duplicate and off-ladder alleles were revealed. Standard single-locus and haplotype-based parameters were calculated and compared between subsets of Y-STR markers established for forensic casework. The PPY23 marker set provides substantially stronger discriminatory power than other available kits but at the same time reveals the same general patterns of population structure as other marker sets. A strong correlation was observed between the number of Y-STRs included in a marker set and some of the forensic parameters under study. Interestingly a weak but consistent trend toward smaller genetic distances resulting from larger numbers of markers became apparent.

A report of the 1995 and 1996 Paternity Testing Workshops of the English Speaking Working Group of the International Society for Forensic Haemogenetics

We report the results of the 1995 and 1996 Paternity Testing Workshops of the English Speaking Working Group of the International Society for Forensic Haemogenetics. In 1995, 18 laboratories participated and in 1996, 21 laboratories participated. Each year, blood samples from three persons (child, mother and alleged father) were sent to participating laboratories which performed paternity testing according to their usual protocols. The results and answers to questionnaires concerning methods were compiled and are presented in this report. From the questionnaires, a general tendency to a more frequent use of polymerase chain reaction (PCR) based methods was seen. In 1996, 62% of the laboratories used PCR based methods. Ten per cent of the laboratories used only PCR based methods. The remaining 90% of the laboratories performed restriction fragment length polymorphism (RFLP) investigations of variable numbers of tandem repeat (VNTR) loci with single locus probes (SLPs) either alone or in combination with PCR based typing, multi locus probing, classical systems (ABO etc.), or serological HLA typing. In 1996, typing with classical systems was used in 29% of the laboratories. The majority of the laboratories performed RFLP typing of VNTR loci using very similar methods. The results and the inter-laboratory variations of the measured lengths of the DNA-fragments of the VNTR regions D2S44, D7S21, D7S22, and D12S11 of the trios were analysed. The overall coefficient of variation was 2.15% in 1995 and 1.43% in 1996. During the period 1991-1996, the inter-laboratory variation has decreased, most probably due to the fact that the methods have now been optimised and the majority of the participating laboratories have adopted the standardised method for RFLP typing with SLPs which was agreed upon for investigations in crime cases by the European DNA Profiling Group. In 1996, eight laboratories reported the results of PCR based typing of the short tandem repeat (STR) locus HumTH01, six laboratories reported results of HumVWA31A typing, and five laboratories reported the results of typing of the STR locus HumF13A1 and the VNTR locus D1S80. The results were concordant although the nomenclature was slightly inconsistent concerning the classification of an irregular repeat of the HumTH01 system.

An investigation of the HUMVWA31A locus in British Caucasians

A number of short tandem repeat (STR) loci are currently being examined for their usefulness as markers of identity; HUMVWA31A is one such locus. We used a high-sieving agarose technique to type 200 British Caucasians for this locus. Comparison of the resultant allele frequencies with other published databases showed that their distributions were similar. Observed heterozygosity was similar to that reported in other population studies but significantly lower than expected. A goodness of fit analysis of observed and expected genotypes was highly significant, suggesting a deviation from Hardy-Weinberg equilibrium, although reports from other populations appeared not to show such differences. The possibility that mistyping had led to an erroneous statistic was investigated by detailed examination of samples within our own laboratory and in two other laboratories, involving verification of the agarose typing by both automated fluorescent detection techniques, and sequencing. Other reported deviations from Hardy-Weinberg equilibrium in STR systems are discussed.