Global transcriptional analysis of Methanosarcina mazei strain Gö1 under different nitrogen availabilities

Quantitative analysis of amino acid excretion by Methanothermobacter marburgensis under N2-fixing conditions

Methanogenic archaea (methanogens) possess fascinating metabolic characteristics, such as the ability to fix molecular nitrogen (N2). Methanogens are of biotechnological importance due to the ability to produce methane (CH4) from molecular hydrogen (H2) and carbon dioxide (CO2) and to excrete proteinogenic amino acids. This study focuses on analyzing the link between biological methanogenesis and amino acid excretion under N2-fixing conditions. Among five hydrogenotrophic, autotrophic methanogens, Methanothermobacter marburgensis was prioritized and further cultivated in closed batch cultivation mode under N2-fixing conditions. M. marburgensis was grown on chemically defined minimal medium with different concentrations of ammonium in a H2/CO2/N2 atmosphere. This enabled the quantification of ammonia uptake, N2-fixation, amino acid excretion and the conversion of H2/CO2 to CH4. To quantify N2-fixation rates in a mass balance setting a novel method has been established. The method utilizes the pressure drop below a certain threshold pressure in closed batch cultivation mode - the threshold pressure for N2-fixation (THpN2fix). Using the THpN2fix method, volumetric N2-fixation rates of M. marburgensis as high as 0.91 mmol L-1 h-1 were determined. Excretion of amino acids was found with highest detected values of glutamic acid, alanine, glycine and asparagine. The highest total amino acid excretion of 7.5 µmol L-1 h-1 was detected with H2/CO2/N2 at an ammonium concentration of 40 mmol L-1. This study sheds light on the link between methanogenesis, biological N2-fixation, and proteinogenic amino acid excretion. The concomitant production of amino acids and CH4 could become of biotechnological relevance in an integrated approach coupling biomethanation and N2-fixation in a biorefinery concept.

Uncovering the small proteome of Methanosarcina mazei using Ribo-seq and peptidomics under different nitrogen conditions

The mesophilic methanogenic archaeal model organism Methanosarcina mazei strain Gö1 is crucial for climate and environmental research due to its ability to produce methane. Here, we establish a Ribo-seq protocol for M. mazei strain Gö1 under two growth conditions (nitrogen sufficiency and limitation). The translation of 93 previously annotated and 314 unannotated small ORFs, coding for proteins ≤ 70 amino acids, is predicted with high confidence based on Ribo-seq data. LC-MS analysis validates the translation for 62 annotated small ORFs and 26 unannotated small ORFs. Epitope tagging followed by immunoblotting analysis confirms the translation of 13 out of 16 selected unannotated small ORFs. A comprehensive differential transcription and translation analysis reveals that 29 of 314 unannotated small ORFs are differentially regulated in response to nitrogen availability at the transcriptional and 49 at the translational level. A high number of reported small RNAs are emerging as dual-function RNAs, including sRNA154, the central regulatory small RNA of nitrogen metabolism. Several unannotated small ORFs are conserved in Methanosarcina species and overproducing several (small ORF encoded) small proteins suggests key physiological functions. Overall, the comprehensive analysis opens an avenue to elucidate the function(s) of multitudinous small proteins and dual-function RNAs in M. mazei.

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.

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.

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.

Archaeal Connectase is a specific and efficient protein ligase related to proteasome β subunits

Sequence-specific protein ligations are widely used to produce customized proteins "on demand." Such chimeric, immobilized, fluorophore-conjugated or segmentally labeled proteins are generated using a range of chemical, (split) intein, split domain, or enzymatic methods. Where short ligation motifs and good chemoselectivity are required, ligase enzymes are often chosen, although they have a number of disadvantages, for example poor catalytic efficiency, low substrate specificity, and side reactions. Here, we describe a sequence-specific protein ligase with more favorable characteristics. This ligase, Connectase, is a monomeric homolog of 20S proteasome subunits in methanogenic archaea. In pulldown experiments with Methanosarcina mazei cell extract, we identify a physiological substrate in methyltransferase A (MtrA), a key enzyme of archaeal methanogenesis. Using microscale thermophoresis and X-ray crystallography, we show that only a short sequence of about 20 residues derived from MtrA and containing a highly conserved KDPGA motif is required for this high-affinity interaction. Finally, in quantitative activity assays, we demonstrate that this recognition tag can be repurposed to allow the ligation of two unrelated proteins. Connectase catalyzes such ligations at substantially higher rates, with higher yields, but without detectable side reactions when compared with a reference enzyme. It thus presents an attractive tool for the development of new methods, for example in the preparation of selectively labeled proteins for NMR, the covalent and geometrically defined attachment of proteins on surfaces for cryo-electron microscopy, or the generation of multispecific antibodies.

Comparative genomic analysis reveals starvation survival systems in Methanothermobacter thermautotrophicus ΔH

Methanothermobacter thermautotrophicus ΔH (MTH) is a thermophilic hydrogenotrophic methanogenic archaeon capable of reducing CO₂ with H₂ to produce methane gas. It is the potential candidate in the biomethanation of CO₂ and CO in anaerobic reactors and biogas upgrading process. However, systematic studies addressing its genome conservation and function remain scant in this genome. In this study, we have evaluated its evolutionary resemblance and metabolic discrepancy, particularly in starvation survival systems by comparing the genomic contexts with Methanothermobacter marburgensis str. Marburg (MMG) and Methanobacterium formicicum DSM 1535 (MFO). The phylogenomic analysis of this study indicated that there was a strong phylogenomic signal among MTH, MMG, and MFO in the whole-genome tree. DNA replication machinery was conserved in the MTH genome and might have evolved at different evolution rates. Genome synteny analysis observed collinearity of either gene orders or gene families has to be maintained with syntenic blocks located in the syntenic out-paralogs. A genome-wide metabolic analysis identified some unique putative metabolic subsystems in MTH, which are proposed to determine its growth characteristics in diverse environments. MTH genome comprised of 93 unique genes-coding for starvation survival and stress-response proteins. These proteins confer its adaptation to nutritional deprivation and other abiotic stresses. MTH has a typical system to withstand its growth and cell viability during stable operation and recovery after prolonged starvation. Thus, the present work will provide an insight to improve the genome refinement and metabolic reconstruction in parallel to other closely related species.

The potential for polyphosphate metabolism in Archaea and anaerobic polyphosphate formation in Methanosarcina mazei

