The Non-Coding Regulatory RNA Revolution in Archaea

TargetRNA3: predicting prokaryotic RNA regulatory targets with machine learning

Small regulatory RNAs pervade prokaryotes, with the best-studied family of these non-coding genes corresponding to trans-acting regulators that bind via base pairing to their message targets. Given the increasing frequency with which these genes are being identified, it is important that methods for illuminating their regulatory targets keep pace. Using a machine learning approach, we investigate thousands of interactions between small RNAs and their targets, and we interrogate more than a hundred features indicative of these interactions. We present a new method, TargetRNA3, for predicting targets of small RNA regulators and show that it outperforms existing approaches. TargetRNA3 is available at https://cs.wellesley.edu/~btjaden/TargetRNA3 .

Little reason to call them small noncoding RNAs

Hundreds of different species of small RNAs can populate a bacterial cell. This small transcriptome contains important information for the adaptation of cellular physiology to environmental changes. Underlying cellular networks involving small RNAs are RNA-RNA and RNA-protein interactions, which are often intertwined. In addition, small RNAs can function as mRNAs. In general, small RNAs are referred to as noncoding because very few are known to contain translated open reading frames. In this article, we intend to highlight that the number of small RNAs that fall within the set of translated RNAs is bound to increase. In addition, we aim to emphasize that the dynamics of the small transcriptome involve different functional codes, not just the genetic code. Therefore, since the role of small RNAs is always code-driven, we believe that there is little reason to continue calling them small noncoding RNAs.

Functional diversity of nanohaloarchaea within xylan-degrading consortia

Extremely halophilic representatives of the phylum Candidatus Nanohaloarchaeota (members of the DPANN superphyla) are obligately associated with extremely halophilic archaea of the phylum Halobacteriota (according to the GTDB taxonomy). Using culture-independent molecular techniques, their presence in various hypersaline ecosystems around the world has been confirmed over the past decade. However, the vast majority of nanohaloarchaea remain uncultivated, and thus their metabolic capabilities and ecophysiology are currently poorly understood. Using the (meta)genomic, transcriptomic, and DNA methylome platforms, the metabolism and functional prediction of the ecophysiology of two novel extremely halophilic symbiotic nanohaloarchaea (Ca. Nanohalococcus occultus and Ca. Nanohalovita haloferacivicina) stably cultivated in the laboratory as members of a xylose-degrading binary culture with a haloarchaeal host, Haloferax lucentense, was determined. Like all known DPANN superphylum nanoorganisms, these new sugar-fermenting nanohaloarchaea lack many fundamental biosynthetic repertoires, making them exclusively dependent on their respective host for survival. In addition, given the cultivability of the new nanohaloarchaea, we managed to discover many unique features in these new organisms that have never been observed in nano-sized archaea both within the phylum Ca. Nanohaloarchaeota and the entire superphylum DPANN. This includes the analysis of the expression of organism-specific non-coding regulatory (nc)RNAs (with an elucidation of their 2D-secondary structures) as well as profiling of DNA methylation. While some ncRNA molecules have been predicted with high confidence as RNAs of an archaeal signal recognition particle involved in delaying protein translation, others resemble the structure of ribosome-associated ncRNAs, although none belong to any known family. Moreover, the new nanohaloarchaea have very complex cellular defense mechanisms. In addition to the defense mechanism provided by the type II restriction-modification system, consisting of Dcm-like DNA methyltransferase and Mrr restriction endonuclease, Ca. Nanohalococcus encodes an active type I-D CRISPR/Cas system, containing 77 spacers divided into two loci. Despite their diminutive genomes and as part of their host interaction mechanism, the genomes of new nanohaloarchaea do encode giant surface proteins, and one of them (9,409 amino acids long) is the largest protein of any sequenced nanohaloarchaea and the largest protein ever discovered in cultivated archaea.

G. de Lomana A, Gomes-Filho JV, Vêncio RZN, Moritz RL, Koide T, Baliga NS. A Genome-Scale Atlas Reveals Complex Interplay of Transcription and Translation in an Archaeon

