Genome-wide analyses in bacteria show small-RNA enrichment for long and conserved intergenic regions

Quantitative morphological analysis of Deinococcus radiodurans elucidates complex dose-dependent nucleoid condensation during recovery from ionizing radiation

The extremophile Deinococcus radiodurans maintains a highly organized and condensed nucleoid as its default state, possibly contributing to its high tolerance to ionizing radiation (IR). Previous studies of the D. radiodurans nucleoid were limited by reliance on manual image annotation and qualitative metrics. Here, we introduce a high-throughput approach to quantify the geometric properties of cells and nucleoids using confocal microscopy, digital reconstructions of cells, and computational modeling. We utilize this novel approach to investigate the dynamic process of nucleoid condensation in response to IR stress. Our quantitative analysis reveals that at the population level, exposure to IR induced nucleoid compaction and decreased the size of D. radiodurans cells. Morphological analysis and clustering identified six distinct sub-populations across all tested experimental conditions. Results indicate that exposure to IR induced fractional redistributions of cells across sub-populations to exhibit morphologies associated with greater nucleoid condensation and decreased the abundance of sub-populations associated with cell division. Nucleoid-associated proteins (NAPs) may link nucleoid compaction and stress tolerance, but their roles in regulating compaction in D. radiodurans are unknown. Imaging of genomic mutants of known and suspected NAPs that contribute to nucleoid condensation found that deletion of nucleic acid-binding proteins, not previously described as NAPs, can remodel the nucleoid by driving condensation or decondensation in the absence of stress and that IR increased the abundance of these morphological states. Thus, our integrated analysis introduces a new methodology for studying environmental influences on bacterial nucleoids and provides an opportunity to further investigate potential regulators of nucleoid condensation.IMPORTANCEDeinococcus radiodurans, an extremophile known for its stress tolerance, constitutively maintains a highly condensed nucleoid. Qualitative studies have described nucleoid behavior under a variety of conditions. However, a lack of quantitative data regarding nucleoid organization and dynamics has limited our understanding of the regulatory mechanisms controlling nucleoid organization in D. radiodurans. Here, we introduce a quantitative approach that enables high-throughput quantitative measurements of subcellular spatial characteristics in bacterial cells. Applying this to wild-type or single-protein-deficient populations of D. radiodurans subjected to ionizing radiation, we identified significant stress-responsive changes in cell shape, nucleoid organization, and morphology. These findings highlight this methodology's adaptability and capacity for quantitatively analyzing the cellular response to stressors for screening cellular proteins involved in bacterial nucleoid organization.

StORF-Reporter: finding genes between genes

Large regions of prokaryotic genomes are currently without any annotation, in part due to well-established limitations of annotation tools. For example, it is routine for genes using alternative start codons to be misreported or completely omitted. Therefore, we present StORF-Reporter, a tool that takes an annotated genome and returns regions that may contain missing CDS genes from unannotated regions. StORF-Reporter consists of two parts. The first begins with the extraction of unannotated regions from an annotated genome. Next, Stop-ORFs (StORFs) are identified in these unannotated regions. StORFs are open reading frames that are delimited by stop codons and thus can capture those genes most often missing in genome annotations. We show this methodology recovers genes missing from canonical genome annotations. We inspect the results of the genomes of model organisms, the pangenome of Escherichia coli, and a set of 5109 prokaryotic genomes of 247 genera from the Ensembl Bacteria database. StORF-Reporter extended the core, soft-core and accessory gene collections, identified novel gene families and extended families into additional genera. The high levels of sequence conservation observed between genera suggest that many of these StORFs are likely to be functional genes that should now be considered for inclusion in canonical annotations.

In silico analysis of Leishmania proteomes and protein-protein interaction network: Prioritizing therapeutic targets and drugs for repurposing to treat leishmaniasis

