Long non coding RNAs (lncRNAs) are dysregulated in Malignant Pleural Mesothelioma (MPM)

Malignant Pleural Mesothelioma (MPM) is an aggressive cancer that is often diagnosed at an advanced stage and is characterized by a long latency period (20-40 years between initial exposure and diagnosis) and prior exposure to asbestos. Currently accurate diagnosis of MPM is difficult due to the lack of sensitive biomarkers and despite minor improvements in treatment, median survival rates do not exceed 12 months. Accumulating evidence suggests that aberrant expression of long non-coding RNAs (lncRNAs) play an important functional role in cancer biology. LncRNAs are a class of recently discovered non-protein coding RNAs >200 nucleotides in length with a role in regulating transcription. Here we used NCode long noncoding microarrays to identify differentially expressed lncRNAs potentially involved in MPM pathogenesis. High priority candidate lncRNAs were selected on the basis of statistical (P<0.05) and biological significance (>3-fold difference). Expression levels of 9 candidate lncRNAs were technically validated using RT-qPCR, and biologically validated in three independent test sets: (1) 57 archived MPM tissues obtained from extrapleural pneumonectomy patients, (2) 15 cryopreserved MPM and 3 benign pleura, and (3) an extended panel of 10 MPM cell lines. RT-qPCR analysis demonstrated consistent up-regulation of these lncRNAs in independent datasets. ROC curve analysis showed that two candidates were able to separate benign pleura and MPM with high sensitivity and specificity, and were associated with nodal metastases and survival following induction chemotherapy. These results suggest that lncRNAs have potential to serve as biomarkers in MPM.

Identification of a long non-coding RNA gene, GHSROS (growth hormone secretagogue receptor opposite strand), which stimulates cell migration in non-small cell lung cancer cell lines

The molecular mechanisms involved in non‑small cell lung cancer tumourigenesis are largely unknown; however, recent studies have suggested that long non-coding RNAs (lncRNAs) are likely to play a role. In this study, we used public databases to identify an mRNA-like, candidate long non-coding RNA, GHSROS (GHSR opposite strand), transcribed from the antisense strand of the ghrelin receptor gene, growth hormone secretagogue receptor (GHSR). Quantitative real-time RT-PCR revealed higher expression of GHSROS in lung cancer tissue compared to adjacent, non-tumour lung tissue. In common with many long non-coding RNAs, GHSROS is 5' capped and 3' polyadenylated (mRNA-like), lacks an extensive open reading frame and harbours a transposable element. Engineered overexpression of GHSROS stimulated cell migration in the A549 and NCI-H1299 non-small cell lung cancer cell lines, but suppressed cell migration in the Beas-2B normal lung-derived bronchoepithelial cell line. This suggests that GHSROS function may be dependent on the oncogenic context. The identification of GHSROS, which is expressed in lung cancer and stimulates cell migration in lung cancer cell lines, contributes to the growing number of non-coding RNAs that play a role in the regulation of tumourigenesis and metastatic cancer progression.

Non-coding RNAs: multi-tasking molecules in the cell

In the last years it has become increasingly clear that the mammalian transcriptome is highly complex and includes a large number of small non-coding RNAs (sncRNAs) and long noncoding RNAs (lncRNAs). Here we review the biogenesis pathways of the three classes of sncRNAs, namely short interfering RNAs (siRNAs), microRNAs (miRNAs) and PIWI-interacting RNAs (piRNAs). These ncRNAs have been extensively studied and are involved in pathways leading to specific gene silencing and the protection of genomes against virus and transposons, for example. Also, lncRNAs have emerged as pivotal molecules for the transcriptional and post-transcriptional regulation of gene expression which is supported by their tissue-specific expression patterns, subcellular distribution, and developmental regulation. Therefore, we also focus our attention on their role in differentiation and development. SncRNAs and lncRNAs play critical roles in defining DNA methylation patterns, as well as chromatin remodeling thus having a substantial effect in epigenetics. The identification of some overlaps in their biogenesis pathways and functional roles raises the hypothesis that these molecules play concerted functions in vivo, creating complex regulatory networks where cooperation with regulatory proteins is necessary. We also highlighted the implications of biogenesis and gene expression deregulation of sncRNAs and lncRNAs in human diseases like cancer.

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.

Transcriptional dynamics of the human DMD locus

