Global mapping transcriptional start sites revealed both transcriptional and post-transcriptional regulation of cold adaptation in the methanogenic archaeon Methanolobus psychrophilus

Promoter engineering for enhanced poly(3-hydroxybutyrate-co-3-hydroxyvalerate) (PHBV) production in Haloferax mediterranei

Polyhydroxyalkanoates (PHAs) are linear polyesters synthesized by various microorganisms. Due to their excellent physicochemical and material properties, PHAs are considered as a promising and sustainable alternative to petroleum-based plastics. Haloferax mediterranei is a promising and cost-effective producer of poly(3-hydroxybutyrate-co-3-hydroxyvalerate) (PHBV). In this study, we aimed to improve PHBV yield of H. mediterranei via promoter engineering. Differential RNA-seq (dRNA-seq) was firstly used to examine the genome-wide transcription start sites (TSSs) of H. mediterranei. Based on the TSS profile, 26 constitutive promoters with expression strength spanning ~103-fold dynamic range were quantitatively characterized using GFP fluorescence reporting system. Furthermore, it was found that the 5' untranslated region (5'-UTR) of phaR encoding PHA synthesis regulatory protein reduced the activity of the promoters derived from phaR, and this regulation occurred at the transcription level. Finally, several suitable promoters were selected to regulate the expression of cimA and bktB to improve the PHBV accumulation and 3HV molar fraction in H. mediterranei, ultimately leading to an increase of 35 % and of 18.4 %, respectively. Taken together, this work established a constitutive promoter library of haloarchaea, which would provide abundant genetic elements for the construction of ideal haloarchaeal chassis cells contributing in biosynthetic biology of extremophiles.

Structures of Saccharolobus solfataricus initiation complexes with leaderless mRNAs highlight archaeal features and eukaryotic proximity

The archaeal ribosome is of the eukaryotic type. TACK and Asgard superphyla, the closest relatives of eukaryotes, have ribosomes containing eukaryotic ribosomal proteins not found in other archaea, eS25, eS26 and eS30. Here, we investigate the case of Saccharolobus solfataricus, a TACK crenarchaeon, using mainly leaderless mRNAs. We characterize the small ribosomal subunit of S. solfataricus bound to SD-leadered or leaderless mRNAs. Cryo-EM structures show eS25, eS26 and eS30 bound to the small subunit. We identify two ribosomal proteins, aS33 and aS34, and an additional domain of eS6. Leaderless mRNAs are bound to the small subunit with contribution of their 5'-triphosphate group. Archaeal eS26 binds to the mRNA exit channel wrapped around the 3' end of rRNA, as in eukaryotes. Its position is not compatible with an SD:antiSD duplex. Our results suggest a positive role of eS26 in leaderless mRNAs translation and possible evolutionary routes from archaeal to eukaryotic translation.

A robust genetic toolbox for fine-tuning gene expression in the CO2-Fixing methanogenic archaeon Methanococcus maripaludis

Libraries of well-characterized genetic elements for fine-tuning gene expression are essential for biological and biotechnological research and applications. The fast-growing and genetically tractable methanogen, Methanococcus maripaludis, is a promising host organism for biotechnological conversion of carbon dioxide and renewable hydrogen into fuels and value-added products, as well as fundamental biological studies of archaea. However, the lack of molecular tools for gene expression has hindered its application as a workhorse to fine-tune gene and metabolic pathway expressions. In this study, we developed a genetic toolbox, including libraries of promoters, ribosome binding sites (RBS), and neutral sites for chromosomal integration, to facilitate precise gene expression in M. maripaludis. We generated a promoter library consisting of 81 constitutive promoters with expression strengths spanning a ∼104-fold dynamic range. Importantly, we identified a base composition rule for strong archaeal promoters and successfully remodeled weak promoters, enhancing their activities by up to 120-fold. We also established an RBS library containing 42 diverse RBS sequences with translation strengths covering a ∼100-fold dynamic range. Additionally, we identified eight neutral sites and developed a one-step, Cas9-based marker-less knock-in approach for chromosomal integration. We successfully applied the characterized promoter and RBS elements to significantly improve recombinant protein expression by 41-fold and modulate essential gene expression to generate corresponding physiological changes in M. maripaludis. Therefore, this work establishes a solid foundation for utilizing this autotrophic methanogen as an ideal workhorse for archaeal biology and biotechnological studies and applications.

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).

Endonucleolytic processing plays a critical role in the maturation of ribosomal RNA in Methanococcus maripaludis

Ribosomal RNA (rRNA) processing and maturation are fundamentally important for ribosome biogenesis, but the mechanisms in archaea, the third form of life, remains largely elusive. This study aimed to investigate the rRNA maturation process in Methanococcus maripaludis, a representative archaeon lacking known 3'-5' exonucleases. Through cleavage site identification and enzymatic assays, the splicing endonuclease EndA was determined to process the bulge-helix-bulge (BHB) motifs in 16S and 23S rRNA precursors. After splicing, the circular processing intermediates were formed and this was confirmed by quantitative RT-PCR and Northern blot. Ribonuclease assay revealed a specific cleavage at a 10-nt A/U-rich motif at the mature 5' end of pre-16S rRNA, which linearized circular pre-16S rRNA intermediate. Further 3'-RACE and ribonuclease assays determined that the endonuclease Nob1 cleaved the 3' extension of pre-16S rRNA, and so generated the mature 3' end. Circularized RT-PCR (cRT-PCR) and 5'-RACE identified two cleavage sites near helix 1 at the 5' end of 23S rRNA, indicating that an RNA structure-based endonucleolytic processing linearized the circular pre-23S rRNA intermediate. In the maturation of pre-5S rRNA, multiple endonucleolytic processing sites were determined at the 10-nt A/U-rich motif in the leader and trailer sequence. This study demonstrates that endonucleolytic processing, particularly at the 10-nt A/U-rich motifs play an essential role in the pre-rRNA maturation of M. maripaludis, indicating diverse pathways of rRNA maturation in archaeal species.

High-Resolution Small RNAs Landscape Provides Insights into Alkane Adaptation in the Marine Alkane-Degrader Alcanivorax dieselolei B-5

