Early evolution of histone mRNA 3' end processing

Evolutionary dynamics of polyadenylation signals and their recognition strategies in protists

The poly(A) signal, together with auxiliary elements, directs cleavage of a pre-mRNA and thus determines the 3' end of the mature transcript. In many species, including humans, the poly(A) signal is an AAUAAA hexamer, but we recently found that the deeply branching eukaryote Giardia lamblia uses a distinct hexamer (AGURAA) and lacks any known auxiliary elements. Our discovery prompted us to explore the evolutionary dynamics of poly(A) signals and auxiliary elements in the eukaryotic kingdom. We use direct RNA sequencing to determine poly(A) signals for four protists within the Metamonada clade (which also contains G. lamblia) and two outgroup protists. These experiments reveal that the AAUAAA hexamer serves as the poly(A) signal in at least four different eukaryotic clades, indicating that it is likely the ancestral signal, whereas the unusual Giardia version is derived. We find that the use and relative strengths of auxiliary elements are also plastic; in fact, within Metamonada, species like G. lamblia make use of a previously unrecognized auxiliary element where nucleotides flanking the poly(A) signal itself specify genuine cleavage sites. Thus, despite the fundamental nature of pre-mRNA cleavage for the expression of all protein-coding genes, the motifs controlling this process are dynamic on evolutionary timescales, providing motivation for future biochemical and structural studies as well as new therapeutic angles to target eukaryotic pathogens.

Translation machinery: the basis of translational control

Messenger RNA (mRNA) translation consists of initiation, elongation, termination, and ribosome recycling, carried out by the translation machinery, primarily including tRNAs, ribosomes, and translation factors (TrFs). Translational regulators transduce signals of growth and development, as well as biotic and abiotic stresses, to the translation machinery, where global or selective translational control occurs to modulate mRNA translation efficiency (TrE). As the basis of translational control, the translation machinery directly determines the quality and quantity of newly synthesized peptides and, ultimately, the cellular adaption. Thus, regulating the availability of diverse machinery components is reviewed as the central strategy of translational control. We provide classical signaling pathways (e.g., integrated stress responses) and cellular behaviors (e.g., liquid-liquid phase separation) to exemplify this strategy within different physiological contexts, particularly during host-microbe interactions. With new technologies developed, further understanding this strategy will speed up translational medicine and translational agriculture.

LeishIF4E2 is a cap-binding protein that plays a role in Leishmania cell cycle progression

Leishmania encode six paralogs of the cap-binding protein eIF4E and five eIF4G candidates, forming unique complexes. Two cap-binding proteins, LeishIF4E1 and LeishIF4E2, do not bind any identified LeishIF4Gs, thus their roles are intriguing. Here, we combine structural prediction, proteomic analysis, and interaction assays to shed light on LeishIF4E2 function. A nonconserved C-terminal extension was identified through structure prediction and sequence alignment. m7 GTP-binding assays involving both recombinant and transgenic LeishIF4E2 with and without the C-terminal extension revealed that this extension functions as a regulatory gate, modulating the cap-binding activity of LeishIF4E2. The interactomes of the two LeishIF4E2 versions were investigated, highlighting the role of the C-terminal extension in binding to SLBP2. SLBP2 is known to interact with a stem-loop structure in the 3' UTRs of histone mRNAs. Consistent with the predicted inhibitory effect of SLBP2 on histone expression in Xenopus laevis, a hemizygous deletion mutant of LeishIF4E2, exhibited an upregulation of several histones. We therefore propose that LeishIF4E2 is involved in histone expression, possibly through its interaction between SLBP2 and LeishIF4E2, thus affecting cell cycle progression. In addition, cell synchronization showed that LeishIF4E2 expression decreased during the S-phase, when histones are known to be synthesized. Previous studies in T. brucei also highlighted an association between TbEIF4E2 and SLBP2, and further reported on an interaction between TbIF4E2 and S-phase-abundant mRNAs. Our results show that overexpression of LeishIF4E2 correlates with upregulation of cell cycle and chromosome maintenance proteins. Along with its effect on histone expression, we propose that LeishIF4E2 is involved in cell cycle progression.

m6A is required for resolving progenitor identity during planarian stem cell differentiation

Regeneration and tissue homeostasis require accurate production of missing cell lineages. Cell production is driven by changes to gene expression, which is shaped by multiple layers of regulation. Here, we find that the ubiquitous mRNA base-modification, m6A, is required for proper cell fate choice and cellular maturation in planarian stem cells (neoblasts). We mapped m6A-enriched regions in 7,600 planarian genes and found that perturbation of the m6A pathway resulted in progressive deterioration of tissues and death. Using single-cell RNA sequencing of >20,000 cells following perturbation of the m6A pathway, we identified an increase in expression of noncanonical histone variants, and that inhibition of the pathway resulted in accumulation of undifferentiated cells throughout the animal in an abnormal transcriptional state. Analysis of >1,000 planarian gene expression datasets revealed that the inhibition of the chromatin modifying complex NuRD had almost indistinguishable consequences, unraveling an unappreciated link between m6A and chromatin modifications. Our findings reveal that m6A is critical for planarian stem cell homeostasis and gene regulation in tissue maintenance and regeneration.

Novel classes and evolutionary turnover of histone H2B variants in the mammalian germline

Histones and their posttranslational modifications facilitate diverse chromatin functions in eukaryotes. Core histones (H2A, H2B, H3, and H4) package genomes after DNA replication. In contrast, variant histones promote specialized chromatin functions, including DNA repair, genome stability, and epigenetic inheritance. Previous studies have identified only a few H2B variants in animals; their roles and evolutionary origins remain largely unknown. Here, using phylogenomic analyses, we reveal the presence of five H2B variants broadly present in mammalian genomes. Three of these variants have been previously described: H2B.1, H2B.L (also called subH2B), and H2B.W. In addition, we identify and describe two new variants: H2B.K and H2B.N. Four of these variants originated in mammals, whereas H2B.K arose prior to the last common ancestor of bony vertebrates. We find that though H2B variants are subject to high gene turnover, most are broadly retained in mammals, including humans. Despite an overall signature of purifying selection, H2B variants evolve more rapidly than core H2B with considerable divergence in sequence and length. All five H2B variants are expressed in the germline. H2B.K and H2B.N are predominantly expressed in oocytes, an atypical expression site for mammalian histone variants. Our findings suggest that H2B variants likely encode potentially redundant but vital functions via unusual chromatin packaging or nonchromatin functions in mammalian germline cells. Our discovery of novel histone variants highlights the advantages of comprehensive phylogenomic analyses and provides unique opportunities to study how innovations in chromatin function evolve.

Genome-Wide Characterization and Analysis of Expression of the Histone Gene Family in Razor Clam, Sinonovacula constricta

