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

Expression of V-nitrogenase and Fe-nitrogenase in Methanosarcina acetivorans is controlled by molybdenum, fixed nitrogen, and the expression of Mo-nitrogenase

All nitrogen-fixing bacteria and archaea (diazotrophs) use molybdenum (Mo) nitrogenase to reduce dinitrogen (N2) to ammonia, with some also containing vanadium (V) and iron-only (Fe) nitrogenases that lack Mo. Among diazotrophs, the regulation and usage of the alternative V-nitrogenase and Fe-nitrogenase in methanogens are largely unknown. Methanosarcina acetivorans contains nif, vnf, and anf gene clusters encoding putative Mo-nitrogenase, V-nitrogenase, and Fe-nitrogenase, respectively. This study investigated nitrogenase expression and growth by M. acetivorans in response to fixed nitrogen, Mo/V availability, and CRISPRi repression of the nif, vnf, and/or anf gene clusters. The availability of Mo and V significantly affected growth of M. acetivorans with N2 but not with NH4Cl. M. acetivorans exhibited the fastest growth rate and highest cell yield during growth with N2 in medium containing Mo, and the slowest growth in medium lacking Mo and V. qPCR analysis revealed the transcription of the nif operon is only moderately affected by depletion of fixed nitrogen and Mo, whereas vnf and anf transcription increased significantly when fixed nitrogen and Mo were depleted, with removal of Mo being key. Immunoblot analysis revealed Mo-nitrogenase is detected when fixed nitrogen is depleted regardless of Mo availability, while V-nitrogenase and Fe-nitrogenase are detected only in the absence of fixed nitrogen and Mo. CRISPRi repression studies revealed that V-nitrogenase and/or Fe-nitrogenase are required for Mo-independent diazotrophy, and unexpectedly that the expression of Mo-nitrogenase is also required. These results reveal that alternative nitrogenase production in M. acetivorans is tightly controlled and dependent on Mo-nitrogenase expression. IMPORTANCE Methanogens and closely related methanotrophs are the only archaea known or predicted to possess nitrogenase. Methanogens play critical roles in both the global biological nitrogen and carbon cycles. Moreover, methanogens are an ancient microbial lineage and nitrogenase likely originated in methanogens. An understanding of the usage and properties of nitrogenases in methanogens can provide new insight into the evolution of nitrogen fixation and aid in the development nitrogenase-based biotechnology. This study provides the first evidence that a methanogen can produce all three forms of nitrogenases, including simultaneously. The results reveal components of Mo-nitrogenase regulate or are needed to produce V-nitrogenase and Fe-nitrogenase in methanogens, a result not seen in bacteria. Overall, this study provides a foundation to understand the assembly, regulation, and activity of the alternative nitrogenases in methanogens.

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.

A CRISPRi-dCas9 System for Archaea and Its Use To Examine Gene Function during Nitrogen Fixation by Methanosarcina acetivorans

