DNA methylation-based organ tissue identification: Marker identification, SNaPshot multiplex assay development, and interlaboratory comparison.

Identifying body fluids and organ tissues is highly significant as they can offer crucial evidence in criminal investigations and aid the court in making informed decisions, primarily through evaluating the biological source and possibly at the activity level up to death or fatal damage. In this study, organ tissue-specific CpG markers were identified from Illumina's methylation EPIC array data of nine organ tissues, including epidermis, dermis, heart, skeletal muscle, blood, kidney, brain, lung, and liver, from autopsies of 10 Koreans. Through the validation test using 43 samples, 18 hypomethylation markers, with two markers for each organ tissue type, were selected to construct a SNaPshot assay. Two multiplex assays involving forward and reverse SBE primers were designed to help investigators accurately determine the organ origin of the analyzed tissue samples through repeated analysis of the same PCR products for markers. The developed multiplex demonstrated high accuracy, achieving 100.0 % correct detection of the presence of nine organ tissue types in 88 samples from autopsies of 10 Asians. However, two lung samples showed additional positive indications of the presence of blood. An interlaboratory comparison using 80 autopsy samples (heart, skeletal muscle, blood, kidney cortex, kidney medulla, brain, lung, and liver) from 10 individuals in Germany revealed overall comparable results with correct detection of the presence of eight organ tissue types in 92.5 % samples (74 of 80 samples). In the case of six samples, it was impossible to determine the correct tissue successfully due to drop-outs of unmethylation signals at target tissue marker loci. One of these lung samples revealed only non-intended off-target signals for blood. The observed differences might be due to differences in sample collection during routine autopsy, technical differences due to the PCR cycler, and the threshold used for signal calling. Indicating the presence of additional tissue type and off-target unmethylation signals seems alleviated by applying more stringent hypomethylation thresholds. Therefore, the developed SNaPshot multiplex assays will be valuable for forensic investigators dealing with organ tissue identification, as well as for prosecutors and defense aiming to establish the circumstances that occurred at the crime scene.

Assessment of mitochondrial DNA copy number variation relative to nuclear DNA quantity between different tissues.

Mitochondrial DNA is a widely tested genetic marker in various fields of research and diagnostics. Nonetheless, there is still little understanding on its abundance and quality within different tissues. Aiming to obtain deeper knowledge about the content and quality of mtDNA, we investigated nine tissues including blood, bone, brain, hair (root and shaft), cardiac muscle, liver, lung, skeletal muscle, and buccal mucosa of 32 deceased individuals using two real-time quantitative PCR-based assays with differently sized mtDNA and nDNA targets. The results revealed that the quantity of nDNA is a weak surrogate to estimate mtDNA quantities among tissues of an individual, as well as tissues across individuals. Especially hair showed extreme variation, depicting a range of multiple magnitudes of mtDNA molecules per hair fragment. Furthermore, degradation can lead to fewer fragments being available for PCR. The results call for parallel determination of the quantity and quality of mtDNA prior to downstream genotyping assays.

Getting the chronological age out of DNA: using insights of age-dependent DNA methylation for forensic DNA applications.

BACKGROUND: DNA analysis for forensic investigations has a long tradition with important developments and optimizations since its first application. Traditionally, short tandem repeats analysis has been the most powerful method for the identification of individuals. However, in addition, epigenetic changes, i.e., DNA methylation, came into focus of forensic DNA research. Chronological age prediction is one promising application to allow for narrowing the pool of possible individuals who caused a trace, as well as to support the identification of unknown bodies and for age verification of living individuals. OBJECTIVE: This review aims to provide an overview of the current knowledge, possibilities, and (current) limitations about DNA methylation-based chronological age prediction with emphasis on forensic application. METHODS: The development, implementation and application of age prediction tools requires a deep understanding about the biological background, the analysis methods, the age-dependent DNA methylation markers, as well as the mathematical models for age prediction and their evaluation. Furthermore, additional influences can have an impact. Therefore, the literature was evaluated in respect to these diverse topics. CONCLUSION: The numerous research efforts in recent years have led to a rapid change in our understanding of the application of DNA methylation for chronological age prediction, which is now on the way to implementation and validation. Knowledge of the various aspects leads to a better understanding and allows a more informed interpretation of DNAm quantification results, as well as the obtained results by the age prediction tools.

