Barcoded cDNA library preparation for small RNA profiling by next-generation sequencing

Epigenetic repression of antiviral genes by SARS-CoV-2 NSP1

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) evades the innate immune machinery through multiple viral proteins, including nonstructural protein 1 (NSP1). While NSP1 is known to suppress translation of host mRNAs, the mechanisms underlying its immune evasion properties remain elusive. By integrating RNA-seq, ribosome footprinting, and ChIP-seq in A549 cells we found that NSP1 predominantly represses transcription of immune-related genes by favoring Histone 3 Lysine 9 dimethylation (H3K9me2). G9a/GLP H3K9 methyltransferase inhibitor UNC0638 restored expression of antiviral genes and restricted SARS-CoV-2 replication. Our multi-omics study unravels an epigenetic mechanism underlying host immune evasion by SARS-CoV-2 NSP1. Elucidating the factors involved in this phenomenon, may have implications for understanding and treating viral infections and other immunomodulatory diseases.

Death Induced by Survival gene Elimination (DISE) correlates with neurotoxicity in Alzheimer's disease and aging

Alzheimer's disease (AD) is characterized by progressive neurodegeneration, but the specific events that cause cell death remain poorly understood. Death Induced by Survival gene Elimination (DISE) is a cell death mechanism mediated by short (s) RNAs acting through the RNA-induced silencing complex (RISC). DISE is thus a form of RNA interference, in which G-rich 6mer seed sequences in the sRNAs (position 2-7) target hundreds of C-rich 6mer seed matches in genes essential for cell survival, resulting in the activation of cell death pathways. Here, using Argonaute precipitation and RNAseq (Ago-RP-Seq), we analyze RISC-bound sRNAs to quantify 6mer seed toxicity in several model systems. In mouse AD models and aging brain, in induced pluripotent stem cell-derived neurons from AD patients, and in cells exposed to Aβ42 oligomers, RISC-bound sRNAs show a shift to more toxic 6mer seeds compared to controls. In contrast, in brains of "SuperAgers", humans over age 80 who have superior memory performance, RISC-bound sRNAs are shifted to more nontoxic 6mer seeds. Cells depleted of nontoxic sRNAs are sensitized to Aβ42-induced cell death, and reintroducing nontoxic RNAs is protective. Altogether, the correlation between DISE and Aβ42 toxicity suggests that increasing the levels of nontoxic miRNAs in the brain or blocking the activity of toxic RISC-bound sRNAs could ameliorate neurodegeneration.

Analysis of the Contribution of 6-mer Seed Toxicity to HIV-1-Induced Cytopathicity

HIV-1 (HIV) infects CD4+ T cells, the gradual depletion of which can lead to AIDS in the absence of antiretroviral therapy (ART). Some cells, however, survive HIV infection and persist as part of the latently infected reservoir that causes recurrent viremia after ART cessation. Improved understanding of the mechanisms of HIV-mediated cell death could lead to a way to clear the latent reservoir. Death induced by survival gene elimination (DISE), an RNA interference (RNAi)-based mechanism, kills cells through short RNAs (sRNAs) with toxic 6-mer seeds (positions 2 to 7 of sRNA). These toxic seeds target the 3' untranslated region (UTR) of mRNAs, decreasing the expression of hundreds of genes critical for cell survival. In most cells under normal conditions, highly expressed cell-encoded nontoxic microRNAs (miRNAs) block access of toxic sRNAs to the RNA-induced silencing complex (RISC) that mediates RNAi, promoting cell survival. HIV has been shown to inhibit the biogenesis of host miRNAs in multiple ways. We now report that HIV infection of cells deficient in miRNA expression or function results in enhanced RISC loading of an HIV-encoded miRNA HIV-miR-TAR-3p, which can kill cells by DISE through a noncanonical (positions 3 to 8) 6-mer seed. In addition, cellular RISC-bound sRNAs shift to lower seed viability. This also occurs after latent HIV provirus reactivation in J-Lat cells, suggesting independence of permissiveness of cells to viral infection. More precise targeting of the balance between protective and cytotoxic sRNAs could provide new avenues to explore novel cell death mechanisms that could be used to kill latent HIV. IMPORTANCE Several mechanisms by which initial HIV infection is cytotoxic to infected cells have been reported and involve various forms of cell death. Characterizing the mechanisms underlying the long-term survival of certain T cells that become persistent provirus reservoirs is critical to developing a cure. We recently discovered death induced by survival gene elimination (DISE), an RNAi-based mechanism of cell death whereby toxic short RNAs (sRNAs) containing 6-mer seed sequences (exerting 6-mer seed toxicity) targeting essential survival genes are loaded into RNA-induced silencing complex (RISC) complexes, resulting in inescapable cell death. We now report that HIV infection in cells with low miRNA expression causes a shift of mostly cellular RISC-bound sRNAs to more toxic seeds. This could prime cells to DISE and is further enhanced by the viral microRNA (miRNA) HIV-miR-TAR-3p, which carries a toxic noncanonical 6-mer seed. Our data provide multiple new avenues to explore novel cell death mechanisms that could be used to kill latent HIV.

CD95/Fas ligand mRNA is toxic to cells through more than one mechanism

CD95/Fas ligand (CD95L) induces apoptosis through protein binding to the CD95 receptor. However, CD95L mRNA also induces toxicity in the absence of CD95 through induction of DISE (Death Induced by Survival Gene Elimination), a form of cell death mediated by RNA interference (RNAi). We now report that CD95L mRNA processing generates a short (s)RNA nearly identical to shL3, a commercial CD95L-targeting shRNA that led to the discovery of DISE. Neither of the miRNA biogenesis proteins Drosha nor Dicer are required for this processing. Interestingly, CD95L toxicity depends on the core component of the RISC, Ago2, in some cell lines, but not in others. In the HCT116 colon cancer cell line, Ago 1-4 appear to function redundantly in RNAi. In fact, Ago 1/2/3 knockout cells retain sensitivity to CD95L mRNA toxicity. Toxicity was only blocked by mutation of all in-frame start codons in the CD95L ORF. Dying cells exhibited an enrichment of RISC bound (R)-sRNAs with toxic 6mer seed sequences, while expression of the non-toxic CD95L mutant enriched for loading of R-sRNAs with nontoxic 6mer seeds. However, CD95L is not the only source of these R-sRNAs. We find that CD95L mRNA may induce DISE directly and indirectly, and that alternate mechanisms may underlie CD95L mRNA processing and toxicity.

Enhanced transcriptome-wide RNA G-quadruplex sequencing for low RNA input samples with rG4-seq 2.0

Background

RNA G-quadruplexes (rG4s) are non-canonical structural motifs that have diverse functional and regulatory roles, for instance in transcription termination, alternative splicing, mRNA localization and stabilization, and translational process. We recently developed the RNA G-quadruplex structure sequencing (rG4-seq) technique and described rG4s in both eukaryotic and prokaryotic transcriptomes. However, rG4-seq suffers from a complicated gel purification step and limited PCR product yield, thus requiring a high amount of RNA input, which limits its applicability in more physiologically or clinically relevant studies often characterized by the limited availability of biological material and low RNA abundance. Here, we redesign and enhance the workflow of rG4-seq to address this issue.

Results

We developed rG4-seq 2.0 by introducing a new ssDNA adapter containing deoxyuridine during library preparation to enhance library quality with no gel purification step, less PCR amplification cycles and higher yield of PCR products. We demonstrate that rG4-seq 2.0 produces high-quality cDNA libraries that support reliable and reproducible rG4 identification at varying RNA inputs, including RNA mounts as low as 10 ng. rG4-seq 2.0 also improved the rG4-seq calling outcome and nucleotide bias in rG4 detection persistent in rG4-seq 1.0. We further provide in vitro mapping of rG4 in the HEK293T cell line, and recommendations for assessing RNA input and sequencing depth for individual rG4 studies based on transcript abundance.

Conclusions

rG4-seq 2.0 can improve the identification and study of rG4s in low abundance transcripts, and our findings can provide insights to optimize cDNA library preparation in other related methods.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01448-3.

Functional editing of endogenous genes through rapid selection of cell pools (Rapid generation of endogenously tagged genes in Drosophila ovarian somatic sheath cells)

