Large-scale genetic characterization of the model sulfate-reducing bacterium, Desulfovibrio vulgaris Hildenborough

Sulfate-reducing bacteria (SRB) are obligate anaerobes that can couple their growth to the reduction of sulfate. Despite the importance of SRB to global nutrient cycles and their damage to the petroleum industry, our molecular understanding of their physiology remains limited. To systematically provide new insights into SRB biology, we generated a randomly barcoded transposon mutant library in the model SRB Desulfovibrio vulgaris Hildenborough (DvH) and used this genome-wide resource to assay the importance of its genes under a range of metabolic and stress conditions. In addition to defining the essential gene set of DvH, we identified a conditional phenotype for 1,137 non-essential genes. Through examination of these conditional phenotypes, we were able to make a number of novel insights into our molecular understanding of DvH, including how this bacterium synthesizes vitamins. For example, we identified DVU0867 as an atypical L-aspartate decarboxylase required for the synthesis of pantothenic acid, provided the first experimental evidence that biotin synthesis in DvH occurs via a specialized acyl carrier protein and without methyl esters, and demonstrated that the uncharacterized dehydrogenase DVU0826:DVU0827 is necessary for the synthesis of pyridoxal phosphate. In addition, we used the mutant fitness data to identify genes involved in the assimilation of diverse nitrogen sources and gained insights into the mechanism of inhibition of chlorate and molybdate. Our large-scale fitness dataset and RB-TnSeq mutant library are community-wide resources that can be used to generate further testable hypotheses into the gene functions of this environmentally and industrially important group of bacteria.

Biofilm Interaction Mapping and Analysis (BIMA) of Interspecific Interactions in Pseudomonas Co-culture Biofilms

Pseudomonas species are ubiquitous in nature and include numerous medically, agriculturally and technologically beneficial strains of which the interspecific interactions are of great interest for biotechnologies. Specifically, co-cultures containing Pseudomonas stutzeri have been used for bioremediation, biocontrol, aquaculture management and wastewater denitrification. Furthermore, the use of P. stutzeri biofilms, in combination with consortia-based approaches, may offer advantages for these processes. Understanding the interspecific interaction within biofilm co-cultures or consortia provides a means for improvement of current technologies. However, the investigation of biofilm-based consortia has been limited. We present an adaptable and scalable method for the analysis of macroscopic interactions (colony morphology, inhibition, and invasion) between colony-forming bacterial strains using an automated printing method followed by analysis of the genes and metabolites involved in the interactions. Using Biofilm Interaction Mapping and Analysis (BIMA), these interactions were investigated between P. stutzeri strain RCH2, a denitrifier isolated from chromium (VI) contaminated soil, and 13 other species of pseudomonas isolated from non-contaminated soil. One interaction partner, Pseudomonas fluorescens N1B4 was selected for mutant fitness profiling of a DNA-barcoded mutant library; with this approach four genes of importance were identified and the effects on interactions were evaluated with deletion mutants and mass spectrometry based metabolomics.

Deletion Mutants, Archived Transposon Library, and Tagged Protein Constructs of the Model Sulfate-Reducing Bacterium Desulfovibrio vulgaris Hildenborough

The dissimilatory sulfate-reducing deltaproteobacterium Desulfovibrio vulgaris Hildenborough (ATCC 29579) was chosen by the research collaboration ENIGMA to explore tools and protocols for bringing this anaerobe to model status. Here, we describe a collection of genetic constructs generated by ENIGMA that are available to the research community.

Novel Mode of Molybdate Inhibition of Desulfovibrio vulgaris Hildenborough

Sulfate-reducing microorganisms (SRM) are found in multiple environments and play a major role in global carbon and sulfur cycling. Because of their growth capabilities and association with metal corrosion, controlling the growth of SRM has become of increased interest. One such mechanism of control has been the use of molybdate (MoO4 2-), which is thought to be a specific inhibitor of SRM. The way in which molybdate inhibits the growth of SRM has been enigmatic. It has been reported that molybdate is involved in a futile energy cycle with the sulfate-activating enzyme, sulfate adenylyl transferase (Sat), which results in loss of cellular ATP. However, we show here that a deletion of this enzyme in the model SRM, Desulfovibrio vulgaris Hildenborough, remained sensitive to molybdate. We performed several subcultures of the ∆sat strain in the presence of increasing concentrations of molybdate and obtained a culture with increased resistance to the inhibitor (up to 3 mM). The culture was re-sequenced and three single nucleotide polymorphisms (SNPs) were identified that were not present in the parental strain. Two of the SNPs seemed unlikely candidates for molybdate resistance due to a lack of conservation of the mutated residues in homologous genes of closely related strains. The remaining SNP was located in DVU2210, a protein containing two domains: a YcaO-like domain and a tetratricopeptide-repeat domain. The SNP resulted in a change of a serine residue to arginine in the ATP-hydrolyzing motif of the YcaO-like domain. Deletion mutants of each of the three genes apparently enriched with SNPs in the presence of inhibitory molybdate and combinations of these genes were generated in the Δsat and wild-type strains. Strains lacking both sat and DVU2210 became more resistant to molybdate. Deletions of the other two genes in which SNPs were observed did not result in increased resistance to molybdate. YcaO-like proteins are distributed across the bacterial and archaeal domains, though the function of these proteins is largely unknown. The role of this protein in D. vulgaris is unknown. Due to the distribution of YcaO-like proteins in prokaryotes, the veracity of molybdate as a specific SRM inhibitor should be reconsidered.

