Comparison of genotyping and weight of evidence results when applying different genotyping strategies on samples from a DNA transfer experiment.

In this study, we assessed to what extent data on the subject of TPPR (transfer, persistence, prevalence, recovery) that are obtained through an older STR typing kit can be used in an activity-level evaluation for a case profiled with a more modern STR kit. Newer kits generally hold more loci and may show higher sensitivity especially when reduced reaction volumes are used, and this could increase the evidential value at the source level. On the other hand, the increased genotyping information may invoke a higher number of contributors in the weight of evidence calculations, which could affect the evidential values as well. An activity scenario well explored in earlier studies [1,2] was redone using volunteers with known DNA profiles. DNA extracts were analyzed with three different approaches, namely using the optimal DNA input for (1) an older and (2) a newer STR typing system, and (3) using a standard, volume-based input combined with replicate PCR analysis with only the newer STR kit. The genotyping results were analyzed for various aspects such as percentage detected alleles and relative peak height contribution for background and the contributors known to be involved in the activity. Next, source-level LRs were calculated and the same trends were observed with regard to inclusionary and exclusionary LRs for persons who had or had not been in direct contact with the sampled areas. We subsequently assessed the impact on the outcome of the activity-level evaluation in an exemplary case by applying the assigned probabilities to a Bayesian network. We infer that data from different STR kits can be combined in the activity-level evaluations.

Inter-laboratory study on standardized MPS libraries: evaluation of performance, concordance, and sensitivity using mixtures and degraded DNA.

We present results from an inter-laboratory massively parallel sequencing (MPS) study in the framework of the SeqForSTRs project to evaluate forensically relevant parameters, such as performance, concordance, and sensitivity, using a standardized sequencing library including reference material, mixtures, and ancient DNA samples. The standardized library was prepared using the ForenSeq DNA Signature Prep Kit (primer mix A). The library was shared between eight European laboratories located in Austria, France, Germany, The Netherlands, and Sweden to perform MPS on their particular MiSeq FGx sequencers. Despite variation in performance between sequencing runs, all laboratories obtained quality metrics that fell within the manufacturer's recommended ranges. Furthermore, differences in locus coverage did not inevitably adversely affect heterozygous balance. Inter-laboratory concordance showed 100% concordant genotypes for the included autosomal and Y-STRs, and still, X-STR concordance exceeded 83%. The exclusive reasons for X-STR discordances were drop-outs at DXS10103. Sensitivity experiments demonstrated that correct allele calling varied between sequencing instruments in particular for lower DNA amounts (≤ 125 pg). The analysis of compromised DNA samples showed the drop-out of one sample (FA10013B01A) while for the remaining three degraded DNA samples MPS was able to successfully type ≥ 87% of all aSTRs, ≥ 78% of all Y-STRs, ≥ 68% of all X-STRs, and ≥ 92% of all iSNPs demonstrating that MPS is a promising tool for human identity testing, which in return, has to undergo rigorous in-house validation before it can be implemented into forensic routine casework.

Extended specificity studies of mRNA assays used to infer human organ tissues and body fluids.

Messenger RNA (mRNA) profiling is a technique increasingly applied for the forensic identification of body fluids and skin. More recently, an mRNA-based organ typing assay was developed which allows for the inference of brain, lung, liver, skeletal muscle, heart, kidney, and skin tissue. When applying this organ typing system in forensic casework for the presence of animal, rather than human, tissue is an alternative scenario to be proposed, for instance that bullets carry cell material from a hunting event. Even though mRNA profiling systems are commonly in silico designed to be primate specific, physical testing against other animal species is generally limited. In this study, human specificity of the organ tissue inferring system was assessed against organ tissue RNAs of various animals. Results confirm human specificity of the system, especially when utilizing interpretation rules considering multiple markers per cell type. Besides, we cross-tested our organ and body fluid mRNA assays against the target types covered by the other assay. Marker expression in the nontarget organ tissues and body fluids was observed to a limited extent, which emphasizes the importance of involving the case-specific context of the forensic samples in deciding which mRNA profiling assay to use and when for interpreting results.

Development of a control region-based mtDNA SNaPshot™ selection tool, integrated into a mini amplicon sequencing method.

