Pervasive transcription: illuminating the dark matter of bacterial transcriptomes

Large-scale investigation of species-specific orphan genes in the human gut microbiome elucidates their evolutionary origins

Species-specific genes, also known as orphans, are ubiquitous across life's domains. In prokaryotes, species-specific orphan genes (SSOGs) are mostly thought to originate in external elements such as viruses followed by horizontal gene transfer, whereas the scenario of native origination, through rapid divergence or de novo, is mostly dismissed. However, quantitative evidence supporting either scenario is lacking. Here, we systematically analyzed genomes from 4644 human gut microbiome species and identified more than 600,000 unique SSOGs, representing an average of 2.6% of a given species' pangenome. These sequences are mostly rare within each species yet show signs of purifying selection. Overall, SSOGs use optimal codons less frequently, and their proteins are more disordered than those of conserved genes (i.e., non-SSOGs). Importantly, across species, the GC content of SSOGs closely matches that of conserved ones. In contrast, the ∼5% of SSOGs that share similarity to known viral sequences have distinct characteristics, including lower GC content. Thus, SSOGs with similarity to viruses differ from the remaining SSOGs, contrasting an external origination scenario for most of them. By examining the orthologous genomic region in closely related species, we show that a small subset of SSOGs likely evolved natively de novo and find that these genes also differ in their properties from the remaining SSOGs. Our results challenge the notion that external elements are the dominant source of prokaryotic genetic novelty and will enable future studies into the biological role and relevance of species-specific genes in the human gut.

Transcriptome analysis of Edwardsiella piscicida during intracellular infection reveals excludons are involved with the activation of a mitochondrion-like energy generation program

Phylogenetic evidence suggests a shared ancestry between mitochondria and modern Proteobacteria, a phylum including several genera of intracellular pathogens. Studying these diverse pathogens, particularly during intracellular infection of their hosts, can reveal characteristics potentially representative of the mitochondrial-Proteobacterial ancestor by identifying traits shared with mitochondria. While transcriptomic approaches can provide global insights into intracellular acclimatization by pathogens, they are often limited by excess host RNAs in extracts. Here, we developed a method employing magnetic nanoparticles to enrich RNA from an intracellular Gammaproteobacterium, Edwardsiella piscicida, within zebrafish, Danio rerio, fin fibroblasts, enabling comprehensive exploration of the bacterial transcriptome. Our findings revealed that the intracellular E. piscicida transcriptome reflects a mitochondrion-like energy generation program characterized by the suppression of glycolysis and sugar transport, coupled with upregulation of the tricarboxylic acid (TCA) cycle and alternative import of simple organic acids that directly flux into TCA cycle intermediates or electron transport chain donors. Additionally, genes predicted to be members of excludons, loci of gene pairs antagonistically co-regulated by overlapping antisense transcription, are significantly enriched in the set of all genes with perturbed sense and antisense transcription, suggesting a general but important involvement of excludons with intracellular acclimatization. Notably, genes involved with the activation of the mitochondrion-like energy generation program, specifically with metabolite import and glycolysis, are also members of predicted excludons. Other intracellular Proteobacterial pathogens appear to employ a similar mitochondrion-like energy generation program, suggesting a potentially conserved mechanism for optimized energy acquisition from hosts centered around the TCA cycle.IMPORTANCEPhylogenetic evidence suggests that mitochondria and Proteobacteria, a phylum encompassing various intracellular pathogens, share a common ancestral lineage. In this study, we developed a novel method employing magnetic nanoparticles to explore the transcriptome of an aquatic Gammaproteobacterium, Edwardsiella piscicida, during intracellular infection of host cells. We show that the strategy E. piscicida uses to generate energy strikingly mirrors the function of mitochondria-energy generators devoid of glycolytic processes. Notably, several implicated genes are members of excludons-gene pairs antagonistically co-regulated by overlapping antisense transcription. Other intracellular Proteobacterial pathogens appear to adopt a similar mitochondrion-like energy generation program, indicating a possibly conserved strategy for optimized energy acquisition from hosts centered around the tricarboxylic acid cycle.

Rho-dependent termination enables cellular pH homeostasis

The termination factor Rho, an ATP-dependent RNA translocase, preempts pervasive transcription processes, thereby rendering genome integrity in bacteria. Here, we show that the loss of Rho function raised the intracellular pH to >8.0 in Escherichia coli. The loss of Rho function upregulates tryptophanase-A (TnaA), an enzyme that catabolizes tryptophan to produce indole, pyruvate, and ammonia. We demonstrate that the enhanced TnaA function had produced the conjugate base ammonia, raising the cellular pH in the Rho-dependent termination defective strains. On the other hand, the constitutively overexpressed Rho lowered the cellular pH to about 6.2, independent of cellular ammonia levels. Since Rho overexpression may increase termination activities, the decrease in cellular pH could result from an excess H+ ion production during ATP hydrolysis by overproduced Rho. Furthermore, we performed in vivo termination assays to show that the efficiency of Rho-dependent termination was increased at both acidic and basic pH ranges. Given that the Rho level remained unchanged, the alkaline pH increases the termination efficiency by stimulating Rho's catalytic activity. We conducted the Rho-mediated RNA release assay from a stalled elongation complex to show an efficient RNA release at alkaline pH, compared to the neutral or acidic pH, that supports our in vivo observation. Whereas acidic pH appeared to increase the termination function by elevating the cellular level of Rho. This study is the first to link Rho function to the cellular pH homeostasis in bacteria. IMPORTANCE The current study shows that the loss or gain of Rho-dependent termination alkalizes or acidifies the cytoplasm, respectively. In the case of loss of Rho function, the tryptophanase-A enzyme is upregulated, and degrades tryptophan, producing ammonia to alkalize cytoplasm. We hypothesize that Rho overproduction by deleting its autoregulatory DNA portion increases termination function, causing excessive ATP hydrolysis to produce H+ ions and cytoplasmic acidification. Therefore, this study is the first to unravel a relationship between Rho function and intrinsic cellular pH homeostasis. Furthermore, the Rho level increases in the absence of autoregulation, causing cytoplasmic acidification. As intracellular pH plays a critical role in enzyme function, such a connection between Rho function and alkalization will have far-reaching implications for bacterial physiology.

Quantification and modeling of turnover dynamics of de novo transcripts in Drosophila melanogaster

Most of the transcribed eukaryotic genomes are composed of non-coding transcripts. Among these transcripts, some are newly transcribed when compared to outgroups and are referred to as de novo transcripts. De novo transcripts have been shown to play a major role in genomic innovations. However, little is known about the rates at which de novo transcripts are gained and lost in individuals of the same species. Here, we address this gap and estimate the de novo transcript turnover rate with an evolutionary model. We use DNA long reads and RNA short reads from seven geographically remote samples of inbred individuals of Drosophila melanogaster to detect de novo transcripts that are gained on a short evolutionary time scale. Overall, each sampled individual contains around 2500 unspliced de novo transcripts, with most of them being sample specific. We estimate that around 0.15 transcripts are gained per year, and that each gained transcript is lost at a rate around 5× 10-5 per year. This high turnover of transcripts suggests frequent exploration of new genomic sequences within species. These rate estimates are essential to comprehend the process and timescale of de novo gene birth.

Metabolic and transcriptional activities underlie stationary-phase Pseudomonas aeruginosa sensitivity to Levofloxacin

IMPORTANCE