The combination of genome-editing and epitope tagging provides a powerful strategy to study proteins with high affinity and specificity while preserving their physiological expression patterns. However, stably modifying endogenous genes in cells that do not allow for clonal selection has been challenging. Here, we present a simple and fast strategy to generate stable, endogenously tagged alleles in a non-transformed cell culture model. At the example of piwi in Drosophila ovarian somatic sheath cells, we show that this strategy enables the generation of an N-terminally tagged protein that emulates the expression level and subcellular localization of the wild type protein and forms functional Piwi-piRNA complexes. We further present a concise workflow to establish endogenously N-terminally and C-terminally tagged proteins, and knockout alleles through rapid selection of cell pools in fly and human models.

Novel MicroRNA-Regulated Transcript Networks Are Associated with Chemotherapy Response in Ovarian Cancer

High-grade serous ovarian cancer (HGSOC) is a highly lethal gynecologic cancer, in part due to resistance to platinum-based chemotherapy reported among 20% of patients. This study aims to generate novel hypotheses of the biological mechanisms underlying chemotherapy resistance, which remain poorly understood. Differential expression analyses of mRNA- and microRNA-sequencing data from HGSOC patients of The Cancer Genome Atlas identified 21 microRNAs associated with angiogenesis and 196 mRNAs enriched for adaptive immunity and translation. Coexpression network analysis identified three microRNA networks associated with chemotherapy response enriched for lipoprotein transport and oncogenic pathways, as well as two mRNA networks enriched for ubiquitination and lipid metabolism. These network modules were replicated in two independent ovarian cancer cohorts. Moreover, integrative analyses of the mRNA/microRNA sequencing and single-nucleotide polymorphisms (SNPs) revealed potential regulation of significant mRNA transcripts by microRNAs and SNPs (expression quantitative trait loci). Thus, we report novel transcriptional networks and biological pathways associated with resistance to platinum-based chemotherapy in HGSOC patients. These results expand our understanding of the effector networks and regulators of chemotherapy response, which will help to improve the management of ovarian cancer.

Serum MicroRNA Transcriptomics and Acute Rejection or Recurrent Hepatitis C Virus in Human Liver Allograft Recipients: A Pilot Study

BACKGROUND

Acute rejection (AR) and recurrent hepatitis C virus (R-HCV) are significant complications in liver allograft recipients. Noninvasive diagnosis of intragraft pathologies may improve their management.

METHODS

We performed small RNA sequencing and microRNA (miRNA) microarray profiling of RNA from sera matched to liver allograft biopsies from patients with nonimmune, nonviral (NINV) native liver disease. Absolute levels of informative miRNAs in 91 sera matched to 91 liver allograft biopsies were quantified using customized real-time quantitative PCR (RT-qPCR) assays: 30 biopsy-matched sera from 26 unique NINV patients and 61 biopsy-matched sera from 41 unique R-HCV patients. The association between biopsy diagnosis and miRNA abundance was analyzed by logistic regression and calculating the area under the receiver operating characteristic curve.

RESULTS

Nine miRNAs-miR-22, miR-34a, miR-122, miR-148a, miR-192, miR-193b, miR-194, miR-210, and miR-885-5p-were identified by both sRNA-seq and TLDA to be associated with NINV-AR. Logistic regression analysis of absolute levels of miRNAs and goodness-of-fit of predictors identified a linear combination of miR-34a + miR-210 (P < 0.0001) as the best statistical model and miR-122 + miR-210 (P < 0.0001) as the best model that included miR-122. A different linear combination of miR-34a + miR-210 (P < 0.0001) was the best model for discriminating NINV-AR from R-HCV with intragraft inflammation, and miR-34a + miR-122 (P < 0.0001) was the best model for discriminating NINV-AR from R-HCV with intragraft fibrosis.

CONCLUSIONS

Circulating levels of miRNAs, quantified using customized RT-qPCR assays, may offer a rapid and noninvasive means of diagnosing AR in human liver allografts and for discriminating AR from intragraft inflammation or fibrosis due to R-HCV.

Topology preserving stratification of tissue neoplasticity using Deep Neural Maps and microRNA signatures

Background

Accurate cancer classification is essential for correct treatment selection and better prognostication. microRNAs (miRNAs) are small RNA molecules that negatively regulate gene expression, and their dyresgulation is a common disease mechanism in many cancers. Through a clearer understanding of miRNA dysregulation in cancer, improved mechanistic knowledge and better treatments can be sought.

Results

We present a topology-preserving deep learning framework to study miRNA dysregulation in cancer. Our study comprises miRNA expression profiles from 3685 cancer and non-cancer tissue samples and hierarchical annotations on organ and neoplasticity status. Using unsupervised learning, a two-dimensional topological map is trained to cluster similar tissue samples. Labelled samples are used after training to identify clustering accuracy in terms of tissue-of-origin and neoplasticity status. In addition, an approach using activation gradients is developed to determine the attention of the networks to miRNAs that drive the clustering. Using this deep learning framework, we classify the neoplasticity status of held-out test samples with an accuracy of 91.07%, the tissue-of-origin with 86.36%, and combined neoplasticity status and tissue-of-origin with an accuracy of 84.28%. The topological maps display the ability of miRNAs to recognize tissue types and neoplasticity status. Importantly, when our approach identifies samples that do not cluster well with their respective classes, activation gradients provide further insight in cancer subtypes or grades.

Conclusions

An unsupervised deep learning approach is developed for cancer classification and interpretation. This work provides an intuitive approach for understanding molecular properties of cancer and has significant potential for cancer classification and treatment selection.

RNA Externalized by Neutrophil Extracellular Traps Promotes Inflammatory Pathways in Endothelial Cells

OBJECTIVE

Neutrophil extracellular traps (NETs) are extracellular lattices composed of nucleic material bound to neutrophil granule proteins. NETs may play pathogenic roles in the development and severity of autoimmune diseases such as systemic lupus erythematosus (SLE), at least in part, through induction of type I interferon (IFN) responses via externalization of oxidized immunostimulatory DNA. A distinct subset of SLE proinflammatory neutrophils (low-density granulocytes [LDGs]) displays enhanced ability to form proinflammatory NETs that damage the vasculature. We undertook this study to assess whether NET-bound RNA can contribute to inflammatory responses in endothelial cells (ECs) and the pathways that mediate this effect.

METHODS

Expression of newly synthesized and total RNA was quantified in NETs from healthy controls and lupus patients. The ability of ECs to take up NET-bound RNA and downstream induction of type I IFN responses were quantified. RNAs present in NETs were sequenced and specific small RNAs were tested for induction of endothelial type I IFN pathways.

RESULTS

NETs extruded RNA that was internalized by ECs, and this was enhanced when NET-bound nucleic acids were oxidized, particularly in lupus LDG-derived NETs. Internalization of NET-bound RNA by ECs was dependent on endosomal Toll-like receptors (TLRs) and the actin cytoskeleton and induced type I IFN-stimulated genes (ISGs). This ISG induction was dependent on NET-associated microRNA let-7b, a small RNA expressed at higher levels in LDG-derived NETs, which acted as a TLR-7 agonist.

CONCLUSION

These findings highlight underappreciated roles for small RNAs externalized in NETs in the induction of proinflammatory responses in vascular cells, with implications for lupus vasculopathy.

Cellular abundance shapes function in piRNA-guided genome defense

Defense against genome invaders universally relies on RNA-guided immunity. Prokaryotic CRISPR-Cas and eukaryotic RNA interference pathways recognize targets by complementary base-pairing, which places the sequences of their guide RNAs at the center of self/nonself discrimination. Here, we explore the sequence space of PIWI-interacting RNAs (piRNAs), the genome defense of animals, and establish functional priority among individual sequences. Our results reveal that only the topmost abundant piRNAs are commonly present in every cell, whereas rare sequences generate cell-to-cell diversity in flies and mice. We identify a skewed distribution of sequence abundance as a hallmark of piRNA populations and show that quantitative differences of more than a 1000-fold are established by conserved mechanisms of biogenesis. Finally, our genomics analyses and direct reporter assays reveal that abundance determines function in piRNA-guided genome defense. Taken together, we identify an effective sequence space and untangle two classes of piRNAs that differ in complexity and function. The first class represents the topmost abundant sequences and drives silencing of genomic parasites. The second class sparsely covers an enormous sequence space. These rare piRNAs cannot function in every cell, every individual, or every generation but create diversity with potential for adaptation in the ongoing arms race with genome invaders.

PAR-CLIP: A Method for Transcriptome-Wide Identification of RNA Binding Protein Interaction Sites

