Dicer ablation affects antibody diversity and cell survival in the B lymphocyte lineage.

To explore the role of Dicer-dependent control mechanisms in B lymphocyte development, we ablated this enzyme in early B cell progenitors. This resulted in a developmental block at the pro- to pre-B cell transition. Gene-expression profiling revealed a miR-17 approximately 92 signature in the 3'UTRs of genes upregulated in Dicer-deficient pro-B cells; a top miR-17 approximately 92 target, the proapoptotic molecule Bim, was highly upregulated. Accordingly, B cell development could be partially rescued by ablation of Bim or transgenic expression of the prosurvival protein Bcl-2. This allowed us to assess the impact of Dicer deficiency on the V(D)J recombination program in developing B cells. We found intact Ig gene rearrangements in immunoglobulin heavy (IgH) and kappa chain loci, but increased sterile transcription and usage of D(H) elements of the DSP family in IgH, and increased N sequence addition in Igkappa due to deregulated transcription of the terminal deoxynucleotidyl transferase gene.

Perfect seed pairing is not a generally reliable predictor for miRNA-target interactions.

We use Caenorhabditis elegans to test proposed general rules for microRNA (miRNA)-target interactions. We show that G.U base pairing is tolerated in the 'seed' region of the lsy-6 miRNA interaction with its in vivo target cog-1, and that 6- to 8-base-pair perfect seed pairing is not a generally reliable predictor for an interaction of lsy-6 with a 3' untranslated region (UTR). Rather, lsy-6 can functionally interact with its target site only in specific 3' UTR contexts. Our findings illustrate the difficulty of establishing generalizable rules of miRNA-target interactions.

Novel structural determinants in human SECIS elements modulate the translational recoding of UGA as selenocysteine.

The selenocysteine insertion sequence (SECIS) element directs the translational recoding of UGA as selenocysteine. In eukaryotes, the SECIS is located downstream of the UGA codon in the 3'-UTR of the selenoprotein mRNA. Despite poor sequence conservation, all SECIS elements form a similar stem-loop structure containing a putative kink-turn motif. We functionally characterized the 26 SECIS elements encoded in the human genome. Surprisingly, the SECIS elements displayed a wide range of UGA recoding activities, spanning several 1000-fold in vivo and several 100-fold in vitro. The difference in activity between a representative strong and weak SECIS element was not explained by differential binding affinity of SECIS binding Protein 2, a limiting factor for selenocysteine incorporation. Using chimeric SECIS molecules, we identified the internal loop and helix 2, which flank the kink-turn motif, as critical determinants of UGA recoding activity. The simultaneous presence of a GC base pair in helix 2 and a U in the 5'-side of the internal loop was a statistically significant predictor of weak recoding activity. Thus, the SECIS contains intrinsic information that modulates selenocysteine incorporation efficiency.

A physical and functional link between splicing factors promotes pre-mRNA 3' end processing.

Polypyrimidine tract-binding protein (PTB) is a splicing regulator that also plays a positive role in pre-mRNA 3' end processing when bound upstream of the polyadenylation signal (pA signal). Here, we address the mechanism of PTB stimulatory function in mRNA 3' end formation. We identify PTB as the protein factor whose binding to the human beta-globin (HBB) 3' UTR is abrogated by a 3' end processing-inactivating mutation. We show that PTB promotes both in vitro 3' end cleavage and polyadenylation and recruits directly the splicing factor hnRNP H to G-rich sequences associated with several pA signals. Increased binding of hnRNP H results in stimulation of polyadenylation through a direct interaction with poly(A) polymerase. Therefore, our results provide evidence of a concerted regulation of pA signal recognition by splicing factors bound to auxiliary polyadenylation sequence elements.

A short motif in Drosophila SECIS Binding Protein 2 provides differential binding affinity to SECIS RNA hairpins.

