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.

Bacteria and fungi in acute cholecystitis. A prospective study comparing next generation sequencing to culture

OBJECTIVES

Guidelines for antibiotic treatment of acute cholecystitis are based on studies using culture techniques for microbial identification. Microbial culture has well described limitations and more comprehensive data on the microbial spectrum may support adjustments of these recommendations. We used next generation sequencing to conduct a thorough microbiological characterization of bile-samples from patients with moderate and severe acute cholecystitis.

METHODS

We prospectively included patients with moderate and severe acute cholecystitis, undergoing percutaneous or perioperative drainage of the gall bladder. Bile samples were analyzed using both culture and deep sequencing of bacterial 16S rRNA and rpoB genes and the fungal ITS2-segment. Clinical details were evaluated by medical record review.

RESULTS

Thirty-six patients with moderate and severe acute cholecystitis were included. Bile from 31 (86%) of these contained bacteria (29) and/or fungi (5) as determined by sequencing. Culture identified only 40 (38%) of the 106 microbes identified by sequencing. In none of the 15 polymicrobial samples did culture detect all present microbes. Frequently identified bacteria often missed by culture included oral streptococci, anaerobic bacteria, enterococci and Enterobacteriaceae other than Klebsiella spp. and Escherichia coli.

CONCLUSIONS

Culture techniques display decreased sensitivity for the microbial diagnostics of acute cholecystitis leaving possible pathogens undetected.

Development and comprehensive evaluation of a noninvasive prenatal paternity testing method through a scaled trial

PURPOSE

To eliminate the miscarriage risks caused by traditional invasive sampling methods, we develop a noninvasive prenatal paternity testing (NIPPT) method and evaluate its efficiency, reliability and sensitivity based on a scaled trial.

METHODS

We use maternal cell-free DNA and massive parallel sequencing to obtain NIPPT genotypes for parents and fetuses based on quality-controlled genome-wide single nucleotide polymorphisms (SNPs). In a preliminary testing, data from 14 pregnant women and 7 negative controls are used for setting threshold of fetal genotyping in reference to postpartum children. After that, those from 349 cases with pregnancies of 6-35 gestational weeks (GW) and 9 negative controls from non-pregnant women who have fertility experience previously are in-depth evaluated.

RESULTS

In all cases, the biological fathers have been successfully identified from unrelated with a combined paternity index (CPI) of 3.58 × 10¹⁸ - 1.46 × 10¹⁶⁵ for the cases versus 1.52 × 10^-22 - 2.30 × 10^-839 for the controls. For negative controls, fetal SNPs originating from previous pregnancies could not be detected. Our NIPPT results completely aligned with the invasive prenatal test results using PCR-CE STR methods.

CONCLUSION

NIPPT can be applied to determine paternity accurately from 6 weeks after conception until birth and may serve as an alternative prenatal paternity test advantageous to the currently-used methods.

Full mtDNA genome sequencing of Brazilian admixed populations: A forensic-focused evaluation of a MPS application as an alternative to Sanger sequencing methods

The use of Massive Parallel Sequencing (MPS) techniques have been proposed by the forensic community as an alternative to Sanger sequencing methods in routine forensic casework analysis regarding mitochondrial DNA (mtDNA). Interesting features of MPS include high throughput, ability to simultaneously genotype a significant number of samples by barcoding techniques, processing automation, reduced time and costs, among others. Advantages include the capability of generating full mtDNA genome sequences versus usual techniques, usually limited to hypervariable or control regions exclusively. In this work, 96 reference single-source samples from three different Brazilian cities were subjected to full mtDNA genome sequencing by MPS techniques using an early-access version of Precision ID mtDNA Whole Genome Panel on an Ion Torrent PGM platform (Thermo Fisher Scientific, Waltham, MA, USA). Complete, high-quality sequences were obtained and sequencing performance was evaluated via four different metrics. As a subset of evaluated samples have been previously submitted for Sanger sequencing of the control region, a comparative analysis of both methods' results was conducted in order to compare technique adequacy within a forensic context. Even though this study is one of the first to report full mtDNA genome sequences for Brazilian admixed populations, the observed haplotypes exhibit a predominance of Native American and African maternal lineages in the studied sample set, reproducing results described in the literature for control regions only. Interpopulation analysis among Brazilian and 26 worldwide populations was also carried out. The results indicate that MPS-generated full mtDNA genome sequences may have great utility in forensic real casework applications, with a pronounced gain of genetic information and discrimination power provided by coding region evaluation and the enhanced capacity of heteroplasmies determination. Database construction and other relevant factors concerning implementation of such techniques in Brazilian forensic laboratories are also discussed.

