Evidence for common short natural trans sense-antisense pairing between transcripts from protein coding genes

Truncating TINF2 p.Tyr312Ter variant and inherited breast cancer susceptibility

TINF2 is a critical subunit of the shelterin complex, which protects and maintains the length of telomeres. Pathogenic missense and truncating TINF2 mutations are causative for dyskeratosis congenita (DC), a rare, dominantly inherited bone marrow failure syndrome characterized by mucocutaneous abnormalities and cancer predisposition. Recent reports indicate that specific TINF2 truncating mutations act as high penetrance cancer predisposition alleles outside DC context, including breast cancer in their tumor spectrum. Here, we have evaluated the role of germline mutations in TINF2 and other shelterin genes in inherited breast cancer susceptibility using exome sequencing data from 98 Northern Finnish breast cancer cases with indication of inherited disease predisposition as a discovery cohort. A single protein truncating variant, TINF2 p.Tyr312Ter, was identified in one of the cases (1/98), and four more carriers were observed in the subsequently genotyped unselected breast cancer cohort (4/1904). None of the carriers were reported to have DC. TINF2 p.Tyr312Ter resulted in stable short form of mRNA transcript, and normal telomere length has been indicated by a recent report. Although recurrent in cases (total of 5/2095), TINF2 p.Tyr312Ter is also present in Finnish population controls (8/12,517), and the observed 4-fold higher frequency in cases falls at most into the range of moderate breast cancer risk alleles (OR 3.74, 95% CI 1.22-11.45, p = 0.029). Current results indicate that not all TINF2 truncating variants are high cancer risk alleles and add further evidence that different TINF2 mutations can have very diverse effects on the disease phenotype.

ASSA: Fast identification of statistically significant interactions between long RNAs

The discovery of thousands of long noncoding RNAs (lncRNAs) in mammals raises a question about their functionality. It has been shown that some of them are involved in post-transcriptional regulation of other RNAs and form inter-molecular duplexes with their targets. Sequence alignment tools have been used for transcriptome-wide prediction of RNA-RNA interactions. However, such approaches have poor prediction accuracy since they ignore RNA's secondary structure. Application of the thermodynamics-based algorithms to long transcripts is not computationally feasible on a large scale. Here, we describe a new computational pipeline ASSA that combines sequence alignment and thermodynamics-based tools for efficient prediction of RNA-RNA interactions between long transcripts. To measure the hybridization strength, the sum energy of all the putative duplexes is computed. The main novelty implemented in ASSA is the ability to quickly estimate the statistical significance of the observed interaction energies. Most of the functional hybridizations between long RNAs were classified as statistically significant. ASSA outperformed 11 other tools in terms of the Area Under the Curve on two out of four test sets. Additionally, our results emphasized a unique property of the [Formula: see text] repeats with respect to the RNA-RNA interactions in the human transcriptome. ASSA is available at https://sourceforge.net/projects/assa/.

Non-coding Transcripts from Enhancers: New Insights into Enhancer Activity and Gene Expression Regulation

Long non-coding RNAs (lncRNAs) have gained widespread interest in the past decade owing to their enormous amount and surprising functions implicated in a variety of biological processes. Some lncRNAs exert function as enhancers, i.e., activating gene transcription by serving as the cis-regulatory molecules. Furthermore, recent studies have demonstrated that many enhancer elements can be transcribed and produce RNA molecules, which are termed as enhancer RNAs (eRNAs). The eRNAs are not merely the by-product of the enhancer transcription. In fact, many of them directly exert or regulate enhancer activity in gene activation through diverse mechanisms. Here, we provide an overview of enhancer activity, transcription of enhancer itself, characteristics of eRNAs, as well as their roles in regulating enhancer activity and gene expression.

Alu-miRNA interactions modulate transcript isoform diversity in stress response and reveal signatures of positive selection

