Search

UNLABELLED: Human age estimation from trace samples may give important leads early in a police investigation by contributing to the description of the perpetrator. Several molecular biomarkers are available for the estimation of chronological age, and currently, DNA methylation patterns are the most promising. In this study, a QIAGEN age protocol for age estimation was tested by five forensic genetic laboratories. The assay comprised bisulfite treatment of the extracted DNA, amplification of five CpG loci (in the genes of ELOVL2, C1orf132, TRIM59, KLF14, and FHL2), and sequencing of the amplicons using the PyroMark Q48 platform. Blood samples from 49 individuals with ages ranging from 18 to 64 years as well as negative and methylation controls were analyzed. An existing age estimation model was applied to display a mean absolute deviation of 3.62 years within the reference data set.KEY POINTS: Age determination as an intelligence tool during investigations can be a powerful tool in forensic genetics.In this study, five laboratories ran 49 samples and obtained a mean absolute deviation of 3.62 years.Five markers were analyzed on a PyroMark Q48 platform.

DNA profiling is central in the field of forensic genetics and allows for the forensic identification of individuals. Biological trace samples consisting of very few cells, however, are a challenge. In Denmark, biological trace samples containing less than 100 pg DNA are termed low template DNA (ltDNA) samples and are not subjected to DNA profiling. In the case of such ltDNA samples, RNA may be used instead of DNA. Each cell only contains two copies of the autosomal DNA, while abundant RNA species may be present and expressed in high copy numbers. As such, much of the genomic information, including sequence variation, should be replicated in multiple copies of messenger RNA (mRNA). Single nucleotide polymorphisms (SNPs) may be used for human identification, and the coding region SNPs in RNA represent a direct link between the DNA and RNA. This thesis describes the development of a prototype mRNA SNP profile for the forensic identification of individuals. Di-allelic SNP markers with high allele frequency were identified in blood stains left to degrade over a two-week time period. Promising markers were retained and further screened in a small subset of the Danish population (41 individuals). Markers in linkage disequilibrium or with skewed allele balances were excluded. The prototype profile of 39 SNPs gives a global mean match probability of 1.88 · 10−12 as well as a mean match probability of 1.86 · 10−10 in the European (non-Finnish) population. Further, the mRNA SNP profile also has the potential for simultaneous identification of the body fluid of a stain and its donor due to the presence of blood-specific mRNA markers. Furthermore, this thesis examined the potential for using RNA markers identified through RNAseq to determine the time since deposition (TSD) of a stain. An accurate method to esti- mate the TSD of forensic trace samples would help investigators evaluate the true re

Estimating an individual’s age can be relevant in several areas primarily related to the clinical and forensic fields. In the latter, estimation of an individual’s chronological age from biological material left by the perpetrator at a crime scene may provide helpful information for police investigation. Estimation of age is also beneficial in immigration cases, where age can affect the person’s protection status under the law, or in disaster victim identification to narrow the list of potential missing persons. In the last decade, research has focused on establishing new approaches for age prediction in the forensic field. From the first forensic age estimations based on morphological inspections of macroscopic changes in bone and teeth, the focus has shifted to molecular methods for age estimation. These methods allow the use of samples from human biological material that does not contain morphological age features and can, in theory, be investigated in traces containing only small amounts of biological material. Molecular methods involving DNA analyses are the primary choice and estimation of DNA methylation levels at specific sites in the genome is the most promising tool. This review aims to provide an overview of the status of forensic age prediction using molecular methods, with particular focus in DNA methylation. The frequent challenges that impact forensic age prediction model development will be addressed, together with the importance of validation efforts within the forensic community.

This report was prepared by Marie-Louise Kampmann, MSc, PhD and Niels Morling, MD, DMSc, Department of Forensic Medicine, University of Copenhagen, Denmark, commissioned by Magistrate Don Mario Carroza Espinosa and Magistrate Ms. Paola Plaza to investigate the death of Pablo Neruda. The objective was to secure the identity of the human remains examined and determine the state of human and bacterial DNA preservation and the presence of contamination and possible pathogens in the remains of Pablo Neruda. This was accomplished with the use of advanced forensic genetic examinations, DNA extraction, sequencing technology, and in-depth bioinformatic analysis.