Alkanes are widespread in the ocean, and Alcanivorax is one of the most ubiquitous alkane-degrading bacteria in the marine ecosystem. Small RNAs (sRNAs) are usually at the heart of regulatory pathways, but sRNA-mediated alkane metabolic adaptability still remains largely unknown due to the difficulties of identification. Here, differential RNA sequencing (dRNA-seq) modified with a size selection (~50-nt to 500-nt) strategy was used to generate high-resolution sRNAs profiling in the model species Alcanivorax dieselolei B-5 under alkane (n-hexadecane) and non-alkane (acetate) conditions. As a result, we identified 549 sRNA candidates at single-nucleotide resolution of 5'-ends, 63.4% of which are with transcription start sites (TSSs), and 36.6% of which are with processing sites (PSSs) at the 5'-ends. These sRNAs originate from almost any location in the genome, regardless of intragenic (65.8%), antisense (20.6%) and intergenic (6.2%) regions, and RNase E may function in the maturation of sRNAs. Most sRNAs locally distribute across the 15 reference genomes of Alcanivorax, and only 7.5% of sRNAs are broadly conserved in this genus. Expression responses to the alkane of several core conserved sRNAs, including 6S RNA, M1 RNA and tmRNA, indicate that they may participate in alkane metabolisms and result in more actively global transcription, RNA processing and stresses mitigation. Two novel CsrA-related sRNAs are identified, which may be involved in the translational activation of alkane metabolism-related genes by sequestering the global repressor CsrA. The relationships of sRNAs with the characterized genes of alkane sensing (ompS), chemotaxis (mcp, cheR, cheW2), transporting (ompT1, ompT2, ompT3) and hydroxylation (alkB1, alkB2, almA) were created based on the genome-wide predicted sRNA-mRNA interactions. Overall, the sRNA landscape lays the ground for uncovering cryptic regulations in critical marine bacterium, among which both the core and species-specific sRNAs are implicated in the alkane adaptive metabolisms.

Genomic remnants of ancestral methanogenesis and hydrogenotrophy in Archaea drive anaerobic carbon cycling

Anaerobic methane metabolism is among the hallmarks of Archaea, originating very early in their evolution. Here, we show that the ancestor of methane metabolizers was an autotrophic CO₂-reducing hydrogenotrophic methanogen that possessed the two main complexes, methyl-CoM reductase (Mcr) and tetrahydromethanopterin-CoM methyltransferase (Mtr), the anaplerotic hydrogenases Eha and Ehb, and a set of other genes collectively called "methanogenesis markers" but could not oxidize alkanes. Overturning recent inferences, we demonstrate that methyl-dependent hydrogenotrophic methanogenesis has emerged multiple times independently, either due to a loss of Mtr while Mcr is inherited vertically or from an ancient lateral acquisition of Mcr. Even if Mcr is lost, Mtr, Eha, Ehb, and the markers can persist, resulting in mixotrophic metabolisms centered around the Wood-Ljungdahl pathway. Through their methanogenesis remnants, Thorarchaeia and two newly reconstructed order-level lineages in Archaeoglobi and Bathyarchaeia act as metabolically versatile players in carbon cycling of anoxic environments across the globe.

Oxidative stress and antioxidant mechanisms of obligate anaerobes involved in biological waste treatment processes: A review

In-depth understanding of the molecular mechanisms and physiological consequences of oxidative stress is still limited for anaerobes. Anaerobic biotechnology has become widely accepted by the wastewater/sludge industry as a better alternative to more conventional but costly aerobic processes. However, the functional anaerobic microorganisms used in anaerobic biotechnology are frequently hampered by reactive oxygen/nitrogen species (ROS/RNS)-mediated oxidative stress caused by exposure to stressful factors (e.g., oxygen and heavy metals), which negatively impact treatment performance. Thus, identifying stressful factors and understanding antioxidative defense mechanisms of functional obligate anaerobes are crucial for the optimization of anaerobic bioprocesses. Herein, we present a comprehensive overview of oxidative stress and antioxidant mechanisms of obligate anaerobes involved in anaerobic bioprocesses; as examples, we focus on anaerobic ammonium oxidation bacteria and methanogenic archaea. We summarize the primary stress factors in anaerobic bioprocesses and the cellular antioxidant defense systems of functional anaerobes, a consortia of enzymatic and nonenzymatic mechanisms. The dual role of ROS/RNS in cellular processes is elaborated; at low concentrations, they have vital cell signaling functions, but at high concentrations, they cause oxidative damage. Finally, we highlight gaps in knowledge and future work to uncover antioxidant and damage repair mechanisms in obligate anaerobes. This review provides in-depth insights and guidance for future research on oxidative stress of obligate anaerobes to boost the accurate regulation of anaerobic bioprocesses in challenging and changing operating conditions.

Lysobacter antarcticus sp. nov., an SUF-system-containing bacterium from Antarctic coastal sediment

A Gram-stain-negative, heterotrophic, aerobic, non-motile, rod-shaped bacterial strain (GW1-59T) belonging to the genus Lysobacter was isolated from coastal sediment collected from the Chinese Great Wall Station, Antarctica. The strain was identified using a polyphasic taxonomic approach. The strain grew well on Reasoner's 2A media and could grow in the presence of 0–4 % (w/v) NaCl (optimum, 1 %), at pH 9.0–11.0 and at 15–37 °C (optimum, 30 °C). Strain GW1-59T possessed ubiquinone-8 as the sole respiratory quinone. The major phospholipids were diphosphatidylglycerol, phosphatidylglycerol and phosphatidylethanolamine. The major fatty acids were summed feature 9 (10-methyl C16 : 0 and/or iso-C17 : 1  ω9c), iso-C15 : 0, iso-C16 : 0, iso-C17 : 0, C16 : 0 and iso-C11 : 0 3-OH. DNA–DNA relatedness with Lysobacter concretionis Ko07T, the nearest phylogenetic relative (98.5 % 16S rRNA gene sequence similarity) was 23.4 % (21.1–25.9 %). The average nucleotide identity value between strain GW1-59T and L. concretionis Ko07T was 80.1 %. The physiological and biochemical results and low level of DNA–DNA relatedness suggested the phenotypic and genotypic differentiation of strain GW1-59T from other Lysobacter species. On the basis of phenotypic, phylogenetic and genotypic data, a novel species, Lysobacter antarcticus sp. nov., is proposed. The type strain is GW1-59T (=CCTCC AB 2019390T=KCTC 72831T).

The Archaeal Small Heat Shock Protein Hsp17.6 Protects Proteins from Oxidative Inactivation

Small heat shock proteins (sHsps) are widely distributed among various types of organisms and function in preventing the irreversible aggregation of thermal denaturing proteins. Here, we report that Hsp17.6 from Methanolobus psychrophilus exhibited protection of proteins from oxidation inactivation. The overexpression of Hsp17.6 in Escherichia coli markedly increased the stationary phase cell density and survivability in HClO and H₂O₂. Treatments with 0.2 mM HClO or 10 mM H₂O₂ reduced malate dehydrogenase (MDH) activity to 57% and 77%, whereas the addition of Hsp17.6 recovered the activity to 70-90% and 86-100%, respectively. A similar effect for superoxide dismutase oxidation was determined for Hsp17.6. Non-reducing sodium dodecyl sulfate polyacrylamide gel electrophoresis assays determined that the Hsp17.6 addition decreased H₂O₂-caused disulfide-linking protein contents and HClO-induced degradation of MDH; meanwhile, Hsp17.6 protein appeared to be oxidized with increased molecular weights. Mass spectrometry identified oxygen atoms introduced into the larger Hsp17.6 molecules, mainly at the aspartate and methionine residues. Substitution of some aspartate residues reduced Hsp17.6 in alleviating H₂O₂- and HClO-caused MDH inactivation and in enhancing the E. coli survivability in H₂O₂ and HClO, suggesting that the archaeal Hsp17.6 oxidation protection might depend on an "oxidant sink" effect, i.e., to consume the oxidants in environments via aspartate oxidation.

