Effects of nitrogen and carbon sources on transcription of soluble methyltransferases in Methanosarcina mazei strain Go1

Amine-recognizing domain in diverse receptors from bacteria and archaea evolved from the universal amino acid sensor

Bacteria possess various receptors that sense different signals and transmit information to enable an optimal adaptation to the environment. A major limitation in microbiology is the lack of information on the signal molecules that activate receptors. Signals recognized by sensor domains are poorly reflected in overall sequence identity, and therefore, the identification of signals from the amino acid sequence of the sensor alone presents a challenge. Biogenic amines are of great physiological importance for microorganisms and humans. They serve as substrates for aerobic and anaerobic growth and play a role of neurotransmitters and osmoprotectants. Here, we report the identification of a sequence motif that is specific for amine-sensing sensor domains that belong to the Cache superfamily of the most abundant extracellular sensors in prokaryotes. We identified approximately 13,000 sensor histidine kinases, chemoreceptors, receptors involved in second messenger homeostasis and Ser/Thr phosphatases from 8,000 bacterial and archaeal species that contain the amine-recognizing motif. The screening of compound libraries and microcalorimetric titrations of selected sensor domains confirmed their ability to specifically bind biogenic amines. Mutants in the amine-binding motif or domains that contain a single mismatch in the binding motif had either no or a largely reduced affinity for amines. We demonstrate that the amine-recognizing domain originated from the universal amino acid-sensing Cache domain, thus providing insight into receptor evolution. Our approach enables precise "wet"-lab experiments to define the function of regulatory systems and therefore holds a strong promise to enable the identification of signals stimulating numerous receptors.

Amine recognizing domain in diverse receptors from bacteria and archaea evolved from the universal amino acid sensor

Bacteria contain many different receptor families that sense different signals permitting an optimal adaptation to the environment. A major limitation in microbiology is the lack of information on the signal molecules that activate receptors. Due to a significant sequence divergence, the signal recognized by sensor domains is only poorly reflected in overall sequence identity. Biogenic amines are of central physiological relevance for microorganisms and serve for example as substrates for aerobic and anaerobic growth, neurotransmitters or osmoprotectants. Based on protein structural information and sequence analysis, we report here the identification of a sequence motif that is specific for amine-sensing dCache sensor domains (dCache_1AM). These domains were identified in more than 13,000 proteins from 8,000 bacterial and archaeal species. dCache_1AM containing receptors were identified in all major receptor families including sensor kinases, chemoreceptors, receptors involved in second messenger homeostasis and Ser/Thr phosphatases. The screening of compound libraries and microcalorimetric titrations of selected dCache_1AM domains confirmed their capacity to specifically bind amines. Mutants in the amine binding motif or domains that contain a single mismatch in the binding motif, had either no or a largely reduced affinity for amines, illustrating the specificity of this motif. We demonstrate that the dCache_1AM domain has evolved from the universal amino acid sensing domain, providing novel insight into receptor evolution. Our approach enables precise "wet"-lab experiments to define the function of regulatory systems and thus holds a strong promise to address an important bottleneck in microbiology: the identification of signals that stimulate numerous receptors.

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

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

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

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

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.

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

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

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.

Community and Proteomic Analysis of Anaerobic Consortia Converting Tetramethylammonium to Methane

Tetramethylammonium-degrading methanogenic consortia from a complete-mixing suspended sludge (CMSS) and an upflow anaerobic sludge blanket (UASB) reactors were studied using multiple PCR-based molecular techniques and shotgun proteomic approach. The prokaryotic 16S rRNA genes of the consortia were analyzed by quantitative PCR, high-throughput sequencing, and DGGE-cloning methods. The results showed that methanogenic archaea were highly predominant in both reactors but differed markedly according to community structure. Community and proteomic analysis revealed that Methanomethylovorans and Methanosarcina were the major players for the demethylation of methylated substrates and methane formation through the reduction pathway of methyl-S-CoM and possibly, acetyl-CoA synthase/decarbonylase-related pathways. Unlike high dominance of one Methanomethylovorans population in the CMSS reactor, diverse methylotrophic Methanosarcina species inhabited in syntrophy-like association with hydrogenotrophic Methanobacterium in the granular sludge of UASB reactor. The overall findings indicated the reactor-dependent community structures of quaternary amines degradation and provided microbial insight for the improved understanding of engineering application.

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.

