Evaluation of reverse transcription yield of RNA standards and forensic samples based on droplet digital PCR

RNA analysis has shown great value in forensic science, such as body fluids and tissue identification, postmortem interval estimation, biological age prediction, etc. Currently, most RNA follow-up experiments involve reverse transcription (RT) procedures. It has been shown that the RT step is variable and has a greater impact on subsequent data analysis, especially for forensic trace samples. However, the pattern of variation between different RNA template inputs and complementary DNA (cDNA) yield is unclear. In this study, a series of 2-fold gradient dilutions of RNA standards (1 μg/μL - 0.24 ng/μL) and forensic samples (including blood samples, saliva samples, bloodstains, and saliva stains) were reverse-transcribed using EasyQuick RT MasterMix. The obtained cDNA was quantified by droplet digital PCR (ddPCR) to assess the RT yield of the ACTB gene. The results showed that the 125 ng RNA template had the highest RT yield in a 10 μL RT reaction system with the selected kit. For all stain samples, the RT yield improved as the amount of RNA template input increased since RNA quantities were below 125 ng. As many commercialized reverse transcription kits using different kinds of enzymes are available for forensic RNA research, we recommend that systematic experiments should be performed in advance to determine the amount of RNA input at the optimum RT yield when using any kit for reverse transcription experiments.

A panel of tetranucleotide STR markers as an alternative approach to forensic DNA identification of wolf and dog

Commercial panels of microsatellite (STR) loci are focused on the use of DNA of the domestic dog (Canis lupus familiaris) and are often inapplicable for genotyping the DNA of the gray wolf (Canis lupus lupus). We propose a CPlex test system, including one hexa- and 12 tetranucleotide autosomal STR loci, as well as two sex loci, that is equally efficient in DNA identification of biological samples of the wolf and the dog. Analysis of molecular variance between samples revealed significant differentiation values (FST = 0.0784, p < 0.001), which allows to use the panel to differentiate wolf and dog samples. Population subdivision coefficients (θ-values) were calculated for each of the 13 STR loci of the developed test system. It was shown that the values of the genotype frequency for dogs and wolves, without and with considering the θ-value, differ by three orders of magnitude (for dogs 8.9 × 10-16 and 2.1 × 10-14 and for wolves 1.9 × 10-15 and 4.5 × 10-14, respectively). The use of population subdivision coefficients will allow to identify the most reliable results of an expert identification study and the power of exclusion provided by the STR loci of the CPlex test system makes it possible to achieve a reliable level of evidence in forensic DNA analysis of both wolves and dogs. The test system has been validated for use in forensic identification of the dog and wolf based on biological traces found at crime scenes, as well as for individual identification and establishing biological relationship of animals of these species.

Cross-species amplification from non-human primate DNAs in commercial human DNA assay kits

Species specificity of commercial human DNA quantification kits and short tandem repeat (STR) profiling kits was examined using primate DNA samples. These samples comprised 33 individuals from eight primate species, each with gender and kinship data, including human (Homo sapiens), chimpanzee (Pan troglodytes), gorilla (Gorilla gorilla), and orangutan (Pongo pygmaeus) of Hominidae family, and Japanese macaque (Macaca fuscata), long-tailed macaque (Macaca fascicularis), hamadryas baboon (Papio hamadryas), and savannah monkey (Chlorocebus sp.) of Cercopithecidae family. The findings revealed varying levels of cross-species amplifications in all non-human DNA samples that correlated with their evolutionary proximity to humans, both kit types. Moreover, cross-species amplification, including female DNA samples, was observed in a Y-chromosomal STR profiling kit. Additionally, species specificity differed among the commercial kits examined. The cross-species amplification data presented in this study offer valuable assistance in interpreting the results of individual human identification in forensic cases involving non-human primates.

Unveiling facial kinship: The BioKinVis dataset for facial kinship verification and genetic association studies

Facial image-based kinship verification represents a burgeoning frontier within the realms of computer vision and biomedicine. Recent genome-wide association studies have underscored the heritability of human facial morphology, revealing its predictability based on genetic information. These revelations form a robust foundation for advancing facial image-based kinship verification. Despite strides in computer vision, there remains a discernible gap between the biomedical and computer vision domains. Notably, the absence of family photo datasets established through biological paternity testing methods poses a significant challenge. This study addresses this gap by introducing the biological kinship visualization dataset, encompassing 5773 individuals from 2412 families with biologically confirmed kinship. Our analysis delves into the distribution and influencing factors of facial similarity among parent-child pairs, probing the potential association between forensic short tandem repeat polymorphisms and facial similarity. Additionally, we have developed a machine learning model for facial image-based kinship verification, achieving an accuracy of 0.80 in the dataset. To facilitate further exploration, we have established an online tool and database, accessible at http://120.55.161.230:88/.

Genetic diversity and forensic statistical support for the 12 X-STR markers in the Malaysian Indian population using Qiagen Investigator® Argus X-12 QS kit

X-chromosome short tandem repeats (X-STRs) are useful for human identification, especially in complex kinship scenarios. Since forensic statistical parameters vary among populations and the X-STRs population data for the diverse population of Peninsular Malaysia's are unavailable, this attempt for Indians (n = 201) appears forensically relevant to support the 12 X-STRs markers' evidential value for human identification in Malaysia. The Qiagen Investigator® Argus X-12 QS kit showed that DXS10135 was the most polymorphic locus with high genetic diversity, polymorphism information richness, heterozygosity, and exclusion power. Based on allele frequencies, the strength of discrimination and mean exclusion chance (MECKrüger, MECKishida, MECDesmarais, and MECDesmaraisDuo) values for the Malaysian Indians were ≥0.999997790686228. As for haplotype frequencies, the overall discrimination power and mean exclusion probability (MECKrüger, MECKishida, MECDesmarais, and MECDesmaraisDuo) were ≥0.9999984801951. The genetic distance, neighbor-joining phylogenetic tree, and principal component analysis also supported the evidential value of the 12 X-STRs markers for forensic practical caseworks in Malaysia.

Strategies to deal with genetic analyzer-specific DNA methylation measurements

Targeted bisulfite sequencing using single-base extension (SBE) can be used to measure DNA methylation via capillary electrophoresis on genetic analyzers in forensic labs. Several accurate age prediction models have been reported using this method. However, using different genetic analyzers with different software settings can generate different methylation values, leading to significant errors in age prediction. To address this issue, the study proposes and compares four methods as follows: (1) adjusting methylation values using numerous actual body fluid DNA samples, (2) adjusting methylation values using control DNAs with varying methylation ratios, (3) constructing new age prediction models for each genetic analyzer type, and (4) constructing new age prediction models that could be applied to all types of genetic analyzers. To test the methods for adjusting values using actual body fluid DNA samples, previously reported adjusting equations were used for blood/saliva DNA age prediction markers (ELOVL2, FHL2, KLF14, MIR29B2CHG/C1orf132, and TRIM59). New equations were generated for semen DNA age prediction markers (TTC7B, LOC401324/cg12837463, and LOC729960/NOX4) by drawing polynomial regression lines between the results of the three types of genetic analyzers (3130, 3500, and SeqStudio). The same method was applied to obtain adjustment equations using 11 control DNA samples. To develop new age prediction models for each genetic analyzer type, linear regression analysis was conducted using DNA methylation data from 150 blood, 150 saliva, and 62 semen samples. For the genetic analyzer-independent models, control DNAs were used to formulate equations for calibrating the bias of the data from each genetic analyzer, and linear regression analysis was performed using calibrated body fluid DNA data. In the comparison results, the genetic analyzer-specific models showed the highest accuracy. However, genetic analyzer-independent models through bias adjustment also provided accurate age prediction results, suggesting its use as an alternative in situations with multiple constraints.

An examination of differences in epigenetic methylation of saliva type samples based on collection method

In the context of forensic casework, it is imperative to both establish a DNA profile from biological specimens and accurately identify the specific bodily fluid source. To achieve this, DNA methylation markers have been developed for the differentiation of blood, semen, vaginal epithelial secretions, and saliva samples. Saliva, alternatively referred to as oral fluid, is recognized for its heterogeneous cellular composition, characterized by a mixture of epithelial, leukocytic, and bacterial cells. Consequently, our research has revealed variations in methylation percentages that correlate with the method employed for collecting saliva samples. To investigate these concepts, we scrutinized four CpG markers situated within or in proximity to the BCAS4, SLC12A8, SOX2OT, and FAM43A genes. Subsequently, we designed primers based on bioinformatically transformed reference sequences for these markers and rigorously assessed their quality by examining dimer and hairpin formation, melting temperature, and specificity. These loci were identified as saliva markers based on either buccal swabs or spit collection. Yet, there has been minimal or no research conducted to explore the variations in methylation between different collection methods. For this study, buccal, lip, tongue, spit, and nasal swabs were collected from 20 individuals (N = 100). Mock forensic samples, which include chewing gum (N = 10) and cigarettes (N = 10), were also tested. DNA was extracted, bisulfite converted, then amplified using in-house designed assays, and pyrosequenced. The methylation levels were compared to other body fluids (semen, blood, vaginal epithelia, and menstrual blood [N = 32]). A total of 608 pyrosequencing results demonstrated that sampling location and collection method can greatly influence the level of methylation, highlighting the importance of examining multiple collection/deposition methods for body fluids when developing epigenetic markers.

