Disentangling genetic effects on transcriptional and post-transcriptional gene regulation through integrating exon and intron expression QTLs.

Expression quantitative trait loci (eQTL) studies typically consider exon expression of genes and discard intronic RNA sequencing reads despite their information on RNA metabolism. Here, we quantify genetic effects on exon and intron levels of genes and their ratio in lymphoblastoid cell lines, revealing thousands of cis-QTLs of each type. While genetic effects are often shared between cis-QTL types, 7814 (47%) are not detected as top cis-QTLs at exon levels. We show that exon levels preferentially capture genetic effects on transcriptional regulation, while exon-intron-ratios better detect those on co- and post-transcriptional processes. Considering all cis-QTL types substantially increases (by 71%) the number of colocalizing variants identified by genome-wide association studies (GWAS). It further allows dissecting the potential gene regulatory processes underlying GWAS associations, suggesting comparable contributions by transcriptional (50%) and co- and post-transcriptional regulation (46%) to complex traits. Overall, integrating intronic RNA sequencing reads in eQTL studies expands our understanding of genetic effects on gene regulatory processes.

Analysis of Eukaryotic lincRNA Sequences Indicates Signatures of Hindered Translation Linked to Selection Pressure.

Long intergenic noncoding RNAs (lincRNAs) represent a large fraction of transcribed loci in eukaryotic genomes. Although classified as noncoding, most lincRNAs contain open reading frames (ORFs), and it remains unclear why cytoplasmic lincRNAs are not or very inefficiently translated. Here, we analyzed signatures of hindered translation in lincRNA sequences from five eukaryotes, covering a range of natural selection pressures. In fission yeast and Caenorhabditis elegans, that is, species under strong selection, we detected significantly shorter ORFs, a suboptimal sequence context around start codons for translation initiation, and trinucleotides ("codons") corresponding to less abundant tRNAs than for neutrally evolving control sequences, likely impeding translation elongation. For human, we detected signatures for cell-type-specific hindrance of lincRNA translation, in particular codons in abundant cytoplasmic lincRNAs corresponding to lower expressed tRNAs than control codons, in three out of five human cell lines. We verified that varying tRNA expression levels between cell lines are reflected in the amount of ribosomes bound to cytoplasmic lincRNAs in each cell line. We further propose that codons at ORF starts are particularly important for reducing ribosome-binding to cytoplasmic lincRNA ORFs. Altogether, our analyses indicate that in species under stronger selection lincRNAs evolved sequence features generally hindering translation and support cell-type-specific hindrance of translation efficiency in human lincRNAs. The sequence signatures we have identified may improve predicting peptide-coding and genuine noncoding lincRNAs in a cell type.

Untargeted Metabolome- and Transcriptome-Wide Association Study Suggests Causal Genes Modulating Metabolite Concentrations in Urine.

Gene products can affect the concentrations of small molecules (aka "metabolites"), and conversely, some metabolites can modulate the concentrations of gene transcripts. While many specific instances of this interplay have been revealed, a global approach to systematically uncover human gene-metabolite interactions is still lacking. We performed a metabolome- and transcriptome-wide association study to identify genes influencing the human metabolome using untargeted metabolome features, extracted from 1H nuclear magnetic resonance spectroscopy (NMR) of urine samples, and gene expression levels, quantified from RNA-Seq of lymphoblastoid cell lines (LCL) from 555 healthy individuals. We identified 20 study-wide significant associations corresponding to 15 genes, of which 5 associations (with 2 genes) were confirmed with follow-up NMR data. Using metabomatching, we identified the metabolites corresponding to metabolome features associated with the genes, namely, N-acetylated compounds with ALMS1 and trimethylamine (TMA) with HPS1. Finally, Mendelian randomization analysis supported a potential causal link between the expression of genes in both the ALMS1- and HPS1-loci and their associated metabolite concentrations. In the case of HPS1, we additionally observed that TMA concentration likely exhibits a reverse causal effect on HPS1 expression levels, indicating a negative feedback loop. Our study highlights how the integration of metabolomics, gene expression, and genetic data can pinpoint causal genes modulating metabolite concentrations.

