Experimental validation of the regulated expression of large numbers of non-coding RNAs from the mouse genome

Recent large-scale analyses of mainly full-length cDNA libraries generated from a variety of mouse tissues indicated that almost half of all representative cloned sequences did flat contain ail apparent protein-coding sequence, and were putatively derived from non-protein-coding RNA (ncRNA) genes. However, many of these clones were singletons and the majority were unspliced, raising the possibility that they may be derived from genomic DNA or unprocessed pre-rnRNA contamination during library construction, or alternatively represent nonspecific transcriptional noise. Here we Show, using reverse transcriptase-dependent PCR, microarray, and Northern blot analyses, that many of these clones were derived from genuine transcripts Of unknown function whose expression appears to be regulated. The ncRNA transcripts have larger exons and fewer introns than protein-coding transcripts. Analysis of the genomic landscape around these sequences indicates that some cDNA clones were produced not from terminal poly(A) tracts but internal priming sites within longer transcripts, only a minority of which is encompassed by known genes. A significant proportion of these transcripts exhibit tissue-specific expression patterns, as well as dynamic changes in their expression in macrophages following lipopolysaccharide Stimulation. Taken together, the data provide strong support for the conclusion that ncRNAs are an important, regulated component of the mammalian transcriptome.

Non-coding RNA

The term non-coding RNA (ncRNA) is commonly employed for RNA that does not encode a protein, but this does not mean that such RNAs do not contain information nor have function. Although it has been generally assumed that most genetic information is transacted by proteins, recent evidence suggests that the majority of the genomes of mammals and other complex organisms is in fact transcribed into ncRNAs, many of which are alternatively spliced and/or processed into smaller products. These ncRNAs include microRNAs and snoRNAs (many if not most of which remain to be identified), as well as likely other classes of yet-to-be-discovered small regulatory RNAs, and tens of thousands of longer transcripts (including complex patterns of interlacing and overlapping sense and antisense transcripts), most of whose functions are unknown. These RNAs (including those derived from introns) appear to comprise a hidden layer of internal signals that control various levels of gene expression in physiology and development, including chromatin architecture/epigenetic memory, transcription, RNA splicing, editing, translation and turnover. RNA regulatory networks may determine most of our complex characteristics, play a significant role in disease and constitute an unexplored world of genetic variation both within and between species.

Non-coding RNAs in the nervous system

Increasing evidence suggests that the development and function of the nervous system is heavily dependent on RNA editing and the intricate spatiotemporal expression of a wide repertoire of non-coding RNAs, including micro RNAs, small nucleolar RNAs and longer non-coding RNAs. Non-coding RNAs may provide the key to understanding the multi-tiered links between neural development, nervous system function, and neurological diseases.

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.

Dissecting the Functional Impacts of Non-Coding Genetic Variation

A large proportion of the variation in traits between individuals can be attributed to variation in the nucleotide sequence of the genome. The most commonly studied traits in human genetics are related to disease and disease susceptibility. Although scientists have identified genetic causes for over 4,000 monogenic diseases, the underlying mechanisms of many highly prevalent multifactorial inheritance disorders such as diabetes, obesity, and cardiovascular disease remain largely unknown. Identifying genetic mechanisms for complex traits has been challenging because most of the variants are located outside of protein-coding regions, and determining the effects of such non-coding variants remains difficult. In this dissertation, I evaluate the hypothesis that such non-coding variants contribute to human traits and diseases by altering the regulation of genes rather than the sequence of those genes. I will specifically focus on studies to determine the functional impacts of genetic variation associated with two related complex traits: gestational hyperglycemia and fetal adiposity. At the genomic locus associated with maternal hyperglycemia, we found that genetic variation in regulatory elements altered the expression of the HKDC1 gene. Furthermore, we demonstrated that HKDC1 phosphorylates glucose in vitro and in vivo, thus demonstrating that HKDC1 is a fifth human hexokinase gene. At the fetal-adiposity associated locus, we identified variants that likely alter VEPH1 expression in preadipocytes during differentiation. To make such studies of regulatory variation high-throughput and routine, we developed POP-STARR, a novel high throughput reporter assay that can empirically measure the effects of regulatory variants directly from patient DNA. By combining targeted genome capture technologies with STARR-seq, we assayed thousands of haplotypes from 760 individuals in a single experiment. We subsequently used POP-STARR to identify three key features of regulatory variants: that regulatory variants typically have weak effects on gene expression; that the effects of regulatory variants are often coordinated with respect to disease-risk, suggesting a general mechanism by which the weak effects can together have phenotypic impact; and that nucleotide transversions have larger impacts on enhancer activity than transitions. Together, the findings presented here demonstrate successful strategies for determining the regulatory mechanisms underlying genetic associations with human traits and diseases, and value of doing so for driving novel biological discovery.

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.

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.

