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.

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.

Footprints of a trypanosomatid RNA world: pre-small subunit rRNA processing by spliced leader addition trans-splicing

The addition of a capped mini-exon [spliced leader (SL)] through trans-splicing is essential for the maturation of RNA polymerase (pol) II-transcribed polycistronic pre-mRNAs in all members of the Trypanosomatidae family. This process is an inter-molecular splicing reaction that follows the same basic rules of cis-splicing reactions. In this study, we demonstrated that mini-exons were added to precursor ribosomal RNA (pre-rRNA) are transcribed by RNA pol I, including the 5' external transcribed spacer (ETS) region. Additionally, we detected the SL-5' ETS molecule using three distinct methods and located the acceptor site between two known 5' ETS rRNA processing sites (A' and A1) in four different trypanosomatids. Moreover, we detected a polyadenylated 5' ETS upstream of the trans-splicing acceptor site, which also occurs in pre-mRNA trans-splicing. After treatment with an indirect trans-splicing inhibitor (sinefungin), we observed SL-5' ETS decay. However, treatment with 5-fluorouracil (a precursor of RNA synthesis that inhibits the degradation of pre-rRNA) led to the accumulation of SL-5' ETS, suggesting that the molecule may play a role in rRNA degradation. The detection of trans-splicing in these molecules may indicate broad RNA-joining properties, regardless of the polymerase used for transcription.

rpb2 is a reliable reference gene for quantitative gene expression analysis in the dermatophyte Trichophyton rubrum

The selection of reference genes used for data normalization to quantify gene expression by real-time PCR amplifications (qRT-PCR) is crucial for the accuracy of this technique. In spite of this, little information regarding such genes for qRT-PCR is available for gene expression analyses in pathogenic fungi. Thus, we investigated the suitability of eight candidate reference genes in isolates of the human dermatophyte Trichophyton rubrum subjected to several environmental challenges, such as drug exposure, interaction with human nail and skin, and heat stress. The stability of these genes was determined by geNorm, NormFinder and Best-Keeper programs. The gene with the most stable expression in the majority of the conditions tested was rpb2 (DNA-dependent RNA polymerase II), which was validated in three T. rubrum strains. Moreover, the combination of rpb2 and chs1 (chitin synthase) genes provided for the most reliable qRT-PCR data normalization in T. rubrum under a broad range of biological conditions. To the best of our knowledge this is the first report on the selection of reference genes for qRT-PCR data normalization in dermatophytes and the results of these studies should permit further analysis of gene expression under several experimental conditions, with improved accuracy and reliability.

Characterization of anti-silencing factor 1 in Leishmania major

Anti-silencing factor 1 (ASF1) is a histone chaperone that contributes to the histone deposition during nucleosome assembly in newly replicated DNA. It is involved in chromatin disassembly, transcription activation and in the cellular response to DNA damage. In Leishmania major the ASF1 gene (LmASF1) is located in chromosome 20 and codes for a protein showing 67% of identity with the Trypanosoma brucei TbASF1a. Compared to orthologous proteins, LmASF1 conserves the main residues relevant for its various biological functions. To study ASF1 in Leishmania we generated a mutant overexpressing LmASF1 in L. major. We observed that the excess of LmASF1 impaired promastigotes growth rates and had no impact on cell cycle progress. Differently from yeast, ASF1 overproduction in Leishmania did not affect expression levels of genes located on telomeres, but led to an upregulation of proteins involved in chromatin remodelling and physiological stress, such as heat shock proteins, oxidoreductase activity and proteolysis. In addition, we observed that LmASF1 mutant is more susceptible to the DNA damaging agent, methyl methane sulphonate, than the control line. Therefore, our study suggests that ASF1 from Leishmania pertains to the chromatin remodelling machinery of the parasite and acts on its response to DNA damage.

Transcriptional profiling of Saccharomyces cerevisiae exposed to propolis

