Highly accessible AU-rich regions in 3' untranslated regions are hotspots for binding of regulatory factors

Revealing miRNAs patterns by employing matrix representations and energy analysis

MicroRNAs (miRNAs) are small, non-coding RNA molecules that regulate gene expression. Despite their relatively short length (about 21 nucleotides), they can regulate thousands of transcripts within a cell. Due to their low complementarity to targets, studying their activity and binding region preferences (3'UTR, 5'UTR, or CDS) is challenging. In this paper, we analyzed a set of human miRNAs to uncover their general patterns. We began with a sequence logo to verify conservation at specific positions. To discover long-range correlations, we employed chaos game representation (CGR) and genomatrix, methods that enable both graphical and analytical analysis of sequence sets and are well-established in bioinformatics. Our results showed that miRNAs exhibit strongly non-random and characteristic patterns. To incorporate physicochemical properties into the analysis, we applied the electron-ion interaction potential (EIIP) parameter. An important part of our study was to validate the division of miRNAs into two parts-seed and puzzle. The seed region is responsible for target binding, while the puzzle region likely interacts with the RISC complex. We estimated duplex binding energy within the 3'UTR, 5'UTR, and CDS regions using the miRanda tool. Based on the median energy distribution, we divided the miRNAs into two subsets, reflecting different patterns in chaos game representation. Interestingly, one subset displayed significant similarity to conserved and highly confidential miRNAs. Our results confirm the low complementarity of miRNA/mRNA interactions and support the functional division of miRNA structure. Additionally, we present findings related to the localization of transcript target sites, which form the basis for further analyses.

pyPAGE: A framework for Addressing biases in gene-set enrichment analysis-A case study on Alzheimer's disease

Inferring the driving regulatory programs from comparative analysis of gene expression data is a cornerstone of systems biology. Many computational frameworks were developed to address this problem, including our iPAGE (information-theoretic Pathway Analysis of Gene Expression) toolset that uses information theory to detect non-random patterns of expression associated with given pathways or regulons. Our recent observations, however, indicate that existing approaches are susceptible to the technical biases that are inherent to most real world annotations. To address this, we have extended our information-theoretic framework to account for specific biases and artifacts in biological networks using the concept of conditional information. To showcase pyPAGE, we performed a comprehensive analysis of regulatory perturbations that underlie the molecular etiology of Alzheimer's disease (AD). pyPAGE successfully recapitulated several known AD-associated gene expression programs. We also discovered several additional regulons whose differential activity is significantly associated with AD. We further explored how these regulators relate to pathological processes in AD through cell-type specific analysis of single cell and spatial gene expression datasets. Our findings showcase the utility of pyPAGE as a precise and reliable biomarker discovery in complex diseases such as Alzheimer's disease.

Heterogenous nuclear ribonucleoprotein D-like controls endothelial cell functions

HnRNPs are ubiquitously expressed RNA-binding proteins, tightly controlling posttranscriptional gene regulation. Consequently, hnRNP networks are essential for cellular homeostasis and their dysregulation is associated with cancer and other diseases. However, the physiological function of hnRNPs in non-cancerous cell systems are poorly understood. We analyzed the importance of HNRNPDL in endothelial cell functions. Knockdown of HNRNPDL led to impaired proliferation, migration and sprouting of spheroids. Transcriptome analysis identified cyclin D1 (CCND1) and tropomyosin 4 (TPM4) as targets of HNRNPDL, reflecting the phenotypic changes after knockdown. Our findings underline the importance of HNRNPDL for the homeostasis of physiological processes in endothelial cells.

Landscape of semi-extractable RNAs across five human cell lines

Phase-separated membraneless organelles often contain RNAs that exhibit unusual semi-extractability using the conventional RNA extraction method, and can be efficiently retrieved by needle shearing or heating during RNA extraction. Semi-extractable RNAs are promising resources for understanding RNA-centric phase separation. However, limited assessments have been performed to systematically identify and characterize semi-extractable RNAs. In this study, 1074 semi-extractable RNAs, including ASAP1, DANT2, EXT1, FTX, IGF1R, LIMS1, NEAT1, PHF21A, PVT1, SCMH1, STRG.3024.1, TBL1X, TCF7L2, TVP23C-CDRT4, UBE2E2, ZCCHC7, ZFAND3 and ZSWIM6, which exhibited consistent semi-extractability were identified across five human cell lines. By integrating publicly available datasets, we found that semi-extractable RNAs tend to be distributed in the nuclear compartments but are dissociated from the chromatin. Long and repeat-containing semi-extractable RNAs act as hubs to provide global RNA-RNA interactions. Semi-extractable RNAs were divided into four groups based on their k-mer content. The NEAT1 group preferred to interact with paraspeckle proteins, such as FUS and NONO, implying that RNAs in this group are potential candidates of architectural RNAs that constitute nuclear bodies.

The emerging roles of CFIm25 (NUDT21/CPSF5) in human biology and disease