During post-transcriptional gene regulation (PTGR), RNA binding proteins (RBPs) interact with all classes of RNA to control RNA maturation, stability, transport, and translation. Here, we describe Photoactivatable-Ribonucleoside-Enhanced Crosslinking and Immunoprecipitation (PAR-CLIP), a transcriptome-scale method for identifying RBP binding sites on target RNAs with nucleotide-level resolution. This method is readily applicable to any protein directly contacting RNA, including RBPs that are predicted to bind in a sequence- or structure-dependent manner at discrete RNA recognition elements (RREs), and those that are thought to bind transiently, such as RNA polymerases or helicases.

Molecular markers of risk of subsequent invasive breast cancer in women with ductal carcinoma in situ: protocol for a population-based cohort study

INTRODUCTION

Ductal carcinoma in situ (DCIS) of the breast is a non-obligate precursor of invasive breast cancer (IBC). Many DCIS patients are either undertreated or overtreated. The overarching goal of the study described here is to facilitate detection of patients with DCIS at risk of IBC development. Here, we propose to use risk factor data and formalin-fixed paraffin-embedded (FFPE) DCIS tissue from a large, ethnically diverse, population-based cohort of 8175 women with a first diagnosis of DCIS and followed for subsequent IBC to: identify/validate miRNA expression changes in DCIS tissue associated with risk of subsequent IBC; evaluate ipsilateral IBC risk in association with two previously identified marker sets (triple immunopositivity for p16, COX-2, Ki67; Oncotype DX Breast DCIS score); examine the association of risk factor data with IBC risk.

METHODS AND ANALYSIS

We are conducting a series of case-control studies nested within the cohort. Cases are women with DCIS who developed subsequent IBC; controls (2/case) are matched to cases on calendar year of and age at DCIS diagnosis. We project 485 cases/970 controls in the aim focused on risk factors. We estimate obtaining FFPE tissue for 320 cases/640 controls for the aim focused on miRNAs; of these, 173 cases/346 controls will be included in the aim focused on p16, COX-2 and Ki67 immunopositivity, and of the latter, 156 case-control pairs will be included in the aim focused on the Oncotype DX Breast DCIS score®. Multivariate conditional logistic regression will be used for statistical analyses.

ETHICS AND DISSEMINATION

Ethics approval was obtained from the Institutional Review Boards of Albert Einstein College of Medicine (IRB 2014-3611), Kaiser Permanente Colorado, Kaiser Permanente Hawaii, Henry Ford Health System, Mayo Clinic, Marshfield Clinic Research Institute and Hackensack Meridian Health, and from Lifespan Research Protection Office. The study results will be presented at meetings and published in peer-reviewed journals.

The Ratio of Toxic-to-Nontoxic miRNAs Predicts Platinum Sensitivity in Ovarian Cancer

Ovarian cancer remains one of the deadliest gynecologic malignancies affecting women, and development of resistance to platinum remains a major barrier to achieving a cure. Multiple mechanisms have been identified to confer platinum resistance. Numerous miRNAs have been linked to platinum sensitivity and resistance in ovarian cancer. miRNA activity occurs mainly when the guide strand of the miRNA, with its seed sequence at position 2-7/8, is loaded into the RNA-induced silencing complex (RISC) and targets complementary short seed matches in the 3' untranslated region of mRNAs. Toxic 6mer seeds, which target genes critical for cancer cell survival, have been found in tumor-suppressive miRNAs. Many siRNAs and short hairpin RNAs (shRNA) can also kill cancer cells via toxic seeds, the most toxic of which carry G-rich 6mer seed sequences. We showed here that treatment of ovarian cancer cells with platinum led to increased RISC-bound miRNAs carrying toxic 6mer seeds and decreased miRNAs with nontoxic seeds. Platinum-tolerant cells did not exhibit this toxicity shift but retained sensitivity to cell death mediated by siRNAs carrying toxic 6mer seeds. Analysis of RISC-bound miRNAs in tumors from patients with ovarian cancer revealed that the ratio between miRNAs with toxic versus nontoxic seeds was predictive of treatment outcome. Application of the 6mer seed toxicity concept to cancer relevant miRNAs provides a new framework for understanding and predicting cancer therapy responses. SIGNIFICANCE: These findings demonstrate that the balance of miRNAs that carry toxic and nontoxic 6mer seeds contributes to platinum resistance in ovarian cancer.

Longitudinal profiling of circulating miRNA during cardiac allograft rejection: a proof-of-concept study

AIMS

Allograft rejection following heart transplantation (HTx) is a serious complication even in the era of modern immunosuppressive regimens and causes up to a third of early deaths after HTx. Allograft rejection is mediated by a cascade of immune mechanisms leading to acute cellular rejection (ACR) and/or antibody-mediated rejection (AMR). The gold standard for monitoring allograft rejection is invasive endomyocardial biopsy that exposes patients to complications. Little is known about the potential of circulating miRNAs as biomarkers to detect cardiac allograft rejection. We here present a systematic analysis of circulating miRNAs as biomarkers and predictors for allograft rejection after HTx using next-generation small RNA sequencing.

METHODS AND RESULTS

We used next-generation small RNA sequencing to investigate circulating miRNAs among HTx recipients (10 healthy controls, 10 heart failure patients, 13 ACR, and 10 AMR). MiRNA profiling was performed at different time points before, during, and after resolution of the rejection episode. We found three miRNAs with significantly increased serum levels in patients with biopsy-proven cardiac rejection when compared with patients without rejection: hsa-miR-139-5p, hsa-miR-151a-5p, and hsa-miR-186-5p. We identified miRNAs that may serve as potential predictors for the subsequent development of ACR: hsa-miR-29c-3p (ACR) and hsa-miR-486-5p (AMR). Overall, hsa-miR-486-5p was most strongly associated with acute rejection episodes.

CONCLUSIONS

Monitoring cardiac allograft rejection using circulating miRNAs might represent an alternative strategy to invasive endomyocardial biopsy.

A miR-375/YAP axis regulates neuroendocrine differentiation and tumorigenesis in lung carcinoid cells

Lung carcinoids are variably aggressive and mechanistically understudied neuroendocrine neoplasms (NENs). Here, we identified and elucidated the function of a miR-375/yes-associated protein (YAP) axis in lung carcinoid (H727) cells. miR-375 and YAP are respectively high and low expressed in wild-type H727 cells. Following lentiviral CRISPR/Cas9-mediated miR-375 depletion, we identified distinct transcriptomic changes including dramatic YAP upregulation. We also observed a significant decrease in neuroendocrine differentiation and substantial reductions in cell proliferation, transformation, and tumor growth in cell culture and xenograft mouse disease models. Similarly, YAP overexpression resulted in distinct and partially overlapping transcriptomic changes, phenocopying the effects of miR-375 depletion in the same models as above. Transient YAP knockdown in miR-375-depleted cells reversed the effects of miR-375 on neuroendocrine differentiation and cell proliferation. Pathways analysis and confirmatory real-time PCR studies of shared dysregulated target genes indicate that this axis controls neuroendocrine related functions such as neural differentiation, exocytosis, and secretion. Taken together, we provide compelling evidence that a miR-375/YAP axis is a critical mediator of neuroendocrine differentiation and tumorigenesis in lung carcinoid cells.

Identification and Quantification of Small RNAs

RNA silencing plays a critical role in diverse biological processes in plants including growth, development, and responses to abiotic and biotic stresses. RNA silencing is guided by small non-coding RNAs (sRNAs) with the length of 21-24 nucleotides (nt) that are loaded into Argonaute (AGO) to repress expression of target loci and transcripts through transcriptional or posttranscriptional gene silencing mechanisms. Identification and quantitative characterization of sRNAs are crucial steps toward appreciation of their functions in biology. Here, we developed a step-by-step protocol to precisely illustrate the process of cloning of sRNA libraries and correspondingly computational analysis of the recovered sRNAs. This protocol can be used in all kinds of organisms, including Arabidopsis, and is compatible with various high-throughput sequence technologies such as Illumina Hiseq. Thus, we wish that this protocol represents an accurate way to identify and quantify sRNAs in vivo.

Classifying Lung Neuroendocrine Neoplasms through MicroRNA Sequence Data Mining

