Anionic G•U pairs in bacterial ribosomal rRNAs.

Wobble GU pairs (or G•U) occur frequently within double-stranded RNA helices interspersed between standard G=C and A-U Watson-Crick pairs. Another type of G•U pair interacting via their Watson-Crick edges has been observed in the A site of ribosome structures between a modified U34 in the tRNA anticodon triplet and G + 3 in the mRNA. In such pairs, the electronic structure of the U is changed with a negative charge on N3(U), resulting in two H-bonds between N1(G)…O4(U) and N2(G)…N3(U). Here, we report that such pairs occur in other highly conserved positions in ribosomal RNAs of bacteria in the absence of U modification. An anionic cis Watson-Crick G•G pair is also observed and well conserved in the small subunit. These pairs are observed in tightly folded regions.

Interactions between terminal ribosomal RNA helices stabilize the E. coli large ribosomal subunit.

The ribosome is a large ribonucleoprotein assembly that uses diverse and complex molecular interactions to maintain proper folding. In vivo assembled ribosomes have been isolated using MS2 tags installed in either the 16S or 23S ribosomal RNAs (rRNAs), to enable studies of ribosome structure and function in vitro. RNA tags in the Escherichia coli 50S subunit have commonly been inserted into an extended helix H98 in 23S rRNA, as this addition does not affect cellular growth or in vitro ribosome activity. Here, we find that E. coli 50S subunits with MS2 tags inserted in H98 are destabilized compared to wild-type (WT) 50S subunits. We identify the loss of RNA-RNA tertiary contacts that bridge helices H1, H94, and H98 as the cause of destabilization. Using cryogenic electron microscopy (cryo-EM), we show that this interaction is disrupted by the addition of the MS2 tag and can be restored through the insertion of a single adenosine in the extended H98 helix. This work establishes ways to improve MS2 tags in the 50S subunit that maintain ribosome stability and investigates a complex RNA tertiary structure that may be important for stability in various bacterial ribosomes.

Identifying ribosome heterogeneity using ribosome profiling.

Recent studies have revealed multiple mechanisms that can lead to heterogeneity in ribosomal composition. This heterogeneity can lead to preferential translation of specific panels of mRNAs, and is defined in large part by the ribosomal protein (RP) content, amongst other things. However, it is currently unknown to what extent ribosomal composition is heterogeneous across tissues, which is compounded by a lack of tools available to study it. Here we present dripARF, a method for detecting differential RP incorporation into the ribosome using Ribosome Profiling (Ribo-seq) data. We combine the 'waste' rRNA fragment data generated in Ribo-seq with the known 3D structure of the human ribosome to predict differences in the composition of ribosomes in the material being studied. We have validated this approach using publicly available data, and have revealed a potential role for eS25/RPS25 in development. Our results indicate that ribosome heterogeneity can be detected in Ribo-seq data, providing a new method to study this phenomenon. Furthermore, with dripARF, previously published Ribo-seq data provides a wealth of new information, allowing the identification of RPs of interest in many disease and normal contexts. dripARF is available as part of the ARF R package and can be accessed through https://github.com/fallerlab/ARF.

Initial study and phylogenetic analysis of hard ticks (Acari: Ixodidae) in Nantong, China along the route of avian migration.

The growing concern about migratory birds potentially spreading ticks due to global warming has become a significant issue. The city of Nantong in this study is situated along the East Asia-Australasian Flyway (EAAF), with numerous wetlands serving as roosting sites for migratory birds. We conducted an investigation of hard ticks and determined the phylogenetic characteristics of tick species in this city. We utilized three different genes for our study: the mitochondrial cytochrome oxidase subunit 1 (COX1) gene, the second internal transcribed spacer (ITS2), and the mitochondrial small subunit rRNA (12 S rRNA) gene. The predominant tick species were Haemaphysalis flava (H. flava) and Haemaphysalis longicornis (H. longicornis). Additionally, specimens of Haemaphysalis campanulata (H. campanulata) and Rhipicephalus sanguineus (R. sanguineus) were collected. The H. flava specimens in this study showed a close genetic relationship with those from inland provinces of China, as well as South Korea and Japan. Furthermore, samples of H. longicornis exhibited a close genetic relationship with those from South Korea, Japan, Australia, and the USA, as well as specific provinces in China. Furthermore, R. sanguineus specimens captured in Nantong showed genetic similarities with specimens from Egypt, Nigeria, and Argentina.

