Deep sequencing analysis of the Methanosarcina mazei Gö1 transcriptome in response to nitrogen availability

Small protein mediates inhibition of ammonium transport in Methanosarcina mazei-an ancient mechanism?

IMPORTANCE

Small proteins containing fewer than 70 amino acids, which were previously disregarded due to computational prediction and biochemical detection challenges, have gained increased attention in the scientific community in recent years. However, the number of functionally characterized small proteins, especially in archaea, is still limited. Here, by using biochemical and genetic approaches, we demonstrate a crucial role of the small protein sP36 in the nitrogen metabolism of M. mazei, which modulates the ammonium transporter AmtB1 according to nitrogen availability. This modulation might represent an ancient archaeal mechanism of AmtB1 inhibition, in contrast to the well-studied uridylylation-dependent regulation in bacteria.

Small proteins in bacteria - Big challenges in prediction and identification

Proteins with up to 100 amino acids have been largely overlooked due to the challenges associated with predicting and identifying them using traditional methods. Recent advances in bioinformatics and machine learning, DNA sequencing, RNA and Ribo-seq technologies, and mass spectrometry (MS) have greatly facilitated the detection and characterisation of these elusive proteins in recent years. This has revealed their crucial role in various cellular processes including regulation, signalling and transport, as toxins and as folding helpers for protein complexes. Consequently, the systematic identification and characterisation of these proteins in bacteria have emerged as a prominent field of interest within the microbial research community. This review provides an overview of different strategies for predicting and identifying these proteins on a large scale, leveraging the power of these advanced technologies. Furthermore, the review offers insights into the future developments that may be expected in this field.

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.

Newly Established Genetic System for Functional Analysis of MetSV

The linear chromosome of the Methanosarcina spherical virus with 10,567 bp exhibits 22 ORFs with mostly unknown functions. Annotation using common tools and databases predicted functions for a few genes like the type B DNA polymerase (MetSVORF07) or the small (MetSVORF15) and major (MetSVORF16) capsid proteins. For verification of assigned functions of additional ORFs, biochemical or genetic approaches were found to be essential. Consequently, we established a genetic system for MetSV by cloning its genome into the E. coli plasmid pCR-XL-2. Comparisons of candidate plasmids with the MetSV reference based on Nanopore sequencing revealed several mutations of yet unknown provenance with an impact on protein-coding sequences. Linear MetSV inserts were generated by BamHI restriction, purified and transformed in Methanosarcina mazei by an optimized liposome-mediated transformation protocol. Analysis of resulting MetSV virions by TEM imaging and infection experiments demonstrated no significant differences between plasmid-born viruses and native MetSV particles regarding their morphology or lytic behavior. The functionality of the genetic system was tested by the generation of a ΔMetSVORF09 mutant that was still infectious. Our genetic system of MetSV, the first functional system for a virus of methanoarchaea, now allows us to obtain deeper insights into MetSV protein functions and virus-host interactions.

Comparative Transcriptomics Sheds Light on Remodeling of Gene Expression during Diazotrophy in the Thermophilic Methanogen Methanothermococcus thermolithotrophicus

Some marine thermophilic methanogens are able to perform energy-consuming nitrogen fixation despite deriving only little energy from hydrogenotrophic methanogenesis. We studied this process in Methanothermococcus thermolithotrophicus DSM 2095, a methanogenic archaeon of the order Methanococcales that contributes to the nitrogen pool in some marine environments. We successfully grew this archaeon under diazotrophic conditions in both batch and fermenter cultures, reaching the highest cell density reported so far. Diazotrophic growth depended strictly on molybdenum and, in contrast to other diazotrophs, was not inhibited by tungstate or vanadium. This suggests an elaborate control of metal uptake and a specific metal recognition system for the insertion into the nitrogenase cofactor. Differential transcriptomics of M. thermolithotrophicus grown under diazotrophic conditions with ammonium-fed cultures as controls revealed upregulation of the nitrogenase machinery, including chaperones, regulators, and molybdate importers, as well as simultaneous upregulation of an ammonium transporter and a putative pathway for nitrate and nitrite utilization. The organism thus employs multiple synergistic strategies for uptake of nitrogen nutrients during the early exponential growth phase without altering transcription levels for genes involved in methanogenesis. As a counterpart, genes coding for transcription and translation processes were downregulated, highlighting the maintenance of an intricate metabolic balance to deal with energy constraints and nutrient limitations imposed by diazotrophy. This switch in the metabolic balance included unexpected processes, such as upregulation of the CRISPR-Cas system, probably caused by drastic changes in transcription levels of putative mobile and virus-like elements. IMPORTANCE The thermophilic anaerobic archaeon M. thermolithotrophicus is a particularly suitable model organism to study the coupling of methanogenesis to diazotrophy. Likewise, its capability of simultaneously reducing N₂ and CO₂ into NH₃ and CH₄ with H₂ makes it a viable target for biofuel production. We optimized M. thermolithotrophicus cultivation, resulting in considerably higher cell yields and enabling the successful establishment of N₂-fixing bioreactors. Improved understanding of the N₂ fixation process would provide novel insights into metabolic adaptations that allow this energy-limited extremophile to thrive under diazotrophy, for instance, by investigating its physiology and uncharacterized nitrogenase. We demonstrated that diazotrophic growth of M. thermolithotrophicus is exclusively dependent on molybdenum, and complementary transcriptomics corroborated the expression of the molybdenum nitrogenase system. Further analyses of differentially expressed genes during diazotrophy across three cultivation time points revealed insights into the response to nitrogen limitation and the coordination of core metabolic processes.