The Chinese razor clam (Sinonovacula constricta), a bivalve species widely distributed in estuaries and mudflats, is often exposed to extreme environmental and microbial stresses. Histones are fundamental components of chromatin and play an important role in innate immunity, as demonstrated by its antimicrobial activities in clams. However, little attention has been paid to histones in bivalves. To fill this gap, we investigated the genomic distribution, structural characteristics, conserved motifs, and phylogenetic relationships of histones in S. constricta. A total of 114 histone genes were detected in the S. constricta genome, which were divided into 25 types in phylogenetic analysis. Among them, partial histones exhibited a tissue-dependent expression pattern, indicating that they may be involved in sustaining the homeostasis of organs/tissues in adult S. constricta. Furthermore, mRNA expression of certain histones changed significantly in S. constricta when infected with Vibrio parahaemolyticus, suggesting that histones play a role in the immune defense of S. constricta. All together, this study on histone genes in S. constricta not only greatly expands our knowledge of histone function in the clam, but also histone evolution in molluscs.

Identification and characterization of histones in Physarum polycephalum evidence a phylogenetic vicinity of Mycetozoans to the animal kingdom

Physarum polycephalum belongs to Mycetozoans, a phylogenetic clade apart from the animal, plant and fungus kingdoms. Histones are nuclear proteins involved in genome organization and regulation and are among the most evolutionary conserved proteins within eukaryotes. Therefore, this raises the question of their conservation in Physarum and the position of this organism within the eukaryotic phylogenic tree based on histone sequences. We carried out a comprehensive study of histones in Physarum polycephalum using genomic, transcriptomic and molecular data. Our results allowed to identify the different isoforms of the core histones H2A, H2B, H3 and H4 which exhibit strong conservation of amino acid residues previously identified as subject to post-translational modifications. Furthermore, we also identified the linker histone H1, the most divergent histone, and characterized a large number of its PTMs by mass spectrometry. We also performed an in-depth investigation of histone genes and transcript structures. Histone proteins are highly conserved in Physarum and their characterization will contribute to a better understanding of the polyphyletic Mycetozoan group. Our data reinforce that P. polycephalum is evolutionary closer to animals than plants and located at the crown of the eukaryotic tree. Our study provides new insights in the evolutionary history of Physarum and eukaryote lineages.

A chemical probe based on the PreQ1 metabolite enables transcriptome-wide mapping of binding sites

The role of metabolite-responsive riboswitches in regulating gene expression in bacteria is well known and makes them useful systems for the study of RNA-small molecule interactions. Here, we study the PreQ₁ riboswitch system, assessing sixteen diverse PreQ₁-derived probes for their ability to selectively modify the class-I PreQ₁ riboswitch aptamer covalently. For the most active probe (11), a diazirine-based photocrosslinking analog of PreQ₁, X-ray crystallography and gel-based competition assays demonstrated the mode of binding of the ligand to the aptamer, and functional assays demonstrated that the probe retains activity against the full riboswitch. Transcriptome-wide mapping using Chem-CLIP revealed a highly selective interaction between the bacterial aptamer and the probe. In addition, a small number of RNA targets in endogenous human transcripts were found to bind specifically to 11, providing evidence for candidate PreQ₁ aptamers in human RNA. This work demonstrates a stark influence of linker chemistry and structure on the ability of molecules to crosslink RNA, reveals that the PreQ₁ aptamer/ligand pair are broadly useful for chemical biology applications, and provides insights into how PreQ₁, which is similar in structure to guanine, interacts with human RNAs.

Novel rRNA-depletion methods for total RNA sequencing and ribosome profiling developed for avian species

Deep sequencing of RNAs has greatly aided the study of the transcriptome, enabling comprehensive gene expression profiling and the identification of novel transcripts. While messenger RNAs (mRNAs) are of the greatest interest in gene expression studies as they encode for proteins, mRNAs make up only 3 to 5% of total RNAs, with the majority comprising ribosomal RNAs (rRNAs). Therefore, applications of deep sequencing to RNA face the challenge of how to efficiently enrich mRNA species prior to library construction. Traditional methods extract mRNAs using oligo-dT primers targeting the poly-A tail on mRNAs; however, this approach is not comprehensive as it does not capture mRNAs lacking the poly-A tail or other long non-coding RNAs that we may be interested in. Alternative mRNA enrichment methods deplete rRNAs, but such approaches require species-specific probes and the commercially available kits are costly and have only been developed for a limited number of model organisms. Here, we describe a quick, cost-effective method for depleting rRNAs using custom-designed oligos, using chickens as an example species for probe design. With this optimized protocol, we have not only removed the rRNAs from total RNAs for RNA-seq library construction but also depleted rRNA fragments from ribosome-protected fragments for ribosome profiling. Currently, this is the only rRNA depletion-based method for avian species; this method thus provides a valuable resource for both the scientific community and the poultry industry.

The RNA structurome in the asexual blood stages of malaria pathogen plasmodium falciparum

Plasmodium falciparum is a deadly human pathogen responsible for the devastating disease called malaria. In this study, we measured the differential accumulation of RNA secondary structures in coding and non-coding transcripts from the asexual developmental cycle in P. falciparum in human red blood cells. Our comprehensive analysis that combined high-throughput nuclease mapping of RNA structures by duplex RNA-seq, SHAPE-directed RNA structure validation, immunoaffinity purification and characterization of antisense RNAs collectively measured differentially base-paired RNA regions throughout the parasite's asexual RBC cycle. Our mapping data not only aligned to a diverse pool of RNAs with known structures but also enabled us to identify new structural RNA regions in the malaria genome. On average, approximately 71% of the genes with secondary structures are found to be protein coding mRNAs. The mapping pattern of these base-paired RNAs corresponded to all regions of mRNAs, including the 5' UTR, CDS and 3' UTR as well as the start and stop codons. Histone family genes which are known to form secondary structures in their mRNAs and transcripts from genes which are important for transcriptional and post-transcriptional control, such as the unique plant-like transcription factor family, ApiAP2, DNA-/RNA-binding protein, Alba3 and proteins important for RBC invasion and malaria cytoadherence also showed strong accumulation of duplex RNA reads in various asexual stages in P. falciparum. Intriguingly, our study determined stage-specific, dynamic relationships between mRNA structural contents and translation efficiency in P. falciparum asexual blood stages, suggesting an essential role of RNA structural changes in malaria gene expression programs.

U7 snRNA: A tool for gene therapy