The Structure of an AAV5-AAVR Complex at 2.5 Å Resolution: Implications for Cellular Entry and Immune Neutralization of AAV Gene Therapy Vectors

Adeno-Associated Virus is the leading vector for gene therapy. Although it is the vector for all in vivo gene therapies approved for clinical use by the US Food and Drug Administration, its biology is still not yet fully understood. It has been shown that different serotypes of AAV bind to their cellular receptor, AAVR, in different ways. Previously we have reported a 2.4Å structure of AAV2 bound to AAVR that shows ordered structure for only one of the two AAVR domains with which AAV2 interacts. In this study we present a 2.5Å resolution structure of AAV5 bound to AAVR. AAV5 binds to the first polycystic kidney disease (PKD) domain of AAVR that was not ordered in the AAV2 structure. Interactions of AAV5 with AAVR are analyzed in detail, and the implications for AAV2 binding are explored through molecular modeling. Moreover, we find that binding sites for the antibodies ADK5a, ADK5b, and 3C5 on AAV5 overlap with the binding site of AAVR. These insights provide a structural foundation for development of gene therapy agents to better evade immune neutralization without disrupting cellular entry.

Characterization of subsurface media from locations up- and down-gradient of a uranium-contaminated aquifer

The processing of sediment to accurately characterize the spatially-resolved depth profiles of geophysical and geochemical properties along with signatures of microbial density and activity remains a challenge especially in complex contaminated areas. This study processed cores from two sediment boreholes from background and contaminated core sediments and surrounding groundwater. Fresh core sediments were compared by depth to capture the changes in sediment structure, sediment minerals, biomass, and pore water geochemistry in terms of major and trace elements including pollutants, cations, anions, and organic acids. Soil porewater samples were matched to groundwater level, flow rate, and preferential flows and compared to homogenized groundwater-only samples from neighboring monitoring wells. Groundwater analysis of nearby wells only revealed high sulfate and nitrate concentrations while the same analysis using sediment pore water samples with depth was able to suggest areas high in sulfate- and nitrate-reducing bacteria based on their decreased concentration and production of reduced by-products that could not be seen in the groundwater samples. Positive correlations among porewater content, total organic carbon, trace metals and clay minerals revealed a more complicated relationship among contaminant, sediment texture, groundwater table, and biomass. The fluctuating capillary interface had high concentrations of Fe and Mn-oxides combined with trace elements including U, Th, Sr, Ba, Cu, and Co. This suggests the mobility of potentially hazardous elements, sediment structure, and biogeochemical factors are all linked together to impact microbial communities, emphasizing that solid interfaces play an important role in determining the abundance of bacteria in the sediments.

Expression and Purification of Adeno-associated Virus Virus-like Particles in a Baculovirus System and AAVR Ectodomain Constructs in E. coli

Adeno-associated virus (AAV) is a promising gene therapy vector and the biophysical characterization of its interactions with host proteins is a critical foundation for engineering tissue targeting and immune escape. Presented here are protocols for the production of: (a) the outer protein shells (virus-like particles or VLPs) for serotype 2 (AAV-2) and (b) two fragments from the binding ectodomain of AAV's cellular receptor, AAVR. His6PKD1-2 comprises the first two polycystic kidney disease (PKD) domains, the minimal required for efficient binding of AAV, expressed with an N-terminal histidine tag. MBP-PKD1-5 is a fusion of the maltose binding protein with all five of the PKD domains of the AAVR receptor. Presented are the expression and purification of milligram quantities, ample for in vitro analyses. For AAV-2, the protocol offers an alternative to the use of (infectious) wild-type virus or transducing vectors. One of the methods for producing transducing vector is in Sf9 cells, and the production of VLPs is based on this. For AAVR, the protocols enable biochemical and biophysical characterization of virus-binding. The minimal two-domain construct allows more saturated binding to symmetry-equivalent sites on the virus, while the larger construct might be better expected to reflect the native receptor.

LurR is a regulator of the central lactate oxidation pathway in sulfate-reducing Desulfovibrio species

The central carbon/lactate utilization pathway in the model sulfate-reducing bacterium, Desulfovibrio vulgaris Hildenborough, is encoded by the highly conserved operon DVU3025-3033. Our earlier in vitro genome-wide study had suggested a network of four two-component system regulators that target this large operon; however, how these four regulators control this operon was not known. Here, we probe the regulation of the lactate utilization operon with mutant strains and DNA-protein binding assays. We show that the LurR response regulator is required for optimal growth and complete lactate utilization, and that it activates the DVU3025-3033 lactate oxidation operon as well as DVU2451, a lactate permease gene, in the presence of lactate. We show by electrophoretic mobility shift assays that LurR binds to three sites in the upstream region of DVU3025, the first gene of the operon. NrfR, a response regulator that is activated under nitrite stress, and LurR share similar binding site motifs and bind the same sites upstream of DVU3025. The DVU3025 promoter also has a binding site motif (Pho box) that is bound by PhoB, a two-component response regulator activated under phosphate limitation. The lactate utilization operon, the regulator LurR, and LurR binding sites are conserved across the order Desulfovibrionales whereas possible modulation of the lactate utilization genes by additional regulators such as NrfR and PhoB appears to be limited to D. vulgaris.