Dual-RNAseq Analysis Unravels Virus-Host Interactions of MetSV and Methanosarcina mazei

Methanosarcina spherical virus (MetSV), infecting Methanosarcina species, encodes 22 genes, but their role in the infection process in combination with host genes has remained unknown. To study the infection process in detail, infected and uninfected M. mazei cultures were compared using dual-RNAseq, qRT-PCRs, and transmission electron microscopy (TEM). The transcriptome analysis strongly indicates a combined role of virus and host genes in replication, virus assembly, and lysis. Thereby, 285 host and virus genes were significantly regulated. Within these 285 regulated genes, a network of the viral polymerase, MetSVORF6, MetSVORF5, MetSVORF2, and the host genes encoding NrdD, NrdG, a CDC48 family protein, and a SSB protein with a role in viral replication was postulated. Ultrastructural analysis at 180 min p.i. revealed many infected cells with virus particles randomly scattered throughout the cytoplasm or attached at the cell surface, and membrane fragments indicating cell lysis. Dual-RNAseq and qRT-PCR analyses suggested a multifactorial lysis reaction in potential connection to the regulation of a cysteine proteinase, a pirin-like protein and a HicB-solo protein. Our study's results led to the first preliminary infection model of MetSV infecting M. mazei, summarizing the key infection steps as follows: replication, assembly, and host cell lysis.

Understanding and application of Bacillus nitrogen regulation: A synthetic biology perspective

BACKGROUND

Nitrogen sources play an essential role in maintaining the physiological and biochemical activity of bacteria. Nitrogen metabolism, which is the core of microorganism metabolism, makes bacteria able to autonomously respond to different external nitrogen environments by exercising complex internal regulatory networks to help them stay in an ideal state. Although various studies have been put forth to better understand this regulation in Bacillus, and many valuable viewpoints have been obtained, these views need to be presented systematically and their possible applications need to be specified.

AIM OF REVIEW

The intention is to provide a deep and comprehensive understanding of nitrogen metabolism in Bacillus, an important industrial microorganism, and thereby apply this regulatory logic to synthetic biology to improve biosynthesis competitiveness. In addition, the potential researches in the future are also discussed.

KEY SCIENTIFIC CONCEPT OF REVIEW

Understanding the meticulous regulation process of nitrogen metabolism in Bacillus not only could facilitate research on metabolic engineering but also could provide constructive insights and inspiration for studies of other microorganisms.

RNAP II produces capped 18S and 25S ribosomal RNAs resistant to 5′-monophosphate dependent processive 5′ to 3′ exonuclease in polymerase switched Saccharomyces cerevisiae

Background

We have previously found that, in the pathogenic yeast Candida albicans, 18S and 25S ribosomal RNA components, containing more than one phosphate on their 5′-end were resistant to 5′-monophosphate requiring 5′ → 3″ exonuclease. Several lines of evidence pointed to RNAP II as the enzyme producing them.

Results

We now show the production of such 18S and 25S rRNAs in Saccharomyces cerevisiae that have been permanently switched to RNAP II (due to deletion of part of RNAP I upstream activator alone, or in combination with deletion of one component of RNAP I itself). They contain more than one phosphate at their 5′-end and an anti-cap specific antibody binds to them indicating capping of these molecules. These molecules are found in RNA isolated from nuclei, therefore are unlikely to have been modified in the cytoplasm.

Conclusions

Our data confirm the existence of such molecules and firmly establish RNAP II playing a role in their production. The fact that we see these molecules in wild type Saccharomyces cerevisiae indicates that they are not only a result of mutations but are part of the cells physiology. This adds another way RNAP II is involved in ribosome production in addition to their role in the production of ribosome associated proteins.

Small Proteins in Archaea, a Mainly Unexplored World

In recent years, increasing numbers of small proteins have moved into the focus of science. Small proteins have been identified and characterized in all three domains of life, but the majority remains functionally uncharacterized, lack secondary structure, and exhibit limited evolutionary conservation. While quite a few have already been described for bacteria and eukaryotic organisms, the amount of known and functionally analyzed archaeal small proteins is still very limited. In this review, we compile the current state of research, show strategies for systematic approaches for global identification of small archaeal proteins, and address selected functionally characterized examples. Besides, we document exemplarily for one archaeon the tool development and optimization to identify small proteins using genome-wide approaches.

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.

Understanding small ORF diversity through a comprehensive transcription feature classification

Small open reading frames (small ORFs/sORFs/smORFs) are potentially coding sequences smaller than 100 codons that have historically been considered junk DNA by gene prediction software and in annotation screening; however, the advent of next-generation sequencing has contributed to the deeper investigation of junk DNA regions and their transcription products, resulting in the emergence of smORFs as a new focus of interest in systems biology. Several smORF peptides were recently reported in noncanonical mRNAs as new players in numerous biological contexts; however, their relevance is still overlooked in coding potential analysis. Hence, this review proposes a smORF classification based on transcriptional features, discussing the most promising approaches to investigate smORFs based on their different characteristics. First, smORFs were divided into nonexpressed (intergenic) and expressed (genic) smORFs. Second, genic smORFs were classified as smORFs located in noncoding RNAs (ncRNAs) or canonical mRNAs. Finally, smORFs in ncRNAs were further subdivided into sequences located in small or long RNAs, whereas smORFs located in canonical mRNAs were subdivided into several specific classes depending on their localization along the gene. We hope that this review provides new insights into large-scale annotations and reinforces the role of smORFs as essential components of a hidden coding DNA world.

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.