The mammalian cleavage factor I subunit CFIm25 (NUDT21) binds to the UGUA sequences of precursor RNAs. Traditionally, CFIm25 is known to facilitate 3' end formation of pre-mRNAs resulting in the formation of polyadenylated transcripts. Recent studies suggest that CFIm25 may be involved in the cyclization and hence generation of circular RNAs (circRNAs) that contain UGUA motifs. These circRNAs act as competing endogenous RNAs (ceRNAs) that disrupt the ceRNA-miRNA-mRNA axis. Other emerging roles of CFIm25 include regulating both alternative splicing and alternative polyadenylation (APA). APA generates different sized transcripts that may code for different proteins, or more commonly transcripts that code for the same protein but differ in the length and sequence content of their 3' UTRs (3' UTR-APA). CFIm25 mediated global changes in 3' UTR-APA affect human physiology including spermatogenesis and the determination of cell fate. Deregulation of CFIm25 and changes in 3' UTR-APA have been implicated in several human diseases including cancer. In many cancers, CFIm25 acts as a tumor suppressor. However, there are some cancers where CFIm25 has the opposite effect. Alterations in CFIm25-driven 3' UTR-APA may also play a role in neural dysfunction and fibrosis. CFIm25 mediated 3' UTR-APA changes can be used to generate specific signatures that can be used as potential biomarkers in development and disease. Due to the emerging role of CFIm25 as a regulator of the aforementioned RNA processing events, modulation of CFIm25 levels may be a novel viable therapeutic approach. This article is categorized under: RNA Processing > 3' End Processing RNA in Disease and Development > RNA in Disease.

MAPKAPK2-centric transcriptome profiling reveals its major role in governing molecular crosstalk of IGFBP2, MUC4, and PRKAR2B during HNSCC pathogenesis

Transcriptome analysis of head and neck squamous cell carcinoma (HNSCC) has been pivotal to comprehending the convoluted biology of HNSCC tumors. MAPKAPK2 or MK2 is a critical modulator of the mRNA turnover of crucial genes involved in HNSCC progression. However, MK2-centric transcriptome profiles of tumors are not well known. This study delves into HNSCC progression with MK2 at the nexus to delineate the biological relevance and intricate crosstalk of MK2 in the tumor milieu. We performed next-generation sequencing-based transcriptome profiling of HNSCC cells and xenograft tumors to ascertain mRNA expression profiles in MK2-wild type and MK2-knockdown conditions. The findings were validated using gene expression assays, immunohistochemistry, and transcript turnover studies. Here, we identified a pool of crucial MK2-regulated candidate genes by annotation and differential gene expression analyses. Regulatory network and pathway enrichment revealed their significance and involvement in the HNSCC pathogenesis. Additionally, 3'-UTR-based filtering recognized important MK2-regulated downstream target genes and validated them by nCounter gene expression assays. Finally, immunohistochemistry and transcript stability studies revealed the putative role of MK2 in regulating the transcript turnover of IGFBP2, MUC4, and PRKAR2B in HNSCC. Conclusively, MK2-regulated candidate genes were identified in this study, and their plausible involvement in HNSCC pathogenesis was elucidated. These genes possess investigative values as targets for diagnosis and therapeutic interventions for HNSCC.

Detection of features predictive of microRNA targets by integration of network data

Gene activity is controlled by multiple molecular mechanisms, for instance through transcription factors or by microRNAs (miRNAs), among others. Established bioinformatics tools for the prediction of miRNA target genes face the challenge of ensuring accuracy, due to high false positive rates. Further, these tools present poor overlap. However, we demonstrated that it is possible to filter good predictions of miRNA targets from the bulk of all predictions by using information from the gene regulatory network. Here, we take advantage of this strategy that selects a large subset of predicted microRNA binding sites as more likely to possess less false-positives because of their over-representation in RE1 silencing transcription factor (REST)-regulated genes from the background of TargetScanHuman 7.2 predictions to identify useful features for the prediction of microRNA targets. These enriched miRNA families would have silencing activity for neural transcripts overlapping the repressive activity on neural genes of REST. We analyze properties of associated microRNA binding sites and contrast the outcome to the background. We found that the selected subset presents significant differences respect to the background: (i) lower GC-content in the vicinity of the predicted miRNA binding site, (ii) more target genes with multiple identical microRNA binding sites and (iii) a higher density of predicted microRNA binding sites close to the 3’ terminal end of the 3’-UTR. These results suggest that network selection of miRNA-mRNA pairs could provide useful features to improve microRNA target prediction.

Leveraging omic features with F3UTER enables identification of unannotated 3'UTRs for synaptic genes

There is growing evidence for the importance of 3' untranslated region (3'UTR) dependent regulatory processes. However, our current human 3'UTR catalogue is incomplete. Here, we develop a machine learning-based framework, leveraging both genomic and tissue-specific transcriptomic features to predict previously unannotated 3'UTRs. We identify unannotated 3'UTRs associated with 1,563 genes across 39 human tissues, with the greatest abundance found in the brain. These unannotated 3'UTRs are significantly enriched for RNA binding protein (RBP) motifs and exhibit high human lineage-specificity. We find that brain-specific unannotated 3'UTRs are enriched for the binding motifs of important neuronal RBPs such as TARDBP and RBFOX1, and their associated genes are involved in synaptic function. Our data is shared through an online resource F3UTER ( https://astx.shinyapps.io/F3UTER/ ). Overall, our data improves 3'UTR annotation and provides additional insights into the mRNA-RBP interactome in the human brain, with implications for our understanding of neurological and neurodevelopmental diseases.

Integrative approaches for studying the role of noncoding RNAs in influencing drug efficacy and toxicity

INTRODUCTION

Drug efficacy and toxicity are important factors for evaluation in drug development. Drug metabolizing enzymes and transporters (DMETs) play an essential role in drug efficacy and toxicity. Noncoding RNAs (ncRNAs) have been implicated to influence inter-individual variations in drug efficacy and safety by regulating DMETs. An efficient strategy is urgently needed to identify and functionally characterize ncRNAs that mediate drug efficacy and toxicity through regulating DMETs.

AREAS COVERED

We outline an integrative strategy to identify ncRNAs that modulate DMETs. We include reliable tools and databases for computational prediction of ncRNA targets with regard to their advantages and limitations. Various biochemical, molecular, and cellular assays are discussed for in vitro experimental verification of the regulatory function of ncRNAs. In vivo approaches for association of ncRNAs with drug treatment and toxicity are also reviewed.

EXPERT OPINION