Primate-specific Alus harbor different regulatory features, including miRNA targets. In this study, we provide evidence for miRNA-mediated modulation of transcript isoform levels during heat-shock response through exaptation of Alu-miRNA sites in mature mRNA. We performed genome-wide expression profiling coupled with functional validation of miRNA target sites within exonized Alus, and analyzed conservation of these targets across primates. We observed that two miRNAs (miR-15a-3p and miR-302d-3p) elevated in stress response, target RAD1, GTSE1, NR2C1, FKBP9 and UBE2I exclusively within Alu. These genes map onto the p53 regulatory network. Ectopic overexpression of miR-15a-3p downregulates GTSE1 and RAD1 at the protein level and enhances cell survival. This Alu-mediated fine-tuning seems to be unique to humans as evident from the absence of orthologous sites in other primate lineages. We further analyzed signatures of selection on Alu-miRNA targets in the genome, using 1000 Genomes Phase-I data. We found that 198 out of 3177 Alu-exonized genes exhibit signatures of selection within Alu-miRNA sites, with 60 of them containing SNPs supported by multiple evidences (global-F_ST > 0.3, pair-wise-F_ST > 0.5, Fay-Wu’s H < −20, iHS > 2.0, high ΔDAF) and implicated in p53 network. We propose that by affecting multiple genes, Alu-miRNA interactions have the potential to facilitate population-level adaptations in response to environmental challenges.

Aberrant expression of miR-153 is associated with overexpression of hypoxia-inducible factor-1α in refractory epilepsy

Evidence suggest that overexpression of hypoxia-inducible factor-1α (HIF-1α) is linked to multidrug resistance of epilepsy. Here we explored whether aberrant expression of HIF-1α is regulated by miRNAs. Genome-wide microRNA expression profiling was performed on temporal cortex resected from mesial temporal lobe epilepsy (mTLE) patients and age-matched controls. miRNAs that are putative regulator of HIF-1α were predicted via target scan and confirmed by real-time quantitative polymerase chain reaction (RT-qPCR). Mimics or miRNA morpholino inhibitors were transfected in astrocytes and luciferase reporter assay was applied to detect HIF-11α expression. Microarray profiling identified down-regulated miR-153 as a putative regulator of HIF-1α in temporal cortex resected from surgical mTLE patients. RT-qPCR confirmed down-regulation of miR-153 in plasma of mTLE patients in an independent validation cohort. Knockdown of miR-153 significantly enhanced expression of HIF-1α while forced expression of miR-153 dramatically inhibited HIF-1α expression in pharmacoresistant astrocyte model. Luciferase assay established that miR-153 might inhibit HIF-1α expression via directly targeting two binding sites in the 3′UTR region of HIF-1α transcript. These data suggest that down-regulation of miR-153 may contribute to enhanced expression of HIF-1α in mTLE and serve as a novel biomarker and treatment target for epilepsy.

Long noncoding RNA as a regulator for transcription

Investigation of noncoding RNAs is in rapid progress, especially regarding translational repression by small (short) noncoding RNAs like microRNAs with 20-25 nucleotide-lengths, while long noncoding RNAs with nucleotide length of more than two hundred are also emerging. Indeed, our analysis has revealed that a long noncoding RNA transcribed from cyclin D1 promoter of 200 and 300 nucleotides exerts transcriptional repression through its binding protein TLS instead of translational repression. Translational repression is executed by short noncoding RNAs, while transcriptional repression is mainly done by long noncoding RNAs. These long noncoding RNAs are heterogeneous molecules and employ divergent molecular mechanisms to exert transcriptional repression. In this review, I overview recent publications regarding the transcription regulation by long noncoding RNAs and explore their biological significance. In addition, the relation between a random transcriptional activity of RNA polymerase II and the origin of long noncoding RNAs is discussed.

Trm9-Catalyzed tRNA Modifications Regulate Global Protein Expression by Codon-Biased Translation

