Modified nucleotides in tRNAs are important determinants of folding, structure and function. Here we identify METTL8 as a mitochondrial matrix protein and active RNA methyltransferase responsible for installing m3C32 in the human mitochondrial (mt-)tRNAThr and mt-tRNASer(UCN). METTL8 crosslinks to the anticodon stem loop (ASL) of many mt-tRNAs in cells, raising the question of how methylation target specificity is achieved. Dissection of mt-tRNA recognition elements revealed U34G35 and t6A37/(ms2)i6A37, present concomitantly only in the ASLs of the two substrate mt-tRNAs, as key determinants for METTL8-mediated methylation of C32. Several lines of evidence demonstrate the influence of U34, G35, and the m3C32 and t6A37/(ms2)i6A37 modifications in mt-tRNAThr/Ser(UCN) on the structure of these mt-tRNAs. Although mt-tRNAThr/Ser(UCN) lacking METTL8-mediated m3C32 are efficiently aminoacylated and associate with mitochondrial ribosomes, mitochondrial translation is mildly impaired by lack of METTL8. Together these results define the cellular targets of METTL8 and shed new light on the role of m3C32 within mt-tRNAs.

The influence of posttranscriptional modification on structural stabilization of tRNA from hyperthermophilic archaea was studied, using Pyrococcus furiosus (growth optimum 100 degrees C) as a primary model. Optical melting temperatures (Tm) of unfractionated tRNA in 20 mM Mg2+ are 97 degrees C for P. furiosus and 101.5 degrees C for Pyrodictium occultum (growth optimum, 105 degrees C). These values are approximately 20 degrees C higher than predicted solely from G-C content and are attributed primarily to posttranscriptional modification. Twenty-three modified nucleosides were determined in total digests of P. furiosus tRNA by combined HPLC-mass spectrometry. From cells cultured at 70, 85, and 100 degrees C, progressively increased levels of modification were observed within three families of nucleosides, the most highly modified forms of which were N4-acetyl-2'-O-methylcytidine (ac4Cm), N2,N2,2'-O-trimethylguanosine (m2(2)Gm), and 5-methyl-2-thiouridine (m5s2U). Nucleosides ac4Cm and m2(2)Gm, which are unique to the archaeal hyperthermophiles, were shown in earlier NMR studies to exhibit unusually high conformational stabilities that favor the C3'-endo form [Kawai, G., et al. (1991) Nucleic Acids Symp. Ser. 21, 49-50; (1992) Nucleosides Nucleotides 11, 759-771]. The sequence location of m5s2U was determined by mass spectrometry to be primarily at tRNA position 54, a site of known thermal stabilization in the bacterial thermophile Thermus thermophilus [Horie, N., et al. (1985) Biochemistry 24, 5711-5715]. It is concluded that selected posttranscriptional modifications in archaeal thermophiles play major stabilizing roles beyond the effects of Mg2+ binding and G-C content, and are proportionally more important and have evolved with greater structural diversity at the nucleoside level in the bacterial thermophiles.

Although the involvement of Ser/Arg-rich (SR) proteins in RNA metabolism is well documented, their role in vertebrate development remains elusive. We, therefore, elected to take advantage of the zebrafish model organism to study the SR genes' functions using the splicing morpholino (sMO) microinjection and the programmable site-specific nucleases. Consistent with previous research, we revealed discrepancies between the mutant and morphant phenotypes and we show that these inconsistencies may result from a large number of unsuspected inadvertent morpholino RNA targets. While microinjection of MOs directed against srsf5a (sMOsrsf5a) led to developmental defects, the corresponding homozygous mutants did not display any phenotypic traits. Furthermore, microinjection of sMOsrsf5a into srsf5a-/- led to the previously observed morphant phenotype. Similar findings were observed for other SR genes. sMOsrsf5a alternative target genes were identified using deep mRNA sequencing. We uncovered that only 11 consecutive bases complementary to sMOsrsf5a are sufficient for binding and subsequent blocking of splice sites. In addition, we observed that sMOsrsf5a secondary targets can be reduced by increasing embryos growth temperature after microinjection. Our data contribute to the debate about MO specificity, efficacy and the number of unknown targeted sequences.