Post-transcriptional regulation is involved in the cold-active methanol-based methanogenic pathway of a psychrophilic methanogen

The methanol-derived methanogenetic pathway contributes to bulk methane production in cold regions, but the cold adaptation mechanisms are obscure. This work investigated the mechanisms using a psychrophilic methylotrophic methanogen Methanolobus psychrophilus R15. R15 possesses two mtaCB operon paralogues-encoding methanol:corrinoid methyltransferase that is key to methanol-based methanogenesis. Molecular combined methanogenic assays determined that MtaC1 is important in methanogenesis at the optimal temperature of 18°C, but MtaC2 can be a cold-adaptive paralogue by highly upregulated at 8°C. The 5'P-seq and 5'RACE all assayed that processing occurred at the 5' untranslated region (5'-UTR) of mtaC2; reporter genes detected higher protein expression, and RNA half-life experiments assayed prolonged lifespan of the processed transcript. Therefore, mtaC2 5'-UTR processing to move the bulged structure elevated both the translation efficiency and transcript stability. 5'P-seq, quantitative RT-PCR and northern blot all identified enhanced mtaC2 5'-UTR processing at 8°C, which could contribute to the upregulation of mtaC2 at cold. The R15 cell extract contains an endoribonuclease cleaving an identified 10 nt-processing motif and the native mtaC2 5'-UTR particularly folded at 8°C. Therefore, this study revealed a 5'-UTR processing mediated post-transcriptional regulation mechanism controlling the cold-adaptive methanol-supported methanogenetic pathway, which may be used by other methylotrophic methanogens.

Integration of absolute multi-omics reveals dynamic protein-to-RNA ratios and metabolic interplay within mixed-domain microbiomes

While the field of microbiology has adapted to the study of complex microbiomes via modern meta-omics techniques, we have not updated our basic knowledge regarding the quantitative levels of DNA, RNA and protein molecules within a microbial cell, which ultimately control cellular function. Here we report the temporal measurements of absolute RNA and protein levels per gene within a mixed bacterial-archaeal consortium. Our analysis of this data reveals an absolute protein-to-RNA ratio of 10²–10⁴ for bacterial populations and 10³–10⁵ for an archaeon, which is more comparable to Eukaryotic representatives’ humans and yeast. Furthermore, we use the linearity between the metaproteome and metatranscriptome over time to identify core functional guilds, hence using a fundamental biological feature (i.e., RNA/protein levels) to highlight phenotypical complementarity. Our findings show that upgrading multi-omic toolkits with traditional absolute measurements unlocks the scaling of core biological questions to dynamic and complex microbiomes, creating a deeper insight into inter-organismal relationships that drive the greater community function.

Here, the authors perform a temporal multi-omic analysis of a minimalistic cellulose-degrading and methane-producing consortium at the strain level and estimate protein-to-RNA ratios and RNA-protein dynamics of the community simultaneously over time.

Recent Advances in Archaeal Translation Initiation

Translation initiation (TI) allows accurate selection of the initiation codon on a messenger RNA (mRNA) and defines the reading frame. In all domains of life, translation initiation generally occurs within a macromolecular complex made up of the small ribosomal subunit, the mRNA, a specialized methionylated initiator tRNA, and translation initiation factors (IFs). Once the start codon is selected at the P site of the ribosome and the large subunit is associated, the IFs are released and a ribosome competent for elongation is formed. However, even if the general principles are the same in the three domains of life, the molecular mechanisms are different in bacteria, eukaryotes, and archaea and may also vary depending on the mRNA. Because TI mechanisms have evolved lately, their studies bring important information about the evolutionary relationships between extant organisms. In this context, recent structural data on ribosomal complexes and genome-wide studies are particularly valuable. This review focuses on archaeal translation initiation highlighting its relationships with either the eukaryotic or the bacterial world. Eukaryotic features of the archaeal small ribosomal subunit are presented. Ribosome evolution and TI mechanisms diversity in archaeal branches are discussed. Next, the use of leaderless mRNAs and that of leadered mRNAs having Shine-Dalgarno sequences is analyzed. Finally, the current knowledge on TI mechanisms of SD-leadered and leaderless mRNAs is detailed.

The Transcriptomic Landscape of Cupriavidus metallidurans CH34 Acutely Exposed to Copper

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

A newly identified duplex RNA unwinding activity of archaeal RNase J depends on processive exoribonucleolysis coupled steric occlusion by its structural archaeal loops

RNase J is a prokaryotic 5'-3' exo/endoribonuclease that functions in mRNA decay and rRNA maturation. Here, we report a novel duplex unwinding activity of mpy-RNase J, an archaeal RNase J from Methanolobus psychrophilus, which enables it to degrade duplex RNAs with hairpins up to 40 bp when linking a 5' single-stranded overhangs of ≥ 7 nt, corresponding to the RNA channel length. A 6-nt RNA-mpy-RNase J-S247A structure reveals the RNA-interacting residues and a steric barrier at the RNA channel entrance comprising two archaeal loops and two helices. Mutagenesis of the residues key to either exoribonucleolysis or RNA translocation diminished the duplex unwinding activity. Substitution of the residues in the steric barrier yielded stalled degradation intermediates at the duplex RNA regions. Thus, an exoribonucleolysis-driven and steric occlusion-based duplex unwinding mechanism was identified. The duplex unwinding activity confers mpy-RNase J the capability of degrading highly structured RNAs, including the bacterial REP RNA, and archaeal mRNAs, rRNAs, tRNAs, SRPs, RNase P and CD-box RNAs, providing an indicative of the potential key roles of mpy-RNase J in pleiotropic RNA metabolisms. Hydrolysis-coupled duplex unwinding activity was also detected in a bacterial RNase J, which may use a shared but slightly different unwinding mechanism from archaeal RNase Js, indicating that duplex unwinding is a common property of the prokaryotic RNase Js.

Comprehensive analysis of the pre-ribosomal RNA maturation pathway in a methanoarchaeon exposes the conserved circularization and linearization mode in archaea

The ribosomal RNA (rRNA) genes are generally organized as an operon and cotranscribed into a polycistronic precursor; therefore, processing and maturation of pre-rRNAs are essential for ribosome biogenesis. However, rRNA maturation pathways of archaea, particularly of methanoarchaea, are scarcely known. Here, we thoroughly elucidated the maturation pathway of the rRNA operon (16S-tRNA^Ala-23S-tRNA^Cys-5S) in Methanolobus psychrophilus, one representative of methanoarchaea. Enzymatic assay demonstrated that EndA, a tRNA splicing endoribonuclease, cleaved bulge-helix-bulge (BHB) motifs buried in the processing stems of pre-16S and pre-23S rRNAs. Northern blot and quantitative PCR detected splicing-coupled circularization of pre-16S and pre-23S rRNAs, which accounted for 2% and 12% of the corresponding rRNAs, respectively. Importantly, endoribonuclease Nob1 was determined to linearize circular pre-16S rRNA at the mature 3' end so to expose the anti-Shine-Dalgarno sequence, while circular pre-23S rRNA was linearized at the mature 5' end by an unknown endoribonuclease. The resultant 5' and 3' extension in linearized pre-16S and pre-23S rRNAs were finally matured through 5'-3' and 3'-5' exoribonucleolytic trimming, respectively. Additionally, a novel processing pathway of endoribonucleolysis coupled with exoribonucleolysis was identified for the pre-5S rRNA maturation in this methanogen, which could be also conserved in most methanogenic euryarchaea. Based on evaluating the phylogenetic conservation of the key elements that are involved in circularization and linearization of pre-rRNA maturation, we predict that the rRNA maturation mode revealed here could be prevalent among archaea.