Reducing the genetic code induces massive rearrangement of the proteome

Expanding the genetic code is an important aim of synthetic biology, but some organisms developed naturally expanded genetic codes long ago over the course of evolution. Less than 1% of all sequenced genomes encode an operon that reassigns the stop codon UAG to pyrrolysine (Pyl), a genetic code variant that results from the biosynthesis of Pyl-tRNA(Pyl). To understand the selective advantage of genetically encoding more than 20 amino acids, we constructed a markerless tRNA(Pyl) deletion strain of Methanosarcina acetivorans (ΔpylT) that cannot decode UAG as Pyl or grow on trimethylamine. Phenotypic defects in the ΔpylT strain were evident in minimal medium containing methanol. Proteomic analyses of wild type (WT) M. acetivorans and ΔpylT cells identified 841 proteins from >7,000 significant peptides detected by MS/MS. Protein production from UAG-containing mRNAs was verified for 19 proteins. Translation of UAG codons was verified by MS/MS for eight proteins, including identification of a Pyl residue in PylB, which catalyzes the first step of Pyl biosynthesis. Deletion of tRNA(Pyl) globally altered the proteome, leading to >300 differentially abundant proteins. Reduction of the genetic code from 21 to 20 amino acids led to significant down-regulation in translation initiation factors, amino acid metabolism, and methanogenesis from methanol, which was offset by a compensatory (100-fold) up-regulation in dimethyl sulfide metabolic enzymes. The data show how a natural proteome adapts to genetic code reduction and indicate that the selective value of an expanded genetic code is related to carbon source range and metabolic efficiency.

Comparative genomics highlights the unique biology of Methanomassiliicoccales, a Thermoplasmatales-related seventh order of methanogenic archaea that encodes pyrrolysine

BACKGROUND

A seventh order of methanogens, the Methanomassiliicoccales, has been identified in diverse anaerobic environments including the gastrointestinal tracts (GIT) of humans and other animals and may contribute significantly to methane emission and global warming. Methanomassiliicoccales are phylogenetically distant from all other orders of methanogens and belong to a large evolutionary branch composed by lineages of non-methanogenic archaea such as Thermoplasmatales, the Deep Hydrothermal Vent Euryarchaeota-2 (DHVE-2, Aciduliprofundum boonei) and the Marine Group-II (MG-II). To better understand this new order and its relationship to other archaea, we manually curated and extensively compared the genome sequences of three Methanomassiliicoccales representatives derived from human GIT microbiota, "Candidatus Methanomethylophilus alvus", "Candidatus Methanomassiliicoccus intestinalis" and Methanomassiliicoccus luminyensis.

RESULTS

Comparative analyses revealed atypical features, such as the scattering of the ribosomal RNA genes in the genome and the absence of eukaryotic-like histone gene otherwise present in most of Euryarchaeota genomes. Previously identified in Thermoplasmatales genomes, these features are presently extended to several completely sequenced genomes of this large evolutionary branch, including MG-II and DHVE2. The three Methanomassiliicoccales genomes share a unique composition of genes involved in energy conservation suggesting an original combination of two main energy conservation processes previously described in other methanogens. They also display substantial differences with each other, such as their codon usage, the nature and origin of their CRISPRs systems and the genes possibly involved in particular environmental adaptations. The genome of M. luminyensis encodes several features to thrive in soil and sediment conditions suggesting its larger environmental distribution than GIT. Conversely, "Ca. M. alvus" and "Ca. M. intestinalis" do not present these features and could be more restricted and specialized on GIT. Prediction of the amber codon usage, either as a termination signal of translation or coding for pyrrolysine revealed contrasted patterns among the three genomes and suggests a different handling of the Pyl-encoding capacity.