Leishmaniasis is a serious world health problem and its current therapies have several limitations demanding to develop novel therapeutics for this disease. The present study aims to prioritize novel broad-spectrum targets using proteomics and protein-protein interaction network (PPIN) data for 11 Leishmania species. Proteome comparison and host non-homology analysis resulted in 3605 pathogen-specific conserved core proteins. Gene ontology analysis indicated their involvement in major molecular functions like DNA binding, transportation, dioxygenase, and catalytic activity. PPIN analysis of these core proteins identified eight hub proteins (viz., vesicle-trafficking protein (LBRM2903_190011800), ribosomal proteins S17 (LBRM2903_34004790) and L2 (LBRM2903_080008100), eukaryotic translation initiation factor 3 (LBRM2903_350086700), replication factor A (LBRM2903_150008000), U3 small nucleolar RNA-associated protein (LBRM2903_340025600), exonuclease (LBRM2903_200021800), and mitochondrial RNA ligase (LBRM2903_200074100)). Among the hub proteins, six were classified as drug targets and two as vaccine candidates. Further, druggability analysis indicated three hub proteins, namely eukaryotic translation initiation factor 3, ribosomal proteins S17 and L2 as druggable. Their three-dimensional structures were modelled and docked with the identified ligands (2-methylthio-N6-isopentenyl-adenosine-5'-monophosphate, artenimol and omacetaxine mepesuccinate). These ligands could be experimentally validated (in vitro and in vivo) and repurposed for the development of novel antileishmanial agents.

Genome-wide high-resolution mapping of DNA methylation reveals epigenetic variation in the offspring of sexual and asexual propagation in Robinia pseudoacacia

We detected the genome-wide pattern of DNA methylation and its association with gene expression in sexual and asexual progenies of mature Robinia pseudoacacia trees. DNA methylation plays an important role in plant reproduction and development. Although some studies on sexual reproduction have been carried out in model plants, little is known about the dynamic changes in DNA methylation and their effect on gene expression in sexual and asexual progeny of woody plants. Here, through whole-genome bisulfite sequencing, we revealed DNA methylation patterns in the sexual and asexual progenies of mature Robinia pseudoacacia to understand the regulation of gene expression by DNA methylation in juvenile seedlings. An average of 53% CG, 34% CHG and 5% CHH contexts was methylated in the leaves of mature and juvenile individuals. The CHH methylation level of asexually propagated seedlings was significantly lower than that of seed-derived seedlings and mature trees. The intergenic regions had the highest methylation level. Analysis of differentially methylated regions (DMRs) showed that most of them were hypermethylated and located in the gene upstream and introns. A total of 24, 108 and 162 differentially expressed genes containing DMRs were identified in root sprouts (RSs), root cuttings (RCs) and seed-derived seedlings (SSs), respectively, and a large proportion of them showed hypermethylation. In addition, DMRs were enriched within GO subcategories including catalytic activity, metabolic process and cellular process. The results reveal widespread DNA methylation changes between mature plants and their progenies through sexual/asexual reproduction, which provides novel insights into DNA methylation reprogramming and the regulation of gene expression in woody plants.

Pan-Genome of the Genus Streptomyces and Prioritization of Biosynthetic Gene Clusters With Potential to Produce Antibiotic Compounds

Species of the genus Streptomyces are known for their ability to produce multiple secondary metabolites; their genomes have been extensively explored to discover new bioactive compounds. The richness of genomic data currently available allows filtering for high quality genomes, which in turn permits reliable comparative genomics studies and an improved prediction of biosynthetic gene clusters (BGCs) through genome mining approaches. In this work, we used 121 genome sequences of the genus Streptomyces in a comparative genomics study with the aim of estimating the genomic diversity by protein domains content, sequence similarity of proteins and conservation of Intergenic Regions (IGRs). We also searched for BGCs but prioritizing those with potential antibiotic activity. Our analysis revealed that the pan-genome of the genus Streptomyces is clearly open, with a high quantity of unique gene families across the different species and that the IGRs are rarely conserved. We also described the phylogenetic relationships of the analyzed genomes using multiple markers, obtaining a trustworthy tree whose relationships were further validated by Average Nucleotide Identity (ANI) calculations. Finally, 33 biosynthetic gene clusters were detected to have potential antibiotic activity and a predicted mode of action, which might serve up as a guide to formulation of related experimental studies.

A Translation-Aborting Small Open Reading Frame in the Intergenic Region Promotes Translation of a Mg2+ Transporter in Salmonella Typhimurium

Translation initiation regions in mRNAs that include the ribosome-binding site (RBS) and the start codon are often sequestered within a secondary structure. Therefore, to initiate protein synthesis, the mRNA secondary structure must be unfolded to allow the RBS to be accessible to the ribosome.

Bacterial mRNAs often harbor upstream open reading frames (uORFs) in the 5′ untranslated regions (UTRs). Translation of the uORF usually affects downstream gene expression at the levels of transcription and/or translation initiation. Unlike other uORFs mostly located in the 5′ UTR, we discovered an 8-amino-acid ORF, designated mgtQ, in the intergenic region between the mgtC virulence gene and the mgtB Mg²⁺ transporter gene in the Salmonella mgtCBRU operon. Translation of mgtQ promotes downstream mgtB Mg²⁺ transporter expression at the level of translation by releasing the ribosome-binding sequence of the mgtB gene that is sequestered in a translation-inhibitory stem-loop structure. Interestingly, mgtQ Asp2 and Glu5 codons that induce ribosome destabilization are required for mgtQ-mediated mgtB translation. Moreover, the mgtQ Asp and Glu codons-mediated mgtB translation is counteracted by the ribosomal subunit L31 that stabilizes ribosome. Substitution of the Asp2 and Glu5 codons in mgtQ decreases MgtB Mg²⁺ transporter production and thus attenuates Salmonella virulence in mice, likely by limiting Mg²⁺ acquisition during infection.

Divergently Transcribed ncRNAs in Escherichia coli: Refinement of the Transcription Starts Assumes Functional Diversification

Non-coding regulatory RNAs (ncRNAs) comprise specialized group of essential genetically encoded biological molecules involved in the wide variety of cellular metabolic processes. The progressive increase in the number of newly identified ncRNAs and the defining of their genome location indicate their predominant nesting in intergenic regions and expression under the control of their own regulatory elements. At the same time, the regulation of ncRNA's transcription cannot be considered in isolation from the processes occurring in the immediate genetic environment. A number of experimental data indicate the notable impact of positional regulation of gene expression mediated by dynamic temporal DNA rearrangements accompanying transcription events in the vicinity of neighboring genes. This issue can be perceived as particularly significant for divergently transcribed ncRNAs being actually subjected to double regulatory pressure. Based on available results of RNAseq experiments for Escherichia coli, we screened out divergent ncRNAs and the adjacent genes for the exact positions of transcription start sites (TSSs) and relative efficiency of RNA production. This analysis revealed extension or shortening of some previously annotated ncRNAs resulting in modified secondary structure, confirmed stable expression of four ncRNAs annotated earlier as putative, and approved the possibility of expression of divergently transcribed ncRNAs containing repetitive extragenic palindromic (REP) elements. The biogenesis of secreted ncRNAs from divergently transcribed ffs, chiX, ralA, and ryhB is discussed taking into account positions of TSSs. Refinement of TSSs for the neighboring genes renders some ncRNAs as true antisense overlapping with 5'UTR of divergently transcribed mRNAs.

Genome-wide identification of Azospirillum brasilense Sp245 small RNAs responsive to nitrogen starvation and likely involvement in plant-microbe interactions

BACKGROUND

Small RNAs (sRNAs) are non-coding RNAs known to regulate various biological functions such as stress adaptation, metabolism, virulence as well as pathogenicity across a wide range of bacteria, mainly by controlling mRNA stabilization or regulating translation. Identification and functional characterization of sRNAs has been carried out in various plant growth-promoting bacteria and they have been shown to help the cells cope up with environmental stress. No study has been carried out to uncover these regulatory molecules in the diazotrophic alpha-proteobacterium Azospirillum brasilense Sp245 to date.

RESULTS

Expression-based sRNA identification (RNA-seq) revealed the first list of ~ 468 sRNA candidate genes in A. brasilense Sp245 that were differentially expressed in nitrogen starvation versus non-starved conditions. In parallel, in silico tools also identified 2 of the above as candidate sRNAs. Altogether, putative candidates were stringently curated from RNA-seq data based on known sRNA parameters (size, location, secondary structure, and abundance). In total, ~ 59 significantly expressed sRNAs were identified in this study of which 53 are potentially novel sRNAs as they have no Rfam and BSRD homologs. Sixteen sRNAs were randomly selected and validated for differential expression, which largely was found to be in congruence with the RNA-seq data.

CONCLUSIONS

Differential expression of 468 A. brasilense sRNAs was indicated by RNA-seq data, a subset of which was confirmed by expression analysis. Four of the significantly expressed sRNAs were not observed in nitrogen starvation while 16 sRNAs were found to be exclusively expressed in nitrogen depletion. Putative candidate sRNAs identified have potential mRNA targets primarily involved in stress (abiotic and biotic) adaptability; regulation of bacterial cellular, biological and molecular pathways such as nitrogen fixation, polyhydroxybutyrate synthesis, chemotaxis, biofilm formation and transcriptional regulation. In addition to directly influencing bacteria, some of these sRNAs also have targets influencing plant-microbe interactions through adhesion of bacteria to plant roots directly, suppressing host response, inducing plant defence and signalling.

Metatranscriptomic analyses of the oral microbiome

Although the composition of the oral human microbiome is now well studied, regulation of genes within oral microbial communities remains mostly uncharacterized. Current concepts of periodontal disease and caries highlight the importance of oral biofilms and their role as etiological agents of those diseases. Currently, there is increased interest in exploring and characterizing changes in the composition and gene-expression profiles of oral microbial communities. These efforts aim to identify changes in functional activities that could explain the transition from health to disease and the reason for the chronicity of those infections. It is now clear that the functions of distinct species within the subgingival microbiota are intimately intertwined with the rest of the microbial community. This point highlights the relevance of examining the expression profile of specific species within the subgingival microbiota in the case of periodontal disease or caries lesions, in the context of the other members of the biofilm in vivo. Metatranscriptomic analysis of the oral community is the starting point for identifying environmental signals that modulate the shift in metabolism of the community from commensal to dysbiotic. These studies give a snapshot of the expression patterns of microbial communities and also allow us to determine triggers to diseases. For example, in the case of caries, studies have unveiled a potential new pathway of sugar metabolism, namely the use of sorbitol as an additional source of carbon by Streptococcus mutans; and in the case of periodontal disease, high levels of extracellular potassium could be a signal of disease. Longitudinal studies are needed to identify the real markers of the initial stages of caries and periodontal disease. More information on the gene-expression profiles of the host, along with the patterns from the microbiome, will lead to a clearer understanding of the modulation of health and disease. This review presents a summary of these initial studies, which have opened the door to a new understanding of the dynamics of the oral community during the dysbiotic process in the oral cavity.

The Transcriptomic Landscape of Cupriavidus metallidurans CH34 Acutely Exposed to Copper

Bacteria are increasingly used for biotechnological applications such as bioremediation, biorecovery, bioproduction, and biosensing. The development of strains suited for such applications requires a thorough understanding of their behavior, with a key role for their transcriptomic landscape. We present a thorough analysis of the transcriptome of Cupriavidus metallidurans CH34 cells acutely exposed to copper by tagRNA-sequencing. C. metallidurans CH34 is a model organism for metal resistance, and its potential as a biosensor and candidate for metal bioremediation has been demonstrated in multiple studies. Several metabolic pathways were impacted by Cu exposure, and a broad spectrum of metal resistance mechanisms, not limited to copper-specific clusters, was overexpressed. In addition, several gene clusters involved in the oxidative stress response and the cysteine-sulfur metabolism were induced. In total, 7500 transcription start sites (TSSs) were annotated and classified with respect to their location relative to coding sequences (CDSs). Predicted TSSs were used to re-annotate 182 CDSs. The TSSs of 2422 CDSs were detected, and consensus promotor logos were derived. Interestingly, many leaderless messenger RNAs (mRNAs) were found. In addition, many mRNAs were transcribed from multiple alternative TSSs. We observed pervasive intragenic TSSs both in sense and antisense to CDSs. Antisense transcripts were enriched near the 5' end of mRNAs, indicating a functional role in post-transcriptional regulation. In total, 578 TSSs were detected in intergenic regions, of which 35 were identified as putative small regulatory RNAs. Finally, we provide a detailed analysis of the main copper resistance clusters in CH34, which include many intragenic and antisense transcripts. These results clearly highlight the ubiquity of noncoding transcripts in the CH34 transcriptome, many of which are putatively involved in the regulation of metal resistance.

Investigation into archaeal extremophilic lifestyles through comparative proteogenomic analysis

Archaea are a group of primary life forms on Earth and could thrive in many unique environments. Their successful colonization of extreme niches requires corresponding adaptations at proteogenomic level in order to maintain stable cellular structures and active physiological functions. Although some studies have already investigated the extremophilic lifestyles of archaeal species based on genomic features and protein structures, there is a lack of comparative proteogenomic analysis in a large scale. In this study, we explored 686 high-quality archaeal genomes (proteomes) sourced from the Pathosystems Resource Integration Center (PATRIC) database. General patterns of genomic features such as genome size, coding capacity (coding genes and non-coding regions), and G + C contents were re-confirmed. Protein domain distribution patterns were then identified across archaeal species. Domains with unknown functions (DUFs) and mini proteins were investigated in terms of their distributions due to their importance in archaeal physiological functions. In addition, physicochemical properties of protein sequences, such as stability, hydrophobicity, isoelectric point, aromaticity and amino acid compositions in corresponding archaeal groups were compared. Unique features associated with extremophilic lifestyles were observed, which suggested that evolutionary adaptations to different extreme environments had intrinsic impacts on archaeal protein features. Taken together, this systematic study facilitates a better understanding of the mechanisms behind the extremophilic lifestyles of archaeal species, which will further contribute to the evolutionary explorations of archaeal adaptations both experimentally and theoretically in the future studies.Communicated by Ramaswamy H. Sarma.

Potential regulation of small RNAs on bacterial function activities in pig farm wastewater treatment plants

Small RNAs (sRNAs) are key players in the regulation of bacterial gene expression. However, the distribution and regulatory functions of sRNA in pig farm wastewater treatment plants (WWTPs) remains unknown. In this study, the wastewaters in anoxic and oxic tanks of the WWTPs were collected. The profiles of the community structure, mRNA expression, and sRNA expression of bacteria in pig farm wastewater were investigated using transcriptome sequencing and qPCR. This study demonstrated that there was a higher abundance of sRNA in the pig farm WWTPs and 52 sRNAs were detected. The sRNAs were mainly present in Proteobacteria and Firmicutes, including the potential human pathogenic bacteria (HPB) (Escherichia, Shigella, Bordetella and Morganella), crop pathogen (Pectobacterium) and denitrifying bacteria (Zobellella). And the sRNAs were involved in the bacterial functional activities such as translation, transcription, drug resistance, membrane transport and amino acid metabolism. In addition, most sRNAs had a higher abundance in anoxic tanks which contained a higher abundance of the genes associated with infectious diseases and drug resistance than that in oxic tanks. The results presented here show that in pig farm WWTPs, sRNA played an important role in bacterial function activities, especially the infectious diseases, drug resistance and denitrification, which can provide a new point of penetration for improving the pig farm WWTPs.

Prevalence of small base-pairing RNAs derived from diverse genomic loci

Small RNAs (sRNAs) that act by base-pairing have been shown to play important roles in fine-tuning the levels and translation of their target transcripts across a variety of model and pathogenic organisms. Work from many different groups in a wide range of bacterial species has provided evidence for the importance and complexity of sRNA regulatory networks, which allow bacteria to quickly respond to changes in their environment. However, despite the expansive literature, much remains to be learned about all aspects of sRNA-mediated regulation, particularly in bacteria beyond the well-characterized Escherichia coli and Salmonella enterica species. Here we discuss what is known, and what remains to be learned, about the identification of regulatory base-pairing RNAs produced from diverse genomic loci including how their expression is regulated. This article is part of a Special Issue entitled: RNA and gene control in bacteria edited by Dr. M. Guillier and F. Repoila.

Multiple Small RNAs Interact to Co-regulate Ethanol Tolerance in Zymomonas mobilis

sRNAs represent a powerful class of regulators that influences multiple mRNA targets in response to environmental changes. However, very few direct sRNA-sRNA interactions have been deeply studied in any organism. Zymomonas mobilis is a bacterium with unique ethanol-producing metabolic pathways in which multiple small RNAs (sRNAs) have recently been identified, some of which show differential expression in ethanol stress. In this study, we show that two sRNAs (Zms4 and Zms6) are upregulated under ethanol stress and have significant impacts on ethanol tolerance and production in Z. mobilis. We conducted multi-omics analysis (combining transcriptomics and sRNA-immunoprecipitation) to map gene networks under the influence of their regulation. We confirmed that Zms4 and Zms6 bind multiple RNA targets and regulate their expressions, influencing many downstream pathways important to ethanol tolerance and production. In particular, Zms4 and Zms6 interact with each other as well as many other sRNAs, forming a novel sRNA-sRNA direct interaction network. This study thus uncovers a sRNA network that co-orchestrates multiple ethanol related pathways through a diverse set of mRNA targets and a large number of sRNAs. To our knowledge, this study represents one of the largest sRNA-sRNA direct interactions uncovered so far.

Identification of small RNAs involved in nitrogen fixation in Anabaena sp. PCC 7120 based on RNA-seq under steady state conditions

Purpose

Anabaena sp. PCC7120 is a genetically tractable model organism for nitrogen fixation and photosynthesis research. The importance of small regulatory RNAs (sRNAs) as mediators of a number of cellular processes in bacteria has begun to be recognized. Bacterial sRNA binds to target genes through base pairing, and play a regulatory role. Many studies have shown that bacterial sRNA can regulate cell stress response, carbon and nitrogen fixation, and so on. However, little is known about sRNAs in Anabaena sp. PCC 7120 regarded to nitrogen fixation under later steady state.

Methods

To provide a comprehensive study of sRNAs in this model organism, the sRNA (< 200 nt) extracted from Anabaena sp. PCC 7120 under nitrogen step-down treatment of 12 days, together with the sRNA from the control, was analyzed using deep RNA sequencing. Possible target genes regulated by all identified putative sRNAs were predicted by IntaRNA and further analyzed for functional categorizations for biological pathways.

Result

Totally, 14,132 transcripts were produced from the de novo assembly. Among them, transcripts that are located either in the intergenic region or antisense strand were kept, which resulted in 1219 sRNA candidates, for further analysis. RPKM-based differential expression analysis showed that 418 sRNAs were significantly differentially expressed between the samples from control (nitrogen addition, N+) and nitrogen depletion, (N−). Among them, 303 sRNAs were significantly upregulated, whereas 115 sRNAs were significantly downregulated. RT-PCR of 18 randomly chosen sRNAs showed a similar pattern as RNA-seq result, which confirmed the reliability of the RNA-seq data. In addition, the possible target genes regulated by unique sRNAs of Anabaena sp. PCC 7120 under nitrogen addition (N+) condition or that under nitrogen depletion (N−) condition were analyzed for functional categorization and biological pathways, which provided the evidences that sRNAs were indeed involved in many different metabolic pathways.

Conclusion

The information from the present study provides a valuable reference for understanding the sRNA-mediated regulation of the nitrogen fixation in Anabaena PCC 7120 under steady state conditions.

Genome-wide identification of the context-dependent sRNA expression in Mycobacterium tuberculosis

BACKGROUND

Tuberculosis remains one of the leading causes of morbidity and mortality worldwide. Therefore, understanding the pathophysiology of Mycobacterium tuberculosis is imperative for developing new drugs. Post-transcriptional regulation plays a significant role in microbial adaptation to different growth conditions. While the proteins associated with gene expression regulation have been extensively studied in the pathogenic strain M. tuberculosis H37Rv, post-transcriptional regulation involving small RNAs (sRNAs) remains poorly understood.

RESULTS

We developed a novel moving-window based approach to detect sRNA expression using RNA-Seq data. Overlaying ChIP-seq data of RNAP (RNA Polymerase) and NusA suggest that these putative sRNA coding regions are significantly bound by the transcription machinery. Besides capturing many experimentally validated sRNAs, we observe the context-dependent expression of novel sRNAs in the intergenic regions of M. tuberculosis genome. For example, ncRv11806 shows expression only in the stationary phase, suggesting its role in mycobacterial latency which is a key attribute to long term pathogenicity. Also, ncRv11875C showed expression in the iron-limited condition, which is prevalent inside the macrophages of the host cells.

CONCLUSION

The systems level analysis of sRNA highlights the condition-specific expression of sRNAs which might enable the pathogen survival by rewiring regulatory circuits.

Applying a New REFINE Approach in Zymomonas mobilis Identifies Novel sRNAs That Confer Improved Stress Tolerance Phenotypes

As global controllers of gene expression, small RNAs represent powerful tools for engineering complex phenotypes. However, a general challenge prevents the more widespread use of sRNA engineering strategies: mechanistic analysis of these regulators in bacteria lags far behind their high-throughput search and discovery. This makes it difficult to understand how to efficiently identify useful sRNAs to engineer a phenotype of interest. To help address this, we developed a forward systems approach to identify naturally occurring sRNAs relevant to a desired phenotype: RNA-seq Examiner for Phenotype-Informed Network Engineering (REFINE). This pipeline uses existing RNA-seq datasets under different growth conditions. It filters the total transcriptome to locate and rank regulatory-RNA-containing regions that can influence a metabolic phenotype of interest, without the need for previous mechanistic characterization. Application of this approach led to the uncovering of six novel sRNAs related to ethanol tolerance in non-model ethanol-producing bacterium Zymomonas mobilis. Furthermore, upon overexpressing multiple sRNA candidates predicted by REFINE, we demonstrate improved ethanol tolerance reflected by up to an approximately twofold increase in relative growth rate compared to controls not expressing these sRNAs in 7% ethanol (v/v) RMG-supplemented media. In this way, the REFINE approach informs strain-engineering strategies that we expect are applicable for general strain engineering.

Comparative Transcriptomic and Functional Assessments of Linezolid-Responsive Small RNA Genes in Staphylococcus aureus

Bacterial small RNAs (sRNAs) are RNA molecules that can have important regulatory roles across gene expression networks. There is a growing understanding of the scope and potential breadth of impact of sRNAs on global gene expression patterns in Staphylococcus aureus, a major human pathogen. Here, transcriptome comparisons were used to examine the roles of sRNA genes with a potential role in the response of S. aureus to antibiotic exposure. Although no measurable impact on key bacterial phenotypes was observed after deleting each of 18 sRNAs identified by these comparisons, this research is significant because it underscores the subtle modes of action of these sometimes abundant molecules within the bacterium.

Staphylococcus aureus contains a repertoire of at least 50 and possibly 500 small RNAs (sRNAs). The functions of most sRNAs are not understood, although some are known to respond to environmental changes, including the presence of antibiotics. Here, in an effort to better understand the roles of sRNAs in the context of antibiotic exposure, we took a clinical methicillin-resistant S. aureus (MRSA) isolate and separately deleted eight sRNAs that were significantly upregulated in response to the last-line antibiotic linezolid as revealed by transcriptome sequencing (RNA-seq) comparisons. We also deleted an additional 10 sRNAs that were either highly expressed or previously found to respond to antibiotic exposure. There were no significant changes for any of the 18 mutants in a variety of phenotypic screens, including MIC screens, growth competition assays in the presence of linezolid, biofilm formation, and resistance to whole-blood killing. These data suggest sRNA functional redundancy, because despite their high expression levels upon antibiotic exposure, individual sRNA genes do not affect readily observable bacterial phenotypes. The sRNA transcriptional changes we measured during antibiotic exposure might also reflect sRNA “indifference,” that is, a general stress response not specifically related to sRNA function. These data underscore the need for sensitive assays and new approaches to try and decipher the functions of sRNA genes in S. aureus.

IMPORTANCE Bacterial small RNAs (sRNAs) are RNA molecules that can have important regulatory roles across gene expression networks. There is a growing understanding of the scope and potential breadth of impact of sRNAs on global gene expression patterns in Staphylococcus aureus, a major human pathogen. Here, transcriptome comparisons were used to examine the roles of sRNA genes with a potential role in the response of S. aureus to antibiotic exposure. Although no measurable impact on key bacterial phenotypes was observed after deleting each of 18 sRNAs identified by these comparisons, this research is significant because it underscores the subtle modes of action of these sometimes abundant molecules within the bacterium.

It's a small, small world: unravelling the role and evolution of small RNAs in organelle and endosymbiont genomes

Organelles and host-restricted bacterial symbionts are characterized by having highly reduced genomes that lack many key regulatory genes and elements. Thus, it has been hypothesized that the eukaryotic nuclear genome is primarily responsible for regulating these symbioses. However, with the discovery of organelle- and symbiont-expressed small RNAs (sRNAs) there is emerging evidence that these sRNAs may play a role in gene regulation as well. Here, we compare the diversity of organelle and bacterial symbiont sRNAs recently identified using genome-enabled '-omic' technologies and discuss their potential role in gene regulation. We also discuss how the genome architecture of small genomes may influence the evolution of these sRNAs and their potential function. Additionally, these new studies suggest that some sRNAs are conserved within organelle and symbiont taxa and respond to changes in the environment and/or their hosts. In summary, these results suggest that organelle and symbiont sRNAs may play a role in gene regulation in addition to nuclear-encoded host mechanisms.

RNA sequencing reveals small RNAs in Bacillus pumilus under different growth phases of the protease fermentation process

Bacillus pumilus, an endospore-forming soil bacterium, produces a wide array of extracellular proteins, such as proteases, which are already applied in the chemical, detergent and leather industries. Small noncoding regulatory RNAs (sRNAs) in bacteria are important RNA regulators that act in response to various environmental signals. Here, an RNA-seq-based transcriptome analysis was applied to B. pumilus SCU11, a strain that produces extracellular alkaline protease, across various growth phases of the protease fermentation process. Through bioinformatics screening of the sequencing data and visual inspection, 84 putative regulatory sRNAs were identified in B. pumilus, including 21 antisense sRNAs and 63 sRNAs in intergenic regions. We experimentally validated the expression of 48 intergenic sRNAs by quantitative RT-PCR (qRT-PCR). Meanwhile, the expression of 6 novel sRNAs was confirmed by northern blotting, and the expression profiles of 5 sRNAs showed close correlation with the growth phase. We revealed that the sRNA Bpsr137 was involved in flagellum and biofilm formation in B. pumilus. The identification of a global set of sRNAs increases the inventory of regulatory sRNAs in Bacillus and implies the important regulatory roles of sRNA in B. pumilus. These findings will contribute another dimension to the optimization of crucial metabolic activities of B. pumilus during a productive fermentation process.

Synthetic Biology of Small RNAs and Riboswitches

In bacteria and archaea, small RNAs (sRNAs) regulate complex networks through antisense interactions with target mRNAs in trans, and riboswitches regulate gene expression in cis based on the ability to bind small-molecule ligands. Although our understanding and characterization of these two important regulatory RNA classes is far from complete, these RNA-based mechanisms have proven useful for a wide variety of synthetic biology applications. Besides classic and contemporary applications in the realm of metabolic engineering and orthogonal gene control, this review also covers newer applications of regulatory RNAs as biosensors, logic gates, and tools to determine RNA-RNA interactions. A separate section focuses on critical insights gained and challenges posed by fundamental studies of sRNAs and riboswitches that should aid future development of synthetic regulatory RNAs.

Functional annotation and distribution overview of RNA families in 27 Streptococcus agalactiae genomes

Background

Streptococcus agalactiae, also known as Group B Streptococcus (GBS), is a Gram-positive bacterium that colonizes the gastrointestinal and genitourinary tract of humans. This bacterium has also been isolated from various animals, such as fish and cattle. Non-coding RNAs (ncRNAs) can act as regulators of gene expression in bacteria, such as Streptococcus pneumoniae and Streptococcus pyogenes. However, little is known about the genomic distribution of ncRNAs and RNA families in S. agalactiae.

Results

Comparative genome analysis of 27 S. agalactiae strains showed more than 5 thousand genomic regions identified and classified as Core, Exclusive, and Shared genome sequences. We identified 27 to 89 RNA families per genome distributed over these regions, from these, 25 were in Core regions while Shared and Exclusive regions showed variations amongst strains. We propose that the amount and type of ncRNA present in each genome can provide a pattern to contribute in the identification of the clonal types.

Conclusions

The identification of RNA families provides an insight over ncRNAs, sRNAs and ribozymes function, that can be further explored as targets for antibiotic development or studied in gene regulation of cellular processes. RNA families could be considered as markers to determine infection capabilities of different strains. Lastly, pan-genome analysis of GBS including the full range of functional transcripts provides a broader approach in the understanding of this pathogen.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4951-z) contains supplementary material, which is available to authorized users.