A streamlined integration of computational prediction and wet-lab validation is important to elucidate mechanisms of ncRNAs in the regulation of DMETs related to drug efficacy and safety. Bioinformatic analyses using open-access tools and databases serve as a powerful booster for ncRNA Research in toxicology. Further refinement of computational algorithms and experimental technologies is needed to improve accuracy and efficiency in ncRNA target identification and characterization.

Comprehensive Identification and Profiling of miRNAs Involved in Terpenoid Synthesis of Gleditsia sinensis Lam

Gleditsia sinensis Lam. is a tree with worldwide distribution and important economic and medicinal values; its pods contain terpenoids including gleditsioside, thiamine, and brassinosteroids. However, thus far, there are few studies on the terpenoid regulation of G. sinensis at the molecular level. microRNA (miRNA) is a class of small RNAs with conserved and crucial roles in the regulation of diverse biological processes during plant growth and development. To identify the miRNAs of G. sinensis and evaluate their involvement in terpenoid synthesis, this investigation quantified the content changes in saponins in pods at three developmental stages: May (pod-setting stage), July (elongation stage), and September (browning stage), and then we performed genome-wide miRNA profiles during the three development stages of the G. sinensis pods. A total of 351 conserved miRNAs belonging to 216 families were identified, among which 36 conserved miRNAs exist specifically in legumes. Through target analysis, 708 unigenes were predicted to be candidate targets of 37 differentially expressed miRNAs. The targets of miR838-3p and miR2093-5p were involved in the derived branches of monoterpenes and gleditsioside, in brassinosteroid biosynthesis (BRB), and in indole alkaloid biosynthesis (IAB). Intriguingly, the targets of miR829-3p.1 were predicted to take part in thiamine biosynthesis, and the targets of miR4414b and miR5037a were involved in the main process of cytokinin synthesis. The corresponding targets participated in BRB, IAB, and terpenoid backbone biosynthesis, which were enriched significantly, suggesting that miR2093-5p, miR4414b, miR5037a, miR829-3p.1, and miR838-3p play indispensable roles in the regulation of triterpenoid saponin and monoterpenoid biosynthesis. To date, this is the first report of miRNA identification in G. sinensis and miRNA expression profiles at different developmental stages of G. sinensis pods, which provides a basis for further uncovering the molecular regulation of terpenoid synthesis in G. sinensis and new insights into the role of miRNAs in legumes.

Cigarette smoking-associated isoform switching and 3' UTR lengthening via alternative polyadenylation

Cigarette smoking induces a profound transcriptomic and systemic inflammatory response. Previous studies have focused on gene level differential expression of smoking, but the genome-wide effects of smoking on alternative isoform regulation have not yet been described. We conducted RNA sequencing in whole-blood samples of 454 current and 767 former smokers in the COPDGene Study, and we analyzed the effects of smoking on differential usage of isoforms and exons. At 10% FDR, we detected 3167 differentially expressed genes, 945 differentially used isoforms and 160 differentially used exons. Isoform switch analysis revealed widespread 3' UTR lengthening associated with cigarette smoking. The lengthening of these 3' UTRs was consistent with alternative usage of distal polyadenylation sites, and these extended 3' UTR regions were significantly enriched with functional sequence elements including microRNA and RNA-protein binding sites. These findings warrant further studies on alternative polyadenylation events as potential biomarkers and novel therapeutic targets for smoking-related diseases.

Differential regulation of mRNA fate by the human Ccr4-Not complex is driven by coding sequence composition and mRNA localization

BACKGROUND

Regulation of protein output at the level of translation allows for a rapid adaptation to dynamic changes to the cell's requirements. This precise control of gene expression is achieved by complex and interlinked biochemical processes that modulate both the protein synthesis rate and stability of each individual mRNA. A major factor coordinating this regulation is the Ccr4-Not complex. Despite playing a role in most stages of the mRNA life cycle, no attempt has been made to take a global integrated view of how the Ccr4-Not complex affects gene expression.

RESULTS

This study has taken a comprehensive approach to investigate post-transcriptional regulation mediated by the Ccr4-Not complex assessing steady-state mRNA levels, ribosome position, mRNA stability, and protein production transcriptome-wide. Depletion of the scaffold protein CNOT1 results in a global upregulation of mRNA stability and the preferential stabilization of mRNAs enriched for G/C-ending codons. We also uncover that mRNAs targeted to the ER for their translation have reduced translational efficiency when CNOT1 is depleted, specifically downstream of the signal sequence cleavage site. In contrast, translationally upregulated mRNAs are normally localized in p-bodies, contain disorder-promoting amino acids, and encode nuclear localized proteins. Finally, we identify ribosome pause sites that are resolved or induced by the depletion of CNOT1.

CONCLUSIONS

We define the key mRNA features that determine how the human Ccr4-Not complex differentially regulates mRNA fate and protein synthesis through a mechanism linked to codon composition, amino acid usage, and mRNA localization.

Alternative polyadenylation: An enigma of transcript length variation in health and disease

Alternative polyadenylation (APA) is a molecular mechanism during a pre-mRNA processing that involves usage of more than one polyadenylation site (PA-site) generating transcripts of varying length from a single gene. The location of a PA-site affects transcript length and coding potential of an mRNA contributing to both mRNA and protein diversification. This variation in the transcript length affects mRNA stability and translation, mRNA subcellular and tissue localization, and protein function. APA is now considered as an important regulatory mechanism in the pathophysiology of human diseases. An important consequence of the changes in the length of 3'-untranslated region (UTR) from disease-induced APA is altered protein expression. Yet, the relationship between 3'-UTR length and protein expression remains a paradox in a majority of diseases. Here, we review occurrence of APA, mechanism of PA-site selection, and consequences of transcript length variation in different diseases. Emerging evidence reveals coordinated involvement of core RNA processing factors including poly(A) polymerases in the PA-site selection in diseases-associated APAs. Targeting such APA regulators will be therapeutically significant in combating drug resistance in cancer and other complex diseases. This article is categorized under: RNA Processing > 3' End Processing RNA in Disease and Development > RNA in Disease Translation > Regulation.