Mid-Infrared Photothermal-Fluorescence In Situ Hybridization for Functional Analysis and Genetic Identification of Single Cells.

Simultaneous identification and metabolic analysis of microbes with single-cell resolution and high throughput are necessary to answer the question of "who eats what, when, and where" in complex microbial communities. Here, we present a mid-infrared photothermal-fluorescence in situ hybridization (MIP-FISH) platform that enables direct bridging of genotype and phenotype. Through multiple improvements of MIP imaging, the sensitive detection of isotopically labeled compounds incorporated into proteins of individual bacterial cells became possible, while simultaneous detection of FISH labeling with rRNA-targeted probes enabled the identification of the analyzed cells. In proof-of-concept experiments, we showed that the clear spectral red shift in the protein amide I region due to incorporation of 13C atoms originating from 13C-labeled glucose can be exploited by MIP-FISH to discriminate and identify 13C-labeled bacterial cells within a complex human gut microbiome sample. The presented methods open new opportunities for single-cell structure-function analyses for microbiology.

Ribo-Pop: simple, cost-effective, and widely applicable ribosomal RNA depletion.

The measurement of RNA abundance derived from massively parallel sequencing experiments is an essential technique. Methods that reduce ribosomal RNA levels are usually required prior to sequencing library construction because ribosomal RNA typically comprises the vast majority of a total RNA sample. For some experiments, ribosomal RNA depletion is favored over poly(A) selection because it offers a more inclusive representation of the transcriptome. However, methods to deplete ribosomal RNA are generally proprietary, complex, inefficient, applicable to only specific species, or compatible with only a narrow range of RNA input levels. Here, we describe Ribo-Pop (ribosomal RNA depletion for popular use), a simple workflow and antisense oligo design strategy that we demonstrate works over a wide input range and can be easily adapted to any organism with a sequenced genome. We provide a computational pipeline for probe selection, a streamlined 20-min protocol, and ready-to-use oligo sequences for several organisms. We anticipate that our simple and generalizable "open source" design strategy would enable virtually any laboratory to pursue full transcriptome sequencing in their organism of interest with minimal time and resource investment.

Isomorphic Fluorescent Nucleosides Facilitate Real-Time Monitoring of RNA Depurination by Ribosome Inactivating Proteins.

Ribosome-inactivating proteins, a family of highly cytotoxic proteins, interfere with protein synthesis by depurinating a specific adenosine residue within the conserved α-sarcin/ricin loop of eukaryotic ribosomal RNA. Besides being biological warfare agents, certain RIPs have been promoted as potential therapeutic tools. Monitoring their deglycosylation activity and their inhibition in real time have remained, however, elusive. Herein, we describe the enzymatic preparation and utility of consensus RIP hairpin substrates in which specific G residues, next to the depurination site, are surgically replaced with tz G and th G, fluorescent G analogs. By strategically modifying key positions with responsive fluorescent surrogate nucleotides, RIP-mediated depurination can be monitored in real time by steady-state fluorescence spectroscopy. Subtle differences observed in preferential depurination sites provide insight into the RNA folding as well as RIPs' substrate recognition features.

Molecular phylogeny and new light microscopic data of Metchnikovella spiralis (Microsporidia: Metchnikovellidae), a hyperparasite of eugregarine Polyrhabdina sp. from the polychaete Pygospio elegans.

Metchnikovellids are a deep-branching group of microsporidia, parasites of gregarines inhabiting the alimentary tract of polychaetes and some other invertebrates. The diversity and phylogeny of these hyperparasites remain poorly studied. Modern descriptions and molecular data are still lacking for many species. The results of a light microscopy study and molecular data for Metchnikovella spiralis Sokolova et al., 2014, a hyperparasite of the eugregarine Polyrhabdina sp., isolated from the polychaete Pygospio elegans, were obtained. The original description of M. spiralis was based primarily on the analysis of stained preparations and transmission electron microscopy images. Here, the species description was complemented with the results of in vivo observations and phylogenetic analysis based on the SSU rRNA gene. It was shown that in this species, free sporogony precedes sac-bound sporogony, as it occurs in the life cycle of most other metchnikovellids. Spore sacs are entwined with spirally wound cords, and possess only one polar plug. Phylogenetic analyses did not group M. spiralis with M. incurvata, another metchnikovellid from the same gregarine species, but placed it as a sister branch to Amphiacantha. The paraphyletic nature of the genus Metchnikovella was discussed. The taxonomic summary for M. spiralis was emended.