A new family of transcriptional regulators of tungstoenzymes and molybdate/tungstate transport

Bacterial genes for molybdenum-containing and tungsten-containing enzymes are often differentially regulated depending on the metal availability in the environment. Here, we describe a new family of transcription factors with an unusual DNA-binding domain related to excisionases of bacteriophages. These transcription factors are associated with genes for various molybdate and tungstate-specific transporting systems as well as molybdo/tungsto-enzymes in a wide range of bacterial genomes. We used a combination of computational and experimental techniques to study a member of the TF family, named TaoR (for tungsten-containing aldehyde oxidoreductase regulator). In Desulfovibrio vulgaris Hildenborough, a model bacterium for sulfate reduction studies, TaoR activates expression of aldehyde oxidoreductase aor and represses tungsten-specific ABC-type transporter tupABC genes under tungsten-replete conditions. TaoR binding sites at aor promoter were identified by electrophoretic mobility shift assay and DNase I footprinting. We also reconstructed TaoR regulons in 45 Deltaproteobacteria by comparative genomics approach and predicted target genes for TaoR family members in other Proteobacteria and Firmicutes.

Deconstructing the Dissimilatory Sulfate Reduction Pathway: Isotope Fractionation of a Mutant Unable of Growth on Sulfate

The sulfur isotope record provides key insight into the history of Earth's redox conditions. A detailed understanding of the metabolisms driving this cycle, and specifically microbial sulfate reduction (MSR), is crucial for accurate paleoenvironmental reconstructions. This includes a precise knowledge of the step-specific sulfur isotope effects during MSR. In this study, we aim at resolving the cellular-level fractionation factor during dissimilatory sulfite reduction to sulfide within MSR, and use this measured isotope effect as a calibration to enhance our understanding of the biochemistry of sulfite reduction. For this, we merge measured isotope effects associated with dissimilatory sulfite reduction with a quantitative model that explicitly links net fractionation, reaction reversibility, and intracellular metabolite levels. The highly targeted experimental aspect of this study was possible by virtue of the availability of a deletion mutant strain of the model sulfate reducer Desulfovibrio vulgaris (strain Hildenborough), in which the sulfite reduction step is isolated from the rest of the metabolic pathway owing to the absence of its QmoABC complex (ΔQmo). This deletion disrupts electron flux and prevents the reduction of adenosine phosphosulfate (APS) to sulfite. When grown in open-system steady-state conditions at 10% maximum growth rate in the presence of sulfite and lactate as electron donor, sulfur isotope fractionation factors averaged -15.9‰ (1 σ = 0.4), which appeared to be statistically indistinguishable from a pure enzyme study with dissimilatory sulfite reductase. We coupled these measurements with an understanding of step-specific equilibrium and kinetic isotope effects, and furthered our mechanistic understanding of the biochemistry of sulfite uptake and ensuing reduction. Our metabolically informed isotope model identifies flavodoxin as the most likely electron carrier performing the transfer of electrons to dissimilatory sulfite reductase. This is in line with previous work on metabolic strategies adopted by sulfate reducers under different energy regimes, and has implications for our understanding of the plasticity of this metabolic pathway at the center of our interpretation of modern and palaeo-environmental records.

Iron- and aluminium-induced depletion of molybdenum in acidic environments impedes the nitrogen cycle

Anthropogenic nitrate contamination is a serious problem in many natural environments. Nitrate removal by microbial action is dependent on the metal molybdenum (Mo), which is required by nitrate reductase for denitrification and dissimilatory nitrate reduction to ammonium. The soluble form of Mo, molybdate (MoO₄ ^2- ), is incorporated into and adsorbed by iron (Fe) and aluminium (Al) (oxy) hydroxide minerals. Herein we used Oak Ridge Reservation (ORR) as a model nitrate-contaminated acidic environment to investigate whether the formation of Fe- and Al-precipitates could impede microbial nitrate removal by depleting Mo. We demonstrate that Fe and Al mineral formation that occurs as the pH of acidic synthetic groundwater is increased, decreases soluble Mo to low picomolar concentrations, a process proposed to mimic environmental diffusion of acidic contaminated groundwater. Analysis of ORR sediments revealed recalcitrant Mo in the contaminated core that co-occurred with Fe and Al, consistent with Mo scavenging by Fe/Al precipitates. Nitrate removal by ORR isolate Pseudomonas fluorescens N2A2 is virtually abolished by Fe/Al precipitate-induced Mo depletion. The depletion of naturally occurring Mo in nitrate- and Fe/Al-contaminated acidic environments like ORR or acid mine drainage sites has the potential to impede microbial-based nitrate reduction thereby extending the duration of nitrate in the environment.