Abstract Background Propolis is a natural product of plant resins collected by honeybees (Apis mellifera) from various plant sources. Our previous studies indicated that propolis sensitivity is dependent on the mitochondrial function and that vacuolar acidification and autophagy are important for yeast cell death caused by propolis. Here, we extended our understanding of propolis-mediated cell death in the yeast Saccharomyces cerevisiae by applying systems biology tools to analyze the transcriptional profiling of cells exposed to propolis. Methods We have used transcriptional profiling of S. cerevisiae exposed to propolis. We validated our findings by using real-time PCR of selected genes. Systems biology tools (physical protein-protein interaction [PPPI] network) were applied to analyse the propolis-induced transcriptional bevavior, aiming to identify which pathways are modulated by propolis in S. cerevisiae and potentially influencing cell death. Results We were able to observe 1,339 genes modulated in at least one time point when compared to the reference time (propolis untreated samples) (t-test, p-value 0.01). Enrichment analysis performed by Gene Ontology (GO) Term finder tool showed enrichment for several biological categories among the genes up-regulated in the microarray hybridization such as transport and transmembrane transport and response to stress. Real-time RT-PCR analysis of selected genes showed by our microarray hybridization approach was capable of providing information about S. cerevisiae gene expression modulation with a considerably high level of confidence. Finally, a physical protein-protein (PPPI) network design and global topological analysis stressed the importance of these pathways in response of S. cerevisiae to propolis and were correlated with the transcriptional data obtained thorough the microarray analysis. Conclusions In summary, our data indicate that propolis is largely affecting several pathways in the eukaryotic cell. However, the most prominent pathways are related to oxidative stress, mitochondrial electron transport chain, vacuolar acidification, regulation of macroautophagy associated with protein target to vacuole, cellular response to starvation, and negative regulation of transcription from RNA polymerase II promoter. Our work emphasizes again the importance of S. cerevisiae as a model system to understand at molecular level the mechanism whereby propolis causes cell death in this organism at the concentration herein tested. Our study is the first one that investigates systematically by using functional genomics how propolis influences and modulates the mRNA abundance of an organism and may stimulate further work on the propolis-mediated cell death mechanisms in fungi.

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.

A distinct tymovirus infecting Cassia hoffmannseggii in Brazil

Leaves of Cassia hoffmannseggii, a wild fabaceous species found in the Atlantic Forest, with a severe mosaic symptom were collected in Pernambuco State, Brazil. By transmission electron microscopy, two types of virus particles were found: the first was recognized as particles of a potyvirus, which was later identified as Cowpea aphid-borne mosaic virus; and the second was isometric and present in high concentration. The observation of vesicles at the periphery of chloroplasts suggested a tymovirus infection, which was confirmed by subsequent assays. A serological assay against several tymovirus antisera resulted in positive reaction of this tymo-like virus with an antiserum of Passion fruit yellow mosaic virus. By means of RT-PCR and using degenerated primers for the conserved region of RNA-dependent RNA polymerase (RdRp) gene of tymoviruses, a specific DNA fragment was amplified and sequenced. Based on this sequence, a specific forward primer was synthesized and successfully used to amplify the 3' terminal genome region, containing the partial RdRp gene and the complete coat protein (CP) sequences. The CP was 188 amino acids (aa) long, and the highest CP aa identity was observed with Kennedya yellow mosaic virus (61 %). Based on the current ICTV demarcation criterion, this isolate was considered as a distinct tymovirus and tentatively named as Cassia yellow mosaic-associated virus.

Aberrant transcript produced by a splice donor site deletion in the TECTA gene is associated with autosomal dominant deafness in a Brazilian family

We ascertained a Brazilian family with nine individuals affected by autosomal dominant nonsyndromic sensorineural hearing loss. The bilateral hearing loss affected mainly mid-high frequencies, was apparently stable with an early onset. Microsatellites close to the DFNA8/DFNA12 locus, which harbors the TECTA gene, showed significant multipoint lod scores (32) close to marker D11S4107. Sequencing of the exons and exon-intron boundaries of the TECTA gene in one affected subject revealed the deletion c.5383 + 5delGTGA in the 5' end of intron 16, that includes the last two bases of the donor splice site consensus sequence. This mutation segregates with deafness within the family. To date, 33 different TECTA mutations associated with autossomal dominant hearing loss have been described. Among them is the mutation reported herein, first described by Hildebrand et al. (2011) in a UK family. The audioprofiles from the UK and Brazilian families were similar. In order to investigate the transcripts produced by the mutated allele, we performed cDNA analysis of a lymphoblastoid cell line from an affected heterozygote with the c.5383 + 5delGTGA and a noncarrier from the same family. The analysis allowed us to identify an aberrant transcript with skipping of exon 16, without affecting the reading frame. One of the dominant TECTA mutations already described, a synonymous substitution in exon 16 (c.5331 G<A), was also shown to affect splicing resulting in an aberrant transcript lacking exon 16. Despite the difference in the DNA level, both the synonymous substitution in exon 16 (c.5331 G<A) and the mutation described herein affect splicing of exon 16, leading to its skipping. At the protein level they would have the same effect, an in-frame deletion of 37 amino-acids (p.S1758Y/G1759_N1795del) probably leading to an impaired function of the ZP domain. Thus, like the TECTA missense mutations associated with dominant hearing loss, the c5383 + 5delGTGA mutation does not have an inactivating effect on the protein. (C) 2012 Elsevier B.V. All rights reserved.

Mechanical loading induces the expression of a Pol I regulon at the onset of skeletal muscle hypertrophy