Transfer RNA[Ser]Sec carries multiple post-transcriptional modifications. The A37G mutation in tRNA[Ser]Sec abrogates isopentenylation of base 37 and has a profound effect on selenoprotein expression in mice. Patients with a homozygous pathogenic p.R323Q variant in tRNA-isopentenyl-transferase (TRIT1) show a severe neurological disorder, and hence we wondered whether selenoprotein expression was impaired. Patient fibroblasts with the homozygous p.R323Q variant did not show a general decrease in selenoprotein expression. However, recombinant human TRIT1R323Q had significantly diminished activities towards several tRNA substrates in vitro. We thus engineered mice conditionally deficient in Trit1 in hepatocytes and neurons. Mass-spectrometry revealed that hypermodification of U34 to mcm5Um occurs independently of isopentenylation of A37 in tRNA[Ser]Sec. Western blotting and 75Se metabolic labeling showed only moderate effects on selenoprotein levels and 75Se incorporation. A detailed analysis of Trit1-deficient liver using ribosomal profiling demonstrated that UGA/Sec re-coding was moderately affected in Selenop, Txnrd1, and Sephs2, but not in Gpx1. 2'O-methylation of U34 in tRNA[Ser]Sec depends on FTSJ1, but does not affect UGA/Sec re-coding in selenoprotein translation. Taken together, our results show that a lack of isopentenylation of tRNA[Ser]Sec affects UGA/Sec read-through but differs from a A37G mutation.

Protein kinases are key signaling enzymes that catalyze the transfer of γ-phosphate from an ATP molecule to a phospho-accepting residue in the substrate. Unraveling the molecular features that govern the preference of kinases for particular residues flanking the phosphoacceptor is important for understanding kinase specificities toward their substrates and for designing substrate-like peptidic inhibitors. We applied ANCHORSmap, a new fragment-based computational approach for mapping amino acid side chains on protein surfaces, to predict and characterize the preference of kinases toward Arginine binding. We focus on positions P-2 and P-5, commonly occupied by Arginine (Arg) in substrates of basophilic Ser/Thr kinases. The method accurately identified all the P-2/P-5 Arg binding sites previously determined by X-ray crystallography and produced Arg preferences that corresponded to those experimentally found by peptide arrays. The predicted Arg-binding positions and their associated pockets were analyzed in terms of shape, physicochemical properties, amino acid composition, and in-silico mutagenesis, providing structural rationalization for previously unexplained trends in kinase preferences toward Arg moieties. This methodology sheds light on several kinases that were described in the literature as having non-trivial preferences for Arg, and provides some surprising departures from the prevailing views regarding residues that determine kinase specificity toward Arg. In particular, we found that the preference for a P-5 Arg is not necessarily governed by the 170/230 acidic pair, as was previously assumed, but by several different pairs of acidic residues, selected from positions 133, 169, and 230 (PKA numbering). The acidic residue at position 230 serves as a pivotal element in recognizing Arg from both the P-2 and P-5 positions.

Nedd4 (Nedd4-1) is an E3 ubiquitin ligase involved in crucial biological processes such as growth factor receptor signaling. While canonical Nedd4-1 comprises a C2-WW(4)-HECT domain architecture, alternative splicing produces non-canonical isoforms that are poorly characterized. Here we characterized Nedd4-1(NE), a primate-specific isoform of Nedd4-1 that contains a large N-terminal Extension (NE) that replaces most of the C2 domain. We show that Nedd4-1(NE) mRNA is ubiquitously expressed in human tissues and cell lines. Moreover, we found that Nedd4-1(NE) is more active than the canonical Nedd4-1 isoform, likely due to the absence of a C2 domain-mediated autoinhibitory mechanism. Additionally, we identified two Thr/Ser phosphoresidues in the NE region that act as binding sites for 14-3-3 proteins, and show that phosphorylation on these sites reduces substrate binding. Finally, we show that the NE region can act as a binding site for the RPB2 subunit of RNA polymerase II, a unique substrate of Nedd4-1(NE) but not the canonical Nedd4-1. Taken together, our results demonstrate that alternative splicing of the ubiquitin ligase Nedd4-1 can produce isoforms that differ in their catalytic activity, binding partners and substrates, and mechanisms of regulation.

