RNAseq analysis of treatment-dependent signaling changes during inflammation in a mouse cutaneous wound healing model.

BACKGROUND: Despite proven therapeutic effects in inflammatory conditions, the specific mechanisms of phytochemical therapies are not well understood. The transcriptome effects of Traumeel (Tr14), a multicomponent natural product, and diclofenac, a non-selective cyclooxygenase (COX) inhibitor, were compared in a mouse cutaneous wound healing model to identify both known and novel pathways for the anti-inflammatory effect of plant-derived natural products. METHODS: Skin samples from abraded mice were analyzed by single-molecule, amplification-free RNAseq transcript profiling at 7 points between 12 and 192 h after injury. Immediately after injury, the wounds were treated with either diclofenac, Tr14, or placebo control (n = 7 per group/time). RNAseq levels were compared between treatment and control at each time point using a systems biology approach. RESULTS: At early time points (12-36 h), both control and Tr14-treated wounds showed marked increase in the inducible COX2 enzyme mRNA, while diclofenac-treated wounds did not. Tr14, in contrast, modulated lipoxygenase transcripts, especially ALOX12/15, and phospholipases involved in arachidonate metabolism. Notably, Tr14 modulated a group of cell-type specific markers, including the T cell receptor, that could be explained by an overarching effect on the type of cells that were recruited into the wound tissue. CONCLUSIONS: Tr14 and diclofenac had very different effects on the COX/LOX synthetic pathway after cutaneous wounding. Tr14 allowed normal autoinduction of COX2 mRNA, but suppressed mRNA levels for key enzymes in the leukotriene synthetic pathway. Tr14 appeared to have a broad 'phytocellular' effect on the wound transcriptome by altering the balance of cell types present in the wound.

RNA sequencing of blood in coronary artery disease: involvement of regulatory T cell imbalance.

BACKGROUND: Cardiovascular disease had a global prevalence of 523 million cases and 18.6 million deaths in 2019. The current standard for diagnosing coronary artery disease (CAD) is coronary angiography. Surprisingly, despite well-established clinical indications, up to 40% of the one million invasive cardiac catheterizations return a result of 'no blockage'. The present studies employed RNA sequencing of whole blood to identify an RNA signature in patients with angiographically confirmed CAD. METHODS: Whole blood RNA was depleted of ribosomal RNA (rRNA) and analyzed by single-molecule sequencing of RNA (RNAseq) to identify transcripts associated with CAD (TRACs) in a discovery group of 96 patients presenting for elective coronary catheterization. The resulting transcript counts were compared between groups to identify differentially expressed genes (DEGs). RESULTS: Surprisingly, 98% of DEGs/TRACs were down-regulated ~ 1.7-fold in patients with mild to severe CAD (> 20% stenosis). The TRACs were independent of comorbid risk factors for CAD, such as sex, hypertension, and smoking. Bioinformatic analysis identified an enrichment in transcripts such as FoxP1, ICOSLG, IKZF4/Eos, SMYD3, TRIM28, and TCF3/E2A that are likely markers of regulatory T cells (Treg), consistent with known reductions in Tregs in CAD. A validation cohort of 80 patients confirmed the overall pattern (92% down-regulation) and supported many of the Treg-related changes. TRACs were enriched for transcripts associated with stress granules, which sequester RNAs, and ciliary and synaptic transcripts, possibly consistent with changes in the immune synapse of developing T cells. CONCLUSIONS: These studies identify a novel mRNA signature of a Treg-like defect in CAD patients and provides a blueprint for a diagnostic test for CAD. The pattern of changes is consistent with stress-related changes in the maturation of T and Treg cells, possibly due to changes in the immune synapse.

Urinary miRNA Biomarkers of Drug-Induced Kidney Injury and Their Site Specificity Within the Nephron.

