The evolution of lncRNA repertoires and expression patterns in tetrapods.

Only a very small fraction of long noncoding RNAs (lncRNAs) are well characterized. The evolutionary history of lncRNAs can provide insights into their functionality, but the absence of lncRNA annotations in non-model organisms has precluded comparative analyses. Here we present a large-scale evolutionary study of lncRNA repertoires and expression patterns, in 11 tetrapod species. We identify approximately 11,000 primate-specific lncRNAs and 2,500 highly conserved lncRNAs, including approximately 400 genes that are likely to have originated more than 300 million years ago. We find that lncRNAs, in particular ancient ones, are in general actively regulated and may function predominantly in embryonic development. Most lncRNAs evolve rapidly in terms of sequence and expression levels, but tissue specificities are often conserved. We compared expression patterns of homologous lncRNA and protein-coding families across tetrapods to reconstruct an evolutionarily conserved co-expression network. This network suggests potential functions for lncRNAs in fundamental processes such as spermatogenesis and synaptic transmission, but also in more specific mechanisms such as placenta development through microRNA production.

H19 lncRNA alters stromal cell growth via IGF signaling in the endometrium of women with endometriosis

Endometriosis affects approximately 15% of reproductive aged women and is associated with chronic pelvic pain and infertility. However, the molecular mechanisms by which endometriosis impacts fertility are poorly understood. The developmentally regulated, imprinted H19 long noncoding RNA (lncRNA) functions to reduce the bioavailability of microRNA let-7 by acting as a molecular sponge. Here we report that H19 expression is significantly decreased in the eutopic endometrium of women with endometriosis as compared to normal controls. We show that decreased H19 increases let-7 activity, which in turn inhibits Igf1r expression at the post-transcriptional level, thereby contributing to reduced proliferation of endometrial stromal cells. We propose that perturbation of this newly identified H19/Let-7/IGF1R regulatory pathway may contribute to impaired endometrial preparation and receptivity for pregnancy in women with endometriosis. Our finding represents the first example of a lncRNA-based mechanism in endometriosis and its associated infertility, thus holding potential in the development of novel therapeutics for women with endometriosis and infertility.

H19 lncRNA alters DNA methylation genome wide by regulating S-adenosylhomocysteine hydrolase

DNA methylation is essential for mammalian development and physiology. Here we report that the developmentally regulated H19 lncRNA binds to and inhibits S-adenosylhomocysteine hydrolase (SAHH), the only mammalian enzyme capable of hydrolysing S-adenosylhomocysteine (SAH). SAH is a potent feedback inhibitor of S-adenosylmethionine (SAM)-dependent methyltransferases that methylate diverse cellular components, including DNA, RNA, proteins, lipids and neurotransmitters. We show that H19 knockdown activates SAHH, leading to increased DNMT3B-mediated methylation of an lncRNA-encoding gene Nctc1 within the Igf2-H19-Nctc1 locus. Genome-wide methylation profiling reveals methylation changes at numerous gene loci consistent with SAHH modulation by H19. Our results uncover an unanticipated regulatory circuit involving broad epigenetic alterations by a single abundantly expressed lncRNA that may underlie gene methylation dynamics of development and diseases and suggest that this mode of regulation may extend to other cellular components.

TAZ role in lung cancer: resistance to BET inhibitors and functional interaction with lncRNA