Most U-rich small nuclear ribonucleoproteins (snRNPs) are complexes that mediate the splicing of pre-mRNAs. U7 snRNP is an exception in that it is not involved in splicing but is a key factor in the unique 3' end processing of replication-dependent histone mRNAs. However, by introducing controlled changes in the U7 snRNA histone binding sequence and in the Sm motif, it can be used as an effective tool for gene therapy. The modified U7 snRNP (U7 Sm OPT) is thus not involved in the processing of replication-dependent histone pre-mRNA but targets splicing by inducing efficient skipping or inclusion of selected exons. U7 Sm OPT is of therapeutic importance in diseases that are an outcome of splicing defects, such as myotonic dystrophy, Duchenne muscular dystrophy, amyotrophic lateral sclerosis, β-thalassemia, HIV-1 infection and spinal muscular atrophy. The benefits of using U7 Sm OPT for gene therapy are its compact size, ability to accumulate in the nucleus without causing any toxic effects in the cells, and no immunoreactivity. The risk of transgene misregulation by using U7 Sm OPT is also low because it is involved in correcting the expression of an endogenous gene controlled by its own regulatory elements. Altogether, using U7 Sm OPT as a tool in gene therapy can ensure lifelong treatment, whereas an oligonucleotide or other drug/compound would require repeated administration. It would thus be strategic to harness these unique properties of U7 snRNP and deploy it as a tool in gene therapy.

High-Resolution Mapping of Transcription Initiation in the Asexual Stages of Toxoplasma gondii

Toxoplasma gondii is a common parasite of humans and animals, causing life-threatening disease in the immunocompromized, fetal abnormalities when contracted during gestation, and recurrent ocular lesions in some patients. Central to the prevalence and pathogenicity of this protozoan is its ability to adapt to a broad range of environments, and to differentiate between acute and chronic stages. These processes are underpinned by a major rewiring of gene expression, yet the mechanisms that regulate transcription in this parasite are only partially characterized. Deciphering these mechanisms requires a precise and comprehensive map of transcription start sites (TSSs); however, Toxoplasma TSSs have remained incompletely defined. To address this challenge, we used 5'-end RNA sequencing to genomically assess transcription initiation in both acute and chronic stages of Toxoplasma. Here, we report an in-depth analysis of transcription initiation at promoters, and provide empirically-defined TSSs for 7603 (91%) protein-coding genes, of which only 1840 concur with existing gene models. Comparing data from acute and chronic stages, we identified instances of stage-specific alternative TSSs that putatively generate mRNA isoforms with distinct 5' termini. Analysis of the nucleotide content and nucleosome occupancy around TSSs allowed us to examine the determinants of TSS choice, and outline features of Toxoplasma promoter architecture. We also found pervasive divergent transcription at Toxoplasma promoters, clustered within the nucleosomes of highly-symmetrical phased arrays, underscoring chromatin contributions to transcription initiation. Corroborating previous observations, we asserted that Toxoplasma 5' leaders are among the longest of any eukaryote studied thus far, displaying a median length of approximately 800 nucleotides. Further highlighting the utility of a precise TSS map, we pinpointed motifs associated with transcription initiation, including the binding sites of the master regulator of chronic-stage differentiation, BFD1, and a novel motif with a similar positional arrangement present at 44% of Toxoplasma promoters. This work provides a critical resource for functional genomics in Toxoplasma, and lays down a foundation to study the interactions between genomic sequences and the regulatory factors that control transcription in this parasite.

Large scale, robust, and accurate whole transcriptome profiling from clinical formalin-fixed paraffin-embedded samples

Transcriptome profiling can provide information of great value in clinical decision-making, yet RNA from readily available formalin-fixed paraffin-embedded (FFPE) tissue is often too degraded for quality sequencing. To assess the clinical utility of FFPE-derived RNA, we performed ribo-deplete RNA extractions on > 3200 FFPE slide samples; 25 of these had direct FFPE vs. fresh frozen (FF) replicates, 57 were sequenced in 2 different labs, 87 underwent multiple library analyses, and 16 had direct microdissected vs. macrodissected replicates. Poly-A versus ribo-depletion RNA extraction methods were compared using transcriptomes of TCGA cohort and 3116 FFPE samples. Compared to FF, FFPE transcripts coding for nuclear/cytoplasmic proteins involved in DNA packaging, replication, and protein synthesis were detected at lower rates and zinc finger family transcripts were of poorer quality. The greatest difference in extraction methods was in histone transcripts which typically lack poly-A tails. Encouragingly, the overall sequencing success rate was 81%. Exome coverage was highly concordant in direct FFPE and FF replicates, with 98% agreement in coding exon coverage and a median correlation of whole transcriptome profiles of 0.95. We provide strong rationale for clinical use of FFPE-derived RNA based on the robustness, reproducibility, and consistency of whole transcriptome profiling.

mRNA Editing, Processing and Quality Control in Caenorhabditis elegans

While DNA serves as the blueprint of life, the distinct functions of each cell are determined by the dynamic expression of genes from the static genome. The amount and specific sequences of RNAs expressed in a given cell involves a number of regulated processes including RNA synthesis (transcription), processing, splicing, modification, polyadenylation, stability, translation, and degradation. As errors during mRNA production can create gene products that are deleterious to the organism, quality control mechanisms exist to survey and remove errors in mRNA expression and processing. Here, we will provide an overview of mRNA processing and quality control mechanisms that occur in Caenorhabditis elegans, with a focus on those that occur on protein-coding genes after transcription initiation. In addition, we will describe the genetic and technical approaches that have allowed studies in C. elegans to reveal important mechanistic insight into these processes.

Single-cell transcription analysis of Plasmodium vivax blood-stage parasites identifies stage- and species-specific profiles of expression

Plasmodium vivax and P. falciparum, the parasites responsible for most human malaria worldwide, exhibit striking biological differences, which have important clinical consequences. Unfortunately, P. vivax, unlike P. falciparum, cannot be cultivated continuously in vitro, which limits our understanding of its biology and, consequently, our ability to effectively control vivax malaria. Here, we describe single-cell gene expression profiles of 9,215 P. vivax parasites from bloodstream infections of Aotus and Saimiri monkeys. Our results show that transcription of most P. vivax genes occurs during short periods of the intraerythrocytic cycle and that this pattern of gene expression is conserved in other Plasmodium species. However, we also identify a strikingly high proportion of species-specific transcripts in late schizonts, possibly associated with the specificity of erythrocyte invasion. Our findings provide new and robust markers of blood-stage parasites, including some that are specific to the elusive P. vivax male gametocytes, and will be useful for analyzing gene expression data from laboratory and field samples.

Analysis of individual Plasmodium vivax parasites reveals the tight control of the expression of most genes during the intra-erythrocytic cycle and the differentiation of male and female gametocytes, and highlights differences between the development of P. vivax and P. falciparum.

Defective Expression of Mitochondrial, Vacuolar H+-ATPase and Histone Genes in a C. elegans Model of SMA