A collaborative exercise on DNA methylation-based age prediction and body fluid typing.

DNA methylation has become one of the most useful biomarkers for age prediction and body fluid identification in the forensic field. Therefore, several assays have been developed to detect age-associated and body fluid-specific DNA methylation changes. Among the many methods developed, SNaPshot-based assays should be particularly useful in forensic laboratories, as they permit multiplex analysis and use the same capillary electrophoresis instrumentation as STR analysis. However, technical validation of any developed assays is crucial for their proper integration into routine forensic workflow. In the present collaborative exercise, two SNaPshot multiplex assays for age prediction and a SNaPshot multiplex for body fluid identification were tested in twelve laboratories. The experimental set-up of the exercise was designed to reflect the entire workflow of SNaPshot-based methylation analysis and involved four increasingly complex tasks designed to detect potential factors influencing methylation measurements. The results of body fluid identification from each laboratory provided sufficient information to determine appropriate age prediction methods in subsequent analysis. In age prediction, systematic measurement differences resulting from the type of genetic analyzer used were identified as the biggest cause of DNA methylation variation between laboratories. Also, the use of a buffer that ensures a high ratio of specific to non-specific primer binding resulted in changes in DNA methylation measurement, especially when using degenerate primers in the PCR reaction. In addition, high input volumes of bisulfite-converted DNA often caused PCR failure, presumably due to carry-over of PCR inhibitors from the bisulfite conversion reaction. The proficiency of the analysts and experimental conditions for efficient SNaPshot reactions were also important for consistent DNA methylation measurement. Several bisulfite conversion kits were used for this study, but differences resulting from the use of any specific kit were not clearly discerned. Even when different experimental settings were used in each laboratory, a positive outcome of the study was a mean absolute age prediction error amongst participant's data of only 2.7 years for semen, 5.0 years for blood and 3.8 years for saliva.

Inference of recent admixture using genotype data.

The inference of biogeographic ancestry (BGA) has become a focus of forensic genetics. Misinference of BGA can have profound unwanted consequences for investigations and society. We show that recent admixture can lead to misclassification and erroneous inference of ancestry proportions, using state of the art analysis tools with (i) simulations, (ii) 1000 genomes project data, and (iii) two individuals analyzed using the ForenSeq DNA Signature Prep Kit. Subsequently, we extend existing tools for estimation of individual ancestry (IA) by allowing for different IA in both parents, leading to estimates of parental individual ancestry (PIA), and a statistical test for recent admixture. Estimation of PIA outperforms IA in most scenarios of recent admixture. Furthermore, additional information about parental ancestry can be acquired with PIA that may guide casework.

Evidence for multi-copy Mega-NUMTs in the human genome.

The maternal mode of mitochondrial DNA (mtDNA) inheritance is central to human genetics. Recently, evidence for bi-parental inheritance of mtDNA was claimed for individuals of three pedigrees that suffered mitochondrial disorders. We sequenced mtDNA using both direct Sanger and Massively Parallel Sequencing in several tissues of eleven maternally related and other affiliated healthy individuals of a family pedigree and observed mixed mitotypes in eight individuals. Cells without nuclear DNA, i.e. thrombocytes and hair shafts, only showed the mitotype of haplogroup (hg) V. Skin biopsies were prepared to generate ρ° cells void of mtDNA, sequencing of which resulted in a hg U4c1 mitotype. The position of the Mega-NUMT sequence was determined by fluorescence in situ hybridization and two different quantitative PCR assays were used to determine the number of contributing mtDNA copies. Thus, evidence for the presence of repetitive, full mitogenome Mega-NUMTs matching haplogroup U4c1 in various tissues of eight maternally related individuals was provided. Multi-copy Mega-NUMTs mimic mixtures of mtDNA that cannot be experimentally avoided and thus may appear in diverse fields of mtDNA research and diagnostics. We demonstrate that hair shaft mtDNA sequencing provides a simple but reliable approach to exclude NUMTs as source of misleading results.