Despite extensive research and introduction of innovative therapy, lung cancer prognosis remains poor, with a five years survival of only 17%. The success of pharmacological treatment is often impaired by drug resistance. Thus, the characterization of response mechanisms to anti-cancer compounds and of the molecular mechanisms supporting lung cancer aggressiveness are crucial for patient’s management. In the first part of this thesis, we characterized the molecular mechanism behind resistance of lung cancer cells to the Inhibitors of the Bromodomain and Extraterminal domain containing Proteins (BETi). Through a CRISPR/Cas9 screening we identified three Hippo Pathway members, LATS2, TAOK1 and NF2 as genes implicated in susceptibility to BETi. These genes confer sensitivity to BETi inhibiting TAZ activity. Conversely, TAZ overexpression increases resistance to BETi. We also displayed that BETi downregulate both YAP, TAZ and TEADs expression in several cancer cell lines, implying a novel BETi-dependent cytotoxic mechanism. In the second part of this work, we attempted to characterize the crosstalk between the TAZ gene and its cognate antisense long-non coding RNA (lncRNA) TAZ-AS202 in lung tumorigenesis. As for TAZ downregulation, TAZ-AS202 silencing impairs NSCLC cells proliferation, migration and invasion, suggesting a pro-tumorigenic function for this lncRNA during lung tumorigenesis. TAZ-AS202 regulates TAZ target genes without altering TAZ expression or localization. This finding implies an uncovered functional cooperation between TAZ and TAZ-AS202. Moreover, we found that the EPH-ephrin signaling receptor EPHB2 is a downstream effector affected by both TAZ and TAZ-AS202 silencing. EPHB2 downregulation significantly attenuates cells proliferation, migration and invasion, suggesting that, at least in part, TAZ-AS202 and TAZ pro-oncogenic activity depends on EPH-ephrin signaling final deregulation. Finally, we started to dissect the mechanism underlying the TAZ-AS202 regulatory activity on EPHB2 in lung cancer, which may involve the existence of an intermediate transcription factor and is the object of our ongoing research.

The lincRNA HOTAIRM1, located in the HOXA genomic region, is expressed in acute myeloid leukemia, impacts prognosis in patients in the intermediate-risk cytogenetic category, and is associated with a distinctive microRNA signature.

Long non-coding RNAs (lncRNAs) are deregulated in several tumors, although their role in acute myeloid leukemia (AML) is mostly unknown.We have examined the expression of the lncRNA HOX antisense intergenic RNA myeloid 1 (HOTAIRM1) in 241 AML patients. We have correlated HOTAIRM1 expression with a miRNA expression profile. We have also analyzed the prognostic value of HOTAIRM1 expression in 215 intermediate-risk AML (IR-AML) patients.The lowest expression level was observed in acute promyelocytic leukemia (P < 0.001) and the highest in t(6;9) AML (P = 0.005). In 215 IR-AML patients, high HOTAIRM1 expression was independently associated with shorter overall survival (OR:2.04;P = 0.001), shorter leukemia-free survival (OR:2.56; P < 0.001) and a higher cumulative incidence of relapse (OR:1.67; P = 0.046). Moreover, HOTAIRM1 maintained its independent prognostic value within the favorable molecular subgroup (OR: 3.43; P = 0.009). Interestingly, HOTAIRM1 was overexpressed in NPM1-mutated AML (P < 0.001) and within this group retained its prognostic value (OR: 2.21; P = 0.01). Moreover, HOTAIRM1 expression was associated with a specific 33-microRNA signature that included miR-196b (P < 0.001). miR-196b is located in the HOX genomic region and has previously been reported to have an independent prognostic value in AML. miR-196b and HOTAIRM1 in combination as a prognostic factor can classify patients as high-, intermediate-, or low-risk (5-year OS: 24% vs 42% vs 70%; P = 0.004).Determination of HOTAIRM1 level at diagnosis provided relevant prognostic information in IR-AML and allowed refinement of risk stratification based on common molecular markers. The prognostic information provided by HOTAIRM1 was strengthened when combined with miR-196b expression. Furthermore, HOTAIRM1 correlated with a 33-miRNA signature.

GENCODE: the reference human genome annotation for The ENCODE Project.