SMARTer single cell total RNA sequencing.

Single cell RNA sequencing methods have been increasingly used to understand cellular heterogeneity. Nevertheless, most of these methods suffer from one or more limitations, such as focusing only on polyadenylated RNA, sequencing of only the 3' end of the transcript, an exuberant fraction of reads mapping to ribosomal RNA, and the unstranded nature of the sequencing data. Here, we developed a novel single cell strand-specific total RNA library preparation method addressing all the aforementioned shortcomings. Our method was validated on a microfluidics system using three different cancer cell lines undergoing a chemical or genetic perturbation and on two other cancer cell lines sorted in microplates. We demonstrate that our total RNA-seq method detects an equal or higher number of genes compared to classic polyA[+] RNA-seq, including novel and non-polyadenylated genes. The obtained RNA expression patterns also recapitulate the expected biological signal. Inherent to total RNA-seq, our method is also able to detect circular RNAs. Taken together, SMARTer single cell total RNA sequencing is very well suited for any single cell sequencing experiment in which transcript level information is needed beyond polyadenylated genes.

Direct Determination of Pseudouridine in RNA by Mass Spectrometry Coupled with Stable Isotope Labeling.

Pseudouridine (Ψ) is the only "mass-silent" nucleoside produced by post-transcriptional RNA modification. We developed a mass spectrometry (MS)-based technique coupled with in vivo deuterium (D) labeling of uridines for direct determination of Ψs in cellular RNA and applied it to the comprehensive analysis of post-transcriptional modifications in human ribosomal RNAs. The method utilizes human TK6/mouse FM3A cells deficient in uridine monophosphate synthase using a CRISPR-Cas9 technique to turn off de novo uridine synthesis and fully labels uridines with D at uracil positions 5 and 6 by cultivating the cells in a medium containing uridine-5,6-D2. The pseudouridylation reaction in those cells results in the exchange of the D at the C5 of uracil with hydrogen from solvent, which produces a -1 Da mass shift, thus allowing MS-based determination of RNA Ψs. We present here the experimental details of this method and show that it allows the identification of all Ψs in human major nuclear and nucleolar RNAs, including several previously unknown Ψs. Because the method allows direct determination of Ψs at the femtomole level of RNA, it will serve as a useful tool for structure/function studies of a wide variety of noncoding RNAs.

MIST: a Multilocus Identification System for Trichoderma.

Due to the rapid expansion in microbial taxonomy, precise identification of common industrially and agriculturally relevant fungi such as Trichoderma species is challenging. In this study, we introduce the online multilocus identification system (MIST) for automated detection of 349 Trichoderma species based on a set of three DNA barcodes. MIST is based on the reference databases of validated sequences of three commonly used phylogenetic markers collected from public databases. The databases consist of 414 complete sequences of the nuclear rRNA internal transcribed spacers (ITS) 1 and 2, 583 sequence fragments of the gene encoding translation elongation factor 1-alpha (tef1), and 534 sequence fragments of the gene encoding RNA polymerase subunit 2 (rpb2). Through MIST, information from different DNA barcodes can be combined and the identification of Trichoderma species can be achieved based on the integrated parametric sequence similarity search (blastn) performed in the manner of a decision tree classifier. In the verification process, MIST provided correct identification for 44 Trichoderma species based on DNA barcodes consisting of tef1 and rpb2 markers. Thus, MIST can be used to obtain an automated species identification as well as to retrieve sequences required for manual identification by means of phylogenetic analysis.IMPORTANCE The genus Trichoderma is important to humankind, with a wide range of applications in industry, agriculture, and bioremediation. Thus, quick and accurate identification of Trichoderma species is paramount, since it is usually the first step in Trichoderma-based research. However, it frequently becomes a limitation, especially for researchers who lack taxonomic knowledge of fungi. Moreover, as the number of Trichoderma-based studies has increased, a growing number of unidentified sequences have been stored in public databases, which has made the species identification more ambiguous. In this study, we provide an easy-to-use tool, MIST, for automated species identification, a list of Trichoderma species, and corresponding sequences of reference DNA barcodes. Therefore, this study will facilitate the research on the biodiversity and applications of the genus Trichoderma.

