Genome mapping and expression analyses of human intronic noncoding RNAs reveal tissue-specific patterns and enrichment in genes related to regulation of transcription

Abstract Background RNAs transcribed from intronic regions of genes are involved in a number of processes related to post-transcriptional control of gene expression. However, the complement of human genes in which introns are transcribed, and the number of intronic transcriptional units and their tissue expression patterns are not known. Results A survey of mRNA and EST public databases revealed more than 55,000 totally intronic noncoding (TIN) RNAs transcribed from the introns of 74% of all unique RefSeq genes. Guided by this information, we designed an oligoarray platform containing sense and antisense probes for each of 7,135 randomly selected TIN transcripts plus the corresponding protein-coding genes. We identified exonic and intronic tissue-specific expression signatures for human liver, prostate and kidney. The most highly expressed antisense TIN RNAs were transcribed from introns of protein-coding genes significantly enriched (p = 0.002 to 0.022) in the 'Regulation of transcription' Gene Ontology category. RNA polymerase II inhibition resulted in increased expression of a fraction of intronic RNAs in cell cultures, suggesting that other RNA polymerases may be involved in their biosynthesis. Members of a subset of intronic and protein-coding signatures transcribed from the same genomic loci have correlated expression patterns, suggesting that intronic RNAs regulate the abundance or the pattern of exon usage in protein-coding messages. Conclusion We have identified diverse intronic RNA expression patterns, pointing to distinct regulatory roles. This gene-oriented approach, using a combined intron-exon oligoarray, should permit further comparative analysis of intronic transcription under various physiological and pathological conditions, thus advancing current knowledge about the biological functions of these noncoding RNAs.

Identification of protein-coding and intronic noncoding RNAs down-regulated in clear cell renal carcinoma

The clear cell subtype of renal cell carcinoma (RCC) is the most lethal and prevalent cancer of the urinary system. To investigate the molecular changes associated with malignant transformation in clear cell RCC, the gene expression profiles of matched samples of tumor and adjacent non-neoplastic tissue were obtained from six patients. A custom-built cDNA microarray platform was used, comprising 2292 probes that map to exons of genes and 822 probes for noncoding RNAs mapping to intronic regions. Intronic transcription was detected in all normal and neoplastic renal tissues. A subset of 55 transcripts was significantly down-regulated in clear cell RCC relative to the matched nontumor tissue as determined by a combination of two statistical tests and leave-one-out patient cross-validation. Among the down-regulated transcripts, 49 mapped to untranslated or coding exons and 6 were intronic relative to known exons of protein-coding genes. Lower levels of expression of SIN3B, TRIP3, SYNJ2BP and NDE1 (P<0.02), and of intronic transcripts derived from SND1 and ACTN4 loci (P<0.05), were confirmed in clear cell RCC by Real-time RT-PCR. A subset of 25 transcripts was deregulated in additional six nonclear cell RCC samples, pointing to common transcriptional alterations in RCC irrespective of the histological subtype or differentiation state of the tumor. Our results indicate a novel set of tumor suppressor gene candidates, including noncoding intronic RNAs, which may play a significant role in malignant transformations of normal renal cells. (C) 2008 Wiley-Liss, Inc.

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.

A highly conserved c-fms gene intronic element controls macrophage-specific and regulated expression

The c fins gene encodes the receptor for macrophage colony-stimulating factor-1. This gene is expressed selectively in the macrophage cell lineage. Previous studies have implicated sequences in intron 2 that control transcript elongation in tissue-specific and regulated expression of c -fms. Four macrophage-specific deoxyribonuclease I (DNase I)-hypersensitive sites (DHSS) were identified within mouse intron 2. Sequences of these DHSS were found to be highly conserved compared with those in the human gene. A 250-bp region we refer to as the fins intronic regulatory element (FIRE), which is even more highly conserved than the c-fins proximal promoter, contains many consensus binding sites for macrophage-expressed transcription factors including Spl, PU.1, and C/EBP. FIRE was found to act as a macrophage-specific enhancer and as a promoter with an antisense orientation preference in transient transfections. In stable transfections of the macrophage line RAW264, as well as in clones selected for high and low-level c -fms mRNA expression, the presence of intron 2 increased the frequency and level of expression of reporter genes compared with those attained using the promoter alone. Removal of FIRE abolished reporter gene expression, revealing a suppressive activity in the remaining intronic sequences. Hence, FIRE is shown to be a key regulatory element in the fins gene.