Transite: A Computational Motif-Based Analysis Platform That Identifies RNA-Binding Proteins Modulating Changes in Gene Expression

RNA-binding proteins (RBPs) play critical roles in regulating gene expression by modulating splicing, RNA stability, and protein translation. Stimulus-induced alterations in RBP function contribute to global changes in gene expression, but identifying which RBPs are responsible for the observed changes remains an unmet need. Here, we present Transite, a computational approach that systematically infers RBPs influencing gene expression through changes in RNA stability and degradation. As a proof of principle, we apply Transite to RNA expression data from human patients with non-small-cell lung cancer whose tumors were sampled at diagnosis or after recurrence following treatment with platinum-based chemotherapy. Transite implicates known RBP regulators of the DNA damage response and identifies hnRNPC as a new modulator of chemotherapeutic resistance, which we subsequently validated experimentally. Transite serves as a framework for the identification of RBPs that drive cell-state transitions and adds additional value to the vast collection of publicly available gene expression datasets.

Krismer et al. present a computational approach to identify RNA-binding proteins (RBPs) that modulate post-transcriptional control of gene expression using RNA expression data as inputs. By applying this approach to publicly available patient datasets, they identify and experimentally confirm that the RBP hnRNPC contributes to chemotherapy resistance in lung cancer.

MiR-124 synergism with ELAVL3 enhances target gene expression to promote neuronal maturity

Neuron-enriched microRNAs (miRNAs), miR-9/9* and miR-124 (miR-9/9*-124), direct cell fate switching of human fibroblasts to neurons when ectopically expressed by repressing antineurogenic genes. How these miRNAs function after the repression of fibroblast genes for neuronal fate remains unclear. Here, we identified targets of miR-9/9*-124 as reprogramming cells activate the neuronal program and reveal the role of miR-124 that directly promotes the expression of its target genes associated with neuronal development and function. The mode of miR-124 as a positive regulator is determined by the binding of both AGO and a neuron-enriched RNA-binding protein, ELAVL3, to target transcripts. Although existing literature indicates that miRNA-ELAVL family protein interaction can result in either target gene up-regulation or down-regulation in a context-dependent manner, we specifically identified neuronal ELAVL3 as the driver for miR-124 target gene up-regulation in neurons. In primary human neurons, repressing miR-124 and ELAVL3 led to the down-regulation of genes involved in neuronal function and process outgrowth and cellular phenotypes of reduced inward currents and neurite outgrowth. Our results highlight the synergistic role between miR-124 and RNA-binding proteins to promote target gene regulation and neuronal function.

Snord116 Post-transcriptionally Increases Nhlh2 mRNA Stability: Implications for Human Prader-Willi Syndrome

The smallest genomic region causing Prader-Willi Syndrome (PWS) deletes the non-coding RNA SNORD116 cluster; however, the function of SNORD116 remains a mystery. Previous work in the field revealed the tantalizing possibility that expression of NHLH2, a gene previously implicated in both obesity and hypogonadism, was downregulated in PWS patients and differentiated stem cells. In silico RNA: RNA modeling identified several potential interaction domains between SNORD116 and NHLH2 mRNA. One of these interaction domains was highly conserved in most vertebrate NHLH2 mRNAs examined. A construct containing the Nhlh2 mRNA, including its 3'-UTR, linked to a c-myc tag was transfected into a hypothalamic neuron cell line in the presence and absence of exogenously-expressed Snord116. Nhlh2 mRNA expression was upregulated in the presence of Snord116 dependent on the length and type of 3'UTR used on the construct. Furthermore, use of actinomycin D to stop new transcription in N29/2 cells demonstrated that the upregulation occurred through increased stability of the Nhlh2 mRNA in the 45 minutes immediately following transcription. In silico modeling also revealed that a single nucleotide variant (SNV) in the NHLH2 mRNA could reduce the predicted interaction strength of the NHLH2:SNORD116 diad. Indeed, use of an Nhlh2 mRNA construct containing this SNV significantly reduces the ability of Snord116 to increase Nhlh2 mRNA levels. For the first time, these data identify a motif and mechanism for SNORD116-mediated regulation of NHLH2, clarifying the mechanism by which deletion of the SNORD116 snoRNAs locus leads to PWS phenotypes.

Translational remodeling by RNA-binding proteins and noncoding RNAs

Responsible for generating the proteome that controls phenotype, translation is the ultimate convergence point for myriad upstream signals that influence gene expression. System-wide adaptive translational reprogramming has recently emerged as a pillar of cellular adaptation. As classic regulators of mRNA stability and translation efficiency, foundational studies established the concept of collaboration and competition between RNA-binding proteins (RBPs) and noncoding RNAs (ncRNAs) on individual mRNAs. Fresh conceptual innovations now highlight stress-activated, evolutionarily conserved RBP networks and ncRNAs that increase the translation efficiency of populations of transcripts encoding proteins that participate in a common cellular process. The discovery of post-transcriptional functions for long noncoding RNAs (lncRNAs) was particularly intriguing given their cell-type-specificity and historical definition as nuclear-functioning epigenetic regulators. The convergence of RBPs, lncRNAs, and microRNAs on functionally related mRNAs to enable adaptive protein synthesis is a newer biological paradigm that highlights their role as "translatome (protein output) remodelers" and reinvigorates the paradigm of "RNA operons." Together, these concepts modernize our understanding of cellular stress adaptation and strategies for therapeutic development. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications Translation > Translation Regulation Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs.

CLIP and Massively Parallel Functional Analysis of CELF6 Reveal a Role in Destabilizing Synaptic Gene mRNAs through Interaction with 3' UTR Elements

