Extent, Causes, and Consequences of Small RNA Expression Variation in Human Adipose Tissue

English version České info

Small RNAs are functional molecules that modulate mRNA transcripts and have been implicated in the aetiology of several common diseases. However, little is known about the extent of their variability within the human population. Here, we characterise the extent, causes, and effects of naturally occurring variation in expression and sequence of small RNAs from adipose tissue in relation to genotype, gene expression, and metabolic traits in the MuTHER reference cohort. We profiled the expression of 15 to 30 base pair RNA molecules in subcutaneous adipose tissue from 131 individuals using high-throughput sequencing, and quantified levels of 591 microRNAs and small nucleolar RNAs. We identified three genetic variants and three RNA editing events. Highly expressed small RNAs are more conserved within mammals than average, as are those with highly variable expression. We identified 14 genetic loci significantly associated with nearby small RNA expression levels, seven of which also regulate an mRNA transcript level in the same region. In addition, these loci are enriched for variants significant in genome-wide association studies for body mass index. Contrary to expectation, we found no evidence for negative correlation between expression level of a microRNA and its target mRNAs. Trunk fat mass, body mass index, and fasting insulin were associated with more than twenty small RNA expression levels each, while fasting glucose had no significant associations. This study highlights the similar genetic complexity and shared genetic control of small RNA and mRNA transcripts, and gives a quantitative picture of small RNA expression variation in the human population.

Vyšlo v časopise: Extent, Causes, and Consequences of Small RNA Expression Variation in Human Adipose Tissue. PLoS Genet 8(5): e32767. doi:10.1371/journal.pgen.1002704
Kategorie: Research Article
prolekare.web.journal.doi_sk: https://doi.org/10.1371/journal.pgen.1002704

Souhrn

Zdroje

1. MattickJSMakuninIV 2006 Non-coding RNA. Hum Mol Genet 15 Spec No 1 R17 29

2. LauNCLimLPWeinsteinEGBartelDP 2001 An abundant class of tiny RNAs with probable regulatory roles in Caenorhabditis elegans. Science 294 858 862

3. Lagos-QuintanaMRauhutRLendeckelWTuschlT 2001 Identification of novel genes coding for small expressed RNAs. Science 294 853 858

4. HeLHannonGJ 2004 MicroRNAs: small RNAs with a big role in gene regulation. Nat Rev Genet 5 522 531

5. CzechBHannonGJ 2011 Small RNA sorting: matchmaking for Argonautes. Nat Rev Genet 12 19 31

6. LeeRCFeinbaumRLAmbrosV 1993 The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell 75 843 854

7. KozomaraAGriffiths-JonesS 2011 miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res 39 D152 157

8. WinterJJungSKellerSGregoryRIDiederichsS 2009 Many roads to maturity: microRNA biogenesis pathways and their regulation. Nat Cell Biol 11 228 234

9. KimVNHanJSiomiMC 2009 Biogenesis of small RNAs in animals. Nat Rev Mol Cell Biol 10 126 139

10. FangZRajewskyN 2011 The impact of miRNA target sites in coding sequences and in 3′UTRs. PLoS ONE 6 e18067 doi:10.1371/journal.pone.0018067

11. KissT 2002 Small nucleolar RNAs: an abundant group of noncoding RNAs with diverse cellular functions. Cell 109 145 148

12. BrameierMHerwigAReinhardtRWalterLGruberJ 2011 Human box C/D snoRNAs with miRNA like functions: expanding the range of regulatory RNAs. Nucleic Acids Res 39 675 686

13. MarzMGruberARHoner Zu SiederdissenCAmmanFBadeltS 2011 Animal snoRNAs and scaRNAs with exceptional structures. RNA Biol 8

14. TaftRJGlazovEALassmannTHayashizakiYCarninciP 2009 Small RNAs derived from snoRNAs. RNA 15 1233 1240

15. LukkMKapusheskyMNikkilaJParkinsonHGoncalvesA 2010 A global map of human gene expression. Nat Biotechnol 28 322 324

16. StrangerBENicaACForrestMSDimasABirdCP 2007 Population genomics of human gene expression. Nat Genet 39 1217 1224

17. EmilssonVThorleifssonGZhangBLeonardsonASZinkF 2008 Genetics of gene expression and its effect on disease. Nature 452 423 428

18. CheungVGNayakRRWangIXElwynSCousinsSM 2010 Polymorphic cis- and trans-regulation of human gene expression. PLoS Biol 8 e1000480 doi:10.1371/journal.pbio.1000480

