Grzegorz Kudla

Synonymous but not the same: the causes and consequences of codon bias

Despite their name, synonymous mutations have significant consequences for cellular processes in all taxa. As a result, an understanding of codon bias is central to fields as diverse as molecular evolution and biotechnology. Although recent advances in sequencing and synthetic biology have helped to resolve longstanding questions about codon bias, they have also uncovered striking patterns that suggest new hypotheses about protein synthesis. Ongoing work to quantify the dynamics of initiation and elongation is as important for understanding natural synonymous variation as it is for designing transgenes in applied contexts.

Mapping the Human miRNA Interactome by CLASH Reveals Frequent Noncanonical Binding

Summary MicroRNAs (miRNAs) play key roles in gene regulation, but reliable bioinformatic or experimental identification of their targets remains difficult. To provide an unbiased view of human miRNA targets, we developed a technique for ligation and sequencing of miRNA-target RNA duplexes associated with human AGO1. Here, we report data sets of more than 18,000 high-confidence miRNA-mRNA interactions. The binding of most miRNAs includes the 5′ seed region, but around 60% of seed interactions are noncanonical, containing bulged or mismatched nucleotides. Moreover, seed interactions are generally accompanied by specific, nonseed base pairing. 18% of miRNA-mRNA interactions involve the miRNA 3′ end, with little evidence for 5′ contacts, and some of these were functionally validated. Analyses of miRNA:mRNA base pairing showed that miRNA species systematically differ in their target RNA interactions, and strongly overrepresented motifs were found in the interaction sites of several miRNAs. We speculate that these affect the response of RISC to miRNA-target binding.

Rate-Limiting Steps in Yeast Protein Translation

Summary Deep sequencing now provides detailed snapshots of ribosome occupancy on mRNAs. We leverage these data to parameterize a computational model of translation, keeping track of every ribosome, tRNA, and mRNA molecule in a yeast cell. We determine the parameter regimes in which fast initiation or high codon bias in a transgene increases protein yield and infer the initiation rates of endogenous Saccharomyces cerevisiae genes, which vary by several orders of magnitude and correlate with 5′ mRNA folding energies. Our model recapitulates the previously reported 5′-to-3′ ramp of decreasing ribosome densities, although our analysis shows that this ramp is caused by rapid initiation of short genes rather than slow codons at the start of transcripts. We conclude that protein production in healthy yeast cells is typically limited by the availability of free ribosomes, whereas protein production under periods of stress can sometimes be rescued by reducing initiation or elongation rates.

Identification of protein binding sites on U3 snoRNA and pre-rRNA by UV cross-linking and high-throughput analysis of cDNAs

The U3 small nucleolar ribonucleoprotein (snoRNP) plays an essential role in ribosome biogenesis but, like many RNA–protein complexes, its architecture is poorly understood. To address this problem, binding sites for the snoRNP proteins Nop1, Nop56, Nop58, and Rrp9 were mapped by UV cross-linking and analysis of cDNAs. Cross-linked protein–RNA complexes were purified under highly-denaturing conditions, ensuring that only direct interactions were detected. Recovered RNA fragments were amplified after linker ligation and cDNA synthesis. Cross-linking was successfully performed either in vitro on purified complexes or in vivo in living cells. Cross-linking sites were precisely mapped either by Sanger sequencing of multiple cloned fragments or direct, high-throughput Solexa sequencing. Analysis of RNAs associated with the snoRNP proteins revealed remarkably high signal-to-noise ratios and identified specific binding sites for each of these proteins on the U3 RNA. The results were consistent with previous data, demonstrating the reliability of the method, but also provided insights into the architecture of the U3 snoRNP. The snoRNP proteins were also cross-linked to pre-rRNA fragments, with preferential association at known sites of box C/D snoRNA function. This finding demonstrates that the snoRNP proteins directly contact the pre-rRNA substrate, suggesting roles in snoRNA recruitment. The techniques reported here should be widely applicable to analyses of RNA–protein interactions.

Cross-linking, ligation, and sequencing of hybrids reveals RNA-RNA interactions in yeast.