Role of non-coding RNA alterations in melanoma

Cutaneous melanoma (CM) is a potentially lethal form of skin cancer and its most important histopathologic factor for staging is Breslow thickness (BT). Its correct determination is fundamental for pathologists. A deeper understanding of the molecular processes guiding CM pathogenesis could improve diagnosis, treatment and prognosis. MicroRNAs (miRNAs) play a key role in CM biology. The firs aim was to investigate miRNA expression in reference to BT assessment. We found that the combined miRNA expression of miR-21-5p and miR-146a-5p above or below 1.5 was significantly associated with overall survival and successfully identified all superficially spreading melanoma (SSM) patients with relapsing suggesting that the combined assessment of these miRNAs expression could aid in SSM staging. Secondly, we focus on multiple primary melanoma (MPM) patients, which develop multiple primary melanomas in their lifetime, and represent a model of high-risk CM occurrence. We explored the miRNome of single CM and MPM: CM and MPM present several dysregulated miRNAs, including key miRNAs involved in epithelial-mesenchymal transition. A different miRNA profile was observed between 1st and 2nd melanoma from the same patient. MiRNA target analysis revealed a more differentiated and less invasive status of MPMs compared to CMs. This characterization of the miRNA regulatory network of MPMs highlights molecular features differentiating this subtype from CM. Recently, NGS experiments revealed the existence of miRNA variants (isomiRs) with different length and sequence. We identified a shorter 3’isoform as tenfold over-represented compared to the canonical form of miR-125a-5p. Target analysis revealed that miRNA shortening could change the pattern of target gene regulation. Finally, we study miRNA and isomiR dysregulation in benign nevi (BN) and CM and in CM and melanoma metastasis. The reported non-random dysregulation of specific isomiRs contributes to the understanding of the complex melanoma pathogenesis and serves as the basis for further functional studies.

Characterisation of RNase Y in H. pylori: a non-essential enzyme involved in the processing of cag-PAI non coding RNA 1 (CncR1) sRNA

The importance of Helicobacter pylori as a human pathogen is underlined by the plethora of diseases it is responsible for. The capacity of H. pylori to adapt to the restricted host-associated environment andto evade the host immune response largely depends on a streamlined signalling network. The peculiar H. pylori small genome size combined with its paucity of transcriptional regulators highlights the relevance of post-transcriptional regulatory mechanisms as small non-coding RNAs (sRNAs). However, among the 8 RNases represented in H. pylori genome, a regulator guiding sRNAs metabolism is still not well studied. We investigated for the first time the physiological role in H. pylori G27 strain of the RNase Y enzyme. In the first line of research we provide a comprehensive characterization of the RNase Y activity by analysing its genomic organization and the factors that orchestrate its expression. Then, based on bioinformatic prediction models, we depict the most relevant determinants of RNase Y function, demonstrating a correlation of both structure and domain organization with orthologues represented in Gram-positive bacteria. To unveil the post-transcriptional regulatory effect exerted by the RNase Y, we compared the transcriptome of an RNase Y knock-out mutant to the parental wild type strain by RNA-seq approach. In the second line of research we characterized the activity of this single strand specific endoribonuclease on cag-PAI non coding RNA 1 (CncR1) sRNA. We found that deletion or inactivation of RNase Y led to the accumulation of a 3’-extended CncR1 (CncR1-L) transcript over time. Moreover, beneath its increased half-life, CncR1-L resembled a CncR1 inactive phenotype. Finally, we focused on the characterization of the in vivo interactome of CncR1. We set up a preliminary MS2-affinity purification coupled with RNA-sequencing (MAPS) approach and we evaluated the enrichment of specific targets, demonstrating the suitability of the technique in the H. pylori G27 strain.