The performance of Verogen’s ForenSeq mtDNA Whole Genome Kit and its possible application in a forensic setting was evaluated. We carried out three sequencing runs that included 1) a dilution series of two different samples (200–6.25 pg input nuclear DNA), 2) casework samples, and 3) hair, blood, and buccal samples from the same individual. For samples described in 1), the method was tested using different reagent volume in the PCR and library building (recommended volume or half-volume). Reference samples were analysed in duplicates, and positive and negative controls were included in every run. All mtDNA profiles were known from previous genotyping with the Precision ID mtDNA Whole Genome Panel or using long range PCR and the MiSeq instrument. The average read depth per sample was > 1200x in all samples, although some regions had read depth 50x. Using half volume of reagents did not compromise the quality of the results and can be a cost-effective alternative for reference samples. For an analytical threshold of 10%, complete concordance was observed between the duplicates and previous results. Heteroplasmy was consistent with prior observations although frequencies varied up to 10% compared to previous results.

In this pilot study, we made dilutions of DNA/RNA co-extractions (1000–12.5 pg gDNA) and whole transcriptome shotgun sequenced RNA on a NovaSeq 6000 platform (Illumina). Whole blood samples from two individuals were sequenced in triplicate for each dilution. More than 450,000 RNA variants were identified to have overlapping positions with the gnomAD 3.0 DNA variant database in the 50 genes with most reads. Between 7000 and 45,000 of these loci were covered with at least 70 reads in the various dilutions. A total of 171 SNPs with global minor allele frequencies > 0.1 were identified in all dilutions with RNA input equivalent to 12.5 pg gDNA. Of these, 56 were genotyped with the Infinium Omni5–4 SNP-chip (Illumina) using DNA from the two individuals and 53 SNP genotypes in 22 autosomal genes were confirmed between RNA and DNA. The mean match probability based on the most polymorphic SNP in the 22 genes was 5 × 10−8. Thus, shotgun RNA sequencing of low template samples may provide both human identification and tissue identification.

Sodium hypochlorite is an efficient reagent for removal of unwanted DNA from laboratory surfaces. Here, we tested two different chlorine wipes and compared their performance to a 0.9–1.8% hypochlorite solution. WipeClean Chlorine Disinfection wipes contain > 0.1 g sodium hypochlorite/kg, whereas WetWipe Chlorine Desinfection wipes contain > 1000 ppm active chlorine. Clean surfaces were contaminated with 10 µL 0.5 ng/µL of massively parallel sequencing libraries. The DNA was dried and left for 45 min before any treatment. The surfaces were cleaned using either 1) a 0.9–1.8% hypochlorite solution and clean wipes, 2) a WipeClean wipe, 3) a WetWipe, or 4) the surface was not cleaned. All experiments were repeated three times. Subsequently, the surfaces were swabbed using cotton swabs. DNA was extracted from the swabs and the DNA concentrations were determined in quadruplicates by real-time PCR. This protocol was repeated after the soft plastic wrapping around the wipes were left open or closed for several weeks. The results showed that the WipeClean wipes efficiently removed DNA for up to four weeks after the box with the wipes were opened, whereas the WetWipe wipes dried faster and gradually lost their cleaning effect.

Background: No biomarkers have been identified that can classify subtypes of hand eczema (HE). Although skin biopsies represent the gold standard for investigations of the skin, the invasive technique is not favorable when investigating skin from sensitive areas. Recent advances in the use of skin-tape strips for molecular investigations enable noninvasive investigations of HE. Objective: By using whole transcriptome sequencing (WTS), the molecular profile of HE according to different localizations on the hands, etiologies, and clinical/morphological subtypes was investigated. Methods: Thirty adult, Danish HE patients, 12 with and 18 without concurrent atopic dermatitis (AD), as well as 16 controls were included. Tape strip samples were collected from lesional, nonlesional, and healthy skin. Total RNA was extracted and WTS was performed. Results: The largest molecular difference of HE patients with and without AD was found in nonlesional skin areas and included a downregulation of CXCL8 for HE patients without AD. Differences between allergic and irritant contact dermatitis included promising epidermal biomarkers such as EPHA1. Conclusion: Skin tape strip samples could be used to assess the gene expression profile of HE on different localizations of the hands. The skin tape strip method identified new molecular markers that showed promising result for the identification of HE subtypes.