Multi-protease Approach for the Improved Identification and Molecular Characterization of Small Proteins and Short Open Reading Frame-Encoded Peptides

The identification of proteins below approximately 70-100 amino acids in bottom-up proteomics is still a challenging task due to the limited number of peptides generated by proteolytic digestion. This includes the short open reading frame-encoded peptides (SEPs), which are a subset of the small proteins that were not previously annotated or that are alternatively encoded. Here, we systematically investigated the use of multiple proteases (trypsin, chymotrypsin, LysC, LysargiNase, and GluC) in GeLC-MS/MS analysis to improve the sequence coverage and the number of identified peptides for small proteins, with a focus on SEPs, in the archaeon Methanosarcina mazei. Combining the data of all proteases, we identified 63 small proteins and additional 28 SEPs with at least two unique peptides, while only 55 small proteins and 22 SEP could be identified using trypsin only. For 27 small proteins and 12 SEPs, a complete sequence coverage was achieved. Moreover, for five SEPs, incorrectly predicted translation start points or potential in vivo proteolytic processing were identified, confirming the data of a previous top-down proteomics study of this organism. The results show clearly that a multi-protease approach allows to improve the identification and molecular characterization of small proteins and SEPs. LC-MS data: ProteomeXchange PXD023921.

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.

Whole-genome Identification of Transcriptional Start Sites by Differential RNA-seq in Bacteria

Gene transcription in bacteria often starts some nucleotides upstream of the start codon. Identifying the specific Transcriptional Start Site (TSS) is essential for genetic manipulation, as in many cases upstream of the start codon there are sequence elements that are involved in gene expression regulation. Taken into account the classical gene structure, we are able to identify two kinds of transcriptional start site: primary and secondary. A primary transcriptional start site is located some nucleotides upstream of the translational start site, while a secondary transcriptional start site is located within the gene encoding sequence. Here, we present a step by step protocol for genome-wide transcriptional start sites determination by differential RNA-sequencing (dRNA-seq) using the enteric pathogen Shigella flexneri serotype 5a strain M90T as model. However, this method can be employed in any other bacterial species of choice. In the first steps, total RNA is purified from bacterial cultures using the hot phenol method. Ribosomal RNA (rRNA) is specifically depleted via hybridization probes using a commercial kit. A 5'-monophosphate-dependent exonuclease (TEX)-treated RNA library enriched in primary transcripts is then prepared for comparison with a library that has not undergone TEX-treatment, followed by ligation of an RNA linker adaptor of known sequence allowing the determination of TSS with single nucleotide precision. Finally, the RNA is processed for Illumina sequencing library preparation and sequenced as purchased service. TSS are identified by in-house bioinformatic analysis. Our protocol is cost-effective as it minimizes the use of commercial kits and employs freely available software.

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.

Multidimensional separation schemes enhance the identification and molecular characterization of low molecular weight proteomes and short open reading frame-encoded peptides in top-down proteomics

Short open reading frame-encoded peptides (SEP) represent a widely undiscovered part of the proteome. The detailed analysis of SEP has, despite inherent limitations such as incomplete sequence coverage, challenges encountered with protein inference, the identification of posttranslational modifications and the assignment of potential N- and C-terminal truncations, predominantly been assessed using bottom-up proteomic workflows. The use of top-down based proteomic workflows is capable of providing an unparalleled level of characterization information, which is of increased importance in the case of alternatively encoded protein products. However, top-down based analysis is not without its own limitations, for which efficient separation prior to MS analysis is a major issue. We established a sample preparation approach for the combined bottom-up and top-down proteomic analysis of SEP. Key improvements were made by the application of solid phase extraction (SPE), which supported enrichment of proteins below ca. 20 kDa, followed by 2D-LC-MS top-down analysis encompassing both HCD and EThcD ion activation. Bottom-up experiments were used to support and confirm top-down data interpretation. This strategy allowed for the top-down characterization of 36 proteoforms mapping to 12 SEP from the archaeon Methanosarcina mazei strain Gö1, with the concurrent detection and identification of several posttranslational modifications in SEP. BIOLOGICAL SIGNIFICANCE: Small or short open reading frames (sORF) have been widely neglected in genome research in the past. With their increasing discovery, the question about the presence and molecular function of their translation products, the short open reading frame-encoded peptides (SEP), arises. As these small proteins are usually below the 10 kDa range, the number of peptides identifiable by bottom-up proteomics is limited which hampers both the identification and the recognition of potential posttranslational modifications. The presented top-down approach allowed for the detection of full length SEP, as well as of terminally truncated proteoforms, and further enabled the identification of disulfide bonds in these small proteins. This demonstrates, that this yet widely undiscovered part of the proteome undergoes the same modifications as classical proteins which is an essential step for future understanding of the biological functions of these molecules.

Genome-Wide Transcription Start Site Mapping and Promoter Assignments to a Sigma Factor in the Human Enteropathogen Clostridioides difficile

