Towards simultaneous individual and tissue identification: A proof-of-principle study on parallel sequencing of STRs, amelogenin, and mRNAs with the Ion Torrent PGM

Forensic Microbiology: When, Where and How

Forensic microbiology is a relatively new discipline, born in part thanks to the development of advanced methodologies for the detection, identification and characterization of microorganisms, and also in relation to the growing impact of infectious diseases of iatrogenic origin. Indeed, the increased application of medical practices, such as transplants, which require immunosuppressive treatments, and the growing demand for prosthetic installations, associated with an increasing threat of antimicrobial resistance, have led to a rise in the number of infections of iatrogenic origin, which entails important medico-legal issues. On the other hand, the possibility of detecting minimal amounts of microorganisms, even in the form of residual traces (e.g., their nucleic acids), and of obtaining gene and genomic sequences at contained costs, has made it possible to ask new questions of whether cases of death or illness might have a microbiological origin, with the possibility of also tracing the origin of the microorganisms involved and reconstructing the chain of contagion. In addition to the more obvious applications, such as those mentioned above related to the origin of iatrogenic infections, or to possible cases of infections not properly diagnosed and treated, a less obvious application of forensic microbiology concerns its use in cases of violence or violent death, where the characterization of the microorganisms can contribute to the reconstruction of the case. Finally, paleomicrobiology, e.g., the reconstruction and characterization of microorganisms in historical or even archaeological remnants, can be considered as a sister discipline of forensic microbiology. In this article, we will review these different aspects and applications of forensic microbiology.

Body Fluid Identification in Samples Collected after Intimate and Social Contact: A Comparison of Two mRNA Profiling Methods and the Additional Information Gained by cSNP Genotypes

The ability to associate a contributor with a specific body fluid in a crime stain can aid casework investigation. The detection of body fluids combined with DNA analyses may supply essential information, but as the two tests are independent, they may not be associated. Recently, the analysis of coding region SNPs (cSNPs) within the RNA transcript has been proven to be a promising method to face this challenge. In this study, we performed targeted RNA sequencing of 158 samples (boxershorts, fingernail swabs and penile swabs) collected from 12 couples at different time points post-intimate contact and after non-intimate contact, using the Ion S5™ System and BFID-cSNP-6F assay. The aim of the study was to compare the performance of the MPS and CE methods in the detection of mRNA markers, and to associate body fluids with contributors by their cSNP genotypes. The results of the study show a lower success rate in the detection of vaginal mucosa by the MPS compared to the CE method. However, the additional information obtained with the cSNP genotypes could successfully associate body fluids with contributors in most cases.

The forensic landscape and the population genetic analyses of Hainan Li based on massively parallel sequencing DNA profiling

Due to the formation of the Qiongzhou Strait by climate change and marine transition, Hainan island was isolated from the mainland southern China during the Last Glacial Maximum. Hainan island, located at the southernmost part of China and separated from the Leizhou Peninsula by the Qiongzhou Strait, laid on one of the modern human northward migration routes from Southeast Asia to East Asia. The Hlai language-speaking Li minority, the second largest population after Han Chinese in Hainan island, is the direct descendants of the initial migrants in Hainan island and has unique ethnic properties and derived characteristics; however, the forensic-associated studies on Hainan Li population are still insufficient. Hence, 136 Hainan Li individuals were genotyped in this study using the MPS-based ForenSeq™ DNA Signature Prep Kit (DNA Primer Set A, DPMA) to characterize the forensic genetic polymorphism landscape, and DNA profiles were obtained from 152 different molecular genetic markers (27 autosomal STRs, 24 Y-STRs, 7 X-STRs, and 94 iiSNPs). A total of 419 distinct length variants and 586 repeat sequence sub-variants, with 31 novel alleles (at 17 loci), were identified across the 58 STR loci from the DNA profiles of Hainan Li population. We evaluated the forensic characteristics and efficiencies of DPMA, demonstrating that the STRs and iiSNPs in DPMA were highly polymorphic in Hainan Li population and could be employed in forensic applications. In addition, we set up three datasets, which included the genetic data of (i) iiSNPs (27 populations, 2640 individuals), (ii) Y-STRs (42 populations, 8281 individuals), and (iii) Y haplogroups (123 populations, 4837 individuals) along with the population ancestries and language families, to perform population genetic analyses separately from different perspectives. In conclusion, the phylogenetic analyses indicated that Hainan Li, with a southern East Asia origin and Tai-Kadai language-speaking language, is an isolated population relatively. But the genetic pool of Hainan Li influenced by the limited gene flows from other Tai-Kadai populations and Hainan populations. Furthermore, the establishment of isolated population models will be beneficial to clarify the exquisite population structures and develop specific genetic markers for subpopulations in forensic genetic fields.