Mutant phenotypes for thousands of bacterial genes of unknown function

One-third of all protein-coding genes from bacterial genomes cannot be annotated with a function. Here, to investigate the functions of these genes, we present genome-wide mutant fitness data from 32 diverse bacteria across dozens of growth conditions. We identified mutant phenotypes for 11,779 protein-coding genes that had not been annotated with a specific function. Many genes could be associated with a specific condition because the gene affected fitness only in that condition, or with another gene in the same bacterium because they had similar mutant phenotypes. Of the poorly annotated genes, 2,316 had associations that have high confidence because they are conserved in other bacteria. By combining these conserved associations with comparative genomics, we identified putative DNA repair proteins; in addition, we propose specific functions for poorly annotated enzymes and transporters and for uncharacterized protein families. Our study demonstrates the scalability of microbial genetics and its utility for improving gene annotations.

Cr(VI) reduction and physiological toxicity are impacted by resource ratio in Desulfovibrio vulgaris

Desulfovibrio spp. are capable of heavy metal reduction and are well-studied systems for understanding metal fate and transport in anaerobic environments. Desulfovibrio vulgaris Hildenborough was grown under environmentally relevant conditions (i.e., temperature, nutrient limitation) to elucidate the impacts on Cr(VI) reduction on cellular physiology. Growth at 20 °C was slower than 30 °C and the presence of 50 μM Cr(VI) caused extended lag times for all conditions, but once growth resumed the growth rate was similar to that without Cr(VI). Cr(VI) reduction rates were greatly diminished at 20 °C for both 50 and 100 μM Cr(VI), particularly for the electron acceptor limited (EAL) condition in which Cr(VI) reduction was much slower, the growth lag much longer (200 h), and viability decreased compared to balanced (BAL) and electron donor limited (EDL) conditions. When sulfate levels were increased in the presence of Cr(VI), cellular responses improved via a shorter lag time to growth. Similar results were observed between the different resource (donor/acceptor) ratio conditions when the sulfate levels were normalized (10 mM), and these results indicated that resource ratio (donor/acceptor) impacted D. vulgaris response to Cr(VI) and not merely sulfate limitation. The results suggest that temperature and resource ratios greatly impacted the extent of Cr(VI) toxicity, Cr(VI) reduction, and the subsequent cellular health via Cr(VI) influx and overall metabolic rate. The results also emphasized the need to perform experiments at lower temperatures with nutrient limitation to make accurate predictions of heavy metal reduction rates as well as physiological states in the environment.

Filling gaps in bacterial amino acid biosynthesis pathways with high-throughput genetics

For many bacteria with sequenced genomes, we do not understand how they synthesize some amino acids. This makes it challenging to reconstruct their metabolism, and has led to speculation that bacteria might be cross-feeding amino acids. We studied heterotrophic bacteria from 10 different genera that grow without added amino acids even though an automated tool predicts that the bacteria have gaps in their amino acid synthesis pathways. Across these bacteria, there were 11 gaps in their amino acid biosynthesis pathways that we could not fill using current knowledge. Using genome-wide mutant fitness data, we identified novel enzymes that fill 9 of the 11 gaps and hence explain the biosynthesis of methionine, threonine, serine, or histidine by bacteria from six genera. We also found that the sulfate-reducing bacterium Desulfovibrio vulgaris synthesizes homocysteine (which is a precursor to methionine) by using DUF39, NIL/ferredoxin, and COG2122 proteins, and that homoserine is not an intermediate in this pathway. Our results suggest that most free-living bacteria can likely make all 20 amino acids and illustrate how high-throughput genetics can uncover previously-unknown amino acid biosynthesis genes.

For a few bacteria, it is well known how they can make all 20 of the standard amino acids (the building blocks of proteins). For many other bacteria, their genome sequence implies that there are gaps in these biosynthetic pathways, so that the bacteria cannot make all of the amino acids and would need to take up some of them from their environment instead. But many bacteria can grow in minimal media (without any amino acids) despite these apparent gaps. We studied 10 bacteria with predicted gaps in amino acid biosynthesis that nevertheless grow in minimal media. Most of these gaps were spurious, but 11 of the gaps were genuine and could not be explained by current knowledge. Using high-throughput genetics, we systematically identified genes that were required for growth in minimal media and identified the biosynthetic genes that fill 9 of the 11 gaps. We hope that this approach can be applied to many more bacteria and will eventually allow us to accurately predict the nutritional requirements of a bacterium from its genome sequence.

Unintended Laboratory-Driven Evolution Reveals Genetic Requirements for Biofilm Formation by Desulfovibrio vulgaris Hildenborough

