Direct RNA motif definition and identification from multiple sequence alignments using secondary structure profiles

RNAMotifProfile: a graph-based approach to build RNA structural motif profiles

RNA structural motifs are the recurrent segments in RNA three-dimensional structures that play a crucial role in the functional diversity of RNAs. Understanding the similarities and variations within these recurrent motif groups is essential for gaining insights into RNA structure and function. While recurrent structural motifs are generally assumed to be composed of the same isosteric base interactions, this consistent pattern is not observed across all examples of these motifs. Existing methods for analyzing and comparing RNA structural motifs may overlook variations in base interactions and associated nucleotides. RNAMotifProfile is a novel profile-to-profile alignment algorithm that generates a comprehensive profile from a group of structural motifs, incorporating all base interactions and associated nucleotides at each position. By structurally aligning input motif instances using a guide-tree-based approach, RNAMotifProfile captures the similarities and variations within recurrent motif groups. Additionally, RNAMotifProfile can function as a motif search tool, enabling the identification of instances of a specific motif family by searching with the corresponding profile. The ability to generate accurate and comprehensive profiles for RNA structural motif families, and to search for these motifs, facilitates a deeper understanding of RNA structure-function relationships and potential applications in RNA engineering and therapeutic design.

Reliable Genomic Integration Sites in Pseudomonas putida Identified by Two-Dimensional Transcriptome Analysis

Genomic integration is commonly used to engineer stable production hosts. However, so far, for many microbial workhorses, only a few integration sites have been characterized, thereby restraining advanced strain engineering that requires multiple insertions. Here, we report on the identification of novel genomic integration sites, so-called landing pads, for Pseudomonas putida KT2440. We identified genomic regions with constant expression patterns under diverse experimental conditions by using RNA-Seq data. Homologous recombination constructs were designed to insert heterologous genes into intergenic sites in these regions, allowing condition-independent gene expression. Ten potential landing pads were characterized using four different msfGFP expression cassettes. An insulated probe sensor was used to study locus-dependent effects on recombinant gene expression, excluding genomic read-through of flanking promoters under changing cultivation conditions. While the reproducibility of expression in the landing pads was very high, the msfGFP signals varied strongly between the different landing pads, confirming a strong influence of the genomic context. To showcase that the identified landing pads are also suitable candidates for heterologous gene expression in other Pseudomonads, four equivalent landing pads were identified and characterized in Pseudomonas taiwanensis VLB120. This study shows that genomic "hot" and "cold" spots exist, causing strong promoter-independent variations in gene expression. This highlights that the genomic context is an additional parameter to consider when designing integrable genomic cassettes for tailored heterologous expression. The set of characterized genomic landing pads presented here further increases the genetic toolbox for deep metabolic engineering in Pseudomonads.

Decoding the Non-coding: Tools and Databases Unveiling the Hidden World of "Junk" RNAs for Innovative Therapeutic Exploration

Non-coding RNAs are pivotal regulators of gene and protein expression, exerting crucial influences on diverse biological processes. Their dysregulation is frequently implicated in the onset and progression of diseases, notably cancer. A profound comprehension of the intricate mechanisms governing ncRNAs is imperative for devising innovative therapeutic interventions against these debilitating conditions. Significantly, nearly 80% of our genome comprises ncRNAs, underscoring their centrality in cellular processes. The elucidation of ncRNA functions is pivotal for grasping the complexities of gene regulation and its implications for human health. Modern genome sequencing techniques yield vast datasets, stored in specialized databases. To harness this wealth of information and to understand the crosstalk of non-coding RNAs, knowledge of available databases is required, and many new sophisticated computational tools have emerged. These tools play a pivotal role in the identification, prediction, and annotation of ncRNAs, thereby facilitating their experimental validation. This Review succinctly outlines the current understanding of ncRNAs, emphasizing their involvement in disease development. It also highlights the databases and tools instrumental in classifying, annotating, and evaluating ncRNAs. By extracting meaningful biological insights from seemingly "junk" data, these tools empower scientists to unravel the intricate roles of ncRNAs in shaping human health.

CA, Harmon B, Espinola JA, Reck Cechinel L, Batabyal R, Freishtat RJ, Hahn A. Bacterial small RNAs may mediate immune response differences seen in respiratory syncytial virus versus rhinovirus bronchiolitis

Bronchiolitis, a viral lower respiratory infection, is the leading cause of infant hospitalization, which is associated with an increased risk for developing asthma later in life. Bronchiolitis can be caused by several respiratory viruses, such as respiratory syncytial virus (RSV), rhinovirus (RV), and others. It can also be caused by a solo infection (e.g., RSV- or RV-only bronchiolitis) or co-infection with two or more viruses. Studies have shown viral etiology-related differences between RSV- and RV-only bronchiolitis in the immune response, human microRNA (miRNA) profiles, and dominance of certain airway microbiome constituents. Here, we identified bacterial small RNAs (sRNAs), the prokaryotic equivalent to eukaryotic miRNAs, that differ between infants of the 35th Multicenter Airway Research Collaboration (MARC-35) cohort with RSV- versus RV-only bronchiolitis. We first derived reference sRNA datasets from cultures of four bacteria known to be associated with bronchiolitis (i.e., Haemophilus influenzae, Moraxella catarrhalis, Moraxella nonliquefaciens, and Streptococcus pneumoniae). Using these reference sRNA datasets, we found several sRNAs associated with RSV- and RV-only bronchiolitis in our human nasal RNA-Seq MARC-35 data. We also determined potential human transcript targets of the bacterial sRNAs and compared expression of the sRNAs between RSV- and RV-only cases. sRNAs are known to downregulate their mRNA target, we found that, compared to those associated with RV-only bronchiolitis, sRNAs associated with RSV-only bronchiolitis may relatively activate the IL-6 and IL-8 pathways and relatively inhibit the IL-17A pathway. These data support that bacteria may be contributing to inflammation differences seen in RSV- and RV-only bronchiolitis, and for the first time indicate that the potential mechanism in doing so may be through bacterial sRNAs.

Comparative RNA Genomics

Over the last quarter of a century it has become clear that RNA is much more than just a boring intermediate in protein expression. Ancient RNAs still appear in the core information metabolism and comprise a surprisingly large component in bacterial gene regulation. A common theme with these types of mostly small RNAs is their reliance of conserved secondary structures. Large-scale sequencing projects, on the other hand, have profoundly changed our understanding of eukaryotic genomes. Pervasively transcribed, they give rise to a plethora of large and evolutionarily extremely flexible non-coding RNAs that exert a vastly diverse array of molecule functions. In this chapter we provide a-necessarily incomplete-overview of the current state of comparative analysis of non-coding RNAs, emphasizing computational approaches as a means to gain a global picture of the modern RNA world.

Antagonistic interactions between phage and host factors control arbitrium lysis-lysogeny decision

Phages can use a small-molecule communication arbitrium system to coordinate lysis-lysogeny decisions, but the underlying mechanism remains unknown. Here we determined that the arbitrium system in Bacillus subtilis phage phi3T modulates the bacterial toxin-antitoxin system MazE-MazF to regulate the phage life cycle. We show that phi3T expresses AimX and YosL, which bind to and inactivate MazF. AimX also inhibits the function of phi3T_93, a protein that promotes lysogeny by binding to MazE and releasing MazF. Overall, these mutually exclusive interactions promote the lytic cycle of the phage. After several rounds of infection, the phage-encoded AimP peptide accumulates intracellularly and inactivates the phage antiterminator AimR, a process that eliminates aimX expression from the aimP promoter. Therefore, when AimP increases, MazF activity promotes reversion back to lysogeny, since AimX is absent. Altogether, our study reveals the evolutionary strategy used by arbitrium to control lysis-lysogeny by domesticating and fine-tuning a phage-defence mechanism.

INTERPIN: A repository for intrinsic transcription termination hairpins in bacteria