The emerging human enteropathogen Clostridioides difficile is the main cause of diarrhea associated with antibiotherapy. Regulatory pathways underlying the adaptive responses remain understudied and the global view of C. difficile promoter structure is still missing. In the genome of C. difficile 630, 22 genes encoding sigma factors are present suggesting a complex pattern of transcription in this bacterium. We present here the first transcriptional map of the C. difficile genome resulting from the identification of transcriptional start sites (TSS), promoter motifs and operon structures. By 5′-end RNA-seq approach, we mapped more than 1000 TSS upstream of genes. In addition to these primary TSS, this analysis revealed complex structure of transcriptional units such as alternative and internal promoters, potential RNA processing events and 5′ untranslated regions. By following an in silico iterative strategy that used as an input previously published consensus sequences and transcriptomic analysis, we identified candidate promoters upstream of most of protein-coding and non-coding RNAs genes. This strategy also led to refine consensus sequences of promoters recognized by major sigma factors of C. difficile. Detailed analysis focuses on the transcription in the pathogenicity locus and regulatory genes, as well as regulons of transition phase and sporulation sigma factors as important components of C. difficile regulatory network governing toxin gene expression and spore formation. Among the still uncharacterized regulons of the major sigma factors of C. difficile, we defined the SigL regulon by combining transcriptome and in silico analyses. We showed that the SigL regulon is largely involved in amino-acid degradation, a metabolism crucial for C. difficile gut colonization. Finally, we combined our TSS mapping, in silico identification of promoters and RNA-seq data to improve gene annotation and to suggest operon organization in C. difficile. These data will considerably improve our knowledge of global regulatory circuits controlling gene expression in C. difficile and will serve as a useful rich resource for scientific community both for the detailed analysis of specific genes and systems biology studies.

The CARF Protein MM_0565 Affects Transcription of the Casposon-Encoded cas1-solo Gene in Methanosarcina mazei Gö1

Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR) loci are found in bacterial and archaeal genomes where they provide the molecular machinery for acquisition of immunity against foreign DNA. In addition to the cas genes fundamentally required for CRISPR activity, a second class of genes is associated with the CRISPR loci, of which many have no reported function in CRISPR-mediated immunity. Here, we characterize MM_0565 associated to the type I-B CRISPR-locus of Methanosarcina mazei Gö1. We show that purified MM_0565 composed of a CRISPR-Cas Associated Rossmann Fold (CARF) and a winged helix-turn-helix domain forms a dimer in solution; in vivo, the dimeric MM_0565 is strongly stabilized under high salt stress. While direct effects on CRISPR-Cas transcription were not detected by genetic approaches, specific binding of MM_0565 to the leader region of both CRISPR-Cas systems was observed by microscale thermophoresis and electromobility shift assays. Moreover, overexpression of MM_0565 strongly induced transcription of the cas1-solo gene located in the recently reported casposon, the gene product of which shows high similarity to classical Cas1 proteins. Based on our findings, and taking the absence of the expressed CRISPR locus-encoded Cas1 protein into account, we hypothesize that MM_0565 might modulate the activity of the CRISPR systems on different levels.

Adaptation of Methanosarcina barkeri 227 as acetate scavenger for succinate fermentation by Actinobacillus succinogenes

Acetate is the main by-product from microbial succinate production. In this study, we performed acetate removal by Methanosarcina barkeri 227 for succinate fermentation by Actinobacillus succinogenes 130Z. The acetoclastic methanogen M. barkeri requires similar environmental factors to A. succinogenes, and the conditions required for co-cultivation were optimized in this study: gas used for anaerobicization, strain adaptation, medium composition, pH adjustment, and inoculation time points. M. barkeri 227 was adapted to acetate for 150 days, which accelerated the acetate consumption to 9-fold (from 190 to 1726 mmol gDW^-1 day^-1). In the acetate-adapted strain, there was a noticeable increase in transcription of genes required for acetoclastic pathway-satP (acetate transporter), ackA (acetate kinase), cdhA (carbon monoxide dehydrogenase/acetyl-CoA synthase complex), and mtrH (methyl-H₄STP:CoM methyltransferase), which was not induced before the adaptation process. The activities of two energy-consuming steps in the pathway-acetate uptake and acetate kinase-increased about 3-fold. This acetate-adapted M. barkeri could be successfully applied to succinate fermentation culture of A. succinogenes, but only after pH adjustment following completion of fermentation. This study suggests the utility of M. barkeri as an acetate scavenger during fermentation for further steps towards genetic and process engineering.

Responses of Methanosarcina barkeri to acetate stress

BACKGROUND

Anaerobic digestion of easily degradable biowaste can lead to the accumulation of volatile fatty acids, which will cause environmental stress to the sensitive methanogens consequently. The metabolic characteristics of methanogens under acetate stress can affect the overall performance of mixed consortia. Nevertheless, there exist huge gaps in understanding the responses of the dominant methanogens to the stress, e.g., Methanosarcinaceae. Such methanogens are resistant to environmental deterioration and able to utilize multiple carbon sources. In this study, transcriptomic and proteomic analyses were conducted to explore the responses of Methanosarcina barkeri strain MS at different acetate concentrations of 10, 25, and 50 mM.

RESULTS