Detection of Heteroplasmic Variants in the Mitochondrial Genome through Massive Parallel Sequencing

Detecting heteroplasmies in the mitochondrial DNA (mtDNA) has been a challenge for many years. In the past, Sanger sequencing was the main option to perform this analysis, however, this method could not detect low frequency heteroplasmies. Massive Parallel Sequencing (MPS) provides the opportunity to study the mtDNA in depth, but a controlled pipeline is necessary to reliably retrieve and quantify the low frequency variants. It has been shown that differences in methods can significantly affect the number and frequency of the retrieved variants. In this protocol, we present a method involving both wet lab and bioinformatics that allows identifying and quantifying single nucleotide variants in the full mtDNA sequence, down to a heteroplasmic load of 1.5%. For this, we set up a PCR-based amplification of the mtDNA, followed by MPS using Illumina chemistry, and variant calling with two different algorithms, mtDNA server and Mutect. The PCR amplification is used to enrich the mitochondrial fraction, while the bioinformatic processing with two algorithms is used to discriminate the true heteroplasmies from background noise. The protocol described here allows for deep sequencing of the mitochondrial DNA in bulk DNA samples as well as single cells (both large cells such as human oocytes, and small-sized single cells such as human embryonic stem cells) with minor modifications to the protocol.

Complex multisystem phenotype associated with the mitochondrial DNA m.5522G>A mutation

Mitochondrial tRNAs are responsible for more than half of pathogenic point mutations in the mitochondrial genome (mtDNA). Different mutations give rise to widely differing phenotypes, ranging from isolated organ-specific diseases to multisystem conditions. Herein, we report a 40-year-old woman presenting with a complex multisystem phenotype including sensorineural hearing loss, retinopathy, severe dilated cardiomyopathy, non-insulin dependent diabetes mellitus, and renal failure. Sequence analysis of mtDNA identified the m.5522G>A mutation in MT-TW, the gene encoding mitochondrial tRNA for tryptophan. The heteroplasmic variant, thus far described once, was almost exclusively confined to skeletal muscle tissue, as shown by massive parallel sequencing and corroborated by an ad hoc designed PCR-based strategy. This patient, presenting a severe, multisystem involvement apparently sparing the brain, contributes to the genetic heterogeneity of mitochondrial diseases caused by mutations in mitochondrial tRNAs.

HLA-G, -E and -F regulatory and coding region variability and haplotypes in the Beninese Toffin population sample

HLA-G/E/F genes exhibit immunomodulatory properties and are expressed in placenta. Little attention has been devoted to the study of these genes in sub-Saharan African populations, which are yet the most diverse. To fill this gap, we evaluated the complete gene variability, approximately 5.1 kb for HLA-G (n = 149), 7.7 kb for HLA-E (n = 150) and 6.2 kb for HLA-F (n = 152) in the remote Beninese Toffin population, using massive parallel sequencing. Overall, 96, 37 and 68 variable sites were detected along the entire HLA-G, -E and -F, respectively, arranged into region-specific haplotypes; i.e., promoter haplotypes (16, 19, and 15 respectively), coding haplotypes (19, 15, and 29 respectively), 3' untranslated region (3'UTR) haplotypes (12, 7 and 2, respectively) and extended haplotypes (33, 31 and 32 respectively). All promoter/coding/3'UTR haplotypes followed the patterns already described in worldwide populations. HLA-E was the most conserved, exhibiting mainly two full-length encoded-molecules (E*01:01 and E*01:03), followed by HLA-F, three full-length proteins (F*01:01, F*01:02 and F*01:03) and HLA-G, four proteins: three full-length (G*01:01, G*01:03 and G*01:04) and one truncated (G*01:05N). Although HLA-G/E/F alleles in the Toffin population were the most frequently observed worldwide, the frequencies of the coding haplotypes were closely similar to those described for other African populations (Guinea-Conakry and Burkina-Faso), when compared to non-African ones (Brazilian), indicating that variable sites along these genes were present in Africa before human dispersion.