High-throughput sequencing identification and characterization of potentially adhesion-related small RNAs in Streptococcus mutans

Purpose

Adherence capacity is one of the principal virulence factors of Streptococcus mutans, and adhesion virulence factors are controlled by small RNAs (sRNAs) at the post-transcriptional level in various bacteria. Here, we aimed to identify and decipher putative adhesion-related sRNAs in clinical strains of S. mutans.

Methodology

RNA deep-sequencing was performed to identify potential sRNAs under different adhesion conditions. The expression of sRNAs was analysed by quantitative real-time PCR (qRT-PCR), and bioinformatic methods were used to predict the functional characteristics of sRNAs.

Results

A total of 736 differentially expressed candidate sRNAs were predicted, and these included 352 sRNAs located on the antisense to mRNA (AM) and 384 sRNAs in intergenic regions (IGRs). The top 7 differentially expressed sRNAs were successfully validated by qRT-PCR in UA159, and 2 of these were further confirmed in 100 clinical isolates. Moreover, the sequences of two sRNAs were conserved in other Streptococcus species, indicating a conserved role in such closely related species. A good correlation between the expression of sRNAs and the adhesion of 100 clinical strains was observed, which, combined with GO and KEGG, provides a perspective for the comprehension of sRNA function annotation.

Conclusion