Selenoproteins contain the amino acid selenocysteine which is encoded by a UGA Sec codon. Recoding UGA Sec requires a complex mechanism, comprising the cis-acting SECIS RNA hairpin in the 3'UTR of selenoprotein mRNAs, and trans-acting factors. Among these, the SECIS Binding Protein 2 (SBP2) is central to the mechanism. SBP2 has been so far functionally characterized only in rats and humans. In this work, we report the characterization of the Drosophila melanogaster SBP2 (dSBP2). Despite its shorter length, it retained the same selenoprotein synthesis-promoting capabilities as the mammalian counterpart. However, a major difference resides in the SECIS recognition pattern: while human SBP2 (hSBP2) binds the distinct form 1 and 2 SECIS RNAs with similar affinities, dSBP2 exhibits high affinity toward form 2 only. In addition, we report the identification of a K (lysine)-rich domain in all SBP2s, essential for SECIS and 60S ribosomal subunit binding, differing from the well-characterized L7Ae RNA-binding domain. Swapping only five amino acids between dSBP2 and hSBP2 in the K-rich domain conferred reversed SECIS-binding properties to the proteins, thus unveiling an important sequence for form 1 binding.

A highly conserved cis-regulatory motif directs differential gonadal synexpression of Dmrt1 transcripts during gonad development.

Differential gene expression largely accounts for the coordinated manifestation of the genetic programme underlying embryonic development and cell differentiation. The 3' untranslated region (3'-UTR) of eukaryotic genes can contain motifs involved in regulation of gene expression at the post-transcriptional level. In the 3'-UTR of dmrt1, a key gene that functions in gonad development and differentiation, an 11-bp protein-binding motif was identified that mediates gonad-specific mRNA localization during embryonic and larval development of fish. Mutations that disrupt the 11-bp motif leading to in vitro protein-binding loss and selective transcript stabilization failure indicate a role for this motif in RNA stabilization through protein binding. The sequence motif was found to be conserved in most of the dmrt1 homologous genes from flies to humans suggesting a widespread conservation of this specific mechanism.

Allegro: analyzing expression and sequence in concert to discover regulatory programs.

A major goal of system biology is the characterization of transcription factors and microRNAs (miRNAs) and the transcriptional programs they regulate. We present Allegro, a method for de-novo discovery of cis-regulatory transcriptional programs through joint analysis of genome-wide expression data and promoter or 3' UTR sequences. The algorithm uses a novel log-likelihood-based, non-parametric model to describe the expression pattern shared by a group of co-regulated genes. We show that Allegro is more accurate and sensitive than existing techniques, and can simultaneously analyze multiple expression datasets with more than 100 conditions. We apply Allegro on datasets from several species and report on the transcriptional modules it uncovers. Our analysis reveals a novel motif over-represented in the promoters of genes highly expressed in murine oocytes, and several new motifs related to fly development. Finally, using stem-cell expression profiles, we identify three miRNA families with pivotal roles in human embryogenesis.

Identification of a signature motif in target mRNAs of RNA-binding protein AUF1.

The ubiquitous RNA-binding protein AUF1 promotes the degradation of some target mRNAs, but increases the stability and translation of other targets. Here, we isolated AUF1-associated mRNAs by immunoprecipitation of (AUF1-RNA) ribonucleoprotein (RNP) complexes from HeLa cells, identified them using microarrays, and used them to elucidate a signature motif shared among AUF1 target transcripts. The predicted AUF1 motif (29-39 nucleotides) contained 79% As and Us, consistent with the AU-rich sequences of reported AUF1 targets. Importantly, 10 out of 15 previously reported AUF1 target mRNAs contained the AUF1 motif. The predicted interactions between AUF1 and target mRNAs were recapitulated in vitro using biotinylated RNAs. Interestingly, further validation of predicted AUF1 target transcripts revealed that AUF1 associates with both the pre-mRNA and the mature mRNA forms. The consequences of AUF1 binding to 10 predicted target mRNAs were tested by silencing AUF1, which elevated the steady-state levels of only four mRNAs, and by overexpressing AUF1, which also lowered the levels of only four mRNAs. In total, we have identified a signature motif in AUF1 target mRNAs, have found that AUF1 also associates with the corresponding pre-mRNAs, and have discovered that altering AUF1 levels alone only modifies the levels of subsets of target mRNAs.

