Mark Gosink

Archival formalin-fixed paraffin-embedded (FFPE) tissue samples offer a vast but largely untapped resource for genomic research. The primary technical issues limiting use of FFPE samples are RNA yield and quality. In this study, we evaluated methods to demodify RNA highly fragmented and crosslinked by formalin fixation. Primary endpoints were RNA recovery, RNA-sequencing quality metrics, and transcriptional responses to a reference chemical (phenobarbital, PB). Frozen mouse liver samples from control and PB groups (n = 6/group) were divided and preserved for 3 months as follows: frozen (FR); 70% ethanol (OH); 10% buffered formalin for 18 h followed by ethanol (18F); or 10% buffered formalin (3F). Samples from OH, 18F, and 3F groups were processed to FFPE blocks and sectioned for RNA isolation. Additional sections from 3F received the following demodification protocols to mitigate RNA damage: short heated incubation with Tris-Acetate-EDTA buffer; overnight heated incubation with an organocatalyst using 2 different isolation kits; or overnight heated incubation without organocatalyst. Ribo-depleted, stranded, total RNA libraries were built and sequenced using the Illumina HiSeq 2500 platform. Overnight incubation (± organocatalyst) increased RNA yield >3-fold and RNA integrity numbers and fragment analysis values by > 1.5- and >3.0-fold, respectively, versus 3F. Postsequencing metrics also showed reduced bias in gene coverage and deletion rates for overnight incubation groups. All demodification groups had increased overlap for differentially expressed genes (77%–84%) and enriched pathways (91%–97%) with FR, with the highest overlap in the organocatalyst groups. These results demonstrate simple changes in RNA isolation methods that can enhance genomic analyses of FFPE samples.

Drug-induced vascular injury (DIVI) in preclinical studies can delay, if not terminate, a drug development program. Clinical detection of DIVI can be very difficult as there are no definitive biomarkers known to reliably detect this disorder in all instances. The preclinical identification of DIVI requires detailed microscopic examination of a wide range of tissues although one of the most commonly affected areas in rats is the mesenteric vasculature. The reason for this predisposition of mesenteric arteries in rats as well as the exact mechanism and cell types involved in the initial development of these lesions have not been fully elucidated. We hypothesized that by using a mixed culture of cells from rat mesenteric tissue, we would be able to identify an RNA expression signature that could predict the in vivo development of DIVI. Five compounds designed to inhibit Phosphodiesterase 4 activity (PDE4i) were chosen as positive controls. PDE4i’s are well known to induce DIVI in the mesenteric vasculature of rats and there is microscopic evidence that this is associated, at least in part, with a proinflammatory mechanism. We surveyed, by qRT-PCR, the expression of 96 genes known to be involved in inflammation and using a Random-Forest model, identified 12 genes predictive of in vivo DIVI outcomes in rats. Using these genes, we were able to cross-validate the ability of the Random-Forest modeling to predict the concentration at which PDE4i caused DIVI in vivo.

MicroRNAs (miRNAs) released into the peripheral circulation upon cellular injury have shown a promise as a new class of tissue-specific biomarkers. We were first to demonstrate that next-generation sequencing analysis of serum from human subjects with acetaminophen-induced liver injury revealed a specific signature of circulating miRNAs. We consequently hypothesized that different types of hepatic liver impairments might feature distinct signatures of circulating miRNAs and that this approach might be useful as minimally invasive diagnostic “liquid biopsies” enabling the interrogation of underlying molecular mechanisms of injury in distant tissues. Therefore we examined serum circulating miRNAs in a total of 72 serum samples from a group of 53 subjects that included patients with accidental acetaminophen overdose, hepatitis B infection, liver cirrhosis and type 2 diabetes as well as gender- and age-matched healthy subjects with no evidence of liver disease. The miRNA signatures were identified using next-generation sequencing that provided analysis for the whole miRNome, including miRNA isoforms. Compared to the healthy subjects, a total of 179 miRNAs showed altered serum levels across the diseased subjects. Although many subjects have elevated alanine aminotransferase suggesting liver impairments, we identified distinct miRNA signatures for different impairments with minimum overlap. Furthermore, the bioinformatics analysis of miRNA signatures revealed relevant molecular pathways associated with the mechanisms of toxicity and or pathogenesis of disease. Interestingly, the high proportion of miRNA isoforms present in the respective signatures indicated a new level of complexity in cellular response to stress or disease. Our study demonstrates for the first time that signatures of circulating miRNAs show specificity for liver injury phenotypes and, once validated, might become useful for diagnosis of organ pathologies as “liquid biopsies”.

A number of chemical compounds have been shown to induce liver tumors in mice but not in other species. While several mechanisms for this species-specific tumorigenicity have been proposed, no definitive mechanism has been established. We examined the effects of the nongenotoxic rodent hepatic carcinogen, WY-14,643, in male mice from a high liver tumor susceptible strain (C3H/HeJ), and from a low tumor susceptible strain (C57BL/6). WY-14,643, a PPARα activator induced widespread increases in the expression of some endogenous retroelements, namely members of LTR and LINE elements in both strains. The expression of a number of known retroviral defense genes was also elevated. We also demonstrated that basal immune-mediated viral defense was elevated in C57BL/6 mice (the resistant strain) and that WY-14,643 further activated those immuno-defense processes. We propose that the previously reported >100X activity of retroelements in mice drives mouse-specific tumorigenicity. We also propose that C57BL/6’s competent immune to retroviral activation allows it to remove cells before the activation of these elements can result in significant chromosomal insertions and mutation. Finally, we showed that WY-14,643 treatment induced gene signatures of DNA recombination in the sensitive C3H/HeJ strain.