The large-scale detection of putative intrinsic transcription terminators is limited to only a few bacteria currently. We discovered a group of hairpins, called cluster hairpins, present within 15 nucleotides from each other. These are expected to work in tandem to cause intrinsic transcription termination (ITT), while the single hairpin can do the same alone. Therefore, exploring these ITT sites and the hairpins across bacterial genomes becomes highly desirable. INTERPIN is the largest archived collection of in silico inferred ITT hairpins in bacteria, covering 12745 bacterial genomes and encompassing ten bacterial phyla for ∼25 million hairpins. Users can obtain details on operons, individual cluster, and single ITT hairpins that were screened therein. Integrated Genome Viewer (IGV) software interactively visualizes hairpin secondary and tertiary structures in the genomic context. We also discuss statistics for the occurrence of cluster or single hairpins and other termination alternatives while showing the validation of predicted hairpins against in vivo detected hairpins. The database is freely available at http://pallab.cds.iisc.ac.in/INTERPIN/. INTERPIN (database and software) can make predictions for both AT and GC-rich genomes, which has not been achieved by any other program so far. It can also be used to improve genome annotation as well as to get predictions to improve the understanding of the ITT pathway by further analysis.

Three genes controlling streptomycin susceptibility in Agrobacterium fabrum

Streptomycin (Sm) is a commonly used antibiotic for its efficacy against diverse bacteria. The plant pathogen Agrobacterium fabrum is a model for studying pathogenesis and interkingdom gene transfer. Streptomycin-resistant variants of A. fabrum are commonly employed in genetic analyses, yet mechanisms of resistance and susceptibility to streptomycin in this organism have not previously been investigated. We observe that resistance to a high concentration of streptomycin arises at high frequency in A. fabrum, and we attribute this trait to the presence of a chromosomal gene (strB) encoding a putative aminoglycoside phosphotransferase. We show how strB, along with rpsL (encoding ribosomal protein S12) and rsmG (encoding a 16S rRNA methyltransferase), modulates streptomycin sensitivity in A. fabrum. IMPORTANCE The plant pathogen Agrobacterium fabrum is a widely used model bacterium for studying biofilms, bacterial motility, pathogenesis, and gene transfer from bacteria to plants. Streptomycin (Sm) is an aminoglycoside antibiotic known for its broad efficacy against gram-negative bacteria. A. fabrum exhibits endogenous resistance to somewhat high levels of streptomycin, but the mechanism underlying this resistance has not been elucidated. Here, we demonstrate that this resistance is caused by a chromosomally encoded streptomycin-inactivating enzyme, StrB, that has not been previously characterized in A. fabrum. Furthermore, we show how the genes rsmG, rpsL, and strB jointly modulate streptomycin susceptibility in A. fabrum.

Small regulatory RNAs are mediators of the Streptococcus mutans SloR regulon

Dental caries is among the most prevalent chronic diseases worldwide. Streptococcus mutans, the chief causative agent of caries, uses a 25-kDa manganese-dependent SloR protein to coordinate the uptake of essential manganese with the transcription of its virulence attributes. Small non-coding RNAs (sRNAs) can either enhance or repress gene expression, and reports in the literature ascribe an emerging role for sRNAs in the environmental stress response. Herein, we focused our attention on 18-50 nt sRNAs as mediators of the S. mutans SloR and manganese regulons. Specifically, the results of RNA sequencing revealed 19 sRNAs in S. mutans, which were differentially transcribed in the SloR-proficient UA159 and SloR-deficient GMS584 strains, and 10 sRNAs that were differentially expressed in UA159 cells grown in the presence of low vs high manganese. We describe SmsR1532 and SmsR1785 as SloR- and manganese-responsive sRNAs that are processed from large transcripts and that bind SloR directly in their promoter regions. The predicted targets of these sRNAs include regulators of metal ion transport, growth management via a toxin-antitoxin operon, and oxidative stress tolerance. These findings support a role for sRNAs in coordinating intracellular metal ion homeostasis with virulence gene control in an important oral cariogen. IMPORTANCE Small regulatory RNAs (sRNAs) are critical mediators of environmental signaling, particularly in bacterial cells under stress, but their role in Streptococcus mutans is poorly understood. S. mutans, the principal causative agent of dental caries, uses a 25-kDa manganese-dependent protein, called SloR, to coordinate the regulated uptake of essential metal ions with the transcription of its virulence genes. In the present study, we identified and characterized sRNAs that are both SloR and manganese responsive. Taken together, this research can elucidate the details of regulatory networks that engage sRNAs in an important oral pathogen and that can enable the development of an effective anti-caries therapeutic.

Extraordinary long-stem confers resistance of intrinsic terminators to processive antitermination

Many prokaryotic operons encode a processive antitermination (P-AT) system. Transcription complexes associated with an antitermination factor can bypass multiple transcription termination signals regardless of their sequences. However, to avoid compromising transcriptional regulation of downstream regions, the terminator at the end of the operon needs to be resistant to antitermination. So far, no studies on the mechanism of resistance to antitermination have been reported. The recently discovered conAn P-AT system is composed of two components that are encoded at the start of many conjugation operons on plasmids of Gram-positive bacteria. Here we report the identification of a conAn-resistant terminator, named TerR, in the conjugation operon of the Bacillus subtilis plasmid pLS20, re-defining the end of the conjugation operon. We investigated the various characteristics of TerR and show that its extraordinary long stem is the determining feature for resistance to antitermination. This is the first P-AT resistance mechanism to be reported.

Mitochondrial genome annotation with MFannot: a critical analysis of gene identification and gene model prediction

