Zahra F. Islam
Paul R. F. Cordero
Chris Greening- School of Biological Sciences, Monash University, Clayton, VIC, Australia
Aerobic soil bacteria persist by scavenging molecular hydrogen (H2) from the atmosphere. This key process is the primary sink in the biogeochemical hydrogen cycle and supports the productivity of oligotrophic ecosystems. In Mycobacterium smegmatis, atmospheric H2 oxidation is catalyzed by two phylogenetically distinct [NiFe]-hydrogenases, Huc (group 2a) and Hhy (group 1h). However, it is currently unresolved how these enzymes transfer electrons derived from H2 oxidation into the aerobic respiratory chain. In this work, we used genetic approaches to confirm that two putative iron-sulfur cluster proteins encoded on the hydrogenase structural operons, HucE and HhyE, are required for H2 consumption in M. smegmatis. Sequence analysis show that these proteins, while homologous, fall into distinct phylogenetic clades and have distinct metal-binding motifs. H2 oxidation was reduced when the genes encoding these proteins were deleted individually and was eliminated when they were deleted in combination. In turn, the growth yield and long-term survival of these deletion strains was modestly but significantly reduced compared to the parent strain. In both biochemical and phenotypic assays, the mutant strains lacking the putative iron-sulfur proteins phenocopied those of hydrogenase structural subunit mutants. We hypothesize that these proteins mediate electron transfer between the catalytic subunits of the hydrogenases and the menaquinone pool of the M. smegmatis respiratory chain; however, other roles (e.g., in maturation) are also plausible and further work is required to resolve their role. The conserved nature of these proteins within most Hhy- or Huc-encoding organisms suggests that these proteins are important determinants of atmospheric H2 oxidation.
Introduction
Over the last decade, various studies have revealed that aerobic bacteria conserve energy during persistence through aerobic respiration of atmospheric hydrogen (H2) (Constant et al., 2010; Greening et al., 2014b, 2015a; Meredith et al., 2014; Liot and Constant, 2016; Islam et al., 2019). This process is now recognized to be important for biogeochemical and ecological reasons. Gas-scavenging soil bacteria serve as the primary sink in the global hydrogen cycle and are responsible for the net consumption of approximately 70 million tonnes of H2 each year (Constant et al., 2009; Ehhalt and Rohrer, 2009; Greening et al., 2014a; Piché-Choquette et al., 2018). More recently, it has been inferred that this process supports the productivity and biodiversity of various ecosystems, especially low-carbon soils (Lynch et al., 2014; Kanno et al., 2015; Khdhiri et al., 2015; Greening et al., 2016; Ji et al., 2017; Bay et al., 2018; Kessler et al., 2019; Piché-Choquette and Constant, 2019). Atmospheric H2 oxidation appears to be a widespread trait among soil bacteria. To date, bacteria from three phyla have been experimentally shown to oxidize atmospheric H2, Actinobacteriota (Constant et al., 2008, 2010; Greening et al., 2014a; Meredith et al., 2014), Acidobacteriota (Greening et al., 2015a; Myers and King, 2016), and Chloroflexota (Islam et al., 2019). However, genomic and metagenomic studies have indicated at least 13 other phyla encode enzymes that can mediate this process (Greening et al., 2016; Carere et al., 2017; Ji et al., 2017; Piché-Choquette et al., 2017).
The genetic basis and physiological role of atmospheric H2 oxidation is now largely understood. This process has been most comprehensively studied in the genetically tractable soil actinobacterium Mycobacterium smegmatis (Greening and Cook, 2014). In this organism, atmospheric H2 oxidation is mediated by two membrane-bound, oxygen-tolerant hydrogenases, Huc (group 2a [NiFe]-hydrogenase, also known as Hyd1 or cyanobacterial-type uptake hydrogenase) and Hhy (group 1h [NiFe]-hydrogenase, also known as Hyd2 or actinobacterial-type uptake hydrogenase) (Berney et al., 2014b). Additionally, M. smegmatis encodes a third [NiFe]-hydrogenase, Hyh (Hyd3), which mediates fermentative H2 production during hypoxia (Berney et al., 2014a). Both H2-oxidizing hydrogenases contain a large subunit containing the [NiFe] active site (HucL, HhyL) and a small subunit containing three iron-sulfur clusters (HucS, HhyS), as well as potential additional subunits (Berney et al., 2014b; Cordero et al., 2019b). These two hydrogenases are upregulated in stationary-phase cells, including in response to organic carbon limitation (Berney and Cook, 2010; Berney et al., 2014b). Consistently, when the structural subunits of these hydrogenases are deleted, strains show reduced growth yield and impaired long-term survival during starvation (Berney and Cook, 2010; Berney et al., 2014a; Greening et al., 2014b). Similar findings have been made in Streptomyces avermitilis; the sole hydrogenase of this organism, Hhy, is exclusively expressed in exospores and strains lacking this enzyme exhibit severe survival defects (Constant et al., 2010; Liot and Constant, 2016). Given these findings, it is proposed that bacteria shift from growing on organic compounds to persisting on atmospheric trace gases. Indeed, theoretical calculations indicate that the energy derived from atmospheric H2 oxidation (0.53 ppmv) can sustain the maintenance of 107 to 108 cells per gram of soil (Conrad, 1999).
Despite this progress, little is currently known about the biochemical basis of atmospheric H2 oxidation. One outstanding question is how electrons derived from H2 oxidation are transferred to the respiratory chain. Most classes of respiratory uptake hydrogenases are predicted to be co-transcribed with a cytochrome b subunit (Greening et al., 2016; Søndergaard et al., 2016). For example, such subunits interact with the prototypical oxygen-tolerant hydrogenases (group 1d [NiFe]-hydrogenases) of Escherichia coli and Ralstonia eutropha; they anchor the hydrogenase to the membrane and transfer electrons from the hydrogenase small subunit to the quinone pool (Frielingsdorf et al., 2011; Volbeda et al., 2013). However, we did not detect equivalent proteins in the operons encoding the structural subunits of Huc (MSMEG_2261–2270) or Hhy (MSMEG_2722 – 2718) in M. smegmatis (Supplementary Figure S1) (Berney et al., 2014b). Putative iron-sulfur proteins, tentatively annotated as HucE (MSMEG_2268) and HhyE (MSMEG_2718), were encoded downstream of the hydrogenase structural subunits and may potentially fulfill this role instead (Berney et al., 2014b; Greening et al., 2015b). In this work, we characterized the effects of deleting these genes on hydrogenase activity, growth, and survival in M. smegmatis. We also investigated their broader conservation in hydrogenase-encoding bacteria.
Materials and Methods
Bacterial Strains and Growth Conditions
All bacterial strains and plasmids used in this study are listed in Supplementary Table S1. Escherichia coli TOP10 was maintained on lysogeny broth (LB) agar plates (10 g L–1 tryptone, 5 g L–1 NaCl, 5 g L–1 yeast extract, 15 g L–1 agar), while Mycobacterium smegmatis mc2155 (Snapper et al., 1990) and derived mutants were maintained on LB agar plates supplemented with 0.05% (w/v) Tween 80 (LBT). For broth culture, E. coli was grown in LB. M. smegmatis was grown in either LBT or in Hartmans de Bont (HdB) minimal medium (Hartmans and De Bont, 1992) supplemented with 0.2% (w/v) glycerol. In all cases, liquid cultures were grown in rotary incubators at 37°C with agitation (200 rpm).
Mutant Strain Construction
Allelic exchange mutagenesis was used to produce markerless deletions of the genes encoding two putative iron-sulfur proteins, hucE (MSMEG_2268) and hhyE (MSMEG_2718) (Supplementary Figure S2). Briefly, a fragment containing fused left and right flanks of the MSMEG_2268 (1800 bp) and MSMEG_2718 (3098 bp) genes were synthesized by GenScript. These fragments were cloned into the SpeI site of the mycobacterial shuttle plasmid pX33 (Gebhard et al., 2006) to yield the constructs pX33-hucE and pX33-hhyE (Supplementary Table S1). These constructs were propagated in E. coli TOP10 and transformed into wild-type M. smegmatis mc2155 cells by electroporation. Gentamycin (5 μg mL–1 for M. smegmatis or 20 μg mL–1 for E. coli) was used in selective solid and liquid medium to propagate pX33. Creation of the double iron-sulfur cluster mutant (ΔhucEΔhhyE) was achieved by transformation of ΔhhyE electrocompetent M. smegmatis mc2155 with the pX33-hucE construct. Briefly, to allow for permissive temperature-sensitive vector replication, transformants were incubated on LBT gentamicin plates at 28°C until colonies were visible (5–7 days). Resultant catechol-positive colonies were subcultured onto fresh LBT gentamicin plates and incubated at 40°C for 3–5 days to facilitate integration of the recombinant plasmid flanks into the chromosome. The second recombination event was facilitated by subculturing catechol-reactive and gentamicin-resistant colonies onto LBT agar plates supplemented with 10% sucrose (w/v) and incubating at 40°C for 3–5 days. Catechol-unreactive colonies were subsequently screened by PCR to discern wild-type revertants from ΔhucE, ΔhhyE and ΔhucEΔhhyE mutants. Primers used for the generation of mutants and for screening are listed in Supplementary Table S2.
Complementation Vector Construction
The genes for the putative iron-sulfur proteins were amplified by PCR and the resulting fragments were cloned into the constitutive expression plasmid pMV261 via PstI/HindIII site for hucE and BamHI/HindIII site for hhyE (Stover et al., 1991) to yield the constructs pMVhucE and pMVhhyE (Supplementary Table S1). Sequence fidelity of the genes was verified through Sanger sequencing and insertion of the genes into the vector was confirmed through restriction-digestion analysis (Supplementary Figure S3). The plasmid constructs were propagated in E. coli DH5α and transformed into M. smegmatis cells by electroporation. Vector pMVhucE was transformed into M. smegmatis wild-type and ΔhucE strains, while pMVhhyE was transformed into wild type and ΔhhyE mutant. In addition, an empty pMV261 was transformed into wild-type, ΔhucE, and ΔhhyE strains. These seven M. smegmatis strains were used for complementation experiments in respirometry and activity staining. Kanamycin (20 μg mL–1 for M. smegmatis or 50 μg mL–1 for E. coli) was used in selective solid and liquid medium to propagate pMV261. Primers used for the generation of the constructs are listed in Supplementary Table S2.
Respirometry Measurements
Cultures of wild-type, derived mutants, and complemented mutant strains of M. smegmatis were grown in 125 mL aerated conical flasks containing 30 mL HdB medium supplemented with 0.2% glycerol. Respirometry measurements were performed with mid-stationary phase cells, i.e., 72 h post ODmax (∼3.0). A Unisense H2 microsensor electrode was polarized at + 800 mV for 1 h using a Unisense multimeter and calibrated against standards of known H2 concentration. Gas-saturated PBS was prepared by bubbling the solution with 100% (v/v) of either H2 or O2 for 5 min. The 1.1 mL microrespiration assay chambers were sequentially amended with stationary-phase cultures (0.9 mL, OD600 = 3.0), H2-saturated PBS (0.1 mL), and O2-saturated PBS (0.1 mL). Chambers were stirred at 250 rpm, 37°C. Changes in H2 concentration were recorded using Unisense Logger Software, and upon observing a linear change in H2 concentration, rates of consumption were calculated over a period of 20 s, which corresponds to the most linear uptake of hydrogen by the cells. Oxidation rates were normalized against total protein concentration, which was determined by the bicinchoninic acid method (Smith et al., 1985) with bovine serum albumin standards.
Activity Staining
Cultures of wild-type, derived mutants, and complemented mutant strains of M. smegmatis were grown in 2.5 L aerated conical flasks containing 500 mL HdB medium supplemented with 0.2% glycerol. For Huc activity staining, cultures of wild-type, ΔhucS, ΔhucE, ΔhucSΔhhyL, ΔhucEΔhhyE, and complemented ΔhucE and wild-type M. smegmatis (either with empty pMV261 or complementation vector pMVhucE) were harvested by centrifugation (10,000 × g, 10 min, 4°C) at early-stationary phase (24 h post ODmax, ∼3.0) (Cordero et al., 2019b). For Hhy activity staining, cultures of wild-type, ΔhhyL, ΔhhyE, ΔhucSΔhhyL, ΔhucEΔhhyE, and complemented ΔhhyE and wild-type M. smegmatis (either with empty pMV261 or complementation vector pMVhhyE) were harvested by centrifugation at mid-stationary phase (72 h post ODmax, ∼3.0) (Cordero et al., 2019b). Harvested cultures were washed in phosphate-buffered saline solution (PBS; 137 mM NaCl, 2.7 mM KCl, 10 mM Na2HPO4, 2 mM KH2PO4, pH 7.4), and resuspended in 16 mL lysis buffer (50 mM Tris-Cl, pH 8.0, 1 mM PMSF, 2 mM MgCl2, 5 mg mL–1 lysozyme, 40 μg mL–1 DNase, 10% glycerol). Resultant cell suspensions were passed through a Constant Systems cell disruptor (40,000 psi, four times), with unbroken cells removed by centrifugation (10,000 × g, 20 min, 4°C) to yield whole-cell lysates. Protein concentration was determined using a bicinchoninic acid assay with bovine serum albumin standards. Next, 20 μg of each whole-cell lysate was loaded onto two native 7.5% (w/v) Bis-Tris polyacrylamide gels prepared as described elsewhere (Walker, 2009) and run alongside a protein standard (NativeMark Unstained Protein Standard, Thermo Fisher Scientific) for 1.5 h at 25 mA. One gel was stained overnight at 4°C with gentle agitation using AcquaStain Protein Gel Stain (Bulldog Bio) for total protein determination. The other gel was incubated for hydrogenase activity staining in 50 mM potassium phosphate buffer (pH 7.0) supplemented with 500 μM nitroblue tetrazolium chloride (NBT) in an anaerobic jar amended with an anaerobic gas mixture (5% H2, 10% CO2, 85% N2 v/v) overnight at room temperature.
Growth and Survival Assays
Cultures of wild type and derived mutants of M. smegmatis were inoculated into 125 mL conical flasks containing 30 mL LBT medium (initial OD600 of 0.001), in six biological replicates. Growth was monitored by measuring optical density at 600 nm (1 cm cuvettes; Eppendorf BioSpectrometer Basic); when OD600 was above 0.5, cultures were diluted ten-fold in LBT before measurement. Specific growth rate during mid-exponential growth was calculated for each replicate using GraphPad Prism (non-linear regression, exponential growth equation, least squares fit). The long-term survival of the cultures was determined by counting colony forming units (CFU mL–1) of cultures 21 days post-ODmax. Cultures were serially diluted in HdB containing no carbon source and spotted on to agar plates in technical quadruplicates. After incubation at 37°C for 3 days, the resultant colonies were counted.
Sequence and Phylogenetic Analysis
Sequences homologous to M. smegmatis HucE (MSMEG_2268) and HhyE (MSMEG_2718) were retrieved by protein BLAST (Altschul et al., 1990) using the National Center for Biotechnology Information (NCBI) Reference Sequence (RefSeq) database (Pruitt et al., 2007). The retrieved hits were cross-referenced with the hydrogenase database (HydDB) (Søndergaard et al., 2016) in order to determine which organisms co-encode HucE with HucL and HhyE with HhyL. For downstream phylogenetic and motif analysis, sequences were filtered to remove truncated HucE/HhyE proteins and retain one protein sequence per genus. This resulted in a representative subset of 52 full-length HucE and 26 full-length HhyE sequences. The retrieved sequences were aligned using ClustalW in MEGA7 (Kumar et al., 2016). The phylogenetic relationships of these sequences were visualized on a maximum-likelihood tree based on the Poisson correction method and bootstrapped with 100 replicates. In addition, WebLogo (Crooks et al., 2004) was used to analyze the conserved motifs containing cysteine and histidine residues predicted to bind iron-sulfur clusters. The web-based software Properon (M. Milton1) was used to generate to-scale genetic organization diagrams of the group 1h and group 2a [NiFe]-hydrogenases, with genes labeled according to the nomenclature in HydDB (Søndergaard et al., 2016).
Results
HucE and HhyE Are Predicted to Be Iron-Sulfur Proteins Associated With Group 2a and Group 1h [NiFe]-Hydrogenases
We investigated the diversity of putative iron-sulfur cluster proteins associated with [NiFe]-hydrogenases by conducting a homology-based search using the amino acid sequences of HucE (MSMEG_2268) and HhyE (MSMEG_2718) from M. smegmatis (Supplementary Figure S1). Homologous sequences were retrieved from 14 phyla and 104 genera of bacteria (Supplementary Figure S4 and Supplementary Table S3).
The evolutionary relationships of these proteins were visualized on a maximum-likelihood phylogenetic tree (Figure 1). All retrieved sequences fall into two robustly supported clades, the HucE proteins associated with group 2a [NiFe]-hydrogenases (Huc) and the HhyE proteins associated with group 1h [NiFe]-hydrogenases (Hhy), that share approximately 27% amino acid identity. HhyE proteins were encoded by various atmospheric H2 oxidizers, including Streptomyces (Berney and Cook, 2010), Rhodococcus (Meredith et al., 2014), Pyrinomonas (Greening et al., 2015a), and Thermogemmatispora (Islam et al., 2019). HucE proteins were encoded by various Cyanobacteria, which are known to recycle H2 produced during the nitrogenase reaction via group 2a [NiFe]-hydrogenases (Houchins and Burris, 1981; Tamagnini et al., 2002), as well as genera capable of aerobic hydrogenotrophic growth such as Nitrospira (Koch et al., 2014), Pseudonocardia (Grostern and Alvarez-Cohen, 2013), and Acidithiobacillus (Schröder et al., 2007). Of the hydrogenase-positive species surveyed, 9.5% lacked HucE and HhyE, including Thermomicrobium (Islam et al., 2019) and Methylacidiphilum (Mohammadi et al., 2017) species known to synthesize mid-affinity group 1h [NiFe]-hydrogenases. In contrast, no HucE or HhyE sequences were retrieved from organisms that lack hydrogenases.
Figure 1. Phylogenetic tree of HucE and HhyE proteins associated with group 2a and 1h [NiFe]-hydrogenases. The tree visualizes the evolutionary relationships between a representative subset of 52 full-length HucE and 26 full-length HhyE sequences. The proteins encoded by Mycobacterium smegmatis are emphasized. The tree was constructed using the maximum-likelihood method (gaps treated with partial deletion), bootstrapped with 100 replicates, and rooted at the mid-point. The sequences used to create this tree are provided in Supplementary Table S3.
Multiple sequence alignments show that HucE and HhyE proteins contain highly conserved motifs potentially involved in binding iron-sulfur clusters (Supplementary Figures S5, S6). Both HucE and HhyE contain a CxxC motif within a domain homologous to NifU proteins (Yuvaniyama et al., 2000). The C-terminus of HhyE proteins contains two CxxC motifs typical of iron-sulfur proteins (e.g., rubredoxins). In contrast, the HucE proteins contain an C-terminal motif CxH(x15–18)CxxC that matches the signature motif of Rieske iron-sulfur clusters (Schmidt and Shaw, 2001) (Supplementary Figure S6). A subset of the species surveyed contain truncated HucE and HhyE proteins that contain the NifU-like domain, but lack the C-terminal domains (Supplementary Figures S4, S5).
HucE and HhyE Are Essential for H2 Oxidation in Mycobacterium smegmatis
We used allelic exchange mutagenesis to generate markerless single and double mutants of the hucE and hhyE genes in M. smegmatis, i.e., ΔhucE, ΔhhyE, and ΔhucEΔhhyE. Gene deletion was confirmed by PCR targeting chromosomal sequences adjacent to the flanking regions used for homologous recombination (Supplementary Figure S2). Assays were used to compare H2 oxidation of these strains with the wild-type strain and strains containing previously generated deletions of the hydrogenase structural subunits, i.e., ΔhucS, ΔhhyL, and ΔhucSΔhhyL, that lack hydrogenase activity (Berney and Cook, 2010; Berney et al., 2014b; Greening et al., 2014a).
We first used a H2 electrode to measure rates of aerobic H2 respiration mediated by whole cells of each strain. There were significant differences in the rate of H2 oxidation for all deletion strains compared to the wild type (Figure 2A). Loss of hucE and hhyE resulted in reductions of 1.8-fold and 8.4-fold, respectively; such reductions were statistically indistinguishable from those observed in the mutants of the hydrogenase structural subunits hucS and hhyL. Deletion of both iron-sulfur proteins (ΔhucEΔhhyE) or both hydrogenase structural subunits (ΔhucSΔhhyL) caused complete cessation of H2 oxidation, highlighting that these two hydrogenases are solely responsible for H2 oxidation and that the putative iron-sulfur proteins are indispensable for this process. The low-level negative rates in ΔhucEΔhhyE and ΔhucSΔhhyL strains most likely reflect drift of the electrode rather than actual H2 production by Hyh (Hyd3), since this hydrogenase is only upregulated during hypoxia (Berney et al., 2014a). We successfully complemented the ΔhucE and ΔhhyE strains by reintroducing the hucE and hhyE genes on the episomal plasmid pMV261 (Stover et al., 1991) (Figure 2B); in contrast, introducing the empty vector caused no effect and neither did introducing the complementation vectors in a wild-type background. This restoration of Huc and Hhy activities in complemented iron-sulfur protein deletion mutants strongly indicate that HucE and HhyE are essential for H2 oxidation. Moreover, the similarity in H2 oxidation rates between the strains containing deletions of the catalytic subunits, compared to the putative iron-sulfur proteins, is consistent with HucE and HhyE being functionally linked with the Huc and Hhy hydrogenases, respectively.
Figure 2. Hydrogen oxidation by wild-type, derived mutants, and complemented mutant strains of M. smegmatis. H2 uptake by whole cells in mid-stationary phase (72 h post ODmax ∼3.0) was measured amperometrically using a Unisense H2 electrode. (A) Comparison of the rates of H2 oxidation between wild-type, single and double mutants of the iron-sulfur proteins (ΔhucE, ΔhhyE, ΔhucEΔhhyE), and single and double mutants of hydrogenase structural subunits (ΔhucS, ΔhhyL, ΔhucSΔhhyL). (B) Rates of H2 oxidation in ΔhucE and ΔhhyE strains complemented with expression of hucE and hhyE, respectively. Controls include wild-type, ΔhucE, and ΔhhyE strains transformed with empty vector pMV261 and wild-type strain transformed with complementation vectors pMVhucE and pMVhhyE. Error bars show standard deviations of three biological replicates and values labeled with different letters are significantly different (p < 0.05) based on a one-way ANOVA.
In an interrelated assay, we performed activity staining of the Huc and Hhy hydrogenases using whole-cell lysates of wild-type and deletion mutant strains, with and without the complementation vectors, in the presence of the artificial electron acceptor nitroblue tetrazolium chloride. In the Huc activity staining gel (Figure 3A), three bands were observed in the whole-cell lysates of wild-type strains, with or without complementation vectors: the top high-MW band, middle mid-MW band, and bottom low-MW band. Both the high-MW and low-MW bands correspond to Huc activity (Cordero et al., 2019b) and these bands were not observed in strains lacking either hucS or hucE. However, Huc activity was restored when the ΔhucE strain was complemented by episomal expression of hucE. For Hhy activity staining (Figure 3B), a mid-sized MW band was be observed in all wild-type strains. This band, which is the same middle band observed in the Huc activity stain (Figure 3A), corresponds to Hhy activity (Greening et al., 2014a; Cordero et al., 2019b). No Hhy staining was detected with the loss of either hhyL or hhyE, but complementation of the ΔhhyE strain with hhyE restored Hhy activity. The similarity in the staining bands observed between ΔhucS and ΔhucE strains or between ΔhhyL and ΔhhyE indicate that the putative iron-sulfur proteins HucE and HhyE, like their respective hydrogenase core subunits HucS and HhyL, are important for hydrogenase activity. The artificial electron acceptor cannot compensate for the loss of HucE/HhyE and neither can HucE for HhyE nor HhyE for HucE. This further supports the model that HucE and HhyE form a functional association with Huc and Hhy, respectively.
Figure 3. Hydrogenase activity staining in wild-type, derived mutants, and complemented mutant strains of M. smegmatis. Whole-cell lysates were used for zymographic staining of H2 uptake in a H2-rich atmosphere with nitroblue tetrazolium as artificial electron acceptor. (A) Huc activity staining of cultures of wild-type, ΔhucS, ΔhucE, ΔhucSΔhhyL, ΔhucEΔhhyE, and complemented ΔhucE and wild-type M. smegmatis (either with empty pMV261 or complementation vector pMVhucE) harvested at early-stationary phase (24 h post ODmax ∼3.0). (B) Hhy activity staining of cultures of wild-type, ΔhhyL, ΔhhyE, ΔhucSΔhhyL, ΔhucEΔhhyE, and complemented ΔhhyE and wild-type M. smegmatis (either with empty pMV261 or complementation vector pMVhhyE) harvested at mid-stationary phase (72 h post ODmax ∼3.0). The original gels and Coomassie stain are shown in Supplementary Figure S7.
HucE and HhyE Mutant Strains Have Significant Growth and Survival Defects
Previous genetic studies have shown that the hydrogenases modestly increase growth yield and long-term survival of M. smegmatis (Berney and Cook, 2010; Greening et al., 2014b). We therefore tested whether these findings extended to the putative iron-sulfur proteins by analyzing the growth rate, growth yield, and long-term survival of the seven aforementioned strains when cultured aerobically on rich media (LBT). In line with previous findings (Berney et al., 2014a; Greening et al., 2014b), no significant differences in specific growth rate were observed between the strains (Figure 4A). However, there was a 10% reduction in the specific growth yield of the HhyE mutant compared to the wild-type strain (ODmaxwt = 4.19 ± 0.21; ODmax ΔhhyE = 3.81 ± 0.09; p = 0.008) (Figure 4B). This phenotype extended to the double mutant strain (ΔhucEΔhhyE) and again phenocopied single and double mutants lacking the hhyL gene.
Figure 4. Comparison of growth and survival between wild-type and mutant strains of M. smegmatis. Seven strains were grown on lysogeny broth supplemented with Tween80 (LBT): wild-type, single and double mutants of the iron-sulfur proteins (ΔhucE, ΔhhyE, ΔhucEΔhhyE), and single and double mutants of hydrogenase structural subunits (ΔhucS, ΔhhyL, ΔhucSΔhhyL). (A) Specific growth rate (μ) during exponential phase. (B) Final growth yield (ODmax) at 24 h post-stationary phase. (C) Long-term survival (CFU mL– 1) at 21 days post-stationary phase. Error bars show standard deviations of six biological replicates. Values labeled with different letters are significantly different (p < 0.05) based on a one-way ANOVA.
We also tested whether the strains were defective in long-term survival by counting colonies of aerobic cultures 21 days following ODmax. There were significant reductions in the survival of most strains compared to the wild-type (Figure 4C). Cell counts were approximately two-fold lower for the ΔhhyE and ΔhhyL strains (p < 0.02), and four-fold lower for the double mutant strains (p < 0.002), relative to the wild-type. These findings agree with previous reports that atmospheric H2 oxidation by the hydrogenases enables M. smegmatis to survive energy starvation (Greening et al., 2014b) and further supports that the putative iron-sulfur proteins contribute to this function. For reasons currently unclear, no phenotypes were observed for the ΔhucE strain.
Discussion
In summary, this study shows that HucE and HhyE are required for the enzymatic activity and physiological function of the mycobacterial uptake hydrogenases. Strains lacking these proteins showed no hydrogenase activity in either amperometric or zymographic assays. Furthermore, they exhibited growth and survival phenotypes similar to those of knockouts of hydrogenase structural subunits (Berney and Cook, 2010; Greening et al., 2014b); as with the structural subunit mutants, these phenotypes are relatively minor, likely reflecting the numerous survival mechanisms present in M. smegmatis such as the ability to persist on carbon monoxide (Cordero et al., 2019a). Despite some sequence similarity between the two proteins, they are non-redundant, as there was no compensation in hydrogenase activity in the single mutant strains. The genomic survey and phylogenetic analysis indicate that hucE and hhyE genes co-evolved with the genes encoding the structural subunits of the group 2a and group 1h [NiFe]-hydrogenases. Their detection in the genomes of most but not all characterized high-affinity H2 oxidizers indicate they are important but overlooked mediators of atmospheric H2 oxidation. They are also associated with the group 2a [NiFe]-hydrogenases of H2-recycling Cyanobacteria and various aerobic hydrogenotrophic bacteria that are not currently known to oxidize atmospheric H2.
This study lends some support to the hypothesis that these proteins serve as the immediate electron acceptors for the group 2a and group 1h [NiFe]-hydrogenases. There are broadly five lines of evidence that support this hypothesis: (i) the presence of highly conserved motifs for binding iron-sulfur clusters, (ii) the essentiality of these proteins for the function of these hydrogenases, (iii) their association with the structural rather than maturation operons of the hydrogenases (Berney et al., 2014b), (iv) co-localization of HhyL, HhyS, and HhyE subunits on native polyacrylamide gels (Cordero et al., 2019b), and (v) their genomic association with hydrogenases that lack known electron transfer subunits (e.g., cytochrome b subunits). With the respect to the latter point, it is interesting that these proteins are conserved in Cyanobacteria, given the immediate electron acceptors of their uptake hydrogenases have long remained enigmatic (Tamagnini et al., 2002). It is also notable that HucE proteins encode the signature motifs of a Rieske iron-sulfur cluster. Given their unusual ligands, these clusters have a higher standard redox potential (Eo’ > −150 mV) than most iron-sulfur clusters (e.g., ferredoxins) (Brown et al., 2008). They would therefore be well-poised to accept the relatively high-potential electrons derived from atmospheric H2 and transfer them to menaquinone. Consistently, zymographic studies suggest that the high-affinity hydrogenases operate at higher redox potential than prototypical hydrogenases, given they are reactive with the nitroblue tetrazolium (Eo’ = −80 mV) but not viologen compounds (Eo’ = −360 mV) (Pinske et al., 2012; Greening et al., 2014a).
While this study demonstrates HucE and HhyE are important for mycobacterial hydrogenase activity, further work is ultimately needed to resolve their respective function. While a role in electron transfer is most plausible, we have not demonstrated that these proteins interact with the hydrogenases and it is notable that the artificial electron acceptor nitroblue tetrazolium chloride cannot compensate for their absence. In this regard, other roles are also possible and compatible with the available evidence, for example as specific assembly factors and/or structural scaffolds for the hydrogenases. For example, it has been demonstrated that a rubredoxin-related protein is important for aerobic maturation of the group 1d [NiFe]-hydrogenase in R. eutropha (Fritsch et al., 2014). Furthermore, it is possible that other hypothetical proteins downstream of HucE and HhyE may also serve as electron acceptor candidates, in particular MSMEG_2717 that shares homology to PHG067, the proposed electron acceptor of R. eutropha (Schäfer et al., 2013). Biochemical studies, including studying the redox chemistry of these proteins and their interactions with the as-yet-unpurified hydrogenases, are now required to distinguish these possibilities and develop a sophisticated understanding of their function.
Data Availability Statement
All datasets generated for this study are included in the article/Supplementary Material.
Author Contributions
CG conceived, designed, and supervised the study. CG and ZI were responsible for the phylogenetic analysis and analyzed the data. ZI and PC were responsible for the knockout generation, hydrogen electrode measurements, and activity staining. PC was responsible for the complementation experiments. ZI was responsible for the phenotypic assays. ZI, PC, and CG wrote the manuscript.
Funding
This work was supported by an ARC DECRA Fellowship (DE170100310; awarded to CG), an Australian Government Research Training Program Stipend (awarded to ZI), and a Monash University Doctoral Scholarship (awarded to PC).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank Dr. Rhys Grinter and the three peer reviewers for their helpful feedback.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.02749/full#supplementary-material
Footnotes
References
Altschul, S. F., Gish, W., Miller, W., Myers, E. W., and Lipman, D. J. (1990). Basic local alignment search tool. J. Mol. Biol. 215, 403–410. doi: 10.1016/S0022-2836(05)80360-2
Bay, S., Ferrari, B., and Greening, C. (2018). Life without water: how do bacteria generate biomass in desert ecosystems? Microbiol. Aust. 39, 28–32. doi: 10.1071/MA18008
Berney, M., and Cook, G. M. (2010). Unique flexibility in energy metabolism allows mycobacteria to combat starvation and hypoxia. PLoS One 5:e8614. doi: 10.1371/journal.pone.0008614
Berney, M., Greening, C., Conrad, R., Jacobs, W. R., and Cook, G. M. (2014a). An obligately aerobic soil bacterium activates fermentative hydrogen production to survive reductive stress during hypoxia. Proc. Natl. Acad. Sci. U.S.A. 111, 11479–11484. doi: 10.1073/pnas.1407034111
Berney, M., Greening, C., Hards, K., Collins, D., and Cook, G. M. (2014b). Three different [NiFe] hydrogenases confer metabolic flexibility in the obligate aerobe Mycobacterium smegmatis. Environ. Microbiol. 16, 318–330. doi: 10.1111/1462-2920.12320
Brown, E. N., Friemann, R., Karlsson, A., Parales, J. V., Couture, M. M.-J., Eltis, L. D., et al. (2008). Determining rieske cluster reduction potentials. JBIC J. Biol. Inorg. Chem. 13:1301. doi: 10.1007/s00775-008-0413-4
Carere, C. R., Hards, K., Houghton, K. M., Power, J. F., McDonald, B., Collet, C., et al. (2017). Mixotrophy drives niche expansion of verrucomicrobial methanotrophs. ISME J. 11, 2599–2610. doi: 10.1038/ismej.2017.112
Conrad, R. (1999). Soil microorganisms oxidizing atmospheric trace gases (CH4, CO, H2, NO). Indian J. Microbiol. 39, 193–203. doi: 10.1128/AEM.00275-17
Constant, P., Chowdhury, S. P., Pratscher, J., and Conrad, R. (2010). Streptomycetes contributing to atmospheric molecular hydrogen soil uptake are widespread and encode a putative high-affinity [NiFe]-hydrogenase. Environ. Microbiol. 12, 821–829. doi: 10.1111/j.1462-2920.2009.02130.x
Constant, P., Poissant, L., and Villemur, R. (2008). Isolation of Streptomyces sp. PCB7, the first microorganism demonstrating high-affinity uptake of tropospheric H2. ISME J. 2, 1066–1076. doi: 10.1038/ismej.2008.59
Constant, P., Poissant, L., and Villemur, R. (2009). Tropospheric H2 budget and the response of its soil uptake under the changing environment. Sci. Total Environ. 407, 1809–1823. doi: 10.1016/j.scitotenv.2008.10.064
Cordero, P. R. F., Bayly, K., Leung, P. M., Huang, C., Islam, Z. F., Schittenhelm, R. B., et al. (2019a). Atmospheric carbon monoxide oxidation is a widespread mechanism supporting microbial survival. ISME J. 13, 2868–2881. doi: 10.1038/s41396-019-0479-8
Cordero, P. R. F., Grinter, R., Hards, K., Cryle, M. J., Warr, C. G., Cook, G. M., et al. (2019b). Two uptake hydrogenases differentially interact with the aerobic respiratory chain during mycobacterial growth and persistence. J. Biol. Chem. (in press). doi: 10.1074/jbc.RA119.011076
Crooks, G. E., Hon, G., Chandonia, J.-M., and Brenner, S. E. (2004). WebLogo: a sequence logo generator. Genome Res. 14, 1188–1190. doi: 10.1101/gr.849004
Ehhalt, D. H., and Rohrer, F. (2009). The tropospheric cycle of H2: a critical review. Tellus B 61, 500–535. doi: 10.1111/j.1600-0889.2009.00416.x
Frielingsdorf, S., Schubert, T., Pohlmann, A., Lenz, O., and Friedrich, B. (2011). A trimeric supercomplex of the oxygen-tolerant membrane-bound [NiFe]-hydrogenase from Ralstonia eutropha H16. Biochemistry 50, 10836–10843. doi: 10.1021/bi201594m
Fritsch, J., Siebert, E., Priebe, J., Zebger, I., Lendzian, F., Teutloff, C., et al. (2014). Rubredoxin-related maturation factor guarantees metal cofactor integrity during aerobic biosynthesis of membrane-bound [NiFe] hydrogenase. J. Biol. Chem. 289, 7982–7993. doi: 10.1074/jbc.M113.544668
Gebhard, S., Tran, S. L., and Cook, G. M. (2006). The Phn system of Mycobacterium smegmatis: a second high-affinity ABC-transporter for phosphate. Microbiology 152, 3453–3465. doi: 10.1099/mic.0.29201-0
Greening, C., Berney, M., Hards, K., Cook, G. M., and Conrad, R. (2014a). A soil actinobacterium scavenges atmospheric H2 using two membrane-associated, oxygen-dependent [NiFe] hydrogenases. Proc. Natl. Acad. Sci. U.S.A. 111, 4257–4261. doi: 10.1073/pnas.1320586111
Greening, C., Villas-Bôas, S. G., Robson, J. R., Berney, M., and Cook, G. M. (2014b). The growth and survival of Mycobacterium smegmatis is enhanced by co-metabolism of atmospheric H2. PLoS One 9:e103034. doi: 10.1371/journal.pone.0103034
Greening, C., Biswas, A., Carere, C. R., Jackson, C. J., Taylor, M. C., Stott, M. B., et al. (2016). Genomic and metagenomic surveys of hydrogenase distribution indicate H2 is a widely utilised energy source for microbial growth and survival. ISME J. 10, 761–777. doi: 10.1038/ismej.2015.153
Greening, C., Carere, C. R., Rushton-Green, R., Harold, L. K., Hards, K., Taylor, M. C., et al. (2015a). Persistence of the dominant soil phylum Acidobacteria by trace gas scavenging. Proc. Natl. Acad. Sci. U.S.A. 112, 10497–10502. doi: 10.1073/pnas.1508385112
Greening, C., Constant, P., Hards, K., Morales, S. E., Oakeshott, J. G., Russell, R. J., et al. (2015b). Atmospheric hydrogen scavenging: from enzymes to ecosystems. Appl. Environ. Microbiol. 81, 1190–1199. doi: 10.1128/AEM.03364-14
Greening, C., and Cook, G. M. (2014). Integration of hydrogenase expression and hydrogen sensing in bacterial cell physiology. Curr. Opin. Microbiol. 18, 30–38. doi: 10.1016/j.mib.2014.02.001
Grostern, A., and Alvarez-Cohen, L. (2013). RubisCO-based CO2 fixation and C1 metabolism in the actinobacterium Pseudonocardia dioxanivorans CB1190. Environ. Microbiol. 15, 3040–3053. doi: 10.1111/1462-2920.12144
Hartmans, S., and De Bont, J. A. (1992). Aerobic vinyl chloride metabolism in Mycobacterium aurum L1. Appl. Environ. Microbiol. 58, 1220–1226.
Houchins, J. P., and Burris, R. H. (1981). Comparative characterization of two distinct hydrogenases from Anabaena sp. strain 7120. J. Bacteriol. 146, 215–221.
Islam, Z. F., Cordero, P. R. F., Feng, J., Chen, Y.-J., Bay, S., Gleadow, R. M., et al. (2019). Two Chloroflexi classes independently evolved the ability to persist on atmospheric hydrogen and carbon monoxide. ISME J. 13, 1801–1813. doi: 10.1038/s41396-019-0393-0
Ji, M., Greening, C., Vanwonterghem, I., Carere, C. R., Bay, S. K., Steen, J. A., et al. (2017). Atmospheric trace gases support primary production in Antarctic desert surface soil. Nature 552, 400–403. doi: 10.1038/nature25014
Kanno, M., Constant, P., Tamaki, H., and Kamagata, Y. (2015). Detection and isolation of plant-associated bacteria scavenging atmospheric molecular hydrogen. Environ. Microbiol. 18, 2495–2506. doi: 10.1111/1462-2920.13162
Kessler, A. J., Chen, Y.-J., Waite, D. W., Hutchinson, T., Koh, S., Popa, M. E., et al. (2019). Bacterial fermentation and respiration processes are uncoupled in permeable sediments. Nat. Microbiol. 4, 1014–1023. doi: 10.1038/s41564-019-0391-z
Khdhiri, M., Hesse, L., Popa, M. E., Quiza, L., Lalonde, I., Meredith, L. K., et al. (2015). Soil carbon content and relative abundance of high affinity H2-oxidizing bacteria predict atmospheric H2 soil uptake activity better than soil microbial community composition. Soil Biol. Biochem. 85, 1–9. doi: 10.1016/j.soilbio.2015.02.030
Koch, H., Galushko, A., Albertsen, M., Schintlmeister, A., Gruber-Dorninger, C., Lucker, S., et al. (2014). Growth of nitrite-oxidizing bacteria by aerobic hydrogen oxidation. Science 345, 1052–1054. doi: 10.1126/science.1256985
Kumar, S., Stecher, G., and Tamura, K. (2016). MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874. doi: 10.1093/molbev/msw054
Liot, Q., and Constant, P. (2016). Breathing air to save energy – new insights into the ecophysiological role of high-affinity [NiFe]-hydrogenase in Streptomyces avermitilis. Microbiologyopen 5, 47–59. doi: 10.1002/mbo3.310
Lynch, R. C., Darcy, J. L., Kane, N. C., Nemergut, D. R., and Schmidt, S. K. (2014). Metagenomic evidence for metabolism of trace atmospheric gases by high-elevation desert Actinobacteria. Front. Microbiol. 5:698. doi: 10.3389/fmicb.2014.00698
Meredith, L. K., Rao, D., Bosak, T., Klepac-Ceraj, V., Tada, K. R., Hansel, C. M., et al. (2014). Consumption of atmospheric hydrogen during the life cycle of soil-dwelling actinobacteria. Environ. Microbiol. Rep. 6, 226–238. doi: 10.1111/1758-2229.12116
Mohammadi, S., Pol, A., van Alen, T. A., Jetten, M. S. M., and Op den Camp, H. J. M. (2017). Methylacidiphilum fumariolicum SolV, a thermoacidophilic “Knallgas” methanotroph with both an oxygen-sensitive and -insensitive hydrogenase. ISME J. 11, 945–958. doi: 10.1038/ismej.2016.171
Myers, M. R., and King, G. M. (2016). Isolation and characterization of Acidobacterium ailaaui sp. nov., a novel member of Acidobacteria subdivision 1, from a geothermally heated Hawaiian microbial mat. Int. J. Syst. Evol. Microbiol. 66, 5328–5335. doi: 10.1099/ijsem.0.001516
Piché-Choquette, S., and Constant, P. (2019). Molecular hydrogen, a neglected key driver of soil biogeochemical processes. Appl. Environ. Microbiol. 85:e02418-18. doi: 10.1128/AEM.02418-18
Piché-Choquette, S., Khdhiri, M., and Constant, P. (2017). Survey of high-affinity H2-oxidizing bacteria in soil reveals their vast diversity yet underrepresentation in genomic databases. Microb. Ecol. 74, 771–775. doi: 10.1007/s00248-017-1011-1
Piché-Choquette, S., Khdhiri, M., and Constant, P. (2018). Dose-response relationships between environmentally-relevant H2 concentrations and the biological sinks of H2, CH4 and CO in soil. Soil Biol. Biochem. 123, 190–199. doi: 10.1016/j.soilbio.2018.05.008
Pinske, C., Jaroschinsky, M., Sargent, F., and Sawers, G. (2012). Zymographic differentiation of [NiFe]-hydrogenases 1, 2 and 3 of Escherichia coli K-12. BMC Microbiol. 12:134. doi: 10.1186/1471-2180-12-134
Pruitt, K. D., Tatusova, T., and Maglott, D. R. (2007). NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 35, D61–D65. doi: 10.1093/nar/gkl842
Schäfer, C., Friedrich, B., and Lenz, O. (2013). Novel, oxygen-insensitive group 5 [NiFe]-hydrogenase in Ralstonia eutropha. Appl. Environ. Microbiol. 79, 5137–5145. doi: 10.1128/AEM.01576-13
Schmidt, C. L., and Shaw, L. (2001). A comprehensive phylogenetic analysis of Rieske and Rieske-type iron-sulfur proteins. J. Bioenerg. Biomembr. 33, 9–26.
Schröder, O., Bleijlevens, B., de Jongh, T. E., Chen, Z., Li, T., Fischer, J., et al. (2007). Characterization of a cyanobacterial-like uptake [NiFe] hydrogenase: EPR and FTIR spectroscopic studies of the enzyme from Acidithiobacillus ferrooxidans. JBIC J. Biol. Inorg. Chem. 12, 212–233. doi: 10.1007/s00775-006-0185-7
Smith, P. K., Krohn, R. II, Hermanson, G. T., Mallia, A. K., Gartner, F. H., and Provenzano, et al. (1985). Measurement of protein using bicinchoninic acid. Anal. Biochem. 150, 76–85. doi: 10.1016/0003-2697(85)90442-7
Snapper, S. B., Melton, R. E., Mustafa, S., Kieser, T., and Jacobs, W. R. J. (1990). Isolation and characterization of efficient plasmid transformation mutants of Mycobacterium smegmatis. Mol. Microbiol. 4, 1911–1919. doi: 10.1111/j.1365-2958.1990.tb02040.x
Søndergaard, D., Pedersen, C. N. S., and Greening, C. (2016). HydDB: a web tool for hydrogenase classification and analysis. Sci. Rep. 6:34212. doi: 10.1038/srep34212
Stover, C. K., De La Cruz, V. F., Fuerst, T. R., Burlein, J. E., Benson, L. A., Bennett, L. T., et al. (1991). New use of BCG for recombinant vaccines. Nature 351:456.
Tamagnini, P., Axelsson, R., Lindberg, P., Oxelfelt, F., Wünschiers, R., and Lindblad, P. (2002). Hydrogenases and hydrogen metabolism of cyanobacteria. Microbiol. Mol. Biol. Rev. 66, 1–20. doi: 10.1128/mmbr.66.1.1-20.2002
Volbeda, A., Darnault, C., Parkin, A., Sargent, F., Armstrong, F. A., and Fontecilla-Camps, J. C. (2013). Crystal structure of the O2-tolerant membrane-bound hydrogenase 1 from Escherichia coli in complex with its cognate cytochrome b. Structure 21, 184–190. doi: 10.1016/j.str.2012.11.010
Walker, J. M. (2009). “Nondenaturing polyacrylamide gel electrophoresis of proteins,” in The Protein Protocols Handbook, ed. J. M. Walker, (Totowa, NJ: Humana Press), 171–176. doi: 10.1007/978-1-59745-198-7_20
Keywords: hydrogenase, Mycobacterium, atmospheric H2, iron-sulfur protein, hydrogen cycle
Citation: Islam ZF, Cordero PRF and Greening C (2019) Putative Iron-Sulfur Proteins Are Required for Hydrogen Consumption and Enhance Survival of Mycobacteria. Front. Microbiol. 10:2749. doi: 10.3389/fmicb.2019.02749
Received: 25 April 2019; Accepted: 12 November 2019;
Published: 22 November 2019.
Edited by:
Patricia Coutinho Dos Santos, Wake Forest University, United StatesReviewed by:
Oliver Lenz, Technische Universität Berlin, GermanyConstanze Pinske, Martin Luther University of Halle-Wittenberg, Germany
Katherine A. Black, Cornell University, United States
Copyright © 2019 Islam, Cordero and Greening. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chris Greening, chris.greening@monash.edu
†These authors have contributed equally to this work