Sudden unexpected death in the young - Value of massive parallel sequencing in postmortem genetic analyses

Cases of sudden cardiac death (SCD) in young and apparently healthy individuals represent a devastating event in affected families. Hereditary arrhythmia syndromes, which include primary electrical heart disorders as well as cardiomyopathies, are known to contribute to a significant number of these sudden death cases. We performed postmortem genetic analyses in young sudden death cases (aged <45years) by means of a defined gene panel using massive parallel sequencing (MPS). The data were evaluated bioinformatically and detected sequence variants were assessed using common databases and applying in silico prediction tools. In this study, we identified variants with likely pathogenic effect in 6 of 9 sudden unexpected death (SUD) cases. Due to the detection of numerous unknown and unclassified variants, interpretation of the results proved to be challenging. However, by means of an appropriate evaluation of the findings, MPS represents an important tool to support the forensic investigation and implies great progress for relatives of young SCD victims facilitating adequate risk stratification and genetic counseling.

Microbial forensics: new breakthroughs and future prospects

Recent advances in genetic data generation, through massive parallel sequencing (MPS), storage and analysis have fostered significant progresses in microbial forensics (or forensic microbiology). Initial applications in circumstances of biocrime, bioterrorism and epidemiology are now accompanied by the prospect of using microorganisms (i) as ancillary evidence in criminal cases; (ii) to clarify causes of death (e.g., drownings, toxicology, hospital-acquired infections, sudden infant death and shaken baby syndromes); (iii) to assist human identification (skin, hair and body fluid microbiomes); (iv) for geolocation (soil microbiome); and (v) to estimate postmortem interval (thanatomicrobiome and epinecrotic microbial community). When compared with classical microbiological methods, MPS offers a diverse range of advantages and alternative possibilities. However, prior to its implementation in the forensic context, critical efforts concerning the elaboration of standards and guidelines consolidated by the creation of robust and comprehensive reference databases must be undertaken.

Optimizing body fluid recognition from microbial taxonomic profiles

In forensics the DNA-profile is used to identify the person who left a biological trace, but information on body fluid can also be essential in the evidence evaluation process. Microbial composition data could potentially be used for body fluid recognition as an improved alternative to the currently used presumptive tests. We have developed a customized workflow for interpretation of bacterial 16S sequence data based on a model composed of Partial Least Squares (PLS) in combination with Linear Discriminant Analysis (LDA). Large data sets from the Human Microbiome Project (HMP) and the American Gut Project (AGP) were used to test different settings in order to optimize performance. From the initial cross-validation of body fluid recognition within the HMP data, the optimal overall accuracy was close to 98%. Sensitivity values for the fecal and oral samples were ≥0.99, followed by the vaginal samples with 0.98 and the skin and nasal samples with 0.96 and 0.81 respectively. Specificity values were high for all 5 categories, mostly >0.99. This optimal performance was achieved by using the following settings: Taxonomic profiles based on operational taxonomic units (OTUs) with 0.98 identity (OTU98), Aitchisons simplex transform with C = 1 pseudo-count and no regularization (r = 1) in the PLS step. Variable selection did not improve the performance further. To test for robustness across sequencing platforms, we also trained the classifier on HMP data and tested on the AGP data set. In this case, the standard OTU based approach showed moderately decline in accuracy. However, by using taxonomic profiles made by direct assignment of reads to a genus, we were able to nearly maintain the high accuracy levels. The optimal combination of settings was still used, except the taxonomic level being genus instead of OTU98. The performance may be improved even further by using higher resolution taxonomic bins.

Improving the diagnosis of cobalamin and related defects by genomic analysis, plus functional and structural assessment of novel variants

Background

Cellular cobalamin defects are a locus and allelic heterogeneous disorder. The gold standard for coming to genetic diagnoses of cobalamin defects has for some time been gene-by-gene Sanger sequencing of individual DNA fragments. Enzymatic and cellular methods are employed before such sequencing to help in the selection of the gene defects to be sought, but this is time-consuming and laborious. Furthermore some cases remain undiagnosed because no biochemical methods have been available to test for cobalamin absorption and transport defects.

Results