The trend of OD600 and the regulation of the specific genes in 50 mM acetate, indicated that high concentration of acetate promoted the acclimation of M. barkeri to acetate stress. Acetate stress hindered the regulation of quorum sensing and thereby eliminated the advantages of cell aggregation, which was beneficial to resist stress. Under acetate stress, M. barkeri allocated more resources to enhance the uptake of iron to maintain the integrities of electron-transport chains and other essential biological processes. Comparing with the initial stages of different acetate concentrations, most of the genes participating in acetoclastic methanogenesis did not show significantly different expressions except hdrB1C1, an electron-bifurcating heterodisulfide reductase participating in energy conversion and improving thermodynamic efficiency. Meanwhile, vnfDGHK and nifDHK participating in nitrogen fixation pathway were upregulated.

CONCLUSION

In this work, transcriptomic and proteomic analyses are combined to reveal the responses of M. barkeri to acetate stress in terms of central metabolic pathways, which provides basic clues for exploring the responses of other specific methanogens under high organics load. Moreover, the results can also be used to gain insights into the complex interactions and geochemical cycles among natural or engineered populations. Furthermore, these findings also provide the potential for designing effective and robust anaerobic digesters with high organic loads.

Gene Expression of Haloferax volcanii on Intermediate and Abundant Sources of Fixed Nitrogen

Haloferax volcanii, a well-developed model archaeon for genomic, transcriptomic, and proteomic analyses, can grow on a defined medium of abundant and intermediate levels of fixed nitrogen. Here we report a global profiling of gene expression of H. volcanii grown on ammonium as an abundant source of fixed nitrogen compared to l-alanine, the latter of which exemplifies an intermediate source of nitrogen that can be obtained from dead cells in natural habitats. By comparing the two growth conditions, 30 genes were found to be differentially expressed, including 16 genes associated with amino acid metabolism and transport. The gene expression profiles contributed to mapping ammonium and l-alanine usage with respect to transporters and metabolic pathways. In addition, conserved DNA motifs were identified in the putative promoter regions and transcription factors were found to be in synteny with the differentially expressed genes, leading us to propose regulons of transcriptionally co-regulated operons. This study provides insight to how H. volcanii responds to and utilizes intermediate vs. abundant sources of fixed nitrogen for growth, with implications for conserved functions in related halophilic archaea.

Global Transcriptional Programs in Archaea Share Features with the Eukaryotic Environmental Stress Response

The environmental stress response (ESR), a global transcriptional program originally identified in yeast, is characterized by a rapid and transient transcriptional response composed of large, oppositely regulated gene clusters. Genes induced during the ESR encode core components of stress tolerance, macromolecular repair, and maintenance of homeostasis. In this review, we investigate the possibility for conservation of the ESR across the eukaryotic and archaeal domains of life. We first re-analyze existing transcriptomics data sets to illustrate that a similar transcriptional response is identifiable in Halobacterium salinarum, an archaeal model organism. To substantiate the archaeal ESR, we calculated gene-by-gene correlations, gene function enrichment, and comparison of temporal dynamics. We note reported examples of variation in the ESR across fungi, then synthesize high-level trends present in expression data of other archaeal species. In particular, we emphasize the need for additional high-throughput time series expression data to further characterize stress-responsive transcriptional programs in the Archaea. Together, this review explores an open question regarding features of global transcriptional stress response programs shared across domains of life.

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 regulatory RNA is involved in RNA duplex formation and biofilm regulation in Sulfolobus acidocaldarius

Non-coding RNAs (ncRNA) are involved in essential biological processes in all three domains of life. The regulatory potential of ncRNAs in Archaea is, however, not fully explored. In this study, RNA-seq analyses identified a set of 29 ncRNA transcripts in the hyperthermophilic archaeon Sulfolobus acidocaldarius that were differentially expressed in response to biofilm formation. The most abundant ncRNA of this set was found to be resistant to RNase R treatment (RNase R resistant RNA, RrrR(+)) due to duplex formation with a reverse complementary RNA (RrrR(−)). The deletion of the RrrR(+) gene resulted in significantly impaired biofilm formation, while its overproduction increased biofilm yield. RrrR(+) was found to act as an antisense RNA against the mRNA of a hypothetical membrane protein. The RrrR(+) transcript was shown to be stabilized by the presence of the RrrR(−) strand in S. acidocaldarius cell extracts. The accumulation of these RrrR duplexes correlates with an apparent absence of dsRNA degrading RNase III domains in archaeal proteins.

Next Generation DNA-Seq and Differential RNA-Seq Allow Re-annotation of the Pyrococcus furiosus DSM 3638 Genome and Provide Insights Into Archaeal Antisense Transcription