Conserved tissue expression signatures of intronic noncoding RNAs transcribed from human and mouse loci

It has been postulated that noncoding RNAs (ncRNAs) are involved in the posttranscriptional control of gene expression, and may have contributed to the emergence of the complex attributes observed in mammalians. We show here that the complement of ncRNAs expressed from intronic regions of the human and mouse genomes comprises at least 78,147 and 39,660 transcriptional units, respectively. To identify conserved intronic sequences expressed in both humans and mice, we used custom-designed human cDNA microarrays to separately interrogate RNA from mouse and human liver, kidney, and prostate tissues. An overlapping tissue expression signature was detected for both species, comprising 198 transcripts; among these, 22 RNAs map to intronic regions with evidence of evolutionary conservation in humans and mice. Transcription of selected human-mouse intronic ncRNAs was confirmed using strand-specific RT-PCR. Altogether, these results support an evolutionarily conserved role of intronic ncRNAs in human and mouse, which are likely to be involved in the fine tuning of gene expression regulation in different mammalian tissues. (C) 2008 Elsevier Inc. All rights reserved.

Androgen responsive intronic non-coding RNAs

Abstract Background Transcription of large numbers of non-coding RNAs originating from intronic regions of human genes has been recently reported, but mechanisms governing their biosynthesis and biological functions are largely unknown. In this work, we evaluated the existence of a common mechanism of transcription regulation shared by protein-coding mRNAs and intronic RNAs by measuring the effect of androgen on the transcriptional profile of a prostate cancer cell line. Results Using a custom-built cDNA microarray enriched in intronic transcribed sequences, we found 39 intronic non-coding RNAs for which levels were significantly regulated by androgen exposure. Orientation-specific reverse transcription-PCR indicated that 10 of the 13 were transcribed in the antisense direction. These transcripts are long (0.5–5 kb), unspliced and apparently do not code for proteins. Interestingly, we found that the relative levels of androgen-regulated intronic transcripts could be correlated with the levels of the corresponding protein-coding gene (asGAS6 and asDNAJC3) or with the alternative usage of exons (asKDELR2 and asITGA6) in the corresponding protein-coding transcripts. Binding of the androgen receptor to a putative regulatory region upstream from asMYO5A, an androgen-regulated antisense intronic transcript, was confirmed by chromatin immunoprecipitation. Conclusion Altogether, these results indicate that at least a fraction of naturally transcribed intronic non-coding RNAs may be regulated by common physiological signals such as hormones, and further corroborate the notion that the intronic complement of the transcriptome play functional roles in the human gene-expression program.

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.

Transcription factor SEF: Purification and functional characterization