Many protein–protein and protein–nucleic acid interactions have been experimentally characterized, whereas RNA–RNA interactions have generally only been predicted computationally. Here, we describe a high-throughput method to identify intramolecular and intermolecular RNA–RNA interactions experimentally by cross-linking, ligation, and sequencing of hybrids (CLASH). As validation, we identified 39 known target sites for box C/D modification-guide small nucleolar RNAs (snoRNAs) on the yeast pre-rRNA. Novel snoRNA-rRNA hybrids were recovered between snR4-5S and U14-25S. These are supported by native electrophoresis and consistent with previously unexplained data. The U3 snoRNA was found to be associated with sequences close to the 3′ side of the central pseudoknot in 18S rRNA, supporting a role in formation of this structure. Applying CLASH to the yeast U2 spliceosomal snRNA led to a revised predicted secondary structure, featuring alternative folding of the 3′ domain and long-range contacts between the 3′ and 5′ domains. CLASH should allow transcriptome-wide analyses of RNA–RNA interactions in many organisms.

Evidence for the biogenesis of more than 1,000 novel human microRNAs

BackgroundMicroRNAs (miRNAs) are established regulators of development, cell identity and disease. Although nearly two thousand human miRNA genes are known and new ones are continuously discovered, no attempt has been made to gauge the total miRNA content of the human genome.ResultsEmploying an innovative computational method on massively pooled small RNA sequencing data, we report 2,469 novel human miRNA candidates of which 1,098 are validated by in-house and published experiments. Almost 300 candidates are robustly expressed in a neuronal cell system and are regulated during differentiation or when biogenesis factors Dicer, Drosha, DGCR8 or Ago2 are silenced. To improve expression profiling, we devised a quantitative miRNA capture system. In a kidney cell system, 400 candidates interact with DGCR8 at transcript positions that suggest miRNA hairpin recognition, and 1,000 of the new miRNA candidates interact with Ago1 or Ago2, indicating that they are directly bound by miRNA effector proteins. From kidney cell CLASH experiments, in which miRNA-target pairs are ligated and sequenced, we observe hundreds of interactions between novel miRNAs and mRNA targets. The novel miRNA candidates are specifically but lowly expressed, raising the possibility that not all may be functional. Interestingly, the majority are evolutionarily young and overrepresented in the human brain.ConclusionsIn summary, we present evidence that the complement of human miRNA genes is substantially larger than anticipated, and that more are likely to be discovered in the future as more tissues and experimental conditions are sequenced to greater depth.

The nuclear RNA polymerase II surveillance system targets polymerase III transcripts: Targets for nuclear surveillance

A key question in nuclear RNA surveillance is how target RNAs are recognized. To address this, we identified in vivo binding sites for nuclear RNA surveillance factors, Nrd1, Nab3 and the Trf4/5–Air1/2–Mtr4 polyadenylation (TRAMP) complex poly(A) polymerase Trf4, by UV crosslinking. Hit clusters were reproducibly found over known binding sites on small nucleolar RNAs (snoRNAs), pre‐mRNAs and cryptic, unstable non‐protein‐coding RNAs (ncRNAs) (‘CUTs’), along with ∼642 predicted long anti‐sense ncRNAs (asRNAs), ∼178 intergenic ncRNAs and, surprisingly, ∼1384 mRNAs. Five putative asRNAs tested were confirmed to exist and were stabilized by loss of Nrd1, Nab3 or Trf4. Mapping of micro‐deletions and substitutions allowed clear definition of preferred, in vivo Nab3 and Nrd1 binding sites. Nrd1 and Nab3 were believed to be Pol II specific but, unexpectedly, bound many oligoadenylated Pol III transcripts, predominately pre‐tRNAs. Depletion of Nrd1 or Nab3 stabilized tested Pol III transcripts and their oligoadenylation was dependent on Nrd1–Nab3 and TRAMP. Surveillance targets were enriched for non‐encoded A‐rich tails. These were generally very short (1–5 nt), potentially explaining why adenylation destabilizes these RNAs while stabilizing mRNAs with long poly(A) tails.

Transcriptome-wide Analysis of Exosome Targets