The bacterial pathogen Pseudomonas aeruginosa is responsible for a variety of chronic human infections. Even in the absence of identifiable resistance mutations, this pathogen can tolerate lethal antibiotic doses through phenotypic strategies like biofilm formation and metabolic quiescence. In this study, we determined that P. aeruginosa maintains greater metabolic activity in the stationary phase compared to the model organism, Escherichia coli, which has traditionally been used to study fluoroquinolone antibiotic tolerance. We demonstrate that hallmarks of E. coli fluoroquinolone tolerance are not conserved in P. aeruginosa, including the timing of cell death and necessity of the SOS DNA damage response for survival. The heightened sensitivity of stationary-phase P. aeruginosa to fluoroquinolones is attributed to maintained transcriptional and reductase activity. Our data suggest that perturbations that suppress transcription and respiration in P. aeruginosa may actually protect the pathogen against this important class of antibiotics.

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.

Comparative Assessment of Pose Prediction Accuracy in RNA-Ligand Docking

Structure-based virtual high-throughput screening is used in early-stage drug discovery. Over the years, docking protocols and scoring functions for protein-ligand complexes have evolved to improve the accuracy in the computation of binding strengths and poses. In the past decade, RNA has also emerged as a target class for new small-molecule drugs. However, most ligand docking programs have been validated and tested for proteins and not RNA. Here, we test the docking power (pose prediction accuracy) of three state-of-the-art docking protocols on 173 RNA-small molecule crystal structures. The programs are AutoDock4 (AD4) and AutoDock Vina (Vina), which were designed for protein targets, and rDock, which was designed for both protein and nucleic acid targets. AD4 performed relatively poorly. For RNA targets for which a crystal structure of a bound ligand used to limit the docking search space is available and for which the goal is to identify new molecules for the same pocket, rDock performs slightly better than Vina, with success rates of 48% and 63%, respectively. However, in the more common type of early-stage drug discovery setting, in which no structure of a ligand-target complex is known and for which a larger search space is defined, rDock performed similarly to Vina, with a low success rate of ∼27%. Vina was found to have bias for ligands with certain physicochemical properties, whereas rDock performs similarly for all ligand properties. Thus, for projects where no ligand-protein structure already exists, Vina and rDock are both applicable. However, the relatively poor performance of all methods relative to protein-target docking illustrates a need for further methods refinement.

DiffSegR: an RNA-seq data driven method for differential expression analysis using changepoint detection

To fully understand gene regulation, it is necessary to have a thorough understanding of both the transcriptome and the enzymatic and RNA-binding activities that shape it. While many RNA-Seq-based tools have been developed to analyze the transcriptome, most only consider the abundance of sequencing reads along annotated patterns (such as genes). These annotations are typically incomplete, leading to errors in the differential expression analysis. To address this issue, we present DiffSegR - an R package that enables the discovery of transcriptome-wide expression differences between two biological conditions using RNA-Seq data. DiffSegR does not require prior annotation and uses a multiple changepoints detection algorithm to identify the boundaries of differentially expressed regions in the per-base log2 fold change. In a few minutes of computation, DiffSegR could rightfully predict the role of chloroplast ribonuclease Mini-III in rRNA maturation and chloroplast ribonuclease PNPase in (3'/5')-degradation of rRNA, mRNA and tRNA precursors as well as intron accumulation. We believe DiffSegR will benefit biologists working on transcriptomics as it allows access to information from a layer of the transcriptome overlooked by the classical differential expression analysis pipelines widely used today. DiffSegR is available at https://aliehrmann.github.io/DiffSegR/index.html.

Non-canonical transcriptional start sites in E. coli O157:H7 EDL933 are regulated and appear in surprisingly high numbers

Analysis of genome wide transcription start sites (TSSs) revealed an unexpected complexity since not only canonical TSS of annotated genes are recognized by RNA polymerase. Non-canonical TSS were detected antisense to, or within, annotated genes as well new intergenic (orphan) TSS, not associated with known genes. Previously, it was hypothesized that many such signals represent noise or pervasive transcription, not associated with a biological function. Here, a modified Cappable-seq protocol allows determining the primary transcriptome of the enterohemorrhagic E. coli O157:H7 EDL933 (EHEC). We used four different growth media, both in exponential and stationary growth phase, replicated each thrice. This yielded 19,975 EHEC canonical and non-canonical TSS, which reproducibly occurring in three biological replicates. This questions the hypothesis of experimental noise or pervasive transcription. Accordingly, conserved promoter motifs were found upstream indicating proper TSSs. More than 50% of 5,567 canonical and between 32% and 47% of 10,355 non-canonical TSS were differentially expressed in different media and growth phases, providing evidence for a potential biological function also of non-canonical TSS. Thus, reproducible and environmentally regulated expression suggests that a substantial number of the non-canonical TSSs may be of unknown function rather than being the result of noise or pervasive transcription.

Comprehensive transcriptional analysis of pig facial skin development

Background

Skin development is a complex process that is influenced by many factors. Pig skin is used as an ideal material for xenografts because it is more anatomically and physiologically similar to human skin. It has been shown that the skin development of different pig breeds is different, and some Chinese pig breeds have the characteristics of skin thickness and facial skin folds, but the specific regulatory mechanism of this skin development is not yet clear.

Methods

In this study, the facial skin of Chenghua sows in the four developmental stages of postnatal Day 3 (D3) , Day 90 (D90) , Day 180 (D180), and Year 3 (Y3) were used as experimental materials, and RNA sequencing (RNA-seq) analysis was used to explore the changes in RNA expression in skin development at the four developmental stages, determine the differentially expressed messenger RNAs (mRNAs), long noncoding RNAs (lncRNAs), microRNAs (miRNAs), and circular RNAs (circRNAs), and perform functional analysis of related genes by Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses.

Results

A pairwise comparison of the four developmental stages identified several differentially expressed genes (DEGs) and found that the number of differentially expressed RNAs (DE RNAs) increased with increasing developmental time intervals. Elastin (ELN) is an important component of the skin. Its content affects the relaxation of the epidermis and dermal connection, and its expression is continuously downregulated during the four developmental stages. The functions of DEGs at different developmental stages were examined by performing GO and KEGG analyses, and the GO terms and enrichment pathways of mRNAs, lncRNAs, miRNAs, and circRNAs highly overlapped, among which the PPAR signaling pathway, a classical pathway for skin development, was enriched by DEGs of D3 vs. D180, D90 vs. D180 and D180 vs. Y3. In addition, we constructed lncRNA-miRNA-mRNA and circRNA-miRNA interaction networks and found genes that may be associated with skin development, but their interactions need further study.

Conclusions

We identified a number of genes associated with skin development, performed functional analyses on some important DEGs and constructed interaction networks that facilitate further studies of skin development.

Genome-Wide Mapping of the Escherichia coli PhoB Regulon Reveals Many Transcriptionally Inert, Intragenic Binding Sites

Genome-scale analyses have revealed many transcription factor binding sites within, rather than upstream of, genes, raising questions as to the function of these binding sites. Here, we use complementary approaches to map the regulon of the Escherichia coli transcription factor PhoB, a response regulator that controls transcription of genes involved in phosphate homeostasis. Strikingly, the majority of PhoB binding sites are located within genes, but these intragenic sites are not associated with detectable transcription regulation and are not evolutionarily conserved. Many intragenic PhoB sites are located in regions bound by H-NS, likely due to shared sequence preferences of PhoB and H-NS. However, these PhoB binding sites are not associated with transcription regulation even in the absence of H-NS. We propose that for many transcription factors, including PhoB, binding sites not associated with promoter sequences are transcriptionally inert and hence are tolerated as genomic "noise." IMPORTANCE Recent studies have revealed large numbers of transcription factor binding sites within the genes of bacteria. The function, if any, of the vast majority of these binding sites has not been investigated. Here, we map the binding of the transcription factor PhoB across the Escherichia coli genome, revealing that the majority of PhoB binding sites are within genes. We show that PhoB binding sites within genes are not associated with regulation of the overlapping genes. Indeed, our data suggest that bacteria tolerate the presence of large numbers of nonregulatory, intragenic binding sites for transcription factors and that these binding sites are not under selective pressure.