Lung neuroendocrine neoplasms (NENs) can be challenging to classify due to subtle histologic differences between pathological types. MicroRNAs (miRNAs) are small RNA molecules that are valuable markers in many neoplastic diseases. To evaluate miRNAs as classificatory markers for lung NENs, we generated comprehensive miRNA expression profiles from 14 typical carcinoid (TC), 15 atypical carcinoid (AC), 11 small cell lung carcinoma (SCLC), and 15 large cell neuroendocrine carcinoma (LCNEC) samples, through barcoded small RNA sequencing. Following sequence annotation and data preprocessing, we randomly assigned these profiles to discovery and validation sets. Through high expression analyses, we found that miR-21 and -375 are abundant in all lung NENs, and that miR-21/miR-375 expression ratios are significantly lower in carcinoids (TC and AC) than in neuroendocrine carcinomas (NECs; SCLC and LCNEC). Subsequently, we ranked and selected miRNAs for use in miRNA-based classification, to discriminate carcinoids from NECs. Using miR-18a and -155 expression, our classifier discriminated these groups in discovery and validation sets, with 93% and 100% accuracy. We also identified miR-17, -103, and -127, and miR-301a, -106b, and -25, as candidate markers for discriminating TC from AC, and SCLC from LCNEC, respectively. However, these promising findings require external validation due to sample size.

Characterizing and classifying neuroendocrine neoplasms through microRNA sequencing and data mining

Neuroendocrine neoplasms (NENs) are clinically diverse and incompletely characterized cancers that are challenging to classify. MicroRNAs (miRNAs) are small regulatory RNAs that can be used to classify cancers. Recently, a morphology-based classification framework for evaluating NENs from different anatomical sites was proposed by experts, with the requirement of improved molecular data integration. Here, we compiled 378 miRNA expression profiles to examine NEN classification through comprehensive miRNA profiling and data mining. Following data preprocessing, our final study cohort included 221 NEN and 114 non-NEN samples, representing 15 NEN pathological types and 5 site-matched non-NEN control groups. Unsupervised hierarchical clustering of miRNA expression profiles clearly separated NENs from non-NENs. Comparative analyses showed that miR-375 and miR-7 expression is substantially higher in NEN cases than non-NEN controls. Correlation analyses showed that NENs from diverse anatomical sites have convergent miRNA expression programs, likely reflecting morphological and functional similarities. Using machine learning approaches, we identified 17 miRNAs to discriminate 15 NEN pathological types and subsequently constructed a multilayer classifier, correctly identifying 217 (98%) of 221 samples and overturning one histological diagnosis. Through our research, we have identified common and type-specific miRNA tissue markers and constructed an accurate miRNA-based classifier, advancing our understanding of NEN diversity.

The Extracellular RNA Communication Consortium: Establishing Foundational Knowledge and Technologies for Extracellular RNA Research

The Extracellular RNA Communication Consortium (ERCC) was launched to accelerate progress in the new field of extracellular RNA (exRNA) biology and to establish whether exRNAs and their carriers, including extracellular vesicles (EVs), can mediate intercellular communication and be utilized for clinical applications. Phase 1 of the ERCC focused on exRNA/EV biogenesis and function, discovery of exRNA biomarkers, development of exRNA/EV-based therapeutics, and construction of a robust set of reference exRNA profiles for a variety of biofluids. Here, we present progress by ERCC investigators in these areas, and we discuss collaborative projects directed at development of robust methods for EV/exRNA isolation and analysis and tools for sharing and computational analysis of exRNA profiling data.

LIN28B Impairs the Transition of hESC-Derived β Cells from the Juvenile to Adult State

Differentiation of human embryonic stem cells into pancreatic β cells holds great promise for the treatment of diabetes. Recent advances have led to the production of glucose-responsive insulin-secreting cells in vitro, but resulting cells remain less mature than their adult primary β cell counterparts. The barrier(s) to in vitro β cell maturation are unclear. Here, we evaluated a potential role for microRNAs. MicroRNA profiling showed high expression of let-7 family microRNAs in vivo, but not in in vitro differentiated β cells. Reduced levels of let-7 in vitro were associated with increased levels of the RNA binding protein LIN28B, a negative regulator of let-7 biogenesis. Ablation of LIN28B during human embryonic stem cell (hESC) differentiation toward β cells led to a more mature glucose-stimulated insulin secretion profile and the suppression of juvenile-specific genes. However, let-7 overexpression had little effect. These results uncover LIN28B as a modulator of β cell maturation in vitro.

Comprehensive transcriptomic analyses of tissue, serum, and serum exosomes from hepatocellular carcinoma patients

BACKGROUND

The expression of microRNAs (miRNAs) is a promising prognostic and diagnostic tool in hepatocellular carcinoma (HCC). Here we performed small RNA sequencing (sRNA-seq) of tissue, serum and serum exosomes to investigate changes in miRNA expression between the different sample types and correlated the expression with clinical parameters. We also performed gene expression arrays on tumor and normal tissue.

RESULTS

Paired tissue, serum and serum exosomes sequencing revealed consistent positive correlation of miR-21 between serum exosomes and tumor tissue, indicating that miR-21 could be exported from tissue to circulation via exosomes. We found that let-7 miRNAs are generally upregulated in serum exosomes compared to whole serum, indicating that these miRNAs could be preferentially loaded into exosomes. Comparing serum from HCC patients with serum from healthy individuals revealed a global increase of miRNAs in serum from HCC patients, including an almost 4-fold increase of several miRNAs, including the liver-specific miR-122. When correlating miRNA expression with clinical parameters we detected significant association between hepatitis B virus (HBV) infection and miR-122 in serum as well as several serum and tissue-miRNAs that correlated with surgery type. We found that miR-141 and miR-146 correlated with cirrhosis in tumor tissue and normal tissue, respectively. Finally, high expression of miR-21 in tumors were associated with poor survival. Focusing on gene expression we found several significant messenger RNAs (mRNAs) between tumor and normal tissue and a Gene Ontology (GO) analysis revealed that these changes were mainly related to cell cycle and metabolism. Further, we detected mRNAs that correlated with cirrhosis and HBV infection in tissue. Finally, GO analysis of predicted targets for miRNAs down-regulated in tumor found that these were enriched for functions related to collagen synthesis.

CONCLUSIONS

Our combined data point to altered miRNA and mRNA expression contributing to both generally impaired lipid metabolism and increased cell proliferation and a miRNA-driven increase in collagen synthesis in HCC. Our results further indicate a correlation in miRNA expression between exosomes, serum, and tissue samples suggesting export from tumors via exosomes. This correlation could provide a basis for a more tumor-specific miRNA profile in serum.

Computational Workflow for Small RNA Profiling in Virus-Infected Plants

In this chapter we describe a series of computational pipelines for the in silico analysis of small RNAs (sRNA) produced in response to viral infections in plants. Our workflow is primarily focused on the analysis of sRNA populations derived from known or previously undescribed viruses infecting host plants. Furthermore, we provide an additional pipeline to examine host-specific endogenous sRNAs activated or specifically expressed during viral infections in plants. We present some key points for a successful and cost-efficient processing of next generation sequencing sRNA libraries, from purification of high quality RNA to guidance for library preparation and sequencing strategies. We report a series of free available tools and programs as well as in-house Perl scripts to perform customized sRNA-seq data mining. Previous bioinformatic background is not required, but experience with basic Unix commands is desirable.

miR-450a Acts as a Tumor Suppressor in Ovarian Cancer by Regulating Energy Metabolism

Dysregulation of miRNA expression is associated with multiple diseases, including cancers, in which small RNAs can have either oncogenic or tumor suppressive functions. Here we investigated the potential tumor suppressive function of miR-450a, one of the most significantly downregulated miRNAs in ovarian cancer. RNA-seq analysis of the ovarian cancer cell line A2780 revealed that overexpression of miR-450a suppressed multiple genes involved in the epithelial-to-mesenchymal transition (EMT). Overexpression of miR-450a reduced tumor migration and invasion and increased anoikis in A2780 and SKOV-3 cell lines and reduced tumor growth in an ovarian tumor xenographic model. Combined AGO-PAR-CLIP and RNA-seq analysis identified a panel of potential miR-450a targets, of which many, including TIMMDC1, MT-ND2, ACO2, and ATP5B, regulate energetic metabolism. Following glutamine withdrawal, miR-450a overexpression decreased mitochondrial membrane potential but increased glucose uptake and viability, characteristics of less invasive ovarian cancer cell lines. In summary, we propose that miR-450a acts as a tumor suppressor in ovarian cancer cells by modulating targets associated with glutaminolysis, which leads to decreased production of lipids, amino acids, and nucleic acids, as well as inhibition of signaling pathways associated with EMT. SIGNIFICANCE: miR-450a limits the metastatic potential of ovarian cancer cells by targeting a set of mitochondrial mRNAs to reduce glycolysis and glutaminolysis.Graphical Abstract: http://cancerres.aacrjournals.org/content/canres/79/13/3294/F1.large.jpg.