Inosine Substitutions in RNA Activate Latent G-Quadruplexes.

It is well-accepted that gene expression is heavily influenced by RNA structure. For instance, stem-loops and G-quadruplexes (rG4s) are dynamic motifs in mRNAs that influence gene expression. Adenosine-to-inosine (A-to-I) editing is a common chemical modification of RNA which introduces a nucleobase that is iso-structural with guanine, thereby changing RNA base-pairing properties. Here, we provide biophysical, chemical, and biological evidence that A-to-I exchange can activate latent rG4s by filling incomplete G-quartets with inosine. We demonstrate the formation of inosine-containing rG4s (GI-quadruplexes) in vitro and verify their activity in cells. GI-quadruplexes adopt parallel topologies, stabilized by potassium ions. They exhibit moderately reduced thermal stability compared to conventional G-quadruplexes. To study inosine-induced structural changes in a naturally occurring RNA, we use a synthetic approach that enables site-specific inosine incorporation in long RNAs. In summary, RNA GI-quadruplexes are a previously unrecognized structural motif that may contribute to the regulation of gene expression in vivo.

Temporal Control of the TGF-ß Signaling Network by Mouse ESC MicroRNA Targets of Different Affinities.

Although microRNAs (miRNAs) function in the control of embryonic stem cell (ESC) pluripotency, a systems-level understanding is still being developed. Through the analysis of progressive Argonaute (Ago)-miRNA depletion and rescue, including stable Ago knockout mouse ESCs, we uncover transforming growth factor beta (TGF-ß) pathway activation as a direct and early response to ESC miRNA reduction. Mechanistically, we link the derepression of weaker miRNA targets, including TGF-ß receptor 1 (Tgfbr1), to the sensitive TGF-ß pathway activation. In contrast, stronger miRNA targets impart a more robust repression, which dampens concurrent transcriptional activation. We verify such dampened induction for TGF-ß antagonist Lefty. We find that TGF-ß pathway activation contributes to the G1 cell-cycle accumulation of miRNA-deficient ESCs. We propose that miRNA target affinity is a determinant of the temporal response to miRNA changes, which enables the coordination of gene network responses.

MINA-1 and WAGO-4 are part of regulatory network coordinating germ cell death and RNAi in C. elegans.

Post-transcriptional control of mRNAs by RNA-binding proteins (RBPs) has a prominent role in the regulation of gene expression. RBPs interact with mRNAs to control their biogenesis, splicing, transport, localization, translation, and stability. Defects in such regulation can lead to a wide range of human diseases from neurological disorders to cancer. Many RBPs are conserved between Caenorhabditis elegans and humans, and several are known to regulate apoptosis in the adult C. elegans germ line. How these RBPs control apoptosis is, however, largely unknown. Here, we identify mina-1(C41G7.3) in a RNA interference-based screen as a novel regulator of apoptosis, which is exclusively expressed in the adult germ line. The absence of MINA-1 causes a dramatic increase in germ cell apoptosis, a reduction in brood size, and an impaired P granules organization and structure. In vivo crosslinking immunoprecipitation experiments revealed that MINA-1 binds a set of mRNAs coding for RBPs associated with germ cell development. Additionally, a system-wide analysis of a mina-1 deletion mutant compared with wild type, including quantitative proteome and transcriptome data, hints to a post-transcriptional regulatory RBP network driven by MINA-1 during germ cell development in C. elegans. In particular, we found that the germline-specific Argonaute WAGO-4 protein levels are increased in mina-1 mutant background. Phenotypic analysis of double mutant mina-1;wago-4 revealed that contemporary loss of MINA-1 and WAGO-4 strongly rescues the phenotypes observed in mina-1 mutant background. To strengthen this functional interaction, we found that upregulation of WAGO-4 in mina-1 mutant animals causes hypersensitivity to exogenous RNAi. Our comprehensive experimental approach allowed us to describe a phenocritical interaction between two RBPs controlling germ cell apoptosis and exogenous RNAi. These findings broaden our understanding of how RBPs can orchestrate different cellular events such as differentiation and death in C. elegans.