A novel small RNA contributes to restrain cellular chain length and anti-phagocytic ability in Streptococcus suis 2

Streptococcus suis serotype 2 (SS2) is an important porcine and human pathogen. Regulatory small non-coding RNAs (sRNAs) play an essential role in diverse physiological processes, although they remain poorly understood in SS2. In this study, we identified eight novel sRNAs through a combination of computational strategies and experimental identification. To explore roles of these novel sRNAs, sRNA34 was preferentially selected to assess phenotypes of the deletion strain in vitro and in vivo. The inactivation of sRNA34 significantly elongated the cellular chain, remarkably increased sensitivity to phagocytosis by RAW264.7, and attenuated virulence in a mouse infection model. Transcriptomic analysis revealed that inactivation of sRNA34 altered expression of multiple genes contributing to cellular chain formation and elongation, indicating a potential mechanism of sRNA34 in maintaining proper bacterial chain length to resist phagocytosis by the host cell. In summary, sRNA34 is a novel sRNA that contributes to cellular chain regulation and the anti-phagocytosis ability of SS2.

The archaeal RNA chaperone TRAM0076 shapes the transcriptome and optimizes the growth of Methanococcus maripaludis

TRAM is a conserved domain among RNA modification proteins that are widely distributed in various organisms. In Archaea, TRAM occurs frequently as a standalone protein with in vitro RNA chaperone activity; however, its biological significance and functional mechanism remain unknown. This work demonstrated that TRAM0076 is an abundant standalone TRAM protein in the genetically tractable methanoarcheaon Methanococcus maripaludis. Deletion of MMP0076, the gene encoding TRAM0076, markedly reduced the growth and altered transcription of 55% of the genome. Substitution mutations of Phe39, Phe42, Phe63, Phe65 and Arg35 in the recombinant TRAM0076 decreased the in vitro duplex RNA unfolding activity. These mutations also prevented complementation of the growth defect of the MMP0076 deletion mutant, indicating that the duplex RNA unfolding activity was essential for its physiological function. A genome-wide mapping of transcription start sites identified many 5' untranslated regions (5'UTRs) of 20-60 nt which could be potential targets of a RNA chaperone. TRAM0076 unfolded three representative 5'UTR structures in vitro and facilitated the in vivo expression of a mCherry reporter system fused to the 5'UTRs, thus behaving like a transcription anti-terminator. Flag-tagged-TRAM0076 co-immunoprecipitated a large number of cellular RNAs, suggesting that TRAM0076 plays multiple roles in addition to unfolding incorrect RNA structures. This work demonstrates that the conserved archaeal RNA chaperone TRAM globally affects gene expression and may represent a transcriptional element in ancient life of the RNA world.

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 computational system for identifying operons based on RNA-seq data

An operon is a set of neighboring genes in a genome that is transcribed as a single polycistronic message. Genes that are part of the same operon often have related functional roles or participate in the same metabolic pathways. The majority of all bacterial genes are co-transcribed with one or more other genes as part of a multi-gene operon. Thus, accurate identification of operons is important in understanding co-regulation of genes and their functional relationships. Here, we present a computational system that uses RNA-seq data to determine operons throughout a genome. The system takes the name of a genome and one or more files of RNA-seq data as input. Our method combines primary genomic sequence information with expression data from the RNA-seq files in a unified probabilistic model in order to identify operons. We assess our method's ability to accurately identify operons in a range of species through comparison to external databases of operons, both experimentally confirmed and computationally predicted, and through focused experiments that confirm new operons identified by our method. Our system is freely available at https://cs.wellesley.edu/~btjaden/Rockhopper/.

Wide diversity of methane and short-chain alkane metabolisms in uncultured archaea

Methanogenesis is an ancient metabolism of key ecological relevance, with direct impact on the evolution of Earth’s climate. Recent results suggest that the diversity of methane metabolisms and their derivations have probably been vastly underestimated. Here, by probing thousands of publicly available metagenomes for homologues of methyl-coenzyme M reductase complex (MCR), we have obtained ten metagenome-assembled genomes (MAGs) belonging to potential methanogenic, anaerobic methanotrophic and short-chain alkane oxidizing archaea. Five of these MAGs represent under-sampled (e.g., Verstraetearchaeota, Methanonatronarchaeia, ANME-1) or previously genomically undescribed (ANME-2c) archaeal lineages. The remaining five MAGs correspond to lineages that are only distantly related to previously known methanogens and span the entire archaeal phylogeny. Comprehensive comparative annotation significantly expands the metabolic diversity and energy conservation systems of MCR-bearing archaea. It also suggests the potential existence of a yet uncharacterized type of methanogenesis linked to short-chain alkane/fatty acid oxidation in a previously undescribed class of archaea (‘Ca. Methanoliparia’). We redefine a common core of marker genes specific to methanogenic, anaerobic methanotrophic and short-chain alkane-oxidizing archaea, and propose a possible scenario for the evolutionary and functional transitions that led to the emergence of such metabolic diversity.

Insights into RNA-processing pathways and associated RNA-degrading enzymes in Archaea

RNA-processing pathways are at the centre of regulation of gene expression. All RNA transcripts undergo multiple maturation steps in addition to covalent chemical modifications to become functional in the cell. This includes destroying unnecessary or defective cellular RNAs. In Archaea, information on mechanisms by which RNA species reach their mature forms and associated RNA-modifying enzymes are still fragmentary. To date, most archaeal actors and pathways have been proposed in light of information gathered from Bacteria and Eukarya. In this context, this review provides a state of the art overview of archaeal endoribonucleases and exoribonucleases that cleave and trim RNA species and also of the key small archaeal proteins that bind RNAs. Furthermore, synthetic up-to-date views of processing and biogenesis pathways of archaeal transfer and ribosomal RNAs as well as of maturation of stable small non-coding RNAs such as CRISPR RNAs, small C/D and H/ACA box guide RNAs, and other emerging classes of small RNAs are described. Finally, prospective post-transcriptional mechanisms to control archaeal messenger RNA quality and quantity are discussed.

Genome-scale analysis of Acetobacterium bakii reveals the cold adaptation of psychrotolerant acetogens by post-transcriptional regulation