CRISPR-based systems are emerging as the premier method to manipulate many cellular processes. In this study, a simple and efficient CRISPR interference (CRISPRi) system for targeted gene repression in archaea was developed. The Methanosarcina acetivorans CRISPR-Cas9 system was repurposed by replacing Cas9 with the catalytically dead Cas9 (dCas9) to generate a CRISPRi-dCas9 system for targeted gene repression. To test the utility of the system, genes involved in nitrogen (N₂) fixation were targeted for dCas9-mediated repression. First, the nif operon (nifHI ₁ I ₂ DKEN) that encodes molybdenum nitrogenase was targeted by separate guide RNAs (gRNAs), one targeting the promoter and the other targeting nifD Remarkably, growth of M. acetivorans with N₂ was abolished by dCas9-mediated repression of the nif operon with each gRNA. The abundance of nif transcripts was >90% reduced in both strains expressing the gRNAs, and NifD was not detected in cell lysate. Next, we targeted NifB, which is required for nitrogenase cofactor biogenesis. Expression of a gRNA targeting the coding sequence of NifB decreased nifB transcript abundance >85% and impaired but did not abolish growth of M. acetivorans with N₂ Finally, to ascertain the ability to study gene regulation using CRISPRi-dCas9, nrpR1, encoding a subunit of the repressor of the nif operon, was targeted. The nrpR1 repression strain grew normally with N₂ but had increased nif operon transcript abundance, consistent with NrpR1 acting as a repressor. These results highlight the utility of the system, whereby a single gRNA when expressed with dCas9 can block transcription of targeted genes and operons in M. acetivorans IMPORTANCE Genetic tools are needed to understand and manipulate the biology of archaea, which serve critical roles in the biosphere. Methanogenic archaea (methanogens) are essential for the biological production of methane, an intermediate in the global carbon cycle, an important greenhouse gas, and a biofuel. The CRISPRi-dCas9 system in the model methanogen Methanosarcina acetivorans is, to our knowledge, the first Cas9-based CRISPR interference system in archaea. Results demonstrate that the system is remarkably efficient in targeted gene repression and provide new insight into nitrogen fixation by methanogens, the only archaea with nitrogenase. Overall, the CRISPRi-dCas9 system provides a simple, yet powerful, genetic tool to control the expression of target genes and operons in methanogens.

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.

Exploring Hydrogenotrophic Methanogenesis: a Genome Scale Metabolic Reconstruction of Methanococcus maripaludis

Hydrogenotrophic methanogenesis occurs in multiple environments, ranging from the intestinal tracts of animals to anaerobic sediments and hot springs. Energy conservation in hydrogenotrophic methanogens was long a mystery; only within the last decade was it reported that net energy conservation for growth depends on electron bifurcation. In this work, we focus on Methanococcus maripaludis, a well-studied hydrogenotrophic marine methanogen. To better understand hydrogenotrophic methanogenesis and compare it with methylotrophic methanogenesis that utilizes oxidative phosphorylation rather than electron bifurcation, we have built iMR539, a genome scale metabolic reconstruction that accounts for 539 of the 1,722 protein-coding genes of M. maripaludis strain S2. Our reconstructed metabolic network uses recent literature to not only represent the central electron bifurcation reaction but also incorporate vital biosynthesis and assimilation pathways, including unique cofactor and coenzyme syntheses. We show that our model accurately predicts experimental growth and gene knockout data, with 93% accuracy and a Matthews correlation coefficient of 0.78. Furthermore, we use our metabolic network reconstruction to probe the implications of electron bifurcation by showing its essentiality, as well as investigating the infeasibility of aceticlastic methanogenesis in the network. Additionally, we demonstrate a method of applying thermodynamic constraints to a metabolic model to quickly estimate overall free-energy changes between what comes in and out of the cell. Finally, we describe a novel reconstruction-specific computational toolbox we created to improve usability. Together, our results provide a computational network for exploring hydrogenotrophic methanogenesis and confirm the importance of electron bifurcation in this process.

IMPORTANCE

Understanding and applying hydrogenotrophic methanogenesis is a promising avenue for developing new bioenergy technologies around methane gas. Although a significant portion of biological methane is generated through this environmentally ubiquitous pathway, existing methanogen models portray the more traditional energy conservation mechanisms that are found in other methanogens. We have constructed a genome scale metabolic network of Methanococcus maripaludis that explicitly accounts for all major reactions involved in hydrogenotrophic methanogenesis. Our reconstruction demonstrates the importance of electron bifurcation in central metabolism, providing both a window into hydrogenotrophic methanogenesis and a hypothesis-generating platform to fuel metabolic engineering efforts.

The Emergence of 2-Oxoglutarate as a Master Regulator Metabolite