The GENCODE Consortium aims to identify all gene features in the human genome using a combination of computational analysis, manual annotation, and experimental validation. Since the first public release of this annotation data set, few new protein-coding loci have been added, yet the number of alternative splicing transcripts annotated has steadily increased. The GENCODE 7 release contains 20,687 protein-coding and 9640 long noncoding RNA loci and has 33,977 coding transcripts not represented in UCSC genes and RefSeq. It also has the most comprehensive annotation of long noncoding RNA (lncRNA) loci publicly available with the predominant transcript form consisting of two exons. We have examined the completeness of the transcript annotation and found that 35% of transcriptional start sites are supported by CAGE clusters and 62% of protein-coding genes have annotated polyA sites. Over one-third of GENCODE protein-coding genes are supported by peptide hits derived from mass spectrometry spectra submitted to Peptide Atlas. New models derived from the Illumina Body Map 2.0 RNA-seq data identify 3689 new loci not currently in GENCODE, of which 3127 consist of two exon models indicating that they are possibly unannotated long noncoding loci. GENCODE 7 is publicly available from gencodegenes.org and via the Ensembl and UCSC Genome Browsers.

Genome-wide profiling of the cardiac transcriptome after myocardial infarction identifies novel heart-specific long non-coding RNAs.

AIM: Heart disease is recognized as a consequence of dysregulation of cardiac gene regulatory networks. Previously, unappreciated components of such networks are the long non-coding RNAs (lncRNAs). Their roles in the heart remain to be elucidated. Thus, this study aimed to systematically characterize the cardiac long non-coding transcriptome post-myocardial infarction and to elucidate their potential roles in cardiac homoeostasis. METHODS AND RESULTS: We annotated the mouse transcriptome after myocardial infarction via RNA sequencing and ab initio transcript reconstruction, and integrated genome-wide approaches to associate specific lncRNAs with developmental processes and physiological parameters. Expression of specific lncRNAs strongly correlated with defined parameters of cardiac dimensions and function. Using chromatin maps to infer lncRNA function, we identified many with potential roles in cardiogenesis and pathological remodelling. The vast majority was associated with active cardiac-specific enhancers. Importantly, oligonucleotide-mediated knockdown implicated novel lncRNAs in controlling expression of key regulatory proteins involved in cardiogenesis. Finally, we identified hundreds of human orthologues and demonstrate that particular candidates were differentially modulated in human heart disease. CONCLUSION: These findings reveal hundreds of novel heart-specific lncRNAs with unique regulatory and functional characteristics relevant to maladaptive remodelling, cardiac function and possibly cardiac regeneration. This new class of molecules represents potential therapeutic targets for cardiac disease. Furthermore, their exquisite correlation with cardiac physiology renders them attractive candidate biomarkers to be used in the clinic.

The striatal long noncoding RNA Abhd11os is neuroprotective against an N-terminal fragment of mutant huntingtin in vivo.

A large number of gene products that are enriched in the striatum have ill-defined functions, although they may have key roles in age-dependent neurodegenerative diseases affecting the striatum, especially Huntington disease (HD). In the present study, we focused on Abhd11os, (called ABHD11-AS1 in human) which is a putative long noncoding RNA (lncRNA) whose expression is enriched in the mouse striatum. We confirm that despite the presence of 2 small open reading frames (ORFs) in its sequence, Abhd11os is not translated into a detectable peptide in living cells. We demonstrate that Abhd11os levels are markedly reduced in different mouse models of HD. We performed in vivo experiments in mice using lentiviral vectors encoding either Abhd11os or a small hairpin RNA targeting Abhd11os. Results show that Abhd11os overexpression produces neuroprotection against an N-terminal fragment of mutant huntingtin, whereas Abhd11os knockdown is protoxic. These novel results indicate that the loss lncRNA Abhd11os likely contribute to striatal vulnerability in HD. Our study emphasizes that lncRNA may play crucial roles in neurodegenerative diseases.

Extensive microRNA-mediated crosstalk between lncRNAs and mRNAs in mouse embryonic stem cells.