Pyrococcus furiosus DSM 3638 is a model organism for hyperthermophilic archaea with an optimal growth temperature near 100°C. The genome was sequenced about 18 years ago. However, some publications suggest that in contrast to other Pyrococcus species, the genome of P. furiosus DSM 3638 is prone to genomic rearrangements. Therefore, we re-sequenced the genome using third generation sequencing techniques. The new de novo assembled genome is 1,889,914 bp in size and exhibits high sequence identity to the published sequence. However, two major deviations were detected: (1) The genome is 18,342 bp smaller than the NCBI reference genome due to a recently described deletion. (2) The region between PF0349 and PF0388 is inverted most likely due an assembly problem for the original sequence. In addition, numerous minor variations, ranging from single nucleotide exchanges, deletions or insertions were identified. The total number of insertion sequence (IS) elements is also reduced from 30 to 24 in the new sequence. Re-sequencing of a 2-year-old “lab culture” using Nanopore sequencing confirmed the overall stability of the P. furiosus DSM 3638 genome even under normal lab conditions without taking any special care. To improve genome annotation, the updated DNA sequence was combined with an RNA sequencing approach. Here, RNAs from eight different growth conditions were pooled to increase the number of detected transcripts. Furthermore, a differential RNA-Seq approach was employed for the identification of transcription start sites (TSSs). In total, 2515 TSSs were detected and classified into 834 primary (pTSS), 797 antisense (aTSS), 739 internal and 145 secondary TSSs. Our analysis of the upstream regions revealed a well conserved archaeal promoter structure. Interrogation of the distances between pTSSs and aTSSs revealed a significant number of antisense transcripts, which are a result of bidirectional transcription from the same TATA box. This mechanism of antisense transcript production could be further confirmed by in vitro transcription experiments. We assume that bidirectional transcription gives rise to non-functional antisense RNAs and that this is a widespread phenomenon in archaea due to the architecture of the TATA element and the symmetric structure of the TATA-binding protein.

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.

Long-term effect of biochar amendment on the biodegradation of petroleum hydrocarbons in soil microbial fuel cells

Biochar is extensively applied in amendment of contaminated soils. However, the effect of biochar on the biodegradation of petroleum hydrocarbons and electricity generation in soil microbial fuel cells (MFCs) remains unclear. Here, three biochars respectively derived from poultry (chicken manure, CB), agriculture (wheat straw, SB) and forestry industries (wood sawdust, WB) were investigated after 223 days of amendment. Consequently, high removal for alkanes was in CB with the mineral nutrition and phosphorus while aromatics were in SB with the most N content and the highest molecular polarity. The lowest removal efficiency of total petroleum hydrocarbons was observed in WB with the highest surface area, whereas the most charge was obtained. The different performance of soil MFCs was due to physicochemical properties of biochar and colonized microbial communities of bacteria and archaea. The abundance of Actinotalea increased by 144-263% in SB and CB while that of Desulfatitalea distinctly increased in WB. Meanwhile, species from Methanosarcina, Methanoculleus, Halovivax and Natronorubrum exerted probably a methanogenic degrading role. This study revealed that the degrader, azotobacter and electricigens exhibited a close relationship in order to degrade hydrocarbons and generate electricity in soil bioelectrochemical remediation systems.

Small RNAs Involved in Regulation of Nitrogen Metabolism

Global (metabolic) regulatory networks allow microorganisms to survive periods of nitrogen starvation or general nutrient stress. Uptake and utilization of various nitrogen sources are thus commonly tightly regulated in Prokarya (Bacteria and Archaea) in response to available nitrogen sources. Those well-studied regulations occur mainly at the transcriptional and posttranslational level. Surprisingly, and in contrast to their involvement in most other stress responses, small RNAs (sRNAs) involved in the response to environmental nitrogen fluctuations are only rarely reported. In addition to sRNAs indirectly affecting nitrogen metabolism, only recently it was demonstrated that three sRNAs were directly involved in regulation of nitrogen metabolism in response to changes in available nitrogen sources. All three trans-acting sRNAs are under direct transcriptional control of global nitrogen regulators and affect expression of components of nitrogen metabolism (glutamine synthetase, nitrogenase, and PII-like proteins) by either masking the ribosome binding site and thus inhibiting translation initiation or stabilizing the respective target mRNAs. Most likely, there are many more sRNAs and other types of noncoding RNAs, e.g., riboswitches, involved in the regulation of nitrogen metabolism in Prokarya that remain to be uncovered. The present review summarizes the current knowledge on sRNAs involved in nitrogen metabolism and their biological functions and targets.

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.

Unique Archaeal Small RNAs

Advances in genome-wide sequence technologies allow for detailed insights into the complexity of RNA landscapes of organisms from all three domains of life. Recent analyses of archaeal transcriptomes identified interaction and regulation networks of noncoding RNAs in this understudied domain. Here, we review current knowledge of small, noncoding RNAs with important functions for the archaeal lifestyle, which often requires adaptation to extreme environments. One focus is RNA metabolism at elevated temperatures in hyperthermophilic archaea, which reveals elevated amounts of RNA-guided RNA modification and virus defense strategies. Genome rearrangement events result in unique fragmentation patterns of noncoding RNA genes that require elaborate maturation pathways to yield functional transcripts. RNA-binding proteins, e.g., L7Ae and LSm, are important for many posttranscriptional control functions of RNA molecules in archaeal cells. We also discuss recent insights into the regulatory potential of their noncoding RNA partners.

Cross-cleavage activity of Cas6b in crRNA processing of two different CRISPR-Cas systems in Methanosarcina mazei Gö1