This paper reports the use of massive parallel sequencing of DNA (exome analysis) for the accurate and rapid genetic diagnosis of cobalamin-related defects in a cohort of affected patients. The method was first validated in an initial cohort with different cobalamin defects. Mendelian segregation, the frequency of mutations, and the comprehensive structural and functional analysis of gene variants, identified disease-causing mutations in 12 genes involved in the absorption and synthesis of active cofactors of vitamin B₁₂ (22 cases), and in the non-cobalamin metabolism-related genes ACSF3 (in four biochemically misdiagnosed patients) and SUCLA2 (in one patient with an unusual presentation). We have identified thirteen new variants all classified as pathogenic according to the ACGM recommendation but four were classified as variant likely pathogenic in MUT and SUCLA2. Functional and structural analysis provided evidences to classify them as pathogenic variants.

Conclusions

The present findings suggest that the technology used is sufficiently sensitive and specific, and the results it provides sufficiently reproducible, to recommend its use as a second-tier test after the biochemical detection of cobalamin disorder markers in the first days of life. However, for accurate diagnoses to be made, biochemical and functional tests that allow comprehensive clinical phenotyping are also needed.

Electronic supplementary material

The online version of this article (10.1186/s13023-018-0862-y) contains supplementary material, which is available to authorized users.

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.

Weekly variations in norovirus genogroup II genotypes in Japanese oysters

Increased levels of norovirus contamination in oysters were reportedly associated with a gastroenteritis epidemic occurring upstream of an oyster farming area. In this study, we monitored the norovirus concentration in oysters weekly between November 2014 and March 2015 and investigated the statistical relationship between norovirus genogroup II (GII) concentrations in oyster and sewage samples and the number of gastroenteritis cases in the area using cross-correlation analysis. A peak correlation coefficient (R = 0.76) at a time lag of +1 week was observed between the number of gastroenteritis cases and norovirus GII concentrations in oysters, indicating that oyster contamination is correlated with the number of gastroenteritis cases with a 1-week delay. Moreover, weekly variations in norovirus GII genotypes in oysters were evaluated using pyrosequencing. Only GII.3 was detected in November and December 2014, whereas GII.17 and GII.4 were present from January to March 2015. GII.17 Kawasaki 2014 strains were detected more frequently than GII.4 Sydney 2012 strains in oyster samples, as previously observed in stool and sewage samples collected during the same study period in Miyagi, Japan. Our observations indicate that there is a time lag between the circulation of norovirus genotypes in the human population and the detection of those genotypes in oysters.

Pancreatic neuroendocrine carcinomas reveal a closer relationship to ductal adenocarcinomas than to neuroendocrine tumors G3

Pancreatic neuroendocrine carcinoma is a rare aggressive tumor commonly harboring TP53 and RB1 alterations and lacking neuroendocrine-related genetic changes such as mutations in MEN1 and ATRX/DAXX. Little is known about its genetic profile with regard to that of pancreatic ductal adenocarcinoma. We therefore conducted a detailed genetic study in 12 pancreatic neuroendocrine carcinomas of large cell (n = 9) and small cell type (n = 3) using massive parallel sequencing applying a 409-gene panel on an Ion Torrent system. The genetic data were compared with known data of pancreatic ductal adenocarcinoma and correlated with exocrine lineage marker expression. A similar analysis was performed in 11 pancreatic neuroendocrine tumors G3. Neuroendocrine carcinomas harbored 63 somatic mutations in 45 different genes, affecting most commonly TP53 (8/12 cases), KRAS (5/12 cases), and RB1 (loss of expression with or without deletion in 4/12 cases). Five carcinomas had both TP53 and KRAS mutations. Neuroendocrine tumors G3 only shared singular mutations in 5 different genes with neuroendocrine carcinomas, including TP53, CDKN2A, ARID1A, LRP1B, and APC, affecting 5 different cases. Most KRAS-positive neuroendocrine carcinomas also expressed MUC1 (4/5) and carcinoembryonic antigen (3/5) as markers of ductal differentiation. Our data indicate that almost half of the pancreatic neuroendocrine carcinomas are genetically and phenotypically related to pancreatic ductal adenocarcinoma, and might therefore respond to chemotherapies targeting the latter carcinomas.

Performance and microbial community structure of a polar Arctic Circle aerobic granular sludge system operating at low temperature