Recently, a handful of intergenic long noncoding RNAs (lncRNAs) have been shown to compete with mRNAs for binding to miRNAs and to contribute to development and disease. Beyond these reports, little is yet known of the extent and functional consequences of miRNA-mediated regulation of mRNA levels by lncRNAs. To gain further insight into lncRNA-mRNA miRNA-mediated crosstalk, we reanalyzed transcriptome-wide changes induced by the targeted knockdown of over 100 lncRNA transcripts in mouse embryonic stem cells (mESCs). We predicted that, on average, almost one-fifth of the transcript level changes induced by lncRNAs are dependent on miRNAs that are highly abundant in mESCs. We validated these findings experimentally by temporally profiling transcriptome-wide changes in gene expression following the loss of miRNA biogenesis in mESCs. Following the depletion of miRNAs, we found that >50% of lncRNAs and their miRNA-dependent mRNA targets were up-regulated coordinately, consistent with their interaction being miRNA-mediated. These lncRNAs are preferentially located in the cytoplasm, and the response elements for miRNAs they share with their targets have been preserved in mammals by purifying selection. Lastly, miRNA-dependent mRNA targets of each lncRNA tended to share common biological functions. Post-transcriptional miRNA-mediated crosstalk between lncRNAs and mRNA, in mESCs, is thus surprisingly prevalent, conserved in mammals, and likely to contribute to critical developmental processes.

Assessing Recent Selection and Functionality at Long Noncoding RNA Loci in the Mouse Genome.

Long noncoding RNAs (lncRNAs) are one of the most intensively studied groups of noncoding elements. Debate continues over what proportion of lncRNAs are functional or merely represent transcriptional noise. Although characterization of individual lncRNAs has identified approximately 200 functional loci across the Eukarya, general surveys have found only modest or no evidence of long-term evolutionary conservation. Although this lack of conservation suggests that most lncRNAs are nonfunctional, the possibility remains that some represent recent evolutionary innovations. We examine recent selection pressures acting on lncRNAs in mouse populations. We compare patterns of within-species nucleotide variation at approximately 10,000 lncRNA loci in a cohort of the wild house mouse, Mus musculus castaneus, with between-species nucleotide divergence from the rat (Rattus norvegicus). Loci under selective constraint are expected to show reduced nucleotide diversity and divergence. We find limited evidence of sequence conservation compared with putatively neutrally evolving ancestral repeats (ARs). Comparisons of sequence diversity and divergence between ARs, protein-coding (PC) exons and lncRNAs, and the associated flanking regions, show weak, but significantly lower levels of sequence diversity and divergence at lncRNAs compared with ARs. lncRNAs conserved deep in the vertebrate phylogeny show lower within-species sequence diversity than lncRNAs in general. A set of 74 functionally characterized lncRNAs show levels of diversity and divergence comparable to PC exons, suggesting that these lncRNAs are under substantial selective constraints. Our results suggest that, in mouse populations, most lncRNA loci evolve at rates similar to ARs, whereas older lncRNAs tend to show signals of selection similar to PC genes.

Involvement of long non-coding RNAs in beta cell failure at the onset of type 1 diabetes in NOD mice.

AIMS/HYPOTHESIS: Exposure of pancreatic beta cells to cytokines released by islet-infiltrating immune cells induces alterations in gene expression, leading to impaired insulin secretion and apoptosis in the initial phases of type 1 diabetes. Long non-coding RNAs (lncRNAs) are a new class of transcripts participating in the development of many diseases. As little is known about their role in insulin-secreting cells, this study aimed to evaluate their contribution to beta cell dysfunction. METHODS: The expression of lncRNAs was determined by microarray in the MIN6 beta cell line exposed to proinflammatory cytokines. The changes induced by cytokines were further assessed by real-time PCR in islets of control and NOD mice. The involvement of selected lncRNAs modified by cytokines was assessed after their overexpression in MIN6 cells and primary islet cells. RESULTS: MIN6 cells were found to express a large number of lncRNAs, many of which were modified by cytokine treatment. The changes in the level of selected lncRNAs were confirmed in mouse islets and an increase in these lncRNAs was also seen in prediabetic NOD mice. Overexpression of these lncRNAs in MIN6 and mouse islet cells, either alone or in combination with cytokines, favoured beta cell apoptosis without affecting insulin production or secretion. Furthermore, overexpression of lncRNA-1 promoted nuclear translocation of nuclear factor of κ light polypeptide gene enhancer in B cells 1 (NF-κB). CONCLUSIONS/INTERPRETATION: Our study shows that lncRNAs are modulated during the development of type 1 diabetes in NOD mice, and that their overexpression sensitises beta cells to apoptosis, probably contributing to their failure during the initial phases of the disease.

