Gene redundancy in the M. acetivorans genome
The M. acetivorans genome contains many seemingly redundant copies of genes annotated with roles in methanogenesis [5]. These include two sets of genes annotated for a molybdate-type formylmethanofuran dehydrogenase (fmd), and two gene sets for a tunsten-type formylmethanofuran dehydrogenase (fwd), five heterodisulfide reductase-like hdrED and hdrABC gene clusters for reduction of Coenzyme M-Coenzyme B heterodisulfide, two sets of vht genes for F420 non-reducing hydrogenase, and two sets of genes for ATP synthesizing complexes [5]. Additional genes include frh hydrogenase-like genes, plus additional genes for rnf- and mrp-type membrane associated bacterial electron transfer complexes, plus genes needed for acetate metabolism (discussed below). Homologous and seemingly "redundant" genes/gene sets are also found in the genomes of M. mazei, and M. barkeri (Table 1). The reason for these genome makeups is currently unknown. M. acetivorans was used as a model microorganism to evaluate expression of over twenty sets of genes using gene specific primer pairs designed to eliminate cross-hybridization when DNA sequence similarity exists (Methods). RT-PCR, pPCR, and 5' analysis was then performed using RNA isolated from M. acetivorans cells grown with either acetate or methanol as the sole source of carbon and energy. In this study, a number of new M. acetivorans gene designations were established to distinguish among homologous orfs (Table 1, and described below).
Formylmethanofuran dehydrogenase (fmd, fwd) gene expression
Two of the four previously annotated sets of genes for formylmethanofuran dehydrogenasethese were designated as molybdenum-type enzymes and are named here as fmdE1F1A1C1D1B1 and fmdF2A2C2D2B2 (Figure 1A). Two additional gene sets were annotated as tungsten-type formylmethanofuran dehydrogenase, and are designated here as fwdD1B1A1C1 and fwdG2B2D2 (Figure 1B). Using qPCR analysis methods (Methods), the molybdenum-type operon reporter genes fmdE1 and fmdA1 (Figure 1A) were shown to be expressed at 14-fold higher levels during methanol growth conditions relative to acetate growth (Figure 1C). The second set of reporter genes (fmdF2, fmdA2, and fmdB2) were expressed about 2-fold higher during these conditions, but the maximal level of expression was less than 5% of that seen for the fmdE1 and fmdA1 genes. Noteworthy, the fmdE1 and fmdA1 gene expression values were within the same range observed for the fpoN and fpoL genes that encode subunits of the F420 H2 dehydrogenase needed for central pathway electron transfer functions (described below). The high transcript abundance of the fmdE1F1A1C1D1B1 gene cluster implies a major role of this gene set during methanogenesis in contrast to that for the fmd2 gene set.
The annotated tungsten containing formylmethanofuran dehydrogenase gene cluster fwdD1B1A1C1 reporter genes designated fwdB1 and fwdA1 (Figure 1B) were also expressed 15-fold higher levels during methanol growth relative to acetate (Figure 1C). Interestingly, this was within the magnitude observed for the fmdE1F1A1C1D1B1 gene cluster. However, the second tungsten-type gene cluster (as reported by the fwdB2 gene), was constitutively expressed and at a level about one-half of that observed for either fwdA1 or fwdB1. These fmd/fwd transcript abundance measurements clearly demonstrate that two of the four fmd and fwd gene clusters (i.e., fmdE1F1A1C1D1B1 and fwdD1B1A1C1) are highly transcribed in response to substrate availability, and furthermore this suggests that two distinct formylmethanofuran dehydrogenase activities are concurrently utilized during methanol growth conditions (discussed below).
Heterodisulfide reductase gene expression
M. acetivorans genome analysis revealed five genes/gene clusters annotated as heterodisulfide reductase, an enzyme essential for electron transfer from methanogenic electron donors to methyl-CoM reductase (Table 1, Figure 2A). These include genes for a membrane-type protein designated here as hdrE1, hdrD1 and hdrD2 similar to those needed for methane formation in M. barkeri[7]. An additional six genes encoding soluble-type heterodisulfide reductase proteins are also present in the genome. They include the hdrA1 gene associated with a poly-ferredoxin-like gene (pfd), an unlinked set of hdrCB genes called hdrC1 and hdrB1, and a third hdr gene cluster designated hdrA2 hdrC2 hdrB2 (Figure 2B).
Quantitative gene expression experiments (Figure 2C) revealed that the membrane-type hdrD1 gene was most highly expressed during acetate cell growth conditions, and where methanol conditions gave slightly lower transcript abundance (ca. 0.7-fold). In contrast, hdrD2 gene expression was very low (i.e., at level of about one twentieth that seen for the hdrD1 gene Figure 2C), suggesting a minor or no direct function in methanogenesis. Interestingly, the abundance of the soluble-type hdrA1pfd and hdrC1B1 gene transcripts were also nearly as high as for the membrane type hdrED1 genes (Figure 2C). Here, acetate growth gave three-fold higher hdrA1 transcript levels versus methanol growth conditions. The participation of a soluble-type hdrABC enzyme in M. acetivorans metabolism is currently unknown but must now be considered. An orf following the hdrA1 gene is annotated as a polyferredoxin (pfd), and this suggests a role for this protein in electron transfer to couple the soluble-type Hdr complex with an appropriate electron donor complex. In contrast, hdrA2 and hdrB2 transcript abundance was about two to twenty-fold lower under the corresponding conditions. This suggests a minor role for the second set of HdrABC-type genes (i.e., hdrA2B2C2) in methanogenesis.
The hdrA1pfd and hdrC1B1 genes for the soluble-type enzyme subunits are located at different chromosomal loci, and are coordinately expressed since their mRNA abundance levels are alike (Figure 2C). Additionally, the PCR-based gene experiments also demonstrate that the hdrA1pfd and the hdrED1 genes are each expressed as operons (data not shown). Taken together, these data are consistent with a need for both a membrane-type and a soluble type Hdr enzyme for electron transfer/energy conservation under acetate and methanol cell growth conditions. This suggests that distinct electron transfer pathways are operating to service the alternative Hdr enzymes.
The vht and frh gene clusters
The M. acetivorans genome lacks an echABCDEF gene cluster encoding an Ech-type hydrogenase with described roles in hydrogen uptake and ion translocation in M. mazei[3, 5]. Since M. acetivorans cells do not exhibit significant hydrogenase activity [8, 9], some other mechanism must provide a means for electron transfer from cellular donor(s) to Hdr. Interestingly, the M. acetivorans genome contains two sets of genes (designated vhtG1A1C1D1, and vhtG2A2C2) for F420-nonreducing hydrogenase-types (Figure 3A, 3B, Table 1). It also contains a set of frhADGB genes for a coenzyme F420-type hydrogenase (Figure 3A). Quantitative RT-PCR assays (Figure 3C) established that the vhtG1 and vhtC1 genes were each expressed at four- to six-fold higher levels during methanol growth conditions, and this is within the range seen for the fpoL and fpoN genes needed for methyl group oxidation for methanol and acetate metabolism. In contrast, expression of the vhtG2 and vhtC2 genes was low under all conditions examined (Ca. about 17-20-fold lower than vhtG1 and vhtC1). Finally, the frhA and frhB gene expression levels were low relative to vhtG1 or fpoL (Figure 3C), and this suggests a minor role for the frhADGB and vhtG2A2C2 gene clusters in either methanol or acetate-dependent cell growth. Since vhtG1 transcript abundance was elevated and about half of that observed for the fpoL and fpoN genes that encode subunits of the F420 H2 dehydrogenase (Figure 3C), this implies a significant physiological role for the vhtG1A1C1D1 gene products during methanol growth. The biochemical and physiological roles for the vhtG1 and vhtC1 hydrogenase-type genes in M. acetivorans are presently unknown.
The rnfXCDGEABY gene cluster is abundantly expressed
M. acetivorans contains a set of six genes (MA0659-0664) annotated as nqr123456[5] that are absent in the M. mazei, and M. barkeri genomes (Table 1). These genes were subsequently re-designated rnfCDGEAB based on sequence comparisons to the rnf and nqr-type genes in other microorganisms, [10]. This gene cluster also contains two additional genes of unknown function that we designate here as rnfX and rnfY (Figure 4A) whereby the first (MA0658) precedes rnfC and the second (MA0665) follows rnfB. We propose that these genes may encode unique input/output modules for membrane associated electron transfer since they are absent in other microbial genomes. During acetate cell growth relative to methanol growth conditions, the rnfX, rnfG, and rnfA reporter genes exhibited elevated transcript abundance (ca. 2.5 to 3.5-fold; Figure 4D). Each gene was also more highly expressed than many reference genes involved in central methanogenesis (e.g., fpoN, and fpoL that encode subunits of the F420 H2 dehydrogenase). Therefore, the rnfXCDGEABY gene expression data support the proposal that the products participate in electron transfer during acetate metabolism as proposed via methanophenazine [10]. In addition, they must also function during methanol culture conditions based on transcript abundance (Figure 4D). Other roles can be envisioned including participation in electron transfer to a soluble-type heterodisulfide reductase via a poly-ferredoxin (e.g., encoded by the hdrA1 pfd and hdrC1B1 gene complex, described below).
The mrpABCDEFG gene cluster is acetate induced
The M. acetivorans genome contains a set of seven genes called mrpABCDEFG (Figure 4B) with similarity to the gene clusters found in a variety of bacterial species but absent in either M. barkeri or M. mazei (Table 1) [5, 11–13]. The mrp- encoded protein complex in Bacillus subtilis was shown to confer a role in multiple drug resistance and/or pH regulation [13, 14]. As revealed by the M. acetivorans transcript analysis studies (Figure 4D), the mrpA and mrpF reporter genes were expressed more highly during acetate cell growth conditions (Ca. 11 to 12-fold) relative to methanol growth. These levels were above the expression levels observed for the ack, pta, and hdr genes needed for acetate utilization, and within the range seen for the rnf gene cluster. These findings imply a major role for the six mrp gene products in acetate metabolism versus methanol metabolism.
Expression of the atp and aha genes encoding ATP synthase complexes
M. acetivorans contains genes for a bacterial-type F0F1 synthase encoded by the MA2441 to MA2433 genes designated here as atpDCIHBEFAG, plus an archaeal-type A0A1 ATP synthase encoded by the ahaHIKECFABD genes (MA4152 to MA4160) (Figure 5). Although prior DNA microarray experiments [6] demonstrated that six of the nine genes in the archaeal-type A0A1 ATP synthase (ahaECFABD) encoding the ATP-hydrolysing/synthesizing domain (A1) were expressed two-fold higher in acetate grown cells relative to methanol, the other genes were not [6]. It is still unknown how their expression varies quantitatively relative to atpDCIHBEFAG gene cluster expression. Corresponding DNA microarray studies with the atpDCIHBEFAG genes that encode a bacterial-like F0F1 complex revealed that only two of the nine genes (atpD and atpC) were expressed significantly higher in acetate by 3.2 and 1.8 fold, respectively: the remaining genes were either not detected or did not exhibit changes. Lastly, relative to central pathway genes for acetate and methanol utilization, it was unresolved how the aha and atp gene sets are expressed since the microarray data did not address this.
From the RT-PCT transcript abundance studies, three representative aha genes representing the archaeal-type A0A1 ATP synthase genes were highly expressed relative to the atp reporter genes (Figure 5C). Acetate cell growth conditions resulted in two-fold higher aha transcript levels relative to methanol cell growth. These genes were the most highly expressed in the cell regardless of the growth condition. In contrast, the bacterial-type F0F1atpD, atpA and atpG genes were expressed at less than 2% of the level seen for the ahaI, ahaC and ahaB genes: this suggests a minor role for the atp genes in methanogenesis in contrast to the aha gene cluster.
Acetate-induced genes
One M. acetivorans gene of unknown function (MA4008) was revealed by our prior DNA microarray studies to be more highly expressed during acetate growth conditions relative to methanol cell growth (Lars Rohlin, personal communication). Inspection of the amino acid sequence revealed six trans-membrane spanning regions reminiscent of a membrane solute uptake system (Additional file 1, Figure S1 and discussed below). To extend these MA4008 gene expression findings, quantitative PCR experiments were performed (Methods, Figure 6A). MA4008 was expressed at a 125-fold higher level during acetate versus methanol cell growth conditions. Interestingly, when methanol was also present in the culture medium in addition to acetate, MA4008 expression was suppressed to a level seen when only methanol was present (ca. by 215 fold). This indicates that the MA4008 gene is expressed only when the energetically superior carbon substrate is absent, consistent with a proposed role in acetate uptake. The M. acetivorans MA4008 orf is designated aceP for its role in an acetate-dependent membrane function. Two other genes required for acetate utilization are ack (MA3606) and pta (MA3607) that encode acetate kinase and phosphoacetyl transferase, respectively ([15] Table 1). Quantitative PCR experiments (Figure 6A) established that both genes were highly expressed and at levels similar to aceP when acetate was the sole substrate. The 11-18-fold differential pta and ack gene expression findings are similar to previous reports in M. acetivorans and M. thermophila[6, 16].
Location of the fpoP, hdrE, hdrA1, mrpA, pta, aceP, and ahaA promoters
The mRNA 5' ends of the fpoPABCDHIJJKLMNO, hdrED1, hdrA1-pfd, mrpABCDEFG, pta ack, aceP and ahaHIKECFABD genes/clusters were determined to locate their corresponding promoter elements. Using primer extension methods (Figure 7A), all but one of the promoter elements were demonstrated to have long un-translated regions (UTR's) that range from 51 to 137 nucleotides in length. For example, the aceP 5' mRNA end is located 104 nucleotides upstream of the translational start site. Similar findings were seen for the mrpA, fpoP, ahaH, hdrE, and hdrA genes. Only the pta gene had a relatively short UTR (i.e., 27 nt). We did not detect mRNA 5' ends for either rnfX or hdrC1. Alignment of all the upstream regions of these promoter elements (Figure 7A) revealed the highly conserved sequence present in other archaeal promoters, the TATA box (Figure 7B) located approximately 20-30 nt upstream of the +1 mRNA start site (discussed below). This site is bound by the TBP protein that aids RNA polymerase binding [17]. In contrast, the BRE box elements were not well conserved. When the UTR elements and the upstream regions were further examined using a suite of bioinformatics tools (Materials), no clearly discernable DNA sequence elements with either dyad symmetry or direct repeats were found. Similarly, no conserved regions within the RNA UTR's were seen for the coordinately expressed hdrA1pfd and hdrC1B1 genes sets.