Bacteria expand their genetic diversity, spread antibiotic resistance genes, and obtain virulence factors through the highly coordinated process of conjugative plasmid transfer (CPT). A plasmid-encoded relaxase enzyme initiates and terminates CPT by nicking and religating the transferred plasmid in a sequence-specific manner. We solved the 2.3 A crystal structure of the relaxase responsible for the spread of the resistance plasmid pCU1 and determined its DNA binding and nicking capabilities. The overall fold of the pCU1 relaxase is similar to that of the F plasmid and plasmid R388 relaxases. However, in the pCU1 structure, the conserved tyrosine residues (Y18,19,26,27) that are required for DNA nicking and religation were displaced up to 14 A out of the relaxase active site, revealing a high degree of mobility in this region of the enzyme. In spite of this flexibility, the tyrosines still cleaved the nic site of the plasmid's origin of transfer, and did so in a sequence-specific, metal-dependent manner. Unexpectedly, the pCU1 relaxase lacked the sequence-specific DNA binding previously reported for the homologous F and R388 relaxase enzymes, despite its high sequence and structural similarity with both proteins. In summary, our work outlines novel structural and functional aspects of the relaxase-mediated conjugative transfer of plasmid pCU1.

Horizontal gene transfer (HGT) has been widely suggested to play a critical role in the environmental adaptation of microbes; however, the number and origin of the genes in microbial genomes obtained through HGT remain unknown as the frequency of detected HGT events is generally underestimated, particularly in the absence of information on donor sequences. As an alternative to phylogeny-based methods that rely on sequence alignments, we have developed an alignment-free clustering method on the basis of an unsupervised neural network "Batch-Learning Self-Organizing Map (BLSOM)" in which sequence fragments are clustered based solely on oligonucleotide similarity without taxonomical information, to detect HGT candidates and their origin in entire genomes. By mapping the microbial genomic sequences on large-scale BLSOMs constructed with nearly all prokaryotic genomes, HGT candidates can be identified, and their origin assigned comprehensively, even for microbial genomes that exhibit high novelty. By focusing on two types of Alphaproteobacteria, specifically psychrotolerant Sphingomonas strains from an Antarctic lake, we detected HGT candidates using BLSOM and found higher proportions of HGT candidates from organisms belonging to Betaproteobacteria in the genomes of these two Antarctic strains compared with those of continental strains. Further, an origin difference was noted in the HGT candidates found in the two Antarctic strains. Although their origins were highly diversified, gene functions related to the cell wall or membrane biogenesis were shared among the HGT candidates. Moreover, analyses of amino acid frequency suggested that housekeeping genes and some HGT candidates of the Antarctic strains exhibited different characteristics to other continental strains. Lys, Ser, Thr, and Val were the amino acids found to be increased in the Antarctic strains, whereas Ala, Arg, Glu, and Leu were decreased. Our findings strongly suggest a low-temperature adaptation process for microbes that may have arisen convergently as an independent evolutionary strategy in each Antarctic strain. Hence, BLSOM analysis could serve as a powerful tool in not only detecting HGT candidates and their origins in entire genomes, but also in providing novel perspectives into the environmental adaptations of microbes.