The aim of this work was to study the performance and microbial community structure of a polar Arctic Circle aerobic granular sludge (AGS) system operating at low temperature. Thus, an AGS bioreactor was operated at 7, 5 and 3 °C of temperature using a cold-adapted sludge from Lapland. At 5 °C, it yielded acceptable conversion rates, in terms of nitrogen, phosphorous, and organic matter. However, under 3 °C a negligible nitrogen and phosphorous removal performance was observed. Below 5 °C, scanning electron microscopy studies showed a wispy, non-dense and irregular granular structure with a strong outgrowth of filamentous. Moreover, Illumina next-generation sequencing showed a heterogeneous microbial population where SM1K20 (Archaea), Trichosporon domesticum (Fungus), and Zooglea, Arcobacter and Acinetobacter (Bacteria) were the dominant phylotypes. Our study suggests that AGS technologies inoculated with North Pole sludge could be operated, in cold regions for a period longer than 3 months (winter season) under 5 °C of water temperature.

MicroRNA Expression Analysis: Next-Generation Sequencing

MicroRNAs are a class of small noncoding RNAs that function as regulators involving in many biological processes. The evaluation of miRNAs and their targets has been aided by miRNA expression profiling studies including multiplex PCR, microarrays, and recent next-generation sequencing tools. Next-generation sequencing has enabled us to profile thousands of genes in a single experiment and overcome the background signal and cross-hybridization issues of microarrays. Next-generation sequencing also allows for the simultaneous confirmation of known miRNAs and discovery of new miRNAs, and significantly reduces costs while providing billions of nucleotide information within a single experiment. Here we describe a detailed procedure of generation of miRNA library for next-generation sequencing to increase the efficiency of adapter ligation and finally construct a more specific cDNA library for sequencing and analyses for miRNA expression profiling.

Analysis of mainland Japanese and Okinawan Japanese populations using the precision ID Ancestry Panel

We typed 165 AIMs in 49 mainland Japanese and 47 Okinawa Japanese using the Precision ID Ancestry Panel (Thermo Fisher Scientific). None of the 165 SNPs showed significant deviation from Hardy-Weinberg equilibrium in the mainland Japanese. One SNP (rs3943253) showed significant deviation from Hardy-Weinberg equilibrium in Okinawa Japanese. Fisher's exact tests showed that the genotype frequencies of 14 loci were significantly different (p<0.05) between the two populations before correction for multiple testing. After Bonferroni correction, only rs671 remained statistically significant (p<0.0003). This SNP is located in the ALDH2 gene. The mutant A allele is associated with increased side effects after alcohol intake. The frequency of the GG genotype (wild type) was higher in the Okinawa Japanese (78.7%) than in mainland Japanese (34.7%; Bonferroni corrected P<0.001). For 31 (63.3%) of the mainland Japanese and 42 (89.4%) of Okinawa Japanese, the highest population likelihood was obtained with the Japanese reference population. However, only in a few individuals, the likelihoods were significantly different from those calculated using reference data from neighboring populations. The likelihoods for mainland Japanese and Okinawa Japanese were not significantly different from each other for any of the investigated individuals. STRUCTURE and PCA analyses showed that mainland Japanese, Okinawa Japanese, and East Asians could not be differentiated with the Precision ID Ancestry Panel.

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.

Postmortem genetic testing should be recommended in sudden cardiac death cases due to thoracic aortic dissection

BACKGROUND

Acute thoracic aortic dissections and ruptures, the main life-threatening complications of the corresponding aneurysms, are an important cause of sudden cardiac death. Despite the usefulness of the molecular diagnosis of these conditions in the clinical setting, the corresponding forensic field remains largely unexplored. The main goal of this study was to explore and validate a new massive parallel sequencing candidate gene​ assay as a diagnostic tool for acute thoracic aortic dissection autopsy cases.

MATERIALS AND METHODS

Massive parallel sequencing of 22 thoracic aortic disease candidate genes performed in 17 cases of thoracic aortic dissection using AmpliSeq and Ion Proton technologies. Genetic variants were filtered by location, type, and frequency at the Exome Aggregation Consortium and an internal database and further classified based on the American College of Medical Genetics and Genomics (ACMG) recommendations published in 2015. All prioritized results were confirmed by traditional sequencing.

RESULTS

From the total of 10 potentially pathogenic genetic variants identified in 7 out of the 17 initial samples, 2 of them were further classified as pathogenic, 2 as likely pathogenic, 1 as possibly benign, and the remaining 5 as variants of uncertain significance, reaching a molecular autopsy yield of 23%, approximately.