Comprehensive Genomic Characterization of Long Non-coding RNAs across Human Cancers.

The discovery of long non-coding RNA (lncRNA) has dramatically altered our understanding of cancer. Here, we describe a comprehensive analysis of lncRNA alterations at transcriptional, genomic, and epigenetic levels in 5,037 human tumor specimens across 13 cancer types from The Cancer Genome Atlas. Our results suggest that the expression and dysregulation of lncRNAs are highly cancer type specific compared with protein-coding genes. Using the integrative data generated by this analysis, we present a clinically guided small interfering RNA screening strategy and a co-expression analysis approach to identify cancer driver lncRNAs and predict their functions. This provides a resource for investigating lncRNAs in cancer and lays the groundwork for the development of new diagnostics and treatments.

Long noncoding intronic RNAs are differentially expressed in primary and metastatic pancreatic cancer

Abstract Background Pancreatic ductal adenocarcinoma (PDAC) is known by its aggressiveness and lack of effective therapeutic options. Thus, improvement in current knowledge of molecular changes associated with pancreatic cancer is urgently needed to explore novel venues of diagnostics and treatment of this dismal disease. While there is mounting evidence that long noncoding RNAs (lncRNAs) transcribed from intronic and intergenic regions of the human genome may play different roles in the regulation of gene expression in normal and cancer cells, their expression pattern and biological relevance in pancreatic cancer is currently unknown. In the present work we investigated the relative abundance of a collection of lncRNAs in patients' pancreatic tissue samples aiming at identifying gene expression profiles correlated to pancreatic cancer and metastasis. Methods Custom 3,355-element spotted cDNA microarray interrogating protein-coding genes and putative lncRNA were used to obtain expression profiles from 38 clinical samples of tumor and non-tumor pancreatic tissues. Bioinformatics analyses were performed to characterize structure and conservation of lncRNAs expressed in pancreatic tissues, as well as to identify expression signatures correlated to tissue histology. Strand-specific reverse transcription followed by PCR and qRT-PCR were employed to determine strandedness of lncRNAs and to validate microarray results, respectively. Results We show that subsets of intronic/intergenic lncRNAs are expressed across tumor and non-tumor pancreatic tissue samples. Enrichment of promoter-associated chromatin marks and over-representation of conserved DNA elements and stable secondary structure predictions suggest that these transcripts are generated from independent transcriptional units and that at least a fraction is under evolutionary selection, and thus potentially functional. Statistically significant expression signatures comprising protein-coding mRNAs and lncRNAs that correlate to PDAC or to pancreatic cancer metastasis were identified. Interestingly, loci harboring intronic lncRNAs differentially expressed in PDAC metastases were enriched in genes associated to the MAPK pathway. Orientation-specific RT-PCR documented that intronic transcripts are expressed in sense, antisense or both orientations relative to protein-coding mRNAs. Differential expression of a subset of intronic lncRNAs (PPP3CB, MAP3K14 and DAPK1 loci) in metastatic samples was confirmed by Real-Time PCR. Conclusion Our findings reveal sets of intronic lncRNAs expressed in pancreatic tissues whose abundance is correlated to PDAC or metastasis, thus pointing to the potential relevance of this class of transcripts in biological processes related to malignant transformation and metastasis in pancreatic cancer.

Expression analysis and in silico characterization of intronic long noncoding RNAs in renal cell carcinoma: emerging functional associations