Drug-induced kidney injury (DIKI) is a major concern in both drug development and clinical practice. There is an unmet need for biomarkers of glomerular damage and more distal renal injury in the loop of Henle and the collecting duct (CD). A cross-laboratory program to identify and characterize urinary microRNA (miRNA) patterns reflecting tissue- or pathology-specific DIKI was conducted. The overall goal was to propose miRNA biomarker candidates for DIKI that could supplement information provided by protein kidney biomarkers in urine. Rats were treated with nephrotoxicants causing injury to distinct nephron segments: the glomerulus, proximal tubule, thick ascending limb (TAL) of the loop of Henle and CD. Meta-analysis identified miR-192-5p as a potential proximal tubule-specific urinary miRNA candidate. This result was supported by data obtained in laser capture microdissection nephron segments showing that miR-192-5p expression was enriched in the proximal tubule. Discriminative miRNAs including miR-221-3p and -222-3p were increased in urine from rats treated with TAL versus proximal tubule toxicants in accordance with their expression localization in the kidney. Urinary miR-210-3p increased up to 40-fold upon treatment with TAL toxicants and was also enriched in laser capture microdissection samples containing TAL and/or CD versus proximal tubule. miR-23a-3p was enriched in the glomerulus and was increased in urine from rats treated with doxorubicin, a glomerular toxicant, but not with toxicants affecting other nephron segments. Taken together these results suggest that urinary miRNA panels sourced from specific nephron regions may be useful to discriminate the pathology of toxicant-induced lesions in the kidney, thereby contributing to DIKI biomarker development needs for industry, clinical, and regulatory use.

The Utility of Novel Urinary Biomarkers in Mice for Drug Development Studies.

Novel urinary protein biomarkers have recently been identified and qualified in rats for the early detection of renal injury in drug development studies. However, there are few reports on the utility of these renal biomarkers in mice, another important and widely used preclinical animal species for drug development studies. The purpose of this study was to assess the value of these recently qualified biomarkers for the early detection of drug-induced kidney injury (DIKI) in different strains of mice using multiple assay panels. To this end, we evaluated biomarker response to kidney injury induced by several nephrotoxic agents including amphotericin B, compound X, and compound Y. Several of the biomarkers were shown to be sensitive to DIKI in mice. When measured, urinary albumin and neutrophil gelatinase-associated lipocalin were highly sensitive to renal tubular injury, regardless of the assay platforms, mouse strain, and nephrotoxic agents. Depending on the type of renal tubular injury, kidney injury molecule-1 was also highly sensitive, regardless of the assay platforms and mouse strain. Osteopontin and cystatin C were modestly to highly sensitive to renal tubular injury, but the assay type and/or the mouse strain should be considered before using these biomarkers. Calbindin D28 was highly sensitive to injury to the distal nephron in mice. To our knowledge, this is the first report that demonstrates the utility of novel urinary biomarkers evaluated across multiple assay platforms and nephrotoxicants in different mice strains with DIKI. These results will help drug developers make informed decisions when selecting urinary biomarkers for monitoring DIKI in mice for toxicology studies.

Gut microbiome and serum metabolome analyses identify molecular biomarkers and altered glutamate metabolism in fibromyalgia.

BACKGROUND: Fibromyalgia is a complex, relatively unknown disease characterised by chronic, widespread musculoskeletal pain. The gut-brain axis connects the gut microbiome with the brain through the enteric nervous system (ENS); its disruption has been associated with psychiatric and gastrointestinal disorders. To gain an insight into the pathogenesis of fibromyalgia and identify diagnostic biomarkers, we combined different omics techniques to analyse microbiome and serum composition. METHODS: We collected faeces and blood samples to study the microbiome, the serum metabolome and circulating cytokines and miRNAs from a cohort of 105 fibromyalgia patients and 54 age- and environment-matched healthy individuals. We sequenced the V3 and V4 regions of the 16S rDNA gene from faeces samples. UPLC-MS metabolomics and custom multiplex cytokine and miRNA analysis (FirePlex™ technology) were used to examine sera samples. Finally, we combined the different data types to search for potential biomarkers. RESULTS: We found that the diversity of bacteria is reduced in fibromyalgia patients. The abundance of the Bifidobacterium and Eubacterium genera (bacteria participating in the metabolism of neurotransmitters in the host) in these patients was significantly reduced. The serum metabolome analysis revealed altered levels of glutamate and serine, suggesting changes in neurotransmitter metabolism. The combined serum metabolomics and gut microbiome datasets showed a certain degree of correlation, reflecting the effect of the microbiome on metabolic activity. We also examined the microbiome and serum metabolites, cytokines and miRNAs as potential sources of molecular biomarkers of fibromyalgia. CONCLUSIONS: Our results show that the microbiome analysis provides more significant biomarkers than the other techniques employed in the work. Gut microbiome analysis combined with serum metabolomics can shed new light onto the pathogenesis of fibromyalgia. We provide a list of bacteria whose abundance changes in this disease and propose several molecules as potential biomarkers that can be used to evaluate the current diagnostic criteria.

Functional genomic analysis identifies miRNA repertoire regulating C. elegans oocyte development.