19. NicaACPartsLGlassDNisbetJBarrettA 2011 The architecture of gene regulatory variation across multiple human tissues: the MuTHER study. PLoS Genet 7 e1002003 doi:10.1371/journal.pgen.1002003

20. SmallKSHedmanAKGrundbergENicaACThorleifssonG 2011 Identification of an imprinted master trans regulator at the KLF14 locus related to multiple metabolic phenotypes. Nat Genet 43 561 564

21. BorelCDeutschSLetourneauAMigliavaccaEMontgomerySB 2011 Identification of cis- and trans-regulatory variation modulating microRNA expression levels in human fibroblasts. Genome Res 21 68 73

22. RantalainenMHerreraBMNicholsonGBowdenRWillsQF 2011 MicroRNA Expression in Abdominal and Gluteal Adipose Tissue Is Associated with mRNA Expression Levels and Partly Genetically Driven. PLoS ONE 6 e27338 doi:10.1371/journal.pone.0027338

23. NicaACMontgomerySBDimasASStrangerBEBeazleyC 2010 Candidate causal regulatory effects by integration of expression QTLs with complex trait genetic associations. PLoS Genet 6 e1000895 doi:10.1371/journal.pgen.1000895

24. FehrmannRSJansenRCVeldinkJHWestraHJArendsD 2011 Trans-eQTLs reveal that independent genetic variants associated with a complex phenotype converge on intermediate genes, with a major role for the HLA. PLoS Genet 7 e1002197 doi:10.1371/journal.pgen.1002197

25. PartsLStegleOWinnJDurbinR 2011 Joint genetic analysis of gene expression data with inferred cellular phenotypes. PLoS Genet 7 e1001276 doi:10.1371/journal.pgen.1001276

26. TakanabeROnoKAbeYTakayaTHorieT 2008 Up-regulated expression of microRNA-143 in association with obesity in adipose tissue of mice fed high-fat diet. Biochem Biophys Res Commun 376 728 732

27. BrestPLapaquettePSouidiMLebrigandKCesaroA 2011 A synonymous variant in IRGM alters a binding site for miR-196 and causes deregulation of IRGM-dependent xenophagy in Crohn's disease. Nat Genet 43 242 245

28. GehrkeSImaiYSokolNLuB 2010 Pathogenic LRRK2 negatively regulates microRNA-mediated translational repression. Nature 466 637 641

29. KumarMSLuJMercerKLGolubTRJacksT 2007 Impaired microRNA processing enhances cellular transformation and tumorigenesis. Nat Genet 39 673 677

30. SahooTdel GaudioDGermanJRShinawiMPetersSU 2008 Prader-Willi phenotype caused by paternal deficiency for the HBII-85 C/D box small nucleolar RNA cluster. Nat Genet 40 719 721

31. EstellerM 2011 Non-coding RNAs in human disease. Nat Rev Genet 12 861 874

32. KlotingNBertholdSKovacsPSchonMRFasshauerM 2009 MicroRNA expression in human omental and subcutaneous adipose tissue. PLoS ONE 4 e4699 doi:10.1371/journal.pone.0004699

33. Grundberg ESKSHedmanÅKNicaACBuilAKeildsonSBellJT 2012 Mapping cis and trans regulatory effects across multiple tissues in twins: the MuTHER study. submitted

34. SpectorTDWilliamsFM 2006 The UK Adult Twin Registry (TwinsUK). Twin Res Hum Genet 9 899 906

35. FlicekPAmodeMRBarrellDBealKBrentS 2011 Ensembl 2011. Nucleic Acids Res 39 D800 806

36. Guerra-AssuncaoJAEnrightAJ 2010 MapMi: automated mapping of microRNA loci. BMC Bioinformatics 11 133

37. KimTKHembergMGrayJMCostaAMBearDM 2010 Widespread transcription at neuronal activity-regulated enhancers. Nature 465 182 187

38. EsauCKangXPeraltaEHansonEMarcussonEG 2004 MicroRNA-143 regulates adipocyte differentiation. J Biol Chem 279 52361 52365

39. van SolingenCSeghersLBijkerkRDuijsJMRoetenMK 2009 Antagomir-mediated silencing of endothelial cell specific microRNA-126 impairs ischemia-induced angiogenesis. J Cell Mol Med 13 1577 1585

40. LeeDYDengZWangCHYangBB 2007 MicroRNA-378 promotes cell survival, tumor growth, and angiogenesis by targeting SuFu and Fus-1 expression. Proc Natl Acad Sci U S A 104 20350 20355