CELF6 is a CELF-RNA-binding protein, and thus part of a protein family with roles in human disease; however, its mRNA targets in the brain are largely unknown. Using cross-linking immunoprecipitation and sequencing (CLIP-seq), we define its CNS targets, which are enriched for 3' UTRs in synaptic protein-coding genes. Using a massively parallel reporter assay framework, we test the consequence of CELF6 expression on target sequences, with and without mutating putative binding motifs. Where CELF6 exerts an effect on sequences, it is largely to decrease RNA abundance, which is reversed by mutating UGU-rich motifs. This is also the case for CELF3-5, with a protein-dependent effect on magnitude. Finally, we demonstrate that targets are derepressed in CELF6-mutant mice, and at least two key CNS proteins, FOS and FGF13, show altered protein expression levels and localization. Our works find, in addition to previously identified roles in splicing, that CELF6 is associated with repression of its CNS targets via the 3' UTR in vivo.

Hearing impairment due to Mir183/96/182 mutations suggests both loss and gain of function effects

The microRNA miR-96 is important for hearing, as point mutations in humans and mice result in dominant progressive hearing loss. Mir96 is expressed in sensory cells along with Mir182 and Mir183, but the roles of these closely-linked microRNAs are as yet unknown. Here we analyse mice carrying null alleles of Mir182, and of Mir183 and Mir96 together to investigate their roles in hearing. We found that Mir183/96 heterozygous mice had normal hearing and homozygotes were completely deaf with abnormal hair cell stereocilia bundles and reduced numbers of inner hair cell synapses at four weeks old. Mir182 knockout mice developed normal hearing then exhibited progressive hearing loss. Our transcriptional analyses revealed significant changes in a range of other genes, but surprisingly there were fewer genes with altered expression in the organ of Corti of Mir183/96 null mice compared with our previous findings in Mir96 Dmdo mutants, which have a point mutation in the miR-96 seed region. This suggests the more severe phenotype of Mir96 Dmdo mutants compared with Mir183/96 mutants, including progressive hearing loss in Mir96 Dmdo heterozygotes, is likely to be mediated by the gain of novel target genes in addition to the loss of its normal targets. We propose three mechanisms of action of mutant miRNAs; loss of targets that are normally completely repressed, loss of targets whose transcription is normally buffered by the miRNA, and gain of novel targets. Any of these mechanisms could lead to a partial loss of a robust cellular identity and consequent dysfunction.

CLIP and Massively Parallel Functional Analysis of CELF6 Reveal a Role in Destabilizing Synaptic Gene mRNAs through Interaction with 3' UTR Elements

CELF6 is a CELF-RNA-binding protein, and thus part of a protein family with roles in human disease; however, its mRNA targets in the brain are largely unknown. Using cross-linking immunoprecipitation and sequencing (CLIP-seq), we define its CNS targets, which are enriched for 3′ UTRs in synaptic protein-coding genes. Using a massively parallel reporter assay framework, we test the consequence of CELF6 expression on target sequences, with and without mutating putative binding motifs. Where CELF6 exerts an effect on sequences, it is largely to decrease RNA abundance, which is reversed by mutating UGU-rich motifs. This is also the case for CELF3–5, with a protein-dependent effect on magnitude. Finally, we demonstrate that targets are derepressed in CELF6-mutant mice, and at least two key CNS proteins, FOS and FGF13, show altered protein expression levels and localization. Our works find, in addition to previously identified roles in splicing, that CELF6 is associated with repression of its CNS targets via the 3′ UTR in vivo.

Rieger et al. assay the function of the RNA-binding protein CELF6 by defining its targets in the brain. They show that CELF6 largely binds 3′ UTRs of synaptic mRNAs. Using a massively parallel reporter assay, they further show that CELF6 and other CELFs are associated with lower mRNA abundance and that targets are derepressed in Celf6-knockout mice in vivo.

Hsa-miR-27a-3p and epidermal growth factor receptor expression analysis in glioblastoma FFPE samples

AIM

Glioblastoma multiforme (GBM) is the most invasive type of glial tumors. MicroRNAs as the small, noncoding RNAs have been revealed that play an important role in tumorigenic processes. So, they may be used as potential biomarkers for detection and prognosis of cancers at the early stages. In addition, they can be applied as therapeutic targets. In the present study, the expression levels of hsa-miR-27a-3p and EGFR were investigated in GBM.

METHODS

Real-time RT-PCR was applied to evaluate hsa-miR-27a-3p and EGFR expressions in the formalin-fixed paraffin-embedded (FFPE) tissue samples obtained from 50 GBM and 50 normal people.

RESULTS

The expression level of hsa-miR-27a-3p and EGFR was significantly different between cases and controls. Positive association was found between gene expressions and immunohistochemistry markers, such as Ki67 and glial fibrillary acidic protein, except for IDH1 status.

CONCLUSION

We showed the association of hsa-miR-27a-3p and EGFR with GBM and it can be concluded that they have a promising potential to use as primary cancer biomarkers.

Structural basis for the recognition of transiently structured AU-rich elements by Roquin

Adenylate/uridylate-rich elements (AREs) are the most common cis-regulatory elements in the 3′-untranslated region (UTR) of mRNAs, where they fine-tune turnover by mediating mRNA decay. They increase plasticity and efficacy of mRNA regulation and are recognized by several ARE-specific RNA-binding proteins (RBPs). Typically, AREs are short linear motifs with a high content of complementary A and U nucleotides and often occur in multiple copies. Although thermodynamically rather unstable, the high AU-content might enable transient secondary structure formation and modify mRNA regulation by RBPs. We have recently suggested that the immunoregulatory RBP Roquin recognizes folded AREs as constitutive decay elements (CDEs), resulting in shape-specific ARE-mediated mRNA degradation. However, the structural evidence for a CDE-like recognition of AREs by Roquin is still lacking. We here present structures of CDE-like folded AREs, both in their free and protein-bound form. Moreover, the AREs in the UCP3 3′-UTR are additionally bound by the canonical ARE-binding protein AUF1 in their linear form, adopting an alternative binding-interface compared to the recognition of their CDE structure by Roquin. Strikingly, our findings thus suggest that AREs can be recognized in multiple ways, allowing control over mRNA regulation by adapting distinct conformational states, thus providing differential accessibility to regulatory RBPs.