(Un)Reliable detection of menstrual blood in forensic casework - evaluation of the Seratec® PMB test with mock samples

The identification of the type of body fluid in crime scene evidence may be crucial, so that the efforts are high to reduce the complexity of these analyses and to minimize time and costs. Reliable immunochromatographic rapid tests for specific and sensitive identification of blood, saliva, urine and sperm secretions are already routinely used in forensic genetics. The recently introduced Seratec® PMB test is said to detect not only hemoglobin, but also differentiate menstrual blood from other secretions containing blood (cells) by detecting D-dimers. In our experimental set-up, menstrual blood could be reliably detected in mock forensic samples. Here, the result was independent of sample age and extraction buffer volume. It was also successfully demonstrated that all secretions without blood cells were negative for both, hemoglobin (P) and D-dimer (M). However, several blood cell-containing secretions/tissues comprising blood (injury), nasal blood, postmortem blood and wound crust also demonstrated positive results for D-dimer (M) and were therefore false positives. For blood (injury) and nasal blood, this result was reproduced for different extraction buffer volumes. The results of this study clearly demonstrate that the Seratec® PMB test is neither useful nor suitable for use in forensic genetics because of the great risk of false positive results which can lead to false conclusions, especially in sexual offense or violent acts.

Effects and considerations of multiplexing ForenSeq Kintelligence libraries with a negative control

The negative template control or negative amplification control has been an essential component of the forensic DNA analysis workflow that helps monitor contamination. As such, the inclusion of a negative control in forensic DNA analysis has been a requirement for all laboratories audited under the FBI's Quality Assurance Standards. As massively parallel sequencing (MPS) becomes more conventional in forensic laboratories, considerations for the inclusion of a negative control in every sequencing run can be evaluated. Although the inclusion of a negative control in library preparation and the first sequencing run has a practical function, there is less utility for its inclusion in all subsequent sequencing runs for that library preparation. Although this is universal to all MPS assays, it is most relevant for an assay that has a low sample multiplexing capacity, such as the ForenSeq Kintelligence Kit (Qiagen/Verogen, Inc.). The ForenSeq Kintelligence Kit is an investigative genetic genealogy (IGG) sequencing-based assay that targets 10,230 forensically relevant single-nucleotide polymorphisms. The manufacturer recommends multiplexing 3 libraries per sequencing run, which includes controls. The purpose of this study was to investigate the effect of the inclusion of a negative control in every Kintelligence sequencing run. We observed that the library generated from a negative amplification control will take 7%-14% of the run output. The loss of sequencing space taken by a negative control decreased the available output for DNA-containing samples, leading in some cases to allele or locus dropout and accompanying higher numbers of sixth to seventh order unknown associations in GEDmatch PRO.

Pilot validation of on-field STR typing and human identity testing by MinION nanopore sequencing

Nanopore sequencing technology has broad application prospects in forensic medicine due to its small size, portability, fast speed, real-time result analysis capabilities, single-molecule sequencing abilities, and simple operation. Here, we demonstrate for the first time that nanopore sequencing platforms can be used to identify individuals in the field. Through scientific and reasonable design, a nanopore MinION MK1B device and other auxiliary devices are integrated into a portable detection box conducive to individual identification at the accident site. Individual identification of 12 samples could be completed within approximately 24 h by jointly detecting 23 short tandem repeat (STR) loci. Through double-blinded experiments, the genotypes of 49 samples were successfully determined, and the accuracy of the STR genotyping was verified by the gold standard. Specifically, the typing success rate for 1150 genotypes was 95.3%, and the accuracy rate was 86.87%. Although this study focused primarily on demonstrating the feasibility of full-process testing, it can be optimistically predicted that further improvements in bioinformatics workflows and nanopore sequencing technology will help enhance the feasibility of Oxford Nanopore Technologies equipment for real-time individual identification at accident sites.

Emerging use of air eDNA and its application to forensic investigations - A review

Biological material is routinely collected at crime scenes and from exhibits and is a key type of evidence during criminal investigations. Improvements in DNA technologies allow collection and profiling of trace samples, comprised of few cells, significantly expanding the types of exhibits targeted for DNA analysis to include touched surfaces. However, success rates from trace and touch DNA samples tend to be poorer compared to other biological materials such as blood. Simultaneously, there have been recent advances in the utility of environmental DNA collection (eDNA) in identification and tracking of different biological organisms and species from bacteria to naked mole rats in different environments, including, soil, ice, snow, air and aquatic. This paper examines the emerging methods and research into eDNA collection, with a special emphasis on the potential forensic applications of human DNA collection from air including challenges and further studies required to progress implementation.

Determining lineages between individuals with high-density mitochondrial and Y-chromosomal single-nucleotide polymorphisms

Genetic genealogy has been more frequently used in forensic investigations in identifying criminals. However, the current genetic genealogy applications usually do not consider lineage markers (including both Y and mitochondrial deoxyribonucleic acid (DNA)), which is probably because not all distant relatives share the same lineage markers. In addition, there is no study to show how to use lineage markers and what methods or thresholds should be applied in genetic genealogy. In this study, we developed a method to quickly determine if two single-nucleotide polymorphism (SNP) profiles are from the same paternal or material lineages by using a mismatch frequency of the SNPs in Y-chromosomal or mitochondrial DNA. For both Y and mitochondrial SNPs, profile pairs from the same or different lineages can be decided with high accuracies (i.e., 0.380% or 0.157% error rates with Y and mitochondrial DNA, respectively). With kinship coefficient filtering based on autosomal SNPs, the accuracies of determining maternal and paternal lineage can be further improved (i.e., 0.120% or 0.057% error rates with Y and mitochondrial DNA, respectively, using a threshold of kinship coefficient >0). This study shows that lineage markers can be very powerful tools with high accuracies to determine lineages, which could help solve cases and reduce costs for genetic genealogy investigations.

Metal-DNA interactions: Exploring the impact of metal ions on key stages of forensic DNA analysis

Forensic DNA analysis continues to be hampered by the complex interactions between metals and DNA. Metal ions may cause direct DNA damage, inhibit DNA extraction and polymerase chain reaction (PCR) amplification or both. This study evaluated the impact of metal ions on DNA extraction, quantitation, and short tandem repeat profiling using cell-free and cellular (saliva) DNA. Of the 11 metals assessed, brass exhibited the strongest PCR inhibitory effects, for both custom and Quantifiler Trio quantitation assays. Metal ion inhibition varied across the two quantitative PCR assays and the amount of DNA template used. The Quantifiler Trio internal PCR control (IPC) only revealed evidence of PCR inhibition at higher metal ion concentrations, limiting the applicability of IPC as an indicator of the presence of metal inhibitor in a sample. Notably, ferrous ions were found to significantly decrease the extraction efficiency of the DNA-IQ DNA extraction system. The amount of DNA degradation and inhibition in saliva samples caused by metal ions increased with a dilution of the sample, suggesting that the saliva matrix provides protection from metal ion effects.

Development of a novel panel for blood identification based on blood-specific CpG-linked SNP markers

Blood-containing mixtures often appear in murder and robbery cases, and their identification plays a significant role in solving crimes. In recent years, the co-detection of DNA methylation markers (CpG) and single nucleotide polymorphism (SNP) markers has been shown to be a promising tool for the identification of semen and its donor. However, similar research on blood stains that are frequently found at crime scenes has not yet been reported. In this study, we employed blood-specific CpG-linked SNP markers (CpG-SNP) for blood-specific genotyping and the linking of blood and its donor. The tissue-specific CpG markers were screened from the literature and further verified by combining bisulfite conversion with amplification-refractory mutation system (ARMS) technology. Meanwhile, adjacent SNP markers with a minor allele frequency (MAF) greater than 0.1 were selected within 400 bp upstream and downstream of the CpG markers. SNP genotyping was performed using SNaPshot technology on a capillary electrophoresis (CE) platform. Finally, a multiplex panel, including 19 blood-specific CpG linked to 23 SNP markers, as well as 1 semen-specific CpG, 1 vaginal secretion-specific CpG, and 1 saliva-specific CpG marker, was constructed successfully. The panel showed good tissue specificity and blood stains stored at room temperature for up to nine months and moderately degraded (4 < DI < 10) could be effectively identified. Moreover, it could also be detected when blood content in the mixed stains was as low as 1%. In addition, 15 ng of DNA used for bisulfite conversion was required for obtaining a complete profile. The cumulative discrimination power of the panel among the Han population of northern China could reach 0.999983. This is the first investigation conducted for the simultaneous identification of blood and its donor regardless of other body fluids included in mixed stains. The successful construction of the panel will play a vital role in the comprehensive analysis of blood-containing mixtures in forensic practice.