Regulation of RNA editing by RNA-binding proteins in human cells.

Adenosine-to-inosine (A-to-I) editing, mediated by the ADAR enzymes, diversifies the transcriptome by altering RNA sequences. Recent studies reported global changes in RNA editing in disease and development. Such widespread editing variations necessitate an improved understanding of the regulatory mechanisms of RNA editing. Here, we study the roles of >200 RNA-binding proteins (RBPs) in mediating RNA editing in two human cell lines. Using RNA-sequencing and global protein-RNA binding data, we identify a number of RBPs as key regulators of A-to-I editing. These RBPs, such as TDP-43, DROSHA, NF45/90 and Ro60, mediate editing through various mechanisms including regulation of ADAR1 expression, interaction with ADAR1, and binding to Alu elements. We highlight that editing regulation by Ro60 is consistent with the global up-regulation of RNA editing in systemic lupus erythematosus. Additionally, most key editing regulators act in a cell type-specific manner. Together, our work provides insights for the regulatory mechanisms of RNA editing.

Structure-mediated modulation of mRNA abundance by A-to-I editing.

RNA editing introduces single nucleotide changes to RNA, thus potentially diversifying gene expression. Recent studies have reported significant changes in RNA editing profiles in disease and development. The functional consequences of these widespread alterations remain elusive because of the unknown function of most RNA editing sites. Here, we carry out a comprehensive analysis of A-to-I editomes in human populations. Surprisingly, we observe highly similar editing profiles across populations despite striking differences in the expression levels of ADAR genes. Striving to explain this discrepancy, we uncover a functional mechanism of A-to-I editing in regulating mRNA abundance. We show that A-to-I editing stabilizes RNA secondary structures and reduces the accessibility of AGO2-miRNA to target sites in mRNAs. The editing-dependent stabilization of mRNAs in turn alters the observed editing levels in the stable RNA repertoire. Our study provides valuable insights into the functional impact of RNA editing in human cells.

Hypothalamic transcriptomes of 99 mouse strains reveal trans eQTL hotspots, splicing QTLs and novel non-coding genes.

Previous studies had shown that the integration of genome wide expression profiles, in metabolic tissues, with genetic and phenotypic variance, provided valuable insight into the underlying molecular mechanisms. We used RNA-Seq to characterize hypothalamic transcriptome in 99 inbred strains of mice from the Hybrid Mouse Diversity Panel (HMDP), a reference resource population for cardiovascular and metabolic traits. We report numerous novel transcripts supported by proteomic analyses, as well as novel non coding RNAs. High resolution genetic mapping of transcript levels in HMDP, reveals both local and trans expression Quantitative Trait Loci (eQTLs) demonstrating 2 trans eQTL 'hotspots' associated with expression of hundreds of genes. We also report thousands of alternative splicing events regulated by genetic variants. Finally, comparison with about 150 metabolic and cardiovascular traits revealed many highly significant associations. Our data provide a rich resource for understanding the many physiologic functions mediated by the hypothalamus and their genetic regulation.

Cooperative target mRNA destabilization and translation inhibition by miR-58 microRNA family in C. elegans.