Spinal muscular atrophy (SMA) is a severe motor neuron degenerative disease caused by loss-of-function mutations in the survival motor neuron gene SMN1. It is widely posited that defective gene expression underlies SMA. However, the identities of these affected genes remain to be elucidated. By analyzing the transcriptome of a Caenorhabditis elegans SMA model at the pre-symptomatic stage, we found that the expression of numerous nuclear encoded mitochondrial genes and vacuolar H⁺-ATPase genes was significantly down-regulated, while that of histone genes was significantly up-regulated. We previously showed that the uaf-1 gene, encoding key splicing factor U2AF large subunit, could affect the behavior and lifespan of smn-1 mutants. Here, we found that smn-1 and uaf-1 interact to affect the recognition of 3′ and 5′ splice sites in a gene-specific manner. Altogether, our results suggest a functional interaction between smn-1 and uaf-1 in affecting RNA splicing and a potential effect of smn-1 on the expression of mitochondrial and histone genes.

How to Recruit the Correct RNA Polymerase? Lessons from snRNA Genes

Nuclear eukaryotic genomes are transcribed by three related RNA polymerases (Pol), which transcribe distinct gene sets. Specific Pol recruitment is achieved through selective core promoter recognition by basal transcription factors (TFs). Transcription by an inappropriate Pol appears to be rare and to generate mostly unstable products. A collection of short noncoding RNA genes [for example, small nuclear RNA (snRNA) or 7SK RNA genes], which play essential roles in processes such as maturation of RNA molecules or control of Pol II transcription elongation, possess highly similar core promoters, and yet are transcribed for some by Pol II and for others by Pol III as a result of small promoter differences. Here we discuss the mechanisms of selective Pol recruitment to such promoters.

Histone supply: Multitiered regulation ensures chromatin dynamics throughout the cell cycle

Mendiratta et al. review the interplay between the different regulatory layers that affect the transcription and dynamics of distinct histone H3 variants along the cell cycle.

As the building blocks of chromatin, histones are central to establish and maintain particular chromatin states associated with given cell fates. Importantly, histones exist as distinct variants whose expression and incorporation into chromatin are tightly regulated during the cell cycle. During S phase, specialized replicative histone variants ensure the bulk of the chromatinization of the duplicating genome. Other non-replicative histone variants deposited throughout the cell cycle at specific loci use pathways uncoupled from DNA synthesis. Here, we review the particular dynamics of expression, cellular transit, assembly, and disassembly of replicative and non-replicative forms of the histone H3. Beyond the role of histone variants in chromatin dynamics, we review our current knowledge concerning their distinct regulation to control their expression at different levels including transcription, posttranscriptional processing, and protein stability. In light of this unique regulation, we highlight situations where perturbations in histone balance may lead to cellular dysfunction and pathologies.

Overexpression of the double homeodomain protein DUX4c interferes with myofibrillogenesis and induces clustering of myonuclei

BACKGROUND

Facioscapulohumeral muscular dystrophy (FSHD) is associated with DNA hypomethylation at the 4q35 D4Z4 repeat array. Both the causal gene DUX4 and its homolog DUX4c are induced. DUX4c is immunodetected in every myonucleus of proliferative cells, while DUX4 is present in only 1/1000 of myonuclei where it initiates a gene deregulation cascade. FSHD primary myoblasts differentiate into either atrophic or disorganized myotubes. DUX4 expression induces atrophic myotubes and associated FSHD markers. Although DUX4 silencing normalizes the FSHD atrophic myotube phenotype, this is not the case for the disorganized phenotype. DUX4c overexpression increases the proliferation rate of human TE671 rhabdomyosarcoma cells and inhibits their differentiation, suggesting a normal role during muscle differentiation.

METHODS

By gain- and loss-of-function experiments in primary human muscle cells, we studied the DUX4c impact on proliferation, differentiation, myotube morphology, and FSHD markers.

RESULTS

In primary myoblasts, DUX4c overexpression increased the staining intensity of KI67 (a proliferation marker) in adjacent cells and delayed differentiation. In differentiating cells, DUX4c overexpression led to the expression of some FSHD markers including β-catenin and to the formation of disorganized myotubes presenting large clusters of nuclei and cytoskeletal defects. These were more severe when DUX4c was expressed before the cytoskeleton reorganized and myofibrils assembled. In addition, endogenous DUX4c was detected at a higher level in FSHD myotubes presenting abnormal clusters of nuclei and cytoskeletal disorganization. We found that the disorganized FSHD myotube phenotype could be rescued by silencing of DUX4c, not DUX4.

CONCLUSION

Excess DUX4c could disturb cytoskeletal organization and nuclear distribution in FSHD myotubes. We suggest that DUX4c up-regulation could contribute to DUX4 toxicity in the muscle fibers by favoring the clustering of myonuclei and therefore facilitating DUX4 diffusion among them. Defining DUX4c functions in the healthy skeletal muscle should help to design new targeted FSHD therapy by DUX4 or DUX4c inhibition without suppressing DUX4c normal function.

Current and Future Methods for mRNA Analysis: A Drive Toward Single Molecule Sequencing

The transcriptome encompasses a range of species including messenger RNA, and other noncoding RNA such as rRNA, tRNA, and short and long noncoding RNAs. Due to the huge role played by mRNA in development and disease, several methods have been developed to sequence and characterize mRNA, with RNA sequencing (RNA-Seq) emerging as the current method of choice particularly for large high-throughput studies. Short-read RNA-Seq which involves sequencing of short cDNA fragments and computationally assembling them to reconstruct the transcriptome, or aligning them to a reference is the most widely used approach. However, due to inherent limitations of this approach in de novo transcriptome assembly and isoform quantification, long-read RNA-Seq approaches, which also happen to be single molecule sequencing approaches, are increasingly becoming the standard for de novo transcriptome assembly and isoform quantification. In this chapter, we review the technical aspects of the current methods of RNA-Seq, both short and long-read approaches, and data analysis methods available. We discuss recent advances in single-cell RNA-Seq and direct RNA-Seq approaches, which perhaps will dominate the future of RNA-Seq.

Histone-derived piRNA biogenesis depends on the ping-pong partners Piwi5 and Ago3 in Aedes aegypti

The piRNA pathway is of key importance in controlling transposable elements in most animal species. In the vector mosquito Aedes aegypti, the presence of eight PIWI proteins and the accumulation of viral piRNAs upon arbovirus infection suggest additional functions of the piRNA pathway beyond genome defense. To better understand the regulatory potential of this pathway, we analyzed in detail host-derived piRNAs in A. aegypti Aag2 cells. We show that a large repertoire of protein-coding genes and non-retroviral integrated RNA virus elements are processed into genic piRNAs by different combinations of PIWI proteins. Among these, we identify a class of genes that produces piRNAs from coding sequences in an Ago3- and Piwi5-dependent fashion. We demonstrate that the replication-dependent histone gene family is a genic source of ping-pong dependent piRNAs and that histone-derived piRNAs are dynamically expressed throughout the cell cycle, suggesting a role for the piRNA pathway in the regulation of histone gene expression. Moreover, our results establish the Aag2 cell line as an accessible experimental model to study gene-derived piRNAs.

Birth and Death of Histone mRNAs