Acetogens synthesize acetyl-CoA via CO2 or CO fixation, producing organic compounds. Despite their ecological and industrial importance, their transcriptional and post-transcriptional regulation has not been systematically studied. With completion of the genome sequence of Acetobacterium bakii (4.28-Mb), we measured changes in the transcriptome of this psychrotolerant acetogen in response to temperature variations under autotrophic and heterotrophic growth conditions. Unexpectedly, acetogenesis genes were highly up-regulated at low temperatures under heterotrophic, as well as autotrophic, growth conditions. To mechanistically understand the transcriptional regulation of acetogenesis genes via changes in RNA secondary structures of 5'-untranslated regions (5'-UTR), the primary transcriptome was experimentally determined, and 1379 transcription start sites (TSS) and 1100 5'-UTR were found. Interestingly, acetogenesis genes contained longer 5'-UTR with lower RNA-folding free energy than other genes, revealing that the 5'-UTRs control the RNA abundance of the acetogenesis genes under low temperature conditions. Our findings suggest that post-transcriptional regulation via RNA conformational changes of 5'-UTRs is necessary for cold-adaptive acetogenesis.

BAC-BROWSER: The Tool for Visualization and Analysis of Prokaryotic Genomes

Prokaryotes are actively studied objects in the scope of genomic regulation. Microbiologists need special tools for complex analysis of data to study and identification of regulatory mechanism in bacteria and archaea. We developed a tool BAC-BROWSER, specifically for visualization and analysis of small prokaryotic genomes. BAC-BROWSER provides tools for different types of analysis to study a wide set of regulatory mechanisms of prokaryotes: -transcriptional regulation by transcription factors (TFs), analysis of TFs, their targets, and binding sites.-other regulatory motifs, promoters, terminators and ribosome binding sites-transcriptional regulation by variation of operon structure, alternative starts or ends of transcription.-non-coding RNAs, antisense RNAs-RNA secondary structure, riboswitches-GC content, GC skew, codon usage BAC-browser incorporated free programs accelerating the verification of obtained results: primer design and oligocalculator, vector visualization, the tool for synthetic gene construction. The program is designed for Windows operating system and freely available for download in http://smdb.rcpcm.org/tools/index.html.

Comparative proteomic analysis of Methanothermobacter thermautotrophicus reveals methane formation from H2 and CO2 under different temperature conditions

The growth of all methanogens is limited to a specific temperature range. However, Methanothermobacter thermautotrophicus can be found in a variety of natural and artificial environments, the temperatures of which sometimes even exceed the temperature growth ranges of thermophiles. As a result, the extent to which methane production and survival are affected by temperature remains unclear. To investigate the mechanisms of methanogenesis that Archaea have evolved to cope with drastic temperature shifts, the responses of Methanothermobacter thermautotrophicus to temperature were investigated under a high temperature growth (71°C) and cold shock (4°C) using Isobaric tags for relative and absolute quantitation (iTRAQ). The results showed that methane formation is decreased and that protein folding and degradation are increased in both high‐ and low‐temperature treatments. In addition, proteins predicted to be involved in processing environmental information processing and in cell membrane/wall/envelope biogenesis may play key roles in affecting methane formation and enhancing the response of M. thermautotrophicus to temperature stress. Analysis of the genomic locations of the genes corresponding to these temperature‐dependent proteins predicted that 77 of the genes likely to form 32 gene clusters. Here, we assess the response of M. thermautotrophicus to different temperatures and provide a new level of understanding of methane formation and cellular putative adaptive responses.

Fine-Tuned Protein Production in Methanosarcina acetivorans C2A

Methanogenic archaea can be integrated into a sustainable, carbon-neutral cycle for producing organic chemicals from C₁ compounds if the rate, yield, and titer of product synthesis can be improved using metabolic engineering. However, metabolic engineering techniques are limited in methanogens by insufficient methods for controlling cellular protein levels. We conducted a systematic approach to tune protein levels in Methanosarcina acetivorans C2A, a model methanogen, by regulating transcription and translation initiation. Rationally designed core promoter and ribosome binding site mutations in M. acetivorans C2A resulted in a predicable change in protein levels over a 60 fold range. The overall range of protein levels was increased an additional 3 fold by introducing the 5' untranslated region of the mcrB transcript. This work demonstrates a wide range of precisely controlled protein levels in M. acetivorans C2A, which will help facilitate systematic metabolic engineering efforts in methanogens.

Toward a Universal Structural and Energetic Model for Prokaryotic Promoters

With almost no consensus promoter sequence in prokaryotes, recruitment of RNA polymerase (RNAP) to precise transcriptional start sites (TSSs) has remained an unsolved puzzle. Uncovering the underlying mechanism is critical for understanding the principle of gene regulation. We attempted to search the hidden code in ∼16,500 promoters of 12 prokaryotes representing two kingdoms in their structure and energetics. Twenty-eight fundamental parameters of DNA structure including backbone angles, basepair axis, and interbasepair and intrabasepair parameters were used, and information was extracted from x-ray crystallography data. Three parameters (solvation energy, hydrogen-bond energy, and stacking energy) were selected for creating energetics profiles using in-house programs. DNA of promoter regions was found to be inherently designed to undergo a change in every parameter undertaken for the study, in all prokaryotes. The change starts from some distance upstream of TSSs and continues past some distance from TSS, hence giving a signature state to promoter regions. These signature states might be the universal hidden codes recognized by RNAP. This observation was reiterated when randomly selected promoter sequences (with little sequence conservation) were subjected to structure generation; all developed into very similar three-dimensional structures quite distinct from those of conventional B-DNA and coding sequences. Fine structural details at important motifs (viz. -11, -35, and -75 positions relative to TSS) of promoters reveal novel to our knowledge and pointed insights for RNAP interaction at these locations; it could be correlated with how some particular structural changes at the -11 region may allow insertion of RNAP amino acids in interbasepair space as well as facilitate the flipping out of bases from the DNA duplex.

Modeling leaderless transcription and atypical genes results in more accurate gene prediction in prokaryotes

In a conventional view of the prokaryotic genome organization, promoters precede operons and ribosome binding sites (RBSs) with Shine-Dalgarno consensus precede genes. However, recent experimental research suggesting a more diverse view motivated us to develop an algorithm with improved gene-finding accuracy. We describe GeneMarkS-2, an ab initio algorithm that uses a model derived by self-training for finding species-specific (native) genes, along with an array of precomputed "heuristic" models designed to identify harder-to-detect genes (likely horizontally transferred). Importantly, we designed GeneMarkS-2 to identify several types of distinct sequence patterns (signals) involved in gene expression control, among them the patterns characteristic for leaderless transcription as well as noncanonical RBS patterns. To assess the accuracy of GeneMarkS-2, we used genes validated by COG (Clusters of Orthologous Groups) annotation, proteomics experiments, and N-terminal protein sequencing. We observed that GeneMarkS-2 performed better on average in all accuracy measures when compared with the current state-of-the-art gene prediction tools. Furthermore, the screening of ∼5000 representative prokaryotic genomes made by GeneMarkS-2 predicted frequent leaderless transcription in both archaea and bacteria. We also observed that the RBS sites in some species with leadered transcription did not necessarily exhibit the Shine-Dalgarno consensus. The modeling of different types of sequence motifs regulating gene expression prompted a division of prokaryotic genomes into five categories with distinct sequence patterns around the gene starts.