The metabolite 2-oxoglutarate (also known as α-ketoglutarate, 2-ketoglutaric acid, or oxoglutaric acid) lies at the intersection between the carbon and nitrogen metabolic pathways. This compound is a key intermediate of one of the most fundamental biochemical pathways in carbon metabolism, the tricarboxylic acid (TCA) cycle. In addition, 2-oxoglutarate also acts as the major carbon skeleton for nitrogen-assimilatory reactions. Experimental data support the conclusion that intracellular levels of 2-oxoglutarate fluctuate according to nitrogen and carbon availability. This review summarizes how nature has capitalized on the ability of 2-oxoglutarate to reflect cellular nutritional status through evolution of a variety of 2-oxoglutarate-sensing regulatory proteins. The number of metabolic pathways known to be regulated by 2-oxoglutarate levels has increased significantly in recent years. The signaling properties of 2-oxoglutarate are highlighted by the fact that this metabolite regulates the synthesis of the well-established master signaling molecule, cyclic AMP (cAMP), in Escherichia coli.

Comparative genomics of three Methanocellales strains reveal novel taxonomic and metabolic features

Methanocellales represents a new order of methanogens, which is widespread in environments and plays specifically the important role in methane emissions from paddy fields. To gain more insights into Methanocellales, comparative genomic studies were performed among three Methanocellales strains through the same annotation pipeline. Genetic relationships among strains revealed by genome alignment, pan-genome reconstruction and comparison of amino average identity suggest that they should be classified in different genera. In addition, multiple copies of cell cycle regulator proteins were identified for the first time in Archaea. Core metabolisms were reconstructed, predicting certain unique and novel features for Methanocellales, including a set of methanogenesis genes potentially organized toward specialization in utilizing low concentrations of H2, a new route of disulfide reduction catalysed by a disulfide-reducing hydrogenase (Drh) complex phylogenetically related to sulfate-reducing prokaryotes, an oxidative tricarboxylic acid (TCA) cycle, a sophisticated nitrogen uptake and regulation system as well as a versatile sulfur utilization system. These core metabolisms are largely conserved among the three strains, but differences in gene copy number and metabolic diversity are evident. The present study thus adds new dimensions to the unique ecophysiology of Methanocellales and offers a road map for further experimental characterization of this methanogen lineage.

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.

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.

Transcription regulation in the third domain

The ability of organisms to sense and respond to their environment is essential to their survival. This is no different for members of the third domain of life, the Archaea. Archaea are found in diverse and often extreme habitats. However, their ability to sense and respond to their environment at the level of gene expression has been understudied when compared to bacteria and eukaryotes. Over the last decade, the field has expanded, and a variety of unique and interesting regulatory schemes have been unraveled. In this review, the current state of knowledge of archaeal transcription regulation is explored.

Random mutagenesis identifies factors involved in formate-dependent growth of the methanogenic archaeon Methanococcus maripaludis

Methane is a key intermediate in the carbon cycle and biologically produced by methanogenic archaea. Most methanogens are able to conserve energy by reducing CO2 to methane using molecular hydrogen as electron donor (hydrogenotrophic methanogenesis), but several hydrogenotrophic methanogens can also use formate as electron donor for methanogenesis. Formate dehydrogenase (Fdh) oxidizes formate to CO2 and is involved in funneling reducing equivalents into the methanogenic pathway, but details on other factors relevant for formate-dependent physiology of methanogens are not available. To learn more about the factors involved in formate-dependent growth of Methanococcus maripaludis strain JJ, we used a recently developed system for random in vitro mutagenesis, which is based on a modified insect transposable element to create 2,865 chromosomal transposon mutants and screened them for impaired growth on formate. Of 12 M. maripaludis transposon-induced mutants exhibiting this phenotype, the transposon insertion sites in the chromosome were mapped. Among the genes, apparently affecting formate-dependent growth were those encoding archaeal transcription factor S, a regulator of ion transport, and carbon monoxide dehydrogenase/acetyl-CoA synthase. Interestingly, in seven of the mutants, transposons were localized in a 10.2 kb region where Fdh1, one of two Fdh isoforms in the organism, is encoded. Two transcription start sites within the 10.2 kb region could be mapped, and quantification of transcripts revealed that transposon insertion in this region diminished fdhA1 expression due to polar effects.

Disruption and complementation of the selenocysteine biosynthesis pathway reveals a hierarchy of selenoprotein gene expression in the archaeon Methanococcus maripaludis