Mitochondrial DNA control region variation in a population sample from Thailand.

Mitochondrial DNA (mtDNA) control region sequences from hair samples of 213 individuals from Thailand were analyzed using Sanger sequencing. A total of 170 different haplotypes were identified, of which 146 occurred only once (unique haplotypes). The dataset showed a random match probability of 0.87% and a haplotype diversity of 0.9960. The samples were assigned to 85 different haplogroups with B5a, F1a1a, and M being the most frequent ones. Pairwise FST-values between this and other Southeast and East Asian populations revealed significant but relatively low differences, indicating a close relation. Heteroplasmic positions were observed in 12.2% of hair samples confirming the frequent appearance of heteroplasmic positions in hairs. This dataset will complement existing data as an mtDNA reference for forensic investigations.

Get it off, but keep it: Efficient cleaning of hair shafts with parallel DNA extraction of the surface stain.

The analysis of hair samples is a common task in forensic investigations. Material transferred to the surface of a hair during a crime challenges the analysis as it has to be removed efficiently. However, the removal of the stain can also lead to a loss of information on stain contributors. DNA analysis of the stain itself might thus be helpful for the forensic investigation. The aim of this study was the examination of different methods to remove common biological surface stains completely from human hair shafts without hampering the parallel DNA extraction of the cleaned hair shaft and the isolated surface stain (blood, saliva, vaginal secretion, semen, and skin flocks). Four different methods of cleaning (water, lysis buffer, swabbing, NaClO) were compared to their cleaning efficiency as well as their success of mtDNA analysis of three hair donors and the original five stains on the hair. In order to test the suitability of this procedure for future analysis methods, a selection of samples were also sequenced with MPS. Additionally, nuclear DNA analysis of the stain DNA was performed using a screening STR assay to test the potential success for detection of a STR profile. The most efficient removal of the stain was achieved using NaClO, however compromising further analysis of the stain DNA. The best results for cleaning and parallel stain analysis were obtained using a swab moistened with 0.5 % SDS for surface cleaning. Especially water failed to remove stains efficiently, leading to a high amount of mixed mtDNA in the DNA extracts. MPS showed an increased sensitivity for detection of minute mixtures.

Proof of concept study of age-dependent DNA methylation markers across different tissues by massive parallel sequencing.

The use of DNA methylation (DNAm) for chronological age determination has been widely investigated within the last few years for its application within the field of forensic genetics. The majority of forensic studies are based on blood, saliva, and buccal cell samples, respectively. Although these types of samples represent an extensive amount of traces found at a crime scene or are readily available from individuals, samples from other tissues can be relevant for forensic investigations. Age determination could be important for cases involving unidentifiable bodies and based on remaining soft tissue e.g. brain and muscle, or completely depend on hard tissue such as bone. However, due to the cell type specificity of DNAm, it is not evident whether cell type specific age-dependent CpG positions are also applicable for age determination in other cell types. Within this pilot study, we investigated whether 13 previously selected age-dependent loci based on whole blood analysis including amongst others ELOVL2, TRIM59, F5, and KLF14 also have predictive value in other forensically relevant tissues. Samples of brain, bone, muscle, buccal swabs, and whole blood of 29 deceased individuals (age range 0-87 years) were analyzed for these 13 age-dependent markers using massive parallel sequencing. Seven of these loci did show age-dependency in all five tissues. The change of DNAm during lifetime was different in the set of tissues analyzed, and sometimes other CpG sites within the loci showed a higher age-dependency. This pilot study shows the potential of existing blood DNAm markers for age-determination to analyze other tissues than blood. We identified seven known blood-based DNAm markers for use in muscle, brain, bone, buccal swabs, and blood. Nevertheless, a different reference set for each tissue is needed to adapt for tissue-specific changes of the DNAm over time.