Cerebral palsy (CP) is the most common cause of movement disorders in children. Next generation sequencing (NGS) studies have previously shown that expression levels are fundamentally different in children with CP compared to typically developing (TD). However, given that children are in full development, we might expect gene expression levels to change once maturity is reached. Therefore, the main purpose of this study was to investigate gene expression levels of 93 target genes in adults with CP using NGS on muscle biopsies of the gastrocnemius, taken from 22 participants (n = 12 adults with CP; n = 10 TD adults). Subsequently, we carried out NGS of the mitochondrial genome to identify mtDNA variants, and additionally we studied the mitochondrial content using transmission electron microscopy images of the gastrocnemius muscle. Finally, we compared systemic ion levels between TD adults and adults with CP. Differential gene expression levels were found in genes involved in muscle contraction (MYH1 and MYBPC2), mitochondrial function kATP5J, CYCS and NDUFB6), calcium handling (CAMK2B and ATP2A), metabolism (LPL), muscle signaling (MYC, CREB1, ACVR2B, LMNA and TRIM54), and ECM (TNC). There was no statistical significant difference between CP and TD for mtDNA variant frequencies and mitochondrial content. The ion levels of Ca2+, Na+ and K+ were statistically significantly reduced while the Cl− levels were significant increased in adults with CP compared to TD adults. These results highlight that most transcriptional differences are related to muscle function in adults with CP and that mitochondrial function might be altered but not mitochondrial content.

Environmental DNA (eDNA) and invertebrate-derived DNA (iDNA) are used to survey biodiversity non-invasively to mitigate difficulties in obtaining wildlife samples, particularly in remote areas or for rare species. Recently, eDNA/iDNA were used to monitor known wildlife pathogens; however, most wildlife pathogens are unknown and often evolutionarily divergent. To detect and identify known and novel mammalian viruses from eDNA/iDNA, we used a curated set of RNA oligonucleotides as viral baits in a hybridization capture system coupled with high-throughput sequencing. We detected multiple known and novel mammalian RNA and DNA viruses from multiple viral families from both waterhole eDNA and leech-derived iDNA. Congruence was found between detected hosts and viruses identified in leeches and waterholes. Our results demonstrate that eDNA/iDNA samples represent an effective non-invasive resource for studying wildlife viral diversity and for detecting novel potentially zoonotic viruses prior to their emergence.

Sudden cardiac death (SCD) is a diagnostic challenge in forensic medicine. In a relatively large proportion of the SCDs, the deaths remain unexplained after autopsy. This challenge is likely caused by unknown disease mechanisms. Changes in DNA methylation have been associated with several heart diseases, but the role of DNA methylation in SCD is unknown. In this study, we investigated DNA methylation in two SCD subtypes, sudden unexplained death (SUD) and sudden unexpected death in epilepsy (SUDEP). We assessed DNA methylation of more than 850,000 positions in cardiac tissue from nine SUD and 14 SUDEP cases using the Illumina Infinium MethylationEPIC BeadChip. In total, six differently methylated regions (DMRs) between the SUD and SUDEP cases were identified. The DMRs were located in proximity to or overlapping genes encoding proteins that are a part of the glutathione S-transferase (GST) superfamily. Whole genome sequencing (WGS) showed that the DNA methylation alterations were not caused by genetic changes, while whole transcriptome sequencing (WTS) showed that DNA methylation was associated with expression levels of the GSTT1 gene. In conclusion, our results indicate that cardiac DNA methylation is similar in SUD and SUDEP, but with regional differential methylation in proximity to GST genes.

Cerebral palsy (CP) is a non-progressive motor disorder that affects posture and gait due to contracture development. The purpose of this study is to analyze a possible relation between muscle stiffness and gene expression levels in muscle tissue of children with CP. Next-generation sequencing (NGS) of gene transcripts was carried out in muscle biopsies from gastrocnemius muscle (n = 13 children with CP and n = 13 typical developed (TD) children). Passive stiffness of the ankle plantarflexors was measured. Structural changes of the basement membranes and the sarcomere length were measured. Twelve pre-defined gene target sub-categories of muscle function, structure and metabolism showed significant differences between muscle tissue of CP and TD children. Passive stiffness was significantly correlated to gene expression levels of HSPG2 (p = 0.02; R2 = 0.67), PRELP (p = 0.002; R2 = 0.84), RYR3 (p = 0.04; R2 = 0.66), C COL5A3 (p = 0.0007; R2 = 0.88), ASPH (p = 0.002; R2 = 0.82) and COL4A6 (p = 0.03; R2 = 0.97). Morphological differences in the basement membrane were observed between children with CP and TD children. The sarcomere length was significantly increased in children with CP when compared with TD (p = 0.04). These findings show that gene targets in the categories: calcium handling, basement membrane and collagens, were significantly correlated to passive muscle stiffness. A Reactome pathway analysis showed that pathways involved in DNA repair, ECM proteoglycans and ion homeostasis were amongst the most upregulated pathways in CP, while pathways involved in collagen fibril crosslinking, collagen fibril assembly and collagen turnover were amongst the most downregulated pathways when compared with TD children. These results underline that contracture formation and motor impairment in CP is an interplay between multiple factors.