In Vivo Sampling of Intracellular Heterogeneity of Pseudomonas putida Enables Multiobjective Optimization of Genetic Devices

The inner physicochemical heterogeneity of bacterial cells generates three-dimensional (3D)-dependent variations of resources for effective expression of given chromosomally located genes. This fact has been exploited for adjusting the most favorable parameters for implanting a complex device for optogenetic control of biofilm formation in the soil bacterium Pseudomonas putida. To this end, a DNA segment encoding a superactive variant of the Caulobacter crescendus diguanylate cyclase PleD expressed under the control of the cyanobacterial light-responsive CcaSR system was placed in a mini-Tn5 transposon vector and randomly inserted through the chromosome of wild-type and biofilm-deficient variants of P. putida lacking the wsp gene cluster. This operation delivered a collection of clones covering a whole range of biofilm-building capacities and dynamic ranges in response to green light. Since the phenotypic output of the device depends on a large number of parameters (multiple promoters, RNA stability, translational efficacy, metabolic precursors, protein folding, etc.), we argue that random chromosomal insertions enable sampling the intracellular milieu for an optimal set of resources that deliver a preset phenotypic specification. Results thus support the notion that the context dependency can be exploited as a tool for multiobjective optimization, rather than a foe to be suppressed in Synthetic Biology constructs.

Population genomics reveals mechanisms and dynamics of de novo expressed open reading frame emergence in Drosophila melanogaster

Novel genes are essential for evolutionary innovations and differ substantially even between closely related species. Recently, multiple studies across many taxa showed that some novel genes arise de novo, that is, from previously noncoding DNA. To characterize the underlying mutations that allowed de novo gene emergence and their order of occurrence, homologous regions must be detected within noncoding sequences in closely related sister genomes. So far, most studies do not detect noncoding homologs of de novo genes because of incomplete assemblies and annotations, and long evolutionary distances separating genomes. Here, we overcome these issues by searching for de novo expressed open reading frames (neORFs), the not-yet fixed precursors of de novo genes that emerged within a single species. We sequenced and assembled genomes with long-read technology and the corresponding transcriptomes from inbred lines of Drosophila melanogaster, derived from seven geographically diverse populations. We found line-specific neORFs in abundance but few neORFs shared by lines, suggesting a rapid turnover. Gain and loss of transcription is more frequent than the creation of ORFs, for example, by forming new start and stop codons. Consequently, the gain of ORFs becomes rate limiting and is frequently the initial step in neORFs emergence. Furthermore, transposable elements (TEs) are major drivers for intragenomic duplications of neORFs, yet TE insertions are less important for the emergence of neORFs. However, highly mutable genomic regions around TEs provide new features that enable gene birth. In conclusion, neORFs have a high birth-death rate, are rapidly purged, but surviving neORFs spread neutrally through populations and within genomes.

RNA polymerase drives ribonucleotide excision DNA repair in E. coli

Ribonuclease HII (RNaseHII) is the principal enzyme that removes misincorporated ribonucleoside monophosphates (rNMPs) from genomic DNA. Here, we present structural, biochemical, and genetic evidence demonstrating that ribonucleotide excision repair (RER) is directly coupled to transcription. Affinity pull-downs and mass-spectrometry-assisted mapping of in cellulo inter-protein cross-linking reveal the majority of RNaseHII molecules interacting with RNA polymerase (RNAP) in E. coli. Cryoelectron microscopy structures of RNaseHII bound to RNAP during elongation, with and without the target rNMP substrate, show specific protein-protein interactions that define the transcription-coupled RER (TC-RER) complex in engaged and unengaged states. The weakening of RNAP-RNaseHII interactions compromises RER in vivo. The structure-functional data support a model where RNaseHII scans DNA in one dimension in search for rNMPs while "riding" the RNAP. We further demonstrate that TC-RER accounts for a significant fraction of repair events, thereby establishing RNAP as a surveillance "vehicle" for detecting the most frequently occurring replication errors.

An unexpected abundance of bidirectional promoters within Salmonella Typhimurium plasmids

Transcription of the DNA template, to generate an RNA message, is the first step in gene expression. The process initiates at DNA sequences called promoters. Conventionally, promoters have been considered to drive transcription in a specific direction. However, in recent work, we showed that many prokaryotic promoters can drive divergent transcription. This is a consequence of key DNA sequences for transcription initiation being inherently symmetrical. Here, we used global transcription start site mapping to determine the prevalence of such bidirectional promoters in Salmonella Typhimurium. Surprisingly, bidirectional promoters occur three times more frequently in plasmid components of the genome compared to chromosomal DNA. Implications for the evolution of promoter sequences are discussed.

Non-coding RNAs and Exosomal Non-coding RNAs in Traumatic Brain Injury: the Small Player with Big Actions

Nowadays, there is an increasing concern regarding traumatic brain injury (TBI) worldwide since substantial morbidity is observed after it, and the long-term consequences that are not yet fully recognized. A number of cellular pathways related to the secondary injury in brain have been identified, including free radical production (owing to mitochondrial dysfunction), excitotoxicity (regulated by excitatory neurotransmitters), apoptosis, and neuroinflammatory responses (as a result of activation of the immune system and central nervous system). In this context, non-coding RNAs (ncRNAs) maintain a fundamental contribution to post-transcriptional regulation. It has been shown that mammalian brains express high levels of ncRNAs that are involved in several brain physiological processes. Furthermore, altered levels of ncRNA expression have been found in those with traumatic as well non-traumatic brain injuries. The current review highlights the primary molecular mechanisms participated in TBI that describes the latest and novel results about changes and role of ncRNAs in TBI in both clinical and experimental research.

Using a whole genome co-expression network to inform the functional characterisation of predicted genomic elements from Mycobacterium tuberculosis transcriptomic data

A whole genome co-expression network was created using Mycobacterium tuberculosis transcriptomic data from publicly available RNA-sequencing experiments covering a wide variety of experimental conditions. The network includes expressed regions with no formal annotation, including putative short RNAs and untranslated regions of expressed transcripts, along with the protein-coding genes. These unannotated expressed transcripts were among the best-connected members of the module sub-networks, making up more than half of the 'hub' elements in modules that include protein-coding genes known to be part of regulatory systems involved in stress response and host adaptation. This data set provides a valuable resource for investigating the role of non-coding RNA, and conserved hypothetical proteins, in transcriptomic remodelling. Based on their connections to genes with known functional groupings and correlations with replicated host conditions, predicted expressed transcripts can be screened as suitable candidates for further experimental validation.

Genome-wide mapping of the Escherichia coli PhoB regulon reveals many transcriptionally inert, intragenic binding sites

Genome-scale analyses have revealed many transcription factor binding sites within, rather than upstream of genes, raising questions as to the function of these binding sites. Here, we use complementary approaches to map the regulon of the Escherichia coli transcription factor PhoB, a response regulator that controls transcription of genes involved in phosphate homeostasis. Strikingly, the majority of PhoB binding sites are located within genes, but these intragenic sites are not associated with detectable transcription regulation and are not evolutionarily conserved. Many intragenic PhoB sites are located in regions bound by H-NS, likely due to shared sequence preferences of PhoB and H-NS. However, these PhoB binding sites are not associated with transcription regulation even in the absence of H-NS. We propose that for many transcription factors, including PhoB, binding sites not associated with promoter sequences are transcriptionally inert, and hence are tolerated as genomic "noise".

IMPORTANCE