Proteins containing selenocysteine are found in members of all three domains of life, Bacteria, Eukarya and Archaea. A dedicated tRNA (tRNA(sec)) serves as a scaffold for selenocysteine synthesis. However, sequence and secondary structures differ in tRNA(sec) from the different domains. An Escherichia coli strain lacking the gene for tRNA(sec) could only be complemented with the homologue from Methanococcus maripaludis when a single base in the anticodon loop was exchanged demonstrating that this base is a crucial determinant for archaeal tRNA(sec) to function in E. coli. Complementation in trans of M. maripaludis JJ mutants lacking tRNA(sec) , O-phosphoseryl-tRNA(sec) kinase or O-phosphoseryl-tRNA(sec) :selenocysteine synthase with the corresponding genes from M. maripaludis S2 restored the mutant's ability to synthesize selenoproteins. However, only partial restoration of the wild-type selenoproteome was observed as only selenocysteine-containing formate dehydrogenase was synthesized. Quantification of transcripts showed that disrupting the pathway of selenocysteine synthesis leads to downregulation of selenoprotein gene expression, concomitant with upregulation of a selenium-independent backup system, which is not re-adjusted upon complementation. This transcriptional arrest was independent of selenophosphate but depended on the 'history' of the mutants and was inheritable, which suggests that a stable genetic switch may cause the resulting hierarchy of selenoproteins synthesized.

Overlapping repressor binding sites regulate expression of the Methanococcus maripaludis glnK(1) operon

The euryarchaeal transcriptional repressor NrpR regulates a variety of nitrogen assimilation genes by 2-oxoglutarate-reversible binding to conserved palindromic operators. The number and positioning of these operators varies among promoter regions of regulated genes, suggesting NrpR can bind in different patterns. Particularly intriguing is the contrast between the nif and glnK(1) promoter regions of Methanococcus maripaludis, where two operators are present but with different configurations. Here we study NrpR binding and regulation at the glnK(1) promoter, where the two operator sequences overlap and occur on opposite faces of the double helix. We find that both operators function in binding, with a dimer of NrpR binding simultaneously to each overlapping operator. We show in vivo that the first operator plays a primary role in regulation and the second operator plays an enhancing role. This is the first demonstration of overlapping operators functioning in Archaea.

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.

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.

Genetic screen for regulatory mutations in Methanococcus maripaludis and its use in identification of induction-deficient mutants of the euryarchaeal repressor NrpR

NrpR is an euryarchaeal transcriptional repressor of nitrogen assimilation genes. Previous studies with Methanococcus maripaludis demonstrated that NrpR binds to palindromic operator sequences, blocking transcription initiation. The metabolite 2-oxoglutarate, an indicator of cellular nitrogen deficiency, induces transcription by lowering the affinity of NrpR for operator DNA. In this report we build on existing genetic tools for M. maripaludis to develop a screen for change-of-function mutations in a transcriptional regulator and demonstrate the use of an X-Gal (5-bromo-4-chloro-3-indolyl-beta-d-galactopyranoside) screen for strict anaerobes. We use the approach to address the primary structural requirements for the response of NrpR to 2-oxoglutarate. nrpR genes from the mesophilic M. maripaludis and the hyperthermophilic Methanopyrus kandleri were targeted for mutagenesis. M. maripaludis nrpR encodes a protein with two homologous NrpR domains while the M. kandleri nrpR homolog encodes a single NrpR domain. Random point mutagenesis and alanine replacement mutagenesis identified two amino acid residues of M. kandleri NrpR involved in induction of gene expression under nitrogen-deficient conditions and thus in the response to 2-oxoglutarate. Mutagenesis of the corresponding regions in either domain of M. maripaludis NrpR resulted in a similar effect, demonstrating a conserved structure-function relationship between the two repressors. The results indicate that in M. maripaludis, both NrpR domains participate in the 2-oxoglutarate response. The approach used here has wide adaptability to other regulatory systems in methanogenic Archaea and other strict anaerobes.