CONCLUSIONS

This study represents the first insights into the genomic organization and metabolic traits of the seventh order of methanogens. It suggests contrasted evolutionary history among the three analyzed Methanomassiliicoccales representatives and provides information on conserved characteristics among the overall methanogens and among Thermoplasmata.

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.

Development of β -lactamase as a tool for monitoring conditional gene expression by a tetracycline-riboswitch in Methanosarcina acetivorans

The use of reporter gene fusions to assess cellular processes such as protein targeting and regulation of transcription or translation is established technology in archaeal, bacterial, and eukaryal genetics. Fluorescent proteins or enzymes resulting in chromogenic substrate turnover, like β-galactosidase, have been particularly useful for microscopic and screening purposes. However, application of such methodology is of limited use for strictly anaerobic organisms due to the requirement of molecular oxygen for chromophore formation or color development. We have developed β-lactamase from Escherichia coli (encoded by bla) in conjunction with the chromogenic substrate nitrocefin into a reporter system usable under anaerobic conditions for the methanogenic archaeon Methanosarcina acetivorans. By using a signal peptide of a putative flagellin from M. acetivorans and different catabolic promoters, we could demonstrate growth substrate-dependent secretion of β-lactamase, facilitating its use in colony screening on agar plates. Furthermore, a series of fusions comprised of a constitutive promoter and sequences encoding variants of the synthetic tetracycline-responsive riboswitch (tc-RS) was created to characterize its influence on translation initiation in M. acetivorans. One tc-RS variant resulted in more than 11-fold tetracycline-dependent regulation of bla expression, which is in the range of regulation by naturally occurring riboswitches. Thus, tc-RS fusions represent the first solely cis-active, that is, factor-independent system for controlled gene expression in Archaea.

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.

Carbon source-dependent expansion of the genetic code in bacteria

Despite the fact that the genetic code is known to vary between organisms in rare cases, it is believed that in the lifetime of a single cell the code is stable. We found Acetohalobium arabaticum cells grown on pyruvate genetically encode 20 amino acids, but in the presence of trimethylamine (TMA), A. arabaticum dynamically expands its genetic code to 21 amino acids including pyrrolysine (Pyl). A. arabaticum is the only known organism that modulates the size of its genetic code in response to its environment and energy source. The gene cassette pylTSBCD, required to biosynthesize and genetically encode UAG codons as Pyl, is present in the genomes of 24 anaerobic archaea and bacteria. Unlike archaeal Pyl-decoding organisms that constitutively encode Pyl, we observed that A. arabaticum controls Pyl encoding by down-regulating transcription of the entire Pyl operon under growth conditions lacking TMA, to the point where no detectable Pyl-tRNA(Pyl) is made in vivo. Pyl-decoding archaea adapted to an expanded genetic code by minimizing TAG codon frequency to typically ~5% of ORFs, whereas Pyl-decoding bacteria (~20% of ORFs contain in-frame TAGs) regulate Pyl-tRNA(Pyl) formation and translation of UAG by transcriptional deactivation of genes in the Pyl operon. We further demonstrate that Pyl encoding occurs in a bacterium that naturally encodes the Pyl operon, and identified Pyl residues by mass spectrometry in A. arabaticum proteins including two methylamine methyltransferases.

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.

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.

Selenocysteine, pyrrolysine, and the unique energy metabolism of methanogenic archaea

Methanogenic archaea are a group of strictly anaerobic microorganisms characterized by their strict dependence on the process of methanogenesis for energy conservation. Among the archaea, they are also the only known group synthesizing proteins containing selenocysteine or pyrrolysine. All but one of the known archaeal pyrrolysine-containing and all but two of the confirmed archaeal selenocysteine-containing protein are involved in methanogenesis. Synthesis of these proteins proceeds through suppression of translational stop codons but otherwise the two systems are fundamentally different. This paper highlights these differences and summarizes the recent developments in selenocysteine- and pyrrolysine-related research on archaea and aims to put this knowledge into the context of their unique energy metabolism.

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.