Recent studies have revealed large numbers of transcription factor binding sites within the genes of bacteria. The function, if any, of the vast majority of these binding sites has not been investigated. Here, we map the binding of the transcription factor PhoB across the Escherichia coli genome, revealing that the majority of PhoB binding sites are within genes. We show that PhoB binding sites within genes are not associated with regulation of the overlapping genes. Indeed, our data suggest that bacteria tolerate the presence of large numbers of non-regulatory, intragenic binding sites for transcription factors, and that these binding sites are not under selective pressure.

Prediction of protein-coding small ORFs in multi-species using integrated sequence-derived features and the random forest model

Proteins encoded by small open reading frames (sORFs) can serve as functional elements playing important roles in vivo. Such sORFs also constitute the potential pool for facilitating the de novo gene birth, driving evolutionary innovation and species diversity. Therefore, their theoretical and experimental identification has become a critical issue. Herein, we proposed a protein-coding sORFs prediction method merely based on integrative sequence-derived features. Our prediction performance is better or comparable compared with other nine prevalent methods, which shows that our method can provide a relatively reliable research tool for the prediction of protein-coding sORFs. Our method allows users to estimate the potential expression of a queried sORF, which has been demonstrated by the correlation analysis between our possibility estimation and codon adaption index (CAI). Based on the features that we used, we demonstrated that the sequence features of the protein-coding sORFs in the two domains have significant differences implying that it might be a relatively hard task in terms of cross-domain prediction, hence domain-specific models were developed, which allowed users to predict protein-coding sORFs both in eukaryotes and prokaryotes. Finally, a web-server was developed and provided to boost and facilitate the study of the related field, which is freely available at http://guolab.whu.edu.cn/codingCapacity/index.html.

The RNA-binding protein RBP33 dampens non-productive transcription in trypanosomes

In-depth analysis of the transcriptomes of several model organisms has revealed that genomes are pervasively transcribed, giving rise to an abundance of non-canonical and mainly antisense RNA polymerase II-derived transcripts that are produced from almost any genomic context. Pervasive RNAs are degraded by surveillance mechanisms, but the repertoire of proteins that control the fate of these non-productive transcripts is still incomplete. Trypanosomes are single-celled eukaryotes that show constitutive RNA polymerase II transcription and in which initiation and termination of transcription occur at a limited number of sites per chromosome. It is not known whether pervasive transcription exists in organisms with unregulated RNA polymerase II activity, and which factors could be involved in the process. We show here that depletion of RBP33 results in overexpression of ∼40% of all annotated genes in the genome, with a marked accumulation of sense and antisense transcripts derived from silenced regions. RBP33 loss does not result in a significant increase in chromatin accessibility. Finally, we have found that transcripts that increase in abundance upon RBP33 knockdown are significantly more stable in RBP33-depleted trypanosomes, and that the exosome complex is responsible for their degradation. Our results provide strong evidence that RBP33 dampens non-productive transcription in trypanosomes.

Faria N, Mwangi M, Miragaia M, de Lencastre H, Tomasz A, Gonçalves Sobral R. Analysis of a Cell Wall Mutant Highlights Rho-Dependent Genome Amplification Events in Staphylococcus aureus

In a study of antibiotic resistance in Staphylococcus aureus, specific cell wall mutants were previously generated for the peptidoglycan biosynthesis gene murF, by the insertion of an integrative plasmid. A collection of 30 independent mutants was obtained, and all harbored a variable number of copies of the inserted plasmid, arranged in tandem in the chromosome. Of the 30 mutants, only 3, F9, F20 and F26, with a lower number of plasmid copies, showed an altered peptidoglycan structure, lower resistance to β-lactams and a different loss-of-function mutation in rho gene, that encodes a transcription termination factor. The rho mutations were found to correlate with the level of oxacillin resistance, since genetic complementation with rho gene reestablished the resistance and cell wall parental profile in F9, F20 and F26 strains. Furthermore, complementation with rho resulted in the amplification of the number of plasmid tandem repeats, suggesting that Rho enabled events of recombination that favored a rearrangement in the chromosome in the region of the impaired murF gene. Although the full mechanism of reversion of the cell wall damage was not fully elucidated, we showed that Rho is involved in the recombination process that mediates the tandem amplification of exogeneous DNA fragments inserted into the chromosome. IMPORTANCE The cell wall of bacteria, namely, peptidoglycan, is the target of several antibiotic classes such as β-lactams. Staphylococcus aureus is well known for its capacity to adapt to antibiotic stress and develop resistance strategies, namely, to β-lactams. In this context, the construction of cell wall mutants provides useful models to study the development of such resistance mechanisms. Here, we characterized a collection of independent mutants, impaired in the same peptidoglycan biosynthetic step, obtained through the insertion of a plasmid in the coding region of murF gene. S. aureus demonstrated the capacity to overcome the cell wall damage by amplifying the copy number of the inserted plasmid, through an undescribed mechanism that involves the Rho transcription termination factor.

Genome-wide analysis suggests the potential role of lncRNAs during seed development and seed size/weight determination in chickpea

The integrated transcriptome data analyses suggested the plausible roles of lncRNAs during seed development in chickpea. The candidate lncRNAs associated with QTLs and those involved in miRNA-mediated seed size/weight determination in chickpea have been identified. Long non-coding RNAs (lncRNAs) are important regulators of various biological processes. Here, we identified lncRNAs at seven successive stages of seed development in small-seeded and large-seeded chickpea cultivars. In total, 4751 lncRNAs implicated in diverse biological processes were identified. Most of lncRNAs were conserved between the two cultivars, whereas only a few of them were conserved in other plants, suggesting their species-specificity. A large number of lncRNAs differentially expressed between the two chickpea cultivars associated with seed development-related processes were identified. The lncRNAs acting as precursors of miRNAs and those mimicking target protein-coding genes of miRNAs involved in seed size/weight determination, including HAIKU1, BIG SEEDS1, and SHB1, were also revealed. Further, lncRNAs located within seed size/weight associated quantitative trait loci were also detected. Overall, we present a comprehensive resource and identified candidate lncRNAs that may play important roles during seed development and seed size/weight determination in chickpea.

Clustered regularly interspaced short palindromic repeats-Cas system: diversity and regulation in Enterobacteriaceae

Insights into the arms race between bacteria and invading mobile genetic elements have revealed the intricacies of the clustered regularly interspaced short palindromic repeats (CRISPR)-Cas system and the counter-defenses of bacteriophages. Incredible spacer diversity but significant spacer conservation among species/subspecies dictates the specificity of the CRISPR-Cas system. Researchers have exploited this feature to type/subtype the bacterial strains, devise targeted antimicrobials and regulate gene expression. This review focuses on the nuances of the CRISPR-Cas systems in Enterobacteriaceae that predominantly harbor type I-E and I-F CRISPR systems. We discuss the systems' regulation by the global regulators, H-NS, LeuO, LRP, cAMP receptor protein and other regulators in response to environmental stress. We further discuss the regulation of noncanonical functions like DNA repair pathways, biofilm formation, quorum sensing and virulence by the CRISPR-Cas system. The review comprehends multiple facets of the CRISPR-Cas system in Enterobacteriaceae including its diverse attributes, association with genetic features, regulation and gene regulatory mechanisms.

Regulatory Interplay between RNase III and Antisense RNAs in E. coli: the Case of AsflhD and FlhD, Component of the Master Regulator of Motility