Detection of circulating extracellular mRNAs by modified small-RNA-sequencing analysis

Extracellular mRNAs (ex-mRNAs) potentially supersede extracellular miRNAs (ex-miRNAs) and other RNA classes as biomarkers. We performed conventional small-RNA-sequencing (sRNA-seq) and sRNA-seq with T4 polynucleotide kinase (PNK) end-treatment of total exRNA isolated from serum and platelet-poor EDTA, ACD, and heparin plasma to study the effect on ex-mRNA capture. Compared to conventional sRNA-seq PNK-treatment increased the detection of informative ex-mRNAs reads up to 50-fold. The exRNA pool was dominated by hematopoietic cells and platelets, with additional contribution from the liver. About 60% of the 15- to 42-nt reads originated from the coding sequences, in a pattern reminiscent of ribosome-profiling. Blood sample type had a considerable influence on the exRNA profile. On average approximately 350 to 1,100 distinct ex-mRNA transcripts were detected depending on plasma type. In serum, additional transcripts from neutrophils and hematopoietic cells increased this number to near 2,300. EDTA and ACD plasma showed a destabilizing effect on ex mRNA and non-coding RNA ribonucleoprotein complexes compared to other plasma types. In a proof-of-concept study, we investigated differences between the exRNA profiles of patients with acute coronary syndrome (ACS) and healthy controls. The improved tissue resolution of ex mRNAs after PNK-treatment enabled us to detect a neutrophil-signature in ACS that escaped detection by ex miRNA analysis.

The splicing factor U2AF1 contributes to cancer progression through a noncanonical role in translation regulation

Palangat et al. uncovered a noncanonical function of U2AF1, showing that it directly binds mature mRNA in the cytoplasm and negatively regulates mRNA translation.

Somatic mutations in the genes encoding components of the spliceosome occur frequently in human neoplasms, including myeloid dysplasias and leukemias, and less often in solid tumors. One of the affected factors, U2AF1, is involved in splice site selection, and the most common change, S34F, alters a conserved nucleic acid-binding domain, recognition of the 3′ splice site, and alternative splicing of many mRNAs. However, the role that this mutation plays in oncogenesis is still unknown. Here, we uncovered a noncanonical function of U2AF1, showing that it directly binds mature mRNA in the cytoplasm and negatively regulates mRNA translation. This splicing-independent role of U2AF1 is altered by the S34F mutation, and polysome profiling indicates that the mutation affects translation of hundreds of mRNA. One functional consequence is increased synthesis of the secreted chemokine interleukin 8, which contributes to metastasis, inflammation, and cancer progression in mice and humans.

Small RNA-seq: The RNA 5'-end adapter ligation problem and how to circumvent it

The preparation of small RNA cDNA sequencing libraries depends on the unbiased ligation of adapters to the RNA ends. Small RNA with 5' recessed ends are poor substrates for enzymatic adapter ligation, but this 5' adapter ligation problem can go undetected if the library preparation steps are not monitored. Here we illustrate the severity of the 5' RNA end ligation problem using several pre-miRNA-like hairpins that allow us to expand the definition of the problem to include 5' ends close to a hairpin stem, whether recessed or in a short extension. The ribosome profiling method can avoid a difficult 5' adapter ligation, but the enzyme typically used to circularize the cDNA has been reported to be biased, calling into question the benefit of this workaround. Using the TS2126 RNA ligase 1 (a.k.a. CircLigase) as the circularizing enzyme, we devised a bias test for the circularization of first strand cDNA. All possible dinucleotides were circle-ligated with similar efficiency. To re-linearize the first strand cDNA in the ribosome profiling approach, we introduce an improved method wherein a single ribonucleotide is placed between the sequencing primer binding sites in the reverse transcriptase primer, which later serves as the point of re-linearization by RNase A. We incorporate this step into the ribosomal profiling method and describe a complete improved library preparation method, Coligo-seq, for the sequencing of small RNA with secondary structure close to the 5' end. This method accepts a variety of 5' modified RNA, including 5' monophosphorylated RNA, as demonstrated by the construction of a HeLa cell microRNA cDNA library.

High-throughput, Efficient, and Unbiased Capture of Small RNAs from Low-input Samples for Sequencing

MicroRNAs hold great promise as biomarkers of disease. However, there are few efficient and robust methods for measuring microRNAs from low input samples. Here, we develop a high-throughput sequencing protocol that efficiently captures small RNAs while minimizing inherent biases associated with library production. The protocol is based on early barcoding such that all downstream manipulations can be performed on a pool of many samples thereby reducing reagent usage and workload. We show that the optimization of adapter concentrations along with the addition of nucleotide modifications and random nucleotides increases the efficiency of small RNA capture. We further show, using unique molecular identifiers, that stochastic capture of low input RNA rather than PCR amplification influences the biased quantitation of intermediately and lowly expressed microRNAs. Our improved method allows the processing of tens to hundreds of samples simultaneously while retaining high efficiency quantitation of microRNAs in low input samples from tissues or bodily fluids.

Evaluating gastroenteropancreatic neuroendocrine tumors through microRNA sequencing

Gastroenteropancreatic neuroendocrine tumors (GEP-NETs) can be challenging to evaluate histologically. MicroRNAs (miRNAs) are small RNA molecules that often are excellent biomarkers due to their abundance, cell-type and disease stage specificity and stability. To evaluate miRNAs as adjunct tissue markers for classifying and grading well-differentiated GEP-NETs, we generated and compared miRNA expression profiles from four pathological types of GEP-NETs. Using quantitative barcoded small RNA sequencing and state-of-the-art sequence annotation, we generated comprehensive miRNA expression profiles from archived pancreatic, ileal, appendiceal and rectal NETs. Following data preprocessing, we randomly assigned sample profiles to discovery (80%) and validation (20%) sets prior to data mining using machine-learning techniques. High expression analyses indicated that miR-375 was the most abundant individual miRNA and miRNA cistron in all samples. Leveraging prior knowledge that GEP-NET behavior is influenced by embryonic derivation, we developed a dual-layer hierarchical classifier for differentiating GEP-NET types. In the first layer, our classifier discriminated midgut (ileum, appendix) from non-midgut (rectum, pancreas) NETs based on miR-615 and -92b expression. In the second layer, our classifier discriminated ileal from appendiceal NETs based on miR-125b, -192 and -149 expression, and rectal from pancreatic NETs based on miR-429 and -487b expression. Our classifier achieved overall accuracies of 98.5% and 94.4% in discovery and validation sets, respectively. We also found provisional evidence that low- and intermediate-grade pancreatic NETs can be discriminated based on miR-328 expression. GEP-NETs can be reliably classified and potentially graded using a limited panel of miRNA markers, complementing morphological and immunohistochemistry-based approaches to histologic evaluation.

MicroRNAs and Transplantation

miRNAs, ∼20 to 22 nucleotide single-stranded RNA species that play a pivotal role in the regulation of protein-coding genes, are emerging as robust biomarkers for assessing allograft status. Herein, the authors briefly review the biogenesis and function of the miRNAs and provide an overview of the tools to quantify miRNAs in tissues and body fluids. They then review their studies of discovery and validation of alterations in miRNA expression within kidney allografts with or without acute rejection, as well as with or without fibrosis, and summarize published data on miRNA expression patterns in kidney transplant recipients.

Proximity-CLIP provides a snapshot of protein-occupied RNA elements in subcellular compartments

Methods to systematically study subcellular RNA localization are limited and lagging behind proteomic tools. Here, we combined APEX2-mediated proximity biotinylation of proteins with photoactivatable ribonucleoside-enhanced crosslinking to simultaneously profile the proteome, as well as the transcriptome bound by RNA-binding proteins in any given subcellular compartment. Our approach is fractionation-independent and enables to study the localization of RNA processing intermediates, as well as the identification of regulatory RNA cis-acting elements occupied by proteins in a cellular compartment-specific manner. We applied Proximity-CLIP to study RNA and protein in the nucleus, cytoplasm and at cell-cell interfaces. Among other insights, we observed frequent transcriptional readthrough continuing for several kilobases downstream of the canonical cleavage and polyadenylation site and a differential RBP occupancy pattern for mRNAs in the nucleus and cytoplasm. Surprisingly, mRNAs localized to cell-cell interfaces often encoded regulatory proteins and contained protein-occupied CUG sequence elements in their 3’ untranslated region.