The clustered regularly interspaced short palindromic repeat (CRISPR) system is a prokaryotic adaptive defense system against foreign nucleic acids. In the methanoarchaeon Methanosarcina mazei Gö1, two types of CRISPR-Cas systems are present (type I-B and type III-C). Both loci encode a Cas6 endonuclease, Cas6b-IB and Cas6b-IIIC, typically responsible for maturation of functional short CRISPR RNAs (crRNAs). To evaluate potential cross cleavage activity, we biochemically characterized both Cas6b proteins regarding their crRNA binding behavior and their ability to process pre-crRNA from the respective CRISPR array in vivo. Maturation of crRNA was studied in the respective single deletion mutants by northern blot and RNA-Seq analysis demonstrating that in vivo primarily Cas6b-IB is responsible for crRNA processing of both CRISPR arrays. Tentative protein level evidence for the translation of both Cas6b proteins under standard growth conditions was detected, arguing for different activities or a potential non-redundant role of Cas6b-IIIC within the cell. Conservation of both Cas6 endonucleases was observed in several other M. mazei isolates, though a wide variety was displayed. In general, repeat and leader sequence conservation revealed a close correlation in the M. mazei strains. The repeat sequences from both CRISPR arrays from M. mazei Gö1 contain the same sequence motif with differences only in two nucleotides. These data stand in contrast to all other analyzed M. mazei isolates, which have at least one additional CRISPR array with repeats belonging to another sequence motif. This conforms to the finding that Cas6b-IB is the crucial and functional endonuclease in M. mazei Gö1. Abbreviations: sRNA: small RNA; crRNA: CRISPR RNA; pre-crRNAs: Precursor CRISPR RNA; CRISPR: clustered regularly interspaced short palindromic repeats; Cas: CRISPR associated; nt: nucleotide; RNP: ribonucleoprotein; RBS: ribosome binding site.

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.

Transcriptional Regulation in Archaea: From Individual Genes to Global Regulatory Networks

Archaea are major contributors to biogeochemical cycles, possess unique metabolic capabilities, and resist extreme stress. To regulate the expression of genes encoding these unique programs, archaeal cells use gene regulatory networks (GRNs) composed of transcription factor proteins and their target genes. Recent developments in genetics, genomics, and computational methods used with archaeal model organisms have enabled the mapping and prediction of global GRN structures. Experimental tests of these predictions have revealed the dynamical function of GRNs in response to environmental variation. Here, we review recent progress made in this area, from investigating the mechanisms of transcriptional regulation of individual genes to small-scale subnetworks and genome-wide global networks. At each level, archaeal GRNs consist of a hybrid of bacterial, eukaryotic, and uniquely archaeal mechanisms. We discuss this theme from the perspective of the role of individual transcription factors in genome-wide regulation, how these proteins interact to compile GRN topological structures, and how these topologies lead to emergent, high-level GRN functions. We conclude by discussing how systems biology approaches are a fruitful avenue for addressing remaining challenges, such as discovering gene function and the evolution of GRNs.

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.

Antisense transcription in Pseudomonas aeruginosa

A large number of antisense transcripts have been detected in diverse microbial genomes and considerable effort has been devoted to elucidating the functional role of antisense transcription. In this study, we reanalysed extensive RNA sequencing data from the opportunistic pathogen Pseudomonas aeruginosa and found that the majority of genes have a propensity for antisense transcription. Although antisense transcripts were found in more than 80 % of the genes of the P. aeruginosa genome, the majority of sequencing reads were mapping sense and only a minority (<2 %) were mapping antisense to genes. Similarly to the sense expression levels, the antisense expression levels varied under different environmental conditions, with the sense and antisense expression levels often being inversely regulated and modulated by the activity of alternative sigma factors. Environment-modulated antisense transcription showed a bias towards being antisense to genes within regions of genomic plasticity and to those encoding small regulatory RNAs. In the future, the validation and functional characterization of antisense transcripts, and novel transcripts that are antisense to small regulatory RNAs in particular, have the potential to contribute to our understanding of the various levels of transcriptional regulation and its dynamics in the bacterial pathogen P. aeruginosa.

Transcriptional Landscape and Regulatory Roles of Small Noncoding RNAs in the Oxidative Stress Response of the Haloarchaeon Haloferax volcanii

Haloarchaea in their natural environment are exposed to hypersalinity, intense solar radiation, and desiccation, all of which generate high levels of oxidative stress. Previous work has shown that haloarchaea are an order of magnitude more resistant to oxidative stress than most mesophilic organisms. Despite this resistance, the pathways haloarchaea use to respond to oxidative stress damage are similar to those of nonresistant organisms, suggesting that regulatory processes might be key to their robustness. Recently, small regulatory noncoding RNAs (sRNAs) were discovered in Archaea under a variety of environmental conditions. We report here the transcriptional landscape and functional roles of sRNAs in the regulation of the oxidative stress response of the model haloarchaeon Haloferax volcanii. Thousands of sRNAs, both intergenic and antisense, were discovered using strand-specific sRNA sequencing (sRNA-seq), comprising 25 to 30% of the total transcriptome under no-challenge and oxidative stress conditions, respectively. We identified hundreds of differentially expressed sRNAs in response to hydrogen peroxide-induced oxidative stress in H. volcanii. The targets of a group of antisense sRNAs decreased in expression when these sRNAs were upregulated, suggesting that sRNAs are potentially playing a negative regulatory role on mRNA targets at the transcript level. Target enrichment of these antisense sRNAs included mRNAs involved in transposon mobility, chemotaxis signaling, peptidase activity, and transcription factors.

IMPORTANCE While a substantial body of experimental work has been done to uncover the functions of small regulatory noncoding RNAs (sRNAs) in gene regulation in Bacteria and Eukarya, the functional roles of sRNAs in Archaea are still poorly understood. This study is the first to establish the regulatory effects of sRNAs on mRNAs during the oxidative stress response in the haloarchaeon Haloferax volcanii. Our work demonstrates that common principles for the response to a major cellular stress exist across the 3 domains of life while uncovering pathways that might be specific to the Archaea. This work also underscores the relevance of sRNAs in adaptation to extreme environmental conditions.

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.