41. WangQHuangZXueHJinCJuXL 2008 MicroRNA miR-24 inhibits erythropoiesis by targeting activin type I receptor ALK4. Blood 111 588 595

42. Bruchova-VotavovaHYoonDPrchalJT 2010 miR-451 enhances erythroid differentiation in K562 cells. Leuk Lymphoma 51 686 693

43. KuchenSReschWYamaneAKuoNLiZ 2010 Regulation of microRNA expression and abundance during lymphopoiesis. Immunity 32 828 839

44. HouJLinLZhouWWangZDingG 2011 Identification of miRNomes in human liver and hepatocellular carcinoma reveals miR-199a/b-3p as therapeutic target for hepatocellular carcinoma. Cancer Cell 19 232 243

45. StarkMSTyagiSNancarrowDJBoyleGMCookAL 2010 Characterization of the Melanoma miRNAome by Deep Sequencing. PLoS ONE 5 e9685 doi:10.1371/journal.pone.0009685

46. SchotteDMoqadamFALange-TurenhoutEAChenCvan IjckenWF 2011 Discovery of new microRNAs by small RNAome deep sequencing in childhood acute lymphoblastic leukemia. Leukemia 25 1389 1399

47. FaghihiMAZhangMHuangJModarresiFVan der BrugMP 2010 Evidence for natural antisense transcript-mediated inhibition of microRNA function. Genome Biol 11 R56

48. SaundersMALiangHLiWH 2007 Human polymorphism at microRNAs and microRNA target sites. Proc Natl Acad Sci U S A 104 3300 3305

49. WymanSKKnoufECParkinRKFritzBRLinDW 2011 Post-transcriptional generation of miRNA variants by multiple nucleotidyl transferases contributes to miRNA transcriptome complexity. Genome Res 21 1450 1461

50. FujitaPARheadBZweigASHinrichsASKarolchikD 2011 The UCSC Genome Browser database: update 2011. Nucleic Acids Res 39 D876 882

51. SchroederAMuellerOStockerSSalowskyRLeiberM 2006 The RIN: an RNA integrity number for assigning integrity values to RNA measurements. BMC Mol Biol 7 3

52. StegleOPartsLDurbinRWinnJ 2010 A Bayesian framework to account for complex non-genetic factors in gene expression levels greatly increases power in eQTL studies. PLoS Comput Biol 6 e1000770 doi:10.1371/journal.pcbi.1000770

53. LeekJTStoreyJD 2007 Capturing heterogeneity in gene expression studies by surrogate variable analysis. PLoS Genet 3 e161 doi:10.1371/journal.pgen.0030161

54. StegleOPartsLPiipariMWinnJDurbinR 2012 Using probabilistic estimation of expression residuals (PEER) to obtain increased power and interpretability of gene expression analyses. Nat Protoc 7 500 507

55. FuJWolfsMGDeelenPWestraHJFehrmannRS 2012 Unraveling the regulatory mechanisms underlying tissue-dependent genetic variation of gene expression. PLoS Genet 8 e1002431 doi:10.1371/journal.pgen.1002431

56. ReimandJArakTViloJ 2011 g:Profiler–a web server for functional interpretation of gene lists (2011 update). Nucleic Acids Res 39 W307 315

57. KidaKNakajimaMMohriTOdaYTakagiS 2011 PPARalpha is regulated by miR-21 and miR-27b in human liver. Pharm Res 28 2467 2476

58. MadhyasthaRMadhyasthaHNakajimaYOmuraSMaruyamaM 2011 MicroRNA signature in diabetic wound healing: promotive role of miR-21 in fibroblast migration. Int Wound J

59. HeAZhuLGuptaNChangYFangF 2007 Overexpression of micro ribonucleic acid 29, highly up-regulated in diabetic rats, leads to insulin resistance in 3T3-L1 adipocytes. Mol Endocrinol 21 2785 2794

60. HerreraBMLockstoneHETaylorJMRiaMBarrettA 2010 Global microRNA expression profiles in insulin target tissues in a spontaneous rat model of type 2 diabetes. Diabetologia 53 1099 1109

61. PullenTJda Silva XavierGKelseyGRutterGA 2011 miR-29a and miR-29b contribute to pancreatic beta-cell-specific silencing of monocarboxylate transporter 1 (Mct1). Mol Cell Biol 31 3182 3194

62. ChenYZhuJLumPYYangXPintoS 2008 Variations in DNA elucidate molecular networks that cause disease. Nature 452 429 435