A developmental validation of the Quick TargSeq 1.0 integrated system for automated DNA genotyping in forensic science for reference samples

Analysis of short tandem repeats (STRs) is a global standard method for human identification. Insertion/Deletion polymorphisms (DIPs) can be used for biogeographical ancestry inference. Current DNA typing involves a trained forensic worker operating several specialized instruments in a controlled laboratory environment, which takes 6-8 h. We developed the Quick TargSeq 1.0 integrated system (hereinafter abbreviated to Quick TargSeq) for automated generation of STR and DIP profiles from buccal swab samples and blood stains. The system fully integrates the processes of DNA extraction, polymerase chain reaction (PCR) amplification, and electrophoresis separation using microfluidic biochip technology. Internal validation studies were performed using RTyper 21 or DIP 38 chip cartridges with single-source reference samples according to the Scientific Working Group for DNA Analysis Methods guidelines. These results indicated that the Quick TargSeq system can process reference samples and generate STR or DIP profiles in approximately 2 h, and the profiles were concordant with those determined using traditional STR or DIP analysis methods. Thus, reproducible and concordant DNA profiles were obtained from reference samples. Throughout the study, no lane-to-lane or run-to-run contamination was observed. The Quick TargSeq system produced full profiles from buccal swabs with at least eight swipes, dried blood spot cards with two 2-mm disks, or 10 ng of purified DNA. Potential PCR inhibitors (i.e., coffee, smoking tobacco, and chewing tobacco) did not appear to affect the amplification reactions of the instrument. The overall success rate and concordance rate of 153 samples were 94.12% and 93.44%, respectively, which is comparable to other commercially available rapid DNA instruments. A blind test initiated by a DNA expert group showed that the system can correctly produce DNA profiles with 97.29% genotype concordance with standard bench-processing methods, and the profiles can be uploaded into the national DNA database. These results demonstrated that the Quick TargSeq system can rapidly generate reliable DNA profiles in an automated manner and has the potential for use in the field and forensic laboratories.

Assessment of the efficiency of DNA isolation and profiling applying a temperature-driven method in human remains

The identification of human remains is of utmost importance in a variety of scenarios. One of the primary identification methods is DNA. DNA extraction from human remains could be difficult, particularly in situations where the remains have been exposed to environmental conditions and other insults. Several studies tried to improve extraction by applying different approaches. ForensicGEM Universal (MicroGem) is a single-tube approach to DNA extraction and a temperature-driven method that could have some advantages with respect to previous techniques, among them, reducing the risk of contamination, not requiring specialized equipment, or several steps to perform. The aim of this study was to assess, for the first time, the efficiency of DNA extraction and quality of STR profiles applying the MicroGem protocol and modifications of this protocol from tooth samples in comparison with automatic extraction (AE). Our results indicated that AE and MicroGem performed similar, though with variability depending on the MicroGem modifications, increasing the DNA yield and STR profile quality when DNA is concentrated with Microcon. These findings demonstrated the efficiency of this methodology for DNA extraction from human remains while also providing a simple and quick technique suitable to apply in a variety of forensic scenarios.

Epigenetic-based age prediction in blood samples: Model development

Aging is a complex process influenced by genetic, epigenetic, and environmental factors that lead to tissue deterioration and frailty. Epigenetic mechanisms, such as DNA methylation, play a significant role in gene expression regulation and aging. This study presents a new age estimation model developed for the Turkish population using blood samples. Eight CpG sites in loci TOM1L1, ELOVL2, ASPA, FHL2, C1orf132, CCDC102B, cg07082267, and RASSF5 were selected based on their correlation with age. Methylation patterns of these sites were analyzed in blood samples from 100 volunteers, grouped into age categories (20-35, 36-55, and ≥56). Sensitivity analysis indicated a reliable performance with DNA inputs ≥1 ng. Statistical modeling, utilizing Multiple Linear Regression, underscores the reliability of the primary 6-CpG model, excluding cg07082267 and TOM1L1. This model demonstrates strong correlations with chronological age (r = 0.941) and explains 88% of the age variance with low error rates (MAE = 4.07, RMSE = 5.73 years). Validation procedures, including a training-test split and fivefold cross-validation, consistently confirm the model's accuracy and consistency. The study indicates minimal variation in error scores across age cohorts and no significant gender differences. The developed model showed strong predictive accuracy, with the ability to estimate age within certain prediction intervals. This study contributes to the age prediction by using DNA methylation patterns, which can have disparate applications, including forensic and clinical assessments.

Developmental validation of a novel 8-dye fluorescent typing system with 59 Y-STRs and 3 Y-indels for forensic applications

An 8-dye fluorescence-labeling forensic Y-chromosomal short tandem repeats (Y-STRs) kit, the 62-plex Y-STR multiplex amplification system, was developed and optimized. The system was validated by testing PCR conditions, stutter ratios (SR) and peak height ratios, sensitivity, mixture samples, precision and accuracy, species-specificity, and inhibition studies according to the Scientific Working Group on DNA Analysis Methods guidelines. PCR-based studies showed that the recommended PCR conditions were optimized for this kit. In the sensitivity study, a full profile was obtained from template DNA with a quantity of u125 pg. Consistent profiles were obtained from three different laboratories. The SRs in all loci were less than 15%, and nice balance and suitable average peak height were shown. No peaks were detected in the profiles of common animal species and microorganisms. In the male-male mixture studies, all loci were observed at a ratio of 1:8, and in the male-female mixture study, all alleles could be profiled at a ratio of 1:500 if the male DNA inputs were ≥0.5 ng/µL. An inhibitor study demonstrated that the kit had varying degrees of resistance to the presence of common inhibitors. Population study demonstrated the 62-plex Y-STR Kit improved the power of discrimination in unrelated Chinese Han males (n = 192). When haplotype diversity was 1, the probability of discrimination power of the 62-plex Y-STR Kit was 0.9948, which is suitable for forensic investigations. The results show that the developed 8-dye fluorescence labeling 62 loci system is sensitive, robust, convenient, and highly informative for forensic applications.

An optimal skeletal element for DNA testing: Evaluation of DNA quantity and quality from various bone types in routine forensic practice

For human identification, the quality and quantity of DNA must be sufficient for amplification and analysis. When DNA extraction from bone tissues and teeth is required, the optimal skeletal elements should be selected as samples for DNA extraction because DNA yield differs among elements. Recently, some studies have reported that a high quantity of high-quality DNA can be extracted from the small cancellous bones of the hands and feet. In this study, we evaluated the effectiveness of small cancellous bones in the human identification of skeletal remains in routine forensic genetic casework. Cancellous bones [phalanges, (meta)carpal bones, and (meta)tarsal bones)] and the cortical bones (femur and petrous bones) and teeth, which have generally been recommended as samples, were collected from the same individuals that needed identifying using DNA analysis in our laboratory. The quantity of DNA from small cancellous bones tended to be higher than that from cortical bones, and the quality from the former was as high as that from the latter. This study showed that in routine forensic casework, the small cancellous bones of the hands and feet should be actively selected as samples for DNA testing.

Development and validation of a new multiplex panel using SNaPshot-based DIP-TriSNP markers for forensic DNA mixtures

Short tandem repeat analysis is challenging when dealing with unbalanced mixtures in forensic cases due to the presence of stutter peaks and large amplicons. In this research, we propose a novel genetic marker called DIP-TriSNP, which combines deletion/insertion polymorphism (DIP) with tri-allelic single nucleotide polymorphism in less than 230 bp length of human genome. Based on multiplex PCR and SNaPShot, a panel, including 14 autosomal DIP-TriSNPs and one Y chromosomal DIP-SNP, had been developed and applied to genotyping 102 unrelated Han Chinese individuals in Sichuan of China and simulated a mixture study. The panel sensitivity can reach as low as 0.1 ng DNA template, and the minor contributor of DNA can be detected with the highest ratio of 19:1, as indicated by the obtained results. In the Sichuan Han population, the cumulative probability of informative genotypes reached 0.997092, with a combined power of discrimination of 0.999999998801. The panel was estimated to detect more than two alleles in at least one locus in 99.69% of mixtures of the Sichuan Han population. In conclusion, DIP-TriSNPs have shown promising as an innovative DNA marker for identifying the minor contributor in unbalanced DNA mixtures, offering advantages such as short amplifications, increased polymorphism, and heightened sensitivity.