Forensic DNA methylation profiling from minimal traces: How low can we go?

Analysis of human DNA methylation (DNAm) can provide additional investigative leads in crime cases, e.g. the type of tissue or body fluid, the chronological age of an individual, and differentiation between identical twins. In contrast to the genetic profile, the DNAm level is not the same in every cell. At the single cell level, DNAm represents a binary event at a defined CpG site (methylated versus non-methylated). The DNAm level from a DNA extract however represents the average level of methylation of the CpG of interest of all molecules in the forensic sample. The variance of DNAm levels between replicates is often attributed to technological issues, i.e. degradation of DNA due to bisulfite treatment, preferential amplification of DNA, and amplification failure. On the other hand, we show that stochastic variations can lead to gross fluctuation in the analysis of methylation levels in samples with low DNA levels. This stochasticity in DNAm results is relevant since low DNA amounts (1pg - 1ng) is rather the norm than the exception when analyzing forensic DNA samples. This study describes a conceptual analysis of DNAm profiling and its dependence on the amount of input DNA. We took a close look at the variation of DNAm analysis due to DNA input and its consequences for different DNAm-based forensic applications. As can be expected, the 95%-confidence interval of measured DNAm becomes narrower with increasing amounts of DNA. We compared this aspect for two different DNAm-based forensic applications: body fluid identification and chronological age determination. Our study shows that DNA amount should be well considered when using DNAm for forensic applications.

Chronological age prediction based on DNA methylation: Massive parallel sequencing and random forest regression.

The use of DNA methylation (DNAm) to obtain additional information in forensic investigations showed to be a promising and increasing field of interest. Prediction of the chronological age based on age-dependent changes in the DNAm of specific CpG sites within the genome is one such potential application. Here we present an age-prediction tool for whole blood based on massive parallel sequencing (MPS) and a random forest machine learning algorithm. MPS allows accurate DNAm determination of pre-selected markers and neighboring CpG-sites to identify the best age-predictive markers for the age-prediction tool. 15 age-dependent markers of different loci were initially chosen based on publicly available 450K microarray data, and 13 finally selected for the age tool based on MPS (DDO, ELOVL2, F5, GRM2, HOXC4, KLF14, LDB2, MEIS1-AS3, NKIRAS2, RPA2, SAMD10, TRIM59, ZYG11A). Whole blood samples of 208 individuals were used for training of the algorithm and a further 104 individuals were used for model evaluation (age 18-69). In the case of KLF14, LDB2, SAMD10, and GRM2, neighboring CpG sites and not the initial 450K sites were chosen for the final model. Cross-validation of the training set leads to a mean absolute deviation (MAD) of 3.21 years and a root-mean square error (RMSE) of 3.97 years. Evaluation of model performance using the test set showed a comparable result (MAD 3.16 years, RMSE 3.93 years). A reduced model based on only the top 4 markers (ELOVL2, F5, KLF14, and TRIM59) resulted in a RMSE of 4.19 years and MAD of 3.24 years for the test set (cross validation training set: RMSE 4.63 years, MAD 3.64 years). The amplified region was additionally investigated for occurrence of SNPs in case of an aberrant DNAm result, which in some cases can be an indication for a deviation in DNAm. Our approach uncovered well-known DNAm age-dependent markers, as well as additional new age-dependent sites for improvement of the model, and allowed the creation of a reliable and accurate epigenetic tool for age-prediction without restriction to a linear change in DNAm with age.

Detection of G1138A Mutation of the FGFR3 Gene in Tooth Material from a 180-Year-Old Museological Achondroplastic Skeleton.

