mRNAs in the methanogenic archaeon Methanococcus vannielii: numbers, half-lives and processing

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

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

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

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

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

Methyl-coenzyme M reductase A as an indicator to estimate methane production from dairy cows

The evaluation of greenhouse gas mitigation strategies requires the quantitative assessment of individual methane production. Because methane measurement in respiration chambers is highly accurate, but also comprises various disadvantages such as limited capacity and high costs, the establishment of an indicator for estimating methane production of individual ruminants would provide an alternative to direct methane measurement. Methyl-coenzyme M reductase is involved in methanogenesis and the subunit α of methyl-coenzyme M reductase is encoded by the mcrA gene of rumen archaea. We therefore examined the relationship between methane emissions of Holstein dairy cows measured in respiration chambers with 2 different diets (high- and medium-concentrate diet) and the mcrA DNA and mcrA cDNA abundance determined from corresponding rumen fluid samples. Whole-body methane production per kilogram of dry matter intake and mcrA DNA normalized to the abundance of the rrs gene coding for 16S rRNA correlated significantly when using qmcrA primers. Use of qmcrA primers also revealed linear correlation between mcrA DNA copy number and methane yield. Regression analyses based on normalized mcrA cDNA abundances revealed no significant linear correlation with methane production per kilogram of dry matter intake. Furthermore, the correlations between normalized mcrA DNA abundance and the rumen fluid concentration of acetic and isobutyric acid were positive, whereas the correlations with propionic and lactic acid were negative. These data suggest that the mcrA DNA approach based on qmcrA primers could potentially be a molecular proxy for methane yield after further refinement.

Process diagnosis using methanogenic Archaea in maize-fed, trace element depleted fermenters

A mesophilic maize-fed pilot-scale fermenter was severely acidified due to trace element (TE) deficiency. Mainly cobalt (0.07 mg * kg(-1) fresh mass (FM)), selenium (0.007 mg * kg(-1) FM) and sodium (13 mg * kg(-1) FM) were depleted. From this inoculum, three lab-scale flow-through fermenters were operated to analyse micronutrient deficiencies and population dynamics in more detail. One fermenter was supplemented with selenium, one with cobalt, and one served as control. After starvation and recovery of the fermenters, the organic loading rate (OLR) was increased. In parallel, the concentration (Real-Time PCR) of methanogens and their population composition (amplicon sequencing) was determined at the DNA and mRNA level. The parameters Metabolic Quotient (MQ) and cDNA/DNA were calculated to assess the activity of the methanogens. The control without TE supplementation acidified first at an OLR of 4.0 kg volatile solids (VS) * m(-3) * d(-1) while the singular addition of selenium and of cobalt positively influenced the fermenter stability up to an OLR of 4.5 or 5.0 kg VS * m(-3) * d(-1), respectively. In the stable process, the methanogenic populations were dominated by probably residual hydrogenotrophic Methanoculleus sp. (DNA-level), but representatives of versatile Methanosarcina sp. were most active (cDNA-level). When the TE supplemented fermenters began to acidify, Methanosarcina spp. were dominant in the whole (DNA-level) and the active (cDNA-level) community. The acidified control fermenter was dominated by Methanobacteriaceae genus IV. Until acidification, the concentration of methanogens increased with higher OLRs. The MQ indicated stress metabolism approximately one month before the TVA/TIC ratio reached a critical level of 0.7, demonstrating its suitability as early warning parameter of process acidification. The development of the cDNA/DNA ratio also reflected the increasing methanogenic activity with higher OLRs. Highest cDNA/DNA values (ca. 2) were obtained at metabolic strain of the methanogens, at the onset of acidification.

Mechanism for stabilizing mRNAs involved in methanol-dependent methanogenesis of cold-adaptive Methanosarcina mazei zm-15

Methylotrophic methanogenesis predominates at low temperatures in the cold Zoige wetland in Tibet. To elucidate the basis of cold-adapted methanogenesis in these habitats, Methanosarcina mazei zm-15 was isolated, and the molecular basis of its cold activity was studied. For this strain, aceticlastic methanogenesis was reduced 7.7-fold during growth at 15°C versus 30°C. Methanol-derived methanogenesis decreased only 3-fold under the same conditions, suggesting that it is more cold adaptive. Reverse transcription-quantitative PCR (RT-qPCR) detected <2-fold difference in the transcript abundances of mtaA1, mtaB1, and mtaC1, the methanol methyltransferase (Mta) genes, in 30°C versus 15°C culture, while ackA and pta mRNAs, encoding acetate kinase (Ack) and phosphotransacetylase (Pta) in aceticlastic methanogenesis, were 4.5- and 6.8-fold higher in 30°C culture than in 15°C culture. The in vivo half-lives of mtaA1 and mtaC1B1 mRNAs were similar in 30°C and 15°C cultures. However, the pta-ackA mRNA half-life was significantly reduced in 15°C culture compared to 30°C culture. Using circularized RNA RT-PCR, large 5' untranslated regions (UTRs) (270 nucleotides [nt] and 238 nt) were identified for mtaA1 and mtaC1B1 mRNAs, while only a 27-nt 5' UTR was present in the pta-ackA transcript. Removal of the 5' UTRs significantly reduced the in vitro half-lives of mtaA1 and mtaC1B1 mRNAs. Remarkably, fusion of the mtaA1 or mtaC1B1 5' UTRs to pta-ackA mRNA increased its in vitro half-life at both 30°C and 15°C. These results demonstrate that the large 5' UTRs significantly enhance the stability of the mRNAs involved in methanol-derived methanogenesis in the cold-adaptive M. mazei zm-15.