This study revealed a multitude of novel putative adhesion-related sRNAs in S. mutans and contributed to a better understanding of information concerning the transcriptional regulation of adhesion in S. mutans.

Fluorescence-Based Methods for Characterizing RNA Interactions In Vivo

Fluorescence-based tools that measure RNA-RNA and RNA-protein interactions in vivo offer useful experimental approaches to probe the complex and dynamic physiological behavior of bacterial RNAs. Here we document the step-by-step design and application of two fluorescence-based methods for studying the regulatory interactions RNAs perform in vivo: (i) the in vivo RNA Structural Sensing System (iRS³) for measuring RNA accessibility and (ii) the trifluorescence complementation (TriFC) assay for measuring RNA-protein interactions.

Identification of small RNAs abundant in Burkholderia cenocepacia biofilms reveal putative regulators with a potential role in carbon and iron metabolism

Small RNAs play a regulatory role in many central metabolic processes of bacteria, as well as in developmental processes such as biofilm formation. Small RNAs of Burkholderia cenocepacia, an opportunistic pathogenic beta-proteobacterium, are to date not well characterised. To address that, we performed genome-wide transcriptome structure analysis of biofilm grown B. cenocepacia J2315. 41 unannotated short transcripts were identified in intergenic regions of the B. cenocepacia genome. 15 of these short transcripts, highly abundant in biofilms, widely conserved in Burkholderia sp. and without known function, were selected for in-depth analysis. Expression profiling showed that most of these sRNAs are more abundant in biofilms than in planktonic cultures. Many are also highly abundant in cells grown in minimal media, suggesting they are involved in adaptation to nutrient limitation and growth arrest. Their computationally predicted targets include a high proportion of genes involved in carbon metabolism. Expression and target genes of one sRNA suggest a potential role in regulating iron homoeostasis. The strategy used for this study to detect sRNAs expressed in B. cenocepacia biofilms has successfully identified sRNAs with a regulatory function.