Throughout the last four centuries, many anatomical museums across the world have collected teratological specimens that became precious objects. These can be regarded as spirits of the past which have captured the morphology of diseases through time. These valuable and irreplaceable specimens can be perfectly used in contemporary dysmorphological or genetic research. Unfortunately, due to the historical nature of these specimens and the regularly used aggressive preservation fluids, DNA degradation is often present. Furthermore, the use of material for DNA extraction is restricted to preserve the appearance of these valuable museological specimens. Thus, the most challenging part in this perspective is to harvest sufficient DNA of good quality for further testing without damaging the specimens. Besides fixated specimens, most teratological collections contain dried skeletal and teeth materials which are an excellent source to extract DNA. We here present a DNA-based method that enables genetic identification of the G1138A mutation of the FGFR3 gene in a 180-year-old achondroplastic skeleton, confirming the previously morphologically determined disease. Nuclear DNA was extracted from a premolar tooth and the mutation was found using Sanger sequencing of a small region of the FGFR3 gene.

Detection of a G>C single nucleotide polymorphism within a repetitive DNA sequence by high-resolution DNA melting.

In standard forensic DNA analysis, single base mutations within short tandem repeats (STR) mostly escape detection. In this study, high-resolution DNA melting (HRM) is compared to minisequencing and Sanger sequencing as to determine the most suitable method for detection of a G to C mutation within a repetitive DNA sequence, the STR system DXS10161. It shows an ATG/ATC polymorphism surrounded by a variable number of (TATC) and (ATCT) motifs. Neutral base changes like G:C to C:G result in very low differences in the melting temperature (T m) of the PCR amplicons. By enhanced resolution of fluorescence vs. temperature in HRM, the technique showed to be suitable for detecting a G to C transversion in this repetitive DNA sequence context. Compared to minisequencing, HRM is more time- and cost-effective. Results were confirmed by Sanger sequencing.

Automated Forensic Animal Family Identification by Nested PCR and Melt Curve Analysis on an Off-the-Shelf Thermocycler Augmented with a Centrifugal Microfluidic Disk Segment.

Nested PCR remains a labor-intensive and error-prone biomolecular analysis. Laboratory workflow automation by precise control of minute liquid volumes in centrifugal microfluidic Lab-on-a-Chip systems holds great potential for such applications. However, the majority of these systems require costly custom-made processing devices. Our idea is to augment a standard laboratory device, here a centrifugal real-time PCR thermocycler, with inbuilt liquid handling capabilities for automation. We have developed a microfluidic disk segment enabling an automated nested real-time PCR assay for identification of common European animal groups adapted to forensic standards. For the first time we utilize a novel combination of fluidic elements, including pre-storage of reagents, to automate the assay at constant rotational frequency of an off-the-shelf thermocycler. It provides a universal duplex pre-amplification of short fragments of the mitochondrial 12S rRNA and cytochrome b genes, animal-group-specific main-amplifications, and melting curve analysis for differentiation. The system was characterized with respect to assay sensitivity, specificity, risk of cross-contamination, and detection of minor components in mixtures. 92.2% of the performed tests were recognized as fluidically failure-free sample handling and used for evaluation. Altogether, augmentation of the standard real-time thermocycler with a self-contained centrifugal microfluidic disk segment resulted in an accelerated and automated analysis reducing hands-on time, and circumventing the risk of contamination associated with regular nested PCR protocols.

Evidence for frequent and tissue-specific sequence heteroplasmy in human mitochondrial DNA.

Mitochondrial point heteroplasmy is a common event observed not only in patients with mitochondrial diseases but also in healthy individuals. We here report a comprehensive investigation of heteroplasmy occurrence in human including the whole mitochondrial control region from nine different tissue types of 100 individuals. Sanger sequencing was used as a standard method and results were supported by cloning, minisequencing, and massively parallel sequencing. Only 12% of all individuals showed no heteroplasmy, whereas 88% showed at least one heteroplasmic position within the investigated tissues. In 66% of individuals up to 8 positions were affected. The highest relative number of heteroplasmies was detected in muscle and liver (79%, 69%), followed by brain, hair, and heart (36.7%-30.2%). Lower percentages were observed in bone, blood, lung, and buccal cells (19.8%-16.2%). Accumulation of position-specific heteroplasmies was found in muscle (positions 64, 72, 73, 189, and 408), liver (position 72) and brain (partial deletion at position 71). Deeper analysis of these specific positions in muscle revealed a non-random appearance and position-specific dependency on age. MtDNA heteroplasmy frequency and its potential functional importance have been underestimated in the past and its occurrence is ubiquitous and dependent at least on age, tissue, and position-specific mutation rates.