BACKGROUND: Skin biopsies represent a gold standard in skin immunology and pathology but can cause pain and induce scarring. Non-invasive techniques will facilitate study recruitment of e.g. patients with paediatric atopic dermatitis (AD), hand eczema or facial dermatitis.OBJECTIVE: By RNA sequencing, we examined whether the stratum corneum transcriptome in AD skin can be assessed by tape stripping, as compared to the epidermal transcriptome of AD in skin biopsies. To make the procedure clinically relevant tape strips were stored and shipped at room temperature for up to 3 days.METHODS: Nine adult Caucasian AD patients and three healthy volunteers were included. Tape samples were collected from non-lesional and lesional skin. Biopsies were collected from lesional skin and were split into epidermis and dermis. Total RNA was extracted, and shotgun sequencing was performed.RESULTS: Shotgun sequencing could be performed on skin cells obtained from two consecutive tape strips which had been stored and shipped at room temperature for up to three days. The most prominent differences between the tape strip and biopsy derived transcriptome were due to structural genes, while established molecular markers of AD, including CCL17, CCL22, IL17A and S100A7-S100A9, were also identified in tape strip samples. Furthermore, the tape strip derived transcriptome showed promise in also analysing the skin microbiome.CONCLUSION: Our study shows that the stratum corneum (SC) transcriptome of AD can be assessed by tape stripping the skin, supporting that this method may be central in future skin biomarker research. NCBI GEO data accession: GSE160501.

Massively parallel sequencing (MPS) has revolutionised clinical genetics and research within human genetics by enabling the detection of variants in multiple genes in several samples at the same time. Today, multiple approaches for MPS of DNA are available, including targeted gene sequencing (TGS) panels, whole exome sequencing (WES), and whole genome sequencing (WGS). As MPS is becoming an integrated part of the work in genetic laboratories, it is important to investigate the variant detection performance of the various MPS methods. We compared the results of single nucleotide variant (SNV) detection of three MPS methods: WGS, WES, and HaloPlex target enrichment sequencing (HES) using matched DNA of 10 individuals. The detection performance was investigated in 100 genes associated with cardiomyopathies and channelopathies. The results showed that WGS overall performed better than those of WES and HES. WGS had a more uniform and widespread coverage of the investigated regions compared to WES and HES, which both had a right-skewed coverage distribution and difficulties in covering regions and genes with high GC-content. WGS and WES showed roughly the same high sensitivities for detection of SNVs, whereas HES showed a lower sensitivity due to a higher number of false negative results.

BACKGROUND: Sudden cardiac death (SCD) is a major public health problem and constitutes a diagnostic and preventive challenge in forensic pathology, especially for cases with structural normal hearts at autopsy, so-called sudden arrhythmic death syndrome (SADS). The identification of new genetic risk factors that predispose to SADS is important, because they may contribute to establish the diagnosis and increase the understanding of disease pathways underlying SADS. Pathogenic mutations in the protein coding regions of cardiac genes were found in relation to SADS. However, much remains unknown about variants in non-coding regions of the genome.METHODS AND RESULTS: In this study, we explored the potential of whole genome sequencing (WGS) and whole transcriptome sequencing (WTS) to find DNA variants in SCD victims with structural normal hearts. With focus on the non-coding regulatory regions, we re-examined a cohort of 13 SADS and sudden unexplained death in infancy (SUDI) victims without disease causing DNA variants in recognized cardiac genes. The genetic re-examination of DNA was carried out using frozen tissue samples and WTS was carried out using five distinct formalin fixed and paraffin embedded (FFPE) cardiac tissue samples from each individual, including anterior and posterior walls of the left ventricle, ventricular papillary muscle, septum, and the right ventricle. We identified 23 candidate variants in regulatory sequences of cardiac genes, including a variant in the promotor region of NEXN, c.-194A>G, that was found to be statistically significantly (p < 0.05) associated with decreased expression of NEXN and cardiac hypertrophy.CONCLUSION: With the use of post-mortem FFPE tissues, we highlight the potential of using WTS investigations and compare gene expression levels with DNA variation in regulatory non-coding regions of the genome for a better understanding of the genetics of cardiac diseases leading to SCD.