Ten years of molecular ballistics-a review and a field guide

Molecular ballistics combines molecular biological, forensic ballistic, and wound ballistic insights and approaches in the description, collection, objective investigation, and contextualization of the complex patterns of biological evidence that are generated by gunshots at biological targets. Setting out in 2010 with two seminal publications proving the principle that DNA from backspatter collected from inside surfaces of firearms can be retreived and successfully be analyzed, molecular ballistics covered a lot of ground until today. In this review, 10 years later, we begin with a comprehensive description and brief history of the field and lay out its intersections with other forensic disciplines like wound ballistics, forensic molecular biology, blood pattern analysis, and crime scene investigation. In an application guide section, we aim to raise consciousness to backspatter traces and the inside surfaces of firearms as sources of forensic evidence. Covering crime scene practical as well as forensic genetic aspects, we introduce operational requirements and lay out possible procedures, including forensic RNA analysis, when searching for, collecting, analyzing, and contextualizing such trace material. We discuss the intricacies and rationales of ballistic model building, employing different tissue, skin, and bone simulants and the advantages of the “triple-contrast” method in molecular ballistics and give advice on how to stage experimental shootings in molecular ballistic research. Finally, we take a look at future applications and prospects of molecular ballistics.

Forensic transcriptome analysis using massively parallel sequencing

The application of transcriptome analyses in forensic genetics has experienced tremendous growth and development in the past decade. The earliest studies and main applications were body fluid and tissue identification, using targeted RNA transcripts and a reverse transcription endpoint PCR method. A number of markers have been identified for the forensically most relevant body fluids and tissues and the method has been successfully used in casework. The introduction of Massively Parallel Sequencing (MPS) opened up new perspectives and opportunities to advance the field. Contrary to genomic DNA where two copies of an autosomal DNA segment are present in a cell, abundant RNA species are expressed in high copy numbers. Even whole transcriptome sequencing (RNA-Seq) of forensically relevant body fluids and of postmortem material was shown to be possible. This review gives an overview on forensic transcriptome analyses and applications. The methods cover whole transcriptome as well as targeted MPS approaches. High resolution forensic transcriptome analyses using MPS are being applied to body fluid/ tissue identification, determination of the age of stains and the age of the donor, the estimation of the post-mortem interval and to post mortem death investigations.

Massive parallel sequencing in forensics: advantages, issues, technicalities, and prospects

In the last decade, next-generation sequencing (NGS) technology, alternatively massive parallel sequencing (MPS), was applied to all fields of biological research. Its introduction to the field of forensics was slower, mainly due to lack of accredited sequencers, kits, and relatively higher sequencing error rates as compared with standardized Sanger sequencing. Currently, a majority of the problematic issues have been solved, which is proven by the body of reports in the literature. Here, we discuss the utility of NGS sequencing in forensics, emphasizing the advantages, issues, the technical aspects of the experiments, commercial solutions, and the potentially interesting applications of MPS.

Highly degraded RNA can still provide molecular information: An in vitro approach