Post-transcriptional modifications of transfer RNAs (tRNAs) have long been recognized to play crucial roles in regulating the rate and fidelity of translation. However, the extent to which they determine global protein production remains poorly understood. Here we use quantitative proteomics to show a direct link between wobble uridine 5-methoxycarbonylmethyl (mcm⁵) and 5-methoxy-carbonyl-methyl-2-thio (mcm⁵s²) modifications catalyzed by tRNA methyltransferase 9 (Trm9) in tRNA^Arg(UCU) and tRNA^Glu(UUC) and selective translation of proteins from genes enriched with their cognate codons. Controlling for bias in protein expression and alternations in mRNA expression, we find that loss of Trm9 selectively impairs expression of proteins from genes enriched with AGA and GAA codons under both normal and stress conditions. Moreover, we show that AGA and GAA codons occur with high frequency in clusters along the transcripts, which may play a role in modulating translation. Consistent with these results, proteins subject to enhanced ribosome pausing in yeast lacking mcm⁵U and mcm⁵s²U are more likely to be down-regulated and contain a larger number of AGA/GAA clusters. Together, these results suggest that Trm9-catalyzed tRNA modifications play a significant role in regulating protein expression within the cell.

Here we present evidence for a more complicated role for transfer RNAs (tRNAs) than as mere adapters that link the genetic code in messenger RNA (mRNA) to the amino acid sequence of a protein during translation. tRNAs have long been known to be modified with dozens of different chemical structures other than the 4 canonical ribonucleosides, though the role of these modifications in controlling translation is poorly understood. By quantifying the expression of thousands of proteins in the yeast S. cerevisiae, we identified a mechanistic link between modified ribonucleosides located at the wobble position of two tRNAs, tRNA^Arg(UCU) and tRNA^Glu(UUC), and the translation of proteins derived from genes enriched with codons read by these tRNAs: AGA and GAA. In cells lacking the enzyme that inserts these modifications, tRNA methyltransferase 9 (Trm9), we found a significant reduction in proteins from genes enriched with AGA and GAA codons and with runs of these codons. Also, mRNAs enriched with runs of AGA and GAA codons are subject to stalled translation on ribosomes in yeast lacking mcm⁵U and mcm⁵s²U. Together, these results reveal a distinct role for Trm9-catalyzed tRNA modifications in selectively regulating the expression of proteins enriched with AGA and GAA codons.

Macrophage infiltration promotes invasiveness of breast cancer cells via activating long non-coding RNA UCA1

There is now considerable evidence supporting the view that macrophage infiltration is playing a critical role in the proliferation and progression of breast cancer but the underlying molecular mechanisms remain largely unknown. To this end, using long non-coding RNA (lncRNA) expression profiling, we examined changes in lncRNA expression in breast cancer cells treated with conditioned medium (CM) from cultured human THP-1 macrophages. We found that treatment with macrophage CM induced the expression of numerous lncRNAs, including urothelial cancer associated 1 (UCA1). Knockdown of UCA1 using shRNA inhibited AKT phosphorylation and abolished invasiveness of tumor cells induced by macrophage CM. Consistent with these results; we further showed that UCA1 level was significantly enhanced in human primary breast tumors and correlated with advanced clinical stage, supporting its role in promoting carcinogenesis and progression of breast cancer. Together, these results suggest that macrophage could promote invasiveness of breast cancer cells by enhancing expression of lncRNA UCA1.

Association between plasma homocysteine and progression of early nephropathy in type 2 diabetic patients

There is now growing evidence supporting the association between renal insufficiency and accumulation of plasma homocysteine (Hcy). However, the role of Hcy in the development of diabetic nephropathy (DN) in type 2 diabetic patients is not clearly elucidated. To this end, we performed a prospective observational study in 208 patients and 49 controls. We show that baseline level of Hcy is significantly enhanced in patients with DN and is associated with the severity of the disease. Focusing on patients at early DN stage (n = 157), after four-year follow-up, we find that increase in plasma Hcy level correlates with greater renal failure characterized by faster decline in estimated glomerular filtration rate (eGFR). Using a multivariate linear regression model, we show that plasma Hcy remains significantly associated with eGFR decline after controlling for other progression promoters. Our results support that plasma Hcy is an independent risk factor as well as an early predictor for DN progression in type 2 diabetic patients.