The sorting of RNA molecules to subcellular locations facilitates the activity of spatially restricted processes. We have analyzed subcellular transcriptomes of FMRP-null mouse neuronal cells to identify transcripts that depend on FMRP for efficient transport to neurites. We found that these transcripts contain an enrichment of G-quadruplex sequences in their 3' UTRs, suggesting that FMRP recognizes them to promote RNA localization. We observed similar results in neurons derived from Fragile X Syndrome patients. We identified the RGG domain of FMRP as important for binding G-quadruplexes and the transport of G-quadruplex-containing transcripts. Finally, we found that the translation and localization targets of FMRP were distinct and that an FMRP mutant that is unable to bind ribosomes still promoted localization of G-quadruplex-containing messages. This suggests that these two regulatory modes of FMRP may be functionally separated. These results provide a framework for the elucidation of similar mechanisms governed by other RNA-binding proteins.

One of the major mechanisms of post-transcriptional gene regulation is achieved by proteins bearing well-defined sequence motifs involved in 'RNA binding'. In eukaryotes, RNA binding proteins (RBPs) are key players of RNA metabolism that includes synthesis, processing, editing, modifying, transport, storage and stability of RNA. In plants, the family of RBPs is vastly expanded compared to other eukaryotes including humans. In this study we identified 363 RBPs in the barley genome. Gene ontology enrichment analysis of barley RBPs indicated these proteins were in all the major cellular compartments and associated with key biological processes including translation, splicing, seed development and stress signaling. Members with the classical RNA binding motifs such as the RNA recognition motif (RRM), KH domain, Helicase, CRM, dsRNA and Pumilio were identified in the repertoire of barley RBPs. Similar to Arabidopsis, the RRM containing RBPs were the most abundant in barley genome. In-depth analysis of the RRM containing proteins - polyA binding proteins, Ser/Arg rich proteins and Glycine-rich RBPs were undertaken. Reanalysis of the proteome dataset of various stages during barley malting identified 38 RBPs suggesting an important role for these proteins during the malting process. This survey provides a systematic analysis of barley RBPs and serves as the basis for the further functional characterization of this important family of proteins.

Horizontal gene transfer is a means by which bacteria, archaea, and eukaryotes are able to trade DNA within and between species. While there are a variety of mechanisms through which this genetic exchange can take place, one means prevalent in the archaeon Haloferax volcanii involves the transient formation of cytoplasmic bridges between cells and is referred to as mating. This process can result in the exchange of very large fragments of DNA between the participating cells. Genes governing the process of mating, including triggers to initiate mating, mechanisms of cell fusion, and DNA exchange, have yet to be characterized. We used a transcriptomic approach to gain a more detailed knowledge of how mating might transpire. By examining the differential expression of genes expressed in cells harvested from mating conditions on a filter over time and comparing them to those expressed in a shaking culture, we were able to identify genes and pathways potentially associated with mating. These analyses provide new insights into both the mechanisms and barriers of mating in Hfx. volcanii.

Herein we present the tRNA core hypothesis, which emphasizes the central role of tRNAs molecules in the origin and evolution of fundamental biological processes. tRNAs gave origin to the first genes (mRNA) and the peptidyl transferase center (rRNA), proto-tRNAs were at the core of a proto-translation system, and the anticodon and operational codes then arose in tRNAs molecules. Metabolic pathways emerged from evolutionary pressures of the decoding systems. The transitions from the RNA world to the ribonucleoprotein world to modern biological systems were driven by three kinds of tRNAs transitions, to wit, tRNAs leading to both mRNA and rRNA.

Prolonged manganese exposure causes manganism, a neurodegenerative movement disorder. The identity of adaptive and nonadaptive cellular processes targeted by manganese remains mostly unexplored. Here we study mechanisms engaged by manganese in genetic cellular models known to increase susceptibility to manganese exposure, the plasma membrane manganese efflux transporter SLC30A10 and the mitochondrial Parkinson's gene PARK2. We found that SLC30A10 and PARK2 mutations as well as manganese exposure compromised the mitochondrial RNA granule composition and function, resulting in disruption of mitochondrial transcript processing. These RNA granule defects led to impaired assembly and function of the mitochondrial respiratory chain. Notably, cells that survived a cytotoxic manganese challenge had impaired RNA granule function, thus suggesting that this granule phenotype was adaptive. CRISPR gene editing of subunits of the mitochondrial RNA granule, FASTKD2 or DHX30, as well as pharmacological inhibition of mitochondrial transcription-translation, were protective rather than deleterious for survival of cells acutely exposed to manganese. Similarly, adult Drosophila mutants with defects in the mitochondrial RNA granule component scully were safeguarded from manganese-induced mortality. We conclude that impairment of the mitochondrial RNA granule function is a protective mechanism for acute manganese toxicity.