Reverse transcriptase activity and untranslated region sharing of a new RTE-like, non-long terminal repeat retrotransposon from the human blood fluke, Schistosoma japonicum

A new RTE-like, non-long terminal repeat retrotransposon, termed SjR2, from the human blood fluke, Schistosoma japonicum, is described. SjR2 is similar to3.9 kb in length and is constituted of a single open reading frame encoding a polyprotein with apurinic/apyrimidinic endonuclease and reverse transcriptase domains. The open reading frame is bounded by 5'- and 3'-terininal untranslated regions and, at its 3-terminus, SjR2 bears a short (TGAC)(3) repeat. Phylogenetic analyses based on conserved domains of reverse transcriptase or endonuclease revealed that SjR2 belonged to the RTE clade of non-long terminal repeat retrotransposons. Further, SjR2 was homologous, but probably not orthologous, to SR2 front the African blood fluke, Schistosoma mansoni; this RTE-like family of non-long terminal repeat retrotransposons appears to have arisen before the divergence of the extant schistosome species. Hybridisation analyses indicated that similar to 10,000 copies of SjR2 were dispersed throughout the S. japonicum chromosomes, accounting for up to 14% of the nuclear genome. Messenger RNAs encoding the reverse transcriptase and endonuclease domains of SjR2 were detected in several developmental stages of the schistosome, indicating that the retrotransposon was actively replicating within the genome of the parasite. Exploration of the coding and non-coding regions of SjR2 revealed two notable characteristics. First, the recombinant reverse transcriptase domain of SjR2 expressed in insect cells primed reverse transcription of SjR2 mRNA in vitro. By contrast, recombinant SjR2-endonuclease did not appear to cleave schistosome or plasmid DNA. Second, the 5'-untranslated region of SjR2 was >80% identical to the 3-untranslated region of a schistosome heat shock protein-70 gene (hsp-70) in the antisense orientation, indicating that SjR2-like elements were probably inserted into the non-coding regions of ancestral S. japonicum HSP-70, probably after the species diverged from S. mansoni. (C) 2002 Australian Society for Parasitology Inc. Published by Elsevier Science Ltd. All rights reserved.

Identification by Real-time PCR of 13 mature microRNAs differentially expressed in colorectal cancer and non-tumoral tissues

MicroRNAs (miRNAs) are short non-coding RNA molecules playing regulatory roles by repressing translation or cleaving RNA transcripts. Although the number of verified human miRNA is still expanding, only few have been functionally described. However, emerging evidences suggest the potential involvement of altered regulation of miRNA in pathogenesis of cancers and these genes are thought to function as both tumours suppressor and oncogenes. In our study, we examined by Real-Time PCR the expression of 156 mature miRNA in colorectal cancer. The analysis by several bioinformatics algorithms of colorectal tumours and adjacent non-neoplastic tissues from patients and colorectal cancer cell lines allowed identifying a group of 13 miRNA whose expression is significantly altered in this tumor. The most significantly deregulated miRNA being miR-31, miR-96, miR-133b, miR-135b, miR-145, and miR-183. In addition, the expression level of miR-31 was correlated with the stage of CRC tumor. Our results suggest that miRNA expression profile could have relevance to the biological and clinical behavior of colorectal neoplasia.

Numerous potentially functional but non-genic conserved sequences on human chromosome 21.

The use of comparative genomics to infer genome function relies on the understanding of how different components of the genome change over evolutionary time. The aim of such comparative analysis is to identify conserved, functionally transcribed sequences such as protein-coding genes and non-coding RNA genes, and other functional sequences such as regulatory regions, as well as other genomic features. Here, we have compared the entire human chromosome 21 with syntenic regions of the mouse genome, and have identified a large number of conserved blocks of unknown function. Although previous studies have made similar observations, it is unknown whether these conserved sequences are genes or not. Here we present an extensive experimental and computational analysis of human chromosome 21 in an effort to assign function to sequences conserved between human chromosome 21 (ref. 8) and the syntenic mouse regions. Our data support the presence of a large number of potentially functional non-genic sequences, probably regulatory and structural. The integration of the properties of the conserved components of human chromosome 21 to the rapidly accumulating functional data for this chromosome will improve considerably our understanding of the role of sequence conservation in mammalian genomes.