The association between MMP-12 82 A/G polymorphism and susceptibility to various malignant tumors: a meta-analysis

Matrix metalloproteinases (MMPs) are a family of zinc-dependent endopeptidases responsible for degrading essentially all components of the extracellular matrix (ECM). Accumulating evidence suggests that MMPs might play a critical role in growth, invasion, and metastasis of malignant tumors. A single nucleotide polymorphism (SNP) in the promoter region of MMP-12, MMP-12 82 A/G (rs2276109), has been recognized to play a critical role in regulating the expression of MMP-12, however, its correlation with tumor susceptibility remains controversial. To address this issue, we performed meta-analysis to investigate the association MMP-12 82 A/G polymorphism and susceptibility of nine malignant tumors from 11 studies, including 6153 cancer patients and 6838 controls. Two reviewers independently screened studies for eligibility and extracted data for included studies. While overall no evident association between MMP-12 82 A/G and tumor susceptibility was observed, subgroup analysis revealed a specific role of G allele in increasing the susceptibility for epithelial ovarian carcinoma (EOC) using the allele model (fixed effects OR = 2.45, 95% CI = 1.46-4.10, P = 0.001) and the dominant model (fixed effects OR = 2.52, 95% CI = 1.49-4.24, P = 0.001). We thus suggest that G allele of MMP-12 82 A/G polymorphism is a genetic risk factor for EOC.

Integrative analysis reveals enhanced regulatory effects of human long intergenic non-coding RNAs in lung adenocarcinoma

Although there is an accumulating appreciation of the key roles that long intergenic non-coding RNAs (lincRNAs) play in diverse cellular processes, our knowledge of how lincRNAs function in cancer remains sparse. Here, we present a comprehensive landscape of RNA-seq transcriptome profiles of lung adenocarcinomas and their paired normal counterparts to unravel gene regulation rules of lincRNAs. Consistent with previous findings of co-expression between neighboring protein-coding genes, lincRNAs were typically co-expressed with their neighboring genes, which was found in both cancerous and normal tissues. By building a mathematical model based on correlated gene expression, we distinguished an additional subset of lincRNAs termed "regulatory lincRNAs", representing their dominant roles in gene regulation. The number of regulatory lincRNAs was significantly higher in cancerous compared to normal tissues, and most of them positively regulated protein-coding genes in trans. Functional validation, using knockdown, determined that regulatory lincRNA, GAS5, affected its predicted protein-coding targets. Moreover, we discovered hundreds of differentially expressed regulatory lincRNAs with inclusion of some cancer-associated lincRNAs. Our integrated analysis reveals enhanced regulatory effects of lincRNAs and provides a resource for the study of regulatory lincRNAs that play critical roles in lung adenocarcinoma.

Enhanced expression of stem cell markers and drug resistance in sphere-forming non-small cell lung cancer cells

There is growing evidence suggesting that cancer stem cells (CSCs) are playing critical roles in tumor progression, metastasis and drug resistance. However, the role of CSCs in non-small cell lung cancer (NSCLC) remains elusive. In this study, we enriched for stem-like cells from tumor spheres derived from NSCLC cell line A549 cultured in serum-free medium. Our results showed that sphere-derived cells expressed various stem cell markers such as CD44, CD133, Sox2 and Oct4. Compared with the corresponding cells in monolayer cultures, sphere-derived cells showed marked morphologic changes and increased expression of the stem cell markers CD133. Furthermore, we found that sphere-derived cells exhibited increased proliferation, cell-cycle progression as well as drug-resistant properties as compared to A549 adherent cells. Consistently, expression of several drug resistance proteins, including lung resistance-related protein (LRP), glutathion-S-transferase-π (GST-π) and multidrug resistance proteins-1 (MRP1) were all significantly enhanced in sphere-derived cells. These results indicate the enrichment of CSCs in sphere cultures and support their role in regulating drug resistance in NSCLC.