In metazoans, histone mRNAs are not polyadenylated but end in a conserved stem-loop. Stem-loop binding protein (SLBP) binds to the stem-loop and is required for all steps in histone mRNA metabolism. The genes for the five histone proteins are linked. A histone locus body (HLB) forms at each histone gene locus. It contains factors essential for transcription and processing of histone mRNAs, and couples transcription and processing. The active form of U7 snRNP contains the HLB component FLASH (FLICE-associated huge protein), the histone cleavage complex (HCC), and a subset of polyadenylation factors including the endonuclease CPSF73. Histone mRNAs are rapidly degraded when DNA replication is inhibited by a 3' to 5' pathway that requires extensive uridylation of mRNA decay intermediates.

Modelling the three-dimensional structure of the right-terminal domain of pospiviroids

Viroids, the smallest know plant pathogens, consist solely of a circular, single-stranded, non-coding RNA. Thus for all of their biological functions, like replication, processing, and transport, they have to present sequence or structural features to exploit host proteins. Viroid binding protein 1 (Virp1) is indispensable for replication of pospiviroids, the largest genus of viroids, in a host plant as well as in protoplasts. Virp1 is known to bind at two sites in the terminal right (TR) domain of pospiviroids; each site consists of a purine- (R-) and a pyrimidine- (Y-)rich motif that are partially base-paired to each other. Here we model the important structural features of the domain and show that it contains an internal loop of two Y · Y cis Watson-Crick/Watson-Crick (cWW) pairs, an asymmetric internal loop including a cWW and a trans Watson/Hoogsteen pair, and a thermodynamically quite stable hairpin loop with several stacking interactions. These features are discussed in connection to the known biological functions of the TR domain.

Genome-Wide Analysis and Functional Characterization of the Polyadenylation Site in Pigs Using RNAseq Data

Polyadenylation, a critical step in the production of mature mRNA for translation in most eukaryotes, involves cleavage and poly(A) tail addition at the 3′ end of mRNAs at the polyadenylation site (PAS). Sometimes, one gene can have more than one PAS, which can produce the alternative polyadenylation (APA) phenomenon and affect the stability, localization and translation of the mRNA. In this study, we discovered 28,363 PASs using pig RNAseq data, with 13,033 located in 7,403 genes. Among the genes, 41% were identified to have more than one PAS. PAS distribution analysis indicated that the PAS position was highly variable in genes. Additionally, the analysis of RNAseq data from the liver and testis showed a difference in their PAS number and usage. RT-PCR and qRT-PCR were performed to confirm our findings by detecting the expression of 3′UTR isoforms for five candidate genes. The analysis of RNAseq data under a different androstenone level and salmonella inoculation indicated that the functional usage of PAS might participate in the immune response and may be related to the androstenone level in pigs. This study provides new insights into pig PAS and facilitates further functional research of PAS.

Genome-Wide Analysis of Polyadenylation Events in Schmidtea mediterranea

In eukaryotes, 3' untranslated regions (UTRs) play important roles in regulating posttranscriptional gene expression. The 3'UTR is defined by regulated cleavage/polyadenylation of the pre-mRNA. The advent of next-generation sequencing technology has now enabled us to identify these events on a genome-wide scale. In this study, we used poly(A)-position profiling by sequencing (3P-Seq) to capture all poly(A) sites across the genome of the freshwater planarian, Schmidtea mediterranea, an ideal model system for exploring the process of regeneration and stem cell function. We identified the 3'UTRs for ∼14,000 transcripts and thus improved the existing gene annotations. We found 97 transcripts, which are polyadenylated within an internal exon, resulting in the shrinking of the ORF and loss of a predicted protein domain. Around 40% of the transcripts in planaria were alternatively polyadenylated (ApA), resulting either in an altered 3'UTR or a change in coding sequence. We identified specific ApA transcript isoforms that were subjected to miRNA mediated gene regulation using degradome sequencing. In this study, we also confirmed a tissue-specific expression pattern for alternate polyadenylated transcripts. The insights from this study highlight the potential role of ApA in regulating the gene expression essential for planarian regeneration.

Genome modification leads to phenotype reversal in human myotonic dystrophy type 1 induced pluripotent stem cell-derived neural stem cells

Myotonic dystrophy type 1 (DM1) is caused by expanded CTG repeats in the 3'-untranslated region (3' UTR) of the DMPK gene. Correcting the mutation in DM1 stem cells would be an important step toward autologous stem cell therapy. The objective of this study is to demonstrate in vitro genome editing to prevent production of toxic mutant transcripts and reverse phenotypes in DM1 stem cells. Genome editing was performed in DM1 neural stem cells (NSCs) derived from human DM1 induced pluripotent stem (iPS) cells. An editing cassette containing SV40/bGH polyA signals was integrated upstream of the CTG repeats by TALEN-mediated homologous recombination (HR). The expression of mutant CUG repeats transcript was monitored by nuclear RNA foci, the molecular hallmarks of DM1, using RNA fluorescence in situ hybridization. Alternative splicing of microtubule-associated protein tau (MAPT) and muscleblind-like (MBNL) proteins were analyzed to further monitor the phenotype reversal after genome modification. The cassette was successfully inserted into DMPK intron 9 and this genomic modification led to complete disappearance of nuclear RNA foci. MAPT and MBNL 1, 2 aberrant splicing in DM1 NSCs were reversed to normal pattern in genome-modified NSCs. Genome modification by integration of exogenous polyA signals upstream of the DMPK CTG repeat expansion prevents the production of toxic RNA and leads to phenotype reversal in human DM1 iPS-cells derived stem cells. Our data provide proof-of-principle evidence that genome modification may be used to generate genetically modified progenitor cells as a first step toward autologous cell transfer therapy for DM1.

A multiprotein occupancy map of the mRNP on the 3' end of histone mRNAs

The animal replication-dependent (RD) histone mRNAs are coordinately regulated with chromosome replication. The RD-histone mRNAs are the only known cellular mRNAs that are not polyadenylated. Instead, the mature transcripts end in a conserved stem-loop (SL) structure. This SL structure interacts with the stem-loop binding protein (SLBP), which is involved in all aspects of RD-histone mRNA metabolism. We used several genomic methods, including high-throughput sequencing of cross-linked immunoprecipitate (HITS-CLIP) to analyze the RNA-binding landscape of SLBP. SLBP was not bound to any RNAs other than histone mRNAs. We performed bioinformatic analyses of the HITS-CLIP data that included (i) clustering genes by sequencing read coverage using CVCA, (ii) mapping the bound RNA fragment termini, and (iii) mapping cross-linking induced mutation sites (CIMS) using CLIP-PyL software. These analyses allowed us to identify specific sites of molecular contact between SLBP and its RD-histone mRNA ligands. We performed in vitro crosslinking assays to refine the CIMS mapping and found that uracils one and three in the loop of the histone mRNA SL preferentially crosslink to SLBP, whereas uracil two in the loop preferentially crosslinks to a separate component, likely the 3'hExo. We also performed a secondary analysis of an iCLIP data set to map UPF1 occupancy across the RD-histone mRNAs and found that UPF1 is bound adjacent to the SLBP-binding site. Multiple proteins likely bind the 3' end of RD-histone mRNAs together with SLBP.