The appearance of pyrrolysine in tRNAHis guanylyltransferase by neutral evolution

tRNA(His) guanylyltransferase (Thg1) post-transcriptionally adds a G (position -1) to the 5'-terminus of tRNA(His). The Methanosarcina acetivorans Thg1 (MaThg1) gene contains an in-frame TAG (amber) codon. Although a UAG codon typically directs translation termination, its presence in Methanosarcina mRNA may lead to pyrrolysine (Pyl) incorporation achieved by Pyl-tRNA(Pyl), the product of pyrrolysyl-tRNA synthetase. Sequencing of the MaThg1 gene and transcript confirmed the amber codon. Translation of MaThg1 mRNA led to a full-length, Pyl-containing, active enzyme as determined by immunoblotting, mass spectrometry, and biochemical analysis. The nature of the inserted amino acid at the position specified by UAG is not critical, as Pyl or Trp insertion yields active MaThg1 variants in M. acetivorans and equal amounts of full-length protein. These data suggest that Pyl insertion is akin to natural suppression and unlike the active stop codon reassignment that is required for selenocysteine insertion. Only three Pyl-containing proteins have been characterized previously, a set of methylamine methyltransferases in which Pyl is assumed to have specifically evolved to be a key active-site constituent. In contrast, Pyl in MaThg1 is a dispensable residue that appears to confer no selective advantage. Phylogenetic analysis suggests that Thg1 is becoming dispensable in the archaea, and furthermore supports the hypothesis that Pyl appeared in MaThg1 as the result of neutral evolution. This indicates that even the most unusual amino acid can play an ordinary role in proteins.

Molecular mechanisms of archaeal RNA polymerase

All cellular life depends on multisubunit RNAPs (RNA polymerases) that are evolutionarily related through the three domains of life. Archaeal RNAPs encompass 12 subunits that contribute in different ways to the assembly and stability of the enzyme, nucleic acid binding, catalysis and specific regulatory interactions with transcription factors. The recent development of methods to reconstitute archaeal RNAP from recombinant materials in conjunction with structural information of multisubunit RNAPs present a potent opportunity to investigate the molecular mechanisms of transcription.

The ether-cleaving methyltransferase system of the strict anaerobe Acetobacterium dehalogenans: analysis and expression of the encoding genes

Anaerobic O-demethylases are inducible multicomponent enzymes which mediate the cleavage of the ether bond of phenyl methyl ethers and the transfer of the methyl group to tetrahydrofolate. The genes of all components (methyltransferases I and II, CP, and activating enzyme [AE]) of the vanillate- and veratrol-O-demethylases of Acetobacterium dehalogenans were sequenced and analyzed. In A. dehalogenans, the genes for methyltransferase I, CP, and methyltransferase II of both O-demethylases are clustered. The single-copy gene for AE is not included in the O-demethylase gene clusters. It was found that AE grouped with COG3894 proteins, the function of which was unknown so far. Genes encoding COG3894 proteins with 20 to 41% amino acid sequence identity with AE are present in numerous genomes of anaerobic microorganisms. Inspection of the domain structure and genetic context of these orthologs predicts that these are also reductive activases for corrinoid enzymes (RACEs), such as carbon monoxide dehydrogenase/acetyl coenzyme A synthases or anaerobic methyltransferases. The genes encoding the O-demethylase components were heterologously expressed with a C-terminal Strep-tag in Escherichia coli, and the recombinant proteins methyltransferase I, CP, and AE were characterized. Gel shift experiments showed that the AE comigrated with the CP. The formation of other protein complexes with the O-demethylase components was not observed under the conditions used. The results point to a strong interaction of the AE with the CP. This is the first report on the functional heterologous expression of acetogenic phenyl methyl ether-cleaving O-demethylases.

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.

Class I and class II lysyl-tRNA synthetase mutants and the genetic encoding of pyrrolysine in Methanosarcina spp