A highly parallel strategy for storage of digital information in living cells

BACKGROUND

Encoding arbitrary digital information in DNA has attracted attention as a potential avenue for large scale and long term data storage. However, in order to enable DNA data storage technologies there needs to be improvements in data storage fidelity (tolerance to mutation), the facility of writing and reading the data (biases and systematic error arising from synthesis and sequencing), and overall scalability.

RESULTS

To this end, we have developed and implemented an encoding scheme that is suitable for detecting and correcting errors that may arise during storage, writing, and reading, such as those arising from nucleotide substitutions, insertions, and deletions. We propose a scheme for parallelized long term storage of encoded sequences that relies on overlaps rather than the address blocks found in previously published work. Using computer simulations, we illustrate the encoding, sequencing, decoding, and recovery of encoded information, ultimately demonstrating the possibility of a successful round-trip read/write. These demonstrations show that in theory a precise control over error tolerance is possible. Even after simulated degradation of DNA, recovery of original data is possible owing to the error correction capabilities built into the encoding strategy. A secondary advantage of our method is that the statistical characteristics (such as repetitiveness and GC-composition) of encoded sequences can also be tailored without sacrificing the overall ability to store large amounts of data. Finally, the combination of the overlap-based partitioning of data with the LZMA compression that is integral to encoding means that the entire sequence must be present for successful decoding. This feature enables inordinately strong encryptions. As a potential application, an encrypted pathogen genome could be distributed and carried by cells without danger of being expressed, and could not even be read out in the absence of the entire DNA consortium.

CONCLUSIONS

We have developed a method for DNA encoding, using a significantly different fundamental approach from existing work, which often performs better than alternatives and allows for a great deal of freedom and flexibility of application.

CD95/Fas ligand mRNA is toxic to cells

CD95/Fas ligand binds to the death receptor CD95 to induce apoptosis in sensitive cells. We previously reported that CD95L mRNA is enriched in sequences that, when converted to si/shRNAs, kill all cancer cells by targeting critical survival genes (Putzbach et al., 2017). We now report expression of full-length CD95L mRNA itself is highly toxic to cells and induces a similar form of cell death. We demonstrate that small (s)RNAs derived from CD95L are loaded into the RNA induced silencing complex (RISC) which is required for the toxicity and processing of CD95L mRNA into sRNAs is independent of both Dicer and Drosha. We provide evidence that in addition to the CD95L transgene a number of endogenous protein coding genes involved in regulating protein translation, particularly under low miRNA conditions, can be processed to sRNAs and loaded into the RISC suggesting a new level of cell fate regulation involving RNAi.

New insights into the cellular temporal response to proteostatic stress

Maintaining a healthy proteome involves all layers of gene expression regulation. By quantifying temporal changes of the transcriptome, translatome, proteome, and RNA-protein interactome in cervical cancer cells, we systematically characterize the molecular landscape in response to proteostatic challenges. We identify shared and specific responses to misfolded proteins and to oxidative stress, two conditions that are tightly linked. We reveal new aspects of the unfolded protein response, including many genes that escape global translation shutdown. A subset of these genes supports rerouting of energy production in the mitochondria. We also find that many genes change at multiple levels, in either the same or opposing directions, and at different time points. We highlight a variety of putative regulatory pathways, including the stress-dependent alternative splicing of aminoacyl-tRNA synthetases, and protein-RNA binding within the 3’ untranslated region of molecular chaperones. These results illustrate the potential of this information-rich resource.

Argonaute-miRNA Complexes Silence Target mRNAs in the Nucleus of Mammalian Stem Cells

In mammals, gene silencing by the RNA-induced silencing complex (RISC) is a well-understood cytoplasmic posttranscriptional gene regulatory mechanism. Here, we show that embryonic stem cells (ESCs) contain high levels of nuclear AGO proteins and that in ESCs nuclear AGO protein activity allows for the onset of differentiation. In the nucleus, AGO proteins interact with core RISC components, including the TNRC6 proteins and the CCR4-NOT deadenylase complex. In contrast to cytoplasmic miRNA-mediated gene silencing that mainly operates on cis-acting elements in mRNA 3' untranslated (UTR) sequences, in the nucleus AGO binding in the coding sequence and potentially introns also contributed to post-transcriptional gene silencing. Thus, nuclear localization of AGO proteins in specific cell types leads to a previously unappreciated expansion of the miRNA-regulated transcriptome.

Roquin targets mRNAs in a 3'-UTR-specific manner by different modes of regulation

The RNA-binding proteins Roquin-1 and Roquin-2 redundantly control gene expression and cell-fate decisions. Here, we show that Roquin not only interacts with stem–loop structures, but also with a linear sequence element present in about half of its targets. Comprehensive analysis of a minimal response element of the Nfkbid 3′-UTR shows that six stem–loop structures cooperate to exert robust and profound post-transcriptional regulation. Only binding of multiple Roquin proteins to several stem–loops exerts full repression, which redundantly involved deadenylation and decapping, but also translational inhibition. Globally, most Roquin targets are regulated by mRNA decay, whereas a small subset, including the Nfat5 mRNA, with more binding sites in their 3′-UTRs, are also subject to translational inhibition. These findings provide insights into how the robustness and magnitude of Roquin-mediated regulation is encoded in complex cis-elements.

Roquin targets are known to contain two types of sequence-structure motifs, the constitutive and the alternative decay elements (CDE and ADE). Here, the authors describe a linear Roquin binding element (LBE) also involved in target recognition, and show that Roquin binding affects the translation of a subset of targeted mRNAs.

Phosphorylation of the Ribosomal Protein RPL12/uL11 Affects Translation during Mitosis

Emerging evidence indicates that heterogeneity in ribosome composition can give rise to specialized functions. Until now, research mainly focused on differences in core ribosomal proteins and associated factors. The effect of posttranslational modifications has not been studied systematically. Analyzing ribosome heterogeneity is challenging because individual proteins can be part of different subcomplexes (40S, 60S, 80S, and polysomes). Here we develop polysome proteome profiling to obtain unbiased proteomic maps across ribosomal subcomplexes. Our method combines extensive fractionation by sucrose gradient centrifugation with quantitative mass spectrometry. The high resolution of the profiles allows us to assign proteins to specific subcomplexes. Phosphoproteomics on the fractions reveals that phosphorylation of serine 38 in RPL12/uL11, a known mitotic CDK1 substrate, is strongly depleted in polysomes. Follow-up experiments confirm that RPL12/uL11 phosphorylation regulates the translation of specific subsets of mRNAs during mitosis. Together, our results show that posttranslational modification of ribosomal proteins can regulate translation.

Identification and removal of sequencing artifacts produced by mispriming during reverse transcription in multiple RNA-seq technologies

The quality of RNA sequencing data relies on specific priming by the primer used for reverse transcription (RT-primer). Nonspecific annealing of the RT-primer to the RNA template can generate reads with incorrect cDNA ends and can cause misinterpretation of data (RT mispriming). This kind of artifact in RNA-seq based technologies is underappreciated and currently no adequate tools exist to computationally remove them from published data sets. We show that mispriming can occur with as little as two bases of complementarity at the 3' end of the primer followed by intermittent regions of complementarity. We also provide a computational pipeline that identifies cDNA reads produced from RT mispriming, allowing users to filter them out from any aligned data set. Using this analysis pipeline, we identify thousands of mispriming events in a dozen published data sets from diverse technologies including short RNA-seq, total/mRNA-seq, HITS-CLIP, and GRO-seq. We further show how RT mispriming can lead to misinterpretation of data. In addition to providing a solution to computationally remove RT-misprimed reads, we also propose an experimental solution to completely avoid RT-mispriming by performing RNA-seq using thermostable group II intron derived reverse transcriptase (TGIRT-seq).

miRge 2.0 for comprehensive analysis of microRNA sequencing data

Background

miRNAs play important roles in the regulation of gene expression. The rapidly developing field of microRNA sequencing (miRNA-seq; small RNA-seq) needs comprehensive, robust, user-friendly and standardized bioinformatics tools to analyze these large datasets. We present miRge 2.0, in which multiple enhancements were made towards these goals.

Results