Mitochondrial DNA (mtDNA) analysis is regularly applied to forensic DNA samples with limited amounts of nuclear DNA (nDNA), such as hair shafts and bones. Generally, this mtDNA analysis involves examination of the hypervariable control region by Sanger sequencing of amplified products. When samples are severely degraded, small-sized amplicons can be applied and an earlier described mini-mtDNA method by Eichmann et al. [1] that accommodates ten mini amplicons in two multiplexes is found to be a very robust approach. However, in cases with large numbers of samples, like when searching for hairs with an mtDNA profile deviant from that of the victim, the method is time (and cost) consuming. Previously, Chemale et al. [2] described a SNaPshot™-based screening tool for a Brazilian population that uses standard-size amplicons for HVS-I and HVS-II. Here, we describe a similar tool adapted to the full control region and compatible with mini-mtDNA amplicons. Eighteen single nucleotide polymorphisms (SNPs) were selected based on their relative frequencies in a European population. They showed a high discriminatory power in a Dutch population (97.2%). The 18 SNPs are assessed in two SNaPshot™ multiplexes that pair to the two mini-mtDNA amplification multiplexes. Degenerate bases are included to limit allele dropout due to SNPs at primer binding site positions. Three SNPs provide haplogroup information. Reliability testing showed no differences with Sanger sequencing results. Since mini-mtSNaPshot screening uses only a small portion of the same PCR products used for Sanger sequencing, no additional DNA extract is consumed, which is forensically advantageous.

Evaluation of samples comprising minute amounts of DNA.

Minute amounts of DNA representing only few diploid cells, may be interrogated using enhanced DNA profiling, which will be accompanied by stochastic amplification effects. Notwithstanding, a weight of evidence statistic may be calculated using current interpretation software. In this study, we profiled single donor, two- and three-person samples having only 3 pg to 12 pg of DNA per contributor using both standard and enhanced capillary electrophoresis (CE) injection settings. Likelihood ratios (LRs) were computed using LRmix Studio, compared for both types of profiles and examined in relation to the amount of DNA, drop-out level, number of detected alleles, peak heights and reproducibility of alleles. Especially for DNA profiles that were generated using enhanced CE, the obtained LRs could indicate strong evidence in favour of the prosecution (log10(LR)>6), also when the amount of DNA represented about half of a diploid cell equivalent in the amplification. These results illustrate that an assessment of the criminalistic relevance of a sample carrying minute amounts of DNA is essential prior to applying enhanced interrogation techniques and/or calculating a weight of evidence statistic.

Oxidation monitoring by fluorescence spectroscopy reveals the age of fingermarks.

No forensic method exists that can reliably estimate the age of fingermarks found at a crime scene. Information on time passed since fingermark deposition is desired as it can be used to distinguish between crime related and unrelated fingermarks and to support or refute statements made by the fingermark donors. We introduce a non-contact method that can estimate the age of fingermarks. Fingermarks were approached as protein-lipid mixtures and an age-estimation model was build based on the expected protein and lipid oxidation reactions. Two measures of oxidation are required from the fingermark to estimate its age: 1) the relative amount of fluorescent oxidation products 2) the rate at which these products are formed. Fluorescence spectroscopy was used to obtain these measures. We tested the method on 44 fingermarks and were able to estimate the age of 55% of the male fingermarks, up to three weeks old with an uncertainty of 1.9 days.

Consensus and pool profiles to assist in the analysis and interpretation of complex low template DNA mixtures.

Forensic analysis of low template (LT) DNA mixtures is particularly complicated when (1) LT components concur with high template components, (2) more than three contributors are present, or (3) contributors are related. In this study, we generated a set of such complex LT mixtures and examined two methods to assist in DNA profile analysis and interpretation: the "n/2" consensus method (Benschop et al. 2011) and the pool profile approach. N/2 consensus profiles include alleles that are reproducibly amplified in at least half of the replications. Pool profiles are generated by injecting a blend of independently amplified PCR products on a capillary electrophoresis instrument. Both approaches resulted in a similar increase in the percentage of detected alleles compared to individual profiles, and both rarely included drop-in alleles in case mixtures of pristine DNAs were used. Interestingly, the consensus and the pool profiles often showed differences for the actual alleles detected for the LT component(s). We estimated the number of contributors using different methods. Better approximations were obtained with data in the consensus and pool profiles compared to the data of the individual profiles. Consensus profiles contain allele calls only, while pool profiles consist of both allele calls and peak height information, which can be of use in (statistical) profile analysis. All advantages and limitations of the various types of profiles were assessed, and based on the results we infer that both consensus and pool profiles (or a combination thereof) are helpful in the interpretation of complex LT DNA mixtures.