Integrative FourD omics approach profiles the target network of the carbon storage regulatory system

Multi-target regulators represent a largely untapped area for metabolic engineering and anti-bacterial development. These regulators are complex to characterize because they often act at multiple levels, affecting proteins, transcripts and metabolites. Therefore, single omics experiments cannot profile their underlying targets and mechanisms. In this work, we used an Integrative FourD omics approach (INFO) that consists of collecting and analyzing systems data throughout multiple time points, using multiple genetic backgrounds, and multiple omics approaches (transcriptomics, proteomics and high throughput sequencing crosslinking immunoprecipitation) to evaluate simultaneous changes in gene expression after imposing an environmental stress that accentuates the regulatory features of a network. Using this approach, we profiled the targets and potential regulatory mechanisms of a global regulatory system, the well-studied carbon storage regulatory (Csr) system of Escherichia coli, which is widespread among bacteria. Using 126 sets of proteomics and transcriptomics data, we identified 136 potential direct CsrA targets, including 50 novel ones, categorized their behaviors into distinct regulatory patterns, and performed in vivo fluorescence-based follow up experiments. The results of this work validate 17 novel mRNAs as authentic direct CsrA targets and demonstrate a generalizable strategy to integrate multiple lines of omics data to identify a core pool of regulator targets.