63. PerryMMMoschosSAWilliamsAEShepherdNJLarner-SvenssonHM 2008 Rapid changes in microRNA-146a expression negatively regulate the IL-1beta-induced inflammatory response in human lung alveolar epithelial cells. J Immunol 180 5689 5698

64. WilliamsAEPerryMMMoschosSALarner-SvenssonHMLindsayMA 2008 Role of miRNA-146a in the regulation of the innate immune response and cancer. Biochem Soc Trans 36 1211 1215

65. LindgrenCMHeidIMRandallJCLaminaCSteinthorsdottirV 2009 Genome-wide association scan meta-analysis identifies three Loci influencing adiposity and fat distribution. PLoS Genet 5 e1000508 doi:10.1371/journal.pgen.1000508

66. TeslovichTMMusunuruKSmithAVEdmondsonACStylianouIM 2010 Biological, clinical and population relevance of 95 loci for blood lipids. Nature 466 707 713

67. DuLSchagemanJJSubausteMCSaberBHammondSM 2009 miR-93, miR-98, and miR-197 regulate expression of tumor suppressor gene FUS1. Mol Cancer Res 7 1234 1243

68. SainiHKGriffiths-JonesSEnrightAJ 2007 Genomic analysis of human microRNA transcripts. Proc Natl Acad Sci U S A 104 17719 17724

69. IsikMKorswagenHCBerezikovE 2010 Expression patterns of intronic microRNAs in Caenorhabditis elegans. Silence 1 5

70. MartinezNJOwMCReece-HoyesJSBarrasaMIAmbrosVR 2008 Genome-scale spatiotemporal analysis of Caenorhabditis elegans microRNA promoter activity. Genome Res 18 2005 2015

71. OzsolakFPolingLLWangZLiuHLiuXS 2008 Chromatin structure analyses identify miRNA promoters. Genes Dev 22 3172 3183

72. BartelDP 2009 MicroRNAs: target recognition and regulatory functions. Cell 136 215 233

73. XiaoFZuoZCaiGKangSGaoX 2009 miRecords: an integrated resource for microRNA-target interactions. Nucleic Acids Res 37 D105 110

74. PapadopoulosGLReczkoMSimossisVASethupathyPHatzigeorgiouAG 2009 The database of experimentally supported targets: a functional update of TarBase. Nucleic Acids Res 37 D155 158

75. van DongenSAbreu-GoodgerCEnrightAJ 2008 Detecting microRNA binding and siRNA off-target effects from expression data. Nat Methods 5 1023 1025

76. RubyJGStarkAJohnstonWKKellisMBartelDP 2007 Evolution, biogenesis, expression, and target predictions of a substantially expanded set of Drosophila microRNAs. Genome Res 17 1850 1864

77. WarfMBJohnsonWEBassBL 2011 Improved annotation of C. elegans microRNAs by deep sequencing reveals structures associated with processing by Drosha and Dicer. RNA 17 563 577

78. Griffiths-JonesSHuiJHMarcoARonshaugenM 2011 MicroRNA evolution by arm switching. EMBO Rep 12 172 177

79. de WitELinsenSECuppenEBerezikovE 2009 Repertoire and evolution of miRNA genes in four divergent nematode species. Genome Res 19 2064 2074

80. OkamuraKPhillipsMDTylerDMDuanHChouYT 2008 The regulatory activity of microRNA* species has substantial influence on microRNA and 3′ UTR evolution. Nat Struct Mol Biol 15 354 363

81. LiangYRidzonDWongLChenC 2007 Characterization of microRNA expression profiles in normal human tissues. BMC Genomics 8 166

82. LandgrafPRusuMSheridanRSewerAIovinoN 2007 A mammalian microRNA expression atlas based on small RNA library sequencing. Cell 129 1401 1414

83. MamanovaLTurnerDJ 2011 Low-bias, strand-specific transcriptome Illumina sequencing by on-flowcell reverse transcription (FRT-seq). Nat Protoc 6 1736 1747

84. AlonSVigneaultFEminagaSChristodoulouDCSeidmanJG 2011 Barcoding bias in high-throughput multiplex sequencing of miRNA. Genome Res 21 1506 1511

85. Van NieuwerburghFSoetaertSPodshivalovaKAy-Lin WangESchafferL 2011 Quantitative Bias in Illumina TruSeq and a Novel Post Amplification Barcoding Strategy for Multiplexed DNA and Small RNA Deep Sequencing. PLoS ONE 6 e26969 doi:10.1371/journal.pone.0026969