Improved analysis of long STR amplicons from degraded single source and mixed DNA.

DNA profiles from degraded samples often suffer from information loss at the longer short tandem repeat (STR) loci. Sensitising the reactions, either by performing additional PCR cycles or increasing the capillary electrophoresis injection settings, carries the risk of over-amplifying or overloading the shorter fragments. We explored whether profiling of degraded DNA can be improved by preferential capturing of the longer amplified fragments. To this aim, a post-PCR purification protocol was developed that is based on AMPure XP beads that have size-selective properties. A comparison was made with an unselective post-PCR purification system (DTR gel filtration) and no purification of the PCR products. Besides a set of differently and serially degraded single source samples, unequal mixtures of degraded DNAs were analysed, in order to extract more genotyping information for the minor contributor without overloading the major component at the shorter amplicons. Purification by the AMPure protocol resulted in higher peak heights especially for the longer amplicons, while DTR gel filtration gave higher peaks for all amplicon sizes. Both purification methods presented more detected alleles, with the AMPure protocol performing slightly better, on average. In conclusion, the in-house developed AMPure protocol can be employed to improve STR profile analysis of degraded single source and (unequally) mixed DNA samples.

Development of a mRNA profiling multiplex for the inference of organ tissues.

Forensic characterisation of organ tissue generally occurs through histological and immunological assays of limited sensitivity. Here, we explore an alternative approach and examine a total of 41 candidate mRNA markers for their ability to differentiate between brain, lung, liver, skeletal muscle, heart, kidney and skin. Various selection rounds are applied involving 85 organ tissues (36 excised autopsy specimens and 49 frozen tissue sections, with at least ten specimens for each organ type), 20 commercially available RNAs from different human tissues and at least two specimens of blood, saliva, semen, vaginal mucosa, menstrual secretion or touch samples. Finally, 14 markers are regarded tissue-specific and included in an endpoint RT-PCR multiplex together with one general muscle, one blood and one housekeeping marker. This 17-plex is successfully used to analyse a blind test set of 20 specimens including mixtures, and samples derived from stabbing of organ tissues. With the blind test set samples, it is shown that an earlier described interpretation strategy for RNA cell typing results [1] is also effective for tissue inference. As organ-typing is embedded in a procedure of combined DNA/RNA extraction and analysis, both donor and organ type information is derived from the same sample. Some autopsy specimens presented DNA profiles characteristic for degraded DNA. Nevertheless, the organ-typing multiplex could generate full RNA profiles, which is probably due to small sizes of the amplicons. This assay provides a novel tool for analysis of samples from violent crimes.

RNA cell typing and DNA profiling of mixed samples: can cell types and donors be associated?

Forensic samples regularly involve mixtures, which are readily recognised in forensic analyses. Combined DNA and mRNA profiling is an upcoming forensic practice to examine donors and cell types from the exact same sample. From DNA profiles individual genotypes may be deconvoluted, but to date no studies have established whether the cell types identified in corresponding RNA profiles can be associated with individual donors. Although RNA expression levels hold many variables from which an association may not be expected, proof of concept is important to forensic experts who may be cross examined about this possible correlation in court settings. Clearly, the gender-specificity of certain body fluids (semen, vaginal mucosa, menstrual secretion) can be instructive. However, when donors of the same gender or gender-neutral cell types are involved, alternatives are needed. Here we analyse basic two-component mixtures (two cell types provided by different donors) composed of six different cell types, and assess whether the heights of DNA and RNA peaks may guide association of donor and cell type. Divergent results were obtained; for some mixtures RNA peak heights followed the DNA results, but for others the major DNA component did not present higher RNA peaks. Also, variation in mixture ratios was observed for RNA profiling replicates and when different donor couples gave the same two body fluids. As sample degradation may affect the two nucleic acids and/or distinct cell types differently (and thus influence donor and cell type association), mixtures were subjected to elevated temperature or UV-light. Variation in DNA and RNA stability was observed both between and within cell types and depended on the method inducing degradation. Taken together, we discourage to associate cell types and donors from peak heights when performing RNA and DNA profiling.