Attempts are frequently made to investigate adverse findings from preclinical toxicology studies in order to better understand underlying toxicity mechanisms. These efforts often begin with limited information, including a description of the adverse finding, knowledge of the structure of the chemical associated with its cause and the intended pharmacological target. ToxEvaluator was developed jointly by Pfizer and the Comparative Toxicogenomics Database (http://ctdbase.org) team at North Carolina State University as an in silico platform to facilitate interpretation of toxicity findings in light of prior knowledge. Through the integration of a diverse set of in silico tools that leverage a number of public and proprietary databases, ToxEvaluator streamlines the process of aggregating and interrogating diverse sources of information. The user enters compound and target identifiers, and selects adverse event descriptors from a safety lexicon and mapped MeSH disease terms. ToxEvaluator provides a summary report with multiple distinct areas organized according to what target or structural aspects have been linked to the adverse finding, including primary pharmacology, structurally similar proprietary compounds, structurally similar public domain compounds, predicted secondary (i.e. off-target) pharmacology and known secondary pharmacology. Similar proprietary compounds and their associated in vivo toxicity findings are reported, along with a link to relevant supporting documents. For similar public domain compounds and interacting targets, ToxEvaluator integrates relationships curated in Comparative Toxicogenomics Database, returning all direct and inferred linkages between them. As an example of its utility, we demonstrate how ToxEvaluator rapidly identified direct (primary pharmacology) and indirect (secondary pharmacology) linkages between cerivastatin and myopathy.

One of the many roles of a toxicologist is to determine if an observed adverse event (AE) is related to a previously unrecognized function of a given gene/protein. Towards that end, he or she will search a variety of public and propriety databases for information linking that protein to the observed AE. However, these databases tend to present all available information about a protein, which can be overwhelming, limiting the ability to find information about the specific toxicity being investigated. ToxReporter compiles information from a broad selection of resources and limits display of the information to user-selected areas of interest. ToxReporter is a PERL-based web-application which utilizes a MySQL database to streamline this process by categorizing public and proprietary domain-derived information into predefined safety categories according to a customizable lexicon. Users can view gene information that is ‘red-flagged’ according to the safety issue under investigation. ToxReporter also uses a scoring system based on relative counts of the red-flags to rank all genes for the amount of information pertaining to each safety issue and to display their scored ranking as an easily interpretable ‘Tox-At-A-Glance’ chart. Although ToxReporter was originally developed to display safety information, its flexible design could easily be adapted to display disease information as well. Database URL: ToxReporter is freely available at https://github.com/mgosink/ToxReporter.

MicroRNAs (miRNAs) are small, non-coding RNAs that regulate protein levels post-transcriptionally. miRNAs play important regulatory roles in many cellular processes and have been implicated in several diseases. Recent studies have reported significant levels of miRNAs in a variety of body fluids, raising the possibility that miRNAs could serve as useful biomarkers. Next-generation sequencing (NGS) is increasingly employed in biomedical investigations. Although concordance between this platform and qRT-PCR based assays has been reported in high quality specimens, information is lacking on comparisons in biofluids especially urine. Here we describe the changes in miRNA expression patterns in a rodent model of renal tubular injury (gentamicin). Our aim is to compare RNA sequencing and qPCR based miRNA profiling in urine specimen from control and rats with confirmed tubular injury. Our preliminary examination of the concordance between miRNA-seq and qRT-PCR in urine specimen suggests minimal agreement between platforms probably due to the differences in sensitivity. Our results suggest that although miRNA-seq has superior specificity, it may not detect low abundant miRNAs in urine samples. Specifically, miRNA-seq did not detect some sequences which were identified by qRT-PCR. On the other hand, the qRT-PCR analysis was not able to detect the miRNA isoforms, which made up the majority of miRNA changes detected by NGS. To our knowledge, this is the first time that miRNA profiling platforms including NGS have been compared in urine specimen. miRNAs identified by both platforms, let-7d, miR-203, and miR-320, may potentially serve as promising novel urinary biomarkers for drug induced renal tubular epithelial injury.

Background: Austism spectrum disorder (ASD) is a heterogeneous behavioral disorder or condition characterized by severe impairment of social engagement and the presence of repetitive activities. The molecular etiology of ASD is still largely unknown despite a strong genetic component. Part of the difficulty in turning genetics into disease mechanisms and potentially new therapeutics is the sheer number and diversity of the genes that have been associated with ASD and ASD symptoms. The goal of this work is to use shRNA-generated models of genetic defects proposed as causative for ASD to identify the common pathways that might explain how they produce a core clinical disability.Methods: Transcript levels ofMecp2,Mef2a,Mef2d,Fmr1,Nlgn1,Nlgn3,Pten, andShank3were knocked-down in mouse primary neuron cultures using shRNA constructs. Whole genome expression analysis was conducted for each of the knockdown cultures as well as a mock-transduced culture and a culture exposed to a lentivirus expressing an anti-luciferase shRNA. Gene set enrichment and a causal reasoning engine was employed to identify pathway level perturbations generated by the transcript knockdown.Results: Quantification of the shRNA targets confirmed the successful knockdown at the transcript and protein levels of at least 75% for each of the genes. After subtracting out potential artifacts caused by viral infection, gene set enrichment and causal reasoning engine analysis showed that a significant number of gene expression changes mapped to pathways associated with neurogenesis, long-term potentiation, and synaptic activity.Conclusions: This work demonstrates that despite the complex genetic nature of ASD, there are common molecular mechanisms that connect many of the best established autism candidate genes. By identifying the key regulatory checkpoints in the interlinking transcriptional networks underlying autism, we are better able to discover the ideal points of intervention that provide the broadest efficacy across the diverse population of autism patients.