MicroRNA-125a is over-expressed in insulin target tissues in a spontaneous rat model of Type 2 Diabetes

Blanca M Herrera, Helen E Lockstone, Jennifer M Taylor, Quin F Wills, Pamela J Kaisaki, Amy Barrett, Carme Camps, Christina Fernandez, Jiannis Ragoussis, Dominique Gauguier, Mark I McCarthy, Cecilia M Lindgren, Blanca M Herrera, Helen E Lockstone, Jennifer M Taylor, Quin F Wills, Pamela J Kaisaki, Amy Barrett, Carme Camps, Christina Fernandez, Jiannis Ragoussis, Dominique Gauguier, Mark I McCarthy, Cecilia M Lindgren

Abstract

Background: MicroRNAs (miRNAs) are non-coding RNA molecules involved in post-transcriptional control of gene expression of a wide number of genes, including those involved in glucose homeostasis. Type 2 diabetes (T2D) is characterized by hyperglycaemia and defects in insulin secretion and action at target tissues. We sought to establish differences in global miRNA expression in two insulin-target tissues from inbred rats of spontaneously diabetic and normoglycaemic strains.

Methods: We used a miRNA microarray platform to measure global miRNA expression in two insulin-target tissues: liver and adipose tissue from inbred rats of spontaneously diabetic (Goto-Kakizaki [GK]) and normoglycaemic (Brown-Norway [BN]) strains which are extensively used in genetic studies of T2D. MiRNA data were integrated with gene expression data from the same rats to investigate how differentially expressed miRNAs affect the expression of predicted target gene transcripts.

Results: The expression of 170 miRNAs was measured in liver and adipose tissue of GK and BN rats. Based on a p-value for differential expression between GK and BN, the most significant change in expression was observed for miR-125a in liver (FC = 5.61, P = 0.001, Padjusted = 0.10); this overexpression was validated using quantitative RT-PCR (FC = 13.15, P = 0.0005). MiR-125a also showed over-expression in the GK vs. BN analysis within adipose tissue (FC = 1.97, P = 0.078, Padjusted = 0.99), as did the previously reported miR-29a (FC = 1.51, P = 0.05, Padjusted = 0.99). In-silico tools assessing the biological role of predicted miR-125a target genes suggest an over-representation of genes involved in the MAPK signaling pathway. Gene expression analysis identified 1308 genes with significantly different expression between GK and BN rats (Padjusted < 0.05): 233 in liver and 1075 in adipose tissue. Pathways related to glucose and lipid metabolism were significantly over-represented among these genes. Enrichment analysis suggested that differentially expressed genes in GK compared to BN included more predicted miR-125a target genes than would be expected by chance in adipose tissue (FDR = 0.006 for up-regulated genes; FDR = 0.036 for down-regulated genes) but not in liver (FDR = 0.074 for up-regulated genes; FDR = 0.248 for down-regulated genes).

Conclusion: MiR-125a is over-expressed in liver in hyperglycaemic GK rats relative to normoglycaemic BN rats, and our array data also suggest miR-125a is over-expressed in adipose tissue. We demonstrate the use of in-silico tools to provide the basis for further investigation of the potential role of miR-125a in T2D. In particular, the enrichment of predicted miR-125a target genes among differentially expressed genes has identified likely target genes and indicates that integrating global miRNA and mRNA expression data may give further insights into miRNA-mediated regulation of gene expression.

Figures

Figure 1
Figure 1
Relative expression of miR-125a in liver from diabetic GK rats compared to BN rats. * P = 0.0005 (miR-125a expression normalized against snoRNA and 4.5S).
Figure 2
Figure 2
Gene expression analysis of adipose tissue and liver reveals differential expression between hyperglycaemic and normoglycaemic rats. The number of genes that are significant (adjusted p < 0.05) in the comparison of 4 hyperglycaemic GK rats and 4 normoglycaemic BN rats in adipose tissue (477 up; 598 down), liver (95 up; 138 down), and in both tissues (41 up; 55 down). Five genes significant in both tissues showed opposite directions of change and are not included in the adipose tissue and liver group.