von Walden F, Casagrande V, Ostlund Farrants AK, Nader GA. Mechanical loading induces the expression of a Pol I regulon at the onset of skeletal muscle hypertrophy. Am J Physiol Cell Physiol 302: C1523-C1530, 2012. First published March 7, 2012; doi:10.1152/ajpcell.00460.2011.-The main goal of the present study was to investigate the regulation of ribosomal DNA (rDNA) gene transcription at the onset of skeletal muscle hypertrophy. Mice were subjected to functional overload of the plantaris by bilateral removal of the synergist muscles. Mechanical loading resulted in muscle hypertrophy with an increase in rRNA content. rDNA transcription, as determined by 45S pre-rRNA abundance, paralleled the increase in rRNA content and was consistent with the onset of the hypertrophic response. Increased transcription and protein expression of c-Myc and its downstream polymerase I (Pol I) regulon (POL1RB, TIF-1A, PAF53, TTF1, TAF1C) was also consistent with the increase in rRNA. Similarly, factors involved in rDNA transcription, such as the upstream binding factor and the Williams syndrome transcription factor, were induced by mechanical loading in a corresponding temporal fashion. Chromatin immunoprecipitation revealed that these factors, together with Pol I, were enriched at the rDNA promoter. This, in addition to an increase in histone H3 lysine 9 acetylation, demonstrates that mechanical loading regulates rRNA synthesis by inducing a gene expression program consisting of a Pol I regulon, together with accessory factors involved in transcription and chromatin remodeling at the rDNA promoter. Altogether, these data indicate that transcriptional and epigenetic mechanisms take place in the regulation of ribosome production at the onset of muscle hypertrophy.

Extracytoplasmic function (ECF) sigma factor sigma(F) is involved in Caulobacter crescentus response to heavy metal stress

Background: The alpha-proteobacterium Caulobacter crescentus inhabits low-nutrient environments and can tolerate certain levels of heavy metals in these sites. It has been reported that C. crescentus responds to exposure to various heavy metals by altering the expression of a large number of genes. Results: In this work, we show that the ECF sigma factor sigma(F) is one of the regulatory proteins involved in the control of the transcriptional response to chromium and cadmium. Microarray experiments indicate that sigma(F) controls eight genes during chromium stress, most of which were previously described as induced by heavy metals. Surprisingly, sigma(F) itself is not strongly auto-regulated under metal stress conditions. Interestingly, sigma(F)-dependent genes are not induced in the presence of agents that generate reactive oxygen species. Promoter analyses revealed that a conserved sigma(F)-dependent sequence is located upstream of all genes of the sigma(F) regulon. In addition, we show that the second gene in the sigF operon acts as a negative regulator of sigma(F) function, and the encoded protein has been named NrsF (Negative regulator of sigma F). Substitution of two conserved cysteine residues (C131 and C181) in NrsF affects its ability to maintain the expression of sigma(F)-dependent genes at basal levels. Furthermore, we show that sigma(F) is released into the cytoplasm during chromium stress and in cells carrying point mutations in both conserved cysteines of the protein NrsF. Conclusion: A possible mechanism for induction of the sigma(F)-dependent genes by chromium and cadmium is the inactivation of the putative anti-sigma factor NrsF, leading to the release of sigma(F) to bind RNA polymerase core and drive transcription of its regulon.

Extracytoplasmic function (ECF) sigma factor σF is involved in Caulobacter crescentus response to heavy metal stress

Background The α-proteobacterium Caulobacter crescentus inhabits low-nutrient environments and can tolerate certain levels of heavy metals in these sites. It has been reported that C. crescentus responds to exposure to various heavy metals by altering the expression of a large number of genes. Results In this work, we show that the ECF sigma factor σF is one of the regulatory proteins involved in the control of the transcriptional response to chromium and cadmium. Microarray experiments indicate that σF controls eight genes during chromium stress, most of which were previously described as induced by heavy metals. Surprisingly, σF itself is not strongly auto-regulated under metal stress conditions. Interestingly, σF-dependent genes are not induced in the presence of agents that generate reactive oxygen species. Promoter analyses revealed that a conserved σF-dependent sequence is located upstream of all genes of the σF regulon. In addition, we show that the second gene in the sigF operon acts as a negative regulator of σF function, and the encoded protein has been named NrsF (Negative regulator of sigma F). Substitution of two conserved cysteine residues (C131 and C181) in NrsF affects its ability to maintain the expression of σF-dependent genes at basal levels. Furthermore, we show that σF is released into the cytoplasm during chromium stress and in cells carrying point mutations in both conserved cysteines of the protein NrsF. Conclusion A possible mechanism for induction of the σF-dependent genes by chromium and cadmium is the inactivation of the putative anti-sigma factor NrsF, leading to the release of σF to bind RNA polymerase core and drive transcription of its regulon.