Translational compensation of genomic instability in neuroblastoma

Cancer-associated gene expression imbalances are conventionally studied at the genomic, epigenomic and transcriptomic levels. Given the relevance of translational control in determining cell phenotypes, we evaluated the translatome, i.e., the transcriptome engaged in translation, as a descriptor of the effects of genetic instability in cancer. We performed this evaluation in high-risk neuroblastomas, which are characterized by a low frequency of point mutations or known cancer-driving genes and by the presence of several segmental chromosomal aberrations that produce gene-copy imbalances that guide aggressiveness. We thus integrated genome, transcriptome, translatome and miRome profiles in a representative panel of high-risk neuroblastoma cell lines. We identified a number of genes whose genomic imbalance was corrected by compensatory adaptations in translational efficiency. The transcriptomic level of these genes was predictive of poor prognosis in more than half of cases, and the genomic imbalances found in their loci were shared by 27 other tumor types. This homeostatic process is also not limited to copy number-altered genes, as we showed the translational stoichiometric rebalance of histone genes. We suggest that the translational buffering of fluctuations in these dose-sensitive transcripts is a potential driving process of neuroblastoma evolution.

Response of yeast cells to high glucose involves molecular and physiological differences when compared to other osmostress conditions

Yeast cells can be affected by several causes of osmotic stress, such as high salt, sorbitol or glucose concentrations. The last condition is particularly interesting during natural processes where this microorganism participates. Response to osmostress requires the HOG (High Osmolarity Glycerol) pathway and several transcription factors, including Hot1, which plays a key role in high glucose concentrations. In this work, we describe how the yeast response to osmotic stress shows differences in accordance with the stress agent responsible for it. Compared with other conditions, under high glucose stress, delocalization of MAPK (Mitogen-Activated Protein Kinase) Hog1 is slower, induction of HOT1 expression is higher and Msn2/4 transcription factors are involved to a lesser extent. The transcriptomic analyses carried out with samples incubated for 30 min in the presence of high glucose or sorbitol reveal the presence of two functional categories with a differential expression between these conditions: glycogen biosynthesis and mobilization, and membrane-anchored proteins. We present data to demonstrate that the cells treated with 20% (w/v) (1.11 M) glucose contain higher chitin levels and are more sensitive to calcofluor white and ethanol.

Conservation and losses of non-coding RNAs in avian genomes

Here we present the results of a large-scale bioinformatics annotation of non-coding RNA loci in 48 avian genomes. Our approach uses probabilistic models of hand-curated families from the Rfam database to infer conserved RNA families within each avian genome. We supplement these annotations with predictions from the tRNA annotation tool, tRNAscan-SE and microRNAs from miRBase. We identify 34 lncRNA-associated loci that are conserved between birds and mammals and validate 12 of these in chicken. We report several intriguing cases where a reported mammalian lncRNA, but not its function, is conserved. We also demonstrate extensive conservation of classical ncRNAs (e.g., tRNAs) and more recently discovered ncRNAs (e.g., snoRNAs and miRNAs) in birds. Furthermore, we describe numerous “losses” of several RNA families, and attribute these to either genuine loss, divergence or missing data. In particular, we show that many of these losses are due to the challenges associated with assembling avian microchromosomes. These combined results illustrate the utility of applying homology-based methods for annotating novel vertebrate genomes.

CstF-64 supports pluripotency and regulates cell cycle progression in embryonic stem cells through histone 3' end processing

Embryonic stem cells (ESCs) exhibit a unique cell cycle with a shortened G1 phase that supports their pluripotency, while apparently buffering them against pro-differentiation stimuli. In ESCs, expression of replication-dependent histones is a main component of this abbreviated G1 phase, although the details of this mechanism are not well understood. Similarly, the role of 3' end processing in regulation of ESC pluripotency and cell cycle is poorly understood. To better understand these processes, we examined mouse ESCs that lack the 3' end-processing factor CstF-64. These ESCs display slower growth, loss of pluripotency and a lengthened G1 phase, correlating with increased polyadenylation of histone mRNAs. Interestingly, these ESCs also express the τCstF-64 paralog of CstF-64. However, τCstF-64 only partially compensates for lost CstF-64 function, despite being recruited to the histone mRNA 3' end-processing complex. Reduction of τCstF-64 in CstF-64-deficient ESCs results in even greater levels of histone mRNA polyadenylation, suggesting that both CstF-64 and τCstF-64 function to inhibit polyadenylation of histone mRNAs. These results suggest that CstF-64 plays a key role in modulating the cell cycle in ESCs while simultaneously controlling histone mRNA 3' end processing.

Structure and function of pre-mRNA 5'-end capping quality control and 3'-end processing

Messenger RNA precursors (pre-mRNAs) are produced as the nascent transcripts of RNA polymerase II (Pol II) in eukaryotes and must undergo extensive maturational processing, including 5′-end capping, splicing, and 3′-end cleavage and polyadenylation. This review will summarize the structural and functional information reported over the past few years on the large machinery required for the 3′-end processing of most pre-mRNAs, as well as the distinct machinery for the 3′-end processing of replication-dependent histone pre-mRNAs, which have provided great insights into the proteins and their subcomplexes in these machineries. Structural and biochemical studies have also led to the identification of a new class of enzymes (the DXO family enzymes) with activity toward intermediates of the 5′-end capping pathway. Functional studies demonstrate that these enzymes are part of a novel quality surveillance mechanism for pre-mRNA 5′-end capping. Incompletely capped pre-mRNAs are produced in yeast and human cells, in contrast to the general belief in the field that capping always proceeds to completion, and incomplete capping leads to defects in splicing and 3′-end cleavage in human cells. The DXO family enzymes are required for the detection and degradation of these defective RNAs.

Towards decrypting cryptobiosis--analyzing anhydrobiosis in the tardigrade Milnesium tardigradum using transcriptome sequencing

Background

Many tardigrade species are capable of anhydrobiosis; however, mechanisms underlying their extreme desiccation resistance remain elusive. This study attempts to quantify the anhydrobiotic transcriptome of the limno-terrestrial tardigrade Milnesium tardigradum.

Results