miRge 2.0 has become more comprehensive with increased functionality including a novel miRNA detection method, A-to-I editing analysis, integrated standardized GFF3 isomiR reporting, and improved alignment to miRNAs. The novel miRNA detection method uniquely uses both miRNA hairpin sequence structure and composition of isomiRs resulting in higher specificity for potential miRNA identification. Using known miRNA data, our support vector machine (SVM) model predicted miRNAs with an average Matthews correlation coefficient (MCC) of 0.939 over 32 human cell datasets and outperformed miRDeep2 and miRAnalyzer regarding phylogenetic conservation. The A-to-I editing detection strongly correlated with a reference dataset with adjusted R² = 0.96. miRge 2.0 is the most up-to-date aligner with custom libraries to both miRBase v22 and MirGeneDB v2.0 for 6 species: human, mouse, rat, fruit fly, nematode and zebrafish; and has a tool to create custom libraries. For user-friendliness, miRge 2.0 is incorporated into bcbio-nextgen and implementable through Bioconda.

Conclusions

miRge 2.0 is a redesigned, leading miRNA RNA-seq aligner with several improvements and novel utilities. miRge 2.0 is freely available at: https://github.com/mhalushka/miRge.

Electronic supplementary material

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

Human plasma and serum extracellular small RNA reference profiles and their clinical utility

Nucleic acids mediate storage and expression of genetic information. Extracellular DNA (exDNA) and exRNA are traces of nucleic acids released from cells into the extracellular environment. Their use as disease biomarkers has been limited by technical challenges in their isolation caused by abundant RNA- and DNA-degrading enzymes in biofluids. Using isolation protocols developed especially for biofluids, we generated plasma and serum exRNA reference profiles from 13 healthy volunteers over time and determined the effect of critical clinical parameters such as gender and fasting. Surprisingly, we encountered one participant with dramatically increased endocrine-origin exRNA contributions stable over 1 year and detectable in all of his samples, thereby demonstrating the robustness of this approach and the clinical potential of circulating RNAs as biomarkers.

Circulating extracellular RNAs (exRNAs) have the potential to serve as biomarkers for a wide range of medical conditions. However, limitations in existing exRNA isolation methods and a lack of knowledge on parameters affecting exRNA variability in human samples may hinder their successful discovery and clinical implementation. Using combinations of denaturants, reducing agents, proteolysis, and revised organic extraction, we developed an automated, high-throughput approach for recovery of exRNAs and exDNA from the same biofluid sample. We applied this method to characterize exRNAs from 312 plasma and serum samples collected from 13 healthy volunteers at 12 time points over a 2-month period. Small RNA cDNA library sequencing identified nearly twofold increased epithelial-, muscle-, and neuroendocrine-cell–specific miRNAs in females, while fasting and hormonal cycle showed little effect. External standardization helped to detect quantitative differences in erythrocyte and platelet-specific miRNA contributions and in miRNA concentrations between biofluids. It also helped to identify a study participant with a unique exRNA phenotype featuring a miRNA signature of up to 20-fold elevated endocrine-cell–specific miRNAs and twofold elevated total miRNA concentrations stable for over 1 year. Collectively, these results demonstrate an efficient and quantitative method to discern exRNA phenotypes and suggest that plasma and serum RNA profiles are stable over months and can be routinely monitored in long-term clinical studies.

Retrospective MicroRNA Sequencing: Complementary DNA Library Preparation Protocol Using Formalin-fixed Paraffin-embedded RNA Specimens

-Archived, clinically classified formalin-fixed paraffin-embedded (FFPE) tissues can provide nucleic acids for retrospective molecular studies of cancer development. By using non-invasive or pre-malignant lesions from patients who later develop invasive disease, gene expression analyses may help identify early molecular alterations that predispose to cancer risk. It has been well described that nucleic acids recovered from FFPE tissues have undergone severe physical damage and chemical modifications, which make their analysis difficult and generally requires adapted assays. MicroRNAs (miRNAs), however, which represent a small class of RNA molecules spanning only up to ~18-24 nucleotides, have been shown to withstand long-term storage and have been successfully analyzed in FFPE samples. Here we present a 3' barcoded complementary DNA (cDNA) library preparation protocol specifically optimized for the analysis of small RNAs extracted from archived tissues, which was recently demonstrated to be robust and highly reproducible when using archived clinical specimens stored for up to 35 years. This library preparation is well adapted to the multiplex analysis of compromised/degraded material where RNA samples (up to 18) are ligated with individual 3' barcoded adapters and then pooled together for subsequent enzymatic and biochemical preparations prior to analysis. All purifications are performed by polyacrylamide gel electrophoresis (PAGE), which allows size-specific selections and enrichments of barcoded small RNA species. This cDNA library preparation is well adapted to minute RNA inputs, as a pilot polymerase chain reaction (PCR) allows determination of a specific amplification cycle to produce optimal amounts of material for next-generation sequencing (NGS). This approach was optimized for the use of degraded FFPE RNA from specimens archived for up to 35 years and provides highly reproducible NGS data.

Modulation of LIN28B/Let-7 Signaling by Propranolol Contributes to Infantile Hemangioma Involution

OBJECTIVE

Infantile hemangiomas (IHs) are the most common benign vascular neoplasms of infancy, characterized by a rapid growth phase followed by a spontaneous involution, or triggered by propranolol treatment by poorly understood mechanisms. LIN28/let-7 axis plays a central role in the regulation of stem cell self-renewal and tumorigenesis. However, the role of LIN28B/let-7 signaling in IH pathogenesis has not yet been elucidated.

APPROACH AND RESULTS

LIN28B is highly expressed in proliferative IH and is less expressed in involuted and in propranolol-treated IH samples as measured by immunofluorescence staining and quantitative RT-PCR. Small RNA sequencing analysis of IH samples revealed a decrease in microRNAs that target LIN28B, including let-7, and an increase in microRNAs in the mir-498(46) cistron. Overexpression of LIN28B in HEK293 cells induced the expression of miR-516b in the mir-498(46) cistron. Propranolol treatment of induced pluripotent stem cells, which express mir-498(46) endogenously, reduced the expression of both LIN28B and mir-498(46) and increased the expression of let-7. Furthermore, propranolol treatment reduced the proliferation of induced pluripotent stem cells and induced epithelial-mesenchymal transition.

CONCLUSIONS

This work uncovers the role of the LIN28B/let-7 switch in IH pathogenesis and provides a novel mechanism by which propranolol induces IH involution. Furthermore, it provides therapeutic implications for cancers in which the LIN28/let-7 pathway is imbalanced.

Defining a microRNA-mRNA interaction map for calcineurin inhibitor induced nephrotoxicity

Calcineurin inhibitors induce nephrotoxicity through poorly understood mechanisms thereby limiting their use in transplantation and other diseases. Here we define a microRNA (miRNA)-messenger RNA (mRNA) interaction map that facilitates exploration into the role of miRNAs in cyclosporine-induced nephrotoxicity (CIN) and the gene pathways they regulate. Using photoactivatable ribonucleoside-enhanced crosslinking and immunoprecipitation (PAR-CLIP), we isolated RNAs associated with Argonaute 2 in the RNA-induced silencing complex (RISC) of cyclosporine A (CsA) treated and control human proximal tubule cells and identified mRNAs undergoing active targeting by miRNAs. CsA causes specific changes in miRNAs and mRNAs associated with RISC, thereby altering post-transcriptional regulation of gene expression. Pathway enrichment analysis identified canonical pathways regulated by miRNAs specifically following CsA treatment. RNA-seq performed on total RNA indicated that only a fraction of total miRNAs and mRNAs are actively targeted in the RISC, indicating that PAR-CLIP more accurately defines meaningful targeting interactions. Our data also revealed a role for miRNAs in calcineurin-independent regulation of JNK and p38 MAPKs caused by targeting of MAP3K1. Together, our data provide a novel resource and unique insights into molecular pathways regulated by miRNAs in CIN. The gene pathways and miRNAs defined may represent novel targets to reduce calcineurin induced nephrotoxicity.

Uncovering cell type-specific complexities of gene expression and RNA metabolism by TU-tagging and EC-tagging