In order to respond to ever-changing environmental cues, bacteria display resilient regulatory mechanisms controlling gene expression. At the post-transcriptional level, this is achieved by a combination of RNA-binding proteins, such as ribonucleases (RNases), and regulatory RNAs, including antisense RNAs (asRNAs). Bound to their complementary mRNA, asRNAs are primary targets for the double-strand-specific endoribonuclease, RNase III. Taking advantage of our own and previously published transcriptomic data sets obtained in strains inactivated for RNase III, we selected several candidate asRNAs and confirmed the existence of RNase III-sensitive asRNAs for crp, ompR, phoP, and flhD genes, all encoding global regulators of gene expression in Escherichia coli. Using FlhD, a component of the master regulator of motility (FlhD₄C₂), as our model, we demonstrate that the asRNA AsflhD, transcribed from the coding sequence of flhD, is involved in the fine-tuning of flhD expression and thus participates in the control of motility.

Pervasive transcription enhances the accessibility of H-NS-silenced promoters and generates bistability in Salmonella virulence gene expression

In Escherichia coli and Salmonella, many genes silenced by the nucleoid structuring protein H-NS are activated upon inhibiting Rho-dependent transcription termination. This response is poorly understood and difficult to reconcile with the view that H-NS acts mainly by blocking transcription initiation. Here we have analyzed the basis for the up-regulation of H-NS-silenced Salmonella pathogenicity island 1 (SPI-1) in cells depleted of Rho-cofactor NusG. Evidence from genetic experiments, semiquantitative 5' rapid amplification of complementary DNA ends sequencing (5' RACE-Seq), and chromatin immunoprecipitation sequencing (ChIP-Seq) shows that transcription originating from spurious antisense promoters, when not stopped by Rho, elongates into a H-NS-bound regulatory region of SPI-1, displacing H-NS and rendering the DNA accessible to the master regulator HilD. In turn, HilD's ability to activate its own transcription triggers a positive feedback loop that results in transcriptional activation of the entire SPI-1. Significantly, single-cell analyses revealed that this mechanism is largely responsible for the coexistence of two subpopulations of cells that either express or do not express SPI-1 genes. We propose that cell-to-cell differences produced by stochastic spurious transcription, combined with feedback loops that perpetuate the activated state, can generate bimodal gene expression patterns in bacterial populations.

Pervasive translation of circular RNAs driven by short IRES-like elements

Some circular RNAs (circRNAs) were found to be translated through IRES-driven mechanism, however the scope and functions of circRNA translation are unclear because endogenous IRESs are rare. To determine the prevalence and mechanism of circRNA translation, we develop a cell-based system to screen random sequences and identify 97 overrepresented hexamers that drive cap-independent circRNA translation. These IRES-like short elements are significantly enriched in endogenous circRNAs and sufficient to drive circRNA translation. We further identify multiple trans-acting factors that bind these IRES-like elements to initiate translation. Using mass-spectrometry data, hundreds of circRNA-coded peptides are identified, most of which have low abundance due to rapid degradation. As judged by mass-spectrometry, 50% of translatable endogenous circRNAs undergo rolling circle translation, several of which are experimentally validated. Consistently, mutations of the IRES-like element in one circRNA reduce its translation. Collectively, our findings suggest a pervasive translation of circRNAs, providing profound implications in translation control.

Unbiased screen of random sequences identified many short IRES-like elements to drive circular RNA translation and hundreds of rolling circle translation events, suggesting a pervasive cap-independent translation in human transcriptome.

In-Depth Temporal Transcriptome Profiling of an Alphaherpesvirus Using Nanopore Sequencing

In this work, a long-read sequencing (LRS) technique based on the Oxford Nanopore Technology MinION platform was used for quantifying and kinetic characterization of the poly(A) fraction of bovine alphaherpesvirus type 1 (BoHV-1) lytic transcriptome across a 12-h infection period. Amplification-based LRS techniques frequently generate artefactual transcription reads and are biased towards the production of shorter amplicons. To avoid these undesired effects, we applied direct cDNA sequencing, an amplification-free technique. Here, we show that a single promoter can produce multiple transcription start sites whose distribution patterns differ among the viral genes but are similar in the same gene at different timepoints. Our investigations revealed that the circ gene is expressed with immediate–early (IE) kinetics by utilizing a special mechanism based on the use of the promoter of another IE gene (bicp4) for the transcriptional control. Furthermore, we detected an overlap between the initiation of DNA replication and the transcription from the bicp22 gene, which suggests an interaction between the two molecular machineries. This study developed a generally applicable LRS-based method for the time-course characterization of transcriptomes of any organism.

The Pet127 protein is a mitochondrial 5'-to-3' exoribonuclease from the PD-(D/E)XK superfamily involved in RNA maturation and intron degradation in yeasts

Pet127 is a mitochondrial protein found in multiple eukaryotic lineages, but absent from several taxa, including plants and animals. Distant homology suggests that it belongs to the divergent PD-(D/E)XK superfamily which includes various nucleases and related proteins. Earlier yeast genetics experiments suggest that it plays a nonessential role in RNA degradation and 5' end processing. Our phylogenetic analysis suggests that it is a primordial eukaryotic invention that was retained in diverse groups, and independently lost several times in the evolution of other organisms. We demonstrate for the first time that the fungal Pet127 protein in vitro is a processive 5'-to-3' exoribonuclease capable of digesting various substrates in a sequence nonspecific manner. Mutations in conserved residues essential in the PD-(D/E)XK superfamily active site abolish the activity of Pet127. Deletion of the PET127 gene in the pathogenic yeast Candida albicans results in a moderate increase in the steady-state levels of several transcripts and in accumulation of unspliced precursors and intronic sequences of three introns. Mutations in the active site residues result in a phenotype identical to that of the deletant, confirming that the exoribonuclease activity is related to the physiological role of the Pet127 protein. Pet127 activity is, however, not essential for maintaining the mitochondrial respiratory activity in C. albicans.

Crucial role and mechanism of transcription-coupled DNA repair in bacteria

Transcription-coupled DNA repair (TCR) is presumed to be a minor sub-pathway of nucleotide excision repair (NER) in bacteria. Global genomic repair is thought to perform the bulk of repair independently of transcription. TCR is also believed to be mediated exclusively by Mfd-a DNA translocase of a marginal NER phenotype^1-3. Here we combined in cellulo cross-linking mass spectrometry with structural, biochemical and genetic approaches to map the interactions within the TCR complex (TCRC) and to determine the actual sequence of events that leads to NER in vivo. We show that RNA polymerase (RNAP) serves as the primary sensor of DNA damage and acts as a platform for the recruitment of NER enzymes. UvrA and UvrD associate with RNAP continuously, forming a surveillance pre-TCRC. In response to DNA damage, pre-TCRC recruits a second UvrD monomer to form a helicase-competent UvrD dimer that promotes backtracking of the TCRC. The weakening of UvrD-RNAP interactions renders cells sensitive to genotoxic stress. TCRC then recruits a second UvrA molecule and UvrB to initiate the repair process. Contrary to the conventional view, we show that TCR accounts for the vast majority of chromosomal repair events; that is, TCR thoroughly dominates over global genomic repair. We also show that TCR is largely independent of Mfd. We propose that Mfd has an indirect role in this process: it participates in removing obstructive RNAPs in front of TCRCs and also in recovering TCRCs from backtracking after repair has been completed.

Pervasive Transcription-coupled DNA repair in E. coli

Global Genomic Repair (GGR) and Transcription-Coupled Repair (TCR) have been viewed, respectively, as major and minor sub-pathways of the nucleotide excision repair (NER) process that removes bulky lesions from the genome. Here we applied a next generation sequencing assay, CPD-seq, in E. coli to measure the levels of cyclobutane pyrimidine dimer (CPD) lesions before, during, and after UV-induced genotoxic stress, and, therefore, to determine the rate of genomic recovery by NER at a single nucleotide resolution. We find that active transcription is necessary for the repair of not only the template strand (TS), but also the non-template strand (NTS), and that the bulk of TCR is independent of Mfd – a DNA translocase that is thought to be necessary and sufficient for TCR in bacteria. We further show that repair of both TS and NTS is enhanced by increased readthrough past Rho-dependent terminators. We demonstrate that UV-induced genotoxic stress promotes global antitermination so that TCR is more accessible to the antisense, intergenic, and other low transcribed regions. Overall, our data suggest that GGR and TCR are essentially the same process required for complete repair of the bacterial genome.