The scale of post-transcriptional regulation and the implications of its interplay with other forms of regulation in environmental acclimation are underexplored for organisms of the domain Archaea. Here, we have investigated the scale of post-transcriptional regulation in the extremely halophilic archaeon Halobacterium salinarum NRC-1 by integrating the transcriptome-wide locations of transcript processing sites (TPSs) and SmAP1 binding, the genome-wide locations of antisense RNAs (asRNAs), and the consequences of RNase_2099C knockout on the differential expression of all genes. This integrated analysis has discovered that 54% of all protein-coding genes in the genome of this haloarchaeon are likely targeted by multiple mechanisms for putative post-transcriptional processing and regulation, with about 20% of genes likely being regulated by combinatorial schemes involving SmAP1, asRNAs, and RNase_2099C. Comparative analysis of mRNA levels (transcriptome sequencing [RNA-Seq]) and protein levels (sequential window acquisition of all theoretical fragment ion spectra mass spectrometry [SWATH-MS]) for 2,579 genes over four phases of batch culture growth in complex medium generated additional evidence for the conditional post-transcriptional regulation of 7% of all protein-coding genes. We demonstrate that post-transcriptional regulation may act to fine-tune specialized and rapid acclimation to stressful environments, e.g., as a switch to turn on gas vesicle biogenesis to promote vertical relocation under anoxic conditions and modulate the frequency of transposition by insertion sequence (IS) elements of the IS200/IS605, IS4, and ISH3 families. Findings from this study are provided as an atlas in a public Web resource (https://halodata.systemsbiology.net). IMPORTANCE While the transcriptional regulation landscape of archaea has been extensively investigated, we currently have limited knowledge about post-transcriptional regulation and its driving mechanisms in this domain of life. In this study, we collected and integrated omics data from multiple sources and technologies to infer post-transcriptionally regulated genes and the putative mechanisms modulating their expression at the protein level in Halobacterium salinarum NRC-1. The results suggest that post-transcriptional regulation may drive environmental acclimation by regulating hallmark biological processes. To foster discoveries by other research groups interested in the topic, we extended our integrated data to the public in the form of an interactive atlas (https://halodata.systemsbiology.net).

Cloning and Sequencing Eukaryotic Small RNAs

Small RNAs are ubiquitous regulators of gene expression that participate in nearly all aspects of physiology in a wide range of organisms. There are many different classes of eukaryotic small RNAs that play regulatory roles at every level of gene expression, including transcription, RNA stability, and translation. While eukaryotic small RNAs display diverse functions across and within classes, they are generally grouped functionally based on the machinery required for their biogenesis, the effector proteins they associate with, and their molecular characteristics. The development of techniques to clone and sequence small RNAs has been critical for their identification, yet the ligation-dependent addition of RNA adapters and the use of reverse transcriptase to generate cDNA in traditional library preparation protocols can be unsuitable to detect certain small RNA subtypes. In particular, 3' or 5' chemical modifications that are characteristic of specific types of small RNAs can impede the ligation-dependent addition of RNA adapters, while internal RNA modifications can interfere with accurate reverse transcription. The inability to clone certain small RNA subtypes with traditional protocols results in an inaccurate assessment of small RNA abundance and diversity, where some RNAs appear over-represented and others are not detected. This overview aims to guide users on how to design small RNA cloning workflows in eukaryotes to more accurately capture specific small RNAs of interest. Hence, we discuss the molecular biology underlying the identification and quantitation of small RNAs, explore the limitations of commonly used protocols, and detail the alternative approaches that can be used to enrich specific small RNA classes. © 2022 Wiley Periodicals LLC.

The origin of RNA interference: Adaptive or neutral evolution?

The origin of RNA interference (RNAi) is usually explained by a defense-based hypothesis, in which RNAi evolved as a defense against transposable elements (TEs) and RNA viruses and was already present in the last eukaryotic common ancestor (LECA). However, since RNA antisense regulation and double-stranded RNAs (dsRNAs) are ancient and widespread phenomena, the origin of defensive RNAi should have occurred in parallel with its regulative functions to avoid imbalances in gene regulation. Thus, we propose a neutral evolutionary hypothesis for the origin of RNAi in which qualitative system drift from a prokaryotic antisense RNA gene regulation mechanism leads to the formation of RNAi through constructive neutral evolution (CNE). We argue that RNAi was already present in the ancestor of LECA before the need for a new defense system arose and that its presence helped to shape eukaryotic genomic architecture and stability.

Where does RNA interference come from? This Essay describes a new step-by-step evolutionary model of how RNA interference might have originated in early eukaryotes through neutral events from the molecular machinery present in prokaryotes.

Small RNA-Sequencing Library Preparation for the Halophilic Archaeon Haloferax volcanii

Posttranscriptional regulation actuated by small RNAs (sRNAs) plays essential roles in a wide variety of cellular processes, especially in stress responses and environmental signaling. Hundreds of sRNAs have recently been discovered in archaea using genome-wide approaches but the molecular mechanisms of only a few have been characterized experimentally. Here, we describe how to build sRNA sequencing libraries using size-selected total RNA in the model archaeon, Haloferax volcanii , to provide a tool to further characterize sRNAs in archaea.

Post-transcriptional regulation of redox homeostasis by the small RNA SHOxi in haloarchaea

While haloarchaea are highly resistant to oxidative stress, a comprehensive understanding of the processes regulating this remarkable response is lacking. Oxidative stress-responsive small non-coding RNAs (sRNAs) have been reported in the model archaeon, Haloferax volc anii, but targets and mechanisms have not been elucidated. Using a combination of high throughput and reverse molecular genetic approaches, we elucidated the functional role of the most up-regulated intergenic sRNA during oxidative stress in H. volcanii, named Small RNA in Haloferax Oxidative Stress (SHOxi). SHOxi was predicted to form a stable secondary structure with a conserved stem-loop region as the potential binding site for trans-targets. NAD-dependent malic enzyme mRNA, identified as a putative target of SHOxi, interacted directly with a putative 'seed' region within the predicted stem loop of SHOxi. Malic enzyme catalyzes the oxidative decarboxylation of malate into pyruvate using NAD+ as a cofactor. The destabilization of malic enzyme mRNA, and the decrease in the NAD+/NADH ratio, resulting from the direct RNA-RNA interaction between SHOxi and its trans-target was essential for the survival of H. volcanii to oxidative stress. These findings indicate that SHOxi likely regulates redox homoeostasis during oxidative stress by the post-transcriptional destabilization of malic enzyme mRNA. SHOxi-mediated regulation provides evidence that the fine-tuning of metabolic cofactors could be a core strategy to mitigate damage from oxidative stress and confer resistance. This study is the first to establish the regulatory effects of sRNAs on mRNAs during the oxidative stress response in Archaea.

Comparative genomics of DNA-binding transcription factors in archaeal and bacterial organisms

Archaea represent a diverse phylogenetic group that includes free-living, extremophile, mesophile, symbiont, and opportunistic organisms. These prokaryotic organisms share a high significant similarity with the basal transcriptional machinery of Eukarya, and they share regulatory mechanisms with Bacteria, such as operonic organization and DNA-binding transcription factors (TFs). In this work, we identified the repertoire of TFs in 415 archaeal genomes and compared them with their counterparts in bacterial genomes. The comparisons of TFs, at a global level and per family, allowed us to identify similarities and differences between the repertoires of regulatory proteins of bacteria and archaea. For example, 11 of 62 families are more highly abundant in archaea than bacteria, and 13 families are abundant in bacteria but not in archaea and 38 families have similar abundances in the two groups. In addition, we found that archaeal TFs have a lower isoelectric point than bacterial proteins, i.e., they contain more acidic amino acids, and are smaller than bacterial TFs. Our findings suggest a divergence occurred for the regulatory proteins, even though they are common to archaea and bacteria. We consider that this analysis contributes to the comprehension of the structure and functionality of regulatory proteins of archaeal organisms.

Halobacterium salinarum and Haloferax volcanii Comparative Transcriptomics Reveals Conserved Transcriptional Processing Sites

Post-transcriptional processing of messenger RNA is an important regulatory strategy that allows relatively fast responses to changes in environmental conditions. In halophile systems biology, the protein perspective of this problem (i.e., ribonucleases which implement the cleavages) is generally more studied than the RNA perspective (i.e., processing sites). In the present in silico work, we mapped genome-wide transcriptional processing sites (TPS) in two halophilic model organisms, Halobacterium salinarum NRC-1 and Haloferax volcanii DS2. TPS were established by reanalysis of publicly available differential RNA-seq (dRNA-seq) data, searching for non-primary (monophosphorylated RNAs) enrichment. We found 2093 TPS in 43% of H. salinarum genes and 3515 TPS in 49% of H. volcanii chromosomal genes. Of the 244 conserved TPS sites found, the majority were located around start and stop codons of orthologous genes. Specific genes are highlighted when discussing antisense, ribosome and insertion sequence associated TPS. Examples include the cell division gene ftsZ2, whose differential processing signal along growth was detected and correlated with post-transcriptional regulation, and biogenesis of sense overlapping transcripts associated with IS200/IS605. We hereby present the comparative, transcriptomics-based processing site maps with a companion browsing interface.

A Small RNA Is Linking CRISPR-Cas and Zinc Transport

The function and mode of action of small regulatory RNAs is currently still understudied in archaea. In the halophilic archaeon Haloferax volcanii, a plethora of sRNAs have been identified; however, in-depth functional analysis is missing for most of them. We selected a small RNA (s479) from Haloferax volcanii for detailed characterization. The sRNA gene is encoded between a CRISPR RNA locus and the Cas protein gene cluster, and the s479 deletion strain is viable and was characterized in detail. Transcriptome studies of wild-type Haloferax cells and the deletion mutant revealed upregulation of six genes in the deletion strain, showing that this sRNA has a clearly defined function. Three of the six upregulated genes encode potential zinc transporter proteins (ZnuA1, ZnuB1, and ZnuC1) suggesting the involvement of s479 in the regulation of zinc transport. Upregulation of these genes in the deletion strain was confirmed by northern blot and proteome analyses. Furthermore, electrophoretic mobility shift assays demonstrate a direct interaction of s479 with the target znuC1 mRNA. Proteome comparison of wild-type and deletion strains further expanded the regulon of s479 deeply rooting this sRNA within the metabolism of H. volcanii especially the regulation of transporter abundance. Interestingly, s479 is not only encoded next to CRISPR-cas genes, but the mature s479 contains a crRNA-like 5' handle, and experiments with Cas protein deletion strains indicate maturation by Cas6 and interaction with Cas proteins. Together, this might suggest that the CRISPR-Cas system is involved in s479 function.

Coupling of Transcription and Translation in Archaea: Cues From the Bacterial World

The lack of a nucleus is the defining cellular feature of bacteria and archaea. Consequently, transcription and translation are occurring in the same compartment, proceed simultaneously and likely in a coupled fashion. Recent cryo-electron microscopy (cryo-EM) and tomography data, also combined with crosslinking-mass spectrometry experiments, have uncovered detailed structural features of the coupling between a transcribing bacterial RNA polymerase (RNAP) and the trailing translating ribosome in Escherichia coli and Mycoplasma pneumoniae. Formation of this supercomplex, called expressome, is mediated by physical interactions between the RNAP-bound transcription elongation factors NusG and/or NusA and the ribosomal proteins including uS10. Based on the structural conservation of the RNAP core enzyme, the ribosome, and the universally conserved elongation factors Spt5 (NusG) and NusA, we discuss requirements and functional implications of transcription-translation coupling in archaea. We furthermore consider additional RNA-mediated and co-transcriptional processes that potentially influence expressome formation in archaea.

RNAmining: A machine learning stand-alone and web server tool for RNA coding potential prediction

Non-coding RNAs (ncRNAs) are important players in the cellular regulation of organisms from different kingdoms. One of the key steps in ncRNAs research is the ability to distinguish coding/non-coding sequences. We applied seven machine learning algorithms (Naive Bayes, Support Vector Machine, K-Nearest Neighbors, Random Forest, Extreme Gradient Boosting, Neural Networks and Deep Learning) through model organisms from different evolutionary branches to create a stand-alone and web server tool (RNAmining) to distinguish coding and non-coding sequences. Firstly, we used coding/non-coding sequences downloaded from Ensembl (April 14th, 2020). Then, coding/non-coding sequences were balanced, had their trinucleotides count analysed (64 features) and we performed a normalization by the sequence length, resulting in total of 180 models. The machine learning algorithms validations were performed using 10-fold cross-validation and we selected the algorithm with the best results (eXtreme Gradient Boosting) to implement at RNAmining. Best F1-scores ranged from 97.56% to 99.57% depending on the organism. Moreover, we produced a benchmarking with other tools already in literature (CPAT, CPC2, RNAcon and TransDecoder) and our results outperformed them. Both stand-alone and web server versions of RNAmining are freely available at https://rnamining.integrativebioinformatics.me/.

Role of Non-coding RNAs in Fungal Pathogenesis and Antifungal Drug Responses

Purpose of Review

Non-coding RNAs (ncRNAs), including regulatory small RNAs (sRNAs) and long non-coding RNAs (lncRNAs), constitute a significant part of eukaryotic genomes; however, their roles in fungi are just starting to emerge. ncRNAs have been shown to regulate gene expression in response to varying environmental conditions (like stress) and response to chemicals, including antifungal drugs. In this review, I highlighted recent studies focusing on the functional roles of ncRNAs in pathogenic fungi.

Recent Findings

Emerging evidence suggests sRNAs (small RNAs) and lncRNAs (long non-coding RNAs) play an important role in fungal pathogenesis and antifungal drug response. Their roles include posttranscriptional gene silencing, histone modification, and chromatin remodeling. Fungal pathogens utilize RNA interference (RNAi) mechanisms to regulate pathogenesis-related genes and can also transfer sRNAs inside the host to suppress host immunity genes to increase virulence. Hosts can also transfer sRNAs to induce RNAi in fungal pathogens to reduce virulence. Additionally, sRNAs and lncRNAs also regulate gene expression in response to antifungal drugs increasing resistance (and possibly tolerance) to drugs.

Summary

Herein, I discuss what is known about ncRNAs in fungal pathogenesis and antifungal drug responses. Advancements in genomic technologies will help identify the ncRNA repertoire in fungal pathogens, and functional studies will elucidate their mechanisms. This will advance our understanding of host-fungal interactions and potentially help develop better treatment strategies.

Screening and verification of long noncoding RNA promoter methylation sites in hepatocellular carcinoma

Background

Long noncoding ribonucleic acid (lncRNA) promoter methylation is closely related to the occurrence and development of hepatocellular carcinoma (HCC). Thus, we aim to screen and verify the lncRNA promoter methylation sites associated with overall survival (OS), vascular invasion, pathological grade, and clinical stage in HCC.

Methods

Methylation-related data including clinical characteristic, transcriptome, methylation, and messenger RNA (mRNA) expression were taken from the Cancer Genome Atlas (TCGA) database. The OS, vascular invasion, pathological grade, and clinical stage-related lncRNA promoter methylation models were developed by the least absolute shrinkage and selection operator (LASSO) algorithm based on the lncRNA promoter methylation sites screened via R software. The Kaplan-Meier analysis, the area under the receiver operating characteristic (ROC) curve (AUC), the calibration curve (C-index) were performed to evaluate the performance of these models. Finally, the methylation-specific polymerase chain reaction (MS-PCR) was performed to verify the accuracy of these models based on 146 HCC tissues from our hospital.

Results

A total of 10 methylation sites were included in the OS-related lncRNA promoter methylation model that could effectively divide HCC patients into high-risk and low-risk groups (P < 0.0001) via survival analysis. COX univariable and multivariable regression analysis found that the OS-related model (P < 0.001, 95% CI 1.378-2.942) and T stage (P < 0.001, 95% CI 1.490-3.418) were independent risk factors affecting OS in HCC patients. The vascular invasion-related model contained 8 methylation sites with its AUC value of 0.657; the pathological grade-related model contained 22 methylation sites with its AUC value of 0.797; the clinical stage-related model contained 13 methylation sites with its AUC of 0.724. Target genes corresponded to vascular invasion-related lncRNA promoter methylation sites were involved in many kinds of biological processes in HCC such as PI3K-Akt signaling pathway. The accuracy of the vascular invasion-related model was consistent with our bioinformatics conclusion after being verified via MS-PCR.

Conclusion

The lncRNA promoter methylation sites are closely correlated with the process of HCC and can be utilized to improve the therapy and prognosis of HCC.

RNA-Binding Proteins Driving the Regulatory Activity of Small Non-coding RNAs in Bacteria

Small non-coding RNAs (sRNAs) are critical post-transcriptional regulators of gene expression. Distinct RNA-binding proteins (RBPs) influence the processing, stability and activity of bacterial small RNAs. The vast majority of bacterial sRNAs interact with mRNA targets, affecting mRNA stability and/or its translation rate. The assistance of RNA-binding proteins facilitates and brings accuracy to sRNA-mRNA basepairing and the RNA chaperones Hfq and ProQ are now recognized as the most prominent RNA matchmakers in bacteria. These RBPs exhibit distinct high affinity RNA-binding surfaces, promoting RNA strand interaction between a trans-encoding sRNA and its mRNA target. Nevertheless, some organisms lack ProQ and/or Hfq homologs, suggesting the existence of other RBPs involved in sRNA function. Along this line of thought, the global regulator CsrA was recently shown to facilitate the access of an sRNA to its target mRNA and may represent an additional factor involved in sRNA function. Ribonucleases (RNases) can be considered a class of RNA-binding proteins with nucleolytic activity that are responsible for RNA maturation and/or degradation. Presently RNase E, RNase III, and PNPase appear to be the main players not only in sRNA turnover but also in sRNA processing. Here we review the current knowledge on the most important bacterial RNA-binding proteins affecting sRNA activity and sRNA-mediated networks.

The rare lncRNA GOLLD is widespread and structurally conserved among Mycobacterium tRNA arrays

Noncoding RNA (ncRNA) genes produce transcripts involved in a wide range of functions, including catalytic and regulatory functions. Besides, some transcripts have highly complex structures that may impact their activities. Among the largest bacterial ncRNAs, there is the rare GOLLD RNA, which is associated with tRNA genes and supposed to be chromosome- and phage-encoded in specialized groups of bacteria, including those from Lactobacillales and Actinomycetales orders. The only GOLLD structure was inferred from a variety of sequences, including many marine metagenomes. To explore GOLLD RNA in bacterial genomes, we mined the GOLLD gene in thousands of Mycobacterium and virus genomes using Infernal software. We identified this gene in 350 mycobacteria, including megaplasmids, and 39 bacteriophages, mainly in the genomic context of tRNA arrays. Mycobacterium GOLLD genes presented a high diversity and were distributed in three phylogenetic groups: (i) Mycobacterium exclusive; (ii) Mycobacterium and mycobacteriophages; and (iii) mycobacteriophage exclusive. We also determined the GOLLD secondary structure of each group using R2 R software based on GOLLD alignments generated by Infernal software. All GOLLD groups displayed a 3' half conserved structure, including utter E-loops pseudoknots substructures, also shared by non-Mycobacterium GOLLD while the 5' half motif was different among the groups. Here, we showed that the lncRNA GOLLD is widespread in Mycobacterium within tRNA arrays and corroborated the previously predicted GOLLD secondary structure.

Screening and verification of microRNA promoter methylation sites in hepatocellular carcinoma

The promoter methylation mode of microribonucleic acid (miRNA) plays a crucial role in the process of hepatocellular carcinoma (HCC). Therefore, the primary purpose of this study was to screen and verify the miRNA methylation sites associated with the overall survival (OS) and clinical characteristics of HCC patients. Methylation-related data were from the Cancer Genome Atlas (TCGA). R software was utilized to screen the methylation sites. The least absolute shrinkage and selection operator algorithm was utilized to develop the miRNA promoter methylation models. Then, methylation-specific polymerase chain reaction was performed with 146 HCC tissues to verify the accuracy of the vascular infiltration-related model. Additionally, we verified the functions of vascular infiltration-related miRNA by utilizing cells transfected with miR-199a-3p mimic. The model for predicting OS of HCC patients contained eight methylation sites. The Kaplan-Meier analysis suggested that the model could divide HCC patients into high- and low-risk groups (P < .0001). COX regression analysis suggested that the model (P < .001; 95% CI, 1.264-2.709) and T category (P < .001; 95% CI, 1.472-3.119) were independent risk factors for affecting OS of HCC patients. The model for predicting vascular infiltration, pathological grade, and clinical stage contained 7, 10, and 9 methylation sites respectively, with their area under the receiver operating characteristic curve (AUC) values 0.667, 0.745, and 0.725, respectively. The functional analysis suggested that miRNA methylation is involved in various biological processes such as WNT, MAPK, and mTOR signaling pathways. The accuracy of the vascular infiltration-related model was consistent with our previous bioinformatics assay. And upregulation of miR-199a-3p decreased migration and invasion abilities. The screened miRNA promoter methylation sites can be served as biomarkers for judging OS, vascular infiltration, pathology grade, and clinical stage. It can also provide new targets for improving the treatment and prognosis of HCC patients.

The Signal and the Noise: Characteristics of Antisense RNA in Complex Microbial Communities

This study systematically evaluated the global patterns of microbial antisense expression across various environments and provides a bird’s-eye view of general patterns observed across data sets, which can provide guidelines in our understanding of antisense expression as well as interpretation of metatranscriptomic data in general. This analysis highlights that in some environments, antisense expression from microbial communities can dominate over regular gene expression. We explored some potential drivers of antisense transcription, but more importantly, this study serves as a starting point, highlighting topics for future research and providing guidelines to include antisense expression in generic bioinformatic pipelines for metatranscriptomic data.

High-throughput sequencing has allowed unprecedented insight into the composition and function of complex microbial communities. With metatranscriptomics, it is possible to interrogate the transcriptomes of multiple organisms simultaneously to get an overview of the gene expression of the entire community. Studies have successfully used metatranscriptomics to identify and describe relationships between gene expression levels and community characteristics. However, metatranscriptomic data sets contain a rich suite of additional information that is just beginning to be explored. Here, we focus on antisense expression in metatranscriptomics, discuss the different computational strategies for handling it, and highlight the strengths but also potentially detrimental effects on downstream analysis and interpretation. We also analyzed the antisense transcriptomes of multiple genomes and metagenome-assembled genomes (MAGs) from five different data sets and found high variability in the levels of antisense transcription for individual species, which were consistent across samples. Importantly, we challenged the conceptual framework that antisense transcription is primarily the product of transcriptional noise and found mixed support, suggesting that the total observed antisense RNA in complex communities arises from the combined effect of unknown biological and technical factors. Antisense transcription can be highly informative, including technical details about data quality and novel insight into the biology of complex microbial communities.

IMPORTANCE This study systematically evaluated the global patterns of microbial antisense expression across various environments and provides a bird’s-eye view of general patterns observed across data sets, which can provide guidelines in our understanding of antisense expression as well as interpretation of metatranscriptomic data in general. This analysis highlights that in some environments, antisense expression from microbial communities can dominate over regular gene expression. We explored some potential drivers of antisense transcription, but more importantly, this study serves as a starting point, highlighting topics for future research and providing guidelines to include antisense expression in generic bioinformatic pipelines for metatranscriptomic data.

Regulatory Noncoding Small RNAs Are Diverse and Abundant in an Extremophilic Microbial Community

Regulatory small RNAs (sRNAs) play large-scale and essential roles in many cellular processes across all domains of life. Microbial sRNAs have been extensively studied in model organisms, but very little is known about the dynamics of sRNA synthesis and their roles in the natural environment. In this study, we discovered hundreds of intergenic (itsRNAs) and antisense (asRNAs) sRNAs expressed in an extremophilic microbial community inhabiting halite nodules (salt rocks) in the Atacama Desert. For this, we built SnapT, a new sRNA annotation pipeline that can be applied to any microbial community. We found asRNAs with expression levels negatively correlated with that of their overlapping putative target and itsRNAs that were conserved and significantly differentially expressed between 2 sampling time points. We demonstrated that we could perform target prediction and correlate expression levels between sRNAs and predicted target mRNAs at the community level. Functions of putative mRNA targets reflected the environmental challenges members of the halite communities were subjected to, including osmotic adjustments to a major rain event and competition for nutrients.IMPORTANCE Microorganisms in the natural world are found in communities, communicating and interacting with each other; therefore, it is essential that microbial regulatory mechanisms, such as gene regulation affected by small RNAs (sRNAs), be investigated at the community level. This work demonstrates that metatranscriptomic field experiments can link environmental variation with changes in RNA pools and have the potential to provide new insights into environmental sensing and responses in natural microbial communities through noncoding RNA-mediated gene regulation.

A Novel Noncoding RNA dsr11 Involved in Heat Stress Tolerance in Deinococcus radiodurans

Deinococcus radiodurans is an extremely resistant bacteria that has evolved masterful strategies to enable survival under various environmental stress conditions. Heat stress is a major environmental stress factor that can cause denaturation of proteins, membrane disruption, and oxidative stress. Previous studies have examined the mechanisms of the heat stress response by analyzing changes in protein levels; however, little is known about the role of small noncoding RNAs (ncRNAs), which are known to play important regulatory functions in bacteria during various environmental stress response. The ncRNA dsr11 of D. radiodurans was previously identified by RNA-seq and Northern blot. In this study, we showed that the transcription level of dsr11 was up-regulated 4.2-fold under heat stress by qRT-PCR analysis. Heat tolerance assay showed that deleting dsr11 significantly inhibited the viability under high temperature conditions. To assess the influence of dsr11 on the D. radiodurans transcriptome, 157 genes were found differentially expressed in the knock-out mutant by RNA-seq experiment. Combining RNA-seq and in silico analysis, we found that trmE (tRNA modification GTPase) and dr_0651 (arginase) were likely to be the direct targets of dsr11. Further microscale thermophoresis results demonstrated that dsr11 can directly bind to the mRNA of trmE and dr_0651. Our results indicated that dsr11 can enhance the tolerance to heat stress of D. radiodurans by binding to trmE and dr_0651 mRNA. Overall, these results extend our understanding of ncRNA regulation and provide new insights into the heat stress response in D. radiodurans.

Application of Transcriptomics to Enhance Early Diagnostics of Mycobacterial Infections, with an Emphasis on Mycobacterium avium ssp. paratuberculosis

Mycobacteria cause a wide variety of disease in human and animals. Species that infect ruminants include M. bovis and M. avium ssp. paratuberculosis (MAP). MAP is the causative agent of Johne’s disease in ruminants, which is a chronic granulomatous enteric infection that leads to severe economic losses worldwide. Characteristic of MAP infection is the long, latent phase in which intermittent shedding can take place, while diagnostic tests are unable to reliably detect an infection in this stage. This leads to unnoticed dissemination within herds and the presence of many undetected, silent carriers, which makes the eradication of Johne’s disease difficult. To improve the control of MAP infection, research is aimed at improving early diagnosis. Transcriptomic approaches can be applied to characterize host-pathogen interactions during infection, and to develop novel biomarkers using transcriptional profiles. Studies have focused on the identification of specific RNAs that are expressed in different infection stages, which will assist in the development and clinical implementation of early diagnostic tests.

A Haloarchaeal Small Regulatory RNA (sRNA) Is Essential for Rapid Adaptation to Phosphate Starvation Conditions

The haloarchaeon Haloferax volcanii contains nearly 2800 small non-coding RNAs (sRNAs). One intergenic sRNA, sRNA132, was chosen for a detailed characterization. A deletion mutant had a growth defect and thus underscored the importance of sRNA132. A microarray analysis identified the transcript of an operon for a phosphate-specific ABC transporter as a putative target of sRNA132. Both the sRNA132 and the operon transcript accumulated under low phosphate concentrations, indicating a positive regulatory role of sRNA132. A kinetic analysis revealed that sRNA132 is essential shortly after the onset of phosphate starvation, while other regulatory processes take over after several hours. Comparison of the transcriptomes of wild-type and the sRNA132 gene deletion mutant 30 min after the onset of phosphate starvation revealed that sRNA132 controls a regulon of about 40 genes. Remarkably, the regulon included a second operon for a phosphate-specific ABC transporter, which also depended on sRNA132 for rapid induction in the absence of phosphate. Competitive growth experiments of the wild-type and ABC transporter operon deletion mutants underscored the importance of both transporters for growth at low phosphate concentrations. Northern blot analyses of four additional members of the sRNA132 regulon verified that all four transcripts depended on sRNA132 for rapid regulation after the onset of phosphate starvation. Importantly, this is the first example for the transient importance of a sRNA for any archaeal and bacterial species. In addition, this study unraveled the first sRNA regulon for haloarchaea.

Characterization of the transcriptome of Haloferax volcanii, grown under four different conditions, with mixed RNA-Seq

Haloferax volcanii is a well-established model species for haloarchaea. Small scale RNomics and bioinformatics predictions were used to identify small non-coding RNAs (sRNAs), and deletion mutants revealed that sRNAs have important regulatory functions. A recent dRNA-Seq study was used to characterize the primary transcriptome. Unexpectedly, it was revealed that, under optimal conditions, H. volcanii contains more non-coding sRNAs than protein-encoding mRNAs. However, the dRNA-Seq approach did not contain any length information. Therefore, a mixed RNA-Seq approach was used to determine transcript length and to identify additional transcripts, which are not present under optimal conditions. In total, 50 million paired end reads of 150 nt length were obtained. 1861 protein-coding RNAs (cdRNAs) were detected, which encoded 3092 proteins. This nearly doubled the coverage of cdRNAs, compared to the previous dRNA-Seq study. About 2/3 of the cdRNAs were monocistronic, and 1/3 covered more than one gene. In addition, 1635 non-coding sRNAs were identified. The highest fraction of non-coding RNAs were cis antisense RNAs (asRNAs). Analysis of the length distribution revealed that sRNAs have a median length of about 150 nt. Based on the RNA-Seq and dRNA-Seq results, genes were chosen to exemplify characteristics of the H. volcanii transcriptome by Northern blot analyses, e.g. 1) the transcript patterns of gene clusters can be straightforward, but also very complex, 2) many transcripts differ in expression level under the four analyzed conditions, 3) some genes are transcribed into RNA isoforms of different length, which can be differentially regulated, 4) transcripts with very long 5'-UTRs and with very long 3'-UTRs exist, and 5) about 30% of all cdRNAs have overlapping 3'-ends, which indicates, together with the asRNAs, that H. volcanii makes ample use of sense-antisense interactions. Taken together, this RNA-Seq study, together with a previous dRNA-Seq study, enabled an unprecedented view on the H. volcanii transcriptome.