A Genome-Wide Prediction and Identification of Intergenic Small RNAs by Comparative Analysis in Mesorhizobium huakuii 7653R

In bacteria, small non-coding RNAs (sRNAs) are critical regulators of cellular adaptation to changes in metabolism, physiology, or the external environment. In the last decade, more than 2000 of sRNA families have been reported in the Rfam database and have been shown to exert various regulatory functions in bacterial transcription and translation. However, little is known about sRNAs and their functions in Mesorhizobium. Here, we predicted putative sRNAs in the intergenic regions (IGRs) of M. huakuii 7653R by genome-wide comparisons with four related Mesorhizobial strains. The expression and transcribed regions of candidate sRNAs were analyzed using a set of high-throughput RNA deep sequencing data. In all, 39 candidate sRNAs were found, with 5 located in the symbiotic megaplasmids and 34 in the chromosome of M. huakuii 7653R. Of these, 24 were annotated as functional sRNAs in the Rfam database and 15 were recognized as putative novel sRNAs. The expression of nine selected sRNAs was confirmed by Northern blotting, and most of the nine selected sRNAs were highly expressed in 28 dpi nodules and under symbiosis-mimicking conditions. For those putative novel sRNAs, functional categorizations of their target genes were performed by analyzing the enriched GO terms. In addition, MH_s15 was shown to be an abundant and conserved sRNA.