In animals, microRNAs frequently form families with related sequences. The functional relevance of miRNA families and the relative contribution of family members to target repression have remained, however, largely unexplored. Here, we used the Caenorhabditis elegans miR-58 miRNA family, composed primarily of the four highly abundant members miR-58.1, miR-80, miR-81, and miR-82, as a model to investigate the redundancy of miRNA family members and their impact on target expression in an in vivo setting. We found that miR-58 family members repress largely overlapping sets of targets in a predominantly additive fashion. Progressive deletions of miR-58 family members lead to cumulative up-regulation of target protein and RNA levels. Phenotypic defects could only be observed in the family quadruple mutant, which also showed the strongest change in target protein levels. Interestingly, although the seed sequences of miR-80 and miR-58.1 differ in a single nucleotide, predicted canonical miR-80 targets were efficiently up-regulated in the mir-58.1 single mutant, indicating functional redundancy of distinct members of this miRNA family. At the aggregate level, target binding leads mainly to mRNA degradation, although we also observed some degree of translational inhibition, particularly in the single miR-58 family mutants. These results provide a framework for understanding how miRNA family members interact to regulate target mRNAs.

miR-CLIP capture of a miRNA targetome uncovers a lincRNA H19-miR-106a interaction.

Identifying the interaction partners of noncoding RNAs is essential for elucidating their functions. We have developed an approach, termed microRNA crosslinking and immunoprecipitation (miR-CLIP), using pre-miRNAs modified with psoralen and biotin to capture their targets in cells. Photo-crosslinking and Argonaute 2 immunopurification followed by streptavidin affinity purification of probe-linked RNAs provided selectivity in the capture of targets, which were identified by deep sequencing. miR-CLIP with pre-miR-106a, a miR-17-5p family member, identified hundreds of putative targets in HeLa cells, many carrying conserved sequences complementary to the miRNA seed but also many that were not predicted computationally. miR-106a overexpression experiments confirmed that miR-CLIP captured functional targets, including H19, a long noncoding RNA that is expressed during skeletal muscle cell differentiation. We showed that miR-17-5p family members bind H19 in HeLa cells and myoblasts. During myoblast differentiation, levels of H19, miR-17-5p family members and mRNA targets changed in a manner suggesting that H19 acts as a 'sponge' for these miRNAs.

Structural and functional implications of the QUA2 domain on RNA recognition by GLD-1.

The STAR family comprises ribonucleic acid (RNA)-binding proteins that play key roles in RNA-regulatory processes. RNA recognition is achieved by a KH domain with an additional α-helix (QUA2) that seems to extend the RNA-binding surface to six nucleotides for SF1 (Homo sapiens) and seven nucleotides for GLD-1 (Caenorhabditis elegans). To understand the structural basis of this probable difference in specificity, we determined the solution structure of GLD-1 KH-QUA2 with the complete consensus sequence identified in the tra-2 gene. Compared to SF1, the GLD-1 KH-QUA2 interface adopts a different conformation resulting indeed in an additional sequence-specific binding pocket for a uracil at the 5'end. The functional relevance of this binding pocket is emphasized by our bioinformatics analysis showing that GLD-1 binding sites with this 5'end uracil are more predictive for the functional response of the messenger RNAs to gld-1 knockout. We further reveal the importance of the KH-QUA2 interface in vitro and that its alteration in vivo affects the level of translational repression dependent on the sequence of the GLD-1 binding motif. In conclusion, we demonstrate that the QUA2 domain distinguishes GLD-1 from other members of the STAR family and contributes more generally to the modulation of RNA-binding affinity and specificity of KH domain containing proteins.

MicroRNA binding sites in the coding region of mRNAs: extending the repertoire of post-transcriptional gene regulation.

It is well established that microRNAs (miRNAs) induce mRNA degradation by binding to 3' untranslated regions (UTRs). The functionality of sites in the coding domain sequence (CDS), on the other hand, remains under discussion. Such sites have limited impact on target mRNA abundance and recent work suggests that miRNAs bind in the CDS to inhibit translation. What then could be the regulatory benefits of translation inhibition through CDS targeting compared to mRNA degradation following 3' UTR binding? We propose that these domain-dependent effects serve to diversify the functional repertoire of post-transcriptional gene expression control. Possible regulatory benefits may include tuning the time-scale and magnitude of post-transcriptional regulation, regulating protein abundance depending on or independently of the cellular state, and regulation of the protein abundance of alternative splice variants. Finally, we review emerging evidence that these ideas may generalize to RNA-binding proteins beyond miRNAs and Argonaute proteins.