Methanosarcina spp. begin methanogenesis from methylamines with methyltransferases made via the translation of UAG as pyrrolysine. In vitro evidence indicates two possible routes to pyrrolysyl-tRNA(Pyl). PylS ligates pyrrolysine to tRNA(Pyl). Alternatively, class I and class II lysyl-tRNA synthetases (LysRS1 and LysRS2) together form lysyl-tRNA(Pyl), a potential intermediate to pyrrolysyl-tRNA(Pyl). The unusual possession of both LysRS1 and LysRS2 by Methanosarcina spp. may also reflect differences in catalytic properties. Here we assessed the in vivo relevance of these hypotheses. The lysK and mtmB transcripts, encoding LysRS1 and monomethylamine methyltransferase, were detectable in Methanosarcina barkeri during early log growth on trimethylamine, but not methanol. In contrast, lysS transcript encoding LysRS2 was detectable during log phase with either substrate. Methanosarcina acetivorans strains bearing deletions of lysK or lysS grew normally on methanol and methylamines with wild-type levels of monomethylamine methyltransferase and aminoacyl-tRNA(Pyl). The lysK and lysS genes could not replace pylS in a recombinant system employing tRNA(Pyl) for UAG suppression. The results support an association of LysRS1 with growth on methylamine, but not an essential role for LysRS1/LysRS2 in the genetic encoding of pyrrolysine. However, decreased lysyl-tRNA(Lys) in the lysS mutant provides a possible rationale for stable transfer of the bacterial lysS gene to methanoarchaea.

In vivo contextual requirements for UAG translation as pyrrolysine

Pyrrolysine and selenocysteine have infiltrated natural genetic codes via the translation of canonical stop codons. UGA translation as selenocysteine is absolutely dependent on message context. Here we describe the first experimental examination of contextual requirements for UAG translation as pyrrolysine. A hexahistidine-tagged Methanosarcina barkeri mtmB1 gene, encoding monomethylamine methyltransferase MtmB1, was introduced into Methanosarcina acetivorans. Host mtmB expression was minimized by growth on methanol and recombinant mtmB1 products monitored by anti-MtmB and anti-hexahistidine immunoblotting. UAG translation was not compromised, as recombinant MtmB1 was 1% of cellular protein with only trace UAG-terminated mtmB1 product detectable. Untranslated regions flanking mtmB1 were not required for UAG translation, but loss of a downstream pyrrolysine insertion sequence (PYLIS) significantly increased the UAG-termination product of mtmB1 and decreased the UAG-translation product, which nonetheless contained pyrrolysine. An in-frame UAG within a bacterial uidA transcript was translated in the methanogen as pyrrolysine with 20% efficiency, suggesting UAG translation in the absence of evolved context. However, predominant UAG-directed termination with enhancement of UAG translation by the PYLIS appears analogous to cis-acting elements for UGA translation as selenocysteine, although different mechanisms may underlie these recoding events.

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

Certain archaeal species can fix molecular nitrogen under nitrogen limiting conditions although little is known about this process at either the genetic or molecular level. To address this on a genome-wide scale, transcriptional analysis was performed on the model methanogen Methanosarcina mazei strain Gö1 using DNA-microarrays. The genomic expression patterns for cells grown under nitrogen fixing conditions versus nitrogen sufficiency (10 mM ammonium) revealed that approximately 5% of all genes are differentially expressed. Besides a small set of genes previously known to be up-regulated under nitrogen limitation, 14 additional genes involved in nitrogen metabolism were identified plus 10 genes encoding potential transcriptional regulators, 13 genes involved in carbon metabolism, 3 genes in general stress response, 8 putative transporter genes, and an additional 21 genes with unknown function. Quantitative reverse transcriptase PCR experiments confirmed the differential expression of a subset of these genes. Promoter analysis revealed a palindromic DNA motif centered nearby the transcriptional start point for several genes up-regulated under nitrogen limitation. A bioinformatics study demonstrated the presence of this motif in the up-stream region of 52 genes genome-wide, the majority of which showed nitrogen dependent differential transcription. We therefore hypothesize that this DNA element is involved in nitrogen control in M. mazei where it may act as a binding site for a regulatory protein.