Summary The exosome plays major roles in RNA processing and surveillance but the in vivo target range and substrate acquisition mechanisms remain unclear. Here we apply in vivo RNA crosslinking (CRAC) to the nucleases (Rrp44, Rrp6), two structural subunits (Rrp41, Csl4) and a cofactor (Trf4) of the yeast exosome. Analysis of wild-type Rrp44 and catalytic mutants showed that both the CUT and SUT classes of non-coding RNA, snoRNAs and, most prominently, pre-tRNAs and other Pol III transcripts are targeted for oligoadenylation and exosome degradation. Unspliced pre-mRNAs were also identified as targets for Rrp44 and Rrp6. CRAC performed using cleavable proteins (split-CRAC) revealed that Rrp44 endonuclease and exonuclease activities cooperate on most substrates. Mapping oligoadenylated reads suggests that the endonuclease activity may release stalled exosome substrates. Rrp6 was preferentially associated with structured targets, which frequently did not associate with the core exosome indicating that substrates follow multiple pathways to the nucleases.

Murine cytomegalovirus encodes a miR-27 inhibitor disguised as a target

Individual microRNAs (miRNAs) are rapidly down-regulated during conditions of cellular activation and infection, but factors mediating miRNA turnover are poorly understood. Infection of mouse cells with murine cytomegalovirus (MCMV) induces the rapid down-regulation of an antiviral cellular miRNA, miR-27. Here, we identify a transcript produced by MCMV that binds to miR-27 and mediates its degradation. UV-crosslinking and high-throughput sequencing [CRAC (UV-crosslinking and analysis of cDNA)] identified MCMV RNA segments associated with the miRNA-binding protein Argonaute 2 (Ago2). A cluster of hits mapped to a predicted miR-27-binding site in the 3′UTR of the previously uncharacterized ORF, m169. The expression kinetics of the m169 transcript correlated with degradation of miR-27 during infection, and m169 expression inhibited miR-27 functional activity in a reporter assay. siRNA knockdown of m169 demonstrated its requirement for miR-27 degradation following infection and did not affect other host miRNAs. Substitution of the miR-27-binding site in m169 to create complementarity to a different cellular miRNA, miR-24, resulted in down-regulation of only miR-24 following infection. The m169 transcript is cytoplasmic, capped, polyadenylated, and interacts with miRNA-27 through seed pairing: characteristic features of the normal messenger RNA (mRNA) targets of miRNAs. This virus–host interaction reveals a mode of miRNA regulation in which a mRNA directs the degradation of a miRNA. We speculate that RNA-mediated miRNA degradation could be a more general viral strategy for manipulating host cells.

PAR-CLIP data indicate that Nrd1-Nab3-dependent transcription termination regulates expression of hundreds of protein coding genes in yeast

BackgroundNrd1 and Nab3 are essential sequence-specific yeast RNA binding proteins that function as a heterodimer in the processing and degradation of diverse classes of RNAs. These proteins also regulate several mRNA coding genes; however, it remains unclear exactly what percentage of the mRNA component of the transcriptome these proteins control. To address this question, we used the pyCRAC software package developed in our laboratory to analyze CRAC and PAR-CLIP data for Nrd1-Nab3-RNA interactions.ResultsWe generated high-resolution maps of Nrd1-Nab3-RNA interactions, from which we have uncovered hundreds of new Nrd1-Nab3 mRNA targets, representing between 20 and 30% of protein-coding transcripts. Although Nrd1 and Nab3 showed a preference for binding near 5′ ends of relatively short transcripts, they bound transcripts throughout coding sequences and 3′ UTRs. Moreover, our data for Nrd1-Nab3 binding to 3′ UTRs was consistent with a role for these proteins in the termination of transcription. Our data also support a tight integration of Nrd1-Nab3 with the nutrient response pathway. Finally, we provide experimental evidence for some of our predictions, using northern blot and RT-PCR assays.ConclusionsCollectively, our data support the notion that Nrd1 and Nab3 function is tightly integrated with the nutrient response and indicate a role for these proteins in the regulation of many mRNA coding genes. Further, we provide evidence to support the hypothesis that Nrd1-Nab3 represents a failsafe termination mechanism in instances of readthrough transcription.