CONCLUSIONS

This massive parallel sequencing candidate gene approach proved useful for the molecular autopsy of aortic dissection sudden cardiac death cases and should therefore be progressively incorporated into the forensic field, being especially beneficial for the anticipated diagnosis and risk stratification of any other family member at risk of developing the same condition.

Statistical method to compare massive parallel sequencing pipelines

Background

Today, sequencing is frequently carried out by Massive Parallel Sequencing (MPS) that cuts drastically sequencing time and expenses. Nevertheless, Sanger sequencing remains the main validation method to confirm the presence of variants. The analysis of MPS data involves the development of several bioinformatic tools, academic or commercial. We present here a statistical method to compare MPS pipelines and test it in a comparison between an academic (BWA-GATK) and a commercial pipeline (TMAP-NextGENe®), with and without reference to a gold standard (here, Sanger sequencing), on a panel of 41 genes in 43 epileptic patients. This method used the number of variants to fit log-linear models for pairwise agreements between pipelines. To assess the heterogeneity of the margins and the odds ratios of agreement, four log-linear models were used: a full model, a homogeneous-margin model, a model with single odds ratio for all patients, and a model with single intercept. Then a log-linear mixed model was fitted considering the biological variability as a random effect.

Results

Among the 390,339 base-pairs sequenced, TMAP-NextGENe® and BWA-GATK found, on average, 2253.49 and 1857.14 variants (single nucleotide variants and indels), respectively. Against the gold standard, the pipelines had similar sensitivities (63.47% vs. 63.42%) and close but significantly different specificities (99.57% vs. 99.65%; p < 0.001). Same-trend results were obtained when only single nucleotide variants were considered (99.98% specificity and 76.81% sensitivity for both pipelines).

Conclusions

The method allows thus pipeline comparison and selection. It is generalizable to all types of MPS data and all pipelines.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1552-9) contains supplementary material, which is available to authorized users.

Intellectual Disability & Rare Disorders: A Diagnostic Challenge

Rare disorders constitute a large and heterogeneous group of diagnoses of which many cause chronic disabilities with significant impact on the lives of affected individuals and their families as well as on the health-care system. Each individual disorder is rare, but when considered as a group, rare disorders are common with a total prevalence of approximately 6-8%. The clinical presentation of these disorders includes a broad diversity of symptoms and signs, often involving the nervous system and resulting in symptoms such as intellectual disability, neuropsychiatric disorders, epilepsy and motor dysfunction. The methods for establishing an etiological diagnosis in patients with rare disorders have improved dramatically during recent years. With the introduction of genomic screening methods, it has been shown that the cause is genetic in the majority of the patients and many will receive an etiological diagnosis in a clinical setting. However, there are a lot of challenges in diagnosing these disorders and despite recent years' advances, a large number of patients with rare disorders still go without an etiological diagnosis. In this chapter we will review the etiology of rare disorders with focus on intellectual disability and what has been learned from massive parallel sequencing studies in deciphering the genetic basis. Furthermore, we will discuss challenges in the etiological diagnostics of these disorders including issues that regard interpretation of the numerous genetic variants detected by genomic screening methods and challenges in the translation of massive parallel sequencing technologies into clinical practice.

Quantitative Genome-Wide Measurements of Meiotic DNA Double-Strand Breaks and Protein Binding in S. pombe

The fission yeast Schizosaccharomyces pombe is especially well suited for studying meiosis in molecular detail. Experiments with S. pombe strains that undergo a nearly synchronous meiosis-at variable temperatures-have elucidated the mechanisms of meiotic progression and the proteins that are involved. For example, studies focused on the initiation of meiotic recombination by programmed DNA double-strand breaks (DSBs) have proven exceptionally informative. In meiosis, some regions of DNA have more frequent DSBs than the surrounding regions. These DSB hotspots can be visualized by Southern blot hybridization of restriction fragments ranging from kilobases (kb) to megabases (Mb) in size. More recently, the benefits of genome-wide analysis to map the distribution and frequency of meiotic DSBs have been attained, with resolution down to the nucleotide level. Infrequent, non-hotspot DSBs previously not detectable have been observed, creating a better understanding of how recombination is regulated. Additional genome-wide analyses have shown proteins that bind specifically to DSB hotspots, providing insight into how the DSB initiation complex functions. We describe here detailed methods for achieving these results.