The long-term survival of RNA in postmortem tissues is a tricky topic. Many aged/forensic specimens show, in fact, high rates of null/inconclusive PCR-based results, while reliable outcomes were sometimes achieved from archaeological samples. On the other hand, several data show that the RNA is a molecule that survives even to several physical-chemical stresses. In the present study, a simple protocol, which was already developed for the prolonged hydrolysis of DNA, was applied to a RNA sample extracted from blood. This protocol is based on the heat-mediated (70°C) hydrolysis for up to 36 h using ultrapure water and di-ethyl-pyro-carbonate-water as hydrolysis medium. Measurable levels of depurination were not found even if microfluidic devices showed a progressive pattern of degradation. The reverse transcription/quantitative PCR analysis of two (60 bp long) housekeeping targets (glyceraldehyde-3-phosphate dehydrogenase and porphobilinogen deaminase) showed that the percentage of amplifiable target (%AT) decreased in relation to the duration of the damaging treatment (r² > 0.973). The comparison of the %AT in the degraded RNA and in the DNA samples that underwent the same damaging treatment showed that the %AT is always higher in RNA, reaching up to three orders of magnitude. Lastly, even the end-point PCR of blood-specific markers gave reliable results, which is in agreement with the body fluid origin of the sample. In conclusion, all the PCR-based results show that RNA maintains the ability to be retro-transcribed in short cDNA fragments even after 36 h of incubation at 70°C in mildly acidic buffers. It is therefore likely that the long-term survival of RNA samples depends mainly on the protection against RNAase attacks rather than on environmental factors (such as humidity and acidity) that are instead of great importance for the stability of DNA. As a final remark, our results suggest that the RNA analysis can be successfully performed even when DNA profiling failed.

Novel taxonomy-independent deep learning microbiome approach allows for accurate classification of different forensically relevant human epithelial materials

Correct identification of different human epithelial materials such as from skin, saliva and vaginal origin is relevant in forensic casework as it provides crucial information for crime reconstruction. However, the overlap in human cell type composition between these three epithelial materials provides challenges for their differentiation and identification when using previously proposed human cell biomarkers, while their microbiota composition largely differs. By using validated 16S rRNA gene massively parallel sequencing data from the Human Microbiome Project of 1636 skin, oral and vaginal samples, 50 taxonomy-independent deep learning networks were trained to classify these three tissues. Validation testing was performed in de-novo generated high-throughput 16S rRNA gene sequencing data using the Ion Torrent^™ Personal Genome Machine from 110 test samples: 56 hand skin, 31 saliva and 23 vaginal secretion specimens. Body-site classification accuracy of these test samples was very high as indicated by AUC values of 0.99 for skin, 0.99 for oral, and 1 for vaginal secretion. Misclassifications were limited to 3 (5%) skin samples. Additional forensic validation testing was performed in mock casework samples by de-novo high-throughput sequencing of 19 freshly-prepared samples and 22 samples aged for 1 up to 7.6 years. All of the 19 fresh and 20 (91%) of the 22 aged mock casework samples were correctly tissue-type classified. Moreover, comparing the microbiome results with outcomes from previous human mRNA-based tissue identification testing in the same 16 aged mock casework samples reveals that our microbiome approach performs better in 12 (75%), similarly in 2 (12.5%), and less good in 2 (12.5%) of the samples. Our results demonstrate that this new microbiome approach allows for accurate tissue-type classification of three human epithelial materials of skin, oral and vaginal origin, which is highly relevant for future forensic investigations.

"The acid test"-validation of the ParaDNA® Body Fluid ID Test for routine forensic casework

The identification of the cellular origin and composition of crime scene-related traces can provide crucial insight into a crime scene reconstruction. In the last decade, especially mRNA-based body fluid and tissue identification (BFI) has been vigorously examined. Besides capillary electrophoretic (CE) and real-time quantitative PCR (RT-qPCR)-based approaches for mRNA detection, melt curve analysis bears potential as a simple-to-use method for BFI. The ParaDNA® Body Fluid ID Test relies on HyBeacon® probes and was developed as a rapid test for mRNA-based BFI of six different body fluids: vaginal fluid, seminal fluid, sperm cells, saliva, menstrual, and peripheral blood. The herein presented work was performed as an "acid test" of the system and should clarify whether the approach matches the requirements of forensic routine casework in German police departments. Tested samples consisted of single source as well as of mixed samples.

Assessment of impact of DNA extraction methods on analysis of human remain samples on massively parallel sequencing success