Modeling the binding specificity of the RNA-binding protein GLD-1 suggests a function of coding region-located sites in translational repression.

To understand the function of the hundreds of RNA-binding proteins (RBPs) that are encoded in animal genomes it is important to identify their target RNAs. Although it is generally accepted that the binding specificity of an RBP is well described in terms of the nucleotide sequence of its binding sites, other factors such as the structural accessibility of binding sites or their clustering, to enable binding of RBP multimers, are also believed to play a role. Here we focus on GLD-1, a translational regulator of Caenorhabditis elegans, whose binding specificity and targets have been studied with a variety of methods such as CLIP (cross-linking and immunoprecipitation), RIP-Chip (microarray measurement of RNAs associated with an immunoprecipitated protein), profiling of polysome-associated mRNAs and biophysical determination of binding affinities of GLD-1 for short nucleotide sequences. We show that a simple biophysical model explains the binding of GLD-1 to mRNA targets to a large extent, and that taking into account the accessibility of putative target sites significantly improves the prediction of GLD-1 binding, particularly due to a more accurate prediction of binding in transcript coding regions. Relating GLD-1 binding to translational repression and stabilization of its target transcripts we find that binding sites along the entire transcripts contribute to functional responses, and that CDS-located sites contribute most to translational repression. Finally, biophysical measurements of GLD-1 affinity for a small number of oligonucleotides appear to allow an accurate reconstruction of the sequence specificity of the protein. This approach can be applied to uncover the specificity and function of other RBPs.

Timing control in regulatory networks by multisite protein modifications.

Computational and experimental studies have yielded quantitative insights into the role for multisite phosphorylation, and other protein modifications, in cell function. This work has emphasized the creation of thresholds and switches for cellular decisions. To date, the dynamics of phosphorylation events have been disregarded yet could be equally relevant for cell function. Here, we discuss theoretical predictions about the kinetic functions of multisite phosphorylation in regulatory networks and how these predictions relate to experimental findings. Using DNA replication as an example, we demonstrate that multisite phosphorylations can support coherent origin firing and robustness against rereplication. We suggest that multisite protein modifications provide a molecular mechanism to robustly time cellular events in the cell cycle, the circadian clock and signal transduction.

Mathematical modelling of DNA replication reveals a trade-off between coherence of origin activation and robustness against rereplication.

Eukaryotic genomes are duplicated from multiple replication origins exactly once per cell cycle. In Saccharomyces cerevisiae, a complex molecular network has been identified that governs the assembly of the replication machinery. Here we develop a mathematical model that links the dynamics of this network to its performance in terms of rate and coherence of origin activation events, number of activated origins, the resulting distribution of replicon sizes and robustness against DNA rereplication. To parameterize the model, we use measured protein expression data and systematically generate kinetic parameter sets by optimizing the coherence of origin firing. While randomly parameterized networks yield unrealistically slow kinetics of replication initiation, networks with optimized parameters account for the experimentally observed distribution of origin firing times. Efficient inhibition of DNA rereplication emerges as a constraint that limits the rate at which replication can be initiated. In addition to the separation between origin licensing and firing, a time delay between the activation of S phase cyclin-dependent kinase (S-Cdk) and the initiation of DNA replication is required for preventing rereplication. Our analysis suggests that distributive multisite phosphorylation of the S-Cdk targets Sld2 and Sld3 can generate both a robust time delay and contribute to switch-like, coherent activation of replication origins. The proposed catalytic function of the complex formed by Dpb11, Sld3 and Sld2 strongly enhances coherence and robustness of origin firing. The model rationalizes how experimentally observed inefficient replication from fewer origins is caused by premature activation of S-Cdk, while premature activity of the S-Cdk targets Sld2 and Sld3 results in DNA rereplication. Thus the model demonstrates how kinetic deregulation of the molecular network governing DNA replication may result in genomic instability.