ProbRank: An efficient DNA database search method for complex mixtures per a quantitative likelihood ratio model.

Searching a DNA Database with a DNA profile from an evidentiary trace can provide investigative leads in a forensic case. Various searching approaches exist such as conventional methods based on matching alleles or more advanced methods computing likelihood ratios (LR) while considering drop-in and drop-out. Here we examine the potential of using a quantitative LR model (EuroForMix model incorporated in ProbRank method) that takes peak heights into account in comparison to a qualitative LR model (LRmix model implemented in SmartRank method). Both methods present DNA database candidates in order of decreasing LR. Especially regarding minor contributors in complex mixtures, the method using the quantitative model outperforms the method using the qualitative model in terms of sensitivity and specificity as more true donors and less adventitious matches are retrieved. ProbRank is to be implemented in DNAStatistX and is sufficiently fast for daily use.

Targeted DNA methylation analysis and prediction of smoking habits in blood based on massively parallel sequencing.

Tobacco smoking is a frequent habit sustained by > 1.3 billion people in 2020 and the leading preventable factor for health risk and premature mortality worldwide. In the forensic context, predicting smoking habits from biological samples may allow broadening DNA phenotyping. In this study, we aimed to implement previously published smoking habit classification models based on blood DNA methylation at 13 CpGs. First, we developed a matching lab tool based on bisulfite conversion and multiplex PCR followed by amplification-free library preparation and targeted paired-end massively parallel sequencing (MPS). Analysis of six technical duplicates revealed high reproducibility of methylation measurements (Pearson correlation of 0.983). Artificially methylated standards uncovered marker-specific amplification bias, which we corrected via bi-exponential models. We then applied our MPS tool to 232 blood samples from Europeans of a wide age range, of which 90 were current, 71 former and 71 never smokers. On average, we obtained 189,000 reads/sample and 15,000 reads/CpG, without marker drop-out. Methylation distributions per smoking category roughly corresponded to previous microarray analysis, showcasing large inter-individual variation but with technology-driven bias. Methylation at 11 out of 13 smoking-CpGs correlated with daily cigarettes in current smokers, while solely one was weakly correlated with time since cessation in former smokers. Interestingly, eight smoking-CpGs correlated with age, and one displayed weak but significant sex-associated methylation differences. Using bias-uncorrected MPS data, smoking habits were relatively accurately predicted using both two- (current/non-current) and three- (never/former/current) category model, but bias correction resulted in worse prediction performance for both models. Finally, to account for technology-driven variation, we built new, joint models with inter-technology corrections, which resulted in improved prediction results for both models, with or without PCR bias correction (e.g. MPS cross-validation F1-score > 0.8; 2-categories). Overall, our novel assay takes us one step closer towards the forensic application of viable smoking habit prediction from blood traces. However, future research is needed towards forensically validating the assay, especially in terms of sensitivity. We also need to further shed light on the employed biomarkers, particularly on the mechanistics, tissue specificity and putative confounders of smoking epigenetic signatures.

Vaginal microbial flora analysis by next generation sequencing and microarrays; can microbes indicate vaginal origin in a forensic context?

