Expression profile analysis of piwi-interacting RNA in forensically relevant biological fluids

Developmental validation of an mRNA kit: A 5-dye multiplex assay designed for body-fluid identification

Identifying the sources of biosamples found at crime scenes is crucial for forensic investigations. Among the markers used for body fluid identification (BFI), mRNA has emerged as a well-studied marker because of its high specificity and remarkable stability. Despite this potential, commercially available mRNA kits specifically designed for BFI are lacking. Therefore, we developed an mRNA kit that includes 21 specific mRNA markers of body fluids, along with three housekeeping genes for BFI, to identify four forensic-relevant fluids (blood, semen, saliva, and vaginal fluids). In this study, we tested 451 single-body-fluid samples, validated the universality of the mRNA kit, and obtained a gene expression profile. We performed the validation studies in triplicates and determined the sensitivity, specificity, stability, precision, and repeatability of the mRNA kit. The sensitivity of the kit was found to be 0.1 ng. Our validation process involved the examination of 59 RNA mixtures, 60 body fluids mixtures, and 20 casework samples, which further established the reliability of the kit. Furthermore, we constructed five classifiers that can handle single-body fluids and mixtures using this kit. The classifiers output possibility values and identify the specific body fluids of interest. Our results showed the reliability and suitability of the BFI kit, and the Random Forest classifier performed the best, with 94% precision. In conclusion, we developed an mRNA kit for BFI which can be a promising tool for forensic practice.

Integrative lncRNA, circRNA, and mRNA analysis reveals expression profiles of six forensic body fluids/tissue

RNAs have attracted much attention in forensic body fluid/tissue identification (BFID) due to their tissue-specific expression characteristics. Among RNAs, long RNAs (e.g., mRNA) have a higher probability of containing more polymorphic sites that can be used to assign the specific donor of the body fluid/tissue. However, few studies have characterized their overall profiles in forensic science. In this study, we sequenced the transcriptomes of 30 samples from venous blood, menstrual blood, semen, saliva, vaginal secretion, and skin tissue, obtaining a comprehensive picture of mRNA, lncRNA, and circRNA profiles. A total of 90,305 mRNAs, 102,906 lncRNAs (including 19,549 novel lncRNAs), and 40,204 circRNAs were detected. RNA type distribution, length distribution, and expression distribution were presented according to their annotation and expression level, and many novel body fluid/tissue-specific RNA markers were identified. Furthermore, the cognate relations among the three RNAs were analyzed according to gene annotations. Finally, SNPs and InDels from RNA transcripts were genotyped, and 21,611 multi-SNP and 4,471 multi-InDel transcriptomic microhaplotypes (tMHs) were identified. These results provide a comprehensive understanding of transcriptome profiles, which could provide new avenues for tracing the origin of the body fluid/tissue and identifying an individual.

RNA quality score evaluation: A preliminary study of RNA integrity number (RIN) and RNA integrity and quality number (RNA IQ)

In the past several years, with the in-depth development of RNA-related research, exploring the application of transcriptome and corresponding RNA biomarkers has become one of the research hotspots in the field of forensic science. High-quality RNA is essential for successful downstream workflows, especially in the steps of screening biomarkers by microarray or RNA sequencing (RNA-seq). Thus, accurately evaluating the quality of RNA samples is a critical step in obtaining meaningful expression data. The RNA integrity number (RIN) generated from the Agilent Bioanalyzer system has been widely used for RNA quality control in the past two decades. Recently, Thermo Fisher Scientific launched a ratiometric fluorescence-based method to quickly check whether an RNA sample has degraded, and the results are presented as RNA integrity and quality number (RNA IQ). Both quality score systems determine RNA quality using a numerical system based on a scale of 1-10, with 1 denoting significantly degraded specimens and 10 representing high-quality, intact RNA samples. In this preliminary study, we evaluated the consistency, reproducibility and linearity of two quality scores in RNA quality determination by analyzing heat- and RNase- artificially degraded samples. Meanwhile, the expression levels of three microRNAs (hsa-let-7 g-5p, hsa-miR-93-5p and hsa-miR-191-5p) in intact and severely degraded RNA samples were estimated by TaqMan-qPCR and droplet digital PCR. Overall, both quality scores showed good repeatability and reproducibility in their respective tests. In the samples subjected to thermal degradation, RIN showed a trend corresponding to heating time, while RNA IQ value showed almost no change on the time gradient. However, in RNase A mediated degradation, RNA IQ value observed better linearity. Furthermore, the expression levels of three microRNAs in the severely degraded samples did not show significant changes compared to the intact RNA samples. RNA degradation is a very complex and highly variable process, which is difficult to comprehensively evaluate through any one index and cannot directly compare these two parameters. Nevertheless, combined with previous research results and the expression levels of three microRNAs in this study, analyzing RNA biomarkers with stable regions or small sizes in challenged samples may be a conservative and reliable approach.