A complete pipeline enables haplotyping and phasing macrohaplotype in long sequencing reads for polyploidy samples and a multi-source DNA mixture

Macrohaplotype combines multiple types of phased DNA variants, increasing forensic discrimination power. High-quality long-sequencing reads, for example, PacBio HiFi reads, provide data to detect macrohaplotypes in multiploidy and DNA mixtures. However, the bioinformatics tools for detecting macrohaplotypes are lacking. In this study, we developed a bioinformatics software, MacroHapCaller, in which targeted loci (i.e., short TRs [STRs], single nucleotide polymorphisms, and insertion and deletions) are genotyped and combined with novel algorithms to call macrohaplotypes from long reads. MacroHapCaller uses physical phasing (i.e., read-backed phasing) to identify macrohaplotypes, and thus it can detect multi-allelic macrohaplotypes for a given sample. MacroHapCaller was validated with data generated from our designed targeted PacBio HiFi sequencing pipeline, which sequenced ∼8-kb amplicon regions harboring 20 core forensic STR loci in human benchmark samples HG002 and HG003. MacroHapCaller also was validated in whole-genome long-read sequencing data. Robust and accurate genotyping and phased macrohaplotypes were obtained with MacroHapCaller compared with the known ground truth. MacroHapCaller achieved a higher or consistent genotyping accuracy and faster speed than existing tools HipSTR and DeepVar. MacroHapCaller enables efficient macrohaplotype analysis from high-throughput sequencing data and supports applications using discriminating macrohaplotypes.

Uncovering Forensic Evidence: A Path to Age Estimation through DNA Methylation

Age estimation is a critical aspect of reconstructing a biological profile in forensic sciences. Diverse biochemical processes have been studied in their correlation with age, and the results have driven DNA methylation to the forefront as a promising biomarker. DNA methylation, an epigenetic modification, has been extensively studied in recent years for developing age estimation models in criminalistics and forensic anthropology. Epigenetic clocks, which analyze DNA sites undergoing hypermethylation or hypomethylation as individuals age, have paved the way for improved prediction models. A wide range of biomarkers and methods for DNA methylation analysis have been proposed, achieving different accuracies across samples and cell types. This review extensively explores literature from the past 5 years, showing scientific efforts toward the ultimate goal: applying age prediction models to assist in human identification.

Biogeographical Ancestry Analyses Using the ForenSeqTM DNA Signature Prep Kit and Multiple Prediction Tools

The inference of biogeographical ancestry (BGA) can assist in police investigations of serious crime cases and help to identify missing people and victims of mass disasters. In this study, we evaluated the typing performance of 56 ancestry-informative SNPs in 177 samples using the ForenSeq™ DNA Signature Prep Kit on the MiSeq FGx system. Furthermore, we compared the prediction accuracy of the tools Universal Analysis Software v1.2 (UAS), the FROG-kb, and GenoGeographer when inferring the ancestry of 503 Europeans, 22 non-Europeans, and 5 individuals with co-ancestry. The kit was highly sensitive with complete aiSNP profiles in samples with as low as 250pg input DNA. However, in line with others, we observed low read depth and occasional drop-out in some SNPs. Therefore, we suggest not using less than the recommended 1ng of input DNA. FROG-kb and GenoGeographer accurately predicted both Europeans (99.6% and 91.8% correct, respectively) and non-Europeans (95.4% and 90.9% correct, respectively). The UAS was highly accurate when predicting Europeans (96.0% correct) but performed poorer when predicting non-Europeans (40.9% correct). None of the tools were able to correctly predict individuals with co-ancestry. Our study demonstrates that the use of multiple prediction tools will increase the prediction accuracy of BGA inference in forensic casework.

Joint application of A-InDels and miniSTRs for forensic personal, full and half sibling identifications, and genetic differentiation analyses in two populations from China

BACKGROUND

Previously, a novel multiplex system of 64 loci was constructed based on capillary electrophoresis platform, including 59 autosomal insertion/deletions (A-InDels), two Y-chromosome InDels, two mini short tandem repeats (miniSTRs), and an Amelogenin gene. The aim of this study is to evaluate the efficiencies of this multiplex system for individual identification, paternity testing and biogeographic ancestry inference in Chinese Hezhou Han (CHH) and Hubei Tujia (CTH) groups, providing valuable insights for forensic anthropology and population genetics research.

RESULTS

The cumulative values of power of discrimination (CDP) and probability of exclusion (CPE) for the 59 A-InDels and two miniSTRs were 0.99999999999999999999999999754, 0.99999905; and 0.99999999999999999999999999998, 0.99999898 in CTH and CHH groups, respectively. When the likelihood ratio thresholds were set to 1 or 10, more than 95% of the full sibling pairs could be identified from unrelated individual pairs, and the false positive rates were less than 1.2% in both CTH and CHH groups. Biogeographic ancestry inference models based on 35 populations were constructed with three algorithms: random forest, adaptive boosting and extreme gradient boosting, and then 10-fold cross-validation analyses were applied to test these three models with the average accuracies of 86.59%, 84.22% and 87.80%, respectively. In addition, we also investigated the genetic relationships between the two studied groups with 33 reference populations using population statistical methods of FST, DA, phylogenetic tree, PCA, STRUCTURE and TreeMix analyses. The present results showed that compared to other continental populations, the CTH and CHH groups had closer genetic affinities to East Asian populations.

CONCLUSIONS

This novel multiplex system has high CDP and CPE in CTH and CHH groups, which can be used as a powerful tool for individual identification and paternity testing. According to various genetic analysis methods, the genetic structures of CTH and CHH groups are relatively similar to the reference East Asian populations.

Review: Wildlife forensic genetics-Biological evidence, DNA markers, analytical approaches, and challenges

Wildlife-related crimes are the second most prevalent lawbreaking offense globally. This illicit trade encompasses hunting, breeding and trafficking. Besides diminishing many species and their habitats and ecosystems, hindering the economic development of local communities that depend on them, undermining the rule of law and financing terrorism, various cross-species transmissions (zoonoses) of pathogens, including COVID-19, can be attributed to wildlife crimes. Wildlife forensics applies interdisciplinary scientific analyses to support law enforcement in investigating wildlife crimes. Its main objectives are to identify the taxonomic species in question, determine if a crime has been committed, link a suspect to the crime and support the conviction and prosecution of the perpetrator. This article reviews wildlife crime and its implications, wildlife forensic science investigation, common forms of wildlife biological evidence, including DNA, wildlife DNA techniques and challenges in wildlife forensic genetics. The article also reviews the contributions of genetic markers such as short tandem repeat (STR) and mitochondrial DNA (mtDNA) markers, which provide the probative genetic data representing the bulk of DNA evidence for solving wildlife crime. This review provides an overview of wildlife DNA databases, which are critical for searching and matching forensic DNA profiles and sequences and establishing how frequent forensic DNA profiles and sequences are in a particular population or geographic region. As such, this review will contain an in-depth analysis of the current status of wildlife forensic genetics, and it will be of general interest to wildlife and conservation biologists, law enforcement officers, and academics interested in combating crimes against wildlife using animal forensic DNA methods.

PCR in Forensic Science: A Critical Review

The polymerase chain reaction (PCR) has played a fundamental role in our understanding of the world, and has applications across a broad range of disciplines. The introduction of PCR into forensic science marked the beginning of a new era of DNA profiling. This era has pushed PCR to its limits and allowed genetic data to be generated from trace DNA. Trace samples contain very small amounts of degraded DNA associated with inhibitory compounds and ions. Despite significant development in the PCR process since it was first introduced, the challenges of profiling inhibited and degraded samples remain. This review examines the evolution of the PCR from its inception in the 1980s, through to its current application in forensic science. The driving factors behind PCR evolution for DNA profiling are discussed along with a critical comparison of cycling conditions used in commercial PCR kits. Newer PCR methods that are currently used in forensic practice and beyond are examined, and possible future directions of PCR for DNA profiling are evaluated.

A novel 193-plex MPS panel integrating STRs and SNPs highlights the application value of forensic genetics in individual identification and paternity testing