Oocyte-specific miRNA function remains unclear in mice and worms because loss of Dgcr8 and Dicer from mouse and worm oocytes, respectively, does not yield oogenic defects. These data lead to several models: (a) miRNAs are not generated in oocytes; (b) miRNAs are generated but do not perform an oogenic function; (c) functional oocyte miRNAs are generated in a manner independent of these enzymes. Here, we test these models using a combination of genomic, expression and functional analyses on the C. elegans germline. We identify a repertoire of at least twenty-three miRNAs that accumulate in four spatial domains in oocytes. Genetic tests demonstrate that oocyte-expressed miRNAs regulate key oogenic processes within their respective expression domains. Unexpectedly, we find that over half of the oocyte-expressed miRNAs are generated through an unknown Drosha independent mechanism. Thus, a functional miRNA repertoire generated via Drosha dependent and independent pathways regulates C. elegans oocyte development.

MicroRNAs Associated With Reverse Left Ventricular Remodeling in Humans Identify Pathways of Heart Failure Progression.

BACKGROUND: Plasma extracellular RNAs have recently garnered interest as biomarkers in heart failure (HF). Most studies in HF focus on single extracellular RNAs related to phenotypes and outcomes, and few describe their functional roles. We hypothesized that clusters of plasma microRNAs (miRNAs) associated with left ventricular (LV) remodeling in human HF would identify novel subsets of genes involved in HF in animal models. METHODS AND RESULTS: We prospectively measured circulating miRNAs in 64 patients with systolic HF (mean age, 64.8 years; 91% men; median LV ejection fraction, 26%) with serial echocardiography (10 months apart) during medical therapy. We defined LV reverse remodeling as a 15% reduction in LV end-systolic volume index. Using principal components analysis, we identified a component associated with LV reverse remodeling (odds ratio=3.99; P=0.01) that provided risk discrimination for LV reverse remodeling superior to a clinical model (C statistic, 0.58 for a clinical model versus 0.71 for RNA-based model). Using network bioinformatics, we uncovered genes not previously widely described in HF regulated simultaneously by >2 miRNAs. We observed increased myocardial expression of these miRNAs during HF development in animals, with downregulation of target gene expression, suggesting coordinate miRNA-mRNA regulation. Target mRNAs were involved in autophagy, metabolism, and inflammation. CONCLUSIONS: Plasma miRNAs associated with LV reverse remodeling in humans are dysregulated in animal HF and target clusters of genes involved in mechanisms implicated in HF. A translational approach integrating human HF, bioinformatics, and model systems may uncover novel pathways involved in HF. CLINICAL TRIAL REGISTRATION: URL: https://www.clinicaltrials.gov. Unique identifier: NCT00351390.

Using FirePlex™ Particle Technology for Multiplex MicroRNA Profiling Without RNA Purification.

Accuracy of miRNA profiling is enhanced when sample processing can be kept to a minimum, avoiding steps such as RNA purification that can introduce bias and inaccuracies. Here we describe a novel multiplex circulating miRNA assay that enables the profiling of up to 65 miRNAs of choice in the same well directly from plasma (including heparin plasma) or serum, with no need for RNA purification. The main component of the assay is FirePlex™ hydrogel particles, which enable the multiplex capture of miRNAs with picomolar sensitivity and high specificity. Results are obtained using conventional flow cytometry and easy to use software, which allows fast analysis and interpretation of the experimental data. This chapter provides methods to profile miRNAs with PCR sensitivity from as little as 10 µL of crude biofluid sample, or from less than 100 pg of purified RNA.

Deep Sequencing Transcriptome Analysis of Murine Wound Healing: Effects of a Multicomponent, Multitarget Natural Product Therapy-Tr14.

Wound healing involves an orchestrated response that engages multiple processes, such as hemostasis, cellular migration, extracellular matrix synthesis, and in particular, inflammation. Using a murine model of cutaneous wound repair, the transcriptome was mapped from 12 h to 8 days post-injury, and in response to a multicomponent, multi-target natural product, Tr14. Using single-molecule RNA sequencing (RNA-seq), there were clear temporal changes in known transcripts related to wound healing pathways, and additional novel transcripts of both coding and non-coding genes. Tr14 treatment modulated >100 transcripts related to key wound repair pathways, such as response to wounding, wound contraction, and cytokine response. The results provide the most precise and comprehensive characterization to date of the transcriptome's response to skin damage, repair, and multicomponent natural product therapy. By understanding the wound repair process, and the effects of natural products, it should be possible to intervene more effectively in diseases involving aberrant repair.

Comparative miRNA Analysis of Urine Extracellular Vesicles Isolated through Five Different Methods.