Cytokine-induced transcription of the serum amyloid A3 (SAA3) gene promoter requires a transcriptional enhancer that contains three functional elements: two C/EBP-binding sites and a third site that interacts with a constitutively expressed transcription factor, SAA3 enhancer factor (SEF). Deletion or site-specific mutations in the SEF-binding site drastically reduced SAA3 promoter activity, strongly suggesting that SEF is important in SAA3 promoter function. To further elucidate its role in the regulation of the SAA3 gene, we purified SEF from HeLa cell nuclear extracts to near homogeneity by using conventional liquid chromatography and DNA-affinity chromatography. Ultraviolet cross-linking and Southwestern experiments indicated that SEF consisted of a single polypeptide with an apparent molecular mass of 65 kDa. Protein sequencing, oligonucleotide competition and antibody supershift experiments identified SEF as transcription factor LBP-1c/CP2/LSF. Cotransfection of SEF expression plasmid with SAA3-luciferase reporter resulted in 3- to 5-fold activation of SAA3 promoter. Interestingly, when SEF-transfected cells were treated with either conditioned medium (CM) or interleukin (IL) 1, the SAA3 promoter was synergistically activated in a dose-dependent manner. Furthermore, when SEF-binding site was mutated, the response of SAA3 promoter to IL-1 or CM stimulation was abolished or drastically decreased, suggesting that SEF may functionally cooperate with an IL-1-inducible transcription factor. Indeed, our functional studies showed that NFκB is a key transcription factor that mediates the IL-1-induced expression of SAA3 gene, and that SEF can synergize with NFκBp65 to activate SAA3 promoter. By coimmunoprecipitation experiments, we found that SEF could specifically interact with NFκBp65, and that the association of these two factors was enhanced upon IL-1 and CM stimulation. This suggests that the molecular basis for the functional synergy between SEF and NFκB may be due to the ability of SEF to physically interact with NPκB. In addition to its interaction with SEF, NFκB-dependent activation also requires the weak κB site in the C element and its interaction with C/EBP. Besides its role in regulating SAA3 gene expression, we provide evidence that SEF could also bind in a sequence-specific manner to the promoters of α2-macroglobulin, Aα fibrinogen, and 6–16 genes and to an intronic enhancer of the human Wilm's tumor 1 gene, suggesting a functional role in the regulation of these genes. By coimmunoprecipitation experiments, we determined that SEF could specifically associate with both Stat3 and Stat2 upon cytokine stimulation. To examine the functional roles of such interactions, we evaluated the effects of SEF on the transcriptional regulation of two reporter genes: Aα fibrinogen and 6–16, which are IL-6- and interferon-α-responsive, respectively. Our results showed that cotransfection of SEF expression plasmid can activate the expression of Aα fibrinogen gene and 6–16 gene. Moreover, SEF can dramatically enhance the interferon-α-induced expression of 6–16 gene and IL-6-induced expression of Aα fibrinogen gene, suggesting that SEF may functionally cooperate with ISGF3 and Stat3 to mediate interferon-α and IL-6 signaling. ^ Our findings that SEF can interact with multiple cytokine-inducible transcription factors to mediate the expression of target genes open a new avenue of investigation of cooperative transcriptional regulation of gene expression, and should further our understanding of differential gene expression in response to a specific stimulus. In summary, our data provide evidence that SEF can mediate the signaling of different cytokines by interacting with various cytokine-inducible transcription factors. ^

An Intronic Enhancer Is Required for Deflagellation-induced Transcriptional Regulation of a Chlamydomonas reinhardtii Dynein Gene

Chlamydomonas reinhardtii flagellar regeneration is accompanied by rapid induction of genes encoding a large set of flagellar structural components and provides a model system to study coordinate gene regulation and organelle assembly. After deflagellation, the abundance of a 70-kDa flagellar dynein intermediate chain (IC70, encoded by ODA6) mRNA increases approximately fourfold within 40 min and returns to predeflagellation levels by ∼90 min. We show by nuclear run-on that this increase results, in part, from increased rates of transcription. To localize cis induction elements, we created an IC70 minigene and measured accumulation, in C. reinhardtii, of transcripts from the endogenous gene and from introduced promoter deletion constructs. Clones containing 416 base pairs (bp) of 5′- and 2 kilobases (kb) of 3′-flanking region retained all sequences necessary for a normal pattern of mRNA abundance change after deflagellation. Extensive 5′- and 3′- flanking region deletions, which removed multiple copies of a proposed deflagellation-response element (the tub box), did not eliminate induction, and the IC70 5′-flanking region alone did not confer deflagellation responsiveness to a promoterless arylsulfatase (ARS) gene. Instead, an intron in the IC70 gene 5′-untranslated region was found to contain the deflagellation response element. These results suggest that the tub box does not play an essential role in deflagellation-induced transcriptional regulation of this dynein gene.