GIGYF1/2-Driven Cooperation between ZNF598 and TTP in Posttranscriptional Regulation of Inflammatory Signaling

Inflammatory signaling is restricted through degradation and the translational repression of cytokine mRNAs. A key factor in this regulation is tristetraprolin (TTP), an RNA-binding protein (RBP) that recruits RNA-destabilizing factors and the translation inhibitory complex 4EHP-GIGYF1/2 to AU-rich element (ARE)-containing mRNAs. Here, we show that the RBP ZNF598 contributes to the same regulatory module in a TTP-like manner. Similar to TTP, ZNF598 harbors three proline-rich motifs that bind the GYF domain of GIGYF1. RNA sequencing experiments showed that ZNF598 is required for the regulation of known TTP targets, including IL-8 and CSF2 mRNA. Furthermore, we demonstrate that ZNF598 binds to IL-8 mRNA, but not TNF mRNA. Collectively, our findings highlight that ZNF598 functions as an RBP that buffers the level of a range of mRNAs. We propose that ZNF598 is a TTP-like factor that can contribute to the regulation of the inflammatory potential of cytokine-producing cells.

Deep neural networks for interpreting RNA-binding protein target preferences

Deep learning has become a powerful paradigm to analyze the binding sites of regulatory factors including RNA-binding proteins (RBPs), owing to its strength to learn complex features from possibly multiple sources of raw data. However, the interpretability of these models, which is crucial to improve our understanding of RBP binding preferences and functions, has not yet been investigated in significant detail. We have designed a multitask and multimodal deep neural network for characterizing in vivo RBP targets. The model incorporates not only the sequence but also the region type of the binding sites as input, which helps the model to boost the prediction performance. To interpret the model, we quantified the contribution of the input features to the predictive score of each RBP. Learning across multiple RBPs at once, we are able to avoid experimental biases and to identify the RNA sequence motifs and transcript context patterns that are the most important for the predictions of each individual RBP. Our findings are consistent with known motifs and binding behaviors and can provide new insights about the regulatory functions of RBPs.

Development of an AAV-Based MicroRNA Gene Therapy to Treat Machado-Joseph Disease

Spinocerebellar ataxia type 3 (SCA3), or Machado-Joseph disease (MJD), is a progressive neurodegenerative disorder caused by a CAG expansion in the ATXN3 gene. The expanded CAG repeat is translated into a prolonged polyglutamine repeat in the ataxin-3 protein and accumulates within inclusions, acquiring toxic properties, which results in degeneration of the cerebellum and brain stem. In the current study, a non-allele-specific ATXN3 silencing approach was investigated using artificial microRNAs engineered to target various regions of the ATXN3 gene (miATXN3). The miATXN3 candidates were screened in vitro based on their silencing efficacy on a luciferase (Luc) reporter co-expressing ATXN3. The three best miATXN3 candidates were further tested for target engagement and potential off-target activity in induced pluripotent stem cells (iPSCs) differentiated into frontal brain-like neurons and in a SCA3 knockin mouse model. Besides a strong reduction of ATXN3 mRNA and protein, small RNA sequencing revealed efficient guide strand processing without passenger strands being produced. We used different methods to predict alteration of off-target genes upon AAV5-miATXN3 treatment and found no evidence for unwanted effects. Furthermore, we demonstrated in a large animal model, the minipig, that intrathecal delivery of AAV5 can transduce the main areas affected in SCA3 patients. These results proved a strong basis to move forward to investigate distribution, efficacy, and safety of AAV5-miATXN3 in large animals.

Alternative polyadenylation of mRNA and its role in cancer

Alternative polyadenylation (APA) is a molecular process that generates diversity at the 3′ end of RNA polymerase II transcripts from over 60% of human genes. APA is derived from the existence of multiple polyadenylation signals (PAS) within the same transcript, and results in the differential inclusion of sequence information at the 3′ end. While APA can occur between two PASs allowing for generation of transcripts with distinct coding potential from a single gene, most APA occurs within the untranslated region (3′UTR) and changes the length and content of these non-coding sequences. APA within the 3′UTR can have tremendous impact on its regulatory potential of the mRNA through a variety of mechanisms, and indeed this layer of gene expression regulation has profound impact on processes vital to cell growth and development. Recent studies have particularly highlighted the importance of APA dysregulation in cancer onset and progression. Here, we review the current knowledge of APA and its impacts on mRNA stability, translation, localization and protein localization. We also discuss the implications of APA dysregulation in cancer research and therapy.

Lessons from studying the AU-rich elements in chronic inflammation and autoimmunity

AU-rich elements (AREs) comprise one of the most widely studied families of regulatory RNA structures met in RNAs engaged in complex immunological reactions. A multitude of genetic, molecular, holistic and functional studies have been utilized for the analyses of the AREs and their interactions to proteins that bind to them. Data stemming from these studies brought forth a world of RNA-related check-points against infection, chronic inflammation, tumor associated immunity, and autoimmunity; and the interest to capitalize the interactions of AREs for clinical management and therapy. They also provided lessons on the cellular capabilities of post-transcriptional control. Originally thought as transcript-restricted regulators of turnover and translation, ARE-binding proteins do in fact harbor great versatility and interactivity across nuclear and cytoplasmic compartments; and act as functional coordinators of immune-cellular programs. Harnessing these deterministic functions requires extensive knowledge of their synergies or antagonisms at a cell-specific level; but holds great promise since it can provide the efficacy of combinatorial therapies with single agents.