Transcription-Coupled DNA repair has been classically defined as the preferential repair of the template strand (TS) over the non-template strand (NTS). Here the authors challenge this classic model of TCR by using a genome-wide repair assay, CPD-seq, as well as RNA-seq, to show that TCR occurs across the entire E. coli genome – including NTS and intergenic regions.

Pervasive translation in Mycobacterium tuberculosis

Most bacterial ORFs are identified by automated prediction algorithms. However, these algorithms often fail to identify ORFs lacking canonical features such as a length of >50 codons or the presence of an upstream Shine-Dalgarno sequence. Here, we use ribosome profiling approaches to identify actively translated ORFs in Mycobacterium tuberculosis. Most of the ORFs we identify have not been previously described, indicating that the M. tuberculosis transcriptome is pervasively translated. The newly described ORFs are predominantly short, with many encoding proteins of ≤50 amino acids. Codon usage of the newly discovered ORFs suggests that most have not been subject to purifying selection, and hence are unlikely to contribute to cell fitness. Nevertheless, we identify 90 new ORFs (median length of 52 codons) that bear the hallmarks of purifying selection. Thus, our data suggest that pervasive translation of short ORFs in Mycobacterium tuberculosis serves as a rich source for the evolution of new functional proteins.

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.

Interplay between non-coding RNA transcription, stringent phenotype and antibiotic production in Streptomyces

While in recent years the key role of non-coding RNAs (ncRNAs) in regulation of gene expression has become increasingly evident, their interaction with the global regulatory circuits is still obscure. Here we analyzed the structure and organization of the transcriptome of Streptomyces ambofaciens, the producer of spiramycin. We identified ncRNAs including 45 small-RNAs (sRNAs) and 119 antisense-RNAs (asRNAs I) that appear transcribed from dedicated promoters. Some sRNAs and asRNAs are unprecedented in Streptomyces, and were predicted to target mRNAs encoding proteins involved in transcription, translation, ribosomal structure and biogenesis, and regulation of morphological and biochemical differentiation. We then compared ncRNA expression in three strains: i.) the wild type strain; ii.) an isogenic pirA-defective mutant with central carbon metabolism imbalance, "relaxed" phenotype, and repression of antibiotic production; iii.) a pirA-derivative strain harboring a "stringent" RNA polymerase that suppresses pirA-associated phenotypes. Data indicated that expression of most ncRNAs was correlated to the stringent/relaxed phenotype suggesting novel effector mechanisms of the stringent response.

Genome-wide profiling of long non-coding RNA of the rice blast fungus Magnaporthe oryzae during infection

Background

Long non-coding RNAs (lncRNAs) play essential roles in developmental processes and disease development at the transcriptional and post-transcriptional levels across diverse taxa. However, only few studies have profiled fungal lncRNAs in a genome-wide manner during host infection.

Results

Infection-associated lncRNAs were identified using lncRNA profiling over six stages of host infection (e.g., vegetative growth, pre-penetration, biotrophic, and necrotrophic stages) in the model pathogenic fungus, Magnaporthe oryzae. We identified 2,601 novel lncRNAs, including 1,286 antisense lncRNAs and 980 intergenic lncRNAs. Among the identified lncRNAs, 755 were expressed in a stage-specific manner and 560 were infection-specifically expressed lncRNAs (ISELs). To decipher the potential roles of lncRNAs during infection, we identified 365 protein-coding genes that were associated with 214 ISELs. Analysis of the predicted functions of these associated genes suggested that lncRNAs regulate pathogenesis-related genes, including xylanases and effectors.

Conclusions

The ISELs and their associated genes provide a comprehensive view of lncRNAs during fungal pathogen-plant interactions. This study expands new insights into the role of lncRNAs in the rice blast fungus, as well as other plant pathogenic fungi.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08380-4.

Genomic inversions in Escherichia coli alter gene expression and are associated with nucleoid protein binding sites

Genomic reorganization, like rearrangements and inversions, influence how genetic information is organized within bacterial genomes. Inversions in particular, facilitate genome evolution through gene gain and loss, and can alter gene expression. Previous studies investigating the impact inversions have on gene expression induced inversions targeting specific genes or examine inversions between distantly related species. This fails to encompass a genome wide perspective on naturally occurring inversions and their post adaptation impact on gene expression. Here we use bioinformatic techniques and multiple RNA-seq datasets to investigate the short- and long-range impact inversions have on genomic gene expression within <i>Escherichia coli</i>. We observed differences in gene expression between homologous inverted and non-inverted genes, even after long term exposure to adaptive selection. In 4% of inversions representing 33 genes, differential gene expression between inverted and non-inverted homologs was detected, with nearly two thirds (71%) of differentially expressed inverted genes having 9.4-85.6 fold higher gene expression. The identified inversions had more overlap than expected with nucleoid associated protein binding sites, which assist in genomic gene expression regulation. Some inversions can drastically impact gene expression even between different strains of <i>E.coli</i>, and could provide a mechanism for the diversification of genetic content through controlled expression changes.

Regulation of Heterogenous LexA Expression in Staphylococcus aureus by an Antisense RNA Originating from Transcriptional Read-Through upon Natural Mispairings in the sbrB Intrinsic Terminator

Bacterial genomes are pervasively transcribed, generating a wide variety of antisense RNAs (asRNAs). Many of them originate from transcriptional read-through events (TREs) during the transcription termination process. Previous transcriptome analyses revealed that the lexA gene from Staphylococcus aureus, which encodes the main SOS response regulator, is affected by the presence of an asRNA. Here, we show that the lexA antisense RNA (lexA-asRNA) is generated by a TRE on the intrinsic terminator (TT_sbrB) of the sbrB gene, which is located downstream of lexA, in the opposite strand. Transcriptional read-through occurs by a natural mutation that destabilizes the TT_sbrB structure and modifies the efficiency of the intrinsic terminator. Restoring the mispairing mutation in the hairpin of TT_sbrB prevented lexA-asRNA transcription. The level of lexA-asRNA directly correlated with cellular stress since the expressions of sbrB and lexA-asRNA depend on the stress transcription factor SigB. Comparative analyses revealed strain-specific nucleotide polymorphisms within TT_sbrB, suggesting that this TT could be prone to accumulating natural mutations. A genome-wide analysis of TREs suggested that mispairings in TT hairpins might provide wider transcriptional connections with downstream genes and, ultimately, transcriptomic variability among S. aureus strains.

Elongating RNA polymerase II and RNA:DNA hybrids hinder fork progression and gene expression at sites of head-on replication-transcription collisions

Uncoordinated clashes between replication forks and transcription cause replication stress and genome instability, which are hallmarks of cancer and neurodegeneration. Here, we investigate the outcomes of head-on replication-transcription collisions, using as a model system budding yeast mutants for the helicase Sen1, the ortholog of human Senataxin. We found that RNA Polymerase II accumulates together with RNA:DNA hybrids at sites of head-on collisions. The replication fork and RNA Polymerase II are both arrested during the clash, leading to DNA damage and, in the long run, the inhibition of gene expression. The inactivation of RNA Polymerase II elongation factors, such as the HMG-like protein Spt2 and the DISF and PAF complexes, but not alterations in chromatin structure, allows replication fork progression through transcribed regions. Attenuation of RNA Polymerase II elongation rescues RNA:DNA hybrid accumulation and DNA damage sensitivity caused by the absence of Sen1, but not of RNase H proteins, suggesting that such enzymes counteract toxic RNA:DNA hybrids at different stages of the cell cycle with Sen1 mainly acting in replication. We suggest that the main obstacle to replication fork progression is the elongating RNA Polymerase II engaged in an R-loop, rather than RNA:DNA hybrids per se or hybrid-associated chromatin modifications.