The Non-Coding Regulatory RNA Revolution in Archaea

Small non-coding RNAs (sRNAs) are ubiquitously found in the three domains of life playing large-scale roles in gene regulation, transposable element silencing and defense against foreign elements. While a substantial body of experimental work has been done to uncover function of sRNAs in Bacteria and Eukarya, the functional roles of sRNAs in Archaea are still poorly understood. Recently, high throughput studies using RNA-sequencing revealed that sRNAs are broadly expressed in the Archaea, comprising thousands of transcripts within the transcriptome during non-challenged and stressed conditions. Antisense sRNAs, which overlap a portion of a gene on the opposite strand (cis-acting), are the most abundantly expressed non-coding RNAs and they can be classified based on their binding patterns to mRNAs (3' untranslated region (UTR), 5' UTR, CDS-binding). These antisense sRNAs target many genes and pathways, suggesting extensive roles in gene regulation. Intergenic sRNAs are less abundantly expressed and their targets are difficult to find because of a lack of complete overlap between sRNAs and target mRNAs (trans-acting). While many sRNAs have been validated experimentally, a regulatory role has only been reported for very few of them. Further work is needed to elucidate sRNA-RNA binding mechanisms, the molecular determinants of sRNA-mediated regulation, whether protein components are involved and how sRNAs integrate with complex regulatory networks.

Genome-scale analysis of Methicillin-resistant Staphylococcus aureus USA300 reveals a tradeoff between pathogenesis and drug resistance

Staphylococcus aureus infection is a rising public health care threat. S. aureus is believed to have elaborate regulatory networks that orchestrate its virulence. Despite its importance, the systematic understanding of the transcriptional landscape of S. aureus is limited. Here, we describe the primary transcriptome landscape of an epidemic USA300 isolate of community-acquired methicillin-resistant S. aureus. We experimentally determined 1,861 transcription start sites with their principal promoter elements, including well-conserved -35 and -10 elements and weakly conserved -16 element and 5' untranslated regions containing AG-rich Shine-Dalgarno sequence. In addition, we identified 225 genes whose transcription was initiated from multiple transcription start sites, suggesting potential regulatory functions at transcription level. Along with the transcription unit architecture derived by integrating the primary transcriptome analysis with operon prediction, the measurement of differential gene expression revealed the regulatory framework of the virulence regulator Agr, the SarA-family transcriptional regulators, and β-lactam resistance regulators. Interestingly, we observed a complex interplay between virulence regulation, β-lactam resistance, and metabolism, suggesting a possible tradeoff between pathogenesis and drug resistance in the USA300 strain. Our results provide platform resource for the location of transcription initiation and an in-depth understanding of transcriptional regulation of pathogenesis, virulence, and antibiotic resistance in S. aureus.

Benefit from decline: the primary transcriptome of Alteromonas macleodii str. Te101 during Trichodesmium demise

Interactions between co-existing microorganisms deeply affect the physiology of the involved organisms and, ultimately, the function of the ecosystem as a whole. Copiotrophic Alteromonas are marine gammaproteobacteria that thrive during the late stages of phytoplankton blooms in the marine environment and in laboratory co-cultures with cyanobacteria such as Trichodesmium. The response of this heterotroph to the sometimes rapid and transient changes in nutrient supply when the phototroph crashes is not well understood. Here, we isolated and sequenced the strain Alteromonas macleodii str. Te101 from a laboratory culture of Trichodesmium erythraeum IMS101, yielding a chromosome of 4.63 Mb and a single plasmid of 237 kb. Increasing salinities to ≥43 ppt inhibited the growth of Trichodesmium but stimulated growth of the associated Alteromonas. We characterized the transcriptomic responses of both microorganisms and identified the complement of active transcriptional start sites in Alteromonas at single-nucleotide resolution. In replicate cultures, a similar set of genes became activated in Alteromonas when growth rates of Trichodesmium declined and mortality was high. The parallel activation of fliA, rpoS and of flagellar assembly and growth-related genes indicated that Alteromonas might have increased cell motility, growth, and multiple biosynthetic activities. Genes with the highest expression in the data set were three small RNAs (Aln1a-c) that were identified as analogs of the small RNAs CsrB-C in E. coli or RsmX-Z in pathogenic bacteria. Together with the carbon storage protein A (CsrA) homolog Te101_05290, these RNAs likely control the expression of numerous genes in responding to changes in the environment.

Genome-wide mRNA processing in methanogenic archaea reveals post-transcriptional regulation of ribosomal protein synthesis

Unlike stable RNAs that require processing for maturation, prokaryotic cellular mRNAs generally follow an 'all-or-none' pattern. Herein, we used a 5΄ monophosphate transcript sequencing (5΄P-seq) that specifically captured the 5΄-end of processed transcripts and mapped the genome-wide RNA processing sites (PSSs) in a methanogenic archaeon. Following statistical analysis and stringent filtration, we identified 1429 PSSs, among which 23.5% and 5.4% were located in 5΄ untranslated region (uPSS) and intergenic region (iPSS), respectively. A predominant uridine downstream PSSs served as a processing signature. Remarkably, 5΄P-seq detected overrepresented uPSS and iPSS in the polycistronic operons encoding ribosomal proteins, and the majority upstream and proximal ribosome binding sites, suggesting a regulatory role of processing on translation initiation. The processed transcripts showed increased stability and translation efficiency. Particularly, processing within the tricistronic transcript of rplA-rplJ-rplL enhanced the translation of rplL, which can provide a driving force for the 1:4 stoichiometry of L10 to L12 in the ribosome. Growth-associated mRNA processing intensities were also correlated with the cellular ribosomal protein levels, thereby suggesting that mRNA processing is involved in tuning growth-dependent ribosome synthesis. In conclusion, our findings suggest that mRNA processing-mediated post-transcriptional regulation is a potential mechanism of ribosomal protein synthesis and stoichiometry.

Conserved TRAM Domain Functions as an Archaeal Cold Shock Protein via RNA Chaperone Activity

Cold shock proteins (Csps) enable organisms to acclimate to and survive in cold environments and the bacterial CspA family exerts the cold protection via its RNA chaperone activity. However, most Archaea do not contain orthologs to the bacterial csp. TRAM, a conserved domain among RNA modification proteins ubiquitously distributed in organisms, occurs as an individual protein in most archaeal phyla and has a structural similarity to Csp proteins, yet its biological functions remain unknown. Through physiological and biochemical studies on four TRAM proteins from a cold adaptive archaeon Methanolobus psychrophilus R15, this work demonstrated that TRAM is an archaeal Csp and exhibits RNA chaperone activity. Three TRAM encoding genes (Mpsy_0643, Mpsy_3043, and Mpsy_3066) exhibited remarkable cold-shock induced transcription and were preferentially translated at lower temperature (18°C), while the fourth (Mpsy_2002) was constitutively expressed. They were all able to complement the cspABGE mutant of Escherichia coli BX04 that does not grow in cold temperatures and showed transcriptional antitermination. TRAM3066 (gene product of Mpsy_3066) and TRAM2002 (gene product of Mpsy_2002) displayed sequence-non-specific RNA but not DNA binding activity, and TRAM3066 assisted RNases in degradation of structured RNA, thus validating the RNA chaperone activity of TRAMs. Given the chaperone activity, TRAM is predicted to function beyond a Csp.