Massively parallel sequencing (MPS) has emerged as a promising technology for targeting multiple genetic loci simultaneously in forensic genetics. Here, a novel 193-plex panel was designed to target 28 A-STRs, 41 Y-STRs, 21 X-STRs, 3 sex-identified loci, and 100 A-SNPs by employing a single-end 400 bp sequencing strategy on the MGISEQ-2000™ platform. In the present study, a series of validations and sequencing of 1642 population samples were performed to evaluate the overall performance of the MPS-based panel and its practicality in forensic application according to the SWGDAM guidelines. In general, the 193-plex markers in our panel showed good performance in terms of species specificity, stability, and repeatability. Compared to commercial kits, this panel achieved 100% concordance for standard gDNA and 99.87% concordance for 14,560 population genotypes. Moreover, this panel detected 100% of the loci from 0.5 ng of DNA template and all unique alleles at a 1:4 DNA mixture ratio (0.2 ng minor contributor), and the applicability of the proposed approach for tracing and degrading DNA was further supported by case samples. In addition, several forensic parameters of STRs and SNPs were calculated in a population study. High CPE and CPD values greater than 0.9999999 were clearly demonstrated and these results could be useful references for the application of this panel in individual identification and paternity testing. Overall, this 193-plex MPS panel has been shown to be a reliable, repeatable, robust, inexpensive, and powerful tool sufficient for forensic practice.

Developmental validation of a custom-designed Multi-InDel panel: A five-dye multiplex amplification system for challenging DNA samples

The continual investigation of novel genetic markers has yielded promising solutions for addressing the challenges encountered in forensic DNA analysis. In this study, we have introduced a custom-designed panel capable of simultaneously amplifying 41 novel Multi-insertion/deletion (Multi-InDel) markers and an amelogenin locus using the capillary electrophoresis platform. Through a developmental validation study conducted in accordance with guidelines recommended by the Scientific Working Group on DNA Analysis Methods, we demonstrated that the new Multi-InDel system exhibited the sensitivity to produce reliable genotyping profiles with as little as 62.5 pg of template DNA. Accurate and complete genotyping profiles could be obtained even in the presence of specific concentrations of PCR inhibitors. Furthermore, the maximum amplicon size for this system was limited to under 220 bp in the genotyping profile, resulting in its superior efficiency compared to commercially available short tandem repeat kits for both naturally and artificially degraded samples. In the context of mixed DNA analysis, the Multi-InDel system was proved informative in the identification of two-person DNA mixture, even when the template DNA of the minor contributor was as low as 50 pg. In conclusion, a series of performance evaluation studies have provided compelling evidence that the new Multi-InDel system holds promise as a valuable tool for forensic DNA analysis.

Development of a 39 MM-InDel multiplex assay for the forensic application

Insertion/deletion polymorphisms (InDels) are a category of highly prevalent markers in the human genome, characterized by their distinctive attributes, including short amplicon sizes and low mutation rates, which have shown great potential in forensic applications. Multi-allelic InDel and multi-InDel markers, collectively abbreviated as MM-InDels, were developed to enhance polymorphism by the introduction of novel alleles. Nevertheless, the relatively low mutation rates of InDels, coupled with the founder effect, result in distinct allele frequency distributions among populations. The divergent characteristics of InDels in different populations also pose challenges to the establishment of universally efficient InDel multiplex assays. To enhance the system efficiency of the InDel assay and its applicability across diverse populations, 39 MM-InDels with high polymorphism in five different ancestry superpopulations were selected from the 1000 Genomes Project dataset and combined with an amelogenin gender marker to construct a multiplex assay (named MMIDplex). The combined power of discrimination and the cumulative probability of exclusion of 39 MM-InDels were 1 - 1.3 × 10-23 and 1 - 9.83 × 10-6 in the Chinese Han population, and larger than 1-10-19 and 1-10-4 in the reference populations, relatively. These results demonstrate that the MMIDplex assay has the potential to obtain sufficient power for individual identification and paternity test in global populations.

Predicting age from blood by droplet digital PCR using a set of three DNA methylation markers

Evaluation of DNA methylation (DNAm) patterns is a promising tool for age estimation. The duplex droplet digital PCR (ddPCR) method has been recently investigated for DNAm evaluation, revealing to be a potential methodology for DNAm evaluation and molecular age estimation. In this study, we evaluated DNAm levels of CpGs located at the three age-associated genes ELOVL2, FHL2 and PDE4C using ddPCR to develop an age prediction model. Blood-derived DNA samples from 58 healthy individuals (42 women and 16 men; aged 1-93 years old) were submitted to bisulfite conversion followed by ddPCR using dual-labeled probes targeting methylated and unmethylated DNA sequences. Simple linear regression statistics revealed a strong correlation between DNAm levels and chronological age for FHL2 (R = 0.948; P = 1.472 × 10-29) and PDE4C (R = 0.819; P = 3.917 × 10-15), addressing only one CpG for each gene. For the ELOVL2 gene, evaluating five CpG sites in simultaneous, revealed a strong age correlation (R = 0.887; P = 2.099 × 10-20) in a simple linear regression statistics and very strong age correlation (R = 0.926; P = 2.202 × 10-25) when using quadratic regression statistics. The multivariable regression analysis, using methylation information captured on ELOVL2 (squared), FHL2 and PDE4C genes, revealed a very strong age correlation (R = 0.970; P = 5.356 ×10-33), explaining 93.7 % of age variance, displaying a mean absolute deviation (MAD) between chronological and predicted age of 4.657 years (RMSE = 6.044). We postulate that the ddPCR method should be further investigated for DNAm-based age prediction, because it is a relatively simple and an accurate method that can be routinely used in forensic laboratories for testing a few numbers of markers.

Sequence-based population structure, relatedness, and inbreeding estimates for forensic autosomal STR markers

Population data have become available for sequence data to aid forensic investigations and prepare the forensic community in the move towards implementing NGS methods. This comes with a need for updated population genetic parameters estimates to allow DNA evidence evaluations using sequence data. Initial work has been done on a small sample and here we expand this work by providing estimates of population structure and relatedness for autosomal STR data generated by sequencing technologies. We also discuss the effect of inbreeding on forensic calculations and discuss why the use of genotypic-based estimates may be preferred over allelic-based estimates.

An mRNA profiling assay incorporating coding region InDels for body fluid identification and the inference of the donor in mixed samples

Biological traces discovered at crime scenes hold significant significance in forensic investigations. In cases involving mixed body fluid stains, the evidentiary value of DNA profiles depends on the type of body fluid from which the DNA was obtained. Recently, coding region polymorphism analysis has proved to be a promising method for directly linking specific body fluids to their respective DNA contributors in mixtures, which may help to avoid "association fallacy" between separate DNA and RNA evidence. In this study, we present an update on previously reported coding region Single Nucleotide Polymorphisms (cSNPs) by exploring the potential application of coding region Insertion/Deletion polymorphisms (cInDels). Nine promising cInDels, selected from 70 mRNA markers based on stringent screening criteria, were integrated into an existing mRNA profiling assay. Subsequently, the body fluid specificity of our cInDel assay and the genotyping consistency between complementary DNA (cDNA) and genomic DNA (gDNA) were examined. Our study demonstrates that cInDels can function as important multifunctional genetic markers, as they provide not only the ability to confirm the presence of forensically relevant body fluids, but also the ability to associate/dissociate specific body fluids with particular donors.

DNA methylation-based age estimation and quantification of the degradation levels of bisulfite-converted DNA

DNA methylation modifications are known to influence epigenetic phenomena and have been a focus of forensic science research for some time. Degraded DNA after bisulfite treatment is widely used in DNA methylation analysis. In this study, we analyzed methylation levels at 12 CpG sites of four selected genomic regions by pyrosequencing after bisulfite treatment. DNA was extracted from buccal swab samples collected from 102 Japanese individuals who were 21-77 years old. We also developed a simple method to quantify the degradation levels of bisulfite-converted DNA by real-time PCR, and evaluated the effect of DNA degradation on age estimation. We found that the methylation levels and chronological ages were highly correlated in the four selected regions, and the mean absolute deviation (MAD) between chronological and estimated ages was low at 3.88 years. These results indicated that pyrosequencing analysis at the 12 CpGs was useful for age estimation in the Japanese population. To develop a sensitive quantification method, we analyzed the amplification efficiency of short and long fragments from 10 regions by real-time PCR. The amplification efficiency was highest for CCDC102B, and the degradation levels of bisulfite-converted DNA for the 102 samples were categorized as moderately or heavily degraded. For the younger age groups (20-49 years), the MADs were lower for moderately degraded DNA than they were for heavily degraded DNA. This finding indicates that degradation levels affected the accuracy of age estimation in most of the samples; the exception was the samples from the 50-77 years age group.

Development and comparison of forensic interval age prediction models by statistical and machine learning methods based on the methylation rates of ELOVL2 in blood DNA