A computational strategy for the search of regulatory small RNAs in Actinobacillus pleuropneumoniae

Bacterial regulatory small RNAs (sRNAs) play important roles in gene regulation and are frequently connected to the expression of virulence factors in diverse bacteria. Only a few sRNAs have been described for Pasteurellaceae pathogens and no in-depth analysis of sRNAs has been described for Actinobacillus pleuropneumoniae, the causative agent of porcine pleuropneumonia, responsible for considerable losses in the swine industry. To search for sRNAs in A. pleuropneumoniae, we developed a strategy for the computational analysis of the bacterial genome by using four algorithms with different approaches, followed by experimental validation. The coding strand and expression of 17 out of 23 RNA candidates were confirmed by Northern blotting, RT-PCR, and RNA sequencing. Among them, two are likely riboswitches, three are housekeeping regulatory RNAs, two are the widely studied GcvB and 6S sRNAs, and 10 are putative novel trans-acting sRNAs, never before described for any bacteria. The latter group has several potential mRNA targets, many of which are involved with virulence, stress resistance, or metabolism, and connect the sRNAs in a complex gene regulatory network. The sRNAs identified are well conserved among the Pasteurellaceae that are evolutionarily closer to A. pleuropneumoniae and/or share the same host. Our results show that the combination of newly developed computational programs can be successfully utilized for the discovery of novel sRNAs and indicate an intricate system of gene regulation through sRNAs in A. pleuropneumoniae and in other Pasteurellaceae, thus providing clues for novel aspects of virulence that will be explored in further studies.

Effect of Shear Stress on Pseudomonas aeruginosa Isolated from the Cystic Fibrosis Lung

Chronic colonization of the lungs by Pseudomonas aeruginosa is one of the major causes of morbidity and mortality in cystic fibrosis (CF) patients. To gain insights into the characteristic biofilm phenotype of P. aeruginosa in the CF lungs, mimicking the CF lung environment is critical. We previously showed that growth of the non-CF-adapted P. aeruginosa PAO1 strain in a rotating wall vessel, a device that simulates the low fluid shear (LS) conditions present in the CF lung, leads to the formation of in-suspension, self-aggregating biofilms. In the present study, we determined the phenotypic and transcriptomic changes associated with the growth of a highly adapted, transmissible P. aeruginosa CF strain in artificial sputum medium under LS conditions. Robust self-aggregating biofilms were observed only under LS conditions. Growth under LS conditions resulted in the upregulation of genes involved in stress response, alginate biosynthesis, denitrification, glycine betaine biosynthesis, glycerol metabolism, and cell shape maintenance, while genes involved in phenazine biosynthesis, type VI secretion, and multidrug efflux were downregulated. In addition, a number of small RNAs appeared to be involved in the response to shear stress. Finally, quorum sensing was found to be slightly but significantly affected by shear stress, resulting in higher production of autoinducer molecules during growth under high fluid shear (HS) conditions. In summary, our study revealed a way to modulate the behavior of a highly adapted P. aeruginosa CF strain by means of introducing shear stress, driving it from a biofilm lifestyle to a more planktonic lifestyle.

Biofilm formation by Pseudomonas aeruginosa is one of the hallmarks of chronic cystic fibrosis (CF) lung infections. The biofilm matrix protects this bacterium from antibiotics as well as from the immune system. Hence, the prevention or reversion of biofilm formation is believed to have a great impact on treatment of chronic P. aeruginosa CF lung infections. In the present study, we showed that it is possible to modulate the behavior of a highly adapted transmissible P. aeruginosa CF isolate at both the transcriptomic and phenotypic levels by introducing shear stress in a CF-like environment, driving it from a biofilm to a planktonic lifestyle. Consequently, the results obtained in this study are of great importance with regard to therapeutic applications that introduce shear stress in the lungs of CF patients.