Redescription of Rigidohymena inquieta (Stokes, 1887) Berger, 2011 as Metahymena inquieta gen. nov., comb. nov. (Ciliophora, Hypotricha) Based on Morphology, Morphogenesis, and Molecular Phylogeny.

The morphology and morphogenesis of Rigidohymena inquieta (Stokes, 1887) Berger, 2011, isolated from a lawn soil in the campus of the University of Ulsan, Korea, was studied, using live observation and protargol impregnation. The molecular phylogeny was studied based on the SSU rRNA gene sequences. The morphology of the Korean population of R. inquieta matches the previously known populations; however, the morphogenetic pattern shows differences to the species R. candens in the involvement of cirrus V/3 in the anlagen formation. A novel genus namely Metahymena gen. nov. has been erected for the present species based on the ontogenetic difference, and the new combination Metahymena inquieta gen. nov., comb. nov. is proposed. The morphology, morphogenesis, distribution, and phylogeny of M. inquieta are presented. The morphologic and morphogenetic data corroborate the phylogenetic analyses as M. inquieta clusters among the stylonychid ciliates in a clade distant from Rigidohymena candens.

Improved RNA modification mapping of cellular non-coding RNAs using C- and U-specific RNases.

Locating ribonucleoside modifications within an RNA sequence requires digestion of the RNA into oligoribonucleotides of amenable size for subsequent analysis by LC-MS (liquid chromatography-mass spectrometry). This approach, widely referred to as RNA modification mapping, is facilitated through ribonucleases (RNases) such as T1 (guanosine-specific), U2 (purine-selective) and A (pyrimidine-specific) among others. Sequence coverage by these enzymes depends on positioning of the recognized nucleobase (such as guanine or purine or pyrimidine) in the sequence and its ribonucleotide composition. Using E. coli transfer RNA (tRNA) and ribosomal RNA (rRNA) as model samples, we demonstrate the ability of complementary nucleobase-specific ribonucleases cusativin (C-specific) and MC1 (U-specific) to generate digestion products that facilitate confident mapping of modifications in regions such as G-rich and pyrimidine-rich segments of RNA, and to distinguish C to U sequence differences. These enzymes also increase the number of oligonucleotide digestion products that are unique to a specific RNA sequence. Further, with these additional RNases, multiple modifications can be localized with high confidence in a single set of experiments with minimal dependence on the individual tRNA abundance in a mixture. The sequence overlaps observed with these complementary digestion products and that of RNase T1 improved sequence coverage to 75% or above. A similar level of sequence coverage was also observed for the 2904 nt long 23S rRNA indicating their utility has no dependence on RNA size. Wide-scale adoption of these additional modification mapping tools could help expedite the characterization of modified RNA sequences to understand their structural and functional role in various living systems.

A formal redefinition of the genera Nosema and Vairimorpha (Microsporidia: Nosematidae) and reassignment of species based on molecular phylogenetics.

The microsporidian genera Nosema and Vairimorpha comprise a clade described from insects. Currently the genus Nosema is defined as having a dimorphic life cycle characterized by diplokaryotic stages and diplosporoblastic sporogony with two functionally and morphologically distinct spore types ("early" or "primary" and "environmental"). The Vairimorpha life cycle, in addition to a Nosema-type diplokaryotic sporogony, includes an octosporoblastic sporogony producing eight uninucleate spores (octospores) within a sporophorous vesicle. Molecular phylogeny, however, has clearly demonstrated that the genera Nosema and Vairimorpha, characterized by the absence or presence of uninucleate octospores, respectively, represent two polyphyletic taxa, and that octosporogony is turned on and off frequently within taxa, depending on environmental factors such as host species and rearing temperature. In addition, recent studies have shown that both branches of the Vairimorpha-Nosema clade contain species that are uninucleate throughout their life cycle. The SSU rRNA gene sequence data reveal two distinct clades, those closely related to Vairimorpha necatrix, the type species for the genus Vairimorpha, and those closely related to Nosema bombycis, the type species for the genus Nosema. Here, we redefine the two genera, giving priority to molecular character states over those observed at the developmental, structural or ultrastructural levels and present a list of revised species designations. Using this approach, a series of species are renamed (combination novum) and members of two genera, Rugispora and Oligosporidium, are reassigned to Vairimorpha because of their phylogenetic position. Moreover, the family Nosematidae is redefined and includes the genera Nosema and Vairimorpha comprising a monophyletic lineage of Microsporidia.