Forensic analysis of biological traces generally encompasses the investigation of both the person who contributed to the trace and the body site(s) from which the trace originates. For instance, for sexual assault cases, it can be beneficial to distinguish vaginal samples from skin or saliva samples. In this study, we explored the use of microbial flora to indicate vaginal origin. First, we explored the vaginal microbiome for a large set of clinical vaginal samples (n = 240) by next generation sequencing (n = 338,184 sequence reads) and found 1,619 different sequences. Next, we selected 389 candidate probes targeting genera or species and designed a microarray, with which we analysed a diverse set of samples; 43 DNA extracts from vaginal samples and 25 DNA extracts from samples from other body sites, including sites in close proximity of or in contact with the vagina. Finally, we used the microarray results and next generation sequencing dataset to assess the potential for a future approach that uses microbial markers to indicate vaginal origin. Since no candidate genera/species were found to positively identify all vaginal DNA extracts on their own, while excluding all non-vaginal DNA extracts, we deduce that a reliable statement about the cellular origin of a biological trace should be based on the detection of multiple species within various genera. Microarray analysis of a sample will then render a microbial flora pattern that is probably best analysed in a probabilistic approach.

Combining results of forensic STR kits: HDplex validation including allelic association and linkage testing with NGM and Identifiler loci.

The autosomal short tandem repeat (STR) kits that are currently used in forensic science have a high discrimination power. However, this discrimination power is sometimes not sufficient for complex kinship analyses or decreases when alleles are missing due to degradation of the DNA. The Investigator HDplex kit contains nine STRs that are additional to the commonly used forensic markers, and we validated this kit to assist human identification. With the increasing number of markers it becomes inevitable that forensic and kinship analyses include two or more STRs present on the same chromosome. To examine whether such markers can be regarded as independent, we evaluated the 30 STRs present in NGM, Identifiler and HDplex. Among these 30 markers, 17 syntenic STR pairs can be formed. Allelic association between these pairs was examined using 335 Dutch reference samples and no linkage disequilibrium was detected, which makes it possible to use the product rule for profile probability calculations in unrelated individuals. Linkage between syntenic STRs was studied by determining the recombination fraction between them in five three-generation CEPH families. The recombination fractions were compared to the physical and genetic distances between the markers. For most types of pedigrees, the kinship analyses can be performed using the product rule, and for those cases that require an alternative calculation method (Gill et al., Forensic Sci Int Genet 6:477-486, 2011), the recombination fractions as determined in this study can be used. Finally, we calculated the (combined) match probabilities, for the supplementary genotyping results of HDplex, NGM and Identifiler.

Structural features of small RNA precursors determine Argonaute loading in Caenorhabditis elegans.

In C. elegans, DCR-1 is required for the maturation of both short interfering RNAs (siRNAs) and microRNAs (miRNAs), which are subsequently loaded into different Argonaute proteins to mediate silencing via distinct mechanisms. We used in vivo analyses to show that precursors of small RNAs contain structural features that direct the small RNAs into the RNA interference (RNAi) pathway or the miRNA-processing pathway. Nucleotide changes in the pre-let-7 miRNA precursor that make its stem fully complementary cause the resulting small RNA to be recognized as siRNA and induce binding to RDE-1, which leads to RNAi. Mismatches of 1 to 3 nucleotides at various positions in the stem of the precursor restore direction into the miRNA pathway, as the largest portion of such small RNA variants is associated with ALG-1. The Argonaute proteins to which the small RNAs are bound determine the silencing mode, and no functional overlap between RDE-1 and ALG-1 was detected.

Advancing FDSTools by integrating STRNaming 1.1.

The introduction of massively parallel sequencing in forensic analysis has been facilitated with typing kits, analysis software and allele naming tools such as the ForenSeq DNA Signature Prep (DSP) kit, FDSTools and STRNaming respectively. Here we describe how FDSTools 2.0 with integrated and refined STRNaming nomenclature was validated for implementation under ISO 17025 accreditation for the ForenSeq DSP kit. Newly-added options result in efficient automatic allele calling for the majority of markers while specific settings are applied for 'novel' sequence variants to avoid the calling of remaining variable noise observed in samples sequenced with the ForenSeq DSP kit that seem to arise in the PCR. Genome-wide built-in reference data allows for greatly simplified configuration of allele naming for human targets.

Experiences in successive forensic analysis of GSR and DNA traces collected on a single tape lift stub.