Small RNAs of Haloferax mediterranei: Identification and Potential Involvement in Nitrogen Metabolism

Small RNAs have been studied in detail in domains Bacteria and Eukarya but, in the case of the domain Archaea, the knowledge is scarce and the physiological function of these small RNAs (sRNAs) is still uncertain. To extend the knowledge of sRNAs in the domain Archaea and their possible role in the regulation of the nitrogen assimilation metabolism in haloarchaea, Haloferax mediterranei has been used as a model microorganism. The bioinformatic approach has allowed for the prediction of 295 putative sRNAs genes in the genome of H. mediterranei, 88 of which have been verified by means of RNA-Sequencing (RNA-Seq). The secondary structure of these sRNAs and their possible targets have been identified. Curiously, some of them present as possible target genes relating to nitrogen assimilation, such as glutamate dehydrogenase and the nitrogen regulatory PII protein. Analysis of RNA-Seq data has also revealed differences in the expression pattern of 16 sRNAs according to the nitrogen source. Consequently, RNomic and bioinformatic approaches used in this work have allowed for the identification of new sRNAs in H. mediterranei, some of which show different expression patterns depending on the nitrogen source. This suggests that these sRNAs could be involved in the regulation of nitrogen assimilation and can constitute an important gene regulatory network.

Nutrient depletion and TOR inhibition induce 18S and 25S ribosomal RNAs resistant to a 5'-phosphate-dependent exonuclease in Candida albicans and other yeasts

Background

Messenger RNA (mRNA) represents a small percentage of RNAs in a cell, with ribosomal RNA (rRNA) making up the bulk of it. To isolate mRNA from eukaryotes, typically poly-A selection is carried out. Recently, a 5´-phosphate-dependent, 5´→3´ processive exonuclease called Terminator has become available. It will digest only RNA that has a 5´-monophosphate end and therefore it is very useful to eliminate most of rRNAs in cell.

Results

We have found that in the pathogenic yeast Candida albicans, while 18S and 25S components isolated from yeast in robust growth phase are easily eliminated by Terminator, those isolated from cells in the nutritionally diminished stationary phase, become resistant to digestion by this enzyme. Additional digestions with alkaline phosphatase, tobacco pyrophosphatase combined with Terminator point toward the 5′-prime end of 18S and 25S as the source of this resistance. Inhibition of TOR by rapamycin also induces resistance by these molecules. We also find that these molecules are incorporated into the ribosome and are not just produced incidentally. Finally, we show that three other yeasts show the same behavior.

Conclusions

Digestion of RNA by Terminator has revealed 18S and 25S rRNA molecules different from the accepted processed ones seen in ribosome generation. The reason for these molecules and the underlying mechanism for their formation is unknown. The preservation of this behavior across these yeasts suggests a useful biological role for it, worthy of further inquiry.

sRNA41 affects ribosome binding sites within polycistronic mRNAs in Methanosarcina mazei Gö1

Several noncoding RNAs potentially involved in nitrogen (N)-regulation have been detected in Methanosarcina mazei, however, targets have been identified only for one of them. Here, we report on the function of sRNA₄₁ , highly expressed under N-sufficiency. Comprising 120 nucleotides, sRNA₄₁ shows high sequence and structural conservation within draft genomes of numerous Methanosarcina species. In silico target prediction revealed several potential targets, including genes of two homologous operons encoding for acetyl-CoA-decarbonylase/synthase complexes (ACDS) representing highly probable target candidates. A highly conserved single stranded region of sRNA₄₁ was predicted to mask six independent ribosome binding sites of these two polycistronic mRNAs and was verified in vitro by microscale thermophoresis. Proteome analysis of the respective sRNA₄₁ -deletion mutant showed increased protein expression of both ACDS complexes in the absence of sRNA₄₁ , whereas no effect on transcript levels was detected, arguing for sRNA₄₁ -mediated post-transcriptional fine-tuning of ACDS expression. We hypothesize that the physiological advantage of downregulating sRNA₄₁ under N-limiting conditions is the resulting increase of ACDS protein levels. This provides sufficient amounts of amino acids for nitrogenase synthesis as well as reducing equivalents and energy for N₂ -fixation, thus linking the carbon and N-metabolism.

RNA Sequencing and Analysis in Microorganisms for Metabolic Network Reconstruction

There is a strict interplay between metabolic networks and transcriptional regulation in bacteria; indeed, the transcriptome regulation, affecting the expression of large gene sets, can be used to predict the likely "on" or "off" state of metabolic genes as a function of environmental factors. Up to date, many bacterial transcriptomes have been studied by RNAseq, hundreds of experiments have been performed, and Giga bases of sequences have been produced. All this transcriptional information could potentially be integrated into metabolic networks in order to obtain a more comprehensive view of their regulation and to increase their prediction power.To get high-quality transcriptomic data, to be integrated into metabolic networks, it is paramount to clearly know how to produce highly informative RNA sequencing libraries and how to manage RNA sequencing data.In this chapter, we will get across the main steps of an RNAseq experiment: from removal of ribosomal RNAs, to strand-specific library preparation, till data analysis and integration. We will try to share our experience and know-how, to give you a precise protocol to follow, and some useful recommendations or tips and tricks to adopt in order to go straightforward toward a successful RNAseq experiment.