Inorganic polyphosphate (polyP) is ubiquitous across all forms of life, but the study of its metabolism has been mainly confined to bacteria and yeasts. Few reports detail the presence and accumulation of polyP in Archaea, and little information is available on its functions and regulation. Here, we report that homologs of bacterial polyP metabolism proteins are present across the major taxa in the Archaea, suggesting that archaeal populations may have a greater contribution to global phosphorus cycling than has previously been recognised. We also demonstrate that polyP accumulation can be induced under strictly anaerobic conditions, in response to changes in phosphate (Pi) availability, i.e. Pi starvation, followed by incubation in Pi replete media (overplus), in cells of the methanogenic archaeon Methanosarcina mazei. Pi-starved M. mazei cells increased transcript abundance of the alkaline phosphatase (phoA) gene and of the high-affinity phosphate transport (pstSCAB-phoU) operon: no increase in polyphosphate kinase 1 (ppk1) transcript abundance was observed. Subsequent incubation of Pi-starved M. mazei cells under Pi replete conditions, led to a 237% increase in intracellular polyphosphate content and a > 5.7-fold increase in ppk1 gene transcripts. Ppk1 expression in M. mazei thus appears not to be under classical phosphate starvation control.

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.

Methylamine-specific methyltransferase paralogs in Methanosarcina are functionally distinct despite frequent gene conversion

Sequenced archaeal genomes are mostly smaller and more streamlined than typical bacterial genomes; however, members of the Methanosarcina genus within the Euryarchaeaota are a significant exception, with M. acetivorans being the largest archaeal genome (5.8 Mbp) sequenced thus far. This finding is partially explained by extensive gene duplication within Methanosarcina spp. Significantly, the evolutionary pressures leading to gene duplication and subsequent genome expansion have not been well investigated, especially with respect to biological methane production (methanogenesis), which is the key biological trait of these environmentally important organisms. In this study, we address this question by specifically probing the functional evolution of two methylamine-specific methyltransferase paralogs in members of the Methanosarcina genus. Using the genetically tractable strain, M. acetivorans, we first show that the two paralogs have distinct cellular functions: one being required for methanogenesis from methylamine, the other for use of methylamine as a nitrogen source. Subsequently, through comparative sequence analyses, we show that functional divergence of paralogs is primarily mediated by divergent evolution of the 5' regulatory region, despite frequent gene conversion within the coding sequence. This unique evolutionary paradigm for functional divergence of genes post-duplication underscores a divergent role for an enzyme singularly associated with methanogenic metabolism in other aspects of cell physiology.

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.

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.

Methanosarcina Spherical Virus, a Novel Archaeal Lytic Virus Targeting Methanosarcina Strains

A novel archaeal lytic virus targeting species of the genus Methanosarcina was isolated using Methanosarcina mazei strain Gö1 as the host. Due to its spherical morphology, the virus was designated Methanosarcina spherical virus (MetSV). Molecular analysis demonstrated that MetSV contains double-stranded linear DNA with a genome size of 10,567 bp containing 22 open reading frames (ORFs), all oriented in the same direction. Functions were predicted for some of these ORFs, i.e., such as DNA polymerase, ATPase, and DNA-binding protein as well as envelope (structural) protein. MetSV-derived spacers in CRISPR loci were detected in several published Methanosarcina draft genomes using bioinformatic tools, revealing a potential protospacer-adjacent motif (PAM) motif (TTA/T). Transcription and expression of several predicted viral ORFs were validated by reverse transcription-PCR (RT-PCR), PAGE analysis, and liquid chromatography-mass spectrometry (LC-MS)-based proteomics. Analysis of core lipids by atmospheric pressure chemical ionization (APCI) mass spectrometry showed that MetSV and Methanosarcina mazei both contain archaeol and glycerol dialkyl glycerol tetraether without a cyclopentane moiety (GDGT-0). The MetSV host range is limited to Methanosarcina strains growing as single cells (M. mazei, Methanosarcina barkeri and Methanosarcina soligelidi). In contrast, strains growing as sarcina-like aggregates were apparently protected from infection. Heterogeneity related to morphology phases in M. mazei cultures allowed acquisition of resistance to MetSV after challenge by growing cultures as sarcina-like aggregates. CRISPR/Cas-mediated resistance was excluded since neither of the two CRISPR arrays showed MetSV-derived spacer acquisition. Based on these findings, we propose that changing the morphology from single cells to sarcina-like aggregates upon rearrangement of the envelope structure prevents infection and subsequent lysis by MetSV.IMPORTANCE Methanoarchaea are among the most abundant organisms on the planet since they are present in high numbers in major anaerobic environments. They convert various carbon sources, e.g., acetate, methylamines, or methanol, to methane and carbon dioxide; thus, they have a significant impact on the emission of major greenhouse gases. Today, very little is known about viruses specifically infecting methanoarchaea that most probably impact the abundance of methanoarchaea in microbial consortia. Here, we characterize the first identified Methanosarcina-infecting virus (MetSV) and show a mechanism for acquiring resistance against MetSV. Based on our results, we propose that growth as sarcina-like aggregates prevents infection and subsequent lysis. These findings allow new insights into the virus-host relationship in methanogenic community structures, their dynamics, and their phase heterogeneity. Moreover, the availability of a specific virus provides new possibilities to deepen our knowledge of the defense mechanisms of potential hosts and offers tools for genetic manipulation.

sRNA154 a newly identified regulator of nitrogen fixation in Methanosarcina mazei strain Gö1

Trans-encoded sRNA₁₅₄ is exclusively expressed under nitrogen (N)-deficiency in Methanosarcina mazei strain Gö1. The sRNA₁₅₄ deletion strain showed a significant decrease in growth under N-limitation, pointing toward a regulatory role of sRNA₁₅₄ in N-metabolism. Aiming to elucidate its regulatory function we characterized sRNA₁₅₄ by means of biochemical and genetic approaches. 24 homologs of sRNA₁₅₄ were identified in recently reported draft genomes of Methanosarcina strains, demonstrating high conservation in sequence and predicted secondary structure with two highly conserved single stranded loops. Transcriptome studies of sRNA₁₅₄ deletion mutants by an RNA-seq approach uncovered nifH- and nrpA-mRNA, encoding the α-subunit of nitrogenase and the transcriptional activator of the nitrogen fixation (nif)-operon, as potential targets besides other components of the N-metabolism. Furthermore, results obtained from stability, complementation and western blot analysis, as well as in silico target predictions combined with electrophoretic mobility shift-assays_, argue for a stabilizing effect of sRNA₁₅₄ on the polycistronic nif-mRNA and nrpA-mRNA by binding with both loops. Further identified N-related targets were studied, which demonstrates that translation initiation of glnA₂-mRNA, encoding glutamine synthetase2, appears to be affected by sRNA₁₅₄ masking the ribosome binding site, whereas glnA₁-mRNA appears to be stabilized by sRNA₁₅₄. Overall, we propose that sRNA₁₅₄ has a crucial regulatory role in N-metabolism in M. mazei by stabilizing the polycistronic mRNA encoding nitrogenase and glnA₁-mRNA, as well as allowing a feed forward regulation of nif-gene expression by stabilizing nrpA-mRNA. Consequently, sRNA₁₅₄ represents the first archaeal sRNA, for which a positive posttranscriptional regulation is demonstrated as well as inhibition of translation initiation.