Target-specific integration of transposable elements for multicopy genes, such as ribosomal RNA and small nuclear RNA (snRNA) genes, is of great interest because of the relatively harmless nature, stable inheritance and possible application for targeted gene delivery of target-specific transposable elements. To date, such strict target specificity has been observed only among non-LTR retrotransposons. We here report a new superfamily of sequence-specific DNA transposons, designated Dada. Dada encodes a DDE-type transposase that shows a distant similarity to transposases encoded by eukaryotic MuDR, hAT, P and Kolobok transposons, as well as the prokaryotic IS256 insertion element. Dada generates 6-7 bp target site duplications upon insertion. One family of Dada DNA transposons targets a specific site inside the U6 snRNA genes and are found in various fish species, water flea, oyster and polycheate worm. Other target sequences of the Dada transposons are U1 snRNA genes and different tRNA genes. The targets are well conserved in multicopy genes, indicating that copy number and sequence conservation are the primary constraints on the target choice of Dada transposons. Dada also opens a new frontier for target-specific gene delivery application.

Tudor-interacting repair regulator (TIRR) is an RNA-binding protein and a negative regulator of the DNA-repair factor p53-binding protein 1 (53BP1). In non-damage conditions, TIRR is bound to 53BP1. After DNA damage, TIRR and 53BP1 dissociate, and 53BP1 binds the chromatin at the double-strand break (DSB) to promote non-homologous end joining (NHEJ)-mediated repair. However, the exact mechanistic details of this dissociation after damage are unknown. Increasing evidence has implicated RNA as a crucial factor in the DNA damage response (DDR). Here, we show that RNA can separate TIRR/53BP1. Specifically, RNA with a hairpin secondary structure, transcribed at the DSB by RNA polymerase II (RNAPII), promotes TIRR/53BP1 complex separation. This hairpin RNA binds to the same residues on TIRR as 53BP1. Our results uncover a role of DNA-damage-derived RNA in modulating a protein-protein interaction and contribute to our understanding of DSB repair.

Parasitoid wasps are important natural enemies of arthropod hosts in natural and agricultural ecosystems and are often associated with viruses or virion-like particles. Here, we report a novel negative-stranded RNA virus from a parasitoid wasp (Pteromalus puparum). The complete viral genome is 12,230 nucleotides in length, containing five non-overlapping, linearly arranged open reading frames. Phylogenetically, the virus clusters with and is a novel member of the mononegaviral family Nyamiviridae, here designated as Pteromalus puparum negative-strand RNA virus 1 (PpNSRV-1). PpNSRV-1 is present in various tissues and life stages of the parasitoid wasp, and is transmitted vertically through infected females and males. Virus infections in field populations of P. puparum wasps ranged from 16.7 to 37.5%, without linearly correlating with temperature. PpNSRV-1 increased adult longevity and impaired several fitness parameters of the wasp, but had no influence on successful parasitism. Strikingly, PpNSRV-1 mediated the offspring sex ratio by decreasing female offspring numbers. RNA interference knockdown of virus open reading frame I eliminated these PpNSRV-1-induced effects. Thus, we infer that PpNSRV-1 has complex effects on its insect host including sex ratio distortion towards males, as well as possible mutualistic benefits through increasing wasp longevity.