Mouse period 2 mRNA circadian oscillation is modulated by PTB-mediated rhythmic mRNA degradation.

Circadian mRNA oscillations are the main feature of core clock genes. Among them, period 2 is a key component in negative-feedback regulation, showing robust diurnal oscillations. Moreover, period 2 has been found to have a physiological role in the cell cycle or the tumor suppression. The present study reports that 3'-untranslated region (UTR)-dependent mRNA decay is involved in the regulation of circadian oscillation of period 2 mRNA. Within the mper2 3'UTR, both the CU-rich region and polypyrimidine tract-binding protein (PTB) are more responsible for mRNA stability and degradation kinetics than are other factors. Depletion of PTB with RNAi results in mper2 mRNA stabilization. During the circadian oscillations of mper2, cytoplasmic PTB showed a reciprocal expression profile compared with mper2 mRNA and its peak amplitude was increased when PTB was depleted. This report on the regulation of mper2 proposes that post-transcriptional mRNA decay mediated by PTB is a fine-tuned regulatory mechanism that includes dampening-down effects during circadian mRNA oscillations.

Long 3'-UTRs target wild-type mRNAs for nonsense-mediated mRNA decay in Saccharomyces cerevisiae.

The nonsense-mediated mRNA decay (NMD) pathway, present in most eukaryotic cells, is a specialized pathway that leads to the recognition and rapid degradation of mRNAs with premature termination codons and, importantly, some wild-type mRNAs. Earlier studies demonstrated that aberrant mRNAs with artificially extended 3'-untranslated regions (3'-UTRs) are degraded by NMD. However, the extent to which wild-type mRNAs with long 3'-UTRs are degraded by NMD is not known. We used a global approach to identify wild-type mRNAs in Saccharomyces cerevisiae that have longer than expected 3'-UTRs, and of these mRNAs tested, 91% were degraded by NMD. We demonstrate for the first time that replacement of the natural, long 3'-UTR from wild-type PGA1 mRNA, which encodes a protein that is important for cell wall biosynthesis, with a short 3'-UTR renders it immune to NMD. The natural PGA1 3'-UTR is sufficient to target a NMD insensitive mRNA for decay by the NMD pathway. Finally, we show that nmd mutants are sensitive to Calcofluor White, which suggests that the regulation of PGA1 and other cell wall biosynthesis proteins by NMD is physiologically significant.

Alternative polyadenylation variants of the RNA binding protein, HuR: abundance, role of AU-rich elements and auto-Regulation.

The RNA-binding protein, HuR, is involved in the stabilization of AU-rich element-containing mRNAs with products that are involved in cell-cycle progression, cell differentiation and inflammation. We show that there are multiple polyadenylation variants of HuR mRNA that differ in their abundance, using both bioinformatics and experimental approaches. A polyadenylation variant with distal poly(A) signal is a rare transcript that harbors functional AU-rich elements (ARE) in the 3'UTR. A minimal 60-nt region, but not a mutant form, fused to reporter-3'UTR constructs was able to downregulate the reporter activity. The most predominant and alternatively polyadenylated mature transcript does not contain the ARE. HuR itself binds HuR mRNA, and upregulated the activity of reporter from constructs fused with ARE-isoform and the HuR ARE. Wild-type tristetraprolin (TTP), but not the zinc finger mutant TTP, competes for HuR binding and upregulation of HuR mRNA. The study shows that the HuR gene codes for several polyadenylation variants differentially regulated by AU-rich elements, and demonstrates an auto-regulatory role of HuR.

Labeled microRNA pull-down assay system: an experimental approach for high-throughput identification of microRNA-target mRNAs.