Spitting in the wind?-The challenges of RNA sequencing for biomarker discovery from saliva

Forensic trace contextualization, i.e., assessing information beyond who deposited a biological stain, has become an issue of great and steadily growing importance in forensic genetic casework and research. The human transcriptome encodes a wide variety of information and thus has received increasing interest for the identification of biomarkers for different aspects of forensic trace contextualization over the past years. Massively parallel sequencing of reverse-transcribed RNA ("RNA sequencing") has emerged as the gold standard technology to characterize the transcriptome in its entirety and identify RNA markers showing significant expression differences not only between different forensically relevant body fluids but also within a single body fluid between forensically relevant conditions of interest. Here, we analyze the quality and composition of four RNA sequencing datasets (whole transcriptome as well as miRNA sequencing) from two different research projects (the RNAgE project and the TrACES project), aiming at identifying contextualizing forensic biomarker from the forensically relevant body fluid saliva. We describe and characterize challenges of RNA sequencing of saliva samples arising from the presence of oral bacteria, the heterogeneity of sample composition, and the confounding factor of degradation. Based on these observations, we formulate recommendations that might help to improve RNA biomarker discovery from the challenging but forensically relevant body fluid saliva.

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 a novel forensic age estimation strategy for aged blood samples by combining piRNA and miRNA markers

In forensic investigations, age estimation is vital for determining whether a suspect is under or over the legally defined adult age. With breakthroughs in RNA sequencing technology, small noncoding RNAs have provided new ways to solve problems related to the age estimation of trace or aged samples, owing to their small molecular weight and better stability. In our previous study, we had applied miRNAs for the age estimation of bloodstains; however, further improvement of the existing model is needed. PIWI-interacting RNAs (PiRNAs), which are 24-32 nt noncoding small RNA molecules involved in the PIWI-piRNA pathway, play an important role in the aging process. In this study, we explored the possibility of simultaneously analyzing piRNAs and miRNAs for better age estimation purpose. Through massively parallel sequencing, five age-related piRNAs were identified in blood samples that had been stored for eight years. Further real-time PCR analysis revealed that two piRNAs (piR-000753 and piR-020548) showed relatively higher efficiency in age estimation. Additionally, two age-related miRNAs (miR-324-3p and miR-330-5p) were used to build the estimation model. Among all algorithms tested, gradient boosting showed the lowest mean absolute error (MAE) and root mean square error (RMSE) values (3.171 and 4.403 years, respectively) for the validation dataset (n = 110). The errors of the model were less than 5 years and 10 years for 81.82% and 96.36% of the samples, respectively. The results suggest that the combined use of piRNA and miRNA markers may increase the accuracy of age estimation, and our new model has great potential for application in forensic casework.

The potential of using non-coding RNAs in forensic science applications

With the continuous development and integration of molecular biology and forensic science, non-coding RNAs (ncRNAs), especially ncRNAs with regulatory functions such as microRNA, long non-coding RNA, and circular RNA, have recently been actively explored by forensic scholars. In this study, we review the literature on these ncRNAs in various fields of forensic science, including postmortem interval determination, wound age estimation, forensic age assessment, cause of death analysis, and body fluid identification, aiming to evaluate the current research and provide a perspective for future applications.

A Comprehensive Characterization of Small RNA Profiles by Massively Parallel Sequencing in Six Forensic Body Fluids/Tissue

Body fluids/tissue identification (BFID) is an essential procedure in forensic practice, and RNA profiling has become one of the most important methods. Small non-coding RNAs, being expressed in high copy numbers and resistant to degradation, have great potential in BFID but have not been comprehensively characterized in common forensic stains. In this study, the miRNA, piRNA, snoRNA, and snRNA were sequenced in 30 forensic relevant samples (menstrual blood, saliva, semen, skin, venous blood, and vaginal secretion) using the BGI platform. Based on small RNA profiles, relative specific markers (RSM) and absolute specific markers (ASM) were defined, which can be used to identify a specific body fluid/tissue out of two or six, respectively. A total of 5204 small RNAs were discovered including 1394 miRNAs (including 236 novel miRNA), 3157 piRNAs, 636 snoRNAs, and 17 snRNAs. RSMs for 15 pairwise body fluid/tissue groups were discovered by differential RNA analysis. In addition, 90 ASMs that were specifically expressed in a certain type of body fluid/tissue were screened, among them, snoRNAs were reported first in forensic genetics. In brief, our study deepened the understanding of small RNA profiles in forensic stains and offered potential BFID markers that can be applied in different forensic scenarios.