Efficient methods for building genomic sequencing libraries from degraded DNA have been in place for Illumina sequencing platforms for some years now, but such methods are still lacking for other sequencing platforms. Here, we provide a protocol for building genomic libraries from degraded DNA (archival or ancient sample material) for sequencing on the Ion Torrent™ high-throughput sequencing platforms. In addition to a reduction in time and cost in comparison to commercial kits, this protocol removes purification steps prior to library amplification, an important consideration for work involving historical samples. Libraries prepared using this method are appropriate for either shotgun sequencing or enrichment-based downstream approaches.

Prenatal paternity testing often relies on invasive procedures that cause risk to both the mother and the foetus. Non-invasive, prenatal paternity testing by investigating paternally inherited single nucleotide polymorphisms (SNPs) in cell-free foetal DNA (cffDNA) in maternal plasma was performed at consecutive time points during early gestation. Plasma from 15 pregnant women was investigated at consecutive time points from gestational weeks (GWs) 4–20. The Precision ID Identity Panel and an Ion S5 Sequencer was used to analyse the cffDNA. Paternally inherited foetal SNP alleles were detected from GW7. The median foetal fractions were 0%, 3.9%, 5.1%, 5.2%, and 4.7% at GWs 4, 7, 12, 16, and 20, respectively. The corresponding median numbers of detected paternally inherited foetal autosomal SNP alleles were 0, 3, 9, 10, and 12, respectively. The typical (i.e. geometric mean) paternity indices at GW12 and GW20 were 24 (range 0.0035–8389) and 199 (range 5.1–30,137), respectively. The method is very promising. However, the method can be improved by shortening the lengths of the PCR amplicons and increasing the number of SNPs. To our knowledge, this is the first study to successfully identify paternally inherited foetal SNP alleles at consecutive time points in early gestation independently of the foetal gender.

Quantification of massively parallel sequencing libraries is important for acquisition of monoclonal beads or clusters prior to clonal amplification and to avoid large variations in library coverage when multiple samples are included in one sequencing analysis. No gold standard for quantification of libraries exists. We assessed eight methods of quantification of libraries by quantifying 54 amplicon, six capture, and six shotgun fragment libraries. Chemically synthesized double-stranded DNA was also quantified. Light spectrophotometry, i.e. NanoDrop, was found to give the highest concentration estimates followed by Qubit and electrophoresis-based instruments (Bioanalyzer, TapeStation, GX Touch, and Fragment Analyzer), while SYBR Green and TaqMan based qPCR assays gave the lowest estimates. qPCR gave more accurate predictions of sequencing coverage than Qubit and TapeStation did. Costs, time-consumption, workflow simplicity, and ability to quantify multiple samples are discussed. Technical specifications, advantages, and disadvantages of the various methods are pointed out.

DNA methylation patterns may be used for identification of body fluids and for age estimation of human individuals. We evaluated some of the challenges and pitfalls of studying methylated CpG sites. We compared the methylated CpG analysis of two different methods 1) massively parallel sequencing (MPS) using the Illumina. MiSeq and 2) the iPLEX assay on the MassARRAY® System. On the Illumina MiSeq, the standard deviation of the fraction of methylation was under 3% between replicates, whereas the reproducibility of the MassARRAY® was very difficult to achieve. We tested the denaturation and conversion times of the Premium Bisulfite kit (Diagenode) in order to optimise the conversion rate and minimise DNA degradation. Finally, we tested the reproducibility of the methylation patterns using the Illumina MiSeq., DNA methylation patterns may be used for identification of body fluids and for age estimation of human individuals. We evaluated some of the challenges and pitfalls of studying methylated CpG sites. We compared the methylated CpG analysis of two different methods 1) massively parallel sequencing (MPS) using the Illumina MiSeq and 2) the iPLEX assay on the MassARRAY® System. On the Illumina MiSeq, the standard deviation of the fraction of methylation was under 3% between replicates, whereas the reproducibility of the MassARRAY® was very difficult to achieve. We tested the denaturation and conversion times of the Premium Bisulfite kit (Diagenode) in order to optimise the conversion rate and minimise DNA degradation. Finally, we tested the reproducibility of the methylation patterns using the Illumina MiSeq