Identification of SL addition trans-splicing acceptor sites in the internal transcribed spacer I region of pre-rRNA in Leishmania (Leishmania) amazonensis

Trypanosomatidae is a family of early branching eukaryotes harbouring a distinctive repertoire of gene expression strategies. Functional mature messenger RNA is generated via the trans-splicing and polyadenylation processing of constitutively transcribed polycistronic units. Recently, trans-splicing of pre-small subunit ribosomal RNA in the 5' external transcribed spacer region and of precursor tRNAsec have been described. Here, we used a previously validated semi-nested reverse transcription-polymerase chain reaction strategy to investigate internal transcribed spacer (ITS) I acceptor sites in total RNA from Leishmania (Leishmania) amazonensis. Two distinct spliced leader-containing RNAs were detected indicating that trans-splicing reactions occur at two AG acceptor sites mapped in this ITS region. These data provide further evidence of the wide spectrum of RNA molecules that act as trans-splicing acceptors in trypanosomatids.

Human immunodeficiency virus prevalence, incidence, and residual risk of transmission by transfusions at Retrovirus Epidemiology Donor Study-II blood centers in Brazil

BACKGROUND: In Brazil nationally representative donor data are limited on human immunodeficiency virus (HIV) prevalence, incidence, and residual transfusion risk. The objective of this study was to analyze HIV data obtained over 24 months by the Retrovirus Epidemiology Donor Study-II program in Brazil. STUDY DESIGN AND METHODS: Donations reactive to third-and fourth-generation immunoassays (IAs) were further confirmed by a less-sensitive (LS) IA algorithm and Western blot (WB). Incidence was calculated for first-time (FT) donors using the LS-EIA results and for repeat donors with a model developed to include all donors with a previous negative donation. Residual risk was projected by multiplying composite FT and repeat donor incidence rates by HIV marker-negative infectious window periods. RESULTS: HIV prevalence among FT donors was 92.2/ 105 donations. FT and repeat donor and composite incidences were 38.5 (95% confidence interval [CI], 25.651.4), 22.5 (95% CI, 17.6-28.0), and 27.5 (95% CI, 22.0-33.0) per 100,000 person-years, respectively. Male and community donors had higher prevalence and incidence rates than female and replacement donors. The estimated residual risk of HIV transfusion transmission was 11.3 per 106 donations (95% CI, 8.4-14.2), which could be reduced to 4.2 per 106 donations (95% CI, 3.2-5.2) by use of individual-donation nucleic acid testing (NAT). CONCLUSION: The incidence and residual transfusion risk of HIV infection are relatively high in Brazil. Implementation of NAT will not be sufficient to decrease transmission rates to levels seen in the United States or Europe; therefore, other measures focused on decreasing donations by at-risk individuals are also necessary.

Correlates of human immunodeficiency virus cervicovaginal shedding among postmenopausal and fertile-aged women

Objective: The aims of this study were to compare the intensity of human immunodeficiency virus (HIV)-RNA genital shedding among postmenopausal (PM) and fertile-aged (F) women and to investigate the association between viral shedding and gynecological features, HIV plasma viral loads, and other markers of HIV disease progression. Methods: We interviewed 146 HIV-infected women (73 PM/73 F) in search of gynecological complaints and signs and symptoms of HIV disease and obtained additional information concerning HIV infection by medical chart review. Cervicovaginal lavages (CVLs) were collected for assessment of HIV shedding. Laboratory analyses included CD4(+) cell counts, HIV-RNA quantitation in plasma and CVL, and screening for concurrent genital infections. Results: HIV-RNA genital shedding was detected in 16.4% of PM and 21.9% of F women (P = 0.400), and the intensity of HIV shedding did not differ between both groups (means-PM: 1.4log/mL; F: 1.4log/mL; P = 0.587). Three women (2 PM/1 F) exhibited viral shedding in the absence of detectable viremia. HIV plasma viral loads correlated with HIV shedding in both groups. In multivariable analysis, HIV plasma viral loads were independently associated with HIV shedding in both groups. Moreover, the intensity of shedding was independently associated with vaginal pH, tumor necrosis factor a concentrations in CVL, and HIV plasma viral loads. Conclusions: Despite significant changes that occur in the vaginal mucosa of PM women, HIV cervicovaginal shedding was not significantly influenced by this state in our cohort. In contrast, increased vaginal pH and genital inflammation, evidenced by increased tumor necrosis factor alpha concentrations in CVL and HIV plasma viral loads, were independently associated with the intensity of HIV shedding in PM and F women.