Ultraconserved elements in insect genomes: A highly conserved intronic sequence implicated in the control of homothorax mRNA splicing

Recently, we identified a large number of ultraconserved (uc) sequences in noncoding regions of human, mouse, and rat genomes that appear to be essential for vertebrate and amniote ontogeny. Here, we used similar methods to identify ultraconserved genomic regions between the insect species Drosophila melanogaster and Drosophila pseudoobscura, as well as the more distantly related Anopheles gambiae. As with vertebrates, ultraconserved sequences in insects appear to Occur primarily in intergenic and intronic sequences, and at intron-exon junctions. The sequences are significantly associated with genes encoding developmental regulators and transcription factors, but are less frequent and are smaller in size than in vertebrates. The longest identical, nongapped orthologous match between the three genomes was found within the homothorax (hth) gene. This sequence spans an internal exon-intron junction, with the majority located within the intron, and is predicted to form a highly stable stem-loop RNA structure. Real-time quantitative PCR analysis of different hth splice isoforms and Northern blotting showed that the conserved element is associated with a high incidence of intron retention in hth pre-mRNA, suggesting that the conserved intronic element is critically important in the post-transcriptional regulation of hth expression in Diptera.

Tshr Intronic Polymorphisms (rs179247 And Rs12885526) And Their Role In The Susceptibility Of The Brazilian Population To Graves' Disease And Graves' Ophthalmopathy.

Intronic thyroid-stimulating hormone receptor polymorphisms have been associated with the risk for both Graves' disease and Graves' ophthalmopathy, but results have been inconsistent among different populations. We aimed to investigate the influence of thyroid-stimulating hormone receptor intronic polymorphisms in a large well-characterized population of GD patients. We studied 279 Graves' disease patients (231 females and 48 males, 39.80 ± 11.69 years old), including 144 with Graves' ophthalmopathy, matched to 296 healthy control individuals. Thyroid-stimulating hormone receptor genotypes of rs179247 and rs12885526 were determined by Real Time PCR TaqMan(®) SNP Genotyping. A multivariate analysis showed that the inheritance of the thyroid-stimulating hormone receptor AA genotype for rs179247 increased the risk for Graves' disease (OR = 2.821; 95 % CI 1.595-4.990; p = 0.0004), whereas the thyroid-stimulating hormone receptor GG genotype for rs12885526 increased the risk for Graves' ophthalmopathy (OR = 2.940; 95 % CI 1.320-6.548; p = 0.0083). Individuals with Graves' ophthalmopathy also presented lower mean thyrotropin receptor antibodies levels (96.3 ± 143.9 U/L) than individuals without Graves' ophthalmopathy (98.3 ± 201.9 U/L). We did not find any association between the investigated polymorphisms and patients clinical features or outcome. We demonstrate that thyroid-stimulating hormone receptor intronic polymorphisms are associated with the susceptibility to Graves' disease and Graves' ophthalmopathy in the Brazilian population, but do not appear to influence the disease course.

Essential regulatory elements for NHE3 gene transcription in renal proximal tubule cells

The objectives of the present study were to identify the cis-elements of the promoter absolutely required for the efficient rat NHE3 gene transcription and to locate positive and negative regulatory elements in the 5’-flanking sequence (5’FS), which might modulate the gene expression in proximal tubules, and to compare this result to those reported for intestinal cell lines. We analyzed the promoter activity of different 5’FS segments of the rat NHE3 gene, in the OKP renal proximal tubule cell line by measuring the activity of the reporter gene luciferase. Because the segment spanning the first 157 bp of 5’FS was the most active it was studied in more detail by sequential deletions, point mutations, and gel shift assays. The essential elements for gene transcription are in the region -85 to -33, where we can identify consensual binding sites for Sp1 and EGR-1, which are relevant to NHE3 gene basal transcription. Although a low level of transcription is still possible when the first 25 bp of the 5’FS are used as promoter, efficient transcription only occurs with 44 bp of 5’FS. There are negative regulatory elements in the segments spanning -1196 to -889 and -467 to -152, and positive enhancers between -889 and -479 bp of 5’FS. Transcription factors in the OKP cell nuclear extract efficiently bound to DNA elements of rat NHE3 promoter as demonstrated by gel shift assays, suggesting a high level of similarity between transcription factors of both species, including Sp1 and EGR-1.