Riboregulation in bacteria: From general principles to novel mechanisms of the trp attenuator and its sRNA and peptide products

Gene expression strategies ensuring bacterial survival and competitiveness rely on cis- and trans-acting RNA-regulators (riboregulators). Among the cis-acting riboregulators are transcriptional and translational attenuators, and antisense RNAs (asRNAs). The trans-acting riboregulators are small RNAs (sRNAs) that bind proteins or base pairs with other RNAs. This classification is artificial since some regulatory RNAs act both in cis and in trans, or function in addition as small mRNAs. A prominent example is the archetypical, ribosome-dependent attenuator of tryptophan (Trp) biosynthesis genes. It responds by transcription attenuation to two signals, Trp availability and inhibition of translation, and gives rise to two trans-acting products, the attenuator sRNA rnTrpL and the leader peptide peTrpL. In Escherichia coli, rnTrpL links Trp availability to initiation of chromosome replication and in Sinorhizobium meliloti, it coordinates regulation of split tryptophan biosynthesis operons. Furthermore, in S. meliloti, peTrpL is involved in mRNA destabilization in response to antibiotic exposure. It forms two types of asRNA-containing, antibiotic-dependent ribonucleoprotein complexes (ARNPs), one of them changing the target specificity of rnTrpL. The posttranscriptional role of peTrpL indicates two emerging paradigms: (1) sRNA reprograming by small molecules and (2) direct involvement of antibiotics in regulatory RNPs. They broaden our view on RNA-based mechanisms and may inspire new approaches for studying, detecting, and using antibacterial compounds. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Small Molecule-RNA Interactions RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs.

A vast pool of lineage-specific microproteins encoded by long non-coding RNAs in plants

Pervasive transcription of eukaryotic genomes results in expression of long non-coding RNAs (lncRNAs) most of which are poorly conserved in evolution and appear to be non-functional. However, some lncRNAs have been shown to perform specific functions, in particular, transcription regulation. Thousands of small open reading frames (smORFs, <100 codons) located on lncRNAs potentially might be translated into peptides or microproteins. We report a comprehensive analysis of the conservation and evolutionary trajectories of lncRNAs-smORFs from the moss Physcomitrium patens across transcriptomes of 479 plant species. Although thousands of smORFs are subject to substantial purifying selection, the majority of the smORFs appear to be evolutionary young and could represent a major pool for functional innovation. Using nanopore RNA sequencing, we show that, on average, the transcriptional level of conserved smORFs is higher than that of non-conserved smORFs. Proteomic analysis confirmed translation of 82 novel species-specific smORFs. Numerous conserved smORFs containing low complexity regions (LCRs) or transmembrane domains were identified, the biological functions of a selected LCR-smORF were demonstrated experimentally. Thus, microproteins encoded by smORFs are a major, functionally diverse component of the plant proteome.

Non-Coding RNAs: Emerging Therapeutic Targets in Spinal Cord Ischemia-Reperfusion Injury

Paralysis or paraplegia caused by transient or permanent spinal cord ischemia–reperfusion injury (SCIRI) remains one of the most devastating post-operative complications after thoracoabdominal aortic surgery, even though perioperative strategies and surgical techniques continue to improve. Uncovering the molecular and cellular pathophysiological processes in SCIRI has become a top priority. Recently, the expression, function, and mechanism of non-coding RNAs (ncRNAs) in various diseases have drawn wide attention. Non-coding RNAs contain a variety of biological functions but do not code for proteins. Previous studies have shown that ncRNAs play a critical role in SCIRI. However, the character of ncRNAs in attenuating SCIRI has not been systematically summarized. This review article will be the first time to assemble the knowledge of ncRNAs regulating apoptosis, inflammation, autophagy, and oxidative stress to attenuate SCIRI. A better understanding of the functional significance of ncRNAs following SCIRI could help us to identify novel therapeutic targets and develop potential therapeutic strategies. All the current research about the function of nRNAs in SCIRI will be summarized one by one in this review.

Dynamics of the compartmentalized Streptomyces chromosome during metabolic differentiation

Bacteria of the genus Streptomyces are prolific producers of specialized metabolites, including antibiotics. The linear chromosome includes a central region harboring core genes, as well as extremities enriched in specialized metabolite biosynthetic gene clusters. Here, we show that chromosome structure in Streptomyces ambofaciens correlates with genetic compartmentalization during exponential phase. Conserved, large and highly transcribed genes form boundaries that segment the central part of the chromosome into domains, whereas the terminal ends tend to be transcriptionally quiescent compartments with different structural features. The onset of metabolic differentiation is accompanied by a rearrangement of chromosome architecture, from a rather 'open' to a 'closed' conformation, in which highly expressed specialized metabolite biosynthetic genes form new boundaries. Thus, our results indicate that the linear chromosome of S. ambofaciens is partitioned into structurally distinct entities, suggesting a link between chromosome folding, gene expression and genome evolution.

The emerging role of circular RNAs in spinal cord injury

Spinal cord injury (SCI) is one kind of severe diseases with high mortality and morbidity worldwide, and lacks effective therapeutic interventions currently, which leads to not only permanent neurological impairments but also heavy social and economic burden. Recent studies have proved that circRNAs are highly expressed in neural tissues, regulating the neuronal and synaptic functions. What's more, significantly altered circRNAs expression profiles are closely associated with the pathophysiology of SCI. In this review, we summarize the current advance on the role of circRNAs in SCI, which may provide a better understanding of pathogenesis and therapeutic strategies of SCI.

THE TRANSLATIONAL POTENTIAL OF THIS ARTICLE

The Translational potential of this article is that A further understanding of circRNAs in the pathogenesis of SCI will promote the circRNA-based clinical applications.

Burning the Candle at Both Ends: Have Exoribonucleases Driven Divergence of Regulatory RNA Mechanisms in Bacteria?

Regulatory RNAs have emerged as ubiquitous gene regulators in all bacterial species studied to date. The combination of sequence-specific RNA interactions and malleable RNA structure has allowed regulatory RNA to adopt different mechanisms of gene regulation in a diversity of genetic backgrounds. In the model GammaproteobacteriaEscherichia coli and Salmonella, the regulatory RNA chaperone Hfq appears to play a global role in gene regulation, directly controlling ∼20 to 25% of the entire transcriptome. While the model FirmicutesBacillus subtilis and Staphylococcus aureus encode a Hfq homologue, its role has been significantly depreciated. These bacteria also have marked differences in RNA turnover. E. coli and Salmonella degrade RNA through internal endonucleolytic and 3′→5′ exonucleolytic cleavage that appears to allow transient accumulation of mRNA 3′ UTR cleavage fragments that contain stabilizing 3′ structures. In contrast, B. subtilis and S. aureus are able to exonucleolytically attack internally cleaved RNA from both the 5′ and 3′ ends, efficiently degrading mRNA 3′ UTR fragments. Here, we propose that the lack of 5′→3′ exoribonuclease activity in Gammaproteobacteria has allowed the accumulation of mRNA 3′ UTR ends as the “default” setting. This in turn may have provided a larger pool of unconstrained RNA sequences that has fueled the expansion of Hfq function and small RNA (sRNA) regulation in E. coli and Salmonella. Conversely, the exoribonuclease RNase J may be a significant barrier to the evolution of 3′ UTR sRNAs in B. subtilis and S. aureus that has limited the pool of RNA ligands available to Hfq and other sRNA chaperones, depreciating their function in these model Firmicutes.