Skeletal remains recovered from missing persons' cases are often exposed to harsh environmental conditions resulting in the DNA being damaged, degraded, and/or the samples containing PCR inhibitors. In this study, the efficacy of common extraction methods was evaluated to remove high levels of PCR inhibitors commonly encountered with human remains, and their downstream compatibility with the two leading sequencing chemistries and platforms for human identification purposes. Blood, hair, and bone samples were spiked with high levels of inhibitors commonly identified in each particular substrate in order to test the efficiency of various DNA extraction methods prior to sequencing. Samples were extracted using three commercial extraction kits (DNA IQ™, DNA Investigator, and PrepFiler® BTA), organic (blood and hair only), and two total demineralization protocols (bone only)). Massively parallel sequencing (MPS) was performed using two different systems: Precision ID chemistry and a custom AmpliSeq™ STR and iiSNP panel on the Ion S5™ System and the ForenSeq DNA Signature Prep Kit on the MiSeq FGx™. The overall results showed that all DNA extraction methods were efficient and are fully compatible with both MPS systems. Key performance indicators such as STR and SNP reportable alleles, read depth, and heterozygote balance were comparable for each extraction method. In samples where CE-based STRs yielded partial profiles (bone), MPS-based STRs generated more complete or full profiles. Moreover, MPS panels contain more STR loci than current CE-based STR kits and also include SNPs, which can further increase the power of discrimination obtained from these samples, making MPS a desirable choice for the forensic analysis of such challenging samples.

Massively Parallel Sequencing of Forensic STRs Using the Ion Chef™ and the Ion S5™ XL Systems

Next-generation sequencing (NGS) has been used to genotype forensic short tandem repeat (STR) markers for individual identification and kinship analysis. STR data from several NGS platforms have been published, but forensic application trials using the Ion S5™ XL system have not been reported. In this work, we report sensitivity, reproducibility, mixture, simulated degradation, and casework sample data on the Ion Chef™ and S5™ XL systems using an early access 25-plex panel. Sensitivity experiments showed that over 97% of the alleles were detectable with down to 62 pg input of genomic DNA. In mixture studies, alleles from minor contributors were correctly assigned at 1:9 and 9:1 ratios. NGS successfully gave 12 full genotype results from 13 challenging casework samples, compared with five full results using the CE platform. In conclusion, the Ion Chef™ and the Ion S5™ XL systems provided an alternative and promising approach for forensic STR genotyping.

A review of bioinformatic methods for forensic DNA analyses

Short tandem repeats, single nucleotide polymorphisms, and whole mitochondrial analyses are three classes of markers which will play an important role in the future of forensic DNA typing. The arrival of massively parallel sequencing platforms in forensic science reveals new information such as insights into the complexity and variability of the markers that were previously unseen, along with amounts of data too immense for analyses by manual means. Along with the sequencing chemistries employed, bioinformatic methods are required to process and interpret this new and extensive data. As more is learnt about the use of these new technologies for forensic applications, development and standardization of efficient, favourable tools for each stage of data processing is being carried out, and faster, more accurate methods that improve on the original approaches have been developed. As forensic laboratories search for the optimal pipeline of tools, sequencer manufacturers have incorporated pipelines into sequencer software to make analyses convenient. This review explores the current state of bioinformatic methods and tools used for the analyses of forensic markers sequenced on the massively parallel sequencing (MPS) platforms currently most widely used.

A technique for setting analytical thresholds in massively parallel sequencing-based forensic DNA analysis

Amplicon (targeted) sequencing by massively parallel sequencing (PCR-MPS) is a potential method for use in forensic DNA analyses. In this application, PCR-MPS may supplement or replace other instrumental analysis methods such as capillary electrophoresis and Sanger sequencing for STR and mitochondrial DNA typing, respectively. PCR-MPS also may enable the expansion of forensic DNA analysis methods to include new marker systems such as single nucleotide polymorphisms (SNPs) and insertion/deletions (indels) that currently are assayable using various instrumental analysis methods including microarray and quantitative PCR. Acceptance of PCR-MPS as a forensic method will depend in part upon developing protocols and criteria that define the limitations of a method, including a defensible analytical threshold or method detection limit. This paper describes an approach to establish objective analytical thresholds suitable for multiplexed PCR-MPS methods. A definition is proposed for PCR-MPS method background noise, and an analytical threshold based on background noise is described.