Transcriptional profiling of Rickettsia prowazekii coding and non-coding transcripts during in vitro host-pathogen and vector-pathogen interactions

Natural pathogen transmission of Rickettsia prowazekii, the etiologic agent of epidemic typhus, to humans is associated with arthropods, including human body lice, ticks, and ectoparasites of eastern flying squirrel. Recently, we have documented the presence of small RNAs in Rickettsia species and expression of R. prowazekii sRNAs during infection of cultured human microvascular endothelial cells (HMECs), which represent the primary target cells during human infections. Bacterial noncoding transcripts are now well established as critical post-transcriptional regulators of virulence and adaptation mechanisms in varying host environments. Despite their importance, little is known about the expression profile and regulatory activities of R. prowazekii sRNAs (Rp_sRs) in different host cells encountered as part of the natural life-cycle. To investigate the sRNA expression profile of R. prowazekii during infection of arthropod host cells, we employed an approach combining in vitro infection, bioinformatics, RNA sequencing, and PCR-based quantitation. Global analysis of R. prowazekii transcriptome by strand-specific RNA sequencing enabled us to identify 67 cis-acting (antisense) and 26 trans-acting (intergenic) Rp_sRs expressed during the infection of Amblyomma americanum (AAE2) cells. Comparative evaluation of expression during R. prowazekii infection of HMECs and AAE2 cells by quantitative RT-PCR demonstrated significantly higher expression of four selected Rp_sRs in tick AAE2 cells. Examination of the coding transcriptome revealed differential up-regulation of >150 rickettsial genes in either HMECs or AAE2 cells and yielded evidence for host cell-dependent utilization of alternative transcription start sites by 18 rickettsial genes. Our results thus suggest noticeable differences in the expression of both Rp_sRs as well as the coding transcriptome and the exploitation of multiple transcription initiation sites for select genes during the infection of human endothelium and tick vector cells as the host and yield new insights into rickettsial virulence and transmission mechanisms.

A global analysis of transcription reveals two modes of Spt4/5 recruitment to archaeal RNA polymerase

The archaeal transcription apparatus is closely related to the eukaryotic RNA polymerase (RNAP) II system, while archaeal genomes are more similar to bacteria with densely packed genes organized in operons. This makes understanding transcription in archaea vital, both in terms of molecular mechanisms and evolution. Very little is known about how archaeal cells orchestrate transcription on a systems level. We have characterized the genome-wide occupancy of the Methanocaldococcus jannaschii transcription machinery and its transcriptome. Our data reveal how the TATA and BRE promoter elements facilitate recruitment of the essential initiation factors TATA-binding protein and transcription factor B, respectively, which in turn are responsible for the loading of RNAP into the transcription units. The occupancies of RNAP and Spt4/5 strongly correlate with each other and with RNA levels. Our results show that Spt4/5 is a general elongation factor in archaea as its presence on all genes matches RNAP. Spt4/5 is recruited proximal to the transcription start site on the majority of transcription units, while on a subset of genes, including rRNA and CRISPR loci, Spt4/5 is recruited to the transcription elongation complex during early elongation within 500 base pairs of the transcription start site and akin to its bacterial homologue NusG.

Genome-wide primary transcriptome analysis of H2-producing archaeon Thermococcus onnurineus NA1

In spite of their pivotal roles in transcriptional and post-transcriptional processes, the regulatory elements of archaeal genomes are not yet fully understood. Here, we determine the primary transcriptome of the H2-producing archaeon Thermococcus onnurineus NA1. We identified 1,082 purine-rich transcription initiation sites along with well-conserved TATA box, A-rich B recognition element (BRE), and promoter proximal element (PPE) motif in promoter regions, a high pyrimidine nucleotide content (T/C) at the -1 position, and Shine-Dalgarno (SD) motifs (GGDGRD) in 5' untranslated regions (5' UTRs). Along with differential transcript levels, 117 leaderless genes and 86 non-coding RNAs (ncRNAs) were identified, representing diverse cellular functions and potential regulatory functions under the different growth conditions. Interestingly, we observed low GC content in ncRNAs for RNA-based regulation via unstructured forms or interaction with other cellular components. Further comparative analysis of T. onnurineus upstream regulatory sequences with those of closely related archaeal genomes demonstrated that transcription of orthologous genes are initiated by highly conserved promoter sequences, however their upstream sequences for transcriptional and translational regulation are largely diverse. These results provide the genetic information of T. onnurineus for its future application in metabolic engineering.

ReadXplorer 2-detailed read mapping analysis and visualization from one single source

MOTIVATION

The vast amount of already available and currently generated read mapping data requires comprehensive visualization, and should benefit from bioinformatics tools offering a wide spectrum of analysis functionality from just one source. Appropriate handling of multiple mapped reads during mapping analyses remains an issue that demands improvement.

RESULTS

The capabilities of the read mapping analysis and visualization tool ReadXplorer were vastly enhanced. Here, we present an even finer granulated read mapping classification, improving the level of detail for analyses and visualizations. The spectrum of automatic analysis functions has been broadened to include genome rearrangement detection as well as correlation analysis between two mapping data sets. Existing functions were refined and enhanced, namely the computation of differentially expressed genes, the read count and normalization analysis and the transcription start site detection. Additionally, ReadXplorer 2 features a highly improved support for large eukaryotic data sets and a command line version, enabling its integration into workflows. Finally, the new version is now able to display any kind of tabular results from other bioinformatics tools.

AVAILABILITY AND IMPLEMENTATION

http://www.readxplorer.org CONTACT: readxplorer@computational.bio.uni-giessen.deSupplementary information: Supplementary data are available at Bioinformatics online.

Negative Regulation of Ectoine Uptake and Catabolism in Sinorhizobium meliloti: Characterization of the EhuR Gene

Ectoine has osmoprotective effects on Sinorhizobium meliloti that differ from its effects in other bacteria. Ectoine does not accumulate in S. meliloti cells; instead, it is degraded. The products of the ehuABCD-eutABCDE operon were previously discovered to be responsible for the uptake and catabolism of ectoine in S. meliloti However, the mechanism by which ectoine is involved in the regulation of the ehuABCD-eutABCDE operon remains unclear. The ehuR gene, which is upstream of and oriented in the same direction as the ehuABCD-eutABCDE operon, encodes a member of the MocR/GntR family of transcriptional regulators. Quantitative reverse transcription-PCR and promoter-lacZ reporter fusion experiments revealed that EhuR represses transcription of the ehuABCD-eutABCDE operon, but this repression is inhibited in the presence of ectoine. Electrophoretic mobility shift assays and DNase I footprinting assays revealed that EhuR bound specifically to the DNA regions overlapping the -35 region of the ehuA promoter and the +1 region of the ehuR promoter. Surface plasmon resonance assays further demonstrated direct interactions between EhuR and the two promoters, although EhuR was found to have higher affinity for the ehuA promoter than for the ehuR promoter. In vitro, DNA binding by EhuR could be directly inhibited by a degradation product of ectoine. Our work demonstrates that EhuR is an important negative transcriptional regulator involved in the regulation of ectoine uptake and catabolism and is likely regulated by one or more end products of ectoine catabolism.