In this study, we tested various approaches to remove DNA from hard laboratory surfaces. We contaminated clean surfaces with four different concentrations of massively parallel sequencing libraries. The DNA was dried and left for 45 min and for 24 h, respectively, before any treatment. The surfaces were cleaned with six different methods water, 96% ethanol, water followed by 96% ethanol, 3–6% hypochlorite solution, 0.9–1.8% hypochlorite solution, and no treatment. Subsequently, the surfaces were swabbed using cotton. DNA was extracted from the swabs and the DNA concentrations were determined by real-time PCR. The results showed that leaving the DNA traces for 24 h decreased the amount of amplifiable DNA approximately four times. The concentrations of amplifiable DNA were approximately five times lower after cleaning with 96% ethanol. Cleaning with water and water followed by 96% ethanol reduced the amount of amplifiable DNA 100–200 times, whereas cleaning with hypochlorite removed all traces of amplifiable DNA. This indicated that the ‘mechanical’ cleaning was efficient but cleaning with hypochlorite was superior. In conclusion, it is recommended to clean laboratory surfaces with 0.9–1.8% hypochlorite in order to eliminate laboratory contamination and simultaneously minimise any poisonous gases., In this study, we tested various approaches to remove DNA from hard laboratory surfaces. We contaminated clean surfaces with four different concentrations of massively parallel sequencing libraries. The DNA was dried and left for 45 min and for 24 h, respectively, before any treatment. The surfaces were cleaned with six different methods water, 96% ethanol, water followed by 96% ethanol, 3–6% hypochlorite solution, 0.9–1.8% hypochlorite solution, and no treatment. Subsequently, the surfaces were swabbed using cotton. DNA was extracted from the swabs and the DNA concentrations were determined by real-time PCR. The results showed that leaving the DNA traces for 24 h decreased the amount of amplifiable DNA approximately four times. The concentrations of amplifiable DNA were approximately five times lower after cleaning with 96% ethanol. Cleaning with water and water followed by 96% ethanol reduced the amount of amplifiable DNA 100–200 times, whereas cleaning with hypochlorite removed all traces of amplifiable DNA. This indicated that the ‘mechanical’ cleaning was efficient but cleaning with hypochlorite was superior. In conclusion, it is recommended to clean laboratory surfaces with 0.9–1.8% hypochlorite in order to eliminate laboratory contamination and simultaneously minimise any poisonous gases.

The Precision ID Identity Panel was used to type 109 Somali individuals in order to obtain allele frequencies for the Somali population. These frequencies were used to establish a Somali HID-SNP database, which will be used for the biostatistic calculations in family and immigration cases. Genotypes obtained with the Precision ID Identity Panel were found to be almost in complete concordance with genotypes obtained with the SNPforID PCR-SBE-CE assay. In seven SNP loci, silent alleles were identified, of which most were previously described in the literature. The project also set out to compare different AmpliSeq™ workflows to investigate the possibility of using automated library building in forensic genetic case work. In order to do so, the SNP typing of the Somalis was performed using three different workflows: 1) manual library building and sequencing on the Ion PGM™, 2) automated library building using the Biomek(®)3000 and sequencing on the Ion PGM™, and 3) automated library building using the Ion Chef™ and sequencing on the Ion S5™. AmpliSeq™ workflows were compared based on coverage, locus balance, noise, and heterozygote balance. Overall, the Ion Chef™/Ion S5™ workflow was found to give the best results and required least hands-on time in the laboratory. However, the Ion Chef™/Ion S5™ workflow was also the most expensive. The number of libraries that may be constructed in one Ion Chef™ library building run was limited to eight, which is too little for high throughput workflows. The Biomek(®)3000/Ion PGM™ workflow was found to perform similarly to the manual/Ion PGM™ workflow. This argues for the use of automated library building in forensic genetic case work. Automated library building decreases the workload of the laboratory staff, decreases the risk of pipetting errors, and simplifies the daily workflow in forensic genetic laboratories.

Surveillance of wild vertebrates can be challenging in remote and inaccessible areas such as tropical rainforests. Blood-feeding parasites, such as leeches, can facilitate wild vertebrate monitoring by targeting residual DNA from the animals the leeches feed on. Successes in detecting host DNA from leeches suggest that host viruses may also be detectable. To systematically test this hypothesis, we performed a proof of concept study using quantitative PCR (qPCR) to detect DNA viruses (bovine herpesvirus [BHV], human adenovirus [HAdV]) and RNA viruses (influenza A [InfA] and measles morbillivirus [MeV]) from nucleic acids extracted from medicinal leeches fed with blood spiked with each virus. All viruses except BHV showed a gradual decline in concentration from day 1 to 50, and all except BHV were detectable in at least half of the samples even after 50 days. BHV exhibited a rapid decline at day 27 and was undetectable at day 50. Our findings in medicinal leeches indicate that leeches collected in the wild might be an untapped resource for detecting vertebrate viruses and could provide new opportunities to study wildlife viral diseases of rare species in challenging environments, where capturing and handling of animals is difficult.