High-resolution melting of 12S rRNA and cytochrome b DNA sequences for discrimination of species within distinct European animal families.

The cheap and easy identification of species is necessary within multiple fields of molecular biology. The use of high-resolution melting (HRM) of DNA provides a fast closed-tube method for analysis of the sequence composition of the mitochondrial genes 12S rRNA and cytochrome b. We investigated the potential use of HRM for species identification within eleven different animal groups commonly found in Europe by animal-group-specific DNA amplification followed by DNA melting. Influence factors as DNA amount, additional single base alterations, and the existence of mixed samples were taken into consideration. Visual inspection combined with mathematical evaluation of the curve shapes did resolve nearly all species within an animal group. The assay can therefore not only be used for identification of animal groups and mixture analysis but also for species identification within the respective groups. The use of a universal 12S rRNA system additionally revealed a possible approach for species discrimination, mostly by exclusion. The use of the HRM assay showed to be a reliable, fast, and cheap method for species discrimination within a broad range of different animal species and can be used in a flexible "modular" manner depending on the question to be solved.

Bite through the tent.

The authors report on a young boy who was bitten into his face by an unknown animal while being asleep in a tent. Given the bite marks and the location of the scene, members of the mustelidae and canidae families were the first "suspects." Deoxyribunucleic acid (DNA) recovered from the tent's wall was analyzed with regard to parts of the mitochondrial 12S ribosomal ribunucleic acid (12S rRNA) and cytochrome b (cytb) genes as well as nuclear short tandem repeats (STRs). Since Sanger sequencing revealed a mixed sequence with a strong human component overlying the nonhuman contributor, an animal screening using a duplex real-time polymerase chain reaction (PCR) with an intercalating dye and melt curve analysis was employed. The results were later confirmed by cloning. The applied commercial canine STR kit verified the animal family (canidae) but did not help in discriminating the species due to cross-species amplification. In the presented case, the real-time PCR assay offered the cheapest and fastest method for animal family determination, which then allowed for an appropriate and sample-saving strategy to characterize the causative animal species.

Cocktail approach for in vivo phenotyping of 5 major CYP450 isoenzymes: development of an effective sampling, extraction, and analytical procedure and pilot study with comparative genotyping.

In this study, the authors developed a phenotyping method for CYP1A2, 2C9, 2C19, 2D6, and 3A4 using a cocktail of 100 mg caffeine, 125 mg tolbutamide, 20 mg omeprazole, 30 mg dextromethorphan, and 2 mg midazolam. A simple sampling scheme was established collecting 3 blood samples at 0, 4, and 24 hours followed by solid-phase extraction and liquid chromatography/tandem mass spectrometry analysis. After addition of 8 deuterated internal standards and extraction, the analytes were separated using gradient elution with ammonium acetate and methanol. Data acquisition was performed on a triple quadrupole linear ion trap mass spectrometer in multiple-reaction monitoring mode with positive electrospray ionization. The assay was validated according to international guidelines: limits of quantification (LOQs) were between 0.25 and 1.0 ng/mL for all analytes, except for paraxanthine and caffeine (20 ng/mL). Extraction efficiencies ranged between 77% and 103% and matrix effects between 23% and 95%; precision and accuracy data fulfilled accepted criteria. Calibration curves from LOQ to 1000 ng/mL were established for undiluted and 1:10 diluted plasma (r > 0.998). The method was tested in a pilot study with 14 volunteers. Additional genotyping of the probands generally demonstrated good accordance with the measured phenotyping indices but also disclosed certain contradictory results.