IMPORTANCE

Sinorhizobium meliloti is an important soil bacterium that displays symbiotic interactions with legume hosts. Ectoine serves as a key osmoprotectant for S. meliloti However, ectoine does not accumulate in the cells; rather, it is degraded. In this study, we characterized the transcriptional regulation of the operon responsible for ectoine uptake and catabolism in S. meliloti We identified and characterized the transcription repressor EhuR, which is the first MocR/GntR family member found to be involved in the regulation of compatible solute uptake and catabolism. More importantly, we demonstrated for the first time that an ectoine catabolic end product could modulate EhuR DNA-binding activity. Therefore, this work provides new insights into the unique mechanism of ectoine-induced osmoprotection in S. meliloti.

Genome-wide gene expression and RNA half-life measurements allow predictions of regulation and metabolic behavior in Methanosarcina acetivorans

BACKGROUND

While a few studies on the variations in mRNA expression and half-lives measured under different growth conditions have been used to predict patterns of regulation in bacterial organisms, the extent to which this information can also play a role in defining metabolic phenotypes has yet to be examined systematically. Here we present the first comprehensive study for a model methanogen.

RESULTS

We use expression and half-life data for the methanogen Methanosarcina acetivorans growing on fast- and slow-growth substrates to examine the regulation of its genes. Unlike Escherichia coli where only small shifts in half-lives were observed, we found that most mRNA have significantly longer half-lives for slow growth on acetate compared to fast growth on methanol or trimethylamine. Interestingly, half-life shifts are not uniform across functional classes of enzymes, suggesting the existence of a selective stabilization mechanism for mRNAs. Using the transcriptomics data we determined whether transcription or degradation rate controls the change in transcript abundance. Degradation was found to control abundance for about half of the metabolic genes underscoring its role in regulating metabolism. Genes involved in half of the metabolic reactions were found to be differentially expressed among the substrates suggesting the existence of drastically different metabolic phenotypes that extend beyond just the methanogenesis pathways. By integrating expression data with an updated metabolic model of the organism (iST807) significant differences in pathway flux and production of metabolites were predicted for the three growth substrates.

CONCLUSIONS

This study provides the first global picture of differential expression and half-lives for a class II methanogen, as well as provides the first evidence in a single organism that drastic genome-wide shifts in RNA half-lives can be modulated by growth substrate. We determined which genes in each metabolic pathway control the flux and classified them as regulated by transcription (e.g. transcription factor) or degradation (e.g. post-transcriptional modification). We found that more than half of genes in metabolism were controlled by degradation. Our results suggest that M. acetivorans employs extensive post-transcriptional regulation to optimize key metabolic steps, and more generally that degradation could play a much greater role in optimizing an organism's metabolism than previously thought.

Genome-wide identification of transcriptional start sites in the haloarchaeon Haloferax volcanii based on differential RNA-Seq (dRNA-Seq)

BACKGROUND

Differential RNA-Seq (dRNA-Seq) is a recently developed method of performing primary transcriptome analyses that allows for the genome-wide mapping of transcriptional start sites (TSSs) and the identification of novel transcripts. Although the transcriptomes of diverse bacterial species have been characterized by dRNA-Seq, the transcriptome analysis of archaeal species is still rather limited. Therefore, we used dRNA-Seq to characterize the primary transcriptome of the model archaeon Haloferax volcanii.

RESULTS

Three independent cultures of Hfx. volcanii grown under optimal conditions to the mid-exponential growth phase were used to determine the primary transcriptome and map the 5'-ends of the transcripts. In total, 4749 potential TSSs were detected. A position weight matrix (PWM) was derived for the promoter predictions, and the results showed that 64 % of the TSSs were preceded by stringent or relaxed basal promoters. Of the identified TSSs, 1851 belonged to protein-coding genes. Thus, fewer than half (46 %) of the 4040 protein-coding genes were expressed under optimal growth conditions. Seventy-two percent of all protein-coding transcripts were leaderless, which emphasized that this pathway is the major pathway for translation initiation in haloarchaea. A total of 2898 of the TSSs belonged to potential non-coding RNAs, which accounted for an unexpectedly high fraction (61 %) of all transcripts. Most of the non-coding TSSs had not been previously described (2792) and represented novel sequences (59 % of all TSSs). A large fraction of the potential novel non-coding transcripts were cis-antisense RNAs (1244 aTSSs). A strong negative correlation between the levels of antisense transcripts and cognate sense mRNAs was found, which suggested that the negative regulation of gene expression via antisense RNAs may play an important role in haloarchaea. The other types of novel non-coding transcripts corresponded to internal transcripts overlapping with mRNAs (1153 iTSSs) and intergenic small RNA (sRNA) candidates (395 TSSs).

CONCLUSION

This study provides a comprehensive map of the primary transcriptome of Hfx. volcanii grown under optimal conditions. Fewer than half of all protein-coding genes have been transcribed under these conditions. Unexpectedly, more than half of the detected TSSs belonged to several classes of non-coding RNAs. Thus, RNA-based regulation appears to play a more important role in haloarchaea than previously anticipated.

Molecular Mechanisms of Transcription Initiation-Structure, Function, and Evolution of TFE/TFIIE-Like Factors and Open Complex Formation

Transcription initiation requires that the promoter DNA is melted and the template strand is loaded into the active site of the RNA polymerase (RNAP), forming the open complex (OC). The archaeal initiation factor TFE and its eukaryotic counterpart TFIIE facilitate this process. Recent structural and biophysical studies have revealed the position of TFE/TFIIE within the pre-initiation complex (PIC) and illuminated its role in OC formation. TFE operates via allosteric and direct mechanisms. Firstly, it interacts with the RNAP and induces the opening of the flexible RNAP clamp domain, concomitant with DNA melting and template loading. Secondly, TFE binds physically to single-stranded DNA in the transcription bubble of the OC and increases its stability. The identification of the β-subunit of archaeal TFE enabled us to reconstruct the evolutionary history of TFE/TFIIE-like factors, which is characterised by winged helix (WH) domain expansion in eukaryotes and loss of metal centres including iron-sulfur clusters and Zinc ribbons. OC formation is an important target for the regulation of transcription in all domains of life. We propose that TFE and the bacterial general transcription factor CarD, although structurally and evolutionary unrelated, show interesting parallels in their mechanism to enhance OC formation. We argue that OC formation is used as a way to regulate transcription in all domains of life, and these regulatory mechanisms coevolved with the basal transcription machinery.