TRIM29 functions as an oncogene in gastric cancer and is regulated by miR-185

Tripartite motif-containing 29 (TRIM29) belongs to TRIM family of transcription factors and may function as an oncogene or a tumor suppressor depending on the tumor types. Overexpression of TRIM29 is frequently observed in gastric cancer but the underlying mechanisms remain largely unknown. In the present study, we investigated the function of TRIM29 in gastric cancer-derived cell line MGC803. RNAi-mediated silencing of TRIM29 resulted in significantly reduced cell proliferation and colony formation, as well as G1-S cell cycle arrest and apoptosis. Interestingly, expression levels of β-catenin, cyclin D1 and c-Myc were all downregulated in TRIM29 knockdown cells, indicating that TRIM29 is involved in regulating the activity of Wnt/β-catenin signaling pathway. Furthermore, based on target prediction and luciferase assay, we identified TRIM29 as a potential target of miR-185, which is frequently downregulated in gastric cancer. Over-expression of miR-185 in MGC803 cells inhibited TRIM29 expression and activity of Wnt/β-catenin signaling. Taken together, our results suggest that TRIM29 functions as an oncogene in gastric cancer and is regulated by miR-185.

Regulatory Roles for Long ncRNA and mRNA

Recent advances in high-throughput sequencing technology have identified the transcription of a much larger portion of the genome than previously anticipated. Especially in the context of cancer it has become clear that aberrant transcription of both protein-coding and long non-coding RNAs (lncRNAs) are frequent events. The current dogma of RNA function describes mRNA to be responsible for the synthesis of proteins, whereas non-coding RNA can have regulatory or epigenetic functions. However, this distinction between protein coding and regulatory ability of transcripts may not be that strict. Here, we review the increasing body of evidence for the existence of multifunctional RNAs that have both protein-coding and trans-regulatory roles. Moreover, we demonstrate that coding transcripts bind to components of the Polycomb Repressor Complex 2 (PRC2) with similar affinities as non-coding transcripts, revealing potential epigenetic regulation by mRNAs. We hypothesize that studies on the regulatory ability of disease-associated mRNAs will form an important new field of research.

Heat shock factor binding in Alu repeats expands its involvement in stress through an antisense mechanism

Background

Alu RNAs are present at elevated levels in stress conditions and, consequently, Alu repeats are increasingly being associated with the physiological stress response. Alu repeats are known to harbor transcription factor binding sites that modulate RNA pol II transcription and Alu RNAs act as transcriptional co-repressors through pol II binding in the promoter regions of heat shock responsive genes. An observation of a putative heat shock factor (HSF) binding site in Alu led us to explore whether, through HSF binding, these elements could further contribute to the heat shock response repertoire.

Results

Alu density was significantly enriched in transcripts that are down-regulated following heat shock recovery in HeLa cells. ChIP analysis confirmed HSF binding to a consensus motif exhibiting positional conservation across various Alu subfamilies, and reporter constructs demonstrated a sequence-specific two-fold induction of these sites in response to heat shock. These motifs were over-represented in the genic regions of down-regulated transcripts in antisense oriented Alus. Affymetrix Exon arrays detected antisense signals in a significant fraction of the down-regulated transcripts, 50% of which harbored HSF sites within 5 kb. siRNA knockdown of the selected antisense transcripts led to the over-expression, following heat shock, of their corresponding down-regulated transcripts. The antisense transcripts were significantly enriched in processes related to RNA pol III transcription and the TFIIIC complex.

Conclusions

We demonstrate a non-random presence of Alu repeats harboring HSF sites in heat shock responsive transcripts. This presence underlies an antisense-mediated mechanism that represents a novel component of Alu and HSF involvement in the heat shock response.