The HID-Ion AmpliSeq™ Identity Panel is a next-generation sequencing assay with 90 autosomal and 34 Y-chromosome SNPs that are amplified in one PCR step and subsequently sequenced using the Ion Personal Genome Machine (Ion PGM™) System. This assay was validated for relationship testing in our ISO 17025 accredited laboratory in 2015. Here, the essential parts of the validation report submitted to the Danish Accreditation Fund are presented. A total of 100 unrelated Danes were typed in duplicates and the locus balance, heterozygote balance (Hb) and noise levels were analysed in detail. Two loci were disregarded for casework because genotyping was uncertain. Hb for rs7520386 was skewed and high levels of noise were observed in rs576261. Three general acceptance criteria for analysis of single-source samples were defined: (i) sequencing depth > 200 reads, (ii) noise level < 3% and (iii) Hb > 0.3. A Python script named SNPonPGM was developed to assist the analyst by highlighting loci that do not fulfil the general acceptance criteria. Furthermore, SNPonPGM has functions that reduce the hands-on time of the reporting officer to a few minutes per case. Mixtures with DNA from two individuals in a 1:24 ratio were readily identified using the three criteria and the SNPonPGM script.

Here, we demonstrate that punches from buccal swab samples preserved on FTA cards can be used for high-throughput DNA sequencing, also known as massively parallel sequencing (MPS). We typed 44 reference samples with the HID-Ion AmpliSeq Identity Panel using washed 1.2 mm punches from FTA cards with buccal swabs and compared the results with those obtained with DNA extracted using the EZ1 DNA Investigator Kit. Concordant profiles were obtained for all samples. Our protocol includes simple punch, wash, and PCR steps, reducing cost and hands-on time in the laboratory. Furthermore, it facilitates automation of DNA sequencing.

To ensure efficient sequencing, the DNA of next-generation sequencing (NGS) libraries must be quantified correctly. Therefore, an accurate, sensitive and stable method for DNA quantification is crucial. In this study, seven different methods for DNA quantification were compared to each other by quantifying NGS libraries for the Ion TorrentTM and Illumina1 platforms as well as dsDNA oligos with known DNA concentrations. Rather large variations in library concentration estimates were observed. The differences between the highest and lowest concentration estimates varied with a factor of 5–100 depending on the library concentration. The Bioanalyzer, TapeStation and Qubit1 instruments gave concentrations closest to the expected when quantifying dsDNA oligos. At very low concentrations (2–4 pg/ul) only the Bioanalyzer could reliably quantify the dsDNA oligos., To ensure efficient sequencing, the DNA of next-generation sequencing (NGS) libraries must be quantified correctly. Therefore, an accurate, sensitive and stable method for DNA quantification is crucial. In this study, seven different methods for DNA quantification were compared to each other by quantifying NGS libraries for the Ion TorrentTM and Illumina1 platforms as well as dsDNA oligos with known DNA concentrations. Rather large variations in library concentration estimates were observed. The differences between the highest and lowest concentration estimates varied with a factor of 5–100 depending on the library concentration. The Bioanalyzer, TapeStation and Qubit1 instruments gave concentrations closest to the expected when quantifying dsDNA oligos. At very low concentrations (2–4 pg/ul) only the Bioanalyzer could reliably quantify the dsDNA oligos.

Solution hybridization capture methods utilize biotinylated oligonucleotides as baits to enrich homologous sequences from next generation sequencing (NGS) libraries. Coupled with NGS, the method generates kilo to gigabases of high confidence consensus targeted sequence. However, in many experiments, a non-negligible fraction of the resulting sequence reads are not homologous to the bait. We demonstrate that during capture, the bait-hybridized library molecules add additional flanking library sequences iteratively, such that baits limited to targeting relatively short regions (e.g. few hundred nucleotides) can result in enrichment across entire mitochondrial and bacterial genomes. Our findings suggest that some of the off-target sequences derived in capture experiments are non-randomly enriched, and that CapFlank will facilitate targeted enrichment of large contiguous sequences with minimal prior target sequence information.