Changes in the Small Noncoding RNAome During M1 and M2 Macrophage Polarization

Macrophages belong to a special phagocytic subgroup of human leukocytes and are one of the important cells of the human immune system. Small noncoding RNAs are a group of small RNA molecules that can be transcribed without the ability to encode proteins but could play a specific function in cells. SncRNAs mainly include microRNAs (miRNAs) and piwi-interacting RNAs (piRNAs), small nucleolar RNAs (snoRNAs), small nuclear RNAs (snRNAs) and repeat RNAs. We used high-throughput sequencing analysis and qPCR to detect the expression changes of the small noncoding RNAome during macrophage polarization. Our results showed that 84 miRNAs and 47 miRNAs with were downregulated during M1 macrophage polarization and that 11 miRNAs were upregulated and 19 miRNAs were downregulated during M2 macrophage polarization. MiR-novel-3-nature and miR-27b-5p could promote expression of TNF-α which was marker gene of M1 macrophages. The piRNA analysis results showed that 69 piRNAs were upregulated and 61 piRNAs were downregulated during M1 macrophage polarization and that 3 piRNAs were upregulated and 10 piRNAs were downregulated during M2 macrophage polarization. DQ551351 and DQ551308 could promote the mRNA expression of TNF-α and DQ551351overexpression promoted the antitumor activity of M1 macrophages. SnoRNA results showed that 62 snoRNAs were upregulated and 59 snoRNAs were downregulated during M1 macrophage polarization, whereas 6 snoRNAs were upregulated and 10 snoRNAs were downregulated during M2 macrophage polarization. Overexpression of snoRNA ENSMUST00000158683.2 could inhibit expression of TNF-α. For snRNA, we found that 12 snRNAs were upregulated and 15 snRNAs were downregulated during M1 macrophage polarization and that 2 snRNAs were upregulated during M2 macrophage polarization. ENSMUSG00000096786 could promote expression of IL-1 and iNOS and ENSMUSG00000096786 overexpression promoted the antitumor activity of M1 macrophages. Analysis of repeat RNAs showed that 7 repeat RNAs were upregulated and 9 repeat RNAs were downregulated during M1 macrophage polarization and that 2 repeat RNAs were downregulated during M2 macrophage polarization. We first reported the expression changes of piRNA, snoRNA, snRNA and repeat RNA during macrophage polarization, and preliminarily confirmed that piRNA, snoRNA and snRNA can regulate the function of macrophages.

Forensic nanopore sequencing of microhaplotype markers using QitanTech's QNome

In recent years, extraordinary progress has been made in genome sequencing technologies, which has led to a decrease in cost and an increase in the diversity of sequenced genomes. Nanopore sequencing is one of the latest genome sequencing technologies. It aims to sequence longer contiguous pieces of DNA, which are essential for resolving structurally complex regions, and provides a new approach for forensic genetics to detect longer markers in real time. To date, multiple studies have been conducted to sequence forensic markers using MinION from Oxford Nanopore Technologies (ONT), and the results indicate that nanopore sequencing holds promise for forensic applications. Qitan Technology (QitanTech) recently launched its first commercial nanopore genome sequencer, QNome. It could achieve a read length of more than 150 kbp, and could generate approximately 500 Mb of data in 8 h. In this pilot study, we explored and validated this alternative nanopore sequencing device for microhaplotype (MH) profiling using a custom set of 15 MH loci. Seventy single-contributor samples were divided into 7 batches, each of which included 10 samples and control DNA 9947A and was sequenced by QNome. MH genotypes generated from QNome were compared to those from Ion Torrent sequencing (Ion S5XL system) to evaluate the accuracy and stability. Twelve samples randomly selected from the last three batches and Control DNA 9947A were also subjected to ONT MinION sequencing (with R9.4 flow cell) for parallel comparison. Based on MHtyper, a bioinformatics workflow developed for automated MH designation, all MH loci can be genotyped and reliably phased using the QNome data, with an overall accuracy of 99.83% (4 errors among 2310 genotypes). Three occurred near or in the region of homopolymer sequences, and one existed within 50 bp of the start of the sequencing reaction. In the last 15 samples (12 individual samples and 3 replicates of control DNA 9947A), two SNPs located at 4-mer homopolymers failed to obtain reliable genotypes on the MinION data. This study shows the potential of state-of-the-art nanopore sequencing methods to analyze forensic MH markers. Given the rapid pace of change, sporadic and nonrepetitive errors presented in this study are expected to be resolved by further developments of nanopore technologies and analysis tools.

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.