Novel small RNA (sRNA) landscape of the starvation-stress response transcriptome of Salmonella enterica serovar typhimurium

Small RNAs (sRNAs) are short (∼50-200 nucleotides) noncoding RNAs that regulate cellular activities across bacteria. Salmonella enterica starved of a carbon-energy (C) source experience a host of genetic and physiological changes broadly referred to as the starvation-stress response (SSR). In an attempt to identify novel sRNAs contributing to SSR control, we grew log-phase, 5-h C-starved and 24-h C-starved cultures of the virulent Salmonella enterica subspecies enterica serovar Typhimurium strain SL1344 and comprehensively sequenced their small RNA transcriptomes. Strikingly, after employing a novel strategy for sRNA discovery based on identifying dynamic transcripts arising from "gene-empty" regions, we identify 58 wholly undescribed Salmonella sRNA genes potentially regulating SSR averaging an ∼1,000-fold change in expression between log-phase and C-starved cells. Importantly, the expressions of individual sRNA loci were confirmed by both comprehensive transcriptome analyses and northern blotting of select candidates. Of note, we find 43 candidate sRNAs share significant sequence identity to characterized sRNAs in other bacteria, and ∼70% of our sRNAs likely assume characteristic sRNA structural conformations. In addition, we find 53 of our 58 candidate sRNAs either overlap neighboring mRNA loci or share significant sequence complementarity to mRNAs transcribed elsewhere in the SL1344 genome strongly suggesting they regulate the expression of transcripts via antisense base-pairing. Finally, in addition to this work resulting in the identification of 58 entirely novel Salmonella enterica genes likely participating in the SSR, we also find evidence suggesting that sRNAs are significantly more prevalent than currently appreciated and that Salmonella sRNAs may actually number in the thousands.

Advances in synthetic dynamic circuits design: using novel synthetic parts to engineer new generations of gene oscillations

As bioengineering applications expand, the need to design and implement circuits that exhibit dynamic properties increases. In particular, schemes that control precise patterns of gene expression as a function of time are essential for balancing multiple metabolic objectives in natural and synthetic systems. Given that modularity has been an important component of dynamic circuits, recent efforts to improve dynamic circuits have focused on replacing old parts for new components that increase the robustness, stability, and tunability. In this review, we show that incorporation of novel components such as regulatory noncoding RNAs (ncRNAs), promoter-transcription factor pairs, and metabolite sensors have allowed traditional dynamic circuits to obtain more robust functionality and improved dynamic properties.

Long non-coding RNA: A new paradigm for lung cancer

Lung cancer is the leading cause of cancer-related deaths worldwide. Recent advances in whole genome transcriptome analysis have enabled the identification of numerous members of a novel class of non-coding RNAs, i.e., long non-coding RNAs (lncRNAs), which play important roles in a wide range of biological processes and whose deregulation causes human disease, including cancer. Herein we provide a comprehensive survey of lncRNAs associated with lung cancer, with particular focus on the functions that either facilitate or inhibit the progression of lung cancer and the pathways involved. Emerging data on the use of lncRNAs as biomarkers for the diagnosis and prognosis of cancer are also discussed. We cast this information within the wider perspective of lncRNA biogenesis and molecular functions in the cell. Relationships that exist between lncRNAs, genome-wide transcription, and lung cancer are discussed. Deepening our understanding on these processes is critical not only from a mechanistic standpoint, but also for the development of novel biomarkers and effective therapeutic targets for cancer patients.

Small RNA Transcriptome of the Oral Microbiome during Periodontitis Progression

The oral microbiome is one of the most complex microbial communities in the human body, and due to circumstances not completely understood, the healthy microbial community becomes dysbiotic, giving rise to periodontitis, a polymicrobial inflammatory disease. We previously reported the results of community-wide gene expression changes in the oral microbiome during periodontitis progression and identified signatures associated with increasing severity of the disease. Small noncoding RNAs (sRNAs) are key players in posttranscriptional regulation, especially in fast-changing environments such as the oral cavity. Here, we expanded our analysis to the study of the sRNA metatranscriptome during periodontitis progression on the same samples for which mRNA expression changes were analyzed. We observed differential expression of 12,097 sRNAs, identifying a total of 20 Rfam sRNA families as being overrepresented in progression and 23 at baseline. Gene ontology activities regulated by the differentially expressed (DE) sRNAs included amino acid metabolism, ethanolamine catabolism, signal recognition particle-dependent cotranslational protein targeting to membrane, intron splicing, carbohydrate metabolism, control of plasmid copy number, and response to stress. In integrating patterns of expression of protein coding transcripts and sRNAs, we found that functional activities of genes that correlated positively with profiles of expression of DE sRNAs were involved in pathogenesis, proteolysis, ferrous iron transport, and oligopeptide transport. These findings represent the first integrated sequencing analysis of the community-wide sRNA transcriptome of the oral microbiome during periodontitis progression and show that sRNAs are key regulatory elements of the dysbiotic process leading to disease.

Transcriptional analysis of Deinococcus radiodurans reveals novel small RNAs that are differentially expressed under ionizing radiation

Small noncoding RNAs (sRNAs) are posttranscriptional regulators that have been identified in multiple species and shown to play essential roles in responsive mechanisms to environmental stresses. The natural ability of specific bacteria to resist high levels of radiation has been of high interest to mechanistic studies of DNA repair and biomolecular protection. Deinococcus radiodurans is a model extremophile for radiation studies that can survive doses of ionizing radiation of >12,000 Gy, 3,000 times higher than for most vertebrates. Few studies have investigated posttranscriptional regulatory mechanisms of this organism that could be relevant in its general gene regulatory patterns. In this study, we identified 199 potential sRNA candidates in D. radiodurans by whole-transcriptome deep sequencing analysis and confirmed the expression of 41 sRNAs by Northern blotting and reverse transcriptase PCR (RT-PCR). A total of 8 confirmed sRNAs showed differential expression during recovery after acute ionizing radiation (15 kGy). We have also found and confirmed 7 sRNAs in Deinococcus geothermalis, a closely related radioresistant species. The identification of several novel sRNAs in Deinococcus bacteria raises important questions about the evolution and nature of global gene regulation in radioresistance.