Horizontal gene transfer drives the evolution of Rh50 permeases in prokaryotes

Background

Rh50 proteins belong to the family of ammonia permeases together with their Amt/MEP homologs. Ammonia permeases increase the permeability of NH₃/NH₄⁺ across cell membranes and are believed to be involved in excretion of toxic ammonia and in the maintenance of pH homeostasis. RH50 genes are widespread in eukaryotes but absent in land plants and fungi, and remarkably rare in prokaryotes. The evolutionary history of RH50 genes in prokaryotes is just beginning to be unveiled.

Results

Here, a molecular phylogenetic approach suggests horizontal gene transfer (HGT) as a primary force driving the evolution and spread of RH50 among prokaryotes. In addition, the taxonomic distribution of the RH50 gene among prokaryotes turned out to be very narrow; a single-copy RH50 is present in the genome of only a small proportion of Bacteria, and, first evidence to date, in only three methanogens among Euryarchaea. The coexistence of RH50 and AMT in prokaryotes seems also a rare event. Finally, phylogenetic analyses were used to reconstruct the HGT network along which prokaryotic RH50 evolution has taken place.

Conclusions

The eukaryotic or bacterial “origin” of the RH50 gene remains unsolved. The RH50 prokaryotic HGT network suggests a preferential directionality of transfer from aerobic to anaerobic organisms. The observed HGT events between archaeal methanogens, anaerobic and aerobic ammonia-oxidizing bacteria suggest that syntrophic relationships play a major role in the structuring of the network, and point to oxygen minimum zones as an ecological niche that might be of crucial importance for HGT-driven evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0850-6) contains supplementary material, which is available to authorized users.

Horizontal gene transfer and genome evolution in Methanosarcina

Background

Genomes of Methanosarcina spp. are among the largest archaeal genomes. One suggested reason for that is massive horizontal gene transfer (HGT) from bacteria. Genes of bacterial origin may be involved in the central metabolism and solute transport, in particular sugar synthesis, sulfur metabolism, phosphate metabolism, DNA repair, transport of small molecules etc. Horizontally transferred (HT) genes are considered to play the key role in the ability of Methanosarcina spp. to inhabit diverse environments. At the moment, genomes of three Methanosarcina spp. have been sequenced, and while these genomes vary in length and number of protein-coding genes, they all have been shown to accumulate HT genes. However, previous estimates had been made when fewer archaeal genomes were known. Moreover, several Methanosarcinaceae genomes from other genera have been sequenced recently. Here, we revise the census of genes of bacterial origin in Methanosarcinaceae.

Results

About 5 % of Methanosarcina genes have been shown to be horizontally transferred from various bacterial groups to the last common ancestor either of Methanosarcinaceae, or Methanosarcina, or later in the evolution. Simulation of the composition of the NCBI protein non-redundant database for different years demonstrates that the estimates of the HGT rate have decreased drastically since 2002, the year of publication of the first Methanosarcina genome.

The phylogenetic distribution of HT gene donors is non-uniform. Most HT genes were transferred from Firmicutes and Proteobacteria, while no HGT events from Actinobacteria to the common ancestor of Methanosarcinaceae were found. About 50 % of HT genes are involved in metabolism. Horizontal transfer of transcription factors is not common, while 46 % of horizontally transferred genes have demonstrated differential expression in a variety of conditions. HGT of complete operons is relatively infrequent and half of HT genes do not belong to operons.

Conclusions

While genes of bacterial origin are still more frequent in Methanosarcinaceae than in other Archaea, most HGT events described earlier as Methanosarcina-specific seem to have occurred before the divergence of Methanosarcinaceae. Genes horizontally transferred from bacteria to archaea neither tend to be transferred with their regulators, nor in long operons.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0393-2) contains supplementary material, which is available to authorized users.

The transcriptional activator NrpA is crucial for inducing nitrogen fixation in Methanosarcina mazei Gö1 under nitrogen-limited conditions

With the aim of unraveling their potential involvement in the regulation of nitrogen metabolism in Methanosarcina mazei strain Gö1, we characterized five genes that are differentially transcribed in response to changing nitrogen availability and encoding putative transcriptional regulators. Study of the respective mutant strains under nitrogen-limited conditions revealed a growth delay for M. mazei MM0444::pac and MM1708::pac, and strongly reduced diazotrophic growth for MM0872::pac, whereas the absence of MM2441 or MM2525 did not affect growth behaviour. Transcriptome analyses further demonstrated that only MM1708 - encoding a CxxCG zinc finger protein - plays a regulatory role in nitrogen metabolism, most likely by specifically enhancing transcription of the N2 fixation (nif) operon under nitrogen-limited conditions. In agreement with this, a palindromic binding motif was predicted in silico in the nifH promoter region, nine nucleotides upstream of the BRE box, and confirmed to bind purified maltose-binding protein-MM1708 by electromobility shift assays. As MM1708 itself is under the control of the global nitrogen repressor NrpR, this adds a secondary level to the transcriptional regulation of the nif genes, and is most likely crucial for maximal nif induction under nitrogen-limited conditions. This is in accordance with the finding that protein expression of NifH is highly reduced in the absence of MM1708 under nitrogen-limited conditions. On the basis of our findings, we hypothesize that, in M. mazei, nitrogen fixation is controlled by a hierarchical network of two transcriptional regulators, the global nitrogen repressor NrpR, and the newly identified activator NrpA (MM1708), thereby providing tight control of N2 fixation.

Small regulatory RNAs in Archaea