Ecology and characteristics of methanogenic archaea in animals and humans

In this review, the molecular techniques used in animal-based-methanogen studies will be discussed along with how methanogens interact not only with other microorganisms but with their animal hosts as well. These methods not only indicate the diversity and levels of methanogens, but also provide insight on their ecological functions. Most molecular techniques have been based on either 16S rRNA genes or methyl-coenzyme M reductase, a ubiquitous enzyme in methanogens. The most predominant methanogens in animals belong to the genus Methanobrevibacter. Besides methanogens contributing to overall H2 balance, methanogens also have mutual interactions with other bacteria. In addition to shared metabolic synergism, the host animal retrieves additional energy from the diet when methanogens are co-colonized with other normal flora. By comparing genes in methanogens with other bacteria, possible gene transfer between methanogens and other bacteria in the same environments appears to occur. Finally, diets in conjunction with the genetics of methanogens and hosts may represent the biological framework that dictate the extent of methanogen prevalence in these ecosystems. In addition, host evolution including the immune system could serve as an additional selective pressure for methanogen colonization.

Methyl fluoride affects methanogenesis rather than community composition of methanogenic archaea in a rice field soil

The metabolic pathways of methane formation vary with environmental conditions, but whether this can also be linked to changes in the active archaeal community structure remains uncertain. Here, we show that the suppression of aceticlastic methanogenesis by methyl fluoride (CH₃F) caused surprisingly little differences in community composition of active methanogenic archaea from a rice field soil. By measuring the natural abundances of carbon isotopes we found that the effective dose for a 90% inhibition of aceticlastic methanogenesis in anoxic paddy soil incubations was <0.75% CH₃F (v/v). The construction of clone libraries as well as t-RFLP analysis revealed that the active community, as indicated by mcrA transcripts (encoding the α subunit of methyl-coenzyme M reductase, a key enzyme for methanogenesis), remained stable over a wide range of CH₃F concentrations and represented only a subset of the methanogenic community. More precisely, Methanocellaceae were of minor importance, but Methanosarcinaceae dominated the active population, even when CH₃F inhibition only allowed for aceticlastic methanogenesis. In addition, we detected mcrA gene fragments of a so far unrecognised phylogenetic cluster. Transcription of this phylotype at methyl fluoride concentrations suppressing aceticlastic methanogenesis suggests that the respective organisms perform hydrogenotrophic methanogenesis. Hence, the application of CH₃F combined with transcript analysis is not only a useful tool to measure and assign in situ acetate usage, but also to explore substrate usage by as yet uncultivated methanogens.

Genetic manipulation of Methanosarcina spp

The discovery of the third domain of life, the Archaea, is one of the most exciting findings of the last century. These remarkable prokaryotes are well known for their adaptations to extreme environments; however, Archaea have also conquered moderate environments. Many of the archaeal biochemical processes, such as methane production, are unique in nature and therefore of great scientific interest. Although formerly restricted to biochemical and physiological studies, sophisticated systems for genetic manipulation have been developed during the last two decades for methanogenic archaea, halophilic archaea and thermophilic, sulfur-metabolizing archaea. The availability of these tools has allowed for more complete studies of archaeal physiology and metabolism and most importantly provides the basis for the investigation of gene expression, regulation and function. In this review we provide an overview of methods for genetic manipulation of Methanosarcina spp., a group of methanogenic archaea that are key players in the global carbon cycle and which can be found in a variety of anaerobic environments.

Environmental evidence for net methane production and oxidation in putative ANaerobic MEthanotrophic (ANME) archaea

Uncultured ANaerobic MEthanotrophic (ANME) archaea are often assumed to be obligate methanotrophs that are incapable of net methanogenesis, and are therefore used as proxies for anaerobic methane oxidation in many environments in spite of uncertainty regarding their metabolic capabilities. Anaerobic methane oxidation regulates methane emissions in marine sediments and appears to occur through a reversal of a methane-producing metabolism. We tested the assumption that ANME are obligate methanotrophs by detecting and quantifying gene transcription of ANME-1 across zones of methane oxidation versus methane production in sediments from the White Oak River estuary, North Carolina. ANME-1 consistently transcribe 16S rRNA and mRNA of methyl coenzyme M reductase (mcrA), the key gene for methanogenesis, up to 45 cm into methanogenic sediments. CARD-FISH shows that ANME-1 exist as single rod-shaped cells or pairs of cells. Integrating normalized depth distributions of 16S rDNA and rRNA (measured with qPCR and RT-qPCR respectively) shows that 26-77% of the rDNA (a proxy for ANME-1 cell numbers), and 18-76% of the rRNA (a proxy for ANME-1 activity) occurs within methane-producing sediments. These results, along with a re-assessment of the published Iiterature, change the perspective to ANME-1 as methanogens that are also capable of methane oxidation.

Links between methane flux and transcriptional activities of methanogens and methane oxidizers in a blanket peat bog

The relationship between biogeochemical process rates and microbial functional activity was investigated by analysis of the transcriptional dynamics of the key functional genes for methanogenesis (methyl coenzyme M reductase; mcrA) and methane oxidation (particulate methane monooxygenase; pmoA) and in situ methane flux at two peat soil field sites with contrasting net methane-emitting and -oxidizing characteristics. qPCR was used to quantify the abundances of mcrA and pmoA genes and transcripts at two soil depths. Total methanogen and methanotroph transcriptional dynamics, calculated from mcrA and pmoA gene : transcript abundance ratios, were similar at both sites and depths. However, a linear relationship was demonstrated between surface mcrA and pmoA transcript dynamics and surface flux rates at the methane-emitting and methane-oxidizing sites, respectively. Results indicate that methanotroph activity was at least partially substrate-limited at the methane-emitting site and by other factors at the methane-oxidizing site. Soil depth also contributed to the control of surface methane fluxes, but to a lesser extent. Small differences in the soil water content may have contributed to differences in methanogen and methanotroph activities. This study therefore provides a first insight into the regulation of in situ, field-level surface CH(4) flux at the molecular level by an accurate reflection of gene : transcript abundance ratios for the key genes in methane generation and consumption.