The Primary Antisense Transcriptome of Halobacterium salinarum NRC-1

Antisense RNAs (asRNAs) are present in diverse organisms and play important roles in gene regulation. In this work, we mapped the primary antisense transcriptome in the halophilic archaeon Halobacterium salinarum NRC-1. By reanalyzing publicly available data, we mapped antisense transcription start sites (aTSSs) and inferred the probable 3′ ends of these transcripts. We analyzed the resulting asRNAs according to the size, location, function of genes on the opposite strand, expression levels and conservation. We show that at least 21% of the genes contain asRNAs in H. salinarum. Most of these asRNAs are expressed at low levels. They are located antisense to genes related to distinctive characteristics of H. salinarum, such as bacteriorhodopsin, gas vesicles, transposases and other important biological processes such as translation. We provide evidence to support asRNAs in type II toxin–antitoxin systems in archaea. We also analyzed public Ribosome profiling (Ribo-seq) data and found that ~10% of the asRNAs are ribosome-associated non-coding RNAs (rancRNAs), with asRNAs from transposases overrepresented. Using a comparative transcriptomics approach, we found that ~19% of the asRNAs annotated in H. salinarum belong to genes with an ortholog in Haloferax volcanii, in which an aTSS could be identified with positional equivalence. This shows that most asRNAs are not conserved between these halophilic archaea.