A novel dimorphic microsporidian Ameson iseebi sp. nov. infecting muscle of the Japanese spiny lobster, Panulirus japonicus, in western Japan.

Japanese spiny lobsters (Panulirus japonicus) exhibiting white opaque abdominal muscle were found in Mie and Wakayama prefectures, in mid-Western Japan. Microscopically, two types of microsporidian spores, ovoid and rod-shaped, were observed infecting the muscle. Histologically, both types of spore were detected inside myofibers of the abdomen, appendages, and cardiac muscles and were often both observed in a single myofiber simultaneously. Transmission electron microscopy revealed that ovoid spores have villous projections on the surface, and that ovoid and rod-shaped spores have a polar filament with 12 coils and 6 to 8 coils respectively. Merogonic and sporogonic stages were observed around ovoid spores, but rarely around rod-shaped spores. The small subunit ribosomal DNA sequences obtained from both spore types were identical to each other, indicating that this microsporidian exhibits a clear spore dimorphism. Phylogenetic analysis based on the rDNA sequences indicates that this microsporidian is part of a clade consisting of the genera Ameson and Nadelspora, with the most closely related species being A. herrnkindi found in the Caribbean spiny lobster P. argus. Based on ultrastructural features, molecular phylogenetic data, host type and geographical differences among known species in these genera, the species found in whitened abdominal muscles of the Japanese spiny lobster is described as Ameson iseebi sp. nov.

Identification and characterization of three microsporidian genera concurrently infecting a silkworm, Bombyx mori, in Brazil.

Microsporidia are important entomopathogens known for infecting insects such as the silkworm (Bombyx mori) thus impairing global silk production. This study aimed to identify and characterize the microsporidia isolated from a diseased larva of silkworm, collected from a sericulture farm in southern Brazil. Identification was performed by phylogenetic analysis of the nucleotide sequences of the SSU rRNA genes. Characterization was performed by analyzing spore sizes, tissue tropism, internal and external symptoms, and pathogenicity against B. mori. Microsporidia belonging to three different genera were identified, namely, Endoreticulatus, Nosema and Tubulinosema. After inoculation of the mixed spores of the microsporidian isolates into B. mori larvae, a high prevalence of Tubulinosema spp. was observed. This isolate showed high prevalence on the silk glands and a late mortality, initially of around 10% until the 20th day post-inoculation but reaching 91.5% upon pupation. Therefore, we demonstrated that Tubulinosema spp. causes chronic infection with slow pathogenicity. We identified for the first time three different microsporidians concurrently infecting B. mori in Brazil. Tubulinosema is of particular interest because of its potential threat to silk production; it affects the formation of silk glands in B. mori while not presenting distinguishable external symptoms or causing the immediate death of these insects. Further studies focusing on this species, mainly regarding its life cycle within the host and the sublethal effects of surviving individuals, demonstrate the importance of describing it as a new species and improving the characterization of the disease in order to prevent its spread.

Molecular characterization of Urosporidium tapetis sp. nov., a haplosporidian hyperparasite infecting metacercariae of Parvatrema duboisi (Dollfus 1923), a trematode parasite of Manila clam Ruditapes philippinarum on the west coast of Korea.

Recently, a putative new hyperparasitic haplosporidian in the genus Urosporidium was identified from metacercariae of the trematode Parvatrema duboisi infecting Manila clam Ruditapes philippinarum on the west coast of Korea. In this study, we applied small subunit ribosomal DNA (SSU rDNA) sequences as a marker to substantiate the phylogenetic relationship of the unidentified Urosporidium within the Order Haplosporida. In our phylogenetic analysis, the 1890 bp of SSU rDNA sequences obtained were closely related to a haplosporidian parasite forming a sister clade to Urosporidium group, although the gene sequences were only 89.22-89.70% similar to Urosporidium spp. Such molecular phylogenetic distance within the genus suggested that the unidentified Urosporidium is a new member of the genus. Accordingly, we report the unidentified haplosporidian hyperparasite as Urosporidium tapetis sp. nov.

