The objective of this study was to identify the expression profile of long noncoding RNAs (lncRNAs) with concurrent assessment of microRNA (miRNA) and messenger RNA (mRNA) profiles in leiomyomas and paired myometrium using next-generation RNA sequencing and assembly of RNA transcripts. Total RNA was isolated from leiomyoma and paired myometrium (N = 8) and samples from 3 pairs were subjected to RNA sequencing. Normalized assembly of over 48 000 lncRNAs resulted in identification of 45 936 lncRNAs. Of these lncRNAs, 22 148 representing overlapping, intergenic, intronic, and antisense subtypes were expressed in all paired tissues, with 5941 (2813 up- and 3128 downregulated at ≥1.5 fold) differentially expressed in leiomyomas. Concurrent RNA sequencing revealed the expression of 2588 miRNAs and 21 814 mRNAs, of which 392 miRNAs and 16 559 mRNAs were expressed in all paired tissues. Of these transcripts, 56 and 92 miRNAs and 2030 and 1825 mRNAs were up- or downregulated at ≥1.5 fold, respectively, in leiomyoma as compared to myometrium. Using quantitative reverse transcription-polymerase chain reaction (QRT-PCR), we confirmed the expression of hepatocellular carcinoma upregulated (HULC), lnc-maternally expressed 3 (MEG3), long intergenic ncRNA 890 (LINC00890), TSIX, long intergenic ncRNA 473 (LINC00473), lnc-KLF9-1, and lnc-POTEM-3 (lncRNA-ATB) in leiomyoma and matched myometrium (N = 8). Collectively, the results presented here provide a comprehensive expression profile of lncRNAs in leiomyomas with concurrent integrated expression of miRNAs and mRNAs and implicate potential regulatory functions of lncRNAs through interactions with specific miRNAs and mRNAs which are known to be critical in the pathogenesis of leiomyoma.
SegarsJHParrottECNagelJD. Proceedings from the third National Institutes of Health International Congress on Advances in Uterine Leiomyoma Research: comprehensive review, conference summary and future recommendations. Hum Reprod Update. 2014;20(3):309–333.
2.
CheginiNVeralaJLuoXXuJWilliamsRS. Gene expression profile of leiomyoma and myometrium and the effect of gonadotropin releasing hormone analogue therapy. J Soc Gynecol Investig. 2003;10(3):161–171.
3.
GuoHZhangXDongR. Integrated analysis of long noncoding RNAs and mRNAs reveals their potential roles in the pathogenesis of uterine leiomyomas. Oncotarget. 2014;5(8):8625–8636.
4.
LuoXPanQLiuLCheginiN. Genomic and proteomic profiling II: comparative assessment of gene expression profiles in leiomyomas, keloids, and surgically-induced scars. Reprod Biol Endocrinol. 2007;5:35.
5.
PanQLuoXCheginiN. Genomic and proteomic profiling I: leiomyomas in African Americans and Caucasians. Reprod Biol Endocrinol. 2007;5:34.
6.
TsibrisJCSegarsJCoppolaD. Insights from gene arrays on the development and growth regulation of uterine leiomyomata. Fertil Steril. 2002;78(1):114–121.
7.
CheginiN. Proinflammatory and profibrotic mediators: principal effectors of leiomyoma development as a fibrotic disorder. Semin Reprod Med. 2010;28(3):180–203.
8.
YatsenkoSAMittalPWood-TrageserMA. Highly heterogeneous genomic landscape of uterine leiomyomas by whole exome sequencing and genome-wide arrays. Fertil Steril. 2017;107(2):457–466.e9.
9.
MehineMKaasinenEMakinenN. Characterization of uterine leiomyomas by whole-genome sequencing. N Engl J Med. 2013;369(1):43–53.
10.
MehineMMakinenNHeinonenHRAaltonenLAVahteristoP. Genomics of uterine leiomyomas: insights from high-throughput sequencing. Fertil Steril. 2014;102(3):621–629.
11.
HodgeJCKimTMDreyfussJM. Expression profiling of uterine leiomyomata cytogenetic subgroups reveals distinct signatures in matched myometrium: transcriptional profiling of the t(12;14) and evidence in support of predisposing genetic heterogeneity. Hum Mol Genet. 2012;21(10):2312–2329.
12.
GallagherCSMortonCC. Genetic association studies in uterine fibroids: risk alleles presage the path to personalized therapies. Semin Reprod Med. 2016;34(4):235–241.
13.
Liegl-AtzwangerBHeitzerEFlickerK. Exploring chromosomal abnormalities and genetic changes in uterine smooth muscle tumors. Mod Pathol. 2016;29(10):1262–1277.
14.
MehineMKaasinenEHeinonenHR. Integrated data analysis reveals uterine leiomyoma subtypes with distinct driver pathways and biomarkers. Proc Natl Acad Sci U S A. 2016;113(5):1315–1320.
Al-HendyALaknaurADiamondMPIsmailNBoyerTGHalderSK. Silencing Med12 gene reduces proliferation of human leiomyoma cells mediated via Wnt/beta-catenin signaling pathway. Endocrinology. 2017;158(3):592–603.
17.
GeorgievaBMilevIMinkovIDimitrovaIBradfordAPBaevV. Characterization of the uterine leiomyoma microRNAome by deep sequencing. Genomics. 2012;99(5):275–281.
18.
LuoXCheginiN. The expression and potential regulatory function of microRNAs in the pathogenesis of leiomyoma. Semin Reprod Med. 2008;26(6):500–514.
19.
WangTZhangXObijuruL. A micro-RNA signature associated with race, tumor size, and target gene activity in human uterine leiomyomas. Genes Chromosomes Cancer. 2007;46(4):336–347.
20.
ZavadilJYeHLiuZ. Profiling and functional analyses of microRNAs and their target gene products in human uterine leiomyomas. PloS One. 2010;5(8):e12362.
21.
ChuangTDKhorramO. miR-200c regulates IL8 expression by targeting IKBKB: a potential mediator of inflammation in leiomyoma pathogenesis. PloS One. 2014;9(4):e95370.
22.
ChuangTDKhorramO. Mechanisms underlying aberrant expression of miR-29c in uterine leiomyoma. Fertil Steril. 2016;105(1):236–245.e1.
HuntzingerEIzaurraldeE. Gene silencing by microRNAs: contributions of translational repression and mRNA decay. Nat Rev Genet. 2011;12(2):99–110.
25.
SelbachMSchwanhausserBThierfelderNFangZKhaninRRajewskyN. Widespread changes in protein synthesis induced by microRNAs. Nature. 2008;455(7209):58–63.
26.
MaXBecker BuscagliaLEBarkerJRLiY. MicroRNAs in NF-kappaB signaling. Journal of Molecular Cell Biol. 2011;3(3):159–166.
27.
SchickelRBoyerinasBParkSMPeterME. MicroRNAs: key players in the immune system, differentiation, tumorigenesis and cell death. Oncogene. 2008;27(45):5959–5974.
28.
ChuangTDPandaHLuoXCheginiN. miR-200c is aberrantly expressed in leiomyomas in an ethnic-dependent manner and targets ZEBs, VEGFA, TIMP2, and FBLN5. Endocr Relat Cancer. 2012;19(4):541–556.
29.
ChuangTDLuoXPandaHCheginiN. miR-93/106b and their host gene, MCM7, are differentially expressed in leiomyomas and functionally target F3 and IL-8. Mol Endocrinol. 2012;26(6):1028–1042.
30.
FitzgeraldJBChennathukuzhiVKoohestaniFNowakRAChristensonLK. Role of microRNA-21 and programmed cell death 4 in the pathogenesis of human uterine leiomyomas. Fertil Steril. 2012;98(3):726–734.e2.
31.
KlemkeMMeyerAHashemi NezhadMBelgeGBartnitzkeSBullerdiekJ. Loss of let-7 binding sites resulting from truncations of the 3’ untranslated region of HMGA2 mRNA in uterine leiomyomas. Cancer Genet Cytogenet. 2010;196(2):119–123.
32.
QiangWLiuZSernaVA. Down-regulation of miR-29b is essential for pathogenesis of uterine leiomyoma. Endocrinology. 2014;155(3):663–669.
33.
BeermannJPiccoliMTViereckJThumT. Non-coding RNAs in development and disease: background, mechanisms, and therapeutic approaches. Physiol Rev. 2016;96(4):1297–1325.
34.
JiangCLiXZhaoHLiuH. Long non-coding RNAs: potential new biomarkers for predicting tumor invasion and metastasis. Mol Cancer. 2016;15(1):62.
35.
LiJTianHYangJGongZ. Long noncoding RNAs regulate cell Growth, proliferation, and apoptosis. DNA Cell Biol. 2016;35(9):459–470.
36.
UlitskyIBartelDP. lincRNAs: genomics, evolution, and mechanisms. Cell. 2013;154(1):26–46.
37.
HangauerMJVaughnIWMcManusMT. Pervasive transcription of the human genome produces thousands of previously unidentified long intergenic noncoding RNAs. PLoS Genet. 2013;9(6):e1003569.
38.
SigovaAAMullenACMolinieB. Divergent transcription of long noncoding RNA/mRNA gene pairs in embryonic stem cells. Proc Natl Acad Sci U S A. 2013;110(8):2876–2881.
39.
CaoMXJiangYPTangYLLiangXH. The crosstalk between lncRNA and microRNA in cancer metastasis: orchestrating the epithelial-mesenchymal plasticity. Oncotarget. 2017;8(7):12472–12483.
40.
LangmeadBTrapnellCPopMSalzbergSL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):R25.
41.
LangmeadBSalzbergSL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–359.
42.
LiBDeweyCN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323.
43.
SunZEvansJBhagwateA. CAP-miRSeq: a comprehensive analysis pipeline for microRNA sequencing data. BMC Genomics. 2014;15:423.
44.
VoldersPJHelsensKWangX. LNCipedia: a database for annotated human lncRNA transcript sequences and structures. Nucleic Acids Res. 2013;41(Database issue):D246–D251.
45.
VoldersPJVerheggenKMenschaertG. An update on LNCipedia: a database for annotated human lncRNA sequences. Nucleic Acids Res. 2015;43(8):4363–4364.
46.
McCarthyDJChenYSmythGK. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 2012;40(10):4288–4297.
47.
RobinsonMDMcCarthyDJSmythGK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–140.
48.
Huang daWShermanBTLempickiRA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
49.
Huang daWShermanBTLempickiRA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.
50.
ChuangTDKhorramO. Tranilast inhibits genes functionally involved in cell proliferation, fibrosis, and epigenetic regulation and epigenetically induces miR-29c expression in leiomyoma cells [Published online January 1, 2016]. Reprod Sci. 1933719116682878.
51.
ChuangTDHoMKhorramO. The regulatory function of miR-200c on inflammatory and cell-cycle associated genes in SK-LMS-1, A leiomyosarcoma cell line. Reprod Sci. 2015;22(5):563–571.
52.
ChuangTDPearceWJKhorramO. miR-29c induction contributes to downregulation of vascular extracellular matrix proteins by glucocorticoids. Am J Physiol Cell Physiol. 2015;309(2):C117–C125.
53.
KhorramOChuangTDPearceWJ. Long-term effects of maternal undernutrition on offspring carotid artery remodeling: role of miR-29c. J Dev Orig Health Dis. 2015;6(4):342–349.
54.
AlmeidaTAQuispe-RicaldeAMontes de OcaFForondaPHernandezMM. A high-throughput open-array qPCR gene panel to identify housekeeping genes suitable for myometrium and leiomyoma expression analysis. Gynecol Oncol. 2014;134(1):138–143.
55.
BraconiCKogureTValeriN. microRNA-29 can regulate expression of the long non-coding RNA gene MEG3 in hepatocellular cancer. Oncogene. 2011;30(47):4750–4756.
56.
DoJTHanDWGentileLSobek-KlockeIStehlingMScholerHR. Enhanced reprogramming of Xist by induced upregulation of Tsix and Dnmt3a. Stem Cells (Dayton, Ohio). 2008;26(11):2821–2831.
57.
WangZJinninMNakamuraK. Long non-coding RNA TSIX is upregulated in scleroderma dermal fibroblasts and controls collagen mRNA stabilization. Exp Dermatol. 2016;25(2):131–136.
58.
YuanJHYangFWangF. A long noncoding RNA activated by TGF-beta promotes the invasion-metastasis cascade in hepatocellular carcinoma. Cancer Cell. 2014;25(5):666–681.
59.
ShiSJWangLJYuBLiYHJinYBaiXZ. LncRNA-ATB promotes trastuzumab resistance and invasion-metastasis cascade in breast cancer. Oncotarget. 2015;6(13):11652–11663.
60.
TerashimaMTangeSIshimuraASuzukiT. MEG3 long noncoding RNA contributes to the epigenetic regulation of epithelial-mesenchymal transition in lung cancer cell lines. J Biol Chem. 2017;292(1):82–99.
61.
YuXZhengHChanMTWuWK. HULC: an oncogenic long non-coding RNA in human cancer. J Cell Mol Med. 2017;21(2):410–417.
62.
LiSPXuHXYuY. LncRNA HULC enhances epithelial-mesenchymal transition to promote tumorigenesis and metastasis of hepatocellular carcinoma via the miR-200a-3p/ZEB1 signaling pathway. Oncotarget. 2016;7(27):42431–42446.
63.
LiuXChenHKongW. Down-regulated long non-coding RNA-ATB in preeclampsia and its effect on suppressing migration, proliferation, and tube formation of trophoblast cells. Placenta. 2017;49:80–87.
64.
UlitskyI. Evolution to the rescue: using comparative genomics to understand long non-coding RNAs. Nat Rev Genet. 2016;17(10):601–614.
65.
OromUADerrienTGuigoRShiekhattarR. Long noncoding RNAs as enhancers of gene expression. Cold Spring Harb Symp Quant Biol. 2010;75:325–331.
66.
MooreA. Commissioning support. The future will be integrated. Health Serv J. 2013;123(6344):suppl 10–suppl 11.
67.
HuangXYuanTTschannenM. Characterization of human plasma-derived exosomal RNAs by deep sequencing. BMC Genomics. 2013;14:319.
68.
ViereckJThumT. Circulating noncoding RNAs as biomarkers of cardiovascular disease and injury. Circ Res. 2017;120(2):381–399.
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