A journey through the evolutionary diversification of archaeal Lsm and Hfq proteins

Sm-like (Lsm) proteins are found in all three domains of life. They are crucially involved in the RNA metabolism of prokaryotic organisms. To exert their function, they assemble into hexa- or heptameric rings and bind RNA via a conserved binding pocket for uridine stretches in the inner pore of the ring. Despite the conserved secondary structure of Lsm proteins, there are several features that lead to a structural diversification of this protein family that mediates their participation in a variety of processes related to RNA metabolism. Until recently, the cellular function of archaeal Sm-like proteins was not well understood. In this review, we discuss structural features of Lsm proteins with a strong focus on archaeal variants, reflect on the evolutionary development of archaeal Lsm proteins and present recent insights into their biological function.

A Small Regulatory RNA Controls Cell Wall Biosynthesis and Antibiotic Resistance

Small regulatory RNAs play an important role in the adaptation to changing conditions. Here, we describe a differentially expressed small regulatory RNA (sRNA) that affects various cellular processes in the plant pathogen Agrobacterium tumefaciens Using a combination of bioinformatic predictions and comparative proteomics, we identified nine targets, most of which are positively regulated by the sRNA. According to these targets, we named the sRNA PmaR for peptidoglycan biosynthesis, motility, and ampicillin resistance regulator. Agrobacterium spp. are long known to be naturally resistant to high ampicillin concentrations, and we can now explain this phenotype by the positive PmaR-mediated regulation of the beta-lactamase gene ampC Structure probing revealed a spoon-like structure of the sRNA, with a single-stranded loop that is engaged in target interaction in vivo and in vitro Several riboregulators have been implicated in antibiotic resistance mechanisms, such as uptake and efflux transporters, but PmaR represents the first example of an sRNA that directly controls the expression of an antibiotic resistance gene.IMPORTANCE The alphaproteobacterium Agrobacterium tumefaciens is able to infect various eudicots causing crown gall tumor formation. Based on its unique ability of interkingdom gene transfer, Agrobacterium serves as a crucial biotechnological tool for genetic manipulation of plant cells. The presence of hundreds of putative sRNAs in this organism suggests a considerable impact of riboregulation on A. tumefaciens physiology. Here, we characterized the biological function of the sRNA PmaR that controls various processes crucial for growth, motility, and virulence. Among the genes directly targeted by PmaR is ampC coding for a beta-lactamase that confers ampicillin resistance, suggesting that the sRNA is crucial for fitness in the competitive microbial composition of the rhizosphere.