Biofilms of sulfate-reducing bacteria (SRB) are of particular interest as members of this group are culprits in corrosion of industrial metal and concrete pipelines as well as being key players in subsurface metal cycling. Yet the mechanism of biofilm formation by these bacteria has not been determined. Here we show that two supposedly identical wild-type cultures of the SRB Desulfovibrio vulgaris Hildenborough maintained in different laboratories have diverged in biofilm formation. From genome resequencing and subsequent mutant analyses, we discovered that a single nucleotide change within DVU1017, the ABC transporter of a type I secretion system (T1SS), was sufficient to eliminate biofilm formation in D. vulgaris Hildenborough. Two T1SS cargo proteins were identified as likely biofilm structural proteins, and the presence of at least one (with either being sufficient) was shown to be required for biofilm formation. Antibodies specific to these biofilm structural proteins confirmed that DVU1017, and thus the T1SS, is essential for localization of these adhesion proteins on the cell surface. We propose that DVU1017 is a member of the lapB category of microbial surface proteins because of its phenotypic similarity to the adhesin export system described for biofilm formation in the environmental pseudomonads. These findings have led to the identification of two functions required for biofilm formation in D. vulgaris Hildenborough and focus attention on the importance of monitoring laboratory-driven evolution, as phenotypes as fundamental as biofilm formation can be altered.IMPORTANCE The growth of bacteria attached to a surface (i.e., biofilm), specifically biofilms of sulfate-reducing bacteria, has a profound impact on the economy of developed nations due to steel and concrete corrosion in industrial pipelines and processing facilities. Furthermore, the presence of sulfate-reducing bacteria in oil wells causes oil souring from sulfide production, resulting in product loss, a health hazard to workers, and ultimately abandonment of wells. Identification of the required genes is a critical step for determining the mechanism of biofilm formation by sulfate reducers. Here, the transporter by which putative biofilm structural proteins are exported from sulfate-reducing Desulfovibrio vulgaris Hildenborough cells was discovered, and a single nucleotide change within the gene coding for this transporter was found to be sufficient to completely stop formation of biofilm.

Mechanisms of Chromium and Uranium Toxicity in Pseudomonas stutzeri RCH2 Grown under Anaerobic Nitrate-Reducing Conditions

Chromium and uranium are highly toxic metals that contaminate many natural environments. We investigated their mechanisms of toxicity under anaerobic conditions using nitrate-reducing Pseudomonas stutzeri RCH2, which was originally isolated from a chromium-contaminated aquifer. A random barcode transposon site sequencing library of RCH2 was grown in the presence of the chromate oxyanion (Cr[VI]O42−) or uranyl oxycation (U[VI]O22+). Strains lacking genes required for a functional nitrate reductase had decreased fitness as both metals interacted with heme-containing enzymes required for the later steps in the denitrification pathway after nitrate is reduced to nitrite. Cr[VI]-resistance also required genes in the homologous recombination and nucleotide excision DNA repair pathways, showing that DNA is a target of Cr[VI] even under anaerobic conditions. The reduced thiol pool was also identified as a target of Cr[VI] toxicity and psest_2088, a gene of previously unknown function, was shown to have a role in the reduction of sulfite to sulfide. U[VI] resistance mechanisms involved exopolysaccharide synthesis and the universal stress protein UspA. As the first genome-wide fitness analysis of Cr[VI] and U[VI] toxicity under anaerobic conditions, this study provides new insight into the impact of Cr[VI] and U[VI] on an environmental isolate from a chromium contaminated site, as well as into the role of a ubiquitous protein, Psest_2088.

Mechanism for microbial population collapse in a fluctuating resource environment

Managing trade-offs through gene regulation is believed to confer resilience to a microbial community in a fluctuating resource environment. To investigate this hypothesis, we imposed a fluctuating environment that required the sulfate-reducer Desulfovibrio vulgaris to undergo repeated ecologically relevant shifts between retaining metabolic independence (active capacity for sulfate respiration) and becoming metabolically specialized to a mutualistic association with the hydrogen-consuming Methanococcus maripaludis Strikingly, the microbial community became progressively less proficient at restoring the environmentally relevant physiological state after each perturbation and most cultures collapsed within 3-7 shifts. Counterintuitively, the collapse phenomenon was prevented by a single regulatory mutation. We have characterized the mechanism for collapse by conducting RNA-seq analysis, proteomics, microcalorimetry, and single-cell transcriptome analysis. We demonstrate that the collapse was caused by conditional gene regulation, which drove precipitous decline in intracellular abundance of essential transcripts and proteins, imposing greater energetic burden of regulation to restore function in a fluctuating environment.

Novel Metal Cation Resistance Systems from Mutant Fitness Analysis of Denitrifying Pseudomonas stutzeri