Alternative cleavage and polyadenylation in health and disease

Most human genes have multiple sites at which RNA 3' end cleavage and polyadenylation can occur, enabling the expression of distinct transcript isoforms under different conditions. Novel methods to sequence RNA 3' ends have generated comprehensive catalogues of polyadenylation (poly(A)) sites; their analysis using innovative computational methods has revealed how poly(A) site choice is regulated by core RNA 3' end processing factors, such as cleavage factor I and cleavage and polyadenylation specificity factor, as well as by other RNA-binding proteins, particularly splicing factors. Here, we review the experimental and computational methods that have enabled the global mapping of mRNA and of long non-coding RNA 3' ends, quantification of the resulting isoforms and the discovery of regulators of alternative cleavage and polyadenylation (APA). We highlight the different types of APA-derived isoforms and their functional differences, and illustrate how APA contributes to human diseases, including cancer and haematological, immunological and neurological diseases.

Ciphers and Executioners: How 3'-Untranslated Regions Determine the Fate of Messenger RNAs

The sequences and structures of 3'-untranslated regions (3'UTRs) of messenger RNAs govern their stability, localization, and expression. 3'UTR regulatory elements are recognized by a wide variety of trans-acting factors that include microRNAs (miRNAs), their associated machinery, and RNA-binding proteins (RBPs). In turn, these factors instigate common mechanistic strategies to execute the regulatory programs encoded by 3'UTRs. Here, we review classes of factors that recognize 3'UTR regulatory elements and the effector machineries they guide toward mRNAs to dictate their expression and fate. We outline illustrative examples of competitive, cooperative, and coordinated interplay such as mRNA localization and localized translation. We further review the recent advances in the study of mRNP granules and phase transition, and their possible significance for the functions of 3'UTRs. Finally, we highlight some of the most recent strategies aimed at deciphering the complexity of the regulatory codes of 3'UTRs, and identify some of the important remaining challenges.

Antagonistic and cooperative AGO2-PUM interactions in regulating mRNAs

Approximately 1500 RNA-binding proteins (RBPs) profoundly impact mammalian cellular function by controlling distinct sets of transcripts, often using sequence-specific binding to 3′ untranslated regions (UTRs) to regulate mRNA stability and translation. Aside from their individual effects, higher-order combinatorial interactions between RBPs on specific mRNAs have been proposed to underpin the regulatory network. To assess the extent of such co-regulatory control, we took a global experimental approach followed by targeted validation to examine interactions between two well-characterized and highly conserved RBPs, Argonaute2 (AGO2) and Pumilio (PUM1 and PUM2). Transcriptome-wide changes in AGO2-mRNA binding upon PUM knockdown were quantified by CLIP-seq, and the presence of PUM binding on the same 3′UTR corresponded with cooperative and antagonistic effects on AGO2 occupancy. In addition, PUM binding sites that overlap with AGO2 showed differential, weakened binding profiles upon abrogation of AGO2 association, indicative of cooperative interactions. In luciferase reporter validation of candidate 3′UTR sites where AGO2 and PUM colocalized, three sites were identified to host antagonistic interactions, where PUM counteracts miRNA-guided repression. Interestingly, the binding sites for the two proteins are too far for potential antagonism due to steric hindrance, suggesting an alternate mechanism. Our data experimentally confirms the combinatorial regulatory model and indicates that the mostly repressive PUM proteins can change their behavior in a context-dependent manner. Overall, the approach underscores the importance of further elucidation of complex interactions between RBPs and their transcriptome-wide extent.

Terminal exon characterization with TECtool reveals an abundance of cell-specific isoforms