Urine extracellular vesicles are a valuable low-invasive source of information, especially for the cells of the genitourinary tract. In the search for biomarkers, different techniques have been developed to isolate and characterize the cargo of these vesicles. In the present work, we compare five of these different isolation methods (three commercial isolation kits, ultracentrifugation, and lectin-based purification) and perform miRNA profiling using a multiplex miRNA assay. The results showed high correlation through all isolation techniques, and 48 out of 68 miRNAs were detected above the detection limit at least 10 times. The results obtained by multiplex assay were validated through Taqman qPCR. In addition, using this technique combined with a clinically friendly extracellular vesicle (uEV)-enrichment method, we performed the analysis of selected miRNAs in urine from patients affected with bladder cancer, benign prostate hyperplasia, or prostate cancer. Importantly, we found that those miRNAs could be detected in almost 100% of the samples, and no significant differences were observed between groups. Our results support the feasibility of analyzing exosomes-associated miRNAs using a methodology that requires a small volume of urine and is compatible with a clinical environment and high-throughput analysis.

RNA editing regulates transposon-mediated heterochromatic gene silencing.

Heterochromatin formation drives epigenetic mechanisms associated with silenced gene expression. Repressive heterochromatin is established through the RNA interference pathway, triggered by double-stranded RNAs (dsRNAs) that can be modified via RNA editing. However, the biological consequences of such modifications remain enigmatic. Here we show that RNA editing regulates heterochromatic gene silencing in Drosophila. We utilize the binding activity of an RNA-editing enzyme to visualize the in vivo production of a long dsRNA trigger mediated by Hoppel transposable elements. Using homologous recombination, we delete this trigger, dramatically altering heterochromatic gene silencing and chromatin architecture. Furthermore, we show that the trigger RNA is edited and that dADAR serves as a key regulator of chromatin state. Additionally, dADAR auto-editing generates a natural suppressor of gene silencing. Lastly, systemic differences in RNA editing activity generates interindividual variation in silencing state within a population. Our data reveal a global role for RNA editing in regulating gene expression.

Genome-wide analysis of A-to-I RNA editing by single-molecule sequencing in Drosophila.

The accurate and thorough genome-wide detection of adenosine-to-inosine editing, a biologically indispensable process, has proven challenging. Here, we present a discovery pipeline in adult Drosophila, with 3,581 high-confidence editing sites identified with an estimated accuracy of 87%. The target genes and specific sites highlight global biological properties and functions of RNA editing, including hitherto-unknown editing in well-characterized classes of noncoding RNAs and 645 sites that cause amino acid substitutions, usually at conserved positions. The spectrum of functions that these gene targets encompass suggests that editing participates in a diverse set of cellular processes. Editing sites in Drosophila exhibit sequence-motif preferences and tend to be concentrated within a small subset of total RNAs. Finally, editing regulates expression levels of target mRNAs and strongly correlates with alternative splicing.

Assessing the validity and reproducibility of genome-scale predictions.

MOTIVATION: Validation and reproducibility of results is a central and pressing issue in genomics. Several recent embarrassing incidents involving the irreproducibility of high-profile studies have illustrated the importance of this issue and the need for rigorous methods for the assessment of reproducibility. RESULTS: Here, we describe an existing statistical model that is very well suited to this problem. We explain its utility for assessing the reproducibility of validation experiments, and apply it to a genome-scale study of adenosine deaminase acting on RNA (ADAR)-mediated RNA editing in Drosophila. We also introduce a statistical method for planning validation experiments that will obtain the tightest reproducibility confidence limits, which, for a fixed total number of experiments, returns the optimal number of replicates for the study. AVAILABILITY: Downloadable software and a web service for both the analysis of data from a reproducibility study and for the optimal design of these studies is provided at http://ccmbweb.ccv.brown.edu/reproducibility.html .

VlincRNAs controlled by retroviral elements are a hallmark of pluripotency and cancer.