Gunshot residue (GSR) and deoxyribonucleic acid (DNA) analyses are used in the reconstruction of shooting incidents involving firearms. In both fields, adhesive tapes are commonly applied to recover traces. Using a simultaneous sampling approach for retrieving both types of traces from objects related to a shooting can be powerful. We present the results for DNA testing and GSR analysis using "DNA stubs" to recover both types of traces in one sampling. The analyses are successive: stubs are first examined with scanning electron microscopy for the presence of GSR and next subjected to DNA extraction and short tandem repeat (STR) typing. The results show no negative effect on the DNA and GSR results. The risk of DNA contamination during GSR analysis is addressed, showing that it is negligible. Our experiences in casework show that the approach can be successfully implemented.

Evaluation of the VISAGE basic tool for appearance and ancestry inference using ForenSeq® chemistry on the MiSeq FGx® system.

The possibility of providing investigative leads when conventional DNA identification methods fail to solve a case can be of extreme relevance to law enforcement. Therefore, the forensic genetics community has focused research towards the broadened use of DNA, particularly for prediction of appearance traits, bio-geographical ancestry and age. The VISible Attributes through GEnomics (VISAGE) Consortium expanded the use of DNA phenotyping by developing new molecular and statistical tools for appearance, age and ancestry prediction. The VISAGE basic tool for appearance (EVC) and ancestry (BGA) prediction was initially developed using Ampliseq chemistry, but here is being evaluated using ForenSeq chemistry. The VISAGE basic tool offers a total of 41 EVC and 115 BGA SNPs and thus provides more predictions, i.e., skin color, than achieved with the ForenSeq DNA Signature Prep kit that is based on 24 EVC and 56 BGA SNPs. Five VISAGE laboratories participated in collaborative experiments to provide foreground for developmental validation of the assay. Assessment of assay performance and quality metrics, reproducibility, sensitivity, inhibitor tolerance and species specificity are described. Furthermore, the assay was tested using challenging samples such as mock casework samples and artificially degraded DNA. Two different analysis strategies were applied for this study and output on genotype calls and read depth was compared. Overall, inter-laboratory, inter-method and concordance with publicly available data were analysed and compared. Finally, the results showed a reliable and robust tool, which can be easily applied for laboratories already using a MiSeq FGx with ForenSeq reagents.

On the Identification of Body Fluids and Tissues: A Crucial Link in the Investigation and Solution of Crime.

Body fluid and body tissue identification are important in forensic science as they can provide key evidence in a criminal investigation and may assist the court in reaching conclusions. Establishing a link between identifying the fluid or tissue and the DNA profile adds further weight to this evidence. Many forensic laboratories retain techniques for the identification of biological fluids that have been widely used for some time. More recently, many different biomarkers and technologies have been proposed for identification of body fluids and tissues of forensic relevance some of which are now used in forensic casework. Here, we summarize the role of body fluid/ tissue identification in the evaluation of forensic evidence, describe how such evidence is detected at the crime scene and in the laboratory, elaborate different technologies available to do this, and reflect real life experiences. We explain how, by including this information, crucial links can be made to aid in the investigation and solution of crime.

STRNaming: Generating simple, informative names for sequenced STR alleles in a standardised and automated manner.

The introduction of Massively Parallel Sequencing in the forensic domain has exposed the need for comprehensive nomenclature of sequenced Short Tandem Repeat (STR) alleles. In general, three strategies are at hand: 1) the full sequence mapped to the human genome reference sequence, which ensures exact data exchange; 2) shortened, human-readable formats for forensic reporting and data presentation and 3) very short codes that enable compact figures and tables but do not convey any sequence information. Here, we describe an algorithm of the second type: STRNaming, which generates human-readable names for sequenced STR alleles. STRNaming is guided by a reference sequence at each locus and then functions independently to automatically assign a unique, sequence-descriptive name that also includes the capillary electrophoresis allele number. STRNaming settings were established based on preferences that were surveyed internationally in the forensic community. These settings ensure that a small change in the sequence corresponds to a small change in the allele name, which is helpful for recognising for instance stutter products. Sequence variants outside of the repeat units are indicated as simple variant calls. Since the STR name is sequence-descriptive, the sequence can be traced back from the allele name. Because STRNaming is fully guided by an assignable reference sequence, no central coordination or configuration is required and the method will work for any STR locus, be it autosomal, Y-, X-chromosomal in current or future use. The algorithm is publicly available online and offline.