Cell type-specific transcription is a key determinant of cell fate and function. An ongoing challenge in biology is to develop robust and stringent biochemical methods to explore gene expression with cell type specificity. This challenge has become even greater as researchers attempt to apply high-throughput RNA analysis methods under in vivo conditions. TU-tagging and EC-tagging are in vivo biosynthetic RNA tagging techniques that allow spatial and temporal specificity in RNA purification. Spatial specificity is achieved through targeted expression of pyrimidine salvage enzymes (uracil phosphoribosyltransferase and cytosine deaminase) and temporal specificity is achieved by controlling exposure to bioorthogonal substrates of these enzymes (4-thiouracil and 5-ethynylcytosine). Tagged RNAs can be purified from total RNA extracted from an animal or tissue and used in transcriptome profiling analyses. In addition to identifying cell type-specific mRNA profiles, these techniques are applicable to noncoding RNAs and can be used to measure RNA transcription and decay. Potential applications of TU-tagging and EC-tagging also include fluorescent RNA imaging and selective definition of RNA-protein interactions. TU-tagging and EC-tagging hold great promise for supporting research at the intersection of RNA biology and developmental biology. This article is categorized under: Technologies > Analysis of the Transcriptome.

Roquin Suppresses the PI3K-mTOR Signaling Pathway to Inhibit T Helper Cell Differentiation and Conversion of Treg to Tfr Cells

Roquin proteins preclude spontaneous T cell activation and aberrant differentiation of T follicular helper (Tfh) or T helper 17 (Th17) cells. Here we showed that deletion of Roquin-encoding alleles specifically in regulatory T (Treg) cells also caused the activation of conventional T cells. Roquin-deficient Treg cells downregulated CD25, acquired a follicular Treg (Tfr) cell phenotype, and suppressed germinal center reactions but could not protect from colitis. Roquin inhibited the PI3K-mTOR signaling pathway by upregulation of Pten through interfering with miR-17∼92 binding to an overlapping cis-element in the Pten 3' UTR, and downregulated the Foxo1-specific E3 ubiquitin ligase Itch. Loss of Roquin enhanced Akt-mTOR signaling and protein synthesis, whereas inhibition of PI3K or mTOR in Roquin-deficient T cells corrected enhanced Tfh and Th17 or reduced iTreg cell differentiation. Thereby, Roquin-mediated control of PI3K-mTOR signaling prevents autoimmunity by restraining activation and differentiation of conventional T cells and specialization of Treg cells.

The Human CCHC-type Zinc Finger Nucleic Acid-Binding Protein Binds G-Rich Elements in Target mRNA Coding Sequences and Promotes Translation

The CCHC-type zinc finger nucleic acid-binding protein (CNBP/ZNF9) is conserved in eukaryotes and is essential for embryonic development in mammals. It has been implicated in transcriptional, as well as post-transcriptional, gene regulation; however, its nucleic acid ligands and molecular function remain elusive. Here, we use multiple systems-wide approaches to identify CNBP targets and function. We used photoactivatable ribonucleoside-enhanced crosslinking and immunoprecipitation (PAR-CLIP) to identify 8,420 CNBP binding sites on 4,178 mRNAs. CNBP preferentially bound G-rich elements in the target mRNA coding sequences, most of which were previously found to form G-quadruplex and other stable structures in vitro. Functional analyses, including RNA sequencing, ribosome profiling, and quantitative mass spectrometry, revealed that CNBP binding did not influence target mRNA abundance but rather increased their translational efficiency. Considering that CNBP binding prevented G-quadruplex structure formation in vitro, we hypothesize that CNBP is supporting translation by resolving stable structures on mRNAs.

A comprehensive, cell specific microRNA catalogue of human peripheral blood

With this study, we provide a comprehensive reference dataset of detailed miRNA expression profiles from seven types of human peripheral blood cells (NK cells, B lymphocytes, cytotoxic T lymphocytes, T helper cells, monocytes, neutrophils and erythrocytes), serum, exosomes and whole blood. The peripheral blood cells from buffy coats were typed and sorted using FACS/MACS. The overall dataset was generated from 450 small RNA libraries using high-throughput sequencing. By employing a comprehensive bioinformatics and statistical analysis, we show that 3′ trimming modifications as well as composition of 3′ added non-templated nucleotides are distributed in a lineage-specific manner—the closer the hematopoietic progenitors are, the higher their similarities in sequence variation of the 3′ end. Furthermore, we define the blood cell-specific miRNA and isomiR expression patterns and identify novel cell type specific miRNA candidates. The study provides the most comprehensive contribution to date towards a complete miRNA catalogue of human peripheral blood, which can be used as a reference for future studies. The dataset has been deposited in GEO and also can be explored interactively following this link: http://134.245.63.235/ikmb-tools/bloodmiRs.

Identification of metastasis-associated microRNAs in serum from rectal cancer patients

MicroRNAs (miRNAs) are promising prognostic and diagnostic biomarkers due to their high stability in blood. Here we investigate the expression of miRNAs and other noncoding (nc) RNAs in serum of rectal cancer patients. Serum from 96 rectal cancer patients was profiled using small RNA sequencing and expression of small RNAs was correlated with the clinicopathological characteristics of the patients. Multiple classes of RNAs were detected, including miRNAs and fragments of tRNAs, snoRNAs, long ncRNAs, and other classes of RNAs. Several miRNAs, miRNA variants (isomiRs) and other ncRNAs were differentially expressed between Stage IV and Stage I-III rectal cancer patients, including several members of the miR-320 family. Furthermore, we show that high expression of miR-320d as well as one tRNA fragment is associated with poor survival. We also show that several miRNAs and isomiRs are differentially expressed between patients receiving preoperative chemoradiotherapy and patients who did not receive any treatment before serum collection. In summary, our study shows that the expression of miRNAs and other small ncRNAs in serum may be used to predict distant metastasis and survival in rectal cancer.

Alterations of microRNAs throughout the malignant evolution of cutaneous squamous cell carcinoma: the role of miR-497 in epithelial to mesenchymal transition of keratinocytes

Skin carcinogenesis is known to be a multi-step process with several stages along its malignant evolution. We hypothesized that transformation of normal epidermis to cutaneous squamous cell carcinoma (cSCC) is causally linked to alterations in microRNAs (miRNA) expression. For this end we decided to evaluate their alterations in the pathologic states ending in cSCC. Total RNA was extracted from formalin fixed paraffin embedded biopsies of five stages along the malignant evolution of keratinocytes towards cSCC: Normal epidermis, solar elastosis, actinic keratosis KIN1-2, advanced actinic keratosis KIN3 and well-differentiated cSCC. Next-generation small RNA sequencing was performed. We found that 18 miRNAs are overexpressed and 28 miRNAs are underexpressed in cSCC compared to normal epidermis. miR-424, miR-320, miR-222 and miR-15a showed the highest fold change among the overexpressed miRNAs. And miR-100, miR-101 and miR-497 showed the highest fold change among the underexpressed miRNAs. Heat map of hierarchical clustering analysis of significantly changed miRNAs and principle component analysis disclosed that the most prominent change in miRNAs expression occurred in the switch from 'early' stages; normal epidermis, solar elastosis and early actinic keratosis to the 'late' stages of epidermal carcinogenesis; late actinic keratosis and cSCC. We found several miRNAs with 'stage specific' alterations while others display a clear 'gradual', either progressive increase or decrease in expression along the malignant evolution of keratinocytes. The observed alterations focused in miRNAs involved in the regulation of AKT/mTOR or in those involved in epithelial to mesenchymal transition. We chose to concentrate on the evaluation of the molecular role of miR-497. We found that it induces reversion of epithelial to mesenchymal transition. We proved that SERPINE-1 is its biochemical target. The present study allows us to further study the pathways that are regulated by miRNAs along the malignant evolution of keratinocytes towards cSCC.

LARP4 mRNA codon-tRNA match contributes to LARP4 activity for ribosomal protein mRNA poly(A) tail length protection

Messenger RNA function is controlled by the 3' poly(A) tail (PAT) and poly(A)-binding protein (PABP). La-related protein-4 (LARP4) binds poly(A) and PABP. LARP4 mRNA contains a translation-dependent, coding region determinant (CRD) of instability that limits its expression. Although the CRD comprises <10% of LARP4 codons, the mRNA levels vary >20 fold with synonymous CRD substitutions that accommodate tRNA dynamics. Separately, overexpression of the most limiting tRNA increases LARP4 levels and reveals its functional activity, net lengthening of the PATs of heterologous mRNAs with concomitant stabilization, including ribosomal protein (RP) mRNAs. Genetic deletion of cellular LARP4 decreases PAT length and RPmRNA stability. This LARP4 activity requires its PABP-interaction domain and the RNA-binding module which we show is sensitive to poly(A) 3'-termini, consistent with protection from deadenylation. The results indicate that LARP4 is a posttranscriptional regulator of ribosomal protein production in mammalian cells and suggest that this activity can be controlled by tRNA levels.