Feasibility of using probabilistic methods to analyse microRNA quantitative data in forensically relevant body fluids: a proof-of-principle study

Several studies have confirmed that microRNAs (miRNAs) are promising markers for body fluid identification since they were introduced to this field. However, there is no consensus on the choice of reference genes and identification strategies. In this study, 13 potential candidate miRNAs were screened from three forensically relevant body fluid datasets, and the expression of 12 markers in five body fluids was determined using a real-time quantitative method. Two probabilistic approaches, Naive Bayes (NB) and partial least squares discriminant analysis (PLS-DA), were then applied to predict the origin of the samples to determine whether probabilistic methods are helpful in body fluid identification using miRNA quantitative data. Furthermore, 14 reference combinations were used to validate the influence of different reference choices on the predicted results simultaneously. Our results showed that in the NB model, leave-one-out cross-validation (LOOCV) achieved 100% accuracy and the prediction accuracy of the test set was 100% in most reference combinations. In the PLS-DA model, the first two components could interpret about 80% expression variance and LOOCV achieved 100% accuracy when miR-92a-3p was used as the reference. This study preliminarily proved that probabilistic approaches hold huge potential in miRNA-based body fluid identification, and the choice of references influences the prediction results to a certain extent.

Forensic transcriptome analysis using massively parallel sequencing

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

mRNA profiling of mock casework samples: Results of a FoRNAP collaborative exercise

In recent years, forensic mRNA profiling has increasingly been used to identify the origin of human body fluids. By now, several laboratories have implemented mRNA profiling and also use it in criminal casework. In 2018 the FoRNAP (Forensic RNA Profiling) group was established among a number of these laboratories with the aim of sharing experiences, discussing optimization potential, identifying challenges and suggesting solutions with regards to mRNA profiling and casework. To compare mRNA profiling methods and results a collaborative exercise was organized within the FoRNAP group. Seven laboratories from four countries received 16 stains, comprising six pure body fluid / tissue stains and ten mock casework samples. The laboratories were asked to analyze the provided stains with their in-house method (PCR/CE or MPS) and markers of choice. Five laboratories used a DNA/RNA co-extraction strategy. Overall, up to 11 mRNA markers per body fluid were analyzed. We found that mRNA profiling using different extraction and analysis methods as well as different multiplexes can be applied to casework-like samples. In general, high input samples were typed with high accuracy by all laboratories, regardless of the method used. Irrespective of the analysis strategy, samples of low input or mixed stains were more challenging to analyze and interpret since, alike to DNA profiling, a higher number of markers dropped out and/or additional unexpected markers not consistent with the cell type in question were detected. It could be shown that a plethora of different but valid analysis and interpretation strategies exist and are successfully applied in the Forensic Genetics community. Nevertheless, efforts aiming at optimizing and harmonizing interpretation approaches in order to achieve a higher consistency between laboratories might be desirable in the future. The simultaneous extraction of DNA alongside RNA showed to be an effective approach to identify not only the body fluid present but also to identify the donor(s) of the stain. This allows investigators to gain valuable information about the origin of crime scene samples and the course of events in a crime case.

Expression profile analysis and stability evaluation of 18 small RNAs in the Chinese Han population

In recent years, differentially expressed small RNAs have been widely used to identify the compositions of forensically relevant biological samples, and a vast number of such RNA candidates have been proposed. Nevertheless, when assessing the expression levels of target small RNAs using relative quantitative analysis methods, credible internal controls are usually required for reliable data normalization. Therefore, the identification of optimal reference genes is an important task. In this study, the expression profile of 18 small RNA reference genes was characterized in the Chinese Han population using TaqMan real-time quantitative PCR. Systematic evaluations of these candidate genes were performed based on their expression levels and stability in several common types of body fluids (i.e., venous blood, menstrual blood, saliva, semen, and vaginal secretions). Analysis results from the ΔCq method, BestKeeper, NormFinder, and geNorm were integrated by RefFinder for ranking and comparing the candidates in each type of body fluid. Among all the candidates, miR-191 was identified as the most suitable reference gene because it had a favorable ranking value in all tested samples. In addition, miR-423, miR-93, miR-484, and let-7i were also shown to be applicable reference genes. Overall, this study provides detailed assessment results of these candidate genes in different body fluids; thus, it could be used as a guide for the selection of reference genes according to their performance in the sample of choice.