Small regulatory RNAs (sRNAs) are universally distributed in all three domains of life, Archaea, Bacteria, and Eukaryotes. In bacteria, sRNAs typically function by binding near the translation start site of their target mRNAs and thereby inhibit or activate translation. In eukaryotes, miRNAs and siRNAs typically bind to the 3′-untranslated region (3′-UTR) of their target mRNAs and influence translation efficiency and/or mRNA stability. In archaea, sRNAs have been identified in all species investigated using bioinformatic approaches, RNomics, and RNA-Seq. Their size can vary significantly between less than 50 to more than 500 nucleotides. Differential expression of sRNA genes has been studied using northern blot analysis, microarrays, and RNA-Seq. In addition, biological functions have been unraveled by genetic approaches, i.e., by characterization of designed mutants. As in bacteria, it was revealed that archaeal sRNAs are involved in many biological processes, including metabolic regulation, adaptation to extreme conditions, stress responses, and even in regulation of morphology and cellular behavior. Recently, the first target mRNAs were identified in archaea, including one sRNA that binds to the 5′-region of two mRNAs in Methanosarcina mazei Gö1 and a few sRNAs that bind to 3′-UTRs in Sulfolobus solfataricus, three Pyrobaculum species, and Haloferax volcanii, indicating that archaeal sRNAs appear to be able to target both the 5′-UTR or the 3′-UTRs of their respective target mRNAs. In addition, archaea contain tRNA-derived fragments (tRFs), and one tRF has been identified as a major ribosome-binding sRNA in H. volcanii, which downregulates translation in response to stress. Besides regulatory sRNAs, archaea contain further classes of sRNAs, e.g., CRISPR RNAs (crRNAs) and snoRNAs.

Towards a computational model of a methane producing archaeum

Progress towards a complete model of the methanogenic archaeum Methanosarcina acetivorans is reported. We characterized size distribution of the cells using differential interference contrast microscopy, finding them to be ellipsoidal with mean length and width of 2.9 μm and 2.3 μm, respectively, when grown on methanol and 30% smaller when grown on acetate. We used the single molecule pull down (SiMPull) technique to measure average copy number of the Mcr complex and ribosomes. A kinetic model for the methanogenesis pathways based on biochemical studies and recent metabolic reconstructions for several related methanogens is presented. In this model, 26 reactions in the methanogenesis pathways are coupled to a cell mass production reaction that updates enzyme concentrations. RNA expression data (RNA-seq) measured for cell cultures grown on acetate and methanol is used to estimate relative protein production per mole of ATP consumed. The model captures the experimentally observed methane production rates for cells growing on methanol and is most sensitive to the number of methyl-coenzyme-M reductase (Mcr) and methyl-tetrahydromethanopterin:coenzyme-M methyltransferase (Mtr) proteins. A draft transcriptional regulation network based on known interactions is proposed which we intend to integrate with the kinetic model to allow dynamic regulation.

A functional approach to uncover the low-temperature adaptation strategies of the archaeon Methanosarcina barkeri

Low-temperature anaerobic digestion (LTAD) technology is underpinned by a diverse microbial community. The methanogenic archaea represent a key functional group in these consortia, undertaking CO2 reduction as well as acetate and methylated C1 metabolism with subsequent biogas (40 to 60% CH4 and 30 to 50% CO2) formation. However, the cold adaptation strategies, which allow methanogens to function efficiently in LTAD, remain unclear. Here, a pure-culture proteomic approach was employed to study the functional characteristics of Methanosarcina barkeri (optimum growth temperature, 37°C), which has been detected in LTAD bioreactors. Two experimental approaches were undertaken. The first approach aimed to characterize a low-temperature shock response (LTSR) of M. barkeri DSMZ 800(T) grown at 37°C with a temperature drop to 15°C, while the second experimental approach aimed to examine the low-temperature adaptation strategies (LTAS) of the same strain when it was grown at 15°C. The latter experiment employed cell viability and growth measurements (optical density at 600 nm [OD600]), which directly compared M. barkeri cells grown at 15°C with those grown at 37°C. During the LTSR experiment, a total of 127 proteins were detected in 37°C and 15°C samples, with 20 proteins differentially expressed with respect to temperature, while in the LTAS experiment 39% of proteins identified were differentially expressed between phases of growth. Functional categories included methanogenesis, cellular information processing, and chaperones. By applying a polyphasic approach (proteomics and growth studies), insights into the low-temperature adaptation capacity of this mesophilically characterized methanogen were obtained which suggest that the metabolically diverse Methanosarcinaceae could be functionally relevant for LTAD systems.

Two CRISPR-Cas systems in Methanosarcina mazei strain Gö1 display common processing features despite belonging to different types I and III