The human DMD locus encodes dystrophin protein. Absence or reduced levels of dystrophin (DMD or BMD phenotype, respectively) lead to progressive muscle wasting. Little is known about the complex coordination of dystrophin expression and its transcriptional regulation is a field of intense interest. In this work we found that DMD locus harbours multiple long non coding RNAs which orchestrate and control transcription of muscle dystrophin mRNA isoforms. These lncRNAs are tissue-specific and highly expressed during myogenesis, suggesting a possible role in tissue-specific expression of DMD gene isoforms. Their forced ectopic expression in human muscle and neuronal cells leads to a specific and negative regulation of endogenous dystrophin full lenght isoforms. An intriguing aspect regarding the transcription of the DMD locus is the gene size (2.4Mb). The mechanism that ensures the complete synthesis of the primary transcript and the coordinated splicing of 79 exons is still completely unknown. By ChIP-on-chip analyses, we discovered novel regions never been involved before in the transcription regulation of the DMD locus. Specifically, we observed enrichments for Pol II, P-Ser2, P-Ser5, Ac-H3 and 2Me-H3K4 in an intronic region of 3Kb (approximately 21Kb) downstream of the end of DMD exon 52 and in a region of 4Kb spanning the DMD exon 62. Interestingly, this latter region and the TSS of Dp71 are strongly marked by 3Me-H3K36, an histone modification associated with the regulation of splicing process. Furthermore, we also observed strong presence of open chromatin marks (Ac-H3 and 2Me-H3K4) around intron 34 and the exon 45 without presence of RNA pol II. We speculate that these two regions may exert an enhancer-like function on Dp427m promoter, although further investigations are necessary. Finally, we investigated the nuclear-cytoplasmic compartmentalization of the muscular dystrophin mRNA and, specifically, we verified whether the exon skipping therapy could influence its cellular distribution.

A long non-coding RNA targets microRNA miR-34a to regulate colon cancer stem cell asymmetric division.

The roles of long non-coding RNAs (lncRNAs) in regulating cancer and stem cells are being increasingly appreciated. Its diverse mechanisms provide the regulatory network with a bigger repertoire to increase complexity. Here we report a novel LncRNA, Lnc34a, that is enriched in colon cancer stem cells (CCSCs) and initiates asymmetric division by directly targeting the microRNA miR-34a to cause its spatial imbalance. Lnc34a recruits Dnmt3a via PHB2 and HDAC1 to methylate and deacetylate the miR-34a promoter simultaneously, hence epigenetically silencing miR-34a expression independent of its upstream regulator, p53. Lnc34a levels affect CCSC self-renewal and colorectal cancer (CRC) growth in xenograft models. Lnc34a is upregulated in late-stage CRCs, contributing to epigenetic miR-34a silencing and CRC proliferation. The fact that lncRNA targets microRNA highlights the regulatory complexity of non-coding RNAs (ncRNAs), which occupy the bulk of the genome.

Investigation of the role of long non-coding RNAs in oncogene induced cellular senescence

Cellular senescence is a stable arrest of cell proliferation induced by several factors such as activated oncogenes, oxidative stress and shortening of telomeres. Senescence acts as a tumour suppression mechanism to halt the progression of cancer. However, senescence may also impact negatively upon tissue regeneration, thus contributing to the effects of ageing. The eukaryotic genome is controlled by various modes of transcriptional and translational regulation. Focus has therefore centred on the role of long non- coding RNAs (lncRNAs) in regulating the genome. Accordingly, understanding how lncRNAs function to regulate the senescent genome is integral to improving our knowledge and understanding of tumour suppression and ageing. Within this study, I set out to investigate the expression of lncRNAs’ expression within models of senescence. Through a custom expression array, I have shown that expression of multiple different lncRNAs is up-regulated and down regulated in IMR90 replicative senescent fibroblasts and oncogene-induced senescent melanocytes. LncRNA expression was determined to be specific to stable senescence-associated cell arrest and predominantly within the nucleus of senescent cells. In order to examine the function of lncRNA expression in senescence, I selected lncRNA transcript ENST0000430998 (lncRNA_98) to focus my investigations upon. LncRNA_98 was robustly upregulated within multiple models of senescence and efficiently depleted using anti-sense oligonucleotide technology. Characterisation and unbiased RNA-sequencing of lncRNA_98 deficient senescent cells highlighted a list of genes that are regulated by lncRNA_98 expression in senescent cells and may regulate aspects of the senescence program. Specifically, the formation of SAHF was impeded upon depletion of lncRNA_98 expression and levels of total pRB protein expression severely decreased. Validation and recapitulation of consequences of pRB depletion was confirmed through lncRNA_98 knock-out cells generated using CRISPR technology. Surprisingly, inhibition of ATM kinase functions permitted the restoration of pRB protein levels within lncRNA_98 deficient cells. I propose that lncRNA_98 antagonizes the ability of ATM kinase to downregulate pRB expression at a post-transcriptional level, thereby potentiating senescence. Furthermore, lncRNA expression was detected within fibroblasts of old individuals and visualised within senescent melanocytes in human benign nevi, a barrier to melanoma progression. Conversely, mining of 337 TCGA primary melanoma data sets highlighted that the lncRNA_98 gene and its expression was lost from a significant proportion of melanoma samples, consistent with lncRNA_98 having a tumour suppressor functions. The data presented in this study illustrates that lncRNA_98 expression has a regulatory role over pRB expression in senescence and may regulate aspects of tumourigenesis and ageing.