A prerequisite for differential gene expression analysis was the generation of a reference hybrid transcriptome atlas by assembly of Sanger, 454 and Illumina sequence data. The final assembly yielded 79,064 contigs (>100 bp) after removal of ribosomal RNAs. Around 50% of them could be annotated by SwissProt and NCBI non-redundant protein sequences. Analysis using CEGMA predicted 232 (93.5%) out of the 248 highly conserved eukaryotic genes in the assembly. We used this reference transcriptome for mapping and quantifying the expression of transcripts regulated under anhdydrobiosis in a time-series during dehydration and rehydration. 834 of the transcripts were found to be differentially expressed in a single stage (dehydration/inactive tun/rehydration) and 184 were overlapping in two stages while 74 were differentially expressed in all three stages. We have found interesting patterns of differentially expressed transcripts that are in concordance with a common hypothesis of metabolic shutdown during anhydrobiosis. This included down-regulation of several proteins of the DNA replication and translational machinery and protein degradation. Among others, heat shock proteins Hsp27 and Hsp30c were up-regulated in response to dehydration and rehydration. In addition, we observed up-regulation of ployubiquitin-B upon rehydration together with a higher expression level of several DNA repair proteins during rehydration than in the dehydration stage.

Conclusions

Most of the transcripts identified to be differentially expressed had distinct cellular function. Our data suggest a concerted molecular adaptation in M. tardigradum that permits extreme forms of ametabolic states such as anhydrobiosis. It is temping to surmise that the desiccation tolerance of tradigrades can be achieved by a constitutive cellular protection system, probably in conjunction with other mechanisms such as rehydration-induced cellular repair.

Regulation of histone gene transcription in yeast

Histones are the primary protein component of chromatin, the mixture of DNA and proteins that packages the genetic material in eukaryotes. Large amounts of histones are required during the S phase of the cell cycle when genome replication occurs. However, ectopic expression of histones during other cell cycle phases is toxic; thus, histone expression is restricted to the S phase and is tightly regulated at multiple levels, including transcriptional, post-transcriptional, translational, and post-translational. In this review, we discuss mechanisms of regulation of histone gene expression with emphasis on the transcriptional regulation of the replication-dependent histone genes in the model yeast Saccharomyces cerevisiae.

Hyper conserved elements in vertebrate mRNA 3'-UTRs reveal a translational network of RNA-binding proteins controlled by HuR

Little is known regarding the post-transcriptional networks that control gene expression in eukaryotes. Additionally, we still need to understand how these networks evolve, and the relative role played in them by their sequence-dependent regulatory factors, non-coding RNAs (ncRNAs) and RNA-binding proteins (RBPs). Here, we used an approach that relied on both phylogenetic sequence sharing and conservation in the whole mapped 3'-untranslated regions (3'-UTRs) of vertebrate species to gain knowledge on core post-transcriptional networks. The identified human hyper conserved elements (HCEs) were predicted to be preferred binding sites for RBPs and not for ncRNAs, namely microRNAs and long ncRNAs. We found that the HCE map identified a well-known network that post-transcriptionally regulates histone mRNAs. We were then able to discover and experimentally confirm a translational network composed of RNA Recognition Motif (RRM)-type RBP mRNAs that are positively controlled by HuR, another RRM-type RBP. HuR shows a preference for these RBP mRNAs bound in stem-loop motifs, confirming its role as a 'regulator of regulators'. Analysis of the transcriptome-wide HCE distribution revealed a profile of prevalently small clusters separated by unconserved intercluster RNA stretches, which predicts the formation of discrete small ribonucleoprotein complexes in the 3'-UTRs.

Emergence of the β-CASP ribonucleases: highly conserved and ubiquitous metallo-enzymes involved in messenger RNA maturation and degradation

The β-CASP ribonucleases, which are found in the three domains of life, have in common a core of 460 residues containing seven conserved sequence motifs involved in the tight binding of two catalytic zinc ions. A hallmark of these enzymes is their ability to catalyze both endo- and exo-ribonucleolytic degradation. Exo-ribonucleolytic degradation proceeds in the 5' to 3' direction and is sensitive to the phosphorylation state of the 5' end of a transcript. Recent phylogenomic analyses have shown that the β-CASP ribonucleases can be partitioned into two major subdivisions that correspond to orthologs of eukaryal CPSF73 and bacterial RNase J. We discuss the known functions of the CPSF73 and RNase J orthologs, their association into complexes, and their structure as it relates to mechanism of action. Eukaryal CPSF73 is part of a large multiprotein complex that is involved in the maturation of the 3' end of RNA Polymerase II transcripts and the polyadenylation of messenger RNA. RNase J1 and J2 are paralogs in Bacillus subtilis that are involved in the degradation of messenger RNA and the maturation of non-coding RNA. RNase J1 and J2 co-purify as a heteromeric complex and there is recent evidence that they interact with other enzymes to form a bacterial RNA degradosome. Finally, we speculate on the evolutionary origin of β-CASP ribonucleases and on their functions in Archaea. Orthologs of CPSF73 with endo- and exo-ribonuclease activity are strictly conserved throughout the archaea suggesting a role for these enzymes in the maturation and/or degradation of messenger RNA. This article is part of a Special Issue entitled: RNA Decay mechanisms.

Structure of histone mRNA stem-loop, human stem-loop binding protein, and 3'hExo ternary complex

Metazoan replication-dependent histone messenger RNAs (mRNAs) have a conserved stem-loop (SL) at their 3'-end. The stem-loop binding protein (SLBP) specifically recognizes the SL to regulate histone mRNA metabolism, and the 3'-5' exonuclease 3'hExo trims its 3'-end after processing. We report the crystal structure of a ternary complex of human SLBP RNA binding domain, human 3'hExo, and a 26-nucleotide SL RNA. Only one base of the SL is recognized specifically by SLBP, and the two proteins primarily recognize the shape of the RNA. SLBP and 3'hExo have no direct contact with each other, and induced structural changes in the loop of the SL mediate their cooperative binding. The 3' flanking sequence is positioned in the 3'hExo active site, but the ternary complex limits the extent of trimming.

Novel 3' ends that support translation

The 3' ends of two large noncoding RNAs, MALAT1 (metastasis-associated lung adenocarcinoma transcript 1) and MEN β, are formed by cleavage by RNase P and are capped but not polyadenylated. In the November 1, 2012, issue of Genes & Development, Wilusz and colleagues (pp. 2392-2407) show that when these 3' ends are formed on a GFP reporter, the resulting mRNA is exported to the cytoplasm and translated. The 3' end forms a novel triple-helical structure that supports export and translation as well as a poly(A) tail does.

Vita-PAMPs: signatures of microbial viability

Can the innate immune system detect and respond to microbial viability? Using bacteria as a model, we found that indeed the very essence of microbial infectivity, viability itself, can be detected, and notably, in the absence of the activity of virulence factors. The microbial molecule that serves as the signature of viability is bacterial messenger RNA (mRNA), common to all bacteria, and without which bacteria cannot survive. Prokaryotic mRNAs also differ from eukaryotic mRNAs in several ways, and as such, these features all fulfill the criteria, and more, for a pathogen-associated molecular pattern (PAMP) as originally proposed by Charles Janeway. Because these mRNAs are lost from dead bacteria, they belong to a special class of PAMPs, which we call vita-PAMPs. Here we discuss the possible receptors and pathways involved in the detection of bacterial mRNAs, and thus microbial viability. We also consider examples of vita-PAMPs other than bacterial mRNA.