Locating protein-coding sequences under selection for additional, overlapping functions in 29 mammalian genomes

The degeneracy of the genetic code allows protein-coding DNA and RNA sequences to simultaneously encode additional, overlapping functional elements. A sequence in which both protein-coding and additional overlapping functions have evolved under purifying selection should show increased evolutionary conservation compared to typical protein-coding genes--especially at synonymous sites. In this study, we use genome alignments of 29 placental mammals to systematically locate short regions within human ORFs that show conspicuously low estimated rates of synonymous substitution across these species. The 29-species alignment provides statistical power to locate more than 10,000 such regions with resolution down to nine-codon windows, which are found within more than a quarter of all human protein-coding genes and contain ∼2% of their synonymous sites. We collect numerous lines of evidence that the observed synonymous constraint in these regions reflects selection on overlapping functional elements including splicing regulatory elements, dual-coding genes, RNA secondary structures, microRNA target sites, and developmental enhancers. Our results show that overlapping functional elements are common in mammalian genes, despite the vast genomic landscape.

The long arm of long noncoding RNAs: roles as sensors regulating gene transcriptional programs

A major surprise arising from genome-wide analyses has been the observation that the majority of the genome is transcribed, generating noncoding RNAs (ncRNAs). It is still an open question whether some or all of these ncRNAs constitute functional networks regulating gene transcriptional programs. However, in light of recent discoveries and given the diversity and flexibility of long ncRNAs and their abilities to nucleate molecular complexes and to form spatially compact arrays of complexes, it becomes likely that many or most ncRNAs act as sensors and integrators of a wide variety of regulated transcriptional responses and probably epigenetic events. Because many RNA-binding proteins, on binding RNAs, show distinct allosteric conformational alterations, we suggest that a ncRNA/RNA-binding protein-based strategy, perhaps in concert with several other mechanistic strategies, serves to integrate transcriptional, as well as RNA processing, regulatory programs.

Promoter-associated long noncoding RNAs repress transcription through a RNA binding protein TLS

The majority of the human genome is found to be transcribed and generates mostly noncoding (nc) RNAs that do not possess protein information. MicroRNAs are one of the well-identified small ncRNAs, but occupy merely a fraction of ncRNAs. Long (large) ncRNAs are emerging as a novel class of ncRNAs, but knowledge of these ncRNAs is far less accumulated. Long ncRNAs are tentatively classified as an ncRNA species containing more than 200 nucleotides. Recently, a long promoter-associated ncRNA (pncRNA) has been identified to be transcribed from the cyclin D1 promoter upon induction by genotoxic factors like ionizing-irradiation. The cyclin D1 pncRNA is specifically bound with an RNA-binding protein TLS (Translocated in liposarcoma) and exerts transcriptional repression through histone acetyltransferase (HAT) inhibitory activity. Analysis of TLS and the pncRNAs could provide a model for elucidating their roles inregulation of mammalian transcriptional programs. The pncRNA binding to TLS turns out to be an essential event for the HAT inhibitory activity. A key consensus sequence of the pncRNA is composed of GGUG, while not every RNA sequence bearing GGUG is targeted by TLS, suggesting that a secondary structure of the GGUG-bearing RNAs is also involved in recognition by TLS. Taken together, TLS is a unique mediator between signals of the long ncRNAs and transcription, suggesting that RNA networking functions in living cells.(1-3).

Whole transcriptome analysis: what are we still missing?

New technologies such as tag-based sequencing and tiling arrays have provided unique insights into the transcriptional output of cells. Many new RNA classes have been uncovered in the past decade, despite limitations in current technologies. Even as the repertoire of known functional elements of the transcriptome increases and contemporary technologies become mainstream, inadequacies in conventional protocols for library preparation, sequencing and mapping continue to hamper revelation of the entire transcriptome of cells. In this article, we review current protocols and outline their deficiencies. We also provide our view on what we may be overlooking in the transcriptome, despite exhaustive investigations, and indicate future areas of technological development and research.