Metal ion transport systems have been studied extensively, but the specificity of a given transporter is often unclear from amino acid sequence data alone. In this study, predicted Cu(2+) and Zn(2+) resistance systems in Pseudomonas stutzeri strain RCH2 are compared with those experimentally implicated in Cu(2+) and Zn(2+) resistance, as determined by using a DNA-barcoded transposon mutant library. Mutant fitness data obtained under denitrifying conditions are combined with regulon predictions to yield a much more comprehensive picture of Cu(2+) and Zn(2+) resistance in strain RCH2. The results not only considerably expand what is known about well-established metal ion exporters (CzcCBA, CzcD, and CusCBA) and their accessory proteins (CzcI and CusF), they also reveal that isolates with mutations in some predicted Cu(2+) resistance systems do not show decreased fitness relative to the wild type when exposed to Cu(2+) In addition, new genes are identified that have no known connection to Zn(2+) (corB, corC, Psest_3226, Psest_3322, and Psest_0618) or Cu(2+) resistance (Mrp antiporter subunit gene, Psest_2850, and Psest_0584) but are crucial for resistance to these metal cations. Growth of individual deletion mutants lacking corB, corC, Psest_3226, or Psest_3322 confirmed the observed Zn-dependent phenotypes. Notably, to our knowledge, this is the first time a bacterial homolog of TMEM165, a human gene responsible for a congenital glycosylation disorder, has been deleted and the resulting strain characterized. Finally, the fitness values indicate Cu(2+)- and Zn(2+)-based inhibition of nitrite reductase and interference with molybdenum cofactor biosynthesis for nitrate reductase. These results extend the current understanding of Cu(2+) and Zn(2+) efflux and resistance and their effects on denitrifying metabolism.

IMPORTANCE

In this study, genome-wide mutant fitness data in P. stutzeri RCH2 combined with regulon predictions identify several proteins of unknown function that are involved in resisting zinc and copper toxicity. For zinc, these include a member of the UPF0016 protein family that was previously implicated in Ca(2+)/H(+) antiport and a human congenital glycosylation disorder, CorB and CorC, which were previously linked to Mg(2+) transport, and Psest_3322 and Psest_0618, two proteins with no characterized homologs. Experiments using mutants lacking Psest_3226, Psest_3322, corB, corC, or czcI verified their proposed functions, which will enable future studies of these little-characterized zinc resistance determinants. Likewise, Psest_2850, annotated as an ion antiporter subunit, and the conserved hypothetical protein Psest_0584 are implicated in copper resistance. Physiological connections between previous studies and phenotypes presented here are discussed. Functional and mechanistic understanding of transport proteins improves the understanding of systems in which members of the same protein family, including those in humans, can have different functions.

Rapid selective sweep of pre-existing polymorphisms and slow fixation of new mutations in experimental evolution of Desulfovibrio vulgaris

To investigate the genetic basis of microbial evolutionary adaptation to salt (NaCl) stress, populations of Desulfovibrio vulgaris Hildenborough (DvH), a sulfate-reducing bacterium important for the biogeochemical cycling of sulfur, carbon and nitrogen, and potentially the bioremediation of toxic heavy metals and radionuclides, were propagated under salt stress or non-stress conditions for 1200 generations. Whole-genome sequencing revealed 11 mutations in salt stress-evolved clone ES9-11 and 14 mutations in non-stress-evolved clone EC3-10. Whole-population sequencing data suggested the rapid selective sweep of the pre-existing polymorphisms under salt stress within the first 100 generations and the slow fixation of new mutations. Population genotyping data demonstrated that the rapid selective sweep of pre-existing polymorphisms was common in salt stress-evolved populations. In contrast, the selection of pre-existing polymorphisms was largely random in EC populations. Consistently, at 100 generations, stress-evolved population ES9 showed improved salt tolerance, namely increased growth rate (2.0-fold), higher biomass yield (1.8-fold) and shorter lag phase (0.7-fold) under higher salinity conditions. The beneficial nature of several mutations was confirmed by site-directed mutagenesis. All four tested mutations contributed to the shortened lag phases under higher salinity condition. In particular, compared with the salt tolerance improvement in ES9-11, a mutation in a histidine kinase protein gene lytS contributed 27% of the growth rate increase and 23% of the biomass yield increase while a mutation in hypothetical gene DVU2472 contributed 24% of the biomass yield increase. Our results suggested that a few beneficial mutations could lead to dramatic improvements in salt tolerance.

Biofilm growth mode promotes maximum carrying capacity and community stability during product inhibition syntrophy

Sulfate-reducing bacteria (SRB) can interact syntrophically with other community members in the absence of sulfate, and interactions with hydrogen-consuming methanogens are beneficial when these archaea consume potentially inhibitory H₂ produced by the SRB. A dual continuous culture approach was used to characterize population structure within a syntrophic biofilm formed by the SRB Desulfovibrio vulgaris Hildenborough and the methanogenic archaeum Methanococcus maripaludis. Under the tested conditions, monocultures of D. vulgaris formed thin, stable biofilms, but monoculture M. maripaludis did not. Microscopy of intact syntrophic biofilm confirmed that D. vulgaris formed a scaffold for the biofilm, while intermediate and steady-state images revealed that M. maripaludis joined the biofilm later, likely in response to H₂ produced by the SRB. Close interactions in structured biofilm allowed efficient transfer of H₂ to M. maripaludis, and H₂ was only detected in cocultures with a mutant SRB that was deficient in biofilm formation (ΔpilA). M. maripaludis produced more carbohydrate (uronic acid, hexose, and pentose) as a monoculture compared to total coculture biofilm, and this suggested an altered carbon flux during syntrophy. The syntrophic biofilm was structured into ridges (∼300 × 50 μm) and models predicted lactate limitation at ∼50 μm biofilm depth. The biofilm had structure that likely facilitated mass transfer of H₂ and lactate, yet maximized biomass with a more even population composition (number of each organism) when compared to the bulk-phase community. Total biomass protein was equivalent in lactate-limited and lactate-excess conditions when a biofilm was present, but in the absence of biofilm, total biomass protein was significantly reduced. The results suggest that multispecies biofilms create an environment conducive to resource sharing, resulting in increased biomass retention, or carrying capacity, for cooperative populations.