BACKGROUND: The function of the non-coding portion of the human genome remains one of the most important questions of our time. Its vast complexity is exemplified by the recent identification of an unusual and notable component of the transcriptome - very long intergenic non-coding RNAs, termed vlincRNAs. RESULTS: Here we identify 2,147 vlincRNAs covering 10 percent of our genome. We show they are present not only in cancerous cells, but also in primary cells and normal human tissues, and are controlled by canonical promoters. Furthermore, vlincRNA promoters frequently originate from within endogenous retroviral sequences. Strikingly, the number of vlincRNAs expressed from endogenous retroviral promoters strongly correlates with pluripotency or the degree of malignant transformation. These results suggest a previously unknown connection between the pluripotent state and cancer via retroviral repeat-driven expression of vlincRNAs. Finally, we show that vlincRNAs can be syntenically conserved in humans and mouse and their depletion using RNAi can cause apoptosis in cancerous cells. CONCLUSIONS: These intriguing observations suggest that vlincRNAs could create a framework that combines many existing short ESTs and lincRNAs into a landscape of very long transcripts functioning in the regulation of gene expression in the nucleus. Certain types of vlincRNAs participate at specific stages of normal development and, based on analysis of a limited set of cancerous and primary cell lines, they appear to be co-opted by cancer-associated transcriptional programs. This provides additional understanding of transcriptome regulation during the malignant state, and could lead to additional targets and options for its reversal.

On the importance of small changes in RNA expression.

The analysis of the differential expression of genes has been the key goal of many molecular biology methods for decades and will remain with us for decades to come. It constitutes a fundamental resource at our disposal for determining the relationship between products of transcription, biology and disease. The completed genome sequencing of many common species allowed microarrays and RNA sequencing (RNAseq) to become major tools in Systems Biology. However, we estimate that at least half of all experiments ignore transcripts that change less than some subjectively chosen threshold, typically around 2-3 fold. Here we show that a majority of the informative RNAs and differentially expressed transcripts can exhibit fold changes less than 2. We use highly quantitative single-molecule sequencing of total cellular RNA derived from a time course of inflammatory response, a process critical to a large number of diseases. Furthermore, we show that enrichment of biologically-relevant functions occurs even at very low fold changes in RNA levels. In addition, we show that most of the common statistical methods can reliably detect transcripts with low fold change when as few as 3 biological replicates are sequenced using single-molecule based RNAseq. In conclusion, given the prevalence of expression profiling in current research, the loss of data in half of all expression studies results in a significant, yet needless drain on the discovery process.

Intronic RNAs constitute the major fraction of the non-coding RNA in mammalian cells.

BACKGROUND: The function of RNA from the non-coding (the so called "dark matter") regions of the genome has been a subject of considerable recent debate. Perhaps the most controversy is regarding the function of RNAs found in introns of annotated transcripts, where most of the reads that map outside of exons are usually found. However, it has been reported that the levels of RNA in introns are minor relative to those of the corresponding exons, and that changes in the levels of intronic RNAs correlate tightly with that of adjacent exons. This would suggest that RNAs produced from the vast expanse of intronic space are just pieces of pre-mRNAs or excised introns en route to degradation. RESULTS: We present data that challenges the notion that intronic RNAs are mere by-standers in the cell. By performing a highly quantitative RNAseq analysis of transcriptome changes during an inflammation time course, we show that intronic RNAs have a number of features that would be expected from functional, standalone RNA species. We show that there are thousands of introns in the mouse genome that generate RNAs whose overall abundance, which changes throughout the inflammation timecourse, and other properties suggest that they function in yet unknown ways. CONCLUSIONS: So far, the focus of non-coding RNA discovery has shied away from intronic regions as those were believed to simply encode parts of pre-mRNAs. Results presented here suggest a very different situation--the sequences encoded in the introns appear to harbor a yet unexplored reservoir of novel, functional RNAs. As such, they should not be ignored in surveys of functional transcripts or other genomic studies.

Mice cloned from olfactory sensory neurons.

Cloning by nuclear transplantation has been successfully carried out in various mammals, including mice. Until now mice have not been cloned from post-mitotic cells such as neurons. Here, we have generated fertile mouse clones derived by transferring the nuclei of post-mitotic, olfactory sensory neurons into oocytes. These results indicate that the genome of a post-mitotic, terminally differentiated neuron can re-enter the cell cycle and be reprogrammed to a state of totipotency after nuclear transfer. Moreover, the pattern of odorant receptor gene expression and the organization of odorant receptor genes in cloned mice was indistinguishable from wild-type animals, indicating that irreversible changes to the DNA of olfactory neurons do not accompany receptor gene choice.

Coordination of the random asynchronous replication of autosomal loci.

Random monoallelic expression and asynchronous replication define an unusual class of autosomal mammalian genes. We show that every cell has randomly chosen either the maternal or paternal copy of each given autosome pair, such that alleles of these genes scattered across the chosen chromosome replicate earlier than the alleles on the homologous chromosome. Thus, chromosome-pair non-equivalence, rather than being limited to X-chromosome inactivation, is a fundamental property of mouse chromosomes.