References

    1. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–297. doi: 10.1016/S0092-8674(04)00045-5.
    1. Vasudevan S, Tong Y, Steitz JA. Switching from repression to activation: microRNAs can up-regulate translation. Science. 2007;318:1931–1934. doi: 10.1126/science.1149460.
    1. Hutvagner G, Zamore PD. A microRNA in a multiple-turnover RNAi enzyme complex. Science. 2002;297:2056–2060. doi: 10.1126/science.1073827.
    1. Dugas DV, Bartel B. MicroRNA regulation of gene expression in plants. Curr Opin Plant Biol. 2004;7:512–520. doi: 10.1016/j.pbi.2004.07.011.
    1. Lagos-Quintana M, Rauhut R, Lendeckel W, Tuschl T. Identification of novel genes coding for small expressed RNAs. Science. 2001;294:853–858. doi: 10.1126/science.1064921.
    1. Lim LP, Lau NC, Garrett-Engele P, Grimson A, Schelter JM, Castle J, Bartel DP, Linsley PS, Johnson JM. Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature. 2005;433:769–773. doi: 10.1038/nature03315.
    1. Selbach M, Schwanhausser B, Thierfelder N, Fang Z, Khanin R, Rajewsky N. Widespread changes in protein synthesis induced by microRNAs. Nature. 2008;455:58–63. doi: 10.1038/nature07228.
    1. Baek D, Villen J, Shin C, Camargo FD, Gygi SP, Bartel DP. The impact of microRNAs on protein output. Nature. 2008;455:64–71. doi: 10.1038/nature07242.
    1. Krutzfeldt J, Rajewsky N, Braich R, Rajeev KG, Tuschl T, Manoharan M, Stoffel M. Silencing of microRNAs in vivo with 'antagomirs'. Nature. 2005;438:685–689. doi: 10.1038/nature04303.
    1. Sood P, Krek A, Zavolan M, Macino G, Rajewsky N. Cell-type-specific signatures of microRNAs on target mRNA expression. Proc Natl Acad Sci USA. 2006;103:2746–2751. doi: 10.1073/pnas.0511045103.
    1. Wang X. Systematic identification of microRNA functions by combining target prediction and expression profiling. Nucleic Acids Res. 2006;34:1646–1652. doi: 10.1093/nar/gkl068.
    1. Cheng C, Li LM. Inferring microRNA activities by combining gene expression with microRNA target prediction. PLoS ONE. 2008;3:e1989. doi: 10.1371/journal.pone.0001989.
    1. Ambros V. The functions of animal microRNAs. Nature. 2004;431:350–355. doi: 10.1038/nature02871.
    1. Poy MN, Spranger M, Stoffel M. microRNAs and the regulation of glucose and lipid metabolism. Diabetes Obes Metab. 2007;9:67–73. doi: 10.1111/j.1463-1326.2007.00775.x.
    1. Tang X, Tang G, Ozcan S. Role of microRNAs in diabetes. Role of microRNAs in diabetes. Biochim Biophys Acta. 2008;1779:697–701.
    1. Poy MN, Eliasson L, Krutzfeldt J, Kuwajima S, Ma X, Macdonald PE, Pfeffer S, Tuschl T, Rajewsky N, Rorsman P, Stoffel M. A pancreatic islet-specific microRNA regulates insulin secretion. Nature. 2004;432:226–230. doi: 10.1038/nature03076.
    1. Plaisance V, Abderrahmani A, Perret-Menoud V, Jacquemin P, Lemaigre F, Regazzi R. MicroRNA-9 controls the expression of Granuphilin/Slp4 and the secretory response of insulin-producing cells. J Biol Chem. 2006;281:26932–26942. doi: 10.1074/jbc.M601225200.
    1. Baroukh N, Ravier MA, Loder MK, Hill EV, Bounacer A, Scharfmann R, Rutter GA, Van Obberghen E. MicroRNA-124a regulates Foxa2 expression and intracellular signaling in pancreatic beta-cell lines. J Biol Chem. 2007;282:19575–19588. doi: 10.1074/jbc.M611841200.
    1. He A, Zhu L, Gupta N, Chang Y, Fang F. Overexpression of micro ribonucleic acid 29, highly up-regulated in diabetic rats, leads to insulin resistance in 3T3-L1 adipocytes. Mol Endocrinol. 2007;21:2785–2794. doi: 10.1210/me.2007-0167.
    1. Picarel-Blanchot F, Berthelier C, Bailbe D, Portha B. Impaired insulin secretion and excessive hepatic glucose production are both early events in the diabetic GK rat. Am J Physiol. 1996;271:E755–762.
    1. Gauguier D. The rat as a model physiological system. Vol. 3. London: Wiley; 2006.
    1. Solberg LC, Valdar W, Gauguier D, Nunez G, Taylor A, Burnett S, Arboledas-Hita C, Hernandez-Pliego P, Davidson S, Burns P, Bhattacharya S, Hough T, Higgs D, Klenerman P, Cookson WO, Zhang Y, Deacon RM, Rawlins JN, Mott R, Flint J. A protocol for high-throughput phenotyping, suitable for quantitative trait analysis in mice. Mamm Genome. 2006;17:129–146. doi: 10.1007/s00335-005-0112-1.
    1. R: A Language and Environment for Statistical Computing
    1. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5:R80. doi: 10.1186/gb-2004-5-10-r80.
    1. Smyth GK, Speed T. Normalization of cDNA microarray data. Methods. 2003;31:265–273. doi: 10.1016/S1046-2023(03)00155-5.
    1. Smyth GK. Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3
    1. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B. 1995;57:289–300.
    1. Huber W, von Heydebreck A, Sultmann H, Poustka A, Vingron M. Variance stabilization applied to microarray data calibration and to the quantification of differential expression. Bioinformatics. 2002;18:S96–104.
    1. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25:402–408. doi: 10.1006/meth.2001.1262.
    1. Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB. Prediction of mammalian microRNA targets. Cell. 2003;115:787–798. doi: 10.1016/S0092-8674(03)01018-3.
    1. Krek A, Grun D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, MacMenamin P, da Piedade I, Gunsalus KC, Stoffel M, Rajewsky N. Combinatorial microRNA target predictions. Nat Genet. 2005;37:495–500. doi: 10.1038/ng1536.
    1. John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS. Human MicroRNA targets. PLoS Biol. 2004;2:e363. doi: 10.1371/journal.pbio.0020363.
    1. John B, Sander C, Marks DS. Prediction of human microRNA targets. Methods Mol Biol. 2006;342:101–113.
    1. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30. doi: 10.1093/nar/28.1.27.
    1. Carmona-Saez P, Chagoyen M, Tirado F, Carazo JM, Pascual-Montano A. GENECODIS: a web-based tool for finding significant concurrent annotations in gene lists. Genome Biol. 2007;8:R3. doi: 10.1186/gb-2007-8-1-r3.
    1. Hochberg Y, Benjamini Y. More powerful procedures for multiple significance testing. Stat Med. 1990;9:811–818. doi: 10.1002/sim.4780090710.
    1. Wilder SP, Kaisaki PJ, Argoud K, Salhan A, Ragoussis J, Bihoreau MT, Gauguier D. Comparative analysis of methods for gene transcription profiling data derived from different microarray technologies in rat and mouse models of diabetes. BMC Genomics. 2009;10:63. doi: 10.1186/1471-2164-10-63.
    1. Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP. MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell. 2007;27:91–105. doi: 10.1016/j.molcel.2007.06.017.
    1. Yekta S, Shih IH, Bartel DP. MicroRNA-directed cleavage of HOXB8 mRNA. Science. 2004;304:594–596. doi: 10.1126/science.1097434.
    1. Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS. MicroRNA targets in Drosophila. Genome Biol. 2003;5:R1. doi: 10.1186/gb-2003-5-1-r1.
    1. Scott GK, Goga A, Bhaumik D, Berger CE, Sullivan CS, Benz CC. Coordinate suppression of ERBB2 and ERBB3 by enforced expression of micro-RNA miR-125a or miR-125b. J Biol Chem. 2007;282:1479–1486. doi: 10.1074/jbc.M609383200.
    1. Wu L, Belasco JG. Micro-RNA regulation of the mammalian lin-28 gene during neuronal differentiation of embryonal carcinoma cells. Mol Cell Biol. 2005;25:9198–9208. doi: 10.1128/MCB.25.21.9198-9208.2005.
    1. Nam EJ, Yoon H, Kim SW, Kim H, Kim YT, Kim JH, Kim JW, Kim S. MicroRNA expression profiles in serous ovarian carcinoma. Clin Cancer Res. 2008;14:2690–2695. doi: 10.1158/1078-0432.CCR-07-1731.
    1. Mattie MD, Benz CC, Bowers J, Sensinger K, Wong L, Scott GK, Fedele V, Ginzinger D, Getts R, Haqq C. Optimized high-throughput microRNA expression profiling provides novel biomarker assessment of clinical prostate and breast cancer biopsies. Mol Cancer. 2006;5:24. doi: 10.1186/1476-4598-5-24.
    1. Fujishiro M, Gotoh Y, Katagiri H, Sakoda H, Ogihara T, Anai M, Onishi Y, Ono H, Abe M, Shojima N, Fukushima Y, Kikuchi M, Oka Y, Asano T. Three mitogen-activated protein kinases inhibit insulin signaling by different mechanisms in 3T3-L1 adipocytes. Mol Endocrinol. 2003;17:487–497. doi: 10.1210/me.2002-0131.
    1. Engelman JA, Berg AH, Lewis RY, Lisanti MP, Scherer PE. Tumor necrosis factor alpha-mediated insulin resistance, but not dedifferentiation, is abrogated by MEK1/2 inhibitors in 3T3-L1 adipocytes. Mol Endocrinol. 2000;14:1557–1569. doi: 10.1210/me.14.10.1557.
    1. Thorleifsson G, Walters GB, Gudbjartsson DF, Steinthorsdottir V, Sulem P, Helgadottir A, Styrkarsdottir U, Gretarsdottir S, Thorlacius S, Jonsdottir I, Jonsdottir T, Olafsdottir EJ, Olafsdottir GH, Jonsson T, Jonsson F, Borch-Johnsen K, Hansen T, Andersen G, Jorgensen T, Lauritzen T, Aben KK, Verbeek AL, Roeleveld N, Kampman E, Yanek LR, Becker LC, Tryggvadottir L, Rafnar T, Becker DM, Gulcher J, Kiemeney LA, Pedersen O, Kong A, Thorsteinsdottir U, Stefansson K. Genome-wide association yields new sequence variants at seven loci that associate with measures of obesity. Nat Genet. 2009;41:18–24. doi: 10.1038/ng.274.
    1. Gauguier D, Froguel P, Parent V, Bernard C, Bihoreau MT, Portha B, James MR, Penicaud L, Lathrop M, Ktorza A. Chromosomal mapping of genetic loci associated with non-insulin dependent diabetes in the GK rat. Nat Genet. 1996;12:38–43. doi: 10.1038/ng0196-38.
    1. Dumas ME, Wilder SP, Bihoreau MT, Barton RH, Fearnside JF, Argoud K, D'Amato L, Wallis RH, Blancher C, Keun HC, Baunsgaard D, Scott J, Sidelmann UG, Nicholson JK, Gauguier D. Direct quantitative trait locus mapping of mammalian metabolic phenotypes in diabetic and normoglycemic rat models. Nat Genet. 2007;39:666–672. doi: 10.1038/ng2026.
    1. Wallis RH, Collins SC, Kaisaki PJ, Argoud K, Wilder SP, Wallace KJ, Ria M, Ktorza A, Rorsman P, Bihoreau MT, Gauguier D. Pathophysiological, genetic and gene expression features of a novel rodent model of the cardio-metabolic syndrome. PLoS ONE. 2008;3:e2962. doi: 10.1371/journal.pone.0002962.

Source: PubMed

3
구독하다