Clinical diagnostic value of viable Schistosoma japonicum eggs detected in host tissues

Background

Schistosomiasis, one of the neglected tropical diseases, is endemic in more than 70 countries. However, the clinical diagnosis of patients with a low degree of infection is an unsolved technical problem. In areas endemic for schistosomiasis japonica, proctoscopy detection of eggs has been one method used for clinical diagnosis. However, it is often a challenge to find typical live eggs and it is difficult to distinguish live eggs from large numbers of partially degraded and/or completely degraded eggs within colon biopsy tissue. To address this problem, we tested six different morphological and biochemical/molecular markers (ALP; morphological characteristics of egg; CalS (calcified substance); AOS (antioxidase); SDHG (succinic dehydrogenase) and SjR2 mRNA (retrotransposons 2 of S.japonicum genome mRNA)), including four new markers (CalS; AOS; SDHG and SjR2 mRNA.), to determine the viability of S. japonicum eggs deposited in human and mouse colon tissues. Our ultimate aim is to obtain a new method that is more sensitive, practical and accurate to clinically diagnose schistosomiasis.

Methods

Tissue samples were collected from mice at six different time points during S. japonicum infection with or without treatment with praziquantel (PZQ). Four new biochemical or molecular markers were used for the detection of egg viability from mouse liver and intestinal samples: CalS; AOS; SDHG and SjR2 mRNA. Subsequently, all markers were employed for the detection and analysis of eggs deposited in biopsy materials from patients with suspected schistosomiasis japonica for clinical evaluation. Microscopic examination of the egg morphology, worm burden in vivo and ALP (alkaline phosphatase) levels were used as a reference standard to evaluate the sensitivity and reliability of four new markers detecting egg viability.

Results

The results of the study showed that the morphology of S. japonicum eggs deposited in tissues of hosts with schistosomiasis, especially cases with chronic schistosomiasis, is complex and egg viability is difficult to judge morphologically, particularly eggs with a fuzzy structure or partially modified eggs. We found that the majority of the viable schistosome eggs determined by four new markers (CalS, AOS, SDHG and SjR2 mRNA) were morphologically difficult to identify.

Conclusions

Among the markers, the most sensitive and specific method was the detection of SjR2 mRNA and the most simple, rapid and practical method was the detection of SDHG. Therefore, the detection of SDHG is the most practical for clinical application and its use could improve the accuracy in diagnosing active schistosome infection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12879-017-2362-4) contains supplementary material, which is available to authorized users.

Digital fragment analysis of short tandem repeats by high-throughput amplicon sequencing

High‐throughput sequencing has been proposed as a method to genotype microsatellites and overcome the four main technical drawbacks of capillary electrophoresis: amplification artifacts, imprecise sizing, length homoplasy, and limited multiplex capability. The objective of this project was to test a high‐throughput amplicon sequencing approach to fragment analysis of short tandem repeats and characterize its advantages and disadvantages against traditional capillary electrophoresis. We amplified and sequenced 12 muskrat microsatellite loci from 180 muskrat specimens and analyzed the sequencing data for precision of allele calling, propensity for amplification or sequencing artifacts, and for evidence of length homoplasy. Of the 294 total alleles, we detected by sequencing, only 164 alleles would have been detected by capillary electrophoresis as the remaining 130 alleles (44%) would have been hidden by length homoplasy. The ability to detect a greater number of unique alleles resulted in the ability to resolve greater population genetic structure. The primary advantages of fragment analysis by sequencing are the ability to precisely size fragments, resolve length homoplasy, multiplex many individuals and many loci into a single high‐throughput run, and compare data across projects and across laboratories (present and future) with minimal technical calibration. A significant disadvantage of fragment analysis by sequencing is that the method is only practical and cost‐effective when performed on batches of several hundred samples with multiple loci. Future work is needed to optimize throughput while minimizing costs and to update existing microsatellite allele calling and analysis programs to accommodate sequence‐aware microsatellite data.