Erosion of functional independence early in the evolution of a microbial mutualism

Many species have evolved to function as specialized mutualists, often to the detriment of their ability to survive independently. However, there are few, if any, well-controlled observations of the evolutionary processes underlying the genesis of new mutualisms. Here, we show that within the first 1,000 generations of initiating independent syntrophic interactions between a sulfate reducer (Desulfovibrio vulgaris) and a hydrogenotrophic methanogen (Methanococcus maripaludis), D. vulgaris frequently lost the capacity to grow by sulfate respiration, thus losing the primary physiological attribute of the genus. The loss of sulfate respiration was a consequence of mutations in one or more of three key genes in the pathway for sulfate respiration, required for sulfate activation (sat) and sulfate reduction to sulfite (apsA or apsB). Because loss-of-function mutations arose rapidly and independently in replicated experiments, and because these mutations were correlated with enhanced growth rate and productivity, gene loss could be attributed to natural selection, even though these mutations should significantly restrict the independence of the evolved D. vulgaris. Together, these data present an empirical demonstration that specialization for a mutualistic interaction can evolve by natural selection shortly after its origin. They also demonstrate that a sulfate-reducing bacterium can readily evolve to become a specialized syntroph, a situation that may have often occurred in nature.

Identification of a cyclic-di-GMP-modulating response regulator that impacts biofilm formation in a model sulfate reducing bacterium

We surveyed the eight putative cyclic-di-GMP-modulating response regulators (RRs) in Desulfovibrio vulgaris Hildenborough that are predicted to function via two-component signaling. Using purified proteins, we examined cyclic-di-GMP (c-di-GMP) production or turnover in vitro of all eight proteins. The two RRs containing only GGDEF domains (DVU2067, DVU0636) demonstrated c-di-GMP production activity in vitro. Of the remaining proteins, three RRs with HD-GYP domains (DVU0722, DVUA0086, and DVU2933) were confirmed to be Mn²⁺-dependent phosphodiesterases (PDEs) in vitro and converted c-di-GMP to its linear form, pGpG. DVU0408, containing both c-di-GMP production (GGDEF) and degradation domains (EAL), showed c-di-GMP turnover activity in vitro also with production of pGpG. No c-di-GMP related activity could be assigned to the RR DVU0330, containing a metal-dependent phosphohydrolase HD-OD domain, or to the HD-GYP domain RR, DVU1181. Studies included examining the impact of overexpressed cyclic-di-GMP-modulating RRs in the heterologous host E. coli and led to the identification of one RR, DVU0636, with increased cellulose production. Evaluation of a transposon mutant in DVU0636 indicated that the strain was impaired in biofilm formation and demonstrated an altered carbohydrate:protein ratio relative to the D. vulgaris wild type biofilms. However, grown in liquid lactate/sulfate medium, the DVU0636 transposon mutant showed no growth impairment relative to the wild-type strain. Among the eight candidates, only the transposon disruption mutant in the DVU2067 RR presented a growth defect in liquid culture. Our results indicate that, of the two diguanylate cyclases (DGCs) that function as part of two-component signaling, DVU0636 plays an important role in biofilm formation while the function of DVU2067 has pertinence in planktonic growth.

Genetic basis for nitrate resistance in Desulfovibrio strains

Nitrate is an inhibitor of sulfate-reducing bacteria (SRB). In petroleum production sites, amendments of nitrate and nitrite are used to prevent SRB production of sulfide that causes souring of oil wells. A better understanding of nitrate stress responses in the model SRB, Desulfovibrio vulgaris Hildenborough and Desulfovibrio alaskensis G20, will strengthen predictions of environmental outcomes of nitrate application. Nitrate inhibition of SRB has historically been considered to result from the generation of small amounts of nitrite, to which SRB are quite sensitive. Here we explored the possibility that nitrate might inhibit SRB by a mechanism other than through nitrite inhibition. We found that nitrate-stressed D. vulgaris cultures grown in lactate-sulfate conditions eventually grew in the presence of high concentrations of nitrate, and their resistance continued through several subcultures. Nitrate consumption was not detected over the course of the experiment, suggesting adaptation to nitrate. With high-throughput genetic approaches employing TnLE-seq for D. vulgaris and a pooled mutant library of D. alaskensis, we determined the fitness of many transposon mutants of both organisms in nitrate stress conditions. We found that several mutants, including homologs present in both strains, had a greatly increased ability to grow in the presence of nitrate but not nitrite. The mutated genes conferring nitrate resistance included the gene encoding the putative Rex transcriptional regulator (DVU0916/Dde_2702), as well as a cluster of genes (DVU0251-DVU0245/Dde_0597-Dde_0605) that is poorly annotated. Follow-up studies with individual D. vulgaris transposon and deletion mutants confirmed high-throughput results. We conclude that, in D. vulgaris and D. alaskensis, nitrate resistance in wild-type cultures is likely conferred by spontaneous mutations. Furthermore, the mechanisms that confer nitrate resistance may be different from those that confer nitrite resistance.