Compared to nuclear genomes, mitochondrial genomes (mitogenomes) are small and usually code for only a few dozen genes. Still, identifying genes and their structure can be challenging and time-consuming. Even automated tools for mitochondrial genome annotation often require manual analysis and curation by skilled experts. The most difficult steps are (i) the structural modelling of intron-containing genes; (ii) the identification and delineation of Group I and II introns; and (iii) the identification of moderately conserved, non-coding RNA (ncRNA) genes specifying 5S rRNAs, tmRNAs and RNase P RNAs. Additional challenges arise through genetic code evolution which can redefine the translational identity of both start and stop codons, thus obscuring protein-coding genes. Further, RNA editing can render gene identification difficult, if not impossible, without additional RNA sequence data. Current automated mito- and plastid-genome annotators are limited as they are typically tailored to specific eukaryotic groups. The MFannot annotator we developed is unique in its applicability to a broad taxonomic scope, its accuracy in gene model inference, and its capabilities in intron identification and classification. The pipeline leverages curated profile Hidden Markov Models (HMMs), covariance (CMs) and ERPIN models to better capture evolutionarily conserved signatures in the primary sequence (HMMs and CMs) as well as secondary structure (CMs and ERPIN). Here we formally describe MFannot, which has been available as a web-accessible service (https://megasun.bch.umontreal.ca/apps/mfannot/) to the research community for nearly 16 years. Further, we report its performance on particularly intron-rich mitogenomes and describe ongoing and future developments.

Small regulatory RNAs are mediators of the Streptococcus mutans SloR regulon

Dental caries is among the most prevalent chronic infectious diseases worldwide. Streptococcus mutans , the chief causative agent of caries, uses a 25 kDa manganese dependent SloR protein to coordinate the uptake of essential manganese with the transcription of its virulence attributes. Small non-coding RNAs (sRNAs) can either enhance or repress gene expression and reports in the literature ascribe an emerging role for sRNAs in the environmental stress response. Herein, we identify 18-50 nt sRNAs as mediators of the S. mutans SloR and manganese regulons. Specifically, the results of sRNA-seq revealed 56 sRNAs in S. mutans that were differentially transcribed in the SloR-proficient UA159 and SloR-deficient GMS584 strains, and 109 sRNAs that were differentially expressed in UA159 cells grown in the presence of low versus high manganese. We describe SmsR1532 and SmsR1785 as SloR- and/or manganese-responsive sRNAs that are processed from large transcripts, and that bind SloR directly in their promoter regions. The predicted targets of these sRNAs include regulators of metal ion transport, growth management via a toxin-antitoxin operon, and oxidative stress tolerance. These findings support a role for sRNAs in coordinating intracellular metal ion homeostasis with virulence gene control in an important oral cariogen.

IMPORTANCE

Small regulatory RNAs (sRNAs) are critical mediators of environmental signaling, particularly in bacterial cells under stress, but their role in Streptococcus mutans is poorly understood. S. mutans, the principal causative agent of dental caries, uses a 25 kDa manganese-dependent protein, called SloR, to coordinate the regulated uptake of essential metal ions with the transcription of its virulence genes. In the present study, we identified and characterize sRNAs that are both SloR- and manganese-responsive. Taken together, this research can elucidate the details of regulatory networks that engage sRNAs in an important oral pathogen, and that can enable the development of an effective anti-caries therapeutic.

Thirteen dubious ways to detect conserved structural RNAs

Covariation induced by compensatory base substitutions in RNA alignments is a great way to deduce conserved RNA structure, in principle. In practice, success depends on many factors, importantly the quality and depth of the alignment and the choice of covariation statistic. Measuring covariation between pairs of aligned positions is easy. However, using covariation to infer evolutionarily conserved RNA structure is complicated by other extraneous sources of covariation such as that resulting from homologous sequences having evolved from a common ancestor. In order to provide evidence of evolutionarily conserved RNA structure, a method to distinguish covariation due to sources other than RNA structure is necessary. Moreover, there are several sorts of artifactually generated covariation signals that can further confound the analysis. Additionally, some covariation signal is difficult to detect due to incomplete comparative data. Here, we investigate and critically discuss the practice of inferring conserved RNA structure by comparative sequence analysis. We provide new methods on how to approach and decide which of the numerous long non-coding RNAs (lncRNAs) have biologically relevant structures.

Identification and initial characterization of Hfq-associated sRNAs in Histophilus somni strain 2336

Small RNAs (sRNA), in association with the global chaperone regulator Hfq, positively or negatively regulate gene expression in bacteria. For this study, Histophilus somni sRNAs that bind to Hfq were identified and then partially characterized. The Hfq-associated sRNAs in H. somni were isolated and identified by co-immunoprecipitation using anti-Hfq antibody, followed by sRNA sequencing. Sequence analysis of the sRNA samples identified 100 putative sRNAs, out of which 16 were present in pathogenic strain 2336, but not in non-pathogenic strain 129Pt. Bioinformatic analyses suggested that the sRNAs HS9, HS79, and HS97 could bind to many genes putatively involved in virulence/biofilm formation. Furthermore, multi-sequence alignment of the sRNA regions in the genome revealed that HS9 and HS97 could interact with sigma 54, which is a transcription factor linked to important bacterial traits, including motility, virulence, and biofilm formation. Northern blotting was used to determine the approximate size, abundance and any processing events attributed to the sRNAs. Selected sRNA candidates were confirmed to bind Hfq, as determined by electrophoretic mobility shift assays using sRNAs synthesized by in vitro transcription and recombinant Hfq. The exact transcriptional start site of the sRNA candidates was determined by RNA ligase-mediated rapid amplification of cDNA ends, followed by cloning and sequencing. This is the first investigation of H. somni sRNAs that show they may have important regulatory roles in virulence and biofilm formation.

Removal of AMR plasmids using a mobile, broad host-range CRISPR-Cas9 delivery tool

Antimicrobial resistance (AMR) genes are widely disseminated on plasmids. Therefore, interventions aimed at blocking plasmid uptake and transfer may curb the spread of AMR. Previous studies have used CRISPR-Cas-based technology to remove plasmids encoding AMR genes from target bacteria, using either phage- or plasmid-based delivery vehicles that typically have narrow host ranges. To make this technology feasible for removal of AMR plasmids from multiple members of complex microbial communities, an efficient, broad host-range delivery vehicle is needed. We engineered the broad host-range IncP1-plasmid pKJK5 to encode cas9 programmed to target an AMR gene. We demonstrate that the resulting plasmid pKJK5::csg has the ability to block the uptake of AMR plasmids and to remove resident plasmids from Escherichia coli. Furthermore, due to its broad host range, pKJK5::csg successfully blocked AMR plasmid uptake in a range of environmental, pig- and human-associated coliform isolates, as well as in isolates of two species of Pseudomonas. This study firmly establishes pKJK5::csg as a promising broad host-range CRISPR-Cas9 delivery tool for AMR plasmid removal, which has the potential to be applied in complex microbial communities to remove AMR genes from a broad range of bacterial species.

Detection of Streptococcus pyogenes M1UK in Australia and characterization of the mutation driving enhanced expression of superantigen SpeA

A new variant of Streptococcus pyogenes serotype M1 (designated 'M1_UK') has been reported in the United Kingdom, linked with seasonal scarlet fever surges, marked increase in invasive infections, and exhibiting enhanced expression of the superantigen SpeA. The progenitor S. pyogenes 'M1_global' and M1_UK clones can be differentiated by 27 SNPs and 4 indels, yet the mechanism for speA upregulation is unknown. Here we investigate the previously unappreciated expansion of M1_UK in Australia, now isolated from the majority of serious infections caused by serotype M1 S. pyogenes. M1_UK sub-lineages circulating in Australia also contain a novel toxin repertoire associated with epidemic scarlet fever causing S. pyogenes in Asia. A single SNP in the 5' transcriptional leader sequence of the transfer-messenger RNA gene ssrA drives enhanced SpeA superantigen expression as a result of ssrA terminator read-through in the M1_UK lineage. This represents a previously unappreciated mechanism of toxin expression and urges enhanced international surveillance.

Early antitermination in the atypical coliphage mEp021 mediated by the Gp17 protein

The coliphage mEp021 belongs to a phage group with a unique immunity repressor, and its life cycle requires the host factor Nus. mEp021 has been classified as non-lambdoid based on its specific characteristics. The mEp021 genome carries a gene encoding an N_λ-like antiterminator protein, termed Gp17, and three nut sites (nutL, nutR1, and nutR2). Analysis of plasmid constructs containing these nut sites, a transcription terminator, and a GFP reporter gene showed high levels of fluorescence when Gp17 was expressed, but not in its absence. Like lambdoid N proteins, Gp17 has an arginine-rich motif (ARM), and mutations in its arginine codons inhibit its function. In infection assays using the mutant phage mEp021ΔGp17::Kan (where gp17 has been deleted), gene transcripts located downstream of transcription terminators were obtained only when Gp17 was expressed. In contrast to phage lambda, mEp021 virus particle production was partially restored (>1/3 relative to wild type) when nus mutants (nusA1, nusB5, nusC60, and nusE71) were infected with mEp021 and Gp17 was overexpressed. Our results suggest that RNA polymerase reads through the third nut site (nutR2), which is more than 7.9 kbp downstream of nutR1.

Non-coding RNAs in human health and disease: potential function as biomarkers and therapeutic targets

Human diseases have been a critical threat from the beginning of human history. Knowing the origin, course of action and treatment of any disease state is essential. A microscopic approach to the molecular field is a more coherent and accurate way to explore the mechanism, progression, and therapy with the introduction and evolution of technology than a macroscopic approach. Non-coding RNAs (ncRNAs) play increasingly important roles in detecting, developing, and treating all abnormalities related to physiology, pathology, genetics, epigenetics, cancer, and developmental diseases. Noncoding RNAs are becoming increasingly crucial as powerful, multipurpose regulators of all biological processes. Parallel to this, a rising amount of scientific information has revealed links between abnormal noncoding RNA expression and human disorders. Numerous non-coding transcripts with unknown functions have been found in addition to advancements in RNA-sequencing methods. Non-coding linear RNAs come in a variety of forms, including circular RNAs with a continuous closed loop (circRNA), long non-coding RNAs (lncRNA), and microRNAs (miRNA). This comprises specific information on their biogenesis, mode of action, physiological function, and significance concerning disease (such as cancer or cardiovascular diseases and others). This study review focuses on non-coding RNA as specific biomarkers and novel therapeutic targets.

Characterization and genome analysis of Escherichia phage fBC-Eco01, isolated from wastewater in Tunisia

The rise of antibiotic resistance in bacterial strains has led to vigorous exploration for alternative treatments. To this end, phage therapy has been revisited, and it is gaining increasing attention, as it may represent an efficient alternative for treating multiresistant pathogenic bacteria. Phage therapy is considered safe, and phages do not infect eukaryotic cells. There have been many studies investigating phage-host bacteria interactions and the ability of phages to target specific hosts. Escherichia coli is the causative agent of a multitude of infections, ranging from urinary tract infections to sepsis, with growing antibiotic resistance. In this study, we characterized the Escherichia phage fBC-Eco01, which was isolated from a water sample collected at Oued, Tunis. Electron microscopy showed that fBC-Eco01 phage particles have siphovirus morphology, with an icosahedral head of 61 ± 3 nm in diameter and a non-contractile tail of 94 ± 2 nm in length and 12 ± 0.9 nm in width. The genome of fBC-Eco01 is a linear double-stranded DNA of 43.466 bp with a GC content of 50.4%. Comparison to databases allowed annotation of the functions to 39 of the 78 predicted gene products. A single-step growth curve revealed that fBC-Eco01 has a latent period of 30 minutes and a burst size of 175 plaque-forming units (PFU) per infected cell. Genomic analysis indicated that fBC-Eco01 is a member of the subfamily Guernseyvirinae. It is most closely related to a group of phages of the genus Kagunavirus that infect Enterobacter, Raoultella, and Escherichia strains.

Bioinformatic and experimental methods to identify and validate bacterial RNA-human RNA interactions

Ample evidence supports the importance of the microbiota on human health and disease. Recent studies suggest that extracellular vesicles are an important means of bacterial-host communication, in part via the transport of small RNAs (sRNAs). Bacterial sRNAs have been shown to co-precipitate with human and mouse RNA-induced silencing complex, hinting that some may regulate gene expression as eukaryotic microRNAs do. Bioinformatic tools, including those that can incorporate an sRNA's secondary structure, can be used to predict interactions between bacterial sRNAs and human messenger RNAs (mRNAs). Validation of these potential interactions using reproducible experimental methods is essential to move the field forward. This review will cover the evidence of interspecies communication via sRNAs, bioinformatic tools currently available to identify potential bacterial sRNA-host (specifically, human) mRNA interactions, and experimental methods to identify and validate those interactions.

Polar flagellar wrapping and lateral flagella jointly contribute to Shewanella putrefaciens environmental spreading

Flagella enable bacteria to actively spread within the environment. A number of species possess two separate flagellar systems, where in most cases a primary polar flagellar system is supported by distinct secondary lateral flagella under appropriate conditions. Using functional fluorescence tagging on one of these species, Shewanella putrefaciens, as a model system, we explored how two different flagellar systems can exhibit efficient joint function. The S. putrefaciens secondary flagellar filaments are composed as a mixture of two highly homologous non-glycosylated flagellins, FlaA₂ and FlaB₂ . Both are solely sufficient to form a functional filament, however, full spreading motility through soft agar requires both flagellins. During swimming, lateral flagella emerge from the cell surface at angles between 30° and 50°, and only filaments located close to the cell pole may form a bundle. Upon a directional shift from forward to backward swimming initiated by the main polar flagellum, the secondary filaments flip over and thus support propulsion into either direction. Lateral flagella do not inhibit the wrapping of the polar flagellum around the cell body at high load. Accordingly, screw thread-like motility mediated by the primary flagellum and activity of lateral flagella cumulatively supports spreading through constricted environments such as polysaccharide matrices.

Deciphering the mitochondrial genome of Hemerocallis citrina (Asphodelaceae) using a combined assembly and comparative genomic strategy

Hemerocallis citrina is a perennial herbaceous plant that is dedicated to mothers in Chinese culture and is widely distributed across the country. As a popular species with a long history of cultivation and utilization, it is renowned for its remarkable edible and medicinal value. In this study, we integrated Illumina short-read and Oxford Nanopore long-read sequencing to generate a complete mitochondrial genome (mitogenome) assembly of H. citrina. The H. citrina mitogenome has a multiple chromosomal structure consisting of three circular molecules that are 45,607 bp, 239,991 bp, and 182,864 bp long. We correspondingly annotated 66 genes, comprising 45 protein-coding genes (PCGs), 17 tRNA genes, and 4 rRNA genes. Comparative analysis of gene organization indicated that six syntenic gene clusters were conserved in the mitogenomes of the compared plants. The investigation of repeat content revealed repeat-rich nature of the H. citrina mitogenome, for which plentiful dispersed repeats were characterized to correlate with the size of the mitogenome. The codon usage behavior disclosed that Leucine (Leu) and Serine (Ser) were the most preferred amino acids in H. citrina, and nearly all of the codons with relative synonymous codon usage (RSCU) values greater than 1 showed the preference of A or T ending. Moreover, we inferred a total of 679 RNA editing sites in all mitochondrial PCGs, which presented perfect C-to-U types and tended to lead to the alteration of internal codons. Subsequent selective pressure analysis showed that the majority of the PCGs had undergone evolutionary negative selections, with atp9 in particular undergoing strong stabilizing selection, reflecting its indispensable function in mitogenomes. According to the phylogenetic analysis, H. citrina is close to the species Allium cepa (Amaryllidaceae) and Asparagus officinalis (Asparagaceae) in evolutionary terms. Overall, this project presents the first complete mitogenome of H. citrina, which could provide a reference genome for the comprehensive exploration of the Asphodelaceae family and can facilitate further genomic breeding and evolutionary research on this medicine-food homologous plant.

Discovery and Heterologous Expression of Microginins from Microcystis aeruginosa LEGE 91341

Microginins are a large family of cyanobacterial lipopeptide protease inhibitors. A hybrid polyketide synthase/non-ribosomal peptide synthetase biosynthetic gene cluster (BGC) found in several microginin-producing strains─mic─was proposed to encode the production of microginins, based on bioinformatic analysis. Here, we explored a cyanobacterium, Microcystis aeruginosa LEGE 91341, which contains a mic BGC, to discover 12 new microginin variants. The new compounds contain uncommon amino acids, namely, homophenylalanine (Hphe), homotyrosine (Htyr), or methylproline, as well as a 3-aminodecanoic acid (Ada) residue, which in some variants was chlorinated at its terminal methyl group. We have used direct pathway cloning (DiPaC) to heterologously express the mic BGC from M. aeruginosa LEGE 91341 in Escherichia coli, which led to the production of several microginins. This proved that the mic BGC is, in fact, responsible for the biosynthesis of microginins and paves the way to accessing new variants from (meta)genome data or through pathway engineering.

CRISPRtracrRNA: robust approach for CRISPR tracrRNA detection

MOTIVATION

The CRISPR-Cas9 system is a Type II CRISPR system that has rapidly become the most versatile and widespread tool for genome engineering. It consists of two components, the Cas9 effector protein, and a single guide RNA that combines the spacer (for identifying the target) with the tracrRNA, a trans-activating small RNA required for both crRNA maturation and interference. While there are well-established methods for screening Cas effector proteins and CRISPR arrays, the detection of tracrRNA remains the bottleneck in detecting Class 2 CRISPR systems.

RESULTS

We introduce a new pipeline CRISPRtracrRNA for screening and evaluation of tracrRNA candidates in genomes. This pipeline combines evidence from different components of the Cas9-sgRNA complex. The core is a newly developed structural model via covariance models from a sequence-structure alignment of experimentally validated tracrRNAs. As additional evidence, we determine the terminator signal (required for the tracrRNA transcription) and the RNA-RNA interaction between the CRISPR array repeat and the 5'-part of the tracrRNA. Repeats are detected via an ML-based approach (CRISPRidenify). Providing further evidence, we detect the cassette containing the Cas9 (Type II CRISPR systems) and Cas12 (Type V CRISPR systems) effector protein. Our tool is the first for detecting tracrRNA for Type V systems.

AVAILABILITY AND IMPLEMENTATION

The implementation of the CRISPRtracrRNA is available on GitHub upon requesting the access permission, (https://github.com/BackofenLab/CRISPRtracrRNA). Data generated in this study can be obtained upon request to the corresponding person: Rolf Backofen (backofen@informatik.uni-freiburg.de).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Key players in regulatory RNA realm of bacteria

Precise regulation of gene expression is crucial for living cells to adapt for survival in diverse environmental conditions. Among the common cellular regulatory mechanisms, RNA-based regulators play a key role in all domains of life. Discovery of regulatory RNAs have made a paradigm shift in molecular biology as many regulatory functions of RNA have been identified beyond its canonical roles as messenger, ribosomal and transfer RNA. In the complex regulatory RNA network, riboswitches, small RNAs, and RNA thermometers can be identified as some of the key players. Herein, we review the discovery, mechanism, and potential therapeutic use of these classes of regulatory RNAs mainly found in bacteria. Being highly adaptive organisms that inhabit a broad range of ecological niches, bacteria have adopted tight and rapid-responding gene regulation mechanisms. This review aims to highlight how bacteria utilize versatile RNA structures and sequences to build a sophisticated gene regulation network.

•

The three major classes of prokaryotic ncRNAs and their characterized mechanisms of operation in gene regulation.

•

sRNAs emerging as major players in global gene regulatory networks.

•

Riboswitch mediated gene control mechanisms through on/off switches in response to ligand binding.

•

RNA thermo sensors for temperature-dependent gene expression.

•

Therapeutic importance of ncRNAs and computational approaches involved in the discovery of ncRNAs.

The Isolation and Characterization of a Broad Host Range Bcep22-like Podovirus JC1

Bacteriophage JC1 is a Podoviridae phage with a C1 morphotype, isolated on host strain Burkholderia cenocepacia Van1. Phage JC1 is capable of infecting an expansive range of Burkholderia cepacia complex (Bcc) species. The JC1 genome exhibits significant similarity and synteny to Bcep22-like phages and to many Ralstonia phages. The genome of JC1 was determined to be 61,182 bp in length with a 65.4% G + C content and is predicted to encode 76 proteins and 1 tRNA gene. Unlike the other Lessieviruses, JC1 encodes a putative helicase gene in its replication module, and it is in a unique organization not found in previously analyzed phages. The JC1 genome also harbours 3 interesting moron genes, that encode a carbon storage regulator (CsrA), an N-acetyltransferase, and a phosphoadenosine phosphosulfate (PAPS) reductase. JC1 can stably lysogenize its host Van1 and integrates into the 5' end of the gene rimO. This is the first account of stable integration identified for Bcep22-like phages. JC1 has a higher global virulence index at 37 °C than at 30 °C (0.8 and 0.21, respectively); however, infection efficiency and lysogen stability are not affected by a change in temperature, and no observable temperature-sensitive switch between lytic and lysogenic lifestyle appears to exist. Although JC1 can stably lysogenize its host, it possesses some desirable characteristics for use in phage therapy. Phage JC1 has a broad host range and requires the inner core of the bacterial LPS for infection. Bacteria that mutate to evade infection by JC1 may develop a fitness disadvantage as seen in previously characterized LPS mutants lacking inner core.

Crucial role of the chaperonin GroES/EL for heterologous production of the soluble methane monooxygenase from Methylomonas methanica MC09

Methane is a widespread energy source and can serve as an attractive C1 building block for a future bioeconomy. The soluble methane monooxygenase (sMMO) is able to break the strong C−H bond of methane and convert it to methanol. The high structural complexity, multiplex cofactors, and unfamiliar folding or maturation procedures of sMMO have hampered the heterologous production and thus biotechnological applications. Here, we demonstrate the heterologous production of active sMMO from the marine Methylomonas methanica MC09 in Escherichia coli by co‐synthesizing the GroES/EL chaperonin. Iron determination, electron paramagnetic resonance spectroscopy, and native gel immunoblots revealed the incorporation of the non‐heme diiron centre and homodimer formation of active sMMO. The production of recombinant sMMO will enable the expansion of the possibilities of detailed studies, allowing for a variety of novel biotechnological applications.

Refining Mitochondrial Intron Classification With ERPIN: Identification Based on Conservation of Sequence Plus Secondary Structure Motifs

Mitochondrial genomes—in particular those of fungi—often encode genes with a large number of Group I and Group II introns that are conserved at both the sequence and the RNA structure level. They provide a rich resource for the investigation of intron and gene structure, self- and protein-guided splicing mechanisms, and intron evolution. Yet, the degree of sequence conservation of introns is limited, and the primary sequence differs considerably among the distinct intron sub-groups. It makes intron identification, classification, structural modeling, and the inference of gene models a most challenging and error-prone task—frequently passed on to an “expert” for manual intervention. To reduce the need for manual curation of intron structures and mitochondrial gene models, computational methods using ERPIN sequence profiles were initially developed in 2007. Here we present a refinement of search models and alignments using the now abundant publicly available fungal mtDNA sequences. In addition, we have tested in how far members of the originally proposed sub-groups are clearly distinguished and validated by our computational approach. We confirm clearly distinct mitochondrial Group I sub-groups IA1, IA3, IB3, IC1, IC2, and ID. Yet, IB1, IB2, and IB4 ERPIN models are overlapping substantially in predictions, and are therefore combined and reported as IB. We have further explored the conversion of our ERPIN profiles into covariance models (CM). Current limitations and prospects of the CM approach will be discussed.

Conjugation Operons in Gram-Positive Bacteria with and without Antitermination Systems

Genes involved in the same cellular process are often clustered together in an operon whose expression is controlled by an upstream promoter. Generally, the activity of the promoter is strictly controlled. However, spurious transcription undermines this strict regulation, particularly affecting large operons. The negative effects of spurious transcription can be mitigated by the presence of multiple terminators inside the operon, in combination with an antitermination system. Antitermination systems modify the transcription elongation complexes and enable them to bypass terminators. Bacterial conjugation is the process by which a conjugative DNA element is transferred from a donor to a recipient cell. Conjugation involves many genes that are mostly organized in one or a few large operons. It has recently been shown that many conjugation operons present on plasmids replicating in Gram-positive bacteria possess a bipartite antitermination system that allows not only many terminators inside the conjugation operon to be bypassed, but also the differential expression of a subset of genes. Here, we show that some conjugation operons on plasmids belonging to the Inc18 family of Gram-positive broad host-range plasmids do not possess an antitermination system, suggesting that the absence of an antitermination system may have advantages. The possible (dis)advantages of conjugation operons possessing (or not) an antitermination system are discussed.

Genomic Characterization of a Prophage, Smhb1, That Infects Salinivibrio kushneri BNH Isolated from a Namib Desert Saline Spring

Recent years have seen the classification and reclassification of many viruses related to the model enterobacterial phage P2. Here, we report the identification of a prophage (Smhb1) that infects Salinivibrio kushneri BNH isolated from a Namib Desert salt pan (playa). Analysis of the genome revealed that it showed the greatest similarity to P2-like phages that infect Vibrio species and showed no relation to any of the previously described Salinivibrio-infecting phages. Despite being distantly related to these Vibrio infecting phages and sharing the same modular gene arrangement as seen in most P2-like viruses, the nucleotide identity to its closest relatives suggest that, for now, Smhb1 is the lone member of the Peduovirus genus Playavirus. Although host range testing was not extensive and no secondary host could be identified for Smhb1, genomic evidence suggests that the phage is capable of infecting other Salinivibrio species, including Salinivibrio proteolyticus DV isolated from the same playa. Taken together, the analysis presented here demonstrates how adaptable the P2 phage model can be.

In Silico Prediction for ncRNAs in Prokaryotes

The identification and characterization of non-coding RNAs (ncRNAs) in prokaryotes is an important step in the study of the interaction of these molecules with mRNAs-or target proteins, in the post-transcriptional regulation process. Here, we describe one of the main in silico prediction methods in prokaryotes, using the TargetRNA2 tool to predict target mRNAs.

Comparative Transcriptome and Endophytic Bacterial Community Analysis of Morchella conica SH

The precious rare edible fungus Morchella conica is popular worldwide for its rich nutrition, savory flavor, and varieties of bioactive components. Due to its high commercial, nutritional, and medicinal value, it has always been a hot spot. However, the molecular mechanism and endophytic bacterial communities in M. conica were poorly understood. In this study, we sequenced, assembled, and analyzed the genome of M. conica SH. Transcriptome analysis reveals significant differences between the mycelia and fruiting body. As shown in this study, 1,329 and 2,796 genes were specifically expressed in the mycelia and fruiting body, respectively. The Gene Ontology (GO) enrichment showed that RNA polymerase II transcription activity-related genes were enriched in the mycelium-specific gene cluster, and nucleotide binding-related genes were enriched in the fruiting body-specific gene cluster. Further analysis of differentially expressed genes in different development stages resulted in finding two groups with distinct expression patterns. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment displays that glycan degradation and ABC transporters were enriched in the group 1 with low expressed level in the mycelia, while taurine and hypotaurine metabolismand tyrosine metabolism-related genes were significantly enriched in the group 2 with high expressed level in mycelia. Moreover, a dynamic shift of bacterial communities in the developing fruiting body was detected by 16S rRNA sequencing, and co-expression analysis suggested that bacterial communities might play an important role in regulating gene expression. Taken together, our study provided a better understanding of the molecular biology of M. conica SH and direction for future research on artificial cultivation.

Characterization and Genomic Analysis of BUCT549, a Novel Bacteriophage Infecting Vibrio alginolyticus With Flagella as Receptor

Vibrio alginolyticus is one of the most important of pathogens that can infect humans and a variety of aquatic animals, and it can cause food poisoning and septicemia in humans. Widely used antibiotics are gradually losing their usefulness, and phages are gaining more attention as new antibacterial strategies. To have more potential strategies for controlling pathogenic bacteria, we isolated a novel V. alginolyticus phage BUCT549 from seafood market sewage. It was classified as a new member of the family Siphoviridae by transmission electron microscopy and a phylogenetic tree. We propose creating a new genus for BUCT549 based on the intergenomic similarities (maximum is 56%) obtained from VIRIDIC calculations. Phage BUCT549 could be used for phage therapy due to its stability in a wide pH (3.0–11.0) range and high-temperature (up to 60°C) environment. It had a latent period of 30–40 min and a burst size of 141 PFU/infected bacterium. In the phylogenetic tree based on a terminase large subunit, BUCT549 was closely related to eight Vibrio phages with different species of host. Meanwhile, our experiments proved that BUCT549 has the ability to infect a strain of Vibrio parahaemolyticus. A coevolution experiment determined that three strains of tolerant V. alginolyticus evaded phage infestation by mutating the MSHA-related membrane protein expression genes, which caused the loss of flagellum. This research on novel phage identification and the mechanism of infestation will help phages to become an integral part of the strategy for biological control agents.

A database of flavivirus RNA structures with a search algorithm for pseudoknots and triple base interactions

Motivation

The Flavivirus genus includes several important pathogens, such as Zika, dengue and yellow fever virus. Flavivirus RNA genomes contain a number of functionally important structures in their 3′ untranslated regions (3′UTRs). Due to the diversity of sequences and topologies of these structures, their identification is often difficult. In contrast, predictions of such structures are important for understanding of flavivirus replication cycles and development of antiviral strategies.

Results

We have developed an algorithm for structured pattern search in RNA sequences, including secondary structures, pseudoknots and triple base interactions. Using the data on known conserved flavivirus 3′UTR structures, we constructed structural descriptors which covered the diversity of patterns in these motifs. The descriptors and the search algorithm were used for the construction of a database of flavivirus 3′UTR structures. Validating this approach, we identified a number of domains matching a general pattern of exoribonuclease Xrn1-resistant RNAs in the growing group of insect-specific flaviviruses.

Availability and implementation

The Leiden Flavivirus RNA Structure Database is available at https://rna.liacs.nl. The search algorithm is available at https://github.com/LeidenRNA/SRHS.

Supplementary information

Supplementary data are available at Bioinformatics online.

A novel bipartite antitermination system widespread in conjugative elements of Gram-positive bacteria

Transcriptional regulation allows adaptive and coordinated gene expression, and is essential for life. Processive antitermination systems alter the transcription elongation complex to allow the RNA polymerase to read through multiple terminators in an operon. Here, we describe the discovery of a novel bipartite antitermination system that is widespread among conjugative elements from Gram-positive bacteria, which we named conAn. This system is composed of a large RNA element that exerts antitermination, and a protein that functions as a processivity factor. Besides allowing coordinated expression of very long operons, we show that these systems allow differential expression of genes within an operon, and probably contribute to strict regulation of the conjugation genes by minimizing the effects of spurious transcription. Mechanistic features of the conAn system are likely to decisively influence its host range, with important implications for the spread of antibiotic resistance and virulence genes.

Identification and Evolution of Cas9 tracrRNAs

Clustered regularly interspaced palindromic repeats (CRISPR)-associated (Cas)9 transactivating CRISPR RNAs (tracrRNAs) form distinct structures essential for target recognition and cleavage and dictate exchangeability between orthologous proteins. As noncoding RNAs that are often apart from the CRISPR array, their identification can be arduous. In this article, a new bioinformatic method for the detection of Cas9 tracrRNAs is presented. The approach utilizes a covariance model based on both sequence homology and predicted secondary structure to locate tracrRNAs. This method predicts a tracrRNA for 98% of CRISPR-Cas9 systems identified by us. To ensure accuracy, we also benchmark our approach against biochemically vetted tracrRNAs finding false-positive and false-negative rates of 5.5% and 7.1%, respectively. Finally, the association between Cas9 amino acid sequence-based phylogeny and tracrRNA secondary structure is evaluated, revealing strong evidence that secondary structure is evolutionarily conserved among Cas9 lineages. Altogether, our findings provide insight into Cas9 tracrRNA evolution and efforts to characterize the tracrRNA of Cas9 systems.

Chloroplast acquisition without the gene transfer in kleptoplastic sea slugs, Plakobranchus ocellatus

Some sea slugs sequester chloroplasts from algal food in their intestinal cells and photosynthesize for months. This phenomenon, kleptoplasty, poses a question of how the chloroplast retains its activity without the algal nucleus. There have been debates on the horizontal transfer of algal genes to the animal nucleus. To settle the arguments, this study reported the genome of a kleptoplastic sea slug, Plakobranchus ocellatus, and found no evidence of photosynthetic genes encoded on the nucleus. Nevertheless, it was confirmed that light illumination prolongs the life of mollusk under starvation. These data presented a paradigm that a complex adaptive trait, as typified by photosynthesis, can be transferred between eukaryotic kingdoms by a unique organelle transmission without nuclear gene transfer. Our phylogenomic analysis showed that genes for proteolysis and immunity undergo gene expansion and are up-regulated in chloroplast-enriched tissue, suggesting that these molluskan genes are involved in the phenotype acquisition without horizontal gene transfer.

The Mitogenomes of Ophiostoma minus and Ophiostoma piliferum and Comparisons With Other Members of the Ophiostomatales

Fungi assigned to the Ophiostomatales are of economic concern as many are blue-stain fungi and some are plant pathogens. The mitogenomes of two blue-stain fungi, Ophiostoma minus and Ophiostoma piliferum, were sequenced and compared with currently available mitogenomes for other members of the Ophiostomatales. Species representing various genera within the Ophiostomatales have been examined for gene content, gene order, phylogenetic relationships, and the distribution of mobile elements. Gene synteny is conserved among the Ophiostomatales but some members were missing the atp9 gene. A genome wide intron landscape has been prepared to demonstrate the distribution of the mobile genetic elements (group I and II introns and homing endonucleases) and to provide insight into the evolutionary dynamics of introns among members of this group of fungi. Examples of complex introns or nested introns composed of two or three intron modules have been observed in some species. The size variation among the mitogenomes (from 23.7 kb to about 150 kb) is mostly due to the presence and absence of introns. Members of the genus Sporothrix sensu stricto appear to have the smallest mitogenomes due to loss of introns. The taxonomy of the Ophiostomatales has recently undergone considerable revisions; however, some lineages remain unresolved. The data showed that genera such as Raffaelea appear to be polyphyletic and the separation of Sporothrix sensu stricto from Ophiostoma is justified.

Possible role played by the SINE2 element in gene regulation, as demonstrated by differential processing and polyadenylation in avirulent strains of E. histolytica

Entamoeba histolytica represents a useful model in parasitic organisms due to its complex genomic organization and survival mechanisms. To counteract pathogenic organisms, it is necessary to characterize their molecular biology to design new strategies to combat them. In this report, we investigated a less-known genetic element, short interspersed nuclear element 2 (SINE2), that is present in this ameba and is highly transcribed and polyadenylated. In this study, we show that in two different nonvirulent strains of E. histolytica, SINE2 is differentially processed into two transcript fragments, that is, a full-length 560-nt fragment and a shorter 393-nt fragment bearing an approximately 18-nt polyadenylation tail. Sequence analysis of the SINE2 transcript showed that a Musashi-like protein may bind to it. Also, two putative Musashi-like sequences were identified on the transcript. Semiquantitative expression analysis of the two Musashi-like proteins identified in the E. histolytica genome (XP_648918 and XP_649094) showed that XP_64094 is overexpressed in the nonvirulent strains tested. The information available in the literature and the results presented in this report indicate that SINE2 may affect other genes, as observed with the epigenetic silencing of the G3 strain, by an antisense mechanism or via RNA-protein interactions that may ultimately be involved in the phenotype of nonvirulent strains of E. histolytica.

Characterization of the Bacteriophage vB_EfaS-271 Infecting Enterococcus faecalis

A newly isolated bacteriophage infecting Enterococcus faecalis strains has been characterized, including determination of its molecular features. This phage, named vB_EfaS-271, has been classified as a Siphoviridae member, according to electron microscopy characterization of the virions, composed of a 50 nm-diameter head and a long, flexible, noncontractable tail (219 × 12.5 nm). Analysis of the whole dsDNA genome of this phage showed that it consists of 40,197 bp and functional modules containing genes coding for proteins that are involved in DNA replication (including DNA polymerase/primase), morphogenesis, packaging and cell lysis. Mass spectrometry analysis allowed us to identify several phage-encoded proteins. vB_EfaS-271 reveals a relatively narrow host range, as it is able to infect only a few E. faecalis strains. On the other hand, it is a virulent phage (unable to lysogenize host cells), effectively and quickly destroying cultures of sensitive host bacteria, with a latent period as short as 8 min and burst size of approximately 70 phages per cell at 37 °C. This phage was also able to destroy biofilms formed by E. faecalis. These results contribute to our understanding of the biodiversity of bacteriophages, confirming the high variability among these viruses and indicating specific genetic and functional features of vB_EfaS-271.

Genome-Wide Analysis of RNA Decay in the Cyanobacterium Synechococcus sp. Strain PCC 7002

RNA degradation is an important process that affects the final concentration of individual mRNAs, affecting protein expression and cellular physiology. Studies of how RNA is degraded increase our knowledge of this fundamental process as well as enable the creation of genetic tools to manipulate RNA stability. By studying global transcript turnover, we searched for sequence elements that correlated with transcript (in)stability and used these sequences to guide tool design. This study probes global RNA turnover in a cyanobacterium, Synechococcus sp. strain PCC 7002, that both has a unique array of RNases that facilitate RNA degradation and is an industrially relevant strain that could be used to convert CO₂ and sunlight into useful products.

RNA degradation is an important process that influences the ultimate concentration of individual proteins inside cells. While the main enzymes that facilitate this process have been identified, global maps of RNA turnover are available for only a few species. Even in these cases, there are few sequence elements that are known to enhance or destabilize a native transcript; even fewer confer the same effect when added to a heterologous transcript. To address this knowledge gap, we assayed genome-wide RNA degradation in the cyanobacterium Synechococcus sp. strain PCC 7002 by collecting total RNA samples after stopping nascent transcription with rifampin. We quantified the abundance of each position in the transcriptome as a function of time using RNA-sequencing data and later analyzed the global mRNA decay map using machine learning principles. Half-lives, calculated on a per-ORF (open reading frame) basis, were extremely short, with a median half-life of only 0.97 min. Despite extremely rapid turnover of most mRNA, transcripts encoding proteins involved in photosynthesis were both highly expressed and highly stable. Upon inspection of these stable transcripts, we identified an enriched motif in the 3′ untranslated region (UTR) that had similarity to Rho-independent terminators. We built statistical models for half-life prediction and used them to systematically identify sequence motifs in both 5′ and 3′ UTRs that correlate with stabilized transcripts. We found that transcripts linked to a terminator containing a poly(U) tract had a longer half-life than both those without a poly(U) tract and those without a terminator.

IMPORTANCE RNA degradation is an important process that affects the final concentration of individual mRNAs, affecting protein expression and cellular physiology. Studies of how RNA is degraded increase our knowledge of this fundamental process as well as enable the creation of genetic tools to manipulate RNA stability. By studying global transcript turnover, we searched for sequence elements that correlated with transcript (in)stability and used these sequences to guide tool design. This study probes global RNA turnover in a cyanobacterium, Synechococcus sp. strain PCC 7002, that both has a unique array of RNases that facilitate RNA degradation and is an industrially relevant strain that could be used to convert CO₂ and sunlight into useful products.

Molecular Characterization and Comparative Genomic Analysis of vB_PaeP_YA3, a Novel Temperate Bacteriophage of Pseudomonas aeruginosa

It is well known that bacteriophages play crucial roles in many aspects, such as controlling the number and the diversity of bacteria and participating in horizontal gene transfer, which is a key process in the evolution of bacteria. However, so far, the number of temperate bacteriophages is still limited, and their life processes are severely unknown, except for members of the lambdoid family of coliphages. In this study, a novel temperate phage of Pseudomonas aeruginosa, YA3 (vB_PaeP_YA3), was isolated from waste water. The morphology of YA3 suggested that it is a Podoviridae. The YA3 genome is a circular double-stranded DNA of 45,253 bp, with an average G + C content of 57.2%. A total of 65 open reading frames (ORFs) were predicted according to the sequence of YA3’s genome, of which only 32 (49.2%) ORFs were assigned with putative functions and 13 ORFs were confirmed by the structural proteome. Genome and proteome analyses confirmed the lysogenic nature of this phage, which encodes the typical lysogen-related proteins integrase, CI, Cro, and Q protein. The genome of YA3 is most closely related with that of temperate phage vB_PaeP_Tr60_Ab31, whereas the homology coverage is just 48%. There are many critical differences between their genomes, involving promoters, lysis pathways, and regulation patterns. YA3 is capable of stably lysogenizing its host P. aeruginosa PA14, targeting the integration site within the serine tRNA gene (PA14_RS20820), which is similar with phage vB_PaeP_Tr60_Ab31. The phylogenetic analysis is more complicated than we thought. Based on phage terminase large subunit (TerL) and CI proteins, phage YA3 is related with phage lambda, while their genome coverage is extremely low (<1%). Therefore, phage YA3 is a considerably novel lambda-like temperate phage, and a further study of its genome may deepen our understanding of the interaction between lysogenic phages and their bacterial hosts.

An in silico map of the SARS-CoV-2 RNA Structurome

SARS-CoV-2 is a positive-sense single-stranded RNA virus that has exploded throughout the global human population. This pandemic coronavirus strain has taken scientists and public health researchers by surprise and knowledge of its basic biology (e.g. structure/function relationships in its genomic, messenger and template RNAs) and modes for therapeutic intervention lag behind that of other human pathogens. In this report we used a recently-developed bioinformatics approach, ScanFold, to deduce the RNA structural landscape of the SARS-CoV-2 transcriptome. We recapitulate known elements of RNA structure and provide a model for the folding of an essential frameshift signal. Our results find that the SARS-CoV-2 is greatly enriched in unusually stable and likely evolutionarily ordered RNA structure, which provides a huge reservoir of potential drug targets for RNA-binding small molecules. Our results also predict regions that are accessible for intermolecular interactions, which can aid in the design of antisense therapeutics. All results are made available via a public database (the RNAStructuromeDB) where they may hopefully drive drug discovery efforts to inhibit SARS-CoV-2 pathogenesis.

Genome analysis of Lactobacillus plantarum subsp. plantarum KCCP11226 reveals a well-conserved C30 carotenoid biosynthetic pathway

Carotenoids are group of colored terpenoids with antioxidant properties and widespread in nature including in microorganisms. Lactobacillus plantarum subsp. plantarum KCCP11226 was previously isolated from kimchi, while exhibiting the production of 4,4'-diaponeurosporene as a C₃₀ carotenoid. In this study, full genome sequencing of the strain KCCP11226 was performed. Genome analysis revealed that the dehydrosqualene synthase (crtM) and dehydrosqualene desaturase (crtN) genes, which are major genes for biosynthesis of 4,4'-diaponeurosporene, were shown to act as an operon in most L. plantarum strains, but they were uncommon in other Lactobacillus species. In vitro experiments revealed that the production of 4,4'-diaponeurosporene was greatly increased by oxidative stress. In this situation, mRNA expressions of crtN and crtM were also significantly increased. In conclusion, genome analysis of L. plantarum subsp. plantarum KCCP11226 suggested the presence of a well-conserved C₃₀ carotenoid biosynthetic pathway that includes the crtM-crtN operon. The genomic information on L. plantarum subsp. plantarum KCCP11226 could further elucidate the functions of genes involved in isoprenoid biosynthetic pathway, especially in C₃₀ carotenoid biosynthesis.

Inactivation of the Agrobacterium tumefaciens ActSR system affects resistance to multiple stresses with increased H2O2 sensitivity due to reduced expression of hemH

The Agrobacterium tumefaciens ActSR two-component regulatory system is a member of a homologous group of global redox-responsive regulatory systems that adjust the expression of energy-consuming and energy-supplying metabolic pathways in order to maintain cellular redox balance. In this study, the transcriptional organization of the hrpB-actSR locus was determined and the effect of actSR system inactivation on stress resistance was investigated. It was found that hrpB is transcribed as a monocistronic mRNA and actS is transcribed along with actR as a bicistronic mRNA, while actR is also transcribed as a monocistronic message. Each message is initiated from a separate promoter. Inactivation of actR resulted in decreased resistance to membrane stress (sodium dodecyl sulfate), acid stress (pH 5.5), iron starvation (bipyridyl) and iron excess (FeCl₃), and antibiotic stress (tetracycline and ciprofloxacin). Resistance to oxidative stress in the form of organic peroxide (cumene hydroperoxide) increased, while resistance to inorganic peroxide (H₂O₂) decreased. An actR insertion mutant displayed reduced catalase activity, even though transcription of katA and catE remained unchanged. Complementation of the actR inactivation mutant with plasmid-encoded actR or overexpression of hemH, encoding ferrochelatase, restored wild-type catalase activity and H₂O₂ resistance levels. Gel mobility shift and hemH promoter-lacZ fusion results indicated that ActR is a positive regulator of hemH that binds directly to the hemH promoter region. Thus, inactivation of the A. tumefaciens ActSR system affects resistance to multiple stresses, including reduced resistance to H₂O₂ resulting from a reduction in catalase activity due to reduced expression of hemH.

Molecular characterization of RNase III protein of Asaia sp. for developing a robust RNAi-based paratransgensis tool to affect the sexual life-cycle of Plasmodium or Anopheles fitness

Background

According to scientific recommendations, paratransgenesis is one of the solutions for improving the effectiveness of the Global Malaria Eradication Programme. In paratransgenesis, symbiont microorganisms are used for distorting or blocking the parasite life-cycle, affecting the fitness and longevity of vectors or reducing the vectorial competence. It has been revealed recently that bacteria could be used as potent tools for double stranded RNA production and delivery to insects. Moreover, findings showed that RNase III mutant bacteria are more competent for this aim. Asaia spp. have been introduced as potent paratransgenesis candidates for combating malaria and, based on their specific features for this goal, could be considered as effective dsRNA production and delivery tools to Anopheles spp. Therefore, we decided to characterize the rnc gene and its related protein to provide the basic required information for creating an RNase III mutant Asaia bacterium.

Methods

Asaia bacteria were isolated from field-collected Anopheles stephensi mosquitoes. The rnc gene and its surrounding sequences were characterized by rapid amplification of genomic ends. RNase III recombinant protein was expressed in E. coli BL21 and biological activity of the purified recombinant protein was assayed. Furthermore, Asaia RNaseIII amino acid sequence was analyzed by in silico approaches such as homology modeling and docking to determine its structural properties.

Results

In this study, the structure of rnc gene and its related operon from Asaia sp. was determined. In addition, by performing superimposition and docking with specific substrate, the structural features of Asaia RNaseIII protein such as critical residues which are involved and essential for proper folding of active site, binding of magnesium ions and double stranded RNA molecule to protein and cleaving of dsRNA molecules, were determined.

Conclusions

In this study, the basic and essential data for creating an RNase III mutant Asaia sp. strain, which is the first step of developing an efficient RNAi-based paratransgenesis tool, were acquired. Asaia sp. have been found in different medically-important vectors and these data are potentially very helpful for researchers studying paratransgenesis and vector-borne diseases and are interested in applying the RNAi technology in the field.

Computational Tools for Discovering and Engineering Natural Product Biosynthetic Pathways

Natural products (NPs), also known as secondary metabolites, are produced in bacteria, fungi, and plants. NPs represent a rich source of antibacterial, antifungal, and anticancer agents. Recent advances in DNA sequencing technologies and bioinformatics unveiled nature's great potential for synthesizing numerous NPs that may confer unprecedented structural and biological features. However, discovering novel bioactive NPs by genome mining remains a challenge. Moreover, even with interesting bioactivity, the low productivity of many NPs significantly limits their practical applications. Here we discuss the progress in developing bioinformatics tools for efficient discovery of bioactive NPs. In addition, we highlight computational methods for optimizing the productivity of NPs of pharmaceutical importance.

The complete mitochondrial genome of Gruberia lanceolata (Gruber, 1884) Kahl, 1932 (Ciliophora: Heterotrichea)

The ciliate Gruberia lanceolata (Gruber, 1884) Kahl, 1932 belonging to the class Heterotrichea was sampled from the coastal waters of South Korea. The complete mitogenome in its linear form and large size (∼40 kb) was obtained. It consisted of 27 protein-coding genes (PCGs), two ribosomal subunit RNA (rRNA) genes, four transfer RNAs (tRNAs), and ten unclassified open reading frames (ORFs). Their telomeric structures were capped, with repeat regions at both ends. We analyzed its phylogenetic tree using the data of its respiratory chain complex I genes. It can be suggested that the complete mitochondrial genome of G. lanceolata can be recorded as a new class of the mitogenome.

Within-Host Adaptation Mediated by Intergenic Evolution in Pseudomonas aeruginosa

Bacterial pathogens evolve during the course of infection as they adapt to the selective pressures that confront them inside the host. Identification of adaptive mutations and their contributions to pathogen fitness remains a central challenge. Although mutations can either target intergenic or coding regions in the pathogen genome, studies of host adaptation have focused predominantly on molecular evolution within coding regions, whereas the role of intergenic mutations remains unclear. Here, we address this issue and investigate the extent to which intergenic mutations contribute to the evolutionary response of a clinically important bacterial pathogen, Pseudomonas aeruginosa, to the host environment, and whether intergenic mutations have distinct roles in host adaptation. We characterize intergenic evolution in 44 clonal lineages of P. aeruginosa and identify 77 intergenic regions in which parallel evolution occurs. At the genetic level, we find that mutations in regions under selection are located primarily within regulatory elements upstream of transcriptional start sites. At the functional level, we show that some of these mutations both increase or decrease transcription of genes and are directly responsible for evolution of important pathogenic phenotypes including antibiotic sensitivity. Importantly, we find that intergenic mutations facilitate essential genes to become targets of evolution. In summary, our results highlight the evolutionary significance of intergenic mutations in creating host-adapted strains, and that intergenic and coding regions have different qualitative contributions to this process.