Histone Modifications in Major Depressive Disorder and Related Rodent Models

Major depressive disorder (MDD) is a multifactorial disease, weakly linked to multiple genetic risk factors. In contrast to that, environmental factors and "gene × environment" interaction between specific risk genes and environmental factors, such as severe or early stress exposure, have been strongly linked to MDD vulnerability. Stressors can act on the interface between an organism and the environment, the epigenome. The molecular foundation for the impact of stressors on the risk to develop MDD is based on the hormonal stress response itself: the glucocorticoid receptor (GR, encoded by NR3C1). NR3C1 can directly interact with the epigenome in the cell nucleus. Besides DNA methylation, histone modifications have been reported to be crucial targets for the interaction with the stress response system. Here, we review critical findings on the impact of the most relevant histone modifications, i.e. histone acetylation and methylation, in the context of MDD and related animal models. We discuss new treatment options which have been based on these findings, including histone deacetylase inhibitors (HDACis) and drugs targeting specific histone marks, closely linked to psychiatric disease. In this context we talk about contemporary and future approaches required to fully understand (1) the epigenetics of stress-related disease and (2) the mode of action of potential MDD drugs targeting histone modifications. This includes harnessing the unprecedented potentials of genome-wide analysis of the epigenome and transcriptome, in a cell type-specific manner, and the use of epigenome editing technologies to clearly link epigenetic marks on specific genomic loci to functional relevance.

Exome sequencing identifies potential novel candidate genes in patients with unexplained colorectal adenomatous polyposis

In up to 30% of patients with colorectal adenomatous polyposis, no germline mutation in the known genes APC, causing familial adenomatous polyposis, MUTYH, causing MUTYH-associated polyposis, and POLE or POLD1, causing Polymerase-Proofreading-associated polyposis can be identified, although a hereditary etiology is likely. To uncover new causative genes, exome sequencing was performed using DNA from leukocytes and a total of 12 colorectal adenomas from seven unrelated patients with unexplained sporadic adenomatous polyposis. For data analysis and variant filtering, an established bioinformatics pipeline including in-house tools was applied. Variants were filtered for rare truncating point mutations and copy-number variants assuming a dominant, recessive, or tumor suppressor model of inheritance. Subsequently, targeted sequence analysis of the most promising candidate genes was performed in a validation cohort of 191 unrelated patients. All relevant variants were validated by Sanger sequencing. The analysis of exome sequencing data resulted in the identification of rare loss-of-function germline mutations in three promising candidate genes (DSC2, PIEZO1, ZSWIM7). In the validation cohort, further variants predicted to be pathogenic were identified in DSC2 and PIEZO1. According to the somatic mutation spectra, the adenomas in this patient cohort follow the classical pathways of colorectal tumorigenesis. The present study identified three candidate genes which might represent rare causes for a predisposition to colorectal adenoma formation. Especially PIEZO1 (FAM38A) and ZSWIM7 (SWS1) warrant further exploration. To evaluate the clinical relevance of these genes, investigation of larger patient cohorts and functional studies are required.

NGS-QC Generator: A Quality Control System for ChIP-Seq and Related Deep Sequencing-Generated Datasets

The combination of massive parallel sequencing with a variety of modern DNA/RNA enrichment technologies provides means for interrogating functional protein-genome interactions (ChIP-seq), genome-wide transcriptional activity (RNA-seq; GRO-seq), chromatin accessibility (DNase-seq, FAIRE-seq, MNase-seq), and more recently the three-dimensional organization of chromatin (Hi-C, ChIA-PET). In systems biology-based approaches several of these readouts are generally cumulated with the aim of describing living systems through a reconstitution of the genome-regulatory functions. However, an issue that is often underestimated is that conclusions drawn from such multidimensional analyses of NGS-derived datasets critically depend on the quality of the compared datasets. To address this problem, we have developed the NGS-QC Generator, a quality control system that infers quality descriptors for any kind of ChIP-sequencing and related datasets. In this chapter we provide a detailed protocol for (1) assessing quality descriptors with the NGS-QC Generator; (2) to interpret the generated reports; and (3) to explore the database of QC indicators (www.ngs-qc.org) for >21,000 publicly available datasets.