Differential expression in “rogue” paramutation in peas (Pisum sativum L.)

Dissertação de Mestrado, Engenharia Biológica, Faculdade de Ciências e Tecnologia, Universidade do Algarve, 2014

Epigenetics and alternative splicing

Dissertação de Mestrado, Oncobiologia: Mecanismos Moleculares do Cancro, Departamento de Ciências Biomédicas e Medicina, Universidade do Algarve, 2015

The role of non-coding RNA in the development of pulmonary arterial hypertension

Pulmonary arterial hypertension (PAH) is a progressive disease of the small pulmonary arteries, characterised by pulmonary vascular remodelling due to excessive proliferation and resistance to apoptosis of pulmonary artery endothelial cells (PAECs) and pulmonary artery smooth muscle cells (PASMCs). The increased pulmonary vascular resistance and elevated pulmonary artery pressures result in right heart failure and premature death. Germline mutations of the bone morphogenetic protein receptor-2 (bmpr2) gene, a receptor of the transforming growth factor beta (TGF-β) superfamily, account for approximately 75%-80% of the cases of heritable form of PAH (HPAH) and 20% of sporadic cases or idiopathic PAH (IPAH). IPAH patients without known bmpr2 mutations show reduced expression of BMPR2. However only ~ 20% of bmpr2-mutation carriers will develop the disease, due to an incomplete penetrance, thus the need for a ‘second hit’ including other genetic and/or environmental factors is accepted. Diagnosis of PAH occurs most frequently when patients have reached an advanced stage of disease. Although modern PAH therapies can markedly improve a patient’s symptoms and slow the rate of clinical deterioration, the mortality rate from PAH remains unacceptably high. Therefore, the development of novel therapeutic approaches is required for the treatment of this multifaceted disease. Noncoding RNAs (ncRNAs) include microRNAs (miRNAs) and long noncoding RNAs (lncRNAs). MiRNAs are ~ 22 nucleotide long and act as negative regulators of gene ex-pression via degradation or translational inhibition of their target mRNAs. Previous studies showed extensive evidence for the role of miRNAs in the development of PAH. LncRNAs are transcribed RNA molecules greater than 200 nucleotides in length. Similar to classical mRNA, lncRNAs are translated by RNA polymerase II and are generally alternatively spliced and polyadenylated. LncRNAs are highly versatile and function to regulate gene expression by diverse mechanisms. Unlike miRNAs, which exhibit well-defined actions in negatively regulating gene expression via the 3’-UTR of mRNAs, lncRNAs play more diverse and unpredictable regulatory roles. Although a number of lncRNAs have been intensively investigated in the cancer field, studies of the role of lncRNAs in vascular diseases such as PAH are still at a very early stage. The aim of this study was to investigate the involvement of specific ncRNAs in the development of PAH using experimental animal models and cell culture. The first ncRNA we focused on was miR-143, which is up-regulated in the lung and right ventricle tissues of various animal models of PH, as well as in the lungs and PASMCs of PAH patients. We show that genetic ablation of miR-143 is protective against the development of chronic hypoxia induced PH in mice, assessed via measurement of right ventricular systolic pressure (RVSP), right ventricular hypertrophy (RVH) and pulmonary vascular remodelling. We further report that knockdown of miR-143-3p in WT mice via anti-miR-143-3p administration prior to exposure of mice to chronic hypoxia significantly decreases certain indices of PH (RVSP) although no significant changes in RVH and pulmo-nary vascular remodelling were observed. However, a reversal study using antimiR-143-3p treatment to modulate miR-143-3p demonstrated a protective effect on RVSP, RVH, and muscularisation of pulmonary arteries in the mouse chronic hypoxia induced PH model. In vitro experiments showed that miR-143-3p overexpression promotes PASMC migration and inhibits PASMC apoptosis, while knockdown miR-143-3p elicits the opposite effect, with no effects observed on cellular proliferation. Interestingly, miR-143-3p-enriched exosomes derived from PASMCs mediated cell-to-cell communication between PASMCs and PAECs, contributing to the pro-migratory and pro-angiogenic phenotype of PAECs that underlies the pathogenesis of PAH. Previous work has shown that miR-145-5p expression is upregulated in the chronic hypoxia induced mouse model of PH, as well as in PAH patients. Genetic ablation and pharmacological inhibition (subcutaneous injection) of miR-145-5p exert a protective against the de-velopment of PAH. In order to explore the potential for alternative, more lung targeted delivery strategies, miR-145-5p expression was inhibited in WT mice using intranasal-delivered antimiR-145-5p both prior to and post exposure to chronic hypoxia. The decreased expression of miR-145-5p in lung showed no beneficial effect on the development of PH compared with control antimiRNA treated mice exposed to chronic hypoxia. Thus, miR-143-3p modulated both cellular and exosome-mediated responses in pulmonary vascular cells, while the inhibition of miR-143-3p prevented the development of experimental pulmonary hypertension. We focused on two lncRNAs in this project: Myocardin-induced Smooth Muscle Long noncoding RNA, Inducer of Differentiation (MYOSLID) and non-annotated Myolnc16, which were identified from RNA sequencing studies in human coronary artery smooth muscle cells (HCASMCs) that overexpress myocardin. MYOSLID was significantly in-creased in PASMCs from patients with IPAH compared to healthy controls and increased in circulating endothelial progenitor cells (EPCs) from bmpr2 mutant PAH patients. Exposure of PASMCs to hypoxia in vitro led to a significant upregulation in MYOSLID expres-sion. MYOSLID expression was also induced by treatment of PASMC with BMP4, TGF-β and PDGF, which are known to be triggers of PAH in vitro. Small interfering RNA (siR-NA)-mediated knockdown MYOSLID inhibited migration and induced cell apoptosis without affecting cell proliferation and upregulated several genes in the BMP pathway in-cluding bmpr1α, bmpr2, id1, and id3. Modulation of MYOSLID also affected expression of BMPR2 at the protein level. In addition, MYOSLID knockdown affected the BMP-Smad and BMP-non-Smad signalling pathways in PASMCs assessed by phosphorylation of Smad1/5/9 and ERK1/2, respectively. In PAECs, MYOSLID expression was also induced by hypoxia exposure, VEGF and FGF2 treatment. In addition, MYOSLID knockdown sig-nificantly decreased the proliferation of PAECs. Thus, MYOSLID may be a novel modulator in pulmonary vascular cell functions, likely through the BMP-Smad and –non-Smad pathways. Treatment of PASMCs with inflammatory cytokines (IL-1 and TNF-α) significantly in-duced the expression of Myolnc16 at a very early time point. Knockdown of Myolnc16 in vitro decreased the expression of il-6, and upregulated the expression of il-1 and il-8 in PASMCs. Moreover, the expression levels of chemokines (cxcl1, cxcl6 and cxcl8) were sig-nificantly decreased with Myolnc16 knockdown. In addition, Myolnc16 knockdown decreased the MAP kinase signalling pathway assessed by phosphorylation of ERK1/2 and p38 MAPK and inhibited cell migration and proliferation in PASMCs. Thus, Myolnc16 may a novel modulator of PASMCs functions through anti-inflammatory signalling pathways. In summary, in this thesis we have demonstrated how miR-143-3p plays a protective role in the development of PH both in vivo animal models and patients, as well as in vitro cell cul-ture. Moreover, we have showed the role of two novel lncRNAs in pulmonary vascular cells. These ncRNAs represent potential novel therapeutic targets for the treatment of PAH with further work addressing to investigate the target genes, and the pathways modulated by these ncRNAs during the development of PAH.