86. ToedlingJServantNCiaudoCFarinelliLVoinnetO 2012 Deep-Sequencing Protocols Influence the Results Obtained in Small-RNA Sequencing. PLoS ONE 7 e32724 doi:10.1371/journal.pone.0032724

87. DahlmanILinderKArvidsson NordstromEAnderssonILidenJ 2005 Changes in adipose tissue gene expression with energy-restricted diets in obese women. Am J Clin Nutr 81 1275 1285

88. ParraPSerraFPalouA 2010 Expression of adipose microRNAs is sensitive to dietary conjugated linoleic acid treatment in mice. PLoS ONE 5 e13005 doi:10.1371/journal.pone.0013005

89. CushingLKuangPPQianJShaoFWuJ 2011 miR-29 is a major regulator of genes associated with pulmonary fibrosis. Am J Respir Cell Mol Biol 45 287 294

90. GotteMMohrCKooCYStockCVaskeAK 2010 miR-145-dependent targeting of junctional adhesion molecule A and modulation of fascin expression are associated with reduced breast cancer cell motility and invasiveness. Oncogene 29 6569 6580

91. Gaziel-SovranASeguraMFDi MiccoRCollinsMKHannifordD 2011 miR-30b/30d regulation of GalNAc transferases enhances invasion and immunosuppression during metastasis. Cancer Cell 20 104 118

92. WebsterRJGilesKMPriceKJZhangPMMattickJS 2009 Regulation of epidermal growth factor receptor signaling in human cancer cells by microRNA-7. J Biol Chem 284 5731 5741

93. ParkinsonHSarkansUKolesnikovNAbeygunawardenaNBurdettT 2011 ArrayExpress update–an archive of microarray and high-throughput sequencing-based functional genomics experiments. Nucleic Acids Res 39 D1002 1004

94. SkidmorePMCassidyASwaminathanRRichardsJBSpectorTD 2008 Relation of birth weight, body mass index, and change in size from birth to adulthood to insulin resistance in a female twin cohort. J Clin Endocrinol Metab 93 516 520

95. ProkopenkoILangenbergCFlorezJCSaxenaRSoranzoN 2009 Variants in MTNR1B influence fasting glucose levels. Nat Genet 41 77 81

96. FalchiMWilsonSGPaximadasDSwaminathanRSpectorTD 2008 Quantitative linkage analysis for pancreatic B-cell function and insulin resistance in a large twin cohort. Diabetes 57 1120 1124

97. R Core Development Team 2010 R: A language and environment for statistical computing: R Foundation for Statistical Computing, Vienna, Austria

98. GentlemanRCCareyVJBatesDMBolstadBDettlingM 2004 Bioconductor: open software development for computational biology and bioinformatics. Genome Biol 5 R80

99. MorganMAndersSLawrenceMAboyounPPagesH 2009 ShortRead: a bioconductor package for input, quality assessment and exploration of high-throughput sequence data. Bioinformatics 25 2607 2608

100. LiHDurbinR 2009 Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25 1754 1760

101. LangmeadBTrapnellCPopMSalzbergSL 2009 Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol 10 R25

102. AndersSHuberW 2010 Differential expression analysis for sequence count data. Genome Biol 11 R106

103. LiHHandsakerBWysokerAFennellTRuanJ 2009 The Sequence Alignment/Map format and SAMtools. Bioinformatics 25 2078 2079

104. IvoLHofackerWFStadlerPeterFBonhoefferSebastianTackerManfredSchusterPeter 1994 Fast folding and comparison of RNA secondary structures. Mh Chemie 125 167 188

105. LeeMTKimJ 2008 Self containment, a property of modular RNA structures, distinguishes microRNAs. PLoS Comput Biol 4 e1000150 doi:10.1371/journal.pcbi.1000150

106. BarrettTTroupDBWilhiteSELedouxPEvangelistaC 2011 NCBI GEO: archive for functional genomics data sets–10 years on. Nucleic Acids Res 39 D1005 1010

107. SeanDMeltzerPS 2007 GEOquery: a bridge between the Gene Expression Omnibus (GEO) and BioConductor. Bioinformatics 23 1846 1847

108. SmythGK 2005 Limma: linear models for microarray data. Gentleman VCRDudoitSIrizarryRHuberW Bioinformatics and Computational Biology Solutions using R and Bioconductor New York Springer 397 420

109. FriedmanRCFarhKKBurgeCBBartelDP 2009 Most mammalian mRNAs are conserved targets of microRNAs. Genome Res 19 92 105