Interplay between Non-Coding RNA Transcription, Stringent/Relaxed Phenotype and Antibiotic Production in Streptomyces ambofaciens

While in recent years the key role of non-coding RNAs (ncRNAs) in the regulation of gene expression has become increasingly evident, their interaction with the global regulatory circuits is still obscure. Here we analyzed the structure and organization of the transcriptome of Streptomyces ambofaciens, the producer of spiramycin. We identified ncRNAs including 45 small-RNAs (sRNAs) and 119 antisense-RNAs (asRNAs I) that appear transcribed from dedicated promoters. Some sRNAs and asRNAs are unprecedented in Streptomyces and were predicted to target mRNAs encoding proteins involved in transcription, translation, ribosomal structure and biogenesis, and regulation of morphological and biochemical differentiation. We then compared ncRNA expression in three strains: (i) the wild-type strain; (ii) an isogenic pirA-defective mutant with central carbon metabolism imbalance, “relaxed” phenotype, and repression of antibiotic production; and (iii) a pirA-derivative strain harboring a “stringent” RNA polymerase that suppresses pirA-associated phenotypes. Data indicated that the expression of most ncRNAs was correlated to the stringent/relaxed phenotype suggesting novel effector mechanisms of the stringent response.

Pervasive transcription of the mitochondrial genome in Candida albicans is revealed in mutants lacking the mtEXO RNase complex

The mitochondrial genome of the pathogenic yeast Candida albicans displays a typical organization of several (eight) primary transcription units separated by noncoding regions. Presence of genes encoding Complex I subunits and the stability of its mtDNA sequence make it an attractive model to study organellar genome expression using transcriptomic approaches. The main activity responsible for RNA degradation in mitochondria is a two-component complex (mtEXO) consisting of a 3ʹ-5ʹ exoribonuclease, in yeasts encoded by the DSS1 gene, and a conserved Suv3p helicase. In C. albicans, deletion of either DSS1 or SUV3 gene results in multiple defects in mitochondrial genome expression leading to the loss of respiratory competence. Transcriptomic analysis reveals pervasive transcription in mutants lacking the mtEXO activity, with evidence of the entire genome being transcribed, whereas in wild-type strains no RNAs corresponding to a significant fraction of the noncoding genome can be detected. Antisense (‘mirror’) transcripts, absent from normal mitochondria are also prominent in the mutants. The expression of multiple mature transcripts, particularly those translated from bicistronic mRNAs, as well as those that contain introns is affected in the mutants, resulting in a decreased level of proteins and reduced respiratory complex activity. The phenotype is most severe in the case of Complex IV, where a decrease of mature COX1 mRNA level to ~5% results in a complete loss of activity. These results show that RNA degradation by mtEXO is essential for shaping the mitochondrial transcriptome and is required to maintain the functional demarcation between transcription units and non-coding genome segments.

An RNA-centric global view of Clostridioides difficile reveals broad activity of Hfq in a clinically important gram-positive bacterium

The gram-positive human pathogen Clostridioides difficile has emerged as the leading cause of antibiotic-associated diarrhea. However, little is known about the bacterium's transcriptome architecture and mechanisms of posttranscriptional control. Here, we have applied transcription start site and termination mapping to generate a single-nucleotide-resolution RNA map of C. difficile 5' and 3' untranslated regions, operon structures, and noncoding regulators, including 42 sRNAs. Our results indicate functionality of many conserved riboswitches and predict cis-regulatory RNA elements upstream of multidrug resistance (MDR)-type ATP-binding cassette (ABC) transporters and transcriptional regulators. Despite growing evidence for a role of Hfq in RNA-based gene regulation in C. difficile, the functions of Hfq-based posttranscriptional regulatory networks in gram-positive pathogens remain controversial. Using Hfq immunoprecipitation followed by sequencing of bound RNA species (RIP-seq), we identify a large cohort of transcripts bound by Hfq and show that absence of Hfq affects transcript stabilities and steady-state levels. We demonstrate sRNA expression during intestinal colonization by C. difficile and identify infection-related signals impacting its expression. As a proof of concept, we show that the utilization of the abundant intestinal metabolite ethanolamine is regulated by the Hfq-dependent sRNA CDIF630nc_085. Overall, our study lays the foundation for understanding clostridial riboregulation with implications for the infection process and provides evidence for a global role of Hfq in posttranscriptional regulation in a gram-positive bacterium.

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.

Best practices on the differential expression analysis of multi-species RNA-seq

Advances in transcriptome sequencing allow for simultaneous interrogation of differentially expressed genes from multiple species originating from a single RNA sample, termed dual or multi-species transcriptomics. Compared to single-species differential expression analysis, the design of multi-species differential expression experiments must account for the relative abundances of each organism of interest within the sample, often requiring enrichment methods and yielding differences in total read counts across samples. The analysis of multi-species transcriptomics datasets requires modifications to the alignment, quantification, and downstream analysis steps compared to the single-species analysis pipelines. We describe best practices for multi-species transcriptomics and differential gene expression.

Engineering Transcriptional Interference through RNA Polymerase Processivity Control

Antisense transcription is widespread in all kingdoms of life and has been shown to influence gene expression through transcriptional interference (TI), a phenomenon in which one transcriptional process negatively influences another in cis. The processivity, or uninterrupted transcription, of an RNA polymerase (RNAP) is closely tied to levels of antisense transcription in bacterial genomes, but its influence on TI, while likely important, is not well-characterized. Here, we show that TI can be tuned through processivity control via three distinct antitermination strategies: the antibiotic bicyclomycin, phage protein Psu, and ribosome-RNAP coupling. We apply these methods toward TI and tune ribosome-RNAP coupling to produce 38-fold transcription-level gene repression due to both RNAP collisions and antisense RNA interference. We then couple protein roadblock and TI to design minimal genetic NAND and NOR logic gates. Together, these results show the importance of processivity control for strong TI and demonstrate TI's potential for synthetic biology.

Transcription initiation mapping in 31 bovine tissues reveals complex promoter activity, pervasive transcription, and tissue-specific promoter usage

Characterizing transcription start sites is essential for understanding the regulatory mechanisms that control gene expression. Recently, a new bovine genome assembly (ARS-UCD1.2) with high continuity, accuracy, and completeness was released; however, the functional annotation of the bovine genome lacks precise transcription start sites and contains a low number of transcripts in comparison to human and mouse. By using the RAMPAGE approach, this study identified transcription start sites at high resolution in a large collection of bovine tissues. We found several known and novel transcription start sites attributed to promoters of protein-coding and lncRNA genes that were validated through experimental and in silico evidence. With these findings, the annotation of transcription start sites in cattle reached a level comparable to the mouse and human genome annotations. In addition, we identified and characterized transcription start sites for antisense transcripts derived from bidirectional promoters, potential lncRNAs, mRNAs, and pre-miRNAs. We also analyzed the quantitative aspects of RAMPAGE to produce a promoter activity atlas, reaching highly reproducible results comparable to traditional RNA-seq. Coexpression networks revealed considerable use of tissue-specific promoters, especially between brain and testicle, which expressed several genes in common from alternate loci. Furthermore, regions surrounding coexpressed modules were enriched in binding factor motifs representative of each tissue. The comprehensive annotation of promoters in such a large collection of tissues will substantially contribute to our understanding of gene expression in cattle and other mammalian species, shortening the gap between genotypes and phenotypes.