Age estimation can be useful information for narrowing down candidates of unidentified donors in criminal investigations. Various age estimation models based on DNA methylation biomarkers have been developed for forensic usage in the past decade. However, many of these models using ordinary least squares regression cannot generate an appropriate estimation due to the deterioration in prediction accuracy caused by an increased prediction error in older age groups. In the present study, to address this problem, we developed age estimation models that set an appropriate prediction interval for all age groups by two approaches: a statistical method using quantile regression (QR) and a machine learning method using an artificial neural network (ANN). Methylation datasets (n = 1280, age 0-91 years) of the promoter for the gene encoding ELOVL fatty acid elongase 2 were used to develop the QR and ANN models. By validation using several test datasets, both models were shown to enlarge prediction intervals in accordance with aging and have a high level of correct prediction (>90 %) for older age groups. The QR and ANN models also generated a point age prediction with high accuracy. The ANN model enabled a prediction with a mean absolute error (MAE) of 5.3 years and root mean square error (RMSE) of 7.3 years for the test dataset (n = 549), which were comparable to those of the QR model (MAE = 5.6 years, RMSE = 7.8 years). Their applicability to casework was also confirmed using bloodstain samples stored for various periods of time (1-14 years), indicating the stability of the models for aged bloodstain samples. From these results, it was considered that the proposed models can provide more useful and effective age estimation in forensic settings.

A combined molecular approach utilizing microbial DNA and microRNAs in a qPCR multiplex for the classification of five forensically relevant body fluids

Body fluid identification is an essential step in the forensic biology workflow that can assist DNA analysts in determining where to collect DNA evidence. Current presumptive tests lack the specificity that molecular techniques can achieve; therefore, molecular methods, including microRNA (miRNA) and microbial signature characterization, have been extensively researched in the forensic community. Limitations of each method suggest combining molecular markers to increase the discrimination efficiency of multiple body fluids from a single assay. While microbial signatures have been successful in identifying fluids with high bacterial abundances, microRNAs have shown promise in fluids with low microbial abundance (blood and semen). This project synergized the benefits of microRNAs and microbial DNA to identify multiple body fluids using DNA extracts. A reverse transcription (RT)-qPCR duplex targeting miR-891a and let-7g was validated, and miR-891a differential expression was significantly different between blood and semen. The miRNA duplex was incorporated into a previously reported qPCR multiplex targeting 16S rRNA genes of Lactobacillus crispatus, Bacteroides uniformis, and Streptococcus salivarius to presumptively identify vaginal/menstrual secretions, feces, and saliva, respectively. The combined classification regression tree model resulted in the presumptive classification of five body fluids with 94.6% overall accuracy, now including blood and semen identification. These results provide proof of concept that microRNAs and microbial DNA can classify multiple body fluids simultaneously at the quantification step of the current forensic DNA workflow.

Operationalisation of the ForenSeq® Kintelligence Kit for Australian unidentified and missing persons casework

Single nucleotide polymorphism (SNP) genotyping technologies can generate investigative leads for human remains identification, including estimation of biological sex, biogeographical ancestry (BGA), externally visible characteristics (EVCs), identity, uniparental lineage and extended kinship. The ForenSeq® Kintelligence Kit provides forensic laboratories with the ability to apply this suite of genetic tools to forensic samples using one panel targeting 10,230 SNPs (including 56 ancestry-informative, 24 phenotype-informative, 94 identity-informative, 106 X chromosome, 85 Y chromosome and 9867 kinship-informative SNPs) sequenced on the MiSeq FGx® Sequencing System. The ForenSeq® Kintelligence Kit has been internally validated, optimised and operationalised by the Australian Federal Police National DNA Program for Unidentified and Missing Persons (AFP Program) for coronial casework. The internal validation was conducted according to the Scientific Working Group on DNA Analysis Methods guidelines (excluding mixture analysis), focussing on sample types typically encountered in human remains identification casework, such as bones, teeth, nail, blood and hair. The workflow was optimised for a high throughput library preparation and sequencing workflow, and additional analytical thresholds were developed to improve genotyping accuracy for low DNA input samples. Additionally, the genetic intelligence generated from the kit was compared to the self-declared biological sex, EVCs and BGA of the DNA donors to assess concordance. The kit was able to produce high quality SNP profiles from 1.0 ng down to 0.1 ng of DNA, with high repeatability and reproducibility, and minimal background noise. The prediction accuracy for biological sex (95%), hair colour (58%), eye colour (74%) and BGA inferences (consistent: 74%; partially consistent: 10%; inconclusive: 16%) was determined based on self-declared data. Additionally, SNP profiles from a volunteer family group of ten related individuals were uploaded to GEDmatch PRO™ to assess kinship accuracy. The kit was capable of detecting (97%) and accurately classifying (90%) genetic relationships spanning from first to fifth degree. The Kintelligence Kit provides the AFP Program with a robust and reliable genetic intelligence tool for unidentified and missing persons investigations, which has been designed to sequence multiple challenging samples in a single multiplexed assay using existing laboratory instrumentation.

A proof-of-principle study: The potential application of MiniHap biomarkers in ancestry inference based on the QNome nanopore sequencing

Haplotyped SNPs convey forensic-related information, and microhaplotypes (MHs), as the most representative of this kind of marker, have proved the potential value for human forensics. In recent years, nanopore sequencing technology has developed rapidly, with its outstanding ability to sequence long continuous DNA fragments and obtain phase information, making the detection of longer haplotype marker possible. In this proof-of-principle study, we proposed a new type of forensic marker, MiniHap, based on five or more SNPs within a molecular distance less than 800 bp, and investigated the haplotype data of 56 selected MiniHaps in five Chinese populations using the QNome nanopore sequencing. The sequencing performance, allele (haplotype) frequencies, forensic parameters, effective number of alleles (Ae), and informativeness (In) were subsequently calculated. In addition, we performed principal component analysis (PCA), phylogenetic tree, and structure analysis to investigate the population genetic relationships and ancestry components among the five investigated populations and 26 worldwide populations. MiniHap-04 exhibited remarkable forensic efficacy, with 148 haplotypes reported and the Ae was 66.9268. In addition, the power of discrimination (PD) was 0.9934, the probability of exclusion (PE) was 0.9898, and the In value was 0.7893. Of the 56 loci, 85.71% had PD values above 0.85, 66.07% had PE values above 0.54, 67.86% had Ae values over 7.0%, and 55.36% were with In values above 0.2 across all samples, indicating that most of the MiniHaps are suitable for individual identification, paternity testing, mixture deconvolution, and ancestry inference. Moreover, the results of PCA, phylogenetic tree and structure analysis demonstrated that this MiniHap panel had the competency in continental population ancestry inference, but the differentiation within intracontinental/linguistically restricted subpopulations was not ideal. Such findings suggested that the QNome device for MiniHap detection was feasible and this novel marker has the potential in ancestry inference. Yet, the establishment of a more comprehensive database with sufficient reference population data remains necessary to screen more suitable MiniHaps.

Metatranscriptomic characterization of six types of forensic samples and its potential application to body fluid/tissue identification: A pilot study

Microorganisms are potential markers for identifying body fluids (venous and menstrual blood, semen, saliva, and vaginal secretion) and skin tissue in forensic genetics. Existing published studies have mainly focused on investigating microbial DNA by 16 S rRNA gene sequencing or metagenome shotgun sequencing. We rarely find microbial RNA level investigations on common forensic body fluid/tissue. Therefore, the use of metatranscriptomics to characterize common forensic body fluids/tissue has not been explored in detail, and the potential application of metatranscriptomics in forensic science remains unknown. Here, we performed 30 metatranscriptome analyses on six types of common forensic sample from healthy volunteers by massively parallel sequencing. After quality control and host RNA filtering, a total of 345,300 unigenes were assembled from clean reads. Four kingdoms, 137 phyla, 267 classes, 488 orders, 985 families, 2052 genera, and 4690 species were annotated across all samples. Alpha- and beta-diversity and differential analysis were also performed. As a result, the saliva and skin groups demonstrated high alpha diversity (Simpson index), while the venous blood group exhibited the lowest diversity despite a high Chao1 index. Specifically, we discussed potential microorganism contamination and the "core microbiome," which may be of special interest to forensic researchers. In addition, we implemented and evaluated artificial neural network (ANN), random forest (RF), and support vector machine (SVM) models for forensic body fluid/tissue identification (BFID) using genus- and species-level metatranscriptome profiles. The ANN and RF prediction models discriminated six forensic body fluids/tissue, demonstrating that the microbial RNA-based method could be applied to BFID. Unlike metagenomic research, metatranscriptomic analysis can provide information about active microbial communities; thus, it may have greater potential to become a powerful tool in forensic science for microbial-based individual identification. This study represents the first attempt to explore the application potential of metatranscriptome profiles in forensic science. Our findings help deepen our understanding of the microorganism community structure at the RNA level and are beneficial for other forensic applications of metatranscriptomics.