During eukaryotic transcription, RNA polymerase II undergoes dynamic post-translational modifications on the C-terminal domain (CTD) of the largest subunit, generating an information-rich PTM landscape that transcriptional regulators bind. The phosphorylation of Ser5 and Ser2 of CTD heptad occurs spatiotemporally with the transcriptional stages, recruiting different transcriptional regulators to Pol II. To delineate the protein interactomes at different transcriptional stages, we reconstructed phosphorylation patterns of the CTD at Ser5 and Ser2 in vitro. Our results showed that distinct protein interactomes are recruited to RNA polymerase II at different stages of transcription by the phosphorylation of Ser2 and Ser5 of the CTD heptads. In particular, we characterized calcium homeostasis endoplasmic reticulum protein (CHERP) as a regulator bound by phospho-Ser2 heptad. Pol II association with CHERP recruits an accessory splicing complex whose loss results in broad changes in alternative splicing events. Our results shed light on the PTM-coded recruitment process that coordinates transcription.

Circular RNAs (circRNAs) are widely expressed in eukaryotes and are regulated in many biological processes. Although several studies indicate their activity as microRNA (miRNA) and protein sponges, little is known about their ability to directly control mRNA homeostasis. We show that the widely expressed circZNF609 directly interacts with several mRNAs and increases their stability and/or translation by favoring the recruitment of the RNA-binding protein ELAVL1. Particularly, the interaction with CKAP5 mRNA, which interestingly overlaps the back-splicing junction, enhances CKAP5 translation, regulating microtubule function in cancer cells and sustaining cell-cycle progression. Finally, we show that circZNF609 downregulation increases the sensitivity of several cancer cell lines to different microtubule-targeting chemotherapeutic drugs and that locked nucleic acid (LNA) protectors against the pairing region on circZNF609 phenocopy such effects. These data set an example of how the small effects tuned by circZNF609/CKAP5 mRNA interaction might have a potent output in tumor growth and drug response.

SRSF1 governs splicing of over 1,500 mRNA transcripts. SRSF1 contains two RNA-recognition motifs (RRMs) and a C-terminal Arg/Ser-rich region (RS). It has been thought that SRSF1 RRMs exclusively recognize single-stranded exonic splicing enhancers, while RS lacks RNA-binding specificity. With our success in solving the insolubility problem of SRSF1, we can explore the unknown RNA-binding landscape of SRSF1. We find that SRSF1 RS prefers purine over pyrimidine. Moreover, SRSF1 binds to the G-quadruplex (GQ) from the ARPC2 mRNA, with both RRMs and RS being crucial. Our binding assays show that the traditional RNA-binding sites on the RRM tandem and the Arg in RS are responsible for GQ binding. Interestingly, our FRET and circular dichroism data reveal that SRSF1 unfolds the ARPC2 GQ, with RS leading unfolding and RRMs aiding. Our saturation transfer difference NMR results discover that Arg residues in SRSF1 RS interact with the guanine base but other nucleobases, underscoring the uniqueness of the Arg/guanine interaction. Our luciferase assays confirm that SRSF1 can alleviate the inhibitory effect of GQ on gene expression in the cell. Given the prevalence of RNA GQ and SR proteins, our findings unveil unexplored SR protein functions with broad implications in RNA splicing and translation.

The highly conserved, Nxf/Nxt (TAP/p15) RNA nuclear export pathway is important for export of most mRNAs from the nucleus, by interacting with mRNAs and promoting their passage through nuclear pores. Nxt1 is essential for viability; using a partial loss of function allele, we reveal a role for this gene in tissue specific transcription. We show that many Drosophila melanogaster testis-specific mRNAs require Nxt1 for their accumulation. The transcripts that require Nxt1 also depend on a testis-specific transcription complex, tMAC. We show that loss of Nxt1 leads to reduced transcription of tMAC targets. A reporter transcript from a tMAC-dependent promoter is under-expressed in Nxt1 mutants, however the same transcript accumulates in mutants if driven by a tMAC-independent promoter. Thus, in Drosophila primary spermatocytes, the transcription factor used to activate expression of a transcript, rather than the RNA sequence itself or the core transcription machinery, determines whether this expression requires Nxt1. We additionally find that transcripts from intron-less genes are more sensitive to loss of Nxt1 function than those from intron-containing genes and propose a mechanism in which transcript processing feeds back to increase activity of a tissue specific transcription complex.