Deep-Sequencing of Endothelial Cells Exposed to Hypoxia Reveals the Complexity of Known and Novel Micrornas

Total Page:16

File Type:pdf, Size:1020Kb

Deep-Sequencing of Endothelial Cells Exposed to Hypoxia Reveals the Complexity of Known and Novel Micrornas Downloaded from rnajournal.cshlp.org on October 8, 2021 - Published by Cold Spring Harbor Laboratory Press Deep-sequencing of endothelial cells exposed to hypoxia reveals the complexity of known and novel microRNAs CHRISTINE VOELLENKLE,1,6 JEROEN VAN ROOIJ,1,6 ALESSANDRO GUFFANTI,2,6 ELENA BRINI,2 PASQUALE FASANARO,3 ELEONORA ISAIA,1 LARRY CROFT,4 MATEI DAVID,5 MAURIZIO C. CAPOGROSSI,3 ANNA MOLES,2 ARMANDO FELSANI,2 and FABIO MARTELLI3,7 1Laboratorio di Cardiologia Molecolare, IRCCS-Policlinico San Donato, San Donato Milanese, 20097 Milan, Italy 2Genomnia srl, Lainate, 20020 Milan, Italy 3Laboratorio di Patologia Vascolare, Istituto Dermopatico dell’Immacolata-IRCCS, 00167 Rome, Italy 4Institute for Molecular Bioscience, University of Queensland, 4072 Brisbane, Australia 5Computational Biology Laboratory, University of Toronto, M5S 3G4 Toronto, Canada ABSTRACT In order to understand the role of microRNAs (miRNAs) in vascular physiopathology, we took advantage of deep-sequencing techniques to accurately and comprehensively profile the entire miRNA population expressed by endothelial cells exposed to hypoxia. SOLiD sequencing of small RNAs derived from human umbilical vein endothelial cells (HUVECs) exposed to 1%% O2 or normoxia for 24 h yielded more than 22 million reads per library. A customized bioinformatic pipeline identified more than 400 annotated microRNA/microRNA* species with a broad abundance range: miR-21 and miR-126 totaled almost 40%% of all miRNAs. A complex repertoire of isomiRs was found, displaying also 59 variations, potentially affecting target recognition. High- stringency bioinformatic analysis identified microRNA candidates, whose predicted pre-miRNAs folded into a stable hairpin. Validation of a subset by qPCR identified 18 high-confidence novel miRNAs as detectable in independent HUVEC cultures and associated to the RISC complex. The expression of two novel miRNAs was significantly down-modulated by hypoxia, while miR- 210 was significantly induced. Gene ontology analysis of their predicted targets revealed a significant association to hypoxia- inducible factor signaling, cardiovascular diseases, and cancer. Overexpression of the novel miRNAs in hypoxic endothelial cells affected cell growth and confirmed the biological relevance of their down-modulation. In conclusion, deep-sequencing accurately profiled known, variant, and novel microRNAs expressed by endothelial cells in normoxia and hypoxia. Keywords: microRNA; isomiR; deep-sequencing; hypoxia; endothelium INTRODUCTION a ribonucleoprotein complex mediating post-transcriptional gene silencing. Sequence-specific interaction with the func- MicroRNAs (miRNAs) are endogenous, short, 21- to tional miRNA leads to destabilization or translational in- 23-nucleotide (nt), nonprotein-coding RNAs that act as neg- hibition of the target mRNAs by the Ago protein (Winter et al. ative regulators of gene expression, modulating the trans- 2009). For several pre-miRNAs, both hairpin arms are lational efficiency and/or the stability of target messenger processed and give rise to a functional miRNA. In this case, RNAs (mRNAs). These RNA molecules are transcribed as mature miRNAs are designated with the suffix ‘‘-3p’’ or ‘‘-5p,’’ a long primary transcript, termed pri-miRNA, sequentially depending on the originating hairpin arm. Alternatively, if processed by Drosha and Dicer ribonucleases, first to a pre- either one is expressed at much lower levels, it is designated as cursor-miRNA with a characteristic hairpin structure and miRNA* (star form) (Winter et al. 2009; Krol et al. 2010). then to the mature miRNA species. Then, the single-stranded, MiRNAs have essential functions in physiological pro- mature miRNA is loaded together with Argonaute (Ago) cesses, including cellular differentiation, development, and proteins into the RNA-induced silencing complex (RISC), apoptosis, but are also involved in pathological processes, like cancer (Garzon et al. 2009; McCormick et al. 2010) or cardiovascular diseases (Urbich et al. 2008; Dorn 2010; 6These authors contributed equally to this work. Fasanaro et al. 2010). 7 Corresponding author. Hypoxia describes the insufficient oxygen supply to tissues E-mail [email protected]. Article published online ahead of print. Article and publication date are and occurs during several physiopathological circumstances at http://www.rnajournal.org/cgi/doi/10.1261/rna.027615.111. such as rapid tissue growth, acute and chronic ischemia, RNA (2012), 18:00–00. Published by Cold Spring Harbor Laboratory Press. 1 Downloaded from rnajournal.cshlp.org on October 8, 2021 - Published by Cold Spring Harbor Laboratory Press Voellenkle et al. organ and tumor development, and high altitude (Semenza owing to the low probability of being miRNA derived. Then, 2010b). Decreased oxygen tension triggers a distinct pro- all remaining reads were trimmed to 25 nt, the maximum gram of responses aimed at relieving tissue hypoxia and expected length for a miRNA (Winter et al. 2009; Krol et al. removing irreversibly damaged cells. These responses include 2010), and used for mapping against various coding and endothelial cell proliferation, migration, and angiogenesis noncoding transcript and genomic data sets. and also growth arrest and apoptotic cell death (Semenza 2010b). Several reports have demonstrated that specific Identification and quantification of annotated miRNAs are involved in the hypoxic response (Ivan et al. miRNAs 2008). In particular, we and others found that the expression of miR-210 is induced in hypoxic endothelial The customized analysis pipeline represented in Figure 1 cells (Fasanaro et al. 2008; Chan et al. 2009), stimulating was used to quantify previously known (annotated) miR- their vascular endothelial growth factor (VEGF)–driven NAs. The reads corresponding to annotated miRNAs were migration and tubulogenesis (Fasanaro et al. 2008). identified by selecting all the reads mapping against the In order to understand how miRNAs regulate gene human sequences contained in miRBase v14.0, excluding expression in endothelial cells exposed to low-oxygen any read mapping also against a collection of other small tension, it is essential to accurately profile the entire miRNA RNAs databases (transfer RNAs [tRNAs], small nucleolar population. Currently, there are three widely used platforms RNAs [snoRNAs], ribosomal RNAs [rRNAs]). Then, fur- for miRNA profiling: quantitative polymerase chain reac- ther selection criteria based on alignment length, absence of tion (qPCR), microarray hybridization (Shendure 2008), and mismatches in the seed region, and absence of sequencing deep-sequencing (also named next-generation sequencing errors were applied: An average of 1.2 million reads per or massively parallel sequencing). While qPCR and micro- library were identified as annotated miRNAs, representing arrays rely on the availability of miRNA sequences in public 305 mature and 96 miRNA* species. As expected, a very and proprietary databases for designing probes or primers significant correlation between the results obtained analyzing (Griffiths-Jones et al. 2008), deep-sequencing offers the pos- the two normoxic and the two hypoxic libraries was observed sibility to comprehensively assess miRNA expression in a (Supplemental Fig. S1). The counts per miRNA were then timely and cost-effective fashion. Along with unprecedented normalized with respect to the total library size, and all scale and linear range, the sequencing approach does not have thermodynamic biases usually associated with hybridization- based qPCR and microarray platforms (Shendure 2008). More- over, deep-sequencing allows discovery of previously unchar- acterized miRNA species. These include not only new miRNAs that eluded previous cloning but also a panoply of miRNA variants, also known as isomiRs (Landgraf et al. 2007; Morin et al. 2008; Fernandez-Valverde et al. 2010). In this study the power of deep-sequencing technology was adopted to determine the complete repertoire of miRNAs expressed by human endothelial cells and to identify their modulations induced by low-oxygen tension. RESULTS Small RNA library generation, sequencing, and primary analysis Human umbilical vein endothelial cells (HUVECs) were cultured either in normoxic conditions or at 1% oxygen for 24 h. Then, small RNA libraries were prepared and deep- sequenced by using a SOLiD3 Plus Sequencer platform. As customary for these kinds of studies, two independent FIGURE 1. Analysis pipeline of annotated miRNAs. Flowchart de- scribing data analysis pipeline designed for the identification of populations were analyzed, each both in normoxia and in mature and star sequences of annotated miRNAs, their isoforms, hypoxia. The corresponding four cDNA libraries yielded a and novel miRNA star sequences. The number of remaining reads/ total of 114 million sequenced reads, and 78% (89 million, hits/miRNAs after each filtering step, averaged across all four libraries, an average of 22.2 million per library) of these reads passed is indicated in italics. The high-stringency filter of alignment (step 5) consisted in minimum alignment length of 17 nt, with a maximum of the quality control filter (quality value [QV] $ 9). Reads one mismatch and three sequencing errors. Reads with sequencing mapping across genomic sequence by >25 nt were removed errors in the first 15 nt or mismatches in the seed region were excluded. 2 RNA, Vol. 18, No. 3 Downloaded from rnajournal.cshlp.org on October 8, 2021 - Published by Cold Spring Harbor Laboratory Press microRNA
Recommended publications
  • Next Generation DNA Sequencing
    Genes 2010, 1, 385-387; doi:10.3390/genes1030385 OPEN ACCESS genes ISSN 2073-4425 www.mdpi.com/journal/genes Editorial Special Issue: Next Generation DNA Sequencing Paul Richardson PRC LLC, 25 Amber Lane, Lafayette, CA 94549, USA; E-Mail: [email protected]; Tel.: +1 925 301-5460 Received: 21 October 2010 / Accepted: 26 October 2010 / Published: 27 October 2010 Next Generation Sequencing (NGS) refers to technologies that do not rely on traditional dideoxy-nucleotide (Sanger) sequencing where labeled DNA fragments are physically resolved by electrophoresis. These new technologies rely on different strategies, but essentially all of them make use of real-time data collection of a base level incorporation event across a massive number of reactions (on the order of millions versus 96 for capillary electrophoresis for instance). The major commercial NGS platforms available to researchers are the 454 Genome Sequencer (Roche), Illumina (formerly Solexa) Genome analyzer, the SOLiD system (Applied Biosystems/Life Technologies) and the Heliscope (Helicos Corporation). The techniques and different strategies utilized by these platforms are reviewed in a number of the papers in this special issue. These technologies are enabling new applications that take advantage of the massive data produced by this next generation of sequencing instruments. In this special issue, nine papers review and demonstrate the utility and potential of next generation sequencing. One of the biggest consequences with NGS technologies is how to deal with all of the data produced by these platforms. Magi et al. [1] review the software tools available for the multiple functions needed to process and interpret the huge amounts of data produced by these instruments.
    [Show full text]
  • Evaluation of Normalization Methods in Mammalian Microrna-Seq Data
    Downloaded from rnajournal.cshlp.org on September 27, 2021 - Published by Cold Spring Harbor Laboratory Press METHOD Evaluation of normalization methods in mammalian microRNA-Seq data LANA XIA GARMIRE1 and SHANKAR SUBRAMANIAM1 Department of Bioengineering, Jacobs School of Engineering, University of California at San Diego, La Jolla, California 92093-0412, USA ABSTRACT Simple total tag count normalization is inadequate for microRNA sequencing data generated from the next generation sequencing technology. However, so far systematic evaluation of normalization methods on microRNA sequencing data is lacking. We comprehensively evaluate seven commonly used normalization methods including global normalization, Lowess normalization, Trimmed Mean Method (TMM), quantile normalization, scaling normalization, variance stabilization, and invariant method. We assess these methods on two individual experimental data sets with the empirical statistical metrics of mean square error (MSE) and Kolmogorov-Smirnov (K-S) statistic. Additionally, we evaluate the methods with results from quantitative PCR validation. Our results consistently show that Lowess normalization and quantile normalization perform the best, whereas TMM, a method applied to the RNA-Sequencing normalization, performs the worst. The poor performance of TMM normalization is further evidenced by abnormal results from the test of differential expression (DE) of microRNA-Seq data. Comparing with the models used for DE, the choice of normalization method is the primary factor that affects the results of DE. In summary, Lowess normalization and quantile normalization are recommended for normalizing microRNA-Seq data, whereas the TMM method should be used with caution. Keywords: microRNA-Seq; next generation sequencing; statistical normalization; high-throughput data analysis INTRODUCTION much smaller than mRNAs. For example, so far there are <1000 annotated microRNAs in human that are expected to The next generation sequencing (NGS) technology has regulate z30% of genes.
    [Show full text]
  • Ancient Microrna Profiles of a 14,300-Year-Old Canid Samples Confirm
    Downloaded from rnajournal.cshlp.org on October 7, 2021 - Published by Cold Spring Harbor Laboratory Press Fromm et al 1 2 Title: Ancient microRNA profiles of a 14,300-year-old canid samples confirm 3 taxonomic origin and give glimpses into tissue-specific gene regulation from 4 the Pleistocene 5 6 Running title: Sequencing of ancient microRNAs. 7 8 Authors: Bastian Fromm1*, Marcel Tarbier1, Oliver Smith2,3, Emilio Mármol- 9 Sánchez1,4, Love Dalén4,5, M. Tom P. Gilbert2,6, Marc R. Friedländer1* 10 11 Affiliations: 12 1 Science for Life Laboratory, Department of Molecular Biosciences, The Wenner-Gren 13 Institute, Stockholm University, Stockholm, Sweden 14 2 Section for Evolutionary Genomics, The Globe Institute, Faculty of Health and Medical 15 Sciences, University of Copenhagen, Copenhagen, Denmark 16 3 Micropathology Ltd, Venture Centre, University of Warwick Science Park, Coventry, UK 17 4 Centre for Palaeogenetics, Stockholm, Sweden 18 5 Department of Bioinformatics and Genetics, Swedish Museum of Natural History, 19 SE-10405 Stockholm, Sweden 20 6 Norwegian University of Science and Technology, University Museum, Trondheim, 21 Norway 22 23 24 * corresponding authors: [email protected] and [email protected] 25 26 Keywords: microRNAs, paleotranscriptomics, Pleistocene, canid, gene regulation Downloaded from rnajournal.cshlp.org on October 7, 2021 - Published by Cold Spring Harbor Laboratory Press Fromm et al. 27 Abstract 28 DNA sequencing is the current key technology for historic or ancient biological samples and has 29 led to many exciting discoveries in the field of paleogenomics. However, functional insights into 30 tissue identity, cellular composition or gene regulation cannot be gained from DNA.
    [Show full text]
  • Downregulation of Hsa-Microrna-204-5P and Identification of Its Potential Regulatory Network in Non-Small Cell Lung Cancer: RT-Q
    Liang et al. Respiratory Research (2020) 21:60 https://doi.org/10.1186/s12931-020-1274-9 RESEARCH Open Access Downregulation of hsa-microRNA-204-5p and identification of its potential regulatory network in non-small cell lung cancer: RT- qPCR, bioinformatic- and meta-analyses Chang-Yu Liang1†, Zu-Yun Li1†, Ting-Qing Gan2†, Ye-Ying Fang3†, Bin-Liang Gan1†, Wen-Jie Chen1†, Yi-Wu Dang1†, Ke Shi1†, Zhen-Bo Feng1* and Gang Chen1* Abstract Background: Pulmonary malignant neoplasms have a high worldwide morbidity and mortality, so the study of these malignancies using microRNAs (miRNAs) has attracted great interest and enthusiasm. The aim of this study was to determine the clinical effect of hsa-microRNA-204-5p (miR-204-5p) and its underlying molecular mechanisms in non-small cell lung cancer (NSCLC). Methods: Expression of miR-204-5p was investigated by real-time quantitative PCR (RT-qPCR). After data mining from public online repositories, several integrative assessment methods, including receiver operating characteristic (ROC) curves, hazard ratios (HR) with 95% confidence intervals (95% CI), and comprehensive meta-analyses, were conducted to explore the expression and clinical utility of miR-204-5p. The potential objects regulated and controlled by miR-204-5p in the course of NSCLC were identified by estimated target prediction and analysis. The regulatory network of miR-204-5p, with its target genes and transcription factors (TFs), was structured from database evidence and literature references. Results: The expression of miR-204-5p was downregulated in NSCLC, and the downtrend was related to gender, histological type, vascular invasion, tumor size, clinicopathologic grade and lymph node metastasis (P<0.05).
    [Show full text]
  • Analyzing Mirna from Raw NGS Data to Quantification
    G C A T T A C G G C A T genes Article DIANA-mAP: Analyzing miRNA from Raw NGS Data to Quantification Athanasios Alexiou 1,2, Dimitrios Zisis 2, Ioannis Kavakiotis 1, Marios Miliotis 1,2, Antonis Koussounadis 3 , Dimitra Karagkouni 1,2 and Artemis G. Hatzigeorgiou 1,2,3,* 1 DIANA Lab, Department of Computer Science and Biomedical Informatics, University of Thessaly, 35131 Lamia, Greece; [email protected] (A.A.); [email protected] (I.K.); [email protected] (M.M.); [email protected] (D.K.) 2 Hellenic Pasteur Institute, 11521 Athens, Greece; [email protected] 3 Department of Electrical & Computer Engineering, University of Thessaly, 38221 Volos, Greece; [email protected] * Correspondence: [email protected]; Tel.: +30-24210-74758; Fax: +30-24210-74997 Abstract: microRNAs (miRNAs) are small non-coding RNAs (~22 nts) that are considered central post-transcriptional regulators of gene expression and key components in many pathological con- ditions. Next-Generation Sequencing (NGS) technologies have led to inexpensive, massive data production, revolutionizing every research aspect in the fields of biology and medicine. Particularly, small RNA-Seq (sRNA-Seq) enables small non-coding RNA quantification on a high-throughput scale, providing a closer look into the expression profiles of these crucial regulators within the cell. Here, we present DIANA-microRNA-Analysis-Pipeline (DIANA-mAP), a fully automated compu- tational pipeline that allows the user to perform miRNA NGS data analysis from raw sRNA-Seq libraries to quantification and Differential Expression Analysis in an easy, scalable, efficient, and intuitive way. Emphasis has been given to data pre-processing, an early, critical step in the anal- ysis for the robustness of the final results and conclusions.
    [Show full text]
  • Addressing Bias in Small RNA Library Preparation for Sequencing: a New Protocol Recovers Micrornas That Evade Capture by Current Methods
    ORIGINAL RESEARCH published: 22 December 2015 doi: 10.3389/fgene.2015.00352 Addressing Bias in Small RNA Library Preparation for Sequencing: A New Protocol Recovers MicroRNAs that Evade Capture by Current Methods Jeanette Baran-Gale1,2*,C.LisaKurtz2, Michael R. Erdos3, Christina Sison4, Alice Young4, Emily E. Fannin2, Peter S. Chines3 and Praveen Sethupathy1,2* 1 Bioinformatics and Computational Biology Curriculum, School of Medicine, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA, 2 Department of Genetics, School of Medicine, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA, 3 National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA, 4 NIH Intramural Sequencing Center, National Institutes of Health, Rockville, MD, USA Edited by: Ghanshyam Upadhyay, Recent advances in sequencing technology have helped unveil the unexpected City College of New York-CUNY, USA Reviewed by: complexity and diversity of small RNAs. A critical step in small RNA library preparation Sandeep Kumar, for sequencing is the ligation of adapter sequences to both the 5 and 3 ends of small State University of New York College RNAs. Studies have shown that adapter ligation introduces a significant but widely of Optometry, USA Ravi C. Kalathur, unappreciated bias in the results of high-throughput small RNA sequencing. We show New York Structural Biology Center, that due to this bias the two widely used Illumina library preparation protocols produce USA strikingly different microRNA (miRNA) expression profiles in the same batch of cells. Manish Mishra, All Saints University College There are 102 highly expressed miRNAs that are >5-fold differentially detected and of Medicine, Saint Vincent some miRNAs, such as miR-24-3p, are over 30-fold differentially detected.
    [Show full text]
  • Microrna and Circrna Expression Analysis in a Zbtb1 Gene Knockout Monoclonal EL4 Cell Line
    ORIGINAL RESEARCH published: 05 July 2021 doi: 10.3389/fcimb.2021.706919 MicroRNA and circRNA Expression Analysis in a Zbtb1 Gene Knockout Monoclonal EL4 Cell Line † † Jun-Hong Wang 1,2,3 , Chun-Wei Shi 1,2,3 , Yi-Yuan Lu 1,2,3, Yan Zeng 1,2,3, Ming-Yang Cheng 1,2,3, Ru-Yu Wang 1,2,3,YuSun 1,2,3, Yan-Long Jiang 1,2,3, Wen-Tao Yang 1,2,3, Dan-Dan Zhao 1,2,3, Hai-Bin Huang 1,2,3, Li-Ping Ye 1,2,3, Xin Cao 1,2,3*, Gui-Lian Yang 1,2,3* and Chun-Feng Wang 1,2,3* 1 College of Veterinary Medicine, Jilin Agricultural University, Changchun, China, 2 Jilin Provincial Engineering Research Center of Animal Probiotics, Jilin Agricultural University, Changchun, China, 3 Key Laboratory of Animal Production and Product Quality Safety of Ministry of Education, Jilin Agricultural University, Changchun, China Edited by: Yongsheng Liu, Chinese Academy of Agricultural Zinc finger and BTB domain containing 1(Zbtb1) is a transcriptional suppressor protein, Sciences, China and a member of the mammalian Zbtb gene family. Previous studies have shown that Reviewed by: Zbtb1 is essential for T-cell development. However, the role of Zbtb1 in T-cell lymphoma is Peng Tan, Broad Institute, United States undetermined. In this study, an EL4 cell line with Zbtb1 deletion was constructed using the Dongqin Yang, CRISPR-Cas9 technique. The expression profiles of microRNA and circRNA produced by Fudan University, China the control and gene deletion groups were determined by RNA-seq.
    [Show full text]
  • Single-Cell Microrna Sequencing Method Comparison and Application to Cell Lines and Circulating Lung Tumor Cells
    ARTICLE https://doi.org/10.1038/s41467-021-24611-w OPEN Single-cell microRNA sequencing method comparison and application to cell lines and circulating lung tumor cells Sarah M. Hücker 1,8, Tobias Fehlmann 2,8, Christian Werno1, Kathrin Weidele1, Florian Lüke1,3, ✉ Anke Schlenska-Lange4, Christoph A. Klein 1,5, Andreas Keller 2,6,7,9 & Stefan Kirsch 1,9 1234567890():,; Molecular single cell analyses provide insights into physiological and pathological processes. Here, in a stepwise approach, we first evaluate 19 protocols for single cell small RNA sequencing on MCF7 cells spiked with 1 pg of 1,006 miRNAs. Second, we analyze MCF7 single cell equivalents of the eight best protocols. Third, we sequence single cells from eight different cell lines and 67 circulating tumor cells (CTCs) from seven SCLC patients. Altogether, we analyze 244 different samples. We observe high reproducibility within pro- tocols and reads covered a broad spectrum of RNAs. For the 67 CTCs, we detect a median of 68 miRNAs, with 10 miRNAs being expressed in 90% of tested cells. Enrichment analysis suggested the lung as the most likely organ of origin and enrichment of cancer-related categories. Even the identification of non-annotated candidate miRNAs was feasible, underlining the potential of single cell small RNA sequencing. 1 Division of Personalized Tumor Therapy, Fraunhofer-Institute for Toxicology and Experimental Medicine, Regensburg, Germany. 2 Chair for Clinical Bioinformatics, Saarland University, Saarbrücken, Germany. 3 Department of Internal Medicine III Hematology & Oncology, University Medical Center Regensburg, Regensburg, Germany. 4 Department of Oncology and Hematology, Hospital Barmherzige Brüder, Regensburg, Germany.
    [Show full text]
  • Erythrocyte Microrna Sequencing Reveals Differential Expression in Relapsing-Remitting Multiple Sclerosis Kira Groen1,2, Vicki E
    Groen et al. BMC Medical Genomics (2018) 11:48 https://doi.org/10.1186/s12920-018-0365-7 RESEARCH ARTICLE Open Access Erythrocyte microRNA sequencing reveals differential expression in relapsing-remitting multiple sclerosis Kira Groen1,2, Vicki E. Maltby1,2, Rodney A. Lea2,3, Katherine A. Sanders4, J. Lynn Fink5, Rodney J. Scott2,6,7, Lotti Tajouri8 and Jeannette Lechner-Scott1,2,9* Abstract Background: There is a paucity of knowledge concerning erythrocytes in the aetiology of Multiple Sclerosis (MS) despite their potential to contribute to disease through impaired antioxidant capacity and altered haemorheological features. Several studies have identified an abundance of erythrocyte miRNAs and variable profiles associated with disease states, such as sickle cell disease and malaria. The aim of this study was to compare the erythrocyte miRNA profile of relapsing-remitting MS (RRMS) patients to healthy sex- and age-matched controls. Methods: Erythrocytes were purified by density-gradient centrifugation and RNA was extracted. Following library preparation, samples were run on a HiSeq4000 Illumina instrument (paired-end 100 bp sequencing). Sequenced erythrocyte miRNA profiles (9 patients and 9 controls) were analysed by DESeq2. Differentially expressed miRNAs were validated by RT-qPCR using miR-152-3p as an endogenous control and replicated in a larger cohort (20 patients and 18 controls). After logarithmic transformation, differential expression was determined by two-tailed unpaired t-tests. Logistic regression analysis was carried out and receiver operating characteristic (ROC) curves were generated to determine biomarker potential. Results: A total of 236 erythrocyte miRNAs were identified. Of twelve differentially expressed miRNAs in RRMS two showed increased expression (adj.
    [Show full text]
  • Computational Analysis of Microrna Function
    microRNA and smallRNA profiling Anton Enright Group Leader EMBL - European Bioinformatics Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SD, UK http://www.ebi.ac.uk/enright/ [email protected] High Throughput Sequencing Analysis 2013 Thursday, 24 October 13 microRNA Regulation Adapted from He L. and Hannon G.J; Nature reviews genetics 5:552 (2004) EMBL-EBI 2 High throughput RNA Seq Thursday, 24 October 13 microRNA Regulation Adapted from He L. and Hannon G.J; Nature reviews genetics 5:552 (2004) EMBL-EBI 2 High throughput RNA Seq Thursday, 24 October 13 MicroRNA Processing EMBL-EBI 3 High throughput RNA Seq Thursday, 24 October 13 MicroRNA Processing EMBL-EBI 3 High throughput RNA Seq Thursday, 24 October 13 MicroRNA Processing EMBL-EBI 3 High throughput RNA Seq Thursday, 24 October 13 microRNA Profiling EMBL-EBI 4 High throughput RNA Seq Thursday, 24 October 13 microRNA Profiling • MicroRNA microarray • Statistical issues with number of probes • Cross-hybridisation and biases • Cheap and straightforward • Small RNA Sequencing • Quantitation can be difficult due to biases • Ability to detect novel microRNAs, edits and variation • Reproducible • qRT-PCR Based Approaches • quite accurate • lower throughput • Single Molecule profiling • e.g. Nanostring • New technology EMBL-EBI 5 High throughput RNA Seq Thursday, 24 October 13 Nanostring EMBL-EBI 6 High throughput RNA Seq Thursday, 24 October 13 Nanostring EMBL-EBI 6 High throughput RNA Seq Thursday, 24 October 13 microRNA microarrays • MicroRNA array manufacturers love to talk about how many probes they have on their chips. • In reality there are < 1000 microRNAs for your species and likely <30 expressed in your sample.
    [Show full text]
  • Mirge 2.0 for Comprehensive Analysis of Microrna Sequencing Data Yin Lu, Alexander S
    Lu et al. BMC Bioinformatics (2018) 19:275 https://doi.org/10.1186/s12859-018-2287-y SOFTWARE Open Access miRge 2.0 for comprehensive analysis of microRNA sequencing data Yin Lu, Alexander S. Baras and Marc K. Halushka* Abstract 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 R2 = 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.
    [Show full text]
  • Evaluating Gastroenteropancreatic Neuroendocrine Tumors Through Microrna Sequencing
    26 1 Endocrine-Related N Panarelli, K Tyryshkin miRNA-based evaluation of 26:1 47–57 Cancer et al. GEP-NETs RESEARCH Evaluating gastroenteropancreatic neuroendocrine tumors through microRNA sequencing Nicole Panarelli1,*,†, Kathrin Tyryshkin2,*, Justin Jong Mun Wong2, Adrianna Majewski2, Xiaojing Yang2, Theresa Scognamiglio1, Michelle Kang Kim3, Kimberly Bogardus4, Thomas Tuschl4, Yao-Tseng Chen1 and Neil Renwick2,4 1Department of Pathology and Laboratory Medicine, Weill Cornell Medicine, New York, New York, USA 2Laboratory of Translational RNA Biology, Department of Pathology and Molecular Medicine, Queen’s University, Kingston, Ontario, Canada 3Center for Carcinoid and Neuroendocrine Tumors of Mount Sinai, Icahn School of Medicine at Mount Sinai, New York, New York, USA 4HHMI, Laboratory of RNA Molecular Biology, The Rockefeller University, New York, New York, USA Correspondence should be addressed to N Renwick: [email protected] *(N Panarelli and K Tyryshkin contributed equally to this work) †(N Panarelli is now at Department of Pathology Albert Einstein College of Medicine, Bronx, New York, USA) Abstract Gastroenteropancreatic neuroendocrine tumors (GEP-NETs) can be challenging to evaluate Key Words histologically. MicroRNAs (miRNAs) are small RNA molecules that often are excellent f gastroenteropancreatic biomarkers due to their abundance, cell-type and disease stage specificity and stability. To neuroendocrine tumors evaluate miRNAs as adjunct tissue markers for classifying and grading well-differentiated f classification GEP-NETs, we generated and compared miRNA expression profiles from four pathological f biomarkers types of GEP-NETs. Using quantitative barcoded small RNA sequencing and state-of-the- f microRNA art sequence annotation, we generated comprehensive miRNA expression profiles from f small RNA sequencing archived pancreatic, ileal, appendiceal and rectal NETs.
    [Show full text]