Reduced transcription of TCOF1 in adult cells of Treacher Collins syndrome patients

Background: Treacher Collins syndrome (TCS) is an autosomal dominant craniofacial disorder caused by frameshift deletions or duplications in the TCOF1 gene. These mutations cause premature termination codons, which are predicted to lead to mRNA degradation by nonsense mediated mRNA decay (NMD). Haploinsufficiency of the gene product (treacle) during embryonic development is the proposed molecular mechanism underlying TCS. However, it is still unknown if TCOF1 expression levels are decreased in postembryonic human cells. Methods: We have estimated TCOF1 transcript levels through real time PCR in mRNA obtained from leucocytes and mesenchymal cells of TCS patients (n = 23) and controls (n = 18). Mutational screening and analysis of NMD were performed by direct sequencing of gDNA and cDNA, respectively. Results: All the 23 patients had typical clinical features of the syndrome and pathogenic mutations were detected in 19 of them. We demonstrated that the expression level of TCOF1 is 18-31% lower in patients than in controls (p < 0.05), even if we exclude the patients in whom we did not detect the pathogenic mutation. We also observed that the mutant allele is usually less abundant than the wild type one in mesenchymal cells. Conclusions: This is the first study to report decreased expression levels of TCOF1 in TCS adult human cells, but it is still unknown if this finding is associated to any phenotype in adulthood. In addition, as we demonstrated that alleles harboring the pathogenic mutations have lower expression, we herein corroborate the current hypothesis of NMD of the mutant transcript as the explanation for diminished levels of TCOF1 expression. Further, considering that TCOF1 deficiency in adult cells could be associated to pathologic clinical findings, it will be important to verify if TCS patients have an impairment in adult stem cell properties, as this can reduce the efficiency of plastic surgery results during rehabilitation of these patients.

Increased Levels of NOTCH1, NF-kappa B, and Other Interconnected Transcription Factors Characterize Primitive Sets of Hematopoietic Stem Cells

As previously shown, higher levels of NOTCH1 and increased NF-kappa B signaling is a distinctive feature of the more primitive umbilical cord blood (UCB) CD34+ hematopoietic stem cells (HSCs), as compared to bone marrow ( BM). Differences between BM and UCB cell composition also account for this finding. The CD133 marker defines a more primitive cell subset among CD34+ HSC with a proposed hemangioblast potential. To further evaluate the molecular basis related to the more primitive characteristics of UCB and CD133+ HSC, immunomagnetically purified human CD34+ and CD133+ cells from BM and UCB were used on gene expression microarrays studies. UCB CD34+ cells contained a significantly higher proportion of CD133+ cells than BM (70% and 40%, respectively). Cluster analysis showed that BM CD133+ cells grouped with the UCB cells ( CD133+ and CD34+) rather than to BM CD34+ cells. Compared with CD34+ cells, CD133+ had a higher expression of many transcription factors (TFs). Promoter analysis on all these TF genes revealed a significantly higher frequency ( than expected by chance) of NF-kappa B-binding sites (BS), including potentially novel NF-kappa B targets such as RUNX1, GATA3, and USF1. Selected transcripts of TF related to primitive hematopoiesis and self-renewal, such as RUNX1, GATA3, USF1, TAL1, HOXA9, HOXB4, NOTCH1, RELB, and NFKB2 were evaluated by real-time PCR and were all significantly positively correlated. Taken together, our data indicate the existence of an interconnected transcriptional network characterized by higher levels of NOTCH1, NF-kappa B, and other important TFs on more primitive HSC sets.