Pathogens have always had a major interest to humans due to their central role in sickness and death. Influenza A annually kills at least 250,000 humans, and has been the cause of millions of further deaths during pandemic years in the past. Plague (Yersinia pestis) has been the cause of the Black Death that was leading to the desertion of whole cities, and as a result was for centuries one of the most feared events in human life. For both of these organisms we generally have situations with only very small amounts of pathogen nucleic acids available, usually because many interesting samples are degraded, as these samples are historic or generally patient samples. In this thesis we have principally explored methods to extract genetic information from such degraded samples, using the modern ‘high-throughput’ techniques of second generation sequencing (SGS) and mass spectrometry (MS), coupled with target enrichment. Using these techniques we are able to investigate the nucleic acids from these pathogens despite the poor quality or the small amounts. We have developed methods for the SGS of influenza A and employed this technique for patient samples to show intra-host variation together with investigation of the transmission chain. Furthermore, we have developed a new method for the subtyping of influenza A using proteins as these are not as fast degradable as nucleic acids and the method therefore is well-suited for old and degraded samples. Lastly, we have used SGS together with enrichment capture for the detection of Y. pestis in samples from the Justinian plague (600 AD) as an attempt to detect this pathogen as a cause of death in the victims.

Despite its charismatic appeal to both scientists and the general public, remarkably little is known about the giant squid Architeuthis, one of the largest of the invertebrates. Although specimens of Architeuthis are becoming more readily available owing to the advancement of deep-sea fishing techniques, considerable controversy exists with regard to topics as varied as their taxonomy, biology and even behaviour. In this study, we have characterized the mitochondrial genome (mitogenome) diversity of 43 Architeuthis samples collected from across the range of the species, in order to use genetic information to provide new and otherwise difficult to obtain insights into the life of this animal. The results show no detectable phylogenetic structure at the mitochondrial level and, furthermore, that the level of nucleotide diversity is exceptionally low. These observations are consistent with the hypotheses that there is only one global species of giant squid, Architeuthis dux (Steenstrup, 1857), and that it is highly vagile, possibly dispersing through both a drifting paralarval stage and migration of larger individuals. Demographic history analyses of the genetic data suggest that there has been a recent population expansion or selective sweep, which may explain the low level of genetic diversity.

Ribonucleic acids (RNA) are generally considered fragile molecules that are readily degraded. However, there is growing documentation of long-term (from days to centuries) RNA persistence in a variety of contexts and tissue types, and as such a number of academic disciplines are beginning to exploit degraded RNA. While the reasons for its survival are not fully understood, there are several plausible mechanisms that would safeguard this molecule against degradation. However, after examining the literature available on the postmortem instability and decay mechanisms of RNA, it has become clear that limited experimental studies and no reviews offer an overview of these mechanisms. Hence in this review we outline molecular reasons for RNA surviving long-term postmortem, and provide specific examples of RNA survival in forensic, archival and archaeological contexts. A better understanding of the mechanisms of RNA decay will be crucial for developing expectations on its long-term survival.

Influenza viruses such as swine-origin influenza A(H1N1) virus (A(H1N1)pdm09) generate genetic diversity due to the high error rate of their RNA polymerase, often resulting in mixed genotype populations (intra-host variants) within a single infection. This variation helps influenza to rapidly respond to selection pressures, such as those imposed by the immunological host response and antiviral therapy. We have applied deep sequencing to characterize influenza intra-host variation in a transmission chain consisting of three cases due to oseltamivir-sensitive viruses, and one derived oseltamivir-resistant case.

Niche variation owing to individual differences in ecology has been hypothesized to be an early stage of sympatric speciation. Yet to date, no study has tracked niche width over more than a few generations. In this study, we show the presence of isotopic niche variation over millennial timescales and investigate the evolutionary outcomes. Isotopic ratios were measured from tissue samples of sympatric killer whale Orcinus orca lineages from the North Sea, spanning over 10 000 years. Isotopic ratios spanned a range similar to the difference in isotopic values of two known prey items, herring Clupea harengus and harbour seal Phoca vitulina. Two proxies of the stage of speciation, lineage sorting of mitogenomes and genotypic clustering, were both weak to intermediate indicating that speciation has made little progress. Thus, our study confirms that even with the necessary ecological conditions, i.e. among-individual variation in ecology, it is difficult for sympatric speciation to progress in the face of gene flow. In contrast to some theoretical models, our empirical results suggest that sympatric speciation driven by among-individual differences in ecological niche is a slow process and may not reach completion. We argue that sympatric speciation is constrained in this system owing to the plastic nature of the behavioural traits under selection when hunting either mammals or fish.