Grc3 programs the essential endoribonuclease Las1 for specific RNA cleavage.

Las1 is a recently discovered endoribonuclease that collaborates with Grc3-Rat1-Rai1 to process precursor ribosomal RNA (rRNA), yet its mechanism of action remains unknown. Disruption of the mammalian Las1 gene has been linked to congenital lethal motor neuron disease and X-linked intellectual disability disorders, thus highlighting the necessity to understand Las1 regulation and function. Here, we report that the essential Las1 endoribonuclease requires its binding partner, the polynucleotide kinase Grc3, for specific C2 cleavage. Our results establish that Grc3 drives Las1 endoribonuclease cleavage to its targeted C2 site both in vitro and in Saccharomyces cerevisiae. Moreover, we observed Las1-dependent activation of the Grc3 kinase activity exclusively toward single-stranded RNA. Together, Las1 and Grc3 assemble into a tetrameric complex that is required for competent rRNA processing. The tetrameric Grc3/Las1 cross talk draws unexpected parallels to endoribonucleases RNaseL and Ire1, and establishes Grc3/Las1 as a unique member of the RNaseL/Ire1 RNA splicing family. Together, our work provides mechanistic insight for the regulation of the Las1 endoribonuclease and identifies the tetrameric Grc3/Las1 complex as a unique example of a protein-guided programmable endoribonuclease.

Nucleolus-like compartmentalization of the transcription machinery in fast-growing bacterial cells.

We have learned a great deal about RNA polymerase (RNA Pol), transcription factors, and the transcriptional regulation mechanisms in prokaryotes for specific genes, operons, or transcriptomes. However, we have only begun to understand how the transcription machinery is three-dimensionally (3D) organized into bacterial chromosome territories to orchestrate the transcription process and to maintain harmony with the replication machinery in the cell. Much progress has been made recently in our understanding of the spatial organization of the transcription machinery in fast-growing Escherichia coli cells using state-of-the-art superresolution imaging techniques. Co-imaging of RNA polymerase (RNA Pol) with DNA and transcription elongation factors involved in ribosomal RNA (rRNA) synthesis, and ribosome biogenesis has revealed similarities between bacteria and eukaryotes in the spatial organization of the transcription machinery for growth genes, most of which are rRNA genes. Evidence supports the notion that RNA Pol molecules are concentrated, forming foci at the clustering of rRNA operons resembling the eukaryotic nucleolus. RNA Pol foci are proposed to be active transcription factories for both rRNA genes expression and ribosome biogenesis to support maximal growth in optimal growing conditions. Thus, in fast-growing bacterial cells, RNA Pol foci mimic eukaryotic Pol I activity, and transcription factories resemble nucleolus-like compartmentation. In addition, the transcription and replication machineries are mostly segregated in space to avoid the conflict between the two major cellular functions in fast-growing cells.

Method for Direct Mass-Spectrometry-Based Identification of Monomethylated RNA Nucleoside Positional Isomers and Its Application to the Analysis of Leishmania rRNA.

RNA post-transcriptional modifications are common in all kingdoms of life and are predominantly affiliated with methylations at various nucleobase positions. Methylations occur frequently at specific sites on the RNA nucleobases and appear to regulate site-specific intermolecular/intramolecular interactions. Herein, we present a method that utilizes liquid chromatography-mass spectrometry (LC-MS) to identify positional monomethylated RNA nucleoside isomers. The method produces profiles of in-source fragmentation and subsequent tandem mass spectrometry (MS2) (pseudo-MS3) of RNase-digested fragments of an RNA and distinguishes between positional methylated nucleobase isomers by comparing their intranucleobase fragment ion profiles with signature profiles derived from authentic isomers. For method validation, we independently determined the positions of all known monomethylated nucleoside isomers in the Escherichia coli 16S/23S rRNAs. As proof of concept, we further applied this technology to fully characterize the base-modified nucleoside positional isomers, in rRNAs derived from Leishmania donovani, a human blood parasite afflicting millions around the globe. The method described herein will be highly beneficial for the delineation of RNA modification profiles in various cellular RNAs, and as it only requires a subpicomole amount of RNA, it could also be used for the structure-function studies of RNA populations represented in minute amounts in the cell.