Abstract Background Intronic and intergenic long noncoding RNAs (lncRNAs) are emerging gene expression regulators. The molecular pathogenesis of renal cell carcinoma (RCC) is still poorly understood, and in particular, limited studies are available for intronic lncRNAs expressed in RCC Methods Microarray experiments were performed with custom-designed arrays enriched with probes for lncRNAs mapping to intronic genomic regions. Samples from 18 primary RCC tumors and 11 nontumor adjacent matched tissues were analyzed. Meta-analyses were performed with microarray expression data from three additional human tissues (normal liver, prostate tumor and kidney nontumor samples), and with large-scale public data for epigenetic regulatory marks and for evolutionarily conserved sequences. Results A signature of 29 intronic lncRNAs differentially expressed between RCC and nontumor samples was obtained (false discovery rate (FDR) <5%). A signature of 26 intronic lncRNAs significantly correlated with the RCC five-year patient survival outcome was identified (FDR <5%, p-value ≤0.01). We identified 4303 intronic antisense lncRNAs expressed in RCC, of which 22% were significantly (p <0.05) cis correlated with the expression of the mRNA in the same locus across RCC and three other human tissues. Gene Ontology (GO) analysis of those loci pointed to 'regulation of biological processes’ as the main enriched category. A module map analysis of the protein-coding genes significantly (p <0.05) trans correlated with the 20% most abundant lncRNAs, identified 51 enriched GO terms (p <0.05). We determined that 60% of the expressed lncRNAs are evolutionarily conserved. At the genomic loci containing the intronic RCC-expressed lncRNAs, a strong association (p <0.001) was found between their transcription start sites and genomic marks such as CpG islands, RNA Pol II binding and histones methylation and acetylation. Conclusion Intronic antisense lncRNAs are widely expressed in RCC tumors. Some of them are significantly altered in RCC in comparison with nontumor samples. The majority of these lncRNAs is evolutionarily conserved and possibly modulated by epigenetic modifications. Our data suggest that these RCC lncRNAs may contribute to the complex network of regulatory RNAs playing a role in renal cell malignant transformation.

The intronic long noncoding RNA ANRASSF1 recruits PRC2 to the RASSF1A promoter, reducing the expression of RASSF1A and increasing cell proliferation

The down-regulation of the tumor-suppressor gene RASSF1A has been shown to increase cell proliferation in several tumors. RASSF1A expression is regulated through epigenetic events involving the polycomb repressive complex 2 (PRC2); however, the molecular mechanisms modulating the recruitment of this epigenetic modifier to the RASSF1 locus remain largely unknown. Here, we identify and characterize ANRASSF1, an endogenous unspliced long noncoding RNA (lncRNA) that is transcribed from the opposite strand on the RASSF1 gene locus in several cell lines and tissues and binds PRC2. ANRASSF1 is transcribed through RNA polymerase II and is 5'-capped and polyadenylated; it exhibits nuclear localization and has a shorter half-life compared with other lncRNAs that bind PRC2. ANRASSF1 endogenous expression is higher in breast and prostate tumor cell lines compared with non-tumor, and an opposite pattern is observed for RASSF1A. ANRASSF1 ectopic overexpression reduces RASSF1A abundance and increases the proliferation of HeLa cells, whereas ANRASSF1 silencing causes the opposite effects. These changes in ANRASSF1 levels do not affect the RASSF1C isoform abundance. ANRASSF1 overexpression causes a marked increase in both PRC2 occupancy and histone H3K27me3 repressive marks, specifically at the RASSF1A promoter region. No effect of ANRASSF1 overexpression was detected on PRC2 occupancy and histone H3K27me3 at the promoter regions of RASSF1C and the four other neighboring genes, including two well-characterized tumor suppressor genes. Additionally, we demonstrated that ANRASSF1 forms an RNA/DNA hybrid and recruits PRC2 to the RASSF1A promoter. Together, these results demonstrate a novel mechanism of epigenetic repression of the RASSF1A tumor suppressor gene involving antisense unspliced lncRNA, in which ANRASSF1 selectively represses the expression of the RASSF1 isoform overlapping the antisense transcript in a location-specific manner. In a broader perspective, our findings suggest that other non-characterized unspliced intronic lncRNAs transcribed in the human genome might contribute to a location-specific epigenetic modulation of genes.