mRNA-specific translation regulation by a ribosome-associated ncRNA in Haloferax volcanii

Regulation of gene expression at the translational level allows rapid adaptation of cellular proteomes to quickly changing environmental conditions and is thus central for prokaryotic organisms. Small non-coding RNAs (sRNAs) have been reported to effectively orchestrate translation control in bacteria and archaea mainly by targeting mRNAs by partial base complementarity. Here we report an unprecedented mechanism how sRNAs are capable of modulating protein biosynthesis in the halophilic archaeon Haloferax volcanii. By analyzing the ribosome-associated ncRNAs (rancRNAs) under different stress conditions we identified an intergenic sRNA, termed rancRNA_s194, that is primarily expressed during exponential growth under all tested conditions. By interaction with the ribosome rancRNA_s194 inhibits peptide bond formation and protein synthesis in vitro but appears to target a specific mRNA in vivo. The respective knock-out strain shows a reduced lag phase in media containing xylose as sole carbon source and outcompetes the wildtype cells under these conditions. Mass spectrometry, polysome profiling and mRNA binding competition experiments suggest that rancRNA_s194 prevents the cstA mRNA from being efficiently translated by H. volcanii ribosomes. These findings enlarge the regulatory repertoire of archaeal sRNAs in modulating post-transcriptional gene expression.