We developed a simple, direct and cost-effective approach to search for the most likely target genes of a known microRNA (miRNA) in vitro. We term this method 'labeled miRNA pull-down (LAMP)' assay system. Briefly, the pre-miRNA is labeled with digoxigenin (DIG), mixed with cell extracts and immunoprecipitated by anti-DIG antiserum. When the DIG-labeled miRNA and bound mRNA complex are obtained, the total cDNAs are then subcloned and sequenced, or RT-PCR-amplified, to search for the putative target genes of a known miRNA. After successfully identifying the known target genes of Caenorhabditis elegans miRNAs lin-4 and let-7 and zebrafish let-7, we applied LAMP to find the unknown target gene of zebrafish miR-1, which resulted in the identification of hand2. We then confirmed hand2 as a novel target gene of miR-1 by whole-mount in situ hybridization and luciferase reporter gene assay. We further validated this target gene by microarray analysis, and the results showed that hand2 is the top-scoring among 302 predicted putative target genes. We concluded that LAMP is an experimental approach for high-throughput identification of the target gene of known miRNAs from both C. elegans and zebrafish, yielding fewer false positive results than those produced by using only the bioinformatics approach.

The proapoptotic dp5 gene is a direct target of the MLK-JNK-c-Jun pathway in sympathetic neurons.

The death of sympathetic neurons after nerve growth factor (NGF) withdrawal requires de novo gene expression. Dp5 was one of the first NGF withdrawal-induced genes to be identified and it encodes a proapoptotic BH3-only member of the Bcl-2 family. To study how dp5 transcription is regulated by NGF withdrawal we cloned the regulatory regions of the rat dp5 gene and constructed a series of dp5-luciferase reporter plasmids. In microinjection experiments with sympathetic neurons we found that three regions of dp5 contribute to its induction after NGF withdrawal: the promoter, a conserved region in the single intron, and sequences in the 3' untranslated region of the dp5 mRNA. A construct containing all three regions is efficiently activated by NGF withdrawal and, like the endogenous dp5, its induction requires mixed-lineage kinase (MLK) and c-Jun N-terminal kinase (JNK) activity. JNKs phosphorylate the AP-1 transcription factor c-Jun, and thereby increase its activity. We identified a conserved ATF site in the dp5 promoter that binds c-Jun and ATF2, which is critical for dp5 promoter induction after NGF withdrawal. These results suggest that part of the mechanism by which the MLK-JNK-c-Jun pathway promotes neuronal apoptosis is by activating the transcription of the dp5 gene.

Estradiol downregulates miR-21 expression and increases miR-21 target gene expression in MCF-7 breast cancer cells.

Select changes in microRNA (miRNA) expression correlate with estrogen receptor alpha (ER alpha) expression in breast tumors. miR-21 is higher in ER alpha positive than negative tumors, but no one has examined how estradiol (E(2)) regulates miR-21 in breast cancer cells. Here we report that E(2) inhibits miR-21 expression in MCF-7 human breast cancer cells. The E(2)-induced reduction in miR-21 was inhibited by 4-hydroxytamoxifen (4-OHT), ICI 182 780 (Faslodex), and siRNA ER alpha indicating that the suppression is ER alpha-mediated. ER alpha and ER beta agonists PPT and DPN inhibited and 4-OHT increased miR-21 expression. E(2) increased luciferase activity from reporters containing the miR-21 recognition elements from the 3'-UTRs of miR-21 target genes, corroborating that E(2) represses miR-21 expression resulting in a loss of target gene suppression. The E(2)-mediated decrease in miR-21 correlated with increased protein expression of endogenous miR-21-targets Pdcd4, PTEN and Bcl-2. siRNA knockdown of ER alpha blocked the E(2)-induced increase in Pdcd4, PTEN and Bcl-2. Transfection of MCF-7 cells with antisense (AS) to miR-21 mimicked the E(2)-induced increase in Pdcd4, PTEN and Bcl-2. These results are the first to demonstrate that E(2) represses the expression of an oncogenic miRNA, miR-21, by activating estrogen receptor in MCF-7 cells.