Targeted sequencing of high-density SNPs provides an enhanced tool for forensic applications and genetic landscape exploration in Chinese Korean ethnic group

BACKGROUND

In this study, we present a NGS-based panel designed for sequencing 1993 SNP loci for forensic DNA investigation. This panel addresses unique challenges encountered in forensic practice and allows for a comprehensive population genetic study of the Chinese Korean ethnic group. To achieve this, we combine our results with datasets from the 1000 Genomes Project and the Human Genome Diversity Panel.

RESULTS

We demonstrate that this panel is a reliable tool for individual identification and parentage testing, even when dealing with degraded DNA samples featuring exceedingly low SNP detection rates. The performance of this panel for complex kinship determinations, such as half-sibling and grandparent-grandchild scenarios, is also validated by various kinship simulations. Population genetic studies indicate that this panel can uncover population substructures on both global and regional scales. Notably, the Han population can be distinguished from the ethnic minorities in the northern and southern regions of East Asia, suggesting its potential for regional ancestry inference. Furthermore, we highlight that the Chinese Korean ethnic group, along with various Han populations from different regional areas and certain northern ethnic minorities (Daur, Tujia, Japanese, Mongolian, Xibo), exhibit a higher degree of genetic affinities when examined from a genomic perspective.

CONCLUSION

This study provides convincing evidence that the NGS-based panel can serve as a reliable tool for various forensic applications. Moreover, it has helped to enhance our knowledge about the genetic landscape of the Chinese Korean ethnic group.

Development and validation of a custom panel including 114 InDels using massively parallel sequencing for forensic application

Insertion/deletion polymorphisms (InDels) have particular characteristics, such as a relatively low mutation rate, small amplicon size, and no stutter artifacts when genotyped via the capillary electrophoresis platform. It would be an important complementary tool for individual identification and certain kinship analyses. At present, massively parallel sequencing (MPS) has shown excellent application value in forensic studies. Therefore, in this study, we developed a custom MPS InDel panel that contains 114 InDels [77 autosomal InDels (A-InDels), 32 X-chromosomal InDels (X-InDels), and 5 Y-chromosomal InDels) based on previous studies. To assess this panel's performance, several validation experiments were performed, including sensitivity, inhibitor, degraded DNA testing, species specificity, concordance, repeatability, case-type samples, and population studies. The results showed that the lowest DNA input was 0.25 ng. All genotypes were obtained in the presence of 80 ng/µL humic acid, 2000 µmol/L calcium, 3000 µmol/L EDTA and indigo. In degraded DNA testing, 90% of loci could be detected for 16-day-old formalin-fixed hearts. In addition, this panel has good species specificity. The values of combined power of discrimination and the combined power of exclusion for 77 A-InDels were 1-3.9951 × 10-32 and 1-4.2956 × 10-7 , respectively. The combined mean exclusion chance for 32 X-InDels was 0.99999 in trios and 0.99904 in duos. The validation results indicate that this newly developed MPS multiplex system is a robust tool for forensic applications.

Investigation of the impact of biogeographic ancestry on DNA methylation based age predictions comparing a Middle East and a Central European population

DNA methylation based age prediction is a new method in the toolbox of forensic genetics. Typically, the method is applied in the course of police investigation e.g. to predict the age of an unknown person that has left a biological trace at a crime scene. The method can also be used to answer other forensic questions, for example to estimate the age of unknown human bodies in the course of the identification process. In the present study, we tested for a potential impact of biogeographic ancestry (BGA) on age predictions using five age dependent methylated CpG sites within the genetic regions of ELOVL2, MIR29B2CHG, FHL2, KLF14 and TRIM59. We collected 102 blood samples each from donors living in Iraq, Middle East (ME) and Germany, Central Europe (EU). Both sample sets were matched in sex and age ranging from 18 to 68 years with exactly one male and female sample per year of age. All samples were analyzed by bisulfite pyrosequencing applying a multiplex pre-amplification strategy based on a single input of 35 ng converted DNA in the PCR. For the CpGs in MIR29B2CHG, FHL2 and KLF14, we observed significantly different methylation levels between the two populations. While we were able to train two highly accurate prediction models for the respective population with mean absolute deviations between predicted and actual ages (MAD) of 3.34 years for the ME model, and 2.72 years for the EU model, we found an absolute prediction difference between the two population specific models of more than 4 years. A combined model for both populations compensated the methylation difference between the two populations, providing MADs of prediction of only 3.81 years for ME and 3.31 years for EU samples. In total, the results of the present study strongly support the benefit of BGA information for more reliable methylation based age predictions.

Partial validation of multiplexed real-time quantitative PCR assays for forensic body fluid identification

Confirmatory body fluid identification using messenger RNA (mRNA) is a well-established technique to address issues encountered with conventional testing - such as poor sensitivity, specificity, and a lack of available tests for all body fluids of interest. For over a decade, endpoint reverse-transcription polymerase chain reaction (RT-PCR) assays have been used in forensic casework for such purposes. However, in comparison with real-time quantitative RT-PCR (RT-qPCR), endpoint RT-PCR has lower sensitivity, precision, and linear dynamic range. This research details the multiplexing and partial validation of confirmatory RT-qPCR assays. We have previously described novel assays for a range of body fluid targets and identified an optimal commercial kit for their amplification. Here, multiplexing was undertaken to form three assays: circulatory blood (SLC4A1) and menstrual fluid (STC1), saliva (HTN3) and vaginal material (CYP2B7P), and spermatozoa (PRM1) and seminal fluid (KLK2), all including a synthetic internal control RNA. Partial validation of the multiplexed assays incorporated the MIQE guidelines, ISO requirements, and SWGDAM guidelines. Using receiver operating characteristic (ROC) curves, each marker was significantly different from an uninformative assay and optimal cut-offs were all above 35 cycles. All assays showed a wide LDR (ranging from 3 to 5 logs with most R2 > 0.99), and high precision (most mean CV < 1 %). STC1 showed some instances of sporadic expression in blood, semen, and vaginal material at high CT values. CYP2B7P showed off-target expression in semen and blood. The sensitivities were approximated as; saliva: 1 in 1,000 dilution of a whole buccal swab, circulatory blood: 0.01-0.1 µL blood, menstrual fluid: 1 in 10,000 dilution of a whole menstrual swab, spermatozoa: 0.001 µL semen, seminal fluid: 0.01 µL semen, and vaginal material: 1 in 1,000 dilution of a whole vaginal swab. A total of 16 mock body fluid extract mixtures and 18 swab mixtures were tested and had 100% and 99% detection of target markers below each specific cut-off, respectively. Some mixtures containing high volumes of blood and semen showed off-target CYP2B7P expression. The successful application of a probabilistic model to the RT-qPCR data was also demonstrated. Further work will involve full developmental validation.

The screening and validation process of miR-223-3p for saliva identification

More accurate identification of the types of body fluids left at a crime scene is indispensable for improving the judicial chain of evidence. MicroRNAs (miRNAs) have become recognized as ideal molecular markers for the identification of body fluids in forensic science due to their short length, stability and high tissue specificity. In this study, small RNA sequencing was performed on 20 samples of five types of body fluids (peripheral blood, menstrual blood, saliva, semen, and vaginal secretions) with the BGISEQ-500 sequencing platform, and the specific miRNA markers of saliva and vaginal secretions were screened by bioinformatics methods, including differential expression analysis and significant enrichment analysis. Through RT-qPCR validation of 169 samples, we confirmed that miR-223-3p can be used as a saliva-specific marker. In addition, we considered miR-223-3p in combination with four other miRNA molecules (miR-451a, miR-891a-5p, miR-144-5p, miR-203a-3p) that had been previously screened and verified in our laboratory, and seven body fluid prediction models based on machine learning algorithms were constructed and verified. The results showed that a kernel density estimation (KDE) model based on the five miRNA markers for body fluid identification could achieve 100% accuracy in the samples tested in the present study.

Forensic validation of a combined analysis of mRNA and miRNA markers for precise tissue origin inferences of five kinds of body fluids by RT-qPCR

Correctly inferring the tissue origin types of forensic-relevant body fluids left at a crime scene is beneficial for reconstructing a crime scene. However, it is still a challenge to accurately identify different kinds of body fluids at a crime scene. Shorter sequence length and anti-degradation microRNA (miRNA) can be used to infer the tissue sources of biological fluid traces, but a limited number of miRNAs are tissue specific. The application of messenger RNA (mRNA) has been confirmed by different studies based on its high tissue specificity. According to the differential expression features of mRNA or miRNA in forensically relevant body fluids, this study developed a simultaneously reversed mRNA and miRNA system and then used these two types of RNAs for the determinations of five common kinds of body fluids. Compared with previously reported single kind of mRNA or miRNA assay, the combined mRNA and miRNA system showed good advantages for human body fluid identifications, especially it could be applied in mixed samples. In conclusion, the obtained results indicated that this combined mRNA and miRNA system might provide a scientific and accurate reference for body fluid identifications.