Spatial structure and activity of sedimentary microbial communities underlying a Beggiatoa spp. mat in a Gulf of Mexico hydrocarbon seep

BACKGROUND

Subsurface fluids from deep-sea hydrocarbon seeps undergo methane- and sulfur-cycling microbial transformations near the sediment surface. Hydrocarbon seep habitats are naturally patchy, with a mosaic of active seep sediments and non-seep sediments. Microbial community shifts and changing activity patterns on small spatial scales from seep to non-seep sediment remain to be examined in a comprehensive habitat study.

METHODOLOGY/PRINCIPAL FINDINGS

We conducted a transect of biogeochemical measurements and gene expression related to methane- and sulfur-cycling at different sediment depths across a broad Beggiatoa spp. mat at Mississippi Canyon 118 (MC118) in the Gulf of Mexico. High process rates within the mat ( approximately 400 cm and approximately 10 cm from the mat's edge) contrasted with sharply diminished activity at approximately 50 cm outside the mat, as shown by sulfate and methane concentration profiles, radiotracer rates of sulfate reduction and methane oxidation, and stable carbon isotopes. Likewise, 16S ribosomal rRNA, dsrAB (dissimilatory sulfite reductase) and mcrA (methyl coenzyme M reductase) mRNA transcripts of sulfate-reducing bacteria (Desulfobacteraceae and Desulfobulbaceae) and methane-cycling archaea (ANME-1 and ANME-2) were prevalent at the sediment surface under the mat and at its edge. Outside the mat at the surface, 16S rRNA sequences indicated mostly aerobes commonly found in seawater. The seep-related communities persisted at 12-20 cm depth inside and outside the mat. 16S rRNA transcripts and V6-tags reveal that bacterial and archaeal diversity underneath the mat are similar to each other, in contrast to oxic or microoxic habitats that have higher bacterial diversity.

CONCLUSIONS/SIGNIFICANCE

The visual patchiness of microbial mats reflects sharp discontinuities in microbial community structure and activity over sub-meter spatial scales; these discontinuities have to be taken into account in geochemical and microbiological inventories of seep environments. In contrast, 12-20 cm deep in the sediments microbial communities performing methane-cycling and sulfate reduction persist at lower metabolic rates regardless of mat cover, and may increase activity rapidly when subsurface flow changes.

Correlation of methane production and functional gene transcriptional activity in a peat soil

The transcription dynamics of subunit A of the key gene in methanogenesis (methyl coenzyme M reductase; mcrA) was studied to evaluate the relationship between process rate (methanogenesis) and gene transcription dynamics in a peat soil ecosystem. Soil methanogen process rates were determined during incubation of peat slurries at temperatures from 4 to 37 degrees C, and real-time quantitative PCR was applied to quantify the abundances of mcrA genes and transcripts; corresponding transcriptional dynamics were calculated from mcrA transcript/gene ratios. Internal standards suggested unbiased recovery of mRNA abundances in comparison to DNA levels. In comparison to those in pure-culture studies, mcrA transcript/gene ratios indicated underestimation by 1 order of magnitude, possibly due to high proportions of inactive or dead methanogens. Methane production rates were temperature dependent, with maxima at 25 degrees C, but changes in abundance and transcription of the mcrA gene showed no correlation with temperature. However, mcrA transcript/gene ratios correlated weakly (regression coefficient = 0.76) with rates of methanogenesis. Methanogen process rates increased over 3 orders of magnitude, while the corresponding maximum transcript/gene ratio increase was only 18-fold. mcrA transcript dynamics suggested steady-state expression in peat soil after incubation for 24 and 48 h, similar to that in stationary-phase cultures. mcrA transcript/gene ratios are therefore potential in situ indicators of methanogen process rate changes in complex soil systems.

The Sac10b homolog in Methanococcus maripaludis binds DNA at specific sites

The Sac10b protein family, also known as Alba, is widely distributed in Archaea. Sac10b homologs in thermophilic Sulfolobus species are very abundant. They bind both DNA and RNA with high affinity and without sequence specificity, and their physiological functions are still not fully understood. Mma10b from the euryarchaeote Methanococcus maripaludis is a mesophilic member of the Sac10b family. Mma10b is not abundant and constitutes only approximately 0.01% of the total cellular protein. Disruption of mma10b resulted in poor growth of the mutant in minimal medium at near the optimal growth temperature but had no detectable effect on growth in rich medium. Quantitative proteomics, real time reverse transcription-PCR, and enzyme assays revealed that the expression levels of some genes involved in CO(2) assimilation and other activities were changed in the Deltamma10b mutant. Chromatin immunoprecipitation suggested a direct association of Mma10b with an 18-bp DNA binding motif in vivo. Electrophoretic mobility shift assays and DNase I footprinting confirmed that Mma10b preferentially binds specific sequences of DNA with an apparent Kd in the 100 nM range. These results suggested that the physiological role of Mma10b in the mesophilic methanococci is greatly diverged from that of homologs in thermophiles.

Global analysis of mRNA stability in the archaeon Sulfolobus

BACKGROUND

Transcript half-lives differ between organisms, and between groups of genes within the same organism. The mechanisms underlying these differences are not clear, nor are the biochemical properties that determine the stability of a transcript. To address these issues, genome-wide mRNA decay studies have been conducted in eukaryotes and bacteria. In contrast, relatively little is known about RNA stability in the third domain of life, Archaea. Here, we present a microarray-based analysis of mRNA half-lives in the hyperthermophilic crenarchaea Sulfolobus solfataricus and Sulfolobus acidocaldarius, constituting the first genome-wide study of RNA decay in archaea.

RESULTS

The two transcriptomes displayed similar half-life distributions, with medians of about five minutes. Growth-related genes, such as those involved in transcription, translation and energy production, were over-represented among unstable transcripts, whereas uncharacterized genes were over-represented among the most stable. Half-life was negatively correlated with transcript abundance and, unlike the situation in other organisms, also negatively correlated with transcript length.

CONCLUSION

The mRNA half-life distribution of Sulfolobus species is similar to those of much faster growing bacteria, contrasting with the earlier observation that median mRNA half-life is proportional to the minimal length of the cell cycle. Instead, short half-lives may be a general feature of prokaryotic transcriptomes, possibly related to the absence of a nucleus and/or more limited post-transcriptional regulatory mechanisms. The pattern of growth-related transcripts being among the least stable in Sulfolobus may also indicate that the short half-lives reflect a necessity to rapidly reprogram gene expression upon sudden changes in environmental conditions.

Visualization of mcr mRNA in a methanogen by fluorescence in situ hybridization with an oligonucleotide probe and two-pass tyramide signal amplification (two-pass TSA–FISH)

Two-pass tyramide signal amplification-fluorescence in situ hybridization (two-pass TSA-FISH) with a horseradish peroxidase (HRP)-labeled oligonucleotide probe was applied to detect prokaryotic mRNA. In this study, mRNA of a key enzyme for methanogenesis, methyl coenzyme M reductase (mcr), in Methanococcus vannielii was targeted. Applicability of mRNA-targeted probes to in situ hybridization was verified by Clone-FISH. It was observed that sensitivity of two-pass TSA-FISH was significantly higher than that of TSA-FISH, which was further increased by the addition of dextran sulphate in TSA working solution. Signals from two-pass TSA-FISH were more reliable compared to the weak, spotty signals yielded by TSA-FISH.

Translation initiation with GUC codon in the archaeon Halobacterium salinarum: implications for translation of leaderless mRNA and strict correlation between translation initiation and presence of mRNA

We have investigated whether anticodon sequence mutant of an archaeal initiator tRNA can initiate protein synthesis using reporter genes carrying mutations in the initiation codon. Halobacterium salinarum was used as the model organism and the bacterio-opsin gene (bop), which encodes the precursor of the protein component of the purple membrane protein bacterio-opsin (Bop), was chosen as the reporter. We demonstrate that a CAU to GAC anticodon sequence mutant of Haloferax volcanii initiator tRNA can initiate Bop protein synthesis using GUC as the initiation codon in H. salinarum. We generated four mutant bop genes, each carrying the AUG to GUC initiation codon mutation, with or without a compensatory mutation to maintain a predicted stem-loop structure at the 5'-end of the bop mRNA, and with or without mutations to test translation initiation at a site corresponding to the amino terminus of mature bacterio-opsin. H. salinarum chromosomal recombinants containing these mutant genes were phenotypically Pum- (purple membrane negative). Upon transformation with a plasmid carrying the mutant initiator tRNA gene, only strains designed to maintain the bop mRNA stem-loop structure produced Bop and were phenotypically Pum+ as indicated by purple colony colour, and immunoblotting and spectral analysis of cell extracts. Thus GUC can serve as an initiation codon in archaea and the stem-loop structure in the bop mRNA is important for translation. Interestingly, for the same mutant mRNA, only transformants that produce Bop protein contain bop mRNA. These results suggest either a strong coupling between translation and mRNA stability or strong transcriptional polarity in H. salinarum.

Methanol-dependent gene expression demonstrates that methyl-coenzyme M reductase is essential in Methanosarcina acetivorans C2A and allows isolation of mutants with defects in regulation of the methanol utilization pathway

Methanosarcina acetivorans C2A is able to convert several substrates to methane via at least four distinct methanogenic pathways. A common step in each of these pathways is the reduction of methyl-coenzyme M (CoM) to methane catalyzed by methyl-CoM reductase (MCR). Because this enzyme is used in each of the known pathways, the mcrBDCGA operon, which encodes MCR, is expected to be essential. To validate this prediction, a system for conditional gene inactivation was developed. A heterologous copy of the mcrBDCGA operon was placed under the control of the highly regulated mtaC1 promoter, which directs the expression of genes involved in methanol utilization, and recombined onto the M. acetivorans chromosome. This allowed for disruption of the endogenous mcr operon in the presence of methanol. Because the PmtaC1 promoter is transcribed only during growth on methanol, mcrBDCGA was rendered methanol dependent and the strain was unable to grow in trimethylamine media, strongly suggesting that mcrBDCGA is essential. Upon prolonged incubation, suppressed mutants which expressed mcrBDCGA constitutively could be selected. Expression analysis of PmtaC1::uidA gene fusions in several isolated suppressed mutants suggests that they carry trans-active mutations leading to deregulation of all genes under control of this promoter. Subsequently, proteome analysis of one such suppressed mutant revealed that all known proteins derived from mtaC1 promoter-dependent expression were constitutively expressed in this mutant. This genetic system can therefore be employed for the testing of essential genes and for the identification of genes under a common regulatory mechanism by making regulatory mutations phenotypically selectable.

Recoding in archaea

Standard decoding of the genetic information into polypeptides is performed by one of the most sophisticated cell machineries, the translating ribosome, which, by following the genetic code, ensures the correspondence between the mature mRNA and the protein sequence. However, the expression of a minority of genes requires programmed deviations from the standard decoding rules, globally named recoding. This includes ribosome programmed -/+1 frameshifting, ribosome hopping, and stop codon readthrough. Recoding in Archaea was unequivocally demonstrated only for the translation of the UGA stop codon into the amino acid selenocysteine. However, a new recoding event leading to the 22nd amino acid pyrrolysine and the preliminary reports on a gene regulated by programmed -1 frameshifting have been recently described in Archaea. Therefore, it appears that the study of this phenomenon in Archaea is still at its dawn and that most of the genes whose expression is regulated by recoding are still uncharacterized.

Individual gvp transcript segments in Haloferax mediterranei exhibit varying half-lives, which are differentially affected by salt concentration and growth phase

The mc-gvp genes for gas vesicle formation in Haloferax mediterranei are transcribed from two promoters located in front of the mc-gvpA and mc-gvpD genes. The different transcripts originating from both promoters show different abundances dependent on salt concentration in the medium and growth phase. Here we show that the half-lives of these transcripts differ significantly and that the small gvp transcripts exhibit higher stabilities than the larger gvp transcripts. While the stability of most gvp transcripts is independent of the salt concentration in the medium, the gvpA mRNA decays about twice as fast in cultures grown at 18% salt compared to cultures grown at 25% salt. The stability of the 0.45 kb transcript population derived from the 5' part of the gvpD gene depends on the growth phase of the culture. Thus, differences in mRNA stability contribute to the salt-dependent and growth phase-dependent abundance of gvp transcripts. This implies that, like in bacteria and eukarya, mRNA processing contributes to regulated gene expression in archaea.

A comprehensive study into the molecular methodology and molecular biology of methanogenic Archaea

Methanogens belong to the kingdom of Euryarchaeota in the domain of Archaea. The Archaea differ from Bacteria in many aspects important to molecular work. Among these are cell wall composition, their sensitivity to antibiotics, their translation and transcription machinery, and their very strict demands to anaerobic culture conditions. These differences may, at least partly, be responsible for the delay in availability of genetic research tools for methanogens. At present, however, the research within genetics of methanogens and their gene regulation and expression is in rapid progress. Two complete methanogenic genomes have been sequenced and published and more are underway. Besides, sequences are known from a multitude of individual genes from methanogens. Standard methods for simple DNA and RNA work can normally be employed, but permeabilization of the cell wall may demand special procedures. Efficient genetic manipulation systems, including shuttle and integration vector systems, have appeared for mesophilic, but not for thermophilic species within the last few years and will have a major impact on future investigations of methanogenic molecular biology.

Genetic organization of the citCDEF locus and identification of mae and clyR genes from Leuconostoc mesenteroides

In this paper, we describe two open reading frames coding for a NAD-dependent malic enzyme (mae) and a putative regulatory protein (clyR) found in the upstream region of citCDEFG of Leuconostoc mesenteroides subsp. cremoris 195. The transcriptional analysis of the citrate lyase locus revealed one polycistronic mRNA covering the mae and citCDEF genes. This transcript was detected only on RNA prepared from cells grown in the presence of citrate. Primer extension experiments suggest that clyR and the citrate lyase operon are expressed from a bidirectional A-T-rich promoter region located between mae and clyR.

Coordinate transcriptional control in the hyperthermophilic archaeon Sulfolobus solfataricus

The existence of a global gene regulatory system in the hyperthermophilic archaeon Sulfolobus solfataricus is described. The system is responsive to carbon source quality and acts at the level of transcription to coordinate synthesis of three physically unlinked glycosyl hydrolases implicated in carbohydrate utilization. The specific activities of three enzymes, an alpha-glucosidase (malA), a beta-glycosidase (lacS), and an alpha-amylase, were reduced 4-, 20-, and 10-fold, respectively, in response to the addition of supplementary carbon sources to a minimal sucrose medium. Western blot analysis using anti-alpha-glucosidase and anti-beta-glycosidase antibodies indicated that reduced enzyme activities resulted exclusively from decreased enzyme levels. Northern blot analysis of malA and lacS mRNAs revealed that changes in enzyme abundance arose primarily from reductions in transcript concentrations. Culture conditions precipitating rapid changes in lacS gene expression were established to determine the response time of the regulatory system in vivo. Full induction occurred within a single generation whereas full repression occurred more slowly, requiring nearly 38 generations. Since lacS mRNA abundance changed much more rapidly in response to a nutrient down shift than to a nutrient up shift, transcript synthesis rather than degradation likely plays a role in the regulatory response.

Clustered genes encoding the methyltransferases of methanogenesis from monomethylamine

Coenzyme M (CoM) is methylated during methanogenesis from monomethyamine in a reaction catalyzed by three proteins. Using monomethylamine, a 52-kDa polypeptide termed monomethylamine methyltransferase (MMAMT) methylates the corrinoid cofactor bound to a second polypeptide, monomethylamine corrinoid protein (MMCP). Methylated MMCP then serves as a substrate for MT2-A, which methylates CoM. The genes for these proteins are clustered on 6.8 kb of DNA in Methanosarcina barkeri MS. The gene encoding MMCP (mtmC) is located directly upstream of the gene encoding MMAMT (mtmB). The gene encoding MT2-A (mtbA) was found 1.1 kb upstream of mtmC, but no obvious open reading frame was found in the intergenic region between mtbA and mtmC. A single monocistronic transcript was found for mtbA that initiated 76 bp from the translational start. Separate transcripts of 2.4 and 4.7 kb were detected, both of which carried mtmCB. The larger transcript also encoded mtmP, which is homologous to the APC family of cationic amine permeases and may therefore encode a methylamine permease. A single transcriptional start site was found 447 bp upstream of the translational start of mtmC. MtmC possesses the corrinoid binding motif found in corrinoid proteins involved in dimethylsulfide- and methanol-dependent methanogenesis, as well as in methionine synthase. The open reading frame of mtmB was interrupted by a single in-frame, midframe, UAG codon which was also found in mtmB from M. barkeri NIH. A mechanism that circumvents UAG-directed termination of translation must operate during expression of mtmB in this methanogen.

Multiple transcriptional control of the Lactococcus lactis trp operon

The Lactococcus lactis trpEGDCFBA operon is preceded by a noncoding leader region. Transcriptional studies of the trp operon revealed three transcripts with respective sizes of 8 kb (encompassing the entire operon), 290 bases, and 160 bases (corresponding to parts of the leader region). These transcripts most likely result from initiation at the unique Ptrp promoter, transcription termination at either T1 (upstream of the trp operon) or T2 (downstream of the trp operon), and/or processing. Three parameters were shown to differentially affect the amount of these transcripts: (i) following tryptophan depletion, the amount of the 8-kb transcript increases 300- to 500-fold; (ii) depletion in any amino acid increased transcription initiation about fourfold; and (iii) upon entry into stationary phase the amount of the 8-kb transcript decreases abruptly. The tryptophan-dependent transcription control is exerted through transcription antitermination.

Biochemical and serological evidence for an RNase E-like activity in halophilic Archaea

Endoribonuclease RNase E appears to control the rate-limiting step that mediates the degradation of many mRNA species in bacteria. In this work, an RNase E-like activity in Archaea is described. An endoribonucleolytic activity from the extreme halophile Haloarcula marismortui showed the same RNA substrate specificity as the Escherichia coli RNase E and cross-reacted with a monoclonal antibody raised against E. coli RNase E. The archaeal RNase E activity was partially purified from the extreme halophilic cells and shown, contrary to the E. coli enzyme, to require a high salt concentration for cleavage specificity and stability. These data indicate that a halophilic RNA processing enzyme can specifically recognize and cleave mRNA from E. coli in an extremely salty environment (3 M KCI). Having recently been shown in mammalian cells (A. Wennborg, B. Sohlberg, D. Angerer, G. Klein, and A. von Gabain, Proc. Natl. Acad. Sci. USA 92:7322-7326, 1995), RNase E-like activity has now been identified in all three evolutionary domains: Archaea, Bacteria, and Eukarya. This strongly suggests that mRNA decay mechanisms are highly conserved despite quite different environmental conditions.

Sequence and transcript analysis of a novel Methanosarcina barkeri methyltransferase II homolog and its associated corrinoid protein homologous to methionine synthase

The sequence and transcript of the genes encoding a recently discovered coenzyme M methylase in Methanosarcina barkeri were analyzed. This 480-kDa protein is composed of two subunits in equimolar concentrations which bind one corrinoid cofactor per alphabeta dimer. The gene for the alphabeta polypeptide, mtsA, is upstream of that encoding the beta polypeptide, mtsB. The two genes are contiguous and overlap by several nucleotides. A 1.9-kb mRNA species which reacted with probes specific for either mtsA or mtsB was detected. Three possible methanogen consensus BoxA sequences as well as two sets of direct repeats were found upstream of mtsA. The 5' end of the mts transcript was 19 nucleotides upstream of the translational start site of mtsA and was positioned 25 bp from the center of the proximal BoxA sequence. The transcript was most abundant in cells grown to the late log phase on acetate but barely detectable in cells grown on methanol or trimethylamine. The amino acid sequence of MtsB was homologous to the cobalamin-binding fragment of methionine synthase from Escherichia coli and possessed the signature residues involved in binding the corrinoid, including a histidyl residue which ligates cobalt. The sequence of MtsA is homologous to the "A" and "M" isozymes of methylcobamide:coenzyme M methyltransferases (methyltransferase II), indicating that the alpha polypeptide is a new member of the methyltransferase II family of coenzyme M methylases. All three methyltransferase II homolog sequences could be aligned with the sequences of uroporphyrinogen decarboxylase from various sources. The implications of these homologies for the mechanism of corrinoid binding by proteins involved in methylotrophic methanogenesis are discussed.

Carbon monoxide dehydrogenase from Methanosarcina frisia Gö1. Characterization of the enzyme and the regulated expression of two operon-like cdh gene clusters

Carbon monoxide dehydrogenase (Cdh) has been anaerobically purified from Methanosarcina frisia Gö1. The enzyme is a Ni2+-, Fe2+-, and S2--containing alpha2beta2 heterotetramer of 214 kDa with a pI of 5.2 and subunits of 94 and 19 kDa. It has a Vmax of 0.3 mmol of CO min-1 mg-1 and Km values for CO and methyl viologen of approximately 0.9 mM and 0.12 mM, respectively. EPR spectroscopy on the reduced enzyme showed two overlapping signals: one indicative for 2 (4Fe-4S)+ clusters and a second signal that is atypical for standard Fe/S clusters. The latter was, together with high-spin EPR signals of the oxidized enzyme tentatively assigned to an Fe/S cluster of high nuclearity.

Cloning, functional organization, transcript studies, and phylogenetic analysis of the complete nitrogenase structural genes (nifHDK2) and associated genes in the archaeon Methanosarcina barkeri 227

Determination of the nucleotide sequence of the nitrogenase structural genes (nifHDK2) from Methanosarcina barkeri 227 was completed in this study by cloning and sequencing a 2.7-kb BamHI fragment containing the 3' end of nifK2 and 1,390 bp of the nifE2-homologous genes. Open reading frame nifK2 is 1,371 bp long including the stop codon TAA and encodes a polypeptide of 456 amino acids. Phylogenetic analysis of the deduced amino acid sequences of the nifK2 and nifE2 gene products from M. barkeri showed that both genes cluster most closely with the corresponding nif-1 gene products from Clostridium pasteurianum, consistent with our previous analyses of nifH2 and nifD2. The nifE gene product is known to be homologous to that of nifD, and our analysis shows that the branching pattern for the nifE proteins resembles that for the nifD product (with the exception of vnfE from Azotobacter vinelandii), suggesting that a gene duplication occurred before the divergence of nitrogenases. Primer extension showed that nifH2 had a single transcription start site located 34 nucleotides upstream of the ATG translation start site for nifH2, and a sequence resembling the archaeal consensus promoter sequence [TTTA(A/T)ATA] was found 32 nucleotides upstream from that transcription start site. A tract of four T's, previously identified as a transcription termination site in archaea, was found immediately downstream of the nifK2 gene, and a potential promoter was located upstream of the nifE2 gene. Hybridization with nifH2 and nifDK2 probes with M. barkeri RNA revealed a 4.6-kb transcript from N2-grown cells, large enough to harbor nifHDK genes and their internal open reading frames, while no transcript was detected from NH4(+)-grown cells. These results support a model in which the nitrogenase structural genes in M. barkeri are cotranscribed in a single NH4(+)-repressed operon.

Purification of the copper response extracellular proteins secreted by the copper-resistant methanogen Methanobacterium bryantii BKYH and cloning, sequencing, and transcription of the gene encoding these proteins

When the copper-resistant methanogen Methanobacterium bryantii BKYH was exposed to 1 mM Cu(II), it secreted approximately fourfold increased levels of three proteins, copper response extracellular (CRX) proteins. The members of the CRX protein trio had apparent molecular masses of 40.8, 42.3, and 42.9 kDa and were purified together from the culture supernatant and separated from each other by sodium dodecyl sulfate-polyacrylamide gel electrophoresis. The N-terminal amino acid sequences of the three proteins were essentially identical, and antibodies raised against one of the trio reacted with all three proteins and with three other intracellular proteins with slightly higher molecular weights. The N-terminal amino acid sequence of one of these larger proteins was different from that of the secreted CRX proteins. The gene crx, which encodes the CRX proteins, was cloned and sequenced, and crx transcription was characterized. The crx sequence predicts that the encoded polypeptide is synthesized as a precursor with an N-terminal leader peptide, containing 28 amino acid residues, that is removed during the extracellular secretion of the CRX proteins. Transcription was initiated 274 bp upstream from the crx gene, producing an approximately 1.4-kb monocistronic transcript that was present in M. bryantii BKYH cells under all growth conditions but that increased approximately fourfold in vivo in response to Cu addition. The CRX proteins appear to be glycosylated, since they react with concanavalin A and neuraminidase, and to be the products of one gene that have different levels of posttranslational glycosylation. This is supported by very similar chromatographic and electrophoretic properties, identical N-terminal amino acid sequences, immunological cross-reactivities, and the detection of only one crx-related sequence by Southern blotting. Western blots (immunoblots) showed no evidence for CRX proteins in cell lysates of several other Methanobacterium strains.

Organization and growth phase-dependent transcription of methane genes in two regions of the Methanobacterium thermoautotrophicum genome

Two regions of the Methanobacterium thermoautotrophicum genome containing genes that encode enzymes involved in methanogenesis (methane genes) have been cloned and sequenced to determine the extent of methane gene clustering and conservation. One region from the M. thermoautotrophicum strains delta H and Winter, extending approximately 13.5 kb upstream from the adjacent mvhDGAB and mrtBDGA operons that encode the methyl-viologen-reducing hydrogenase (MVH) and the methyl coenzyme M reductase II (MRII), respectively, was sequenced, and 76% sequence identity and very similar gene organizations were demonstrated. Five closely linked open reading frames were located immediately upstream of the mvh operon and were designated flpECBDA. The flpCBD genes encode amino acid sequences that are 31, 47, and 65% identical to the primary sequences of the alpha and beta subunits of formate dehydrogenase and the delta subunit of MVH, respectively. Located immediately upstream of the flp genes was the mth gene, which encodes the H2-dependent methylene-tetrahydromethanopterin dehydrogenase (MTH). In contrast to this mth-flp-mvh-mrt cluster of methane genes, a separate approximately 5.4-kb genomic fragment cloned from M. thermoautotrophicum delta H contained only one methane gene, the mtd gene, which encodes the 8-hydroxy-5-deazaflavin (H2F420)-dependent methylene-tetrahydromethanopterin dehydrogenase (MTD). Northern (RNA) blot experiments demonstrated that mth was transcribed only at early growth stages in fermentor-grown cultures of M. thermoautotrophicum delta H, whereas mtd was transcribed at later growth stages and in the stationary phase. Very similar transcription patterns have been observed by T.D. Pihl, S. Sharma, and J. N. Reeve (J. Bacteriol. 176:6384-6391, 1994) for the MRI- and MRII-encoding operons, mrtBDGA and mcrBDCGA, im M. thermoautotrophicum deltaH, suggesting coordinated regulation of methane gene expression. In contrast to the growth phase-dependent transcription of the mth/mrt and mtd/mcr genes, transcription of the mvhDGAB and frhADGB operons, which encode the two (NiFe) hydrogenases in M. thermoautotrophicum deltaH, was found to occur at all growth stages.

The energy conserving N5-methyltetrahydromethanopterin:coenzyme M methyltransferase complex from Methanobacterium thermoautotrophicum is composed of eight different subunits

N5-Methyltetrahydromethanopterin:coenzyme M methyltransferase (Mtr) from Methanobacterium thermoautotrophicum strain Marburg is a membrane-associated enzyme complex which catalyzes an energy-conserving, sodium-ion-translocating step in methanogenesis from H2 and CO2. We report here that the complex is composed of eight different subunits for which evidence was obtained at the protein, DNA and RNA levels: (a) SDS/PAGE of the purified complex revealed the presence of eight different polypeptides of apparent molecular masses of 34 (MtrH), 28 (MtrE), 24 (MtrC), 23 (MtrA), 21 (MtrD), 13 (MtrG), 12.5 (MtrB) and 12 kDa (MtrF). The N-terminal amino acid sequences of the 12-, 12.5- and 13-kDa polypeptides, which had previously not been accessible, were determined; (b) cloning and sequencing of the corresponding genes revealed the presence of the eight mtr genes organized in a 4.9-kbp gene cluster in the order mtrEDCBAFGH; (c) Northern-blot analysis revealed the presence of a 5-kbp transcript. DNA probes derived from the mtrE and mtrH genes hybridized to the transcript, indicating that the eight mtr genes are organized in a transcription unit. By primer extension, the 5' end of the mtrEDC-BAFGH mRNA was analyzed. The mtr operon was found to be located between the methyl-coenzyme M reductase I operon (mcr) and a downstream open reading frame predicted to encode a Na+/Ca2+, K+ exchanger.

Cloning, sequencing, and transcriptional analysis of the coenzyme F420-dependent methylene-5,6,7,8-tetrahydromethanopterin dehydrogenase gene from Methanobacterium thermoautotrophicum strain Marburg and functional expression in Escherichia coli

Two methylenetetrahydromethanopterin dehydrogenases have been purified from Methanobacterium thermoautotrophicum strain Marburg: one (MTD) is coenzyme F420-dependent and oxygen-stable (Mukhopadhyay, B., and Daniels, L. (1989) Can. J. Microbiol. 35, 499-507), and the other (MTH) is coenzyme F420-independent (or hydrogenase-type) and oxygen-sensitive (Zirngibl, C., Hedderich, R., and Thauer, R. K. (1990) FEBS Lett. 261, 112-116). Based on the NH2-terminal sequence of MTD, a 36-mer oligonucleotide was designed and used to identify and clone a 6.1-kilobase pair EcoRI fragment of M. thermoautotrophicum DNA. Sequencing of this fragment revealed an 825-base pair (bp) MTD encoding gene (mtd), which was expressed in Escherichia coli yielding an enzyme that, like the native enzyme, was oxygen-stable, strictly dependent on coenzyme F420, thermostable, thermophilic, and exhibited maximum activity at an acidic pH. The amino acid sequence predicts that MTD is a hydrophobic and acidic protein with no identifiable homology to MTH (von Bunau, R., Zirngibl, C., Thauer, R. K., and Klein, A. (1991) Eur. J. Biochem. 202, 1205-1208), but comparisons with coenzyme F420 utilizing enzymes revealed a conserved region at the NH2 terminus of MTD that could correspond to the ability to interact with coenzyme F420. The mtd transcript was approximately 900 nucleotides long and initiated 8 bp upstream of the translation initiation codon and 22 bp downstream from an archaeal promoter sequence. The mtd coding sequence was followed by several poly(dT) sequences and an inverted repeat that could be transcription termination signals.

Growth phase-dependent transcription of the genes that encode the two methyl coenzyme M reductase isoenzymes and N5-methyltetrahydromethanopterin:coenzyme M methyltransferase in Methanobacterium thermoautotrophicum delta H

The genes encoding the two isoenzymes of methyl coenzyme M reductase (MRI and MRII) in Methanobacterium thermoautotrophicum delta H have been cloned and sequenced. The MRI-encoding mcr operon (mcrBDCGA) has been located immediately upstream from the mtr operon (mtrEDCBA) that encodes N5-methyltetrahydromethanopterin:coenzyme M methyltransferase, the enzyme that catalyzes the step preceding the MR-catalyzed reaction in methanogenesis. The MRII-encoding mrt operon (mrtBDGA) has been located between the operon that encodes the methyl viologen-reducing hydrogenase and an open reading frame (designated pyrC) predicted to encode dihydroorotase. Surprisingly, the mrt operon has been found to contain only four genes (mrtBDGA), lacking the equivalent of the mcrC gene that is present in all mcr operons. A protocol that isolates transcripts intact from M. thermoautotrophicum delta H cells has been developed and used, with primer extension and Northern (RNA) blot procedures, to identify the sites of transcription initiation upstream of the mcr, mrt, and mtr operons and to determine the relative numbers of these transcripts in cells at different growth stages. Transcription of the mrt operon was found to occur only at early times in batch cultures and was then replaced by transcription of the mcr operon. Transcripts of the mtr operon were detectable at all times; however, at early times, all mtr transcripts were initiated at the mtr promoter, whereas at later times, during mcr transcription, approximately 3% of mcr transcripts were extended to generate mcr plus mtr transcripts that constituted approximately 20% of all mtr transcripts present.