A common variant in the 3'UTR of the GRIK4 glutamate receptor gene affects transcript abundance and protects against bipolar disorder.

Underactivity of the glutamatergic system is an attractive model for the pathophysiology of several major mental illnesses. We previously described a chromosome abnormality disrupting the kainate class ionotropic glutamate receptor gene, GRIK4/KA1, in an individual with schizophrenia and learning disability (mental retardation). We also demonstrated in a case-control study that two physically separated haplotypes within this gene were significantly associated with increased risk of schizophrenia and decreased risk of bipolar disorder, respectively. The latter protective haplotype was located at the 3' end of the gene. We now report the identification from carriers of the protective haplotype of a deletion variant within the 3' untranslated region of the gene. The deletion allele also was found to be negatively associated with bipolar disorder in both initial (P = 0.00000019) and replication (P = 0.0107) case-control studies. Expression studies indicated that deletion-carrying mRNA transcripts were relatively more abundant. We postulate that this may be a direct consequence of the differences in the RNA secondary structures predicted for the insertion and deletion alleles. These data suggest a mechanism whereby the genetic protective effect is mediated through increased kainate receptor expression.

Prediction of human miRNAs using tissue-selective motifs in 3' UTRs.

MicroRNAs (miRNAs) play an important role in posttranscriptional regulation of genes. We developed a method to predict human miRNAs without requiring cross-species conservation. We first identified lowly/moderately expressed tissue-selective genes using EST data and then identified overrepresented motifs of seven nucleotides in the 3' UTRs of these genes. Using these motifs as potential target sites of miRNAs, we recovered more than two-thirds of the known human miRNAs. We then used those motifs that did not match any known human miRNA seed region to infer novel miRNAs. We predicted 36 new human miRNA genes with 44 mature forms and 4 novel alternative mature forms of 2 known miRNA genes when a stringent criterion was used and many more novel miRNAs when a less stringent criterion was used. We tested the expression of 11 predicted miRNAs in three human cell lines and found 5 of them expressed in all three cell lines and 1 expressed in one cell line. We selected 2 of them, P-2 and P-27-5p, to do functional validation, using their mimics and inhibitors and using both luciferase assay and Western blotting. These experiments provided strong evidence that both P-2 and P-27-5p are novel miRNAs and that CREB3L3, which encodes cAMP-responsive element binding protein 3-like 3, is a target gene of P-2, whereas LAMB3, which encodes laminin beta3, is a target gene of P-27-5p.

Molecular architecture of a miRNA-regulated 3' UTR.

Animal genomes contain hundreds of microRNAs (miRNAs), small regulatory RNAs that control gene expression by binding to complementary sites in target mRNAs. Some rules that govern miRNA/target interaction have been elucidated but their general applicability awaits further experimentation on a case-by-case basis. We use here an assay system in transgenic nematodes to analyze the interaction of the Caenorhabditis elegans lsy-6 miRNA with 3' UTR sequences. In contrast to many previously described assay systems used to analyze miRNA/target interactions, our assay system operates within the cellular context in which lsy-6 normally functions, a single neuron in the nervous system of C. elegans. Through extensive mutational analysis, we define features in the known and experimentally validated target of lsy-6, the 3' UTR of the cog-1 homeobox gene, that are required for a functional miRNA/target interaction. We describe that both in the context of the cog-1 3' UTR and in the context of heterologous 3' UTRs, one or more seed matches are not a reliable predictor for a functional miRNA/target interaction. We rather find that two nonsequence specific contextual features beyond miRNA target sites are critical determinants of miRNA-mediated 3' UTR regulation. The contextual features reside 3' of lsy-6 binding sites in the 3' UTR and act in a combinatorial manner; mutation of each results in limited defects in 3' UTR regulation, but a combinatorial deletion results in complete loss of 3' UTR regulation. Together with two lsy-6 sites, these two contextual features are capable of imparting regulation on a heterologous 3' UTR. Moreover, the contextual features need to be present in a specific configuration relative to miRNA binding sites and could either represent protein binding sites or provide an appropriate structural context. We conclude that a given target site resides in a 3' UTR context that evolved beyond target site complementarity to support regulation by a specific miRNA. The large number of 3' UTRs that we analyzed in this study will also be useful to computational biologists in designing the next generation of miRNA/target prediction algorithms.