Biology and Mechanisms of Short RNAs in Caenorhabditis elegans

The significance of noncoding RNAs in animal biology is being increasingly recognized. The nematode Caenorhabditis elegans has an extensive system of short RNAs that includes microRNAs, piRNAs, and endogenous siRNAs, which regulate development, control life span, provide resistance to viruses and transposons, and monitor gene duplications. Progress in our understanding of short RNAs was stimulated by the discovery of RNA interference, a phenomenon of sequence-specific gene silencing induced by exogenous double-stranded RNA, at the turn of the twenty-first century. This chapter provides a broad overview of the exogenous and endogenous RNAi processes in C. elegans and describes recent advances in genetic, genomic, and molecular analyses of nematode's short RNAs and proteins involved in the RNAi-related pathways.

Systematic reconstruction of RNA functional motifs with high-throughput microfluidics

We present RNA-mechanically induced trapping of molecular interactions (RNA-MITOMI), a microfluidic platform that allows integrated synthesis and functional assays for programmable RNA libraries. The interaction of a comprehensive library of RNA mutants with stem-loop-binding protein precisely defined the RNA structural and sequence features that govern affinity. The functional motif reconstructed in a single experiment on our platform uncovers new binding specificities and enriches interpretation of phylogenetic data.

CSR-1 RNAi pathway positively regulates histone expression in C. elegans

Endogenous small interfering RNAs (endo-siRNAs) have been discovered in many organisms, including mammals. In C. elegans, depletion of germline-enriched endo-siRNAs found in complex with the CSR-1 Argonaute protein causes sterility and defects in chromosome segregation in early embryos. We discovered that knockdown of either csr-1, the RNA-dependent RNA polymerase (RdRP) ego-1, or the dicer-related helicase drh-3, leads to defects in histone mRNA processing, resulting in severe depletion of core histone proteins. The maturation of replication-dependent histone mRNAs, unlike that of other mRNAs, requires processing of their 3'UTRs through an endonucleolytic cleavage guided by the U7 snRNA, which is lacking in C. elegans. We found that CSR-1-bound antisense endo-siRNAs match histone mRNAs and mRNA precursors. Consistently, we demonstrate that CSR-1 directly binds to histone mRNA in an ego-1-dependent manner using biotinylated 2'-O-methyl RNA oligonucleotides. Moreover, we demonstrate that increasing the dosage of histone genes rescues the lethality associated with depletion of CSR-1 and EGO-1. These results support a positive and direct effect of RNAi on histone gene expression.

Regulation of leaf morphology by microRNA394 and its target LEAF CURLING RESPONSIVENESS

The present study identified Arabidopsis miR394 and its target, an F-box (SKP1-Cullin/CDC53-F-box) gene At1g27340 (here referred to as LEAF CURLING RESPONSIVENESS, LCR), for regulation of leaf curling-related morphology. The loss-of-function lcr mutants exhibit pleiotropic defects with semi-dwarfism, altered leaf shape and a shorter stem. Overexpression of an miR394-resistant version of LCR under the 35S promoter (35S:m5LCR) and target mimicry MIM394 resulted in a curled-down leaf defect. Conversely, transgenic plants overexpressing 35S:MIR394a/b display a curled-up leaf phenotype. Detailed analyses show that there is a certain level of LCR that is optimal for leaf morphology, but lower or higher levels lead to abnormal leaf development, indicating that expression of miR394 in the leaf lamina is necessary for proper leaf morphology. Because the phytohormone auxin plays a crucial role in leaf morphogenesis and patterning, the DR5-GUS reporter gene was used to monitor the auxin response. We show that DR5 expression patterns in lcr and 35S::m5LCR plants were significantly different from those in the wild type. Also, overexpression of LCR in 35S::m5LCR plants drastically decreased the expression of the auxin-responsive genes IAA3, AXR3 and IAMT1, whereas increased expression of the genes was found in 35S::MIR394a plants. These results indicate that miR394 and its target LCR are involved in the regulation of leaf development.

Interaction of the histone mRNA hairpin with stem-loop binding protein (SLBP) and regulation of the SLBP-RNA complex by phosphorylation and proline isomerization

In metazoans, the majority of histone proteins are generated from replication-dependent histone mRNAs. These mRNAs are unique in that they are not polyadenylated but have a stem-loop structure in their 3' untranslated region. An early event in 3' end formation of histone mRNAs is the binding of stem-loop binding protein (SLBP) to the stem-loop structure. Here we provide insight into the mechanism by which SLBP contacts the histone mRNA. There are two binding sites in the SLBP RNA binding domain for the histone mRNA hairpin. The first binding site (Glu129-Val158) consists of a helix-turn-helix motif that likely recognizes the unpaired uridines in the loop of the histone hairpin and, upon binding, destabilizes the first G-C base pair at the base of the stem. The second binding site lies between residues Arg180 and Pro200, which appears to recognize the second G-C base pair from the base of the stem and possibly regions flanking the stem-loop structure. We show that the SLBP-histone mRNA complex is regulated by threonine phosphorylation and proline isomerization in a conserved TPNK sequence that lies between the two binding sites. Threonine phosphorylation increases the affinity of SLBP for histone mRNA by slowing the off rate for complex dissociation, whereas the adjacent proline acts as a critical hinge that may orient the second binding site for formation of a stable SLBP-histone mRNA complex. The nuclear magnetic resonance and kinetic studies presented here provide a framework for understanding how SLBP recognizes histone mRNA and highlight possible structural roles of phosphorylation and proline isomerization in RNA binding proteins in remodeling ribonucleoprotein complexes.

Identification of the enzyme responsible for N1-methylation of pseudouridine 54 in archaeal tRNAs

tRNAs from all three kingdoms of life contain a variety of modified nucleotides required for their stability, proper folding, and accurate decoding. One prominent example is the eponymous ribothymidine (rT) modification at position 54 in the T-arm of eukaryotic and bacterial tRNAs. In contrast, in most archaea this position is occupied by another hypermodified nucleotide: the isosteric N1-methylated pseudouridine. While the enzyme catalyzing pseudouridine formation at this position is known, the pseudouridine N1-specific methyltransferase responsible for this modification has not yet been experimentally identified. Here, we present biochemical and genetic evidence that the two homologous proteins, Mja_1640 (COG 1901, Pfam DUF358) and Hvo_1989 (Pfam DUF358) from Methanocaldococcus jannaschii and Haloferax volcanii, respectively, are representatives of the methyltransferase responsible for this modification. However, the in-frame deletion of the pseudouridine N1-methyltransferase gene in H. volcanii did not result in a discernable phenotype in line with similar observations for knockouts of other T-arm methylating enzymes.