Sequencing of RNA 3' ends has uncovered numerous sites that do not correspond to the termination sites of known transcripts. Through their 3' untranslated regions, protein-coding RNAs interact with RNA-binding proteins and microRNAs, which regulate many properties, including RNA stability and subcellular localization. We developed the terminal exon characterization (TEC) tool ( http://tectool.unibas.ch ), which can be used with RNA-sequencing data from any species for which a genome annotation that includes sites of RNA cleavage and polyadenylation is available. We discovered hundreds of previously unknown isoforms and cell-type-specific terminal exons in human cells. Ribosome profiling data revealed that many of these isoforms were translated. By applying TECtool to single-cell sequencing data, we found that the newly identified isoforms were expressed in subpopulations of cells. Thus, TECtool enables the identification of previously unknown isoforms in well-studied cell systems and in rare cell types.

Cryptic sequence features in the active postmortem transcriptome

BACKGROUND

Our previous study found that more than 500 transcripts significantly increased in abundance in the zebrafish and mouse several hours to days postmortem relative to live controls. The current literature suggests that most mRNAs are post-transcriptionally regulated in stressful conditions. We rationalized that the postmortem transcripts must contain sequence features (3- to 9- mers) that are unique from those in the rest of the transcriptome and that these features putatively serve as binding sites for proteins and/or non-coding RNAs involved in post-transcriptional regulation.

RESULTS

We identified 5117 and 2245 over-represented sequence features in the mouse and zebrafish, respectively, which represents less than 1.5% of all possible features. Some of these features were disproportionately distributed along the transcripts with high densities in the 3' untranslated regions of the zebrafish (0.3 mers/nt) and the open reading frames of the mouse (0.6 mers/nt). Yet, the highest density (2.3 mers/nt) occurred in the open reading frames of 11 mouse transcripts that lacked 3' or 5' untranslated regions. These results suggest the transcripts with high density of features might serve as 'molecular sponges' that sequester RNA binding proteins and/or microRNAs, and thus indirectly increase the stability and gene expression of other transcripts. In addition, some of the features were identified as binding sites for Rbfox and Hud proteins that are also involved in increasing transcript stability and gene expression.

CONCLUSIONS

Our results are consistent with the hypothesis that transcripts involved in responding to extreme stress, such as organismal death, have sequence features that make them different from the rest of the transcriptome. Some of these features serve as putative binding sites for proteins and non-coding RNAs that determine transcript stability and fate. A small number of the transcripts have high density sequence features, which are presumably involved in sequestering RNA binding proteins and microRNAs and thus preventing regulatory interactions among other transcripts. Our results provide baseline data on post-transcriptional regulation in stressful conditions that has implications for regulation in disease, starvation, and cancer.

Trans-acting translational regulatory RNA binding proteins

The canonical molecular machinery required for global mRNA translation and its control has been well defined, with distinct sets of proteins involved in the processes of translation initiation, elongation and termination. Additionally, noncanonical, trans-acting regulatory RNA-binding proteins (RBPs) are necessary to provide mRNA-specific translation, and these interact with 5' and 3' untranslated regions and coding regions of mRNA to regulate ribosome recruitment and transit. Recently it has also been demonstrated that trans-acting ribosomal proteins direct the translation of specific mRNAs. Importantly, it has been shown that subsets of RBPs often work in concert, forming distinct regulatory complexes upon different cellular perturbation, creating an RBP combinatorial code, which through the translation of specific subsets of mRNAs, dictate cell fate. With the development of new methodologies, a plethora of novel RNA binding proteins have recently been identified, although the function of many of these proteins within mRNA translation is unknown. In this review we will discuss these methodologies and their shortcomings when applied to the study of translation, which need to be addressed to enable a better understanding of trans-acting translational regulatory proteins. Moreover, we discuss the protein domains that are responsible for RNA binding as well as the RNA motifs to which they bind, and the role of trans-acting ribosomal proteins in directing the translation of specific mRNAs. This article is categorized under: RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes Translation > Translation Regulation Translation > Translation Mechanisms.

PTRE-seq reveals mechanism and interactions of RNA binding proteins and miRNAs

RNA binding proteins (RBP) and microRNAs (miRNAs) often bind sequences in 3' untranslated regions (UTRs) of mRNAs, and regulate stability and translation efficiency. With the identification of numerous RBPs and miRNAs, there is an urgent need for new technologies to dissect the function of the cis-acting elements of RBPs and miRNAs. We describe post-transcriptional regulatory element sequencing (PTRE-seq), a massively parallel method for assaying the target sequences of miRNAs and RBPs. We use PTRE-seq to dissect sequence preferences and interactions between miRNAs and RBPs. The binding sites for these effector molecules influenced different aspects of the RNA lifecycle: RNA stability, translation efficiency, and translation initiation. In some cases, post-transcriptional control is modular, with different factors acting independently of each other, while in other cases factors show specific epistatic interactions. The throughput, flexibility, and reproducibility of PTRE-seq make it a valuable tool to study post-transcriptional regulation by 3'UTR elements.

Transcriptome-wide Identification and Validation of Interactions between the miRNA Machinery and HuR on mRNA Targets

The 3' untranslated region (UTR) of mRNAs is the primary regulatory region that mediates post-transcriptional control by microRNAs and RNA-binding proteins in the cytoplasm. Aside from individual sequence-specific binding and regulation, examples of interaction between these factors at particular 3' UTR sites have emerged. However, the whole picture of such higher-order regulatory modules across the transcriptome is lacking. Here, we investigate the interactions between HuR, a ubiquitous RNA-binding protein, and Ago2, a core effector of the miRNA pathway, at the transcriptome-wide level. Using HITS-CLIP, we map HuR and miRNA binding sites on human 3' UTRs and assess their co-occurrence. In addition, we demonstrate global effects of HuR knockdown on Ago2 occupancy, suggesting a co-regulatory relationship. Focusing on sites of Ago2-HuR overlap, 13 candidates were screened in luciferase reporter assays. Eleven sites showed miRNA-dependent repression, as confirmed in Dicer-null cells. To test for HuR's role in co-regulation, we measured the reporters in HuR KO cells. Three of the miRNA sites demonstrated altered activities, indicating that HuR has an effect on miRNA repression at those sites. Our study presents an efficient search and validation system for studying miRNA-HuR interactions, which expands our understanding of the combinatorial post-transcriptional control of gene expression at the 3' UTR.

Handshakes and Fights: The Regulatory Interplay of RNA-Binding Proteins

What drives the flow of signals controlling the outcome of post-transcriptional regulation of gene expression? This regulatory layer, presiding to processes ranging from splicing to mRNA stability and localization, is a key determinant of protein levels and thus cell phenotypes. RNA-binding proteins (RBPs) form a remarkable army of post-transcriptional regulators, strong of more than 1,500 genes implementing this expression fine-tuning plan and implicated in both cell physiology and pathology. RBPs can bind and control a wide array of RNA targets. This sheer amount of interactions form complex regulatory networks (PTRNs) where the action of individual RBPs cannot be easily untangled from each other. While past studies have mostly focused on the action of individual RBPs on their targets, we are now observing an increasing amount of evidence describing the occurrence of interactions between RBPs, defining how common target RNAs are regulated. This suggests that the flow of signals in PTRNs is driven by the intertwined contribution of multiple RBPs, concurrently acting on each of their targets. Understanding how RBPs cooperate and compete is thus of paramount importance to chart the wiring of PTRNs and their impact on cell phenotypes. Here we review the current knowledge about patterns of RBP interaction and attempt at describing their general principles. We also discuss future directions which should be taken to reach a comprehensive understanding of this fundamental aspect of gene expression regulation.