SECISaln, a web-based tool for the creation of structure-based alignments of eukaryotic SECIS elements.

SUMMARY: Selenoproteins contain the 21st amino acid selenocysteine which is encoded by an inframe UGA codon, usually read as a stop. In eukaryotes, its co-translational recoding requires the presence of an RNA stem-loop structure, the SECIS element in the 3 untranslated region of (UTR) selenoprotein mRNAs. Despite little sequence conservation, SECIS elements share the same overall secondary structure. Until recently, the lack of a significantly high number of selenoprotein mRNA sequences hampered the identification of other potential sequence conservation. In this work, the web-based tool SECISaln provides for the first time an extensive structure-based sequence alignment of SECIS elements resulting from the well-defined secondary structure of the SECIS RNA and the increased size of the eukaryotic selenoproteome. We have used SECISaln to improve our knowledge of SECIS secondary structure and to discover novel, conserved nucleotide positions and we believe it will be a useful tool for the selenoprotein and RNA scientific communities. AVAILABILITY: SECISaln is freely available as a web-based tool at http://genome.crg.es/software/secisaln/.

Genome level analysis of rice mRNA 3'-end processing signals and alternative polyadenylation.

The position of a poly(A) site of eukaryotic mRNA is determined by sequence signals in pre-mRNA and a group of polyadenylation factors. To reveal rice poly(A) signals at a genome level, we constructed a dataset of 55 742 authenticated poly(A) sites and characterized the poly(A) signals. This resulted in identifying the typical tripartite cis-elements, including FUE, NUE and CE, as previously observed in Arabidopsis. The average size of the 3'-UTR was 289 nucleotides. When mapped to the genome, however, 15% of these poly(A) sites were found to be located in the currently annotated intergenic regions. Moreover, an extensive alternative polyadenylation profile was evident where 50% of the genes analyzed had more than one unique poly(A) site (excluding microheterogeneity sites), and 13% had four or more poly(A) sites. About 4% of the analyzed genes possessed alternative poly(A) sites at their introns, 5'-UTRs, or protein coding regions. The authenticity of these alternative poly(A) sites was partially confirmed using MPSS data. Analysis of nucleotide profile and signal patterns indicated that there may be a different set of poly(A) signals for those poly(A) sites found in the coding regions. Based on the features of rice poly(A) signals, an updated algorithm termed PASS-Rice was designed to predict poly(A) sites.

UTRome.org: a platform for 3'UTR biology in C. elegans.

Three-prime untranslated regions (3'UTRs) are widely recognized as important post-transcriptional regulatory regions of mRNAs. RNA-binding proteins and small non-coding RNAs such as microRNAs (miRNAs) bind to functional elements within 3'UTRs to influence mRNA stability, translation and localization. These interactions play many important roles in development, metabolism and disease. However, even in the most well-annotated metazoan genomes, 3'UTRs and their functional elements are not well defined. Comprehensive and accurate genome-wide annotation of 3'UTRs and their functional elements is thus critical. We have developed an open-access database, available at http://www.UTRome.org, to provide a rich and comprehensive resource for 3'UTR biology in the well-characterized, experimentally tractable model system Caenorhabditis elegans. UTRome.org combines data from public repositories and a large-scale effort we are undertaking to characterize 3'UTRs and their functional elements in C. elegans, including 3'UTR sequences, graphical displays, predicted and validated functional elements, secondary structure predictions and detailed data from our cloning pipeline. UTRome.org will grow substantially over time to encompass individual 3'UTR isoforms for the majority of genes, new and revised functional elements, and in vivo data on 3'UTR function as they become available. The UTRome database thus represents a powerful tool to better understand the biology of 3'UTRs.