The clustered regularly interspaced short palindromic repeats (CRISPR) system represents a highly adaptive and heritable defense system against foreign nucleic acids in bacteria and archaea. We analyzed the two CRISPR-Cas systems in Methanosarcina mazei strain Gö1. Although belonging to different subtypes (I-B and III-B), the leaders and repeats of both loci are nearly identical. Also, despite many point mutations in each array, a common hairpin motif was identified in the repeats by a bioinformatics analysis and in vitro structural probing. The expression and maturation of CRISPR-derived RNAs (crRNAs) were studied in vitro and in vivo. Both respective potential Cas6b-type endonucleases were purified and their activity tested in vitro. Each protein showed significant activity and could cleave both repeats at the same processing site. Cas6b of subtype III-B, however, was significantly more efficient in its cleavage activity compared with Cas6b of subtype I-B. Northern blot and differential RNAseq analyses were performed to investigate in vivo transcription and maturation of crRNAs, revealing generally very low expression of both systems, whereas significant induction at high NaCl concentrations was observed. crRNAs derived proximal to the leader were generally more abundant than distal ones and in vivo processing sites were clarified for both loci, confirming the previously well-established 8 nt 5' repeat tags. The 3'-ends were more diverse, but generally ended in a prefix of the following repeat sequence (3'-tag). The analysis further revealed a 5'-hydroxy and 3'-phosphate termini architecture of small crRNAs specific for cleavage products of Cas6 endonucleases from type I-E and I-F and type III-B.

Contribution of transcriptomics to systems-level understanding of methanogenic Archaea

Methane-producing Archaea are of interest due to their contribution to atmospheric change and for their roles in technological applications including waste treatment and biofuel production. Although restricted to anaerobic environments, methanogens are found in a wide variety of habitats, where they commonly live in syntrophic relationships with bacterial partners. Owing to tight thermodynamic constraints of methanogenesis alone or in syntrophic metabolism, methanogens must carefully regulate their catabolic pathways including the regulation of RNA transcripts. The transcriptome is a dynamic and important control point in microbial systems. This paper assesses the impact of mRNA (transcriptome) studies on the understanding of methanogenesis with special consideration given to how methanogenesis is regulated to cope with nutrient limitation, environmental variability, and interactions with syntrophic partners. In comparison with traditional microarray-based transcriptome analyses, next-generation high-throughput RNA sequencing is greatly advantageous in assessing transcription start sites, the extent of 5′ untranslated regions, operonic structure, and the presence of small RNAs. We are still in the early stages of understanding RNA regulation but it is already clear that determinants beyond transcript abundance are highly relevant to the lifestyles of methanogens, requiring further study.

An archaeal sRNA targeting cis- and trans-encoded mRNAs via two distinct domains

We report on the characterization and target analysis of the small (s)RNA₁₆₂ in the methanoarchaeon Methanosarcina mazei. Using a combination of genetic approaches, transcriptome analysis and computational predictions, the bicistronic MM2441-MM2440 mRNA encoding the transcription factor MM2441 and a protein of unknown function was identified as a potential target of this sRNA, which due to processing accumulates as three stabile 5′ fragments in late exponential growth. Mobility shift assays using various mutants verified that the non-structured single-stranded linker region of sRNA₁₆₂ (SLR) base-pairs with the MM2440-MM2441 mRNA internally, thereby masking the predicted ribosome binding site of MM2441. This most likely leads to translational repression of the second cistron resulting in dis-coordinated operon expression. Analysis of mutant RNAs in vivo confirmed that the SLR of sRNA₁₆₂ is crucial for target interactions. Furthermore, our results indicate that sRNA₁₆₂-controlled MM2441 is involved in regulating the metabolic switch between the carbon sources methanol and methylamine. Moreover, biochemical studies demonstrated that the 5′ end of sRNA₁₆₂ targets the 5′-untranslated region of the cis-encoded MM2442 mRNA. Overall, this first study of archaeal sRNA/mRNA-target interactions unraveled that sRNA₁₆₂ acts as an antisense (as)RNA on cis- and trans-encoded mRNAs via two distinct domains, indicating that cis-encoded asRNAs can have larger target regulons than previously anticipated.

Global transcriptomic and proteomic responses of Dehalococcoides ethenogenes strain 195 to fixed nitrogen limitation

Bacteria of the genus Dehalococcoides play an important role in the reductive dechlorination of chlorinated ethenes. A systems-level approach was taken in this study to examine the global transcriptomic and proteomic responses of exponentially growing cells of Dehalococcoides ethenogenes strain 195 to fixed nitrogen limitation (FNL), as dechlorination activity and cell yield both decrease during FNL. As expected, the nitrogen-fixing (nif) genes were differentially upregulated in the transcriptome and proteome of strain 195 during FNL. Aside from the nif operon, a putative methylglyoxal synthase-encoding gene (DET1576), the product of which is predicted to catalyze the formation of the toxic electrophile methylglyoxal and is implicated in the uncoupling of anabolism from catabolism in bacteria, was strongly upregulated in the transcriptome and could potentially play a role in the observed growth inhibition during FNL. Carbon catabolism genes were generally downregulated in response to FNL, and a number of transporters were differentially regulated in response to nitrogen limitation, with some playing apparent roles in nitrogen acquisition, while others were associated with general stress responses. A number of genes related to the functions of nucleotide synthesis, replication, transcription, translation, and posttranslational modifications were also differentially expressed. One gene coding for a putative reductive dehalogenase (DET1545) and a number of genes coding for oxidoreductases, which have implications in energy generation and redox reactions, were also differentially regulated. Interestingly, most of the genes within the multiple integrated elements were not differentially expressed. Overall, this study elucidates the molecular responses of strain 195 to FNL and identifies differentially expressed genes that are potential biomarkers to evaluate environmental cellular nitrogen status.

Connection between multimetal(loid) methylation in methanoarchaea and central intermediates of methanogenesis

In spite of the significant impact of biomethylation on the mobility and toxicity of metals and metalloids in the environment, little is known about the biological formation of these methylated metal(loid) compounds. While element-specific methyltransferases have been isolated for arsenic, the striking versatility of methanoarchaea to methylate numerous metal(loid)s, including rare elements like bismuth, is still not understood. Here, we demonstrate that the same metal(loid)s (arsenic, selenium, antimony, tellurium, and bismuth) that are methylated by Methanosarcina mazei in vivo are also methylated by in vitro assays with purified recombinant MtaA, a methyltransferase catalyzing the methyl transfer from methylcobalamin [CH₃Cob(III)] to 2-mercaptoethanesulfonic acid (CoM) in methylotrophic methanogenesis. Detailed studies revealed that cob(I)alamin [Cob(I)], formed by MtaA-catalyzed demethylation of CH₃Cob(III), is the causative agent for the multimetal(loid) methylation observed. Moreover, Cob(I) is also capable of metal(loid) hydride generation. Global transcriptome profiling of M. mazei cultures exposed to bismuth did not reveal induced methyltransferase systems but upregulated regeneration of methanogenic cofactors in the presence of bismuth. Thus, we conclude that the multimetal(loid) methylation in vivo is attributed to side reactions of CH₃Cob(III) with reduced cofactors formed in methanogenesis. The close connection between metal(loid) methylation and methanogenesis explains the general capability of methanoarchaea to methylate metal(loid)s.

Model organisms for genetics in the domain Archaea: methanogens, halophiles, Thermococcales and Sulfolobales

The tree of life is split into three main branches: eukaryotes, bacteria, and archaea. Our knowledge of eukaryotic and bacteria cell biology has been built on a foundation of studies in model organisms, using the complementary approaches of genetics and biochemistry. Archaea have led to some exciting discoveries in the field of biochemistry, but archaeal genetics has been slow to get off the ground, not least because these organisms inhabit some of the more inhospitable places on earth and are therefore believed to be difficult to culture. In fact, many species can be cultivated with relative ease and there has been tremendous progress in the development of genetic tools for both major archaeal phyla, the Euryarchaeota and the Crenarchaeota. There are several model organisms available for methanogens, halophiles, and thermophiles; in the latter group, there are genetic systems for Sulfolobales and Thermococcales. In this review, we present the advantages and disadvantages of working with each archaeal group, give an overview of their different genetic systems, and direct the neophyte archaeologist to the most appropriate model organism.

Identification of diazotrophic microorganisms in marine sediment via fluorescence in situ hybridization coupled to nanoscale secondary ion mass spectrometry (FISH-NanoSIMS)

Growing appreciation for the biogeochemical significance of uncultured microorganisms is changing the focus of environmental microbiology. Techniques designed to investigate microbial metabolism in situ are increasingly popular, from mRNA-targeted fluorescence in situ hybridization (FISH) to the "-omics" revolution, including metagenomics, transcriptomics, and proteomics. Recently, the coupling of FISH with nanometer-scale secondary ion mass spectrometry (NanoSIMS) has taken this movement in a new direction, allowing single-cell metabolic analysis of uncultured microbial phylogenic groups. The main advantage of FISH-NanoSIMS over previous noncultivation-based techniques to probe metabolism is its ability to directly link 16S rRNA phylogenetic identity to metabolic function. In the following chapter, we describe the procedures necessary to identify nitrogen-fixing microbes within marine sediment via FISH-NanoSIMS, using our work on nitrogen fixation by uncultured deep-sea methane-consuming archaea as a case study.

Studying gene regulation in methanogenic archaea

Methanogenic archaea are a unique group of strictly anaerobic microorganisms characterized by their ability, and dependence, to convert simple C1 and C2 compounds to methane for growth. The major models for studying the biology of methanogens are members of the Methanococcus and Methanosarcina species. Recent development of sophisticated tools for molecular analysis and for genetic manipulation allows investigating not only their metabolism but also their cell cycle, and their interaction with the environment in great detail. One aspect of such analyses is assessment and dissection of methanoarchaeal gene regulation, for which, at present, only a handful of cases have been investigated thoroughly, partly due to the great methodological effort required. However, it becomes more and more evident that many new regulatory paradigms can be unraveled in this unique archaeal group. Here, we report both molecular and physiological/genetic methods to assess gene regulation in Methanococcus maripaludis and Methanosarcina acetivorans, which should, however, be applicable for other methanogens as well.

Differentially expressed genes after hyper- and hypo-salt stress in the halophilic archaeonMethanohalophilus portucalensis

Methanohalophilus portucalensis FDF1 can grow over a range of external NaCl concentrations, from 1.2 to 2.9 mol/L. Differential gene expression in response to long-term hyper-salt stress (3.1 mol/L of NaCl) and hypo-salt stress (0.9 mol/L of NaCl) were compared by differential display RT-PCR. Fourteen differentially expressed genes responding to long-term hyper- or hypo-salt stress were detected, cloned, and sequenced. Several of the differentially expressed genes were related to the unique energy-acquiring methanogenesis pathway in this organism, including the transmembrane protein MttP, cobalamin biosynthesis protein, methenyl-H4MPT cyclohydrolase and monomethylamine methyltransferase. One signal transduction histidine kinase was identified from the hyper-salt stress cultures. Moreover, 3 known stress-response gene homologues — the DNA mismatch repair protein, MutS, the universal stress protein, UspA, and a member of the protein-disaggregating multichaperone system, ClpB — were also detected. The transcriptional analysis of these long-term salt stress response and adaptation-related genes for cells immediately after salt stress indicated that the expression of the energy metabolism genes was arrested during hyper-salt shock, while the chaperone clpB gene was stimulated by both hypo- and hyper-salt shock.

Metatranscriptomic array analysis of 'Candidatus Accumulibacter phosphatis'-enriched enhanced biological phosphorus removal sludge

Here we report the first metatranscriptomic analysis of gene expression and regulation of 'Candidatus Accumulibacter'-enriched lab-scale sludge during enhanced biological phosphorus removal (EBPR). Medium density oligonucleotide microarrays were generated with probes targeting most predicted genes hypothesized to be important for the EBPR phenotype. RNA samples were collected at the early stage of anaerobic and aerobic phases (15 min after acetate addition and switching to aeration respectively). We detected the expression of a number of genes involved in the carbon and phosphate metabolisms, as proposed by EBPR models (e.g. polyhydroxyalkanoate synthesis, a split TCA cycle through methylmalonyl-CoA pathway, and polyphosphate formation), as well as novel genes discovered through metagenomic analysis. The comparison between the early stage anaerobic and aerobic gene expression profiles showed that expression levels of most genes were not significantly different between the two stages. The majority of upregulated genes in the aerobic sample are predicted to encode functions such as transcription, translation and protein translocation, reflecting the rapid growth phase of Accumulibacter shortly after being switched to aerobic conditions. Components of the TCA cycle and machinery involved in ATP synthesis were also upregulated during the early aerobic phase. These findings support the predictions of EBPR metabolic models that the oxidation of intracellularly stored carbon polymers through the TCA cycle provides ATP for cell growth when oxygen becomes available. Nitrous oxide reductase was among the very few Accumulibacter genes upregulated in the anaerobic sample, suggesting that its expression is likely induced by the deprivation of oxygen.

Hot Transcriptomics

DNA microarray technology allows for a quick and easy comparison of complete transcriptomes, resulting in improved molecular insight in fluctuations of gene expression. After emergence of the microarray technology about a decade ago, the technique has now matured and has become routine in many molecular biology laboratories. Numerous studies have been performed that have provided global transcription patterns of many organisms under a wide range of conditions. Initially, implementation of this high-throughput technology has lead to high expectations for ground breaking discoveries. Here an evaluation is performed of the insight that transcriptome analysis has brought about in the field of hyperthermophilic archaea. The examples that will be discussed have been selected on the basis of their impact, in terms of either biological insight or technological progress.

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

Methanosarcina mazei and related mesophilic archaea are the only organisms fermenting acetate, methylamines, and methanol to methane and carbon dioxide, contributing significantly to greenhouse gas production. The biochemistry of these metabolic processes is well studied, and genome sequences are available, yet little is known about the overall transcriptional organization and the noncoding regions representing 25% of the 4.01-Mb genome of M. mazei. We present a genome-wide analysis of transcription start sites (TSS) in M. mazei grown under different nitrogen availabilities. Pyrosequencing-based differential analysis of primary vs. processed 5' ends of transcripts discovered 876 TSS across the M. mazei genome. Unlike in other archaea, in which leaderless mRNAs are prevalent, the majority of the detected mRNAs in M. mazei carry long untranslated 5' regions. Our experimental data predict a total of 208 small RNA (sRNA) candidates, mostly from intergenic regions but also antisense to 5' and 3' regions of mRNAs. In addition, 40 new small mRNAs with ORFs of < or = 30 aa were identified, some of which might have dual functions as mRNA and regulatory sRNA. We confirmed differential expression of several sRNA genes in response to nitrogen availability. Inspection of their promoter regions revealed a unique conserved sequence motif associated with nitrogen-responsive regulation, which might serve as a regulator binding site upstream of the common IIB recognition element. Strikingly, several sRNAs antisense to mRNAs encoding transposases indicate nitrogen-dependent transposition events. This global TSS map in archaea will facilitate a better understanding of transcriptional and posttranscriptional control in the third domain of life.

A 5' leader sequence regulates expression of methanosarcinal CO dehydrogenase/acetyl coenzyme A synthase

In vivo expression of CO dehydrogenase/acetyl coenzyme A synthase in Methanosarcina spp. is coordinately regulated in response to substrate by at least two mechanisms: differential transcription initiation and early elongation termination near the 3' end of a 371-bp leader sequence. This is the first report of regulation of transcription elongation in the Archaea.

Expression of ribonuclease A and ribonuclease N1 in the filamentous fungus Neurospora crassa

In this study, we investigated the ability of the fungus Neurospora crassa to produce and secrete two ribonucleases: the heterologous bovine RNase A and the endogenous RNase N(1). A set of expression vectors was constructed, each consisting of an RNase A open reading frame under the control of a specific promoter and each with a specific terminator. N. crassa transformants were analyzed at the transcriptional and protein levels. Irrespective of the promoter used, all transformants showed an RNase A-specific transcript in northern hybridization, but transcriptional strengths differed significantly. The strongest transcription was detected in transformants under the control of the cfp promoter. Western blot analysis and ELISA assays of selected transformants showed an effective secretion up to 356 ng/mL of recombinant RNase A protein. However, the highest ribonuclease activity could be detected in transformants carrying the endogenous RNase N(1) under the control of the ccg1 promoter. Expression and secretion of RNase N(1) thus represent an alternative to recombinant expression of RNase A protein. In conclusion, we have created a viable expression system for expression of homologous and heterologous proteins in N. crassa.

Quantitative proteomics of nutrient limitation in the hydrogenotrophic methanogen Methanococcus maripaludis

Background

Methanogenic Archaea play key metabolic roles in anaerobic ecosystems, where they use H₂ and other substrates to produce methane. Methanococcus maripaludis is a model for studies of the global response to nutrient limitations.

Results

We used high-coverage quantitative proteomics to determine the response of M. maripaludis to growth-limiting levels of H₂, nitrogen, and phosphate. Six to ten percent of the proteome changed significantly with each nutrient limitation. H₂ limitation increased the abundance of a wide variety of proteins involved in methanogenesis. However, one protein involved in methanogenesis decreased: a low-affinity [Fe] hydrogenase, which may dominate over a higher-affinity mechanism when H₂ is abundant. Nitrogen limitation increased known nitrogen assimilation proteins. In addition, the increased abundance of molybdate transport proteins suggested they function for nitrogen fixation. An apparent regulon governed by the euryarchaeal nitrogen regulator NrpR is discussed. Phosphate limitation increased the abundance of three different sets of proteins, suggesting that all three function in phosphate transport.

Conclusion

The global proteomic response of M. maripaludis to each nutrient limitation suggests a wider response than previously appreciated. The results give new insight into the function of several proteins, as well as providing information that should contribute to the formulation of a regulatory network model.

A novel limb in the osmoregulatory network of Methanosarcina mazei Gö1: N(epsilon)-acetyl-beta-lysine can be substituted by glutamate and alanine

N(epsilon)-acetyl-beta-lysine is a unique compatible solute found in methanogenic archaea grown at high salinities. Deletion of the genes that encode the lysine-2,3-aminomutase (ablA) and the beta-lysine acetyltransferase (ablB) abolished the production of N(epsilon)-acetyl-beta-lysine in Methanosarcina mazei Gö1. The mutant grew well at low and intermediate salinities. Interestingly, growth at high salt (800 mM NaCl) was only slowed down but not impaired demonstrating that in M. mazei Gö1 N(epsilon)-acetyl-beta-lysine is not essential for growth at high salinities. Nuclear magnetic resonance (NMR) analysis revealed an increased glutamate pool in the mutant. In addition to alpha-glutamate, a novel solute, alanine, was produced. The intracellular alanine concentration was as high as 0.36 +/- 0.05 micromol (mg protein)-1 representing up to 18% of the total solute pool at 800 mM NaCl. The cellular alanine concentration increased with the salinity of the medium and decreased in the presence of glycine betaine in the medium, indicating that alanine is used as compatible solute by M. mazei Gö1.

Transcriptional profiling of methyltransferase genes during growth of Methanosarcina mazei on trimethylamine

The genomic expression patterns of Methanosarcina mazei growing with trimethylamine were measured in comparison to those of cells grown with methanol. We identified a total of 72 genes with either an increased level (49 genes) or a decreased level (23 genes) of mRNA during growth on trimethylamine with methanol-grown cells as the control. Major differences in transcript levels were observed for the mta, mtb, mtt, and mtm genes, which encode enzymes involved in methane formation from methanol and trimethylamine, respectively. Other differences in mRNA abundance were found for genes encoding enzymes involved in isopentenyl pyrophosphate synthesis and in the formation of aromatic amino acids, as well as a number of proteins with unknown functions. The results were verified by in-depth analysis of methyltransferase genes using specific primers for real-time quantitative reverse transcription-PCR (RT-PCR). The monitored transcript levels of genes encoding corrinoid proteins involved in methyl group transfer from methylated C(1) compounds (mtaC, mtbC, mttC, and mtmC) indicated increased amounts of mRNA from the mtaBC1, mtaBC2, and mtaBC3 operons in methanol-grown cells, whereas mRNA of the mtb1-mtt1 operon was found in high concentrations during trimethylamine consumption. The genes of the mtb1-mtt1 operon encode methyltransferases that are responsible for sequential demethylation of trimethylamine. The analysis of product formation of trimethylamine-grown cells at different optical densities revealed that large amounts of dimethylamine and monomethylamine were excreted into the medium. The intermediate compounds were consumed only in the very late exponential growth phase. RT-PCR analysis of key genes involved in methanogenesis led to the conclusion that M. mazei is able to adapt to changing trimethylamine concentrations and the consumption of intermediate compounds. Hence, we assume that the organism possesses a regulatory network for optimal substrate utilization.

The molecular basis of salt adaptation in Methanosarcina mazei Gö1

The study on the molecular basis of salt adaptation and its regulation in archaea is still in its infancy, but genomics and functional genome analyses combined with classical biochemistry shed light on the processes conferring salt adaptation in the methanogenic archaeon Methanosarcina mazei Gö1. In this article, we will review discoveries made within the last years that will culminate in the description of the overall cellular response of M. mazei Gö1 to elevated salinities. This response includes accumulation of solutes and export of Na+ as well as potential uptake/export of K+ but also a restructuring of the cell surface.

Deletion of the archaeal histone in Methanosarcina mazei Gö1 results in reduced growth and genomic transcription

HMm is the only archaeal histone in Methanosarcina mazei Göl and recombinant HMm, synthesized by expression of MM1825 in Escherichia coli, has been purified and confirmed to have the DNA binding and compaction properties characteristic of an archaeal histone. Insertion of a puromycin resistance conferring cassette (pac) into MM1825 was not lethal but resulted in mutants (M. mazei MM1825::pac) that have impaired ability to grow on methanol and trimethylamine. Loss of HMm also resulted in increased sensitivity to UV light and decreased transcript levels for approximately 25% of all M. mazei genes. For most genes, the transcript decrease was 3- to 10-fold, but transcripts of MM483 (small heat-shock protein), MM1688 (trimethylamine:corrinoid methyl transferase) and MM3195 (transcription regulator), were reduced 100-, 100- and 25-fold, respectively, in M. mazei MM1825::pac cells. Transcripts of only five adjacent genes that appear to constitute an aromatic amino acid biosynthetic operon were elevated in M. mazei MM1825::pac cells. Complementary synthesis of HMm from a plasmid transformed into M. mazei MM1825::pac restored wild-type growth and transcript levels.

Nitrogen Regulation in Bacteria and Archaea

A wide range of Bacteria and Archaea sense cellular 2-oxoglutarate (2OG) as an indicator of nitrogen limitation. 2OG sensor proteins are varied, but most of those studied belong to the PII superfamily. Within the PII superfamily, GlnB and GlnK represent a widespread family of homotrimeric proteins (GlnB-K) that bind and respond to 2OG and ATP. In some bacterial phyla, GlnB-K proteins are covalently modified, depending on enzymes that sense cellular glutamine as an indicator of nitrogen sufficiency. GlnB-K proteins are central clearing houses of nitrogen information and bind and modulate a variety of nitrogen assimilation regulators and enzymes. NifI(1) and NifI(2) comprise a second widespread family of PII proteins (NifI) that are heteromultimeric, respond to 2OG and ATP, and bind and regulate dinitrogenase in Euryarchaeota and many Bacteria.

Diverse homologues of the archaeal repressor NrpR function similarly in nitrogen regulation

NrpR is a transcriptional repressor of nitrogen assimilation genes that was recently discovered and characterized in the methanogenic archaeon Methanococcus maripaludis. NrpR homologues are widely distributed in Euryarchaeota and present in a few bacterial species. They exist in three different domain configurations: a single ORF encoding one NrpR domain following an N-terminal helix-turn-helix (HTH); a single ORF encoding two NrpR domains fused in tandem following an N-terminal HTH; and two separate ORFs, one with a single domain following an N-terminal HTH and one with a single domain without a HTH. Phylogenetic analysis indicated that the NrpR family forms five distinct groups: the single domain HTH type, the two domains of the double domain HTH type and the two separately encoded domains. To determine the function of diverse NrpR homologues, representative genes in were expressed an Methanococcus maripaludis nrpR deletion mutant. Homologues from species that possess a single gene restored regulated repression, regardless of domain structure. In the case of Methanosarcina acetivorans that contains two genes, both were required. The results show that distantly related NrpR homologues that are present in widely dispersed phyla regulate the expression of nitrogen assimilation genes in a similar fashion.

Quantitative proteomic and microarray analysis of the archaeon Methanosarcina acetivorans grown with acetate versus methanol

Methanosarcina acetivorans strain C2A is an acetate- and methanol-utilizing methane-producing organism for which the genome, the largest yet sequenced among the Archaea, reveals extensive physiological diversity. LC linear ion trap-FTICR mass spectrometry was employed to analyze acetate- vs methanol-grown cells metabolically labeled with 14N vs 15N, respectively, to obtain quantitative protein abundance ratios. DNA microarray analyses of acetate- vs methanol-grown cells was also performed to determine gene expression ratios. The combined approaches were highly complementary, extending the physiological understanding of growth and methanogenesis. Of the 1081 proteins detected, 255 were > or =3-fold differentially abundant. DNA microarray analysis revealed 410 genes that were > or =2.5-fold differentially expressed of 1972 genes with detected expression. The ratios of differentially abundant proteins were in good agreement with expression ratios of the encoding genes. Taken together, the results suggest several novel roles for electron transport components specific to acetate-grown cells, including two flavodoxins each specific for growth on acetate or methanol. Protein abundance ratios indicated that duplicate CO dehydrogenase/acetyl-CoA complexes function in the conversion of acetate to methane. Surprisingly, the protein abundance and gene expression ratios indicated a general stress response in acetate- vs methanol-grown cells that included enzymes specific for polyphosphate accumulation and oxidative stress. The microarray analysis identified transcripts of several genes encoding regulatory proteins with identity to the PhoU, MarR, GlnK, and TetR families commonly found in the Bacteria domain. An analysis of neighboring genes suggested roles in controlling phosphate metabolism (PhoU), ammonia assimilation (GlnK), and molybdopterin cofactor biosynthesis (TetR). Finally, the proteomic and microarray results suggested roles for two-component regulatory systems specific for each growth substrate.

Cell envelope stress response in Bacillus licheniformis: integrating comparative genomics, transcriptional profiling, and regulon mining to decipher a complex regulatory network

The envelope is an essential structure of the bacterial cell, and maintaining its integrity is a prerequisite for survival. To ensure proper function, transmembrane signal-transducing systems, such as two-component systems (TCS) and extracytoplasmic function (ECF) sigma factors, closely monitor its condition and respond to harmful perturbations. Both systems consist of a transmembrane sensor protein (histidine kinase or anti-sigma factor, respectively) and a corresponding cytoplasmic transcriptional regulator (response regulator or sigma factor, respectively) that mediates the cellular response through differential gene expression. The regulatory network of the cell envelope stress response is well studied in the gram-positive model organism Bacillus subtilis. It consists of at least two ECF sigma factors and four two-component systems. In this study, we describe the corresponding network in a close relative, Bacillus licheniformis. Based on sequence homology, domain architecture, and genomic context, we identified five TCS and eight ECF sigma factors as potential candidate regulatory systems mediating cell envelope stress response in this organism. We characterized the corresponding regulatory network by comparative transcriptomics and regulon mining as an initial screening tool. Subsequent in-depth transcriptional profiling was applied to define the inducer specificity of each identified cell envelope stress sensor. A total of three TCS and seven ECF sigma factors were shown to be induced by cell envelope stress in B. licheniformis. We noted a number of significant differences, indicative of a regulatory divergence between the two Bacillus species, in addition to the expected overlap in the respective responses.