Exploring the role of CheA3 in Desulfovibrio vulgaris Hildenborough motility

Sulfate-reducing bacteria such as Desulfovibrio vulgaris Hildenborough are often found in environments with limiting growth nutrients. Using lactate as the electron donor and carbon source, and sulfate as the electron acceptor, wild type D. vulgaris shows motility on soft agar plates. We evaluated this phenotype with mutants resulting from insertional inactivation of genes potentially related to motility. Our study revealed that the cheA3 (DVU2072) kinase mutant was impaired in the ability to form motility halos. Insertions in two other cheA loci did not exhibit a loss in this phenotype. The cheA3 mutant was also non-motile in capillary assays. Complementation with a plasmid-borne copy of cheA3 restores wild type phenotypes. The cheA3 mutant displayed a flagellum as observed by electron microscopy, grew normally in liquid medium, and was motile in wet mounts. In the growth conditions used, the D. vulgaris ΔfliA mutant (DVU3229) for FliA, predicted to regulate flagella-related genes including cheA3, was defective both in flagellum formation and in forming the motility halos. In contrast, a deletion of the flp gene (DVU2116) encoding a pilin-related protein was similar to wild type. We conclude that wild type D. vulgaris forms motility halos on solid media that are mediated by flagella-related mechanisms via the CheA3 kinase. The conditions under which the CheA1 (DVU1594) and CheA2 (DVU1960) kinase function remain to be explored.

Fractionation of sulfur isotopes by Desulfovibrio vulgaris mutants lacking hydrogenases or type I tetraheme cytochrome c 3

The sulfur isotope effect produced by sulfate reducing microbes is commonly used to trace biogeochemical cycles of sulfur and carbon in aquatic and sedimentary environments. To test the contribution of intracellular coupling between carbon and sulfur metabolisms to the overall magnitude of the sulfur isotope effect, this study compared sulfur isotope fractionations by mutants of Desulfovibrio vulgaris Hildenborough. We tested mutant strains lacking one or two periplasmic (Hyd, Hyn-1, Hyn-2, and Hys) or cytoplasmic hydrogenases (Ech and CooL), and a mutant lacking type I tetraheme cytochrome (TpI-c 3). In batch culture, wild-type D. vulgaris and its hydrogenase mutants had comparable growth kinetics and produced the same sulfur isotope effects. This is consistent with the reported redundancy of hydrogenases in D. vulgaris. However, the TpI-c 3 mutant (ΔcycA) exhibited slower growth and sulfate reduction rates in batch culture, and produced more H2 and an approximately 50% larger sulfur isotope effect, compared to the wild type. The magnitude of sulfur isotope fractionation in the CycA deletion strain, thus, increased due to the disrupted coupling of the carbon oxidation and sulfate reduction pathways. In continuous culture, wild-type D. vulgaris and the CycA mutant produced similar sulfur isotope effects, underscoring the influence of environmental conditions on the relative contribution of hydrogen cycling to the electron transport. The large sulfur isotope effects associated with the non-ideal stoichiometry of sulfate reduction in this study imply that simultaneous fermentation and sulfate reduction may be responsible for some of the large naturally-occurring sulfur isotope effects. Overall, mutant strains provide a powerful tool to test the effect of specific redox proteins and pathways on sulfur isotope fractionation.

Effect of the deletion of qmoABC and the promoter-distal gene encoding a hypothetical protein on sulfate reduction in Desulfovibrio vulgaris Hildenborough

The pathway of electrons required for the reduction of sulfate in sulfate-reducing bacteria (SRB) is not yet fully characterized. In order to determine the role of a transmembrane protein complex suggested to be involved in this process, a deletion in Desulfovibrio vulgaris Hildenborough was created by marker exchange mutagenesis that eliminated four genes putatively encoding the QmoABC complex and a hypothetical protein (DVU0851). The Qmo (quinone-interacting membrane-bound oxidoreductase) complex is proposed to be responsible for transporting electrons to the dissimilatory adenosine-5'-phosphosulfate reductase in SRB. In support of the predicted role of this complex, the deletion mutant was unable to grow using sulfate as its sole electron acceptor with a range of electron donors. To explore a possible role for the hypothetical protein in sulfate reduction, a second mutant was constructed that had lost only the gene that codes for the DVU0851 protein. The second constructed mutant grew with sulfate as the sole electron acceptor; however, there was a lag that was not present with the wild-type or complemented strain. Neither deletion strain was significantly impaired for growth with sulfite or thiosulfate as the terminal electron acceptor. Complementation of the Delta(qmoABC-DVU0851) mutant with all four genes or only the qmoABC genes restored its ability to grow by sulfate respiration. These results confirmed the prediction that the Qmo complex is in the electron pathway for sulfate reduction and revealed that no other transmembrane complex could compensate when Qmo was lacking.