Human age estimation from blood using mRNA, DNA methylation, DNA rearrangement, and telomere length

Establishing the age of unknown persons, or persons with unknown age, can provide important leads in police investigations, disaster victim identification, fraud cases, and in other legal affairs. Previous methods mostly relied on morphological features available from teeth or skeletal parts. The development of molecular methods for age estimation allowing to use human specimens that possess no morphological age information, such as bloodstains, is extremely valuable as this type of samples is commonly found at crime scenes. Recently, we introduced a DNA-based approach for human age estimation from blood based on the quantification of T-cell specific DNA rearrangements (sjTRECs), which achieves accurate assignment of blood DNA samples to one of four 20-year-interval age categories. Aiming at improving the accuracy of molecular age estimation from blood, we investigated different types of biomarkers. We started out by systematic genome-wide surveys for new age-informative mRNA and DNA methylation markers in blood from the same young and old individuals using microarray technologies. The obtained candidate markers were validated in independent samples covering a wide age range using alternative technologies together with previously proposed DNA methylation, sjTREC, and telomere length markers. Cross-validated multiple regression analysis was applied for estimating and validating the age predictive power of various sets of biomarkers within and across different marker types. We found that DNA methylation markers outperformed mRNA, sjTREC, and telomere length in age predictive power. The best performing model included 8 DNA methylation markers derived from 3 CpG islands reaching a high level of accuracy (cross-validated R(2)=0.88, SE±6.97 years, mean absolute deviation 5.07 years). However, our data also suggest that mRNA markers can provide independent age information: a model using a combined set of 5 DNA methylation markers and one mRNA marker could provide similarly high accuracy (cross-validated R(2)=0.86, SE±7.62 years, mean absolute deviation 4.60 years). Overall, our study provides new and confirms previously suggested molecular biomarkers for age estimation from blood. Moreover, our comparative study design revealed that DNA methylation markers are superior for this purpose over other types of molecular biomarkers tested. While the new and some previous findings are highly promising, before molecular age estimation can eventually meet forensic practice, the proposed biomarkers should be tested further in larger sets of blood samples from both healthy and unhealthy individuals, and markers and genotyping methods shall be validated to meet forensic standards.

Exploring the efficacy of paternity and kinship testing based on single nucleotide polymorphisms

Short tandem repeats (STRs) are conventional genetic markers typically used for paternity and kinship testing. As supplementary markers of STRs, single nucleotide polymorphisms (SNPs) have less discrimination power but broader applicability to degraded samples. The rapid improvement of next-generation sequencing (NGS) and multiplex amplification technologies also make it possible now to simultaneously identify dozens or even hundreds of SNP loci in a single pool. However, few studies have been endeavored to kinship testing based on SNP loci. In this study, we genotyped 90 autosomal human identity SNP loci with NGS, and investigated their testing efficacies based on the likelihood ratio model in eight pedigree scenarios involving paternity, half/full-sibling, uncle/nephew, and first-cousin relationships. We found that these SNPs might be sufficient to discriminate paternity and full-sibling, but impractical for more distant relatives such as uncle and cousin. Furthermore, we conducted an in silico study to obtain the theoretical tendency of how testing efficacy varied with increasing number of SNP loci. For each testing battery in a given pedigree scenario, we obtained distributions of logarithmic likelihood ratio for both simulated relatives and unrelated controls. The proportion of the overlapping area between the two distributions was defined as a false testing level (FTL) to evaluate the testing efficacy. We estimated that 85, 127, 491, and 1,858 putative SNP loci were required to discriminate paternity, full-sibling, half-sibling/uncle-nephew, and first-cousin (FTL, 0.1%), respectively. To test a half-sibling or nephew, an additional uncle relative could be included to decrease the required number of putative SNP loci to ∼320 (FTL, 0.1%). As a systematic computation of paternity and kinship testing based only on SNPs, our results could be informative for further studies and applications on paternity and kinship testing using SNP loci.