Development of an epigenetic age predictor for costal cartilage with a simultaneous somatic tissue differentiation system

Age prediction from DNA has been a topic of interest in recent years due to the promising results obtained when using epigenetic markers. Since DNA methylation gradually changes across the individual's lifetime, prediction models have been developed accordingly for age estimation. The tissue-dependence for this biomarker usually necessitates the development of tissue-specific age prediction models, in this way, multiple models for age inference have been constructed for the most commonly encountered forensic tissues (blood, oral mucosa, semen). The analysis of skeletal remains has also been attempted and prediction models for bone have now been reported. Recently, the VISAGE Enhanced Tool was developed for the simultaneous DNA methylation analysis of 8 age-correlated loci using targeted high-throughput sequencing. It has been shown that this method is compatible with epigenetic age estimation models for blood, buccal cells, and bone. Since when dealing with decomposed cadavers or postmortem samples, cartilage samples are also an important biological source, an age prediction model for cartilage has been generated in the present study based on methylation data collected using the VISAGE Enhanced Tool. In this way, we have developed a forensic cartilage age prediction model using a training set composed of 109 samples (19-74 age range) based on DNA methylation levels from three CpGs in FHL2, TRIM59 and KLF14, using multivariate quantile regression which provides a mean absolute error (MAE) of ± 4.41 years. An independent testing set composed of 72 samples (19-75 age range) was also analyzed and provided an MAE of ± 4.26 years. In addition, we demonstrate that the 8 VISAGE markers, comprising EDARADD, TRIM59, ELOVL2, MIR29B2CHG, PDE4C, ASPA, FHL2 and KLF14, can be used as tissue prediction markers which provide reliable blood, buccal cells, bone, and cartilage differentiation using a developed multinomial logistic regression model. A training set composed of 392 samples (n = 87 blood, n = 86 buccal cells, n = 110 bone and n = 109 cartilage) was used for building the model (correct classifications: 98.72%, sensitivity: 0.988, specificity: 0.996) and validation was performed using a testing set composed of 192 samples (n = 38 blood, n = 36 buccal cells, n = 46 bone and n = 72 cartilage) showing similar predictive success to the training set (correct classifications: 97.4%, sensitivity: 0.968, specificity: 0.991). By developing both a new cartilage age model and a tissue differentiation model, our study significantly expands the use of the VISAGE Enhanced Tool while increasing the amount of DNA methylation-based information obtained from a single sample and a single forensic laboratory analysis. Both models have been placed in the open-access Snipper forensic classification website.

TrACES of time: Transcriptomic analyses for the contextualization of evidential stains - Identification of RNA markers for estimating time-of-day of bloodstain deposition

Obtaining forensically relevant information beyond who deposited a biological stain on how and under which circumstances it was deposited is a question of increasing importance in forensic molecular biology. In the past few years, several studies have been produced on the potential of gene expression analysis to deliver relevant contextualizing information, e.g. on nature and condition of a stain as well as aspects of stain deposition timing. However, previous attempts to predict the time-of-day of sample deposition were all based on and thus limited by previously described diurnal oscillators. Herein, we newly approached this goal by applying current sequencing technologies and statistical methods to identify novel candidate markers for forensic time-of-day predictions from whole transcriptome analyses. To this purpose, we collected whole blood samples from ten individuals at eight different time points throughout the day, performed whole transcriptome sequencing and applied biostatistical algorithms to identify 81 mRNA markers with significantly differential expression as candidates to predict the time of day. In addition, we performed qPCR analysis to assess the characteristics of a subset of 13 candidate predictors in dried and aged blood stains. While we demonstrated the general possibility of using the selected candidate markers to predict time-of-day of sample deposition, we also observed notable variation between different donors and storage conditions, highlighting the relevance of employing accurate quantification methods in combination with robust normalization procedures.This study's results are foundational and may be built upon when developing a targeted assay for time-of-day predictions from forensic blood samples in the future.

Application of SNaPshot Technology in Semen-Specific cSNP Genetic Marker

OBJECTIVES

To explore the feasibility of genetic marker detection of semen-specific coding region single nucleotide polymorphism (cSNP) based on SNaPshot technology in semen stains and mixed body fluid identification.

METHODS

Genomic DNA (gDNA) and total RNA were extracted from 16 semen stains and 11 mixtures composed of semen and venous blood, and the total RNA was reverse transcribed into complementary DNA (cDNA). The cSNP genetic markers were screened on the validated semen-specific mRNA coding genes. The cSNP multiplex detection system based on SNaPshot technology was established, and samples were genotyped by capillary electrophoresis (CE).

RESULTS

A multiplex detection system containing 5 semen-specific cSNPs was successfully established. In 16 semen samples, except the cSNP located in the TGM4 gene showed allele loss in cDNA detection results, the gDNA and cDNA typing results of other cSNPs were highly consistent. When detecting semen-venous blood mixtures, the results of cSNP typing detected were consistent with the genotype of semen donor and were not interfered by the genotype of venous blood donor.

CONCLUSIONS

The method of semen-specific cSNPs detection by SNaPshot technology method can be applied to the genotyping of semen (stains) and provide information for determining the origin of semen in mixed body fluids (stains).

Research Progress on Molecular Biology of Human Height Estimation

As an important anthropometric characteristic, human height not only contributes to the recognition of other anthropological characteristics and genetic risk factors, but also is an important part of forensic DNA phenotyping studies. Accurate estimation of height can provide more complete information about the phenotype of suspects and provide help to solve cases. In recent years, having benefited from the rapid development of molecular biological techniques and bioinformatics, height-related genetics research has made some progress. This paper describes the research progress of human height estimation from the genetic variation and the epigenetic inheritance perspectives and looks into the future research direction.

Construction of an Analysis Model of mRNA Markers in Menstrual Blood Based on Naïve Bayes and Multivariate Logistic Regression Methods

OBJECTIVES

To establish the menstrual blood identification model based on Naïve Bayes and multivariate logistic regression methods by using specific mRNA markers in menstrual blood detection technology combined with statistical methods, and to quantitatively distinguish menstrual blood from other body fluids.

METHODS

Body fluids including 86 menstrual blood, 48 peripheral blood, 48 vaginal secretions, 24 semen and 24 saliva samples were collected. RNA of the samples was extracted and cDNA was obtained by reverse transcription. Five menstrual blood-specific markers including members of the matrix metalloproteinase (MMP) family MMP3, MMP7, MMP11, progestogens associated endometrial protein (PAEP) and stanniocalcin-1 (STC1) were amplified and analyzed by electrophoresis. The results were analyzed by Naïve Bayes and multivariate logistic regression.

RESULTS

The accuracy of the classification model constructed was 88.37% by Naïve Bayes and 91.86% by multivariate logistic regression. In non-menstrual blood samples, the distinguishing accuracy of peripheral blood, saliva and semen was generally higher than 90%, while the distinguishing accuracy of vaginal secretions was lower, which were 16.67% and 33.33%, respectively.

CONCLUSIONS

The mRNA detection technology combined with statistical methods can be used to establish a classification and discrimination model for menstrual blood, which can distignuish the menstrual blood and other body fluids, and quantitative description of analysis results, which has a certain application value in body fluid stain identification.

Cross-Platform Application of DNA Methylation Age Estimation Model in Eastern Chinese Han Population

OBJECTIVES

To evaluate the forensic application value of an age estimation model based on DNA methylation in eastern Chinese Han population, and to provide a theoretical basis for exploring age estimation models suitable for different detection platforms.

METHODS

According to the 6 age-related methylation sites in the published blood DNA methylation age estimation models of Chinese Han population, the DNA methylation level of 48 samples was detected by pyrosequencing and next-generation sequencing (NGS). After submitting DNA methylation levels to the age estimation model, the DNA methylation ages were predicted and compared with their real ages.

RESULTS

The 6 DNA methylation sites in both detection techniques were age-related, with an R2 of 0.85 and a median absolute deviation (MAD) of 4.81 years when using pyrosequencing；with an R2 of 0.84 and MAD of 4.41 years when using NGS.

CONCLUSIONS

The blood DNA methylation age estimation model can be used under pyrosequencing and multi-purpose regional methylation enrichment sequencing technology based on NGS and it can accurately estimate the age.