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.

Rex (encoded by DVU_0916) in Desulfovibrio vulgaris Hildenborough is a repressor of sulfate adenylyl transferase and is regulated by NADH.

Although the enzymes for dissimilatory sulfate reduction by microbes have been studied, the mechanisms for transcriptional regulation of the encoding genes remain unknown. In a number of bacteria the transcriptional regulator Rex has been shown to play a key role as a repressor of genes producing proteins involved in energy conversion. In the model sulfate-reducing microbe Desulfovibrio vulgaris Hildenborough, the gene DVU_0916 was observed to resemble other known Rex proteins. Therefore, the DVU_0916 protein has been predicted to be a transcriptional repressor of genes encoding proteins that function in the process of sulfate reduction in D. vulgaris Hildenborough. Examination of the deduced DVU_0916 protein identified two domains, one a winged helix DNA-binding domain common for transcription factors, and the other a Rossman fold that could potentially interact with pyridine nucleotides. A deletion of the putative rex gene was made in D. vulgaris Hildenborough, and transcript expression studies of sat, encoding sulfate adenylyl transferase, showed increased levels in the D. vulgaris Hildenborough Rex (RexDvH) mutant relative to the parental strain. The RexDvH-binding site upstream of sat was identified, confirming RexDvH to be a repressor of sat. We established in vitro that the presence of elevated NADH disrupted the interaction between RexDvH and DNA. Examination of the 5' transcriptional start site for the sat mRNA revealed two unique start sites, one for respiring cells that correlated with the RexDvH-binding site and a second for fermenting cells. Collectively, these data support the role of RexDvH as a transcription repressor for sat that senses the redox status of the cell.

Generalized schemes for high-throughput manipulation of the Desulfovibrio vulgaris genome.

The ability to conduct advanced functional genomic studies of the thousands of sequenced bacteria has been hampered by the lack of available tools for making high-throughput chromosomal manipulations in a systematic manner that can be applied across diverse species. In this work, we highlight the use of synthetic biological tools to assemble custom suicide vectors with reusable and interchangeable DNA "parts" to facilitate chromosomal modification at designated loci. These constructs enable an array of downstream applications, including gene replacement and the creation of gene fusions with affinity purification or localization tags. We employed this approach to engineer chromosomal modifications in a bacterium that has previously proven difficult to manipulate genetically, Desulfovibrio vulgaris Hildenborough, to generate a library of over 700 strains. Furthermore, we demonstrate how these modifications can be used for examining metabolic pathways, protein-protein interactions, and protein localization. The ubiquity of suicide constructs in gene replacement throughout biology suggests that this approach can be applied to engineer a broad range of species for a diverse array of systems biological applications and is amenable to high-throughput implementation.

Using population mixtures to optimize the utility of genomic databases: linkage disequilibrium and association study design in India.

When performing association studies in populations that have not been the focus of large-scale investigations of haplotype variation, it is often helpful to rely on genomic databases in other populations for study design and analysis - such as in the selection of tag SNPs and in the imputation of missing genotypes. One way of improving the use of these databases is to rely on a mixture of database samples that is similar to the population of interest, rather than using the single most similar database sample. We demonstrate the effectiveness of the mixture approach in the application of African, European, and East Asian HapMap samples for tag SNP selection in populations from India, a genetically intermediate region underrepresented in genomic studies of haplotype variation.

Hexavalent chromium reduction in Desulfovibrio vulgaris Hildenborough causes transitory inhibition of sulfate reduction and cell growth.

Desulfovibrio vulgaris Hildenborough is a well-studied sulfate reducer that can reduce heavy metals and radionuclides [e.g., Cr(VI) and U(VI)]. Cultures grown in a defined medium had a lag period of approximately 30 h when exposed to 0.05 mM Cr(VI). Substrate analyses revealed that although Cr(VI) was reduced within the first 5 h, growth was not observed for an additional 20 h. The growth lag could be explained by a decline in cell viability; however, during this time small amounts of lactate were still utilized without sulfate reduction or acetate formation. Approximately 40 h after Cr exposure (0.05 mM), sulfate reduction occurred concurrently with the accumulation of acetate. Similar amounts of hydrogen were produced by Cr-exposed cells compared to control cells, and lactate was not converted to glycogen during non-growth conditions. D. vulgaris cells treated with a reducing agent and then exposed to Cr(VI) still experienced a growth lag, but the addition of ascorbate at the time of Cr(VI) addition prevented the lag period. In addition, cells grown on pyruvate displayed more tolerance to Cr(VI) compared to lactate-grown cells. These results indicated that D. vulgaris utilized lactate during Cr(VI) exposure without the reduction of sulfate or production of acetate, and that ascorbate and pyruvate could protect D. vulgaris cells from Cr(VI)/Cr(III) toxicity.

Temporal transcriptomic analysis as Desulfovibrio vulgaris Hildenborough transitions into stationary phase during electron donor depletion.

Desulfovibrio vulgaris was cultivated in a defined medium, and biomass was sampled for approximately 70 h to characterize the shifts in gene expression as cells transitioned from the exponential to the stationary phase during electron donor depletion. In addition to temporal transcriptomics, total protein, carbohydrate, lactate, acetate, and sulfate levels were measured. The microarray data were examined for statistically significant expression changes, hierarchical cluster analysis, and promoter element prediction and were validated by quantitative PCR. As the cells transitioned from the exponential phase to the stationary phase, a majority of the down-expressed genes were involved in translation and transcription, and this trend continued at the remaining times. There were general increases in relative expression for intracellular trafficking and secretion, ion transport, and coenzyme metabolism as the cells entered the stationary phase. As expected, the DNA replication machinery was down-expressed, and the expression of genes involved in DNA repair increased during the stationary phase. Genes involved in amino acid acquisition, carbohydrate metabolism, energy production, and cell envelope biogenesis did not exhibit uniform transcriptional responses. Interestingly, most phage-related genes were up-expressed at the onset of the stationary phase. This result suggested that nutrient depletion may affect community dynamics and DNA transfer mechanisms of sulfate-reducing bacteria via the phage cycle. The putative feoAB system (in addition to other presumptive iron metabolism genes) was significantly up-expressed, and this suggested the possible importance of Fe2+ acquisition under metal-reducing conditions. The expression of a large subset of carbohydrate-related genes was altered, and the total cellular carbohydrate levels declined during the growth phase transition. Interestingly, the D. vulgaris genome does not contain a putative rpoS gene, a common attribute of the delta-Proteobacteria genomes sequenced to date, and the transcription profiles of other putative rpo genes were not significantly altered. Our results indicated that in addition to expected changes (e.g., energy conversion, protein turnover, translation, transcription, and DNA replication and repair), genes related to phage, stress response, carbohydrate flux, the outer envelope, and iron homeostasis played important roles as D. vulgaris cells experienced electron donor depletion.

Desulfovibrio desulfuricans G20 tetraheme cytochrome structure at 1.5 Angstrom and cytochrome interaction with metal complexes.

The structure of the type I tetraheme cytochrome c(3) from Desulfovibrio desulfuricans G20 was determined to 1.5 Angstrom by X-ray crystallography. In addition to the oxidized form, the structure of the molybdate-bound form of the protein was determined from oxidized crystals soaked in sodium molybdate. Only small structural shifts were obtained with metal binding, consistent with the remarkable structural stability of this protein. In vitro experiments with pure cytochrome showed that molybdate could oxidize the reduced cytochrome, although not as rapidly as U(VI) present as uranyl acetate. Alterations in the overall conformation and thermostability of the metal-oxidized protein were investigated by circular dichroism studies. Again, only small changes in protein structure were documented. The location of the molybdate ion near heme IV in the crystal structure suggested heme IV as the site of electron exit from the reduced cytochrome and implicated Lys14 and Lys56 in binding. Analysis of structurally conserved water molecules in type I cytochrome c(3) crystal structures identified interactions predicted to be important for protein stability and possibly for intramolecular electron transfer among heme molecules.

Global analysis of heat shock response in Desulfovibrio vulgaris Hildenborough.

Desulfovibrio vulgaris Hildenborough belongs to a class of sulfate-reducing bacteria (SRB) and is found ubiquitously in nature. Given the importance of SRB-mediated reduction for bioremediation of metal ion contaminants, ongoing research on D. vulgaris has been in the direction of elucidating regulatory mechanisms for this organism under a variety of stress conditions. This work presents a global view of this organism's response to elevated growth temperature using whole-cell transcriptomics and proteomics tools. Transcriptional response (1.7-fold change or greater; Z >/= 1.5) ranged from 1,135 genes at 15 min to 1,463 genes at 120 min for a temperature up-shift of 13 degrees C from a growth temperature of 37 degrees C for this organism and suggested both direct and indirect modes of heat sensing. Clusters of orthologous group categories that were significantly affected included posttranslational modifications; protein turnover and chaperones (up-regulated); energy production and conversion (down-regulated), nucleotide transport, metabolism (down-regulated), and translation; ribosomal structure; and biogenesis (down-regulated). Analysis of the genome sequence revealed the presence of features of both negative and positive regulation which included the CIRCE element and promoter sequences corresponding to the alternate sigma factors sigma(32) and sigma(54). While mechanisms of heat shock control for some genes appeared to coincide with those established for Escherichia coli and Bacillus subtilis, the presence of unique control schemes for several other genes was also evident. Analysis of protein expression levels using differential in-gel electrophoresis suggested good agreement with transcriptional profiles of several heat shock proteins, including DnaK (DVU0811), HtpG (DVU2643), HtrA (DVU1468), and AhpC (DVU2247). The proteomics study also suggested the possibility of posttranslational modifications in the chaperones DnaK, AhpC, GroES (DVU1977), and GroEL (DVU1976) and also several periplasmic ABC transporters.

Insights from linked single nucleotide polymorphisms: what we can learn from linkage disequilibrium.

New methods for analyzing sequence polymorphism data have uncovered some striking patterns of linkage disequilibrium in both humans and fruitflies. These methods have revealed examples where the observed amount of linkage disequilibrium is either much more or much less than expected, and have led to advances in our understanding of the forces that affect naturally occurring genetic variation. With the recent explosion of sequence polymorphism data, the prospects for further progress from these methods are quite promising.

Gene conversion and different population histories may explain the contrast between polymorphism and linkage disequilibrium levels.

To characterize linkage disequilibrium (LD) levels in human populations, we have analyzed 10 independent noncoding segments in three population samples from the major ethnic groups--that is, Africans, Asians, and Europeans. Descriptive statistics show that LD decays much faster in the African samples than in the non-African ones. With the assumption of an equilibrium model, we estimated the population crossing-over parameter (4N(e)r(bp), where N(e) is the effective population size and r(bp) is the crossing-over rate per generation between adjacent base pairs) in the presence of gene conversion. In the African sample, LD and polymorphism levels lead to similar estimates of effective population size, as expected under an equilibrium model. Conversely, in both non-African samples, LD levels suggest a smaller effective population size than that implied by polymorphism levels. This observation is paralleled by significant departures from an equilibrium model in the spectrum of allele frequencies of the non-African samples. Besides ruling out the possibility that non-African populations are at equilibrium, these results suggest different demographic history (temporal and spatial) of these groups. Interestingly, the African sample fits the expectations of an equilibrium model based on polymorphism and divergence levels and on frequency spectrum. For this sample, the estimated ratio of gene conversion to crossing-over rates is 7.3 for a mean tract length of 500 bp, suggesting that gene conversion may be more frequent than previously thought. These findings imply that disease-association studies will require a much denser map of polymorphic sites in African than in non-African populations.

Why is there so little intragenic linkage disequilibrium in humans?

The efficient design of association mapping studies relies on a knowledge of the rate of decay of linkage disequilibrium with distance. This rate depends on the population recombination rate, C. An estimate of C for humans is usually obtained from a comparison of physical and genetic maps, assuming an effective population size of approximately 10(4). We demonstrate that under both a constant population size model and a model of long-term exponential growth, there is evidence for more recombination in polymorphism data than is expected from this estimate. An important contribution of gene conversion to meiotic recombination helps to explain our observation, but does not appear to be sufficient. The occurrence of multiple hits at CpG sites and the presence of population structure are not explanations.

Recombination and the frequency spectrum in Drosophila melanogaster and Drosophila simulans.

Most "tests of neutrality" assess whether particular data sets depart from the predictions of a standard neutral model with no recombination. For Drosophila, where nuclear polymorphism data routinely show evidence of genetic exchange, the assumption of no recombination is often unrealistic. In addition, while conservative, this assumption is made at the cost of a great loss in power. Perhaps as a result, tests of the frequency spectrum based on zero recombination suggest an adequate fit of Drosophila polymorphism data to the predictions of the standard neutral model. Here, we analyze the frequency spectrum of a large number of loci in Drosophila melanogaster and D. simulans using two summary statistics. We use an estimate of the population recombination rate based on a laboratory estimate of the rate of crossing over per physical length and an estimate of the species' effective population size. In contrast to previous studies, we find that roughly half of the loci depart from the predictions of the standard neutral model. The extent of the departure depends on the exact recombination rate, but the global pattern that emerges is robust. Interestingly, these departures from neutral expectations are not unidirectional. The large variance in outcomes may be due to a complex demographic history and inconsistent sampling, or to the pervasive action of natural selection.

When did the human population size start increasing?

We analyze the frequency spectra of all available human nuclear sequence data sets by using a model of constant population size followed by exponential growth. Parameters of growth (more extreme than or) comparable to what has been suggested from mtDNA data can be rejected for 6 out of the 10 largest data sets. When the data are separated into African and non-African samples, a constant size no-growth model can be rejected for 4 out of 8 non-African samples. Long-term growth (i.e., starting 50-100 kya) can be rejected for 2 out of 8 African samples and 5 out of 8 non-African ones. Under more complex demographic models, including a bottleneck or population subdivision, more of the data are compatible with long-term growth. One problem with the data used here is that a subset of loci may reflect the action of natural selection as well as of demography. It remains possible that the correct demographic model is one of constant population size followed by long-term growth but that at several loci the demographic signature has been obscured by balancing or diversifying selection. However, it is not clear that the data at these loci are consistent with a simple model of balancing selection; more complicated selective alternatives cannot be tested unless they are made explicit. An alternative explanation is that population size growth is more recent (e.g., upper Paleolithic) and that some of the loci have experienced recent directional selection. Given the available data, the latter hypothesis seems more likely.

Detecting ancient admixture in humans using sequence polymorphism data.

A debate of long-standing interest in human evolution centers around whether archaic human populations (such as the Neanderthals) have contributed to the modern gene pool. A model of ancient population structure with recent mixing is introduced, and it is determined how much information (i.e., sequence data from how many unlinked nuclear loci) would be necessary to distinguish between different demographic scenarios. It is found that approximately 50-100 loci are necessary if plausible parameter estimates are used. There are not enough data available at the present to support either the "single origin" or the "multiregional" model of modern human evolution. However, this information should be available in a few years.

Cytochrome c(3) mutants of Desulfovibrio desulfuricans.

To explore the physiological role of tetraheme cytochrome c(3) in the sulfate-reducing bacterium Desulfovibrio desulfuricans G20, the gene encoding the preapoprotein was cloned, sequenced, and mutated by plasmid insertion. The physical analysis of the DNA from the strain carrying the integrated plasmid showed that the insertion was successful. The growth rate of the mutant on lactate with sulfate was comparable to that of the wild type; however, mutant cultures did not achieve the same cell densities. Pyruvate, the oxidation product of lactate, served as a poor electron source for the mutant. Unexpectedly, the mutant was able to grow on hydrogen-sulfate medium. These data support a role for tetraheme cytochrome c(3) in the electron transport pathway from pyruvate to sulfate or sulfite in D. desulfuricans G20.

A comparison of estimators of the population recombination rate.

Three new estimators of the population recombination rate C = 4Nr are introduced. These estimators summarize the data using the number of distinct haplotypes and the estimated minimum number of recombination events, then calculate the value of C that maximizes the likelihood of obtaining the summarized data. They are compared with a number of previously proposed estimators of the recombination rate. One of the newly proposed estimators is generally better than the others for the parameter values considered here, while the three programs that calculate maximum-likelihood estimates give conflicting results.

Unusual haplotype structure at the proximal breakpoint of In(2L)t in a natural population of Drosophila melanogaster.

The existence of temporally stable frequency clines for In(2L)t in natural populations of Drosophila melanogaster suggests a role for selection in the maintenance of this polymorphism. We have collected nucleotide polymorphism data from the proximal breakpoint junction regions of In(2L)t to infer its evolutionary history. The finding of a novel LINE-like element near the In(2L)t breakpoint junction in sampled inverted chromosomes supports a transposable element-mediated origin for this inversion. An analysis of nucleotide variation in a Costa Rican population sample of standard and inverted chromosomes indicates a unique and relatively recent origin for In(2L)t. Additional In(2L)t alleles from three geographically diverse populations reveal no detectable geographic differentiation. Low levels of In(2L)t nucleotide polymorphism suggest a recent increase in the inversion's frequency in tropical populations. An unusual feature of our sample of standard alleles is a marked heterogeneity in levels of linkage disequilibrium among polymorphic sites across the breakpoint region. We introduce a test of neutral equilibrium haplotype structure that corrects both for multiple tests and for an arbitrarily chosen window size. It reveals that an approximately 1.4-kb region immediately spanning the breakpoint has fewer haplotypes than expected under the neutral model, given the expected level of recombination in this genomic region. Certain features of our data suggest that the unusual pattern in standard chromosomes is the product of selection rather than demography.

Genealogies and weak purifying selection.

The assumption that selection alters the genealogical tree of a sample of alleles from a population relative to the neutral expectation underlies several "tests of neutrality." Two recent papers have studied the effect of purifying selection; their suggestive but incomplete results indicate that, in the single site case, the shape of a gene genealogy for a locus may differ only from the neutral expectation. We verify this finding for weak selection using the "ancestral selection graph." We consider a wider range of models, including both a four-allele single-site model and an infinite-sites model. Our results confirm the previous claim for the symmetric-mutation single site model. We emphasize, however, that a neutral-seeming genealogy is consistent with detectable effects of selection on the distribution of allele frequences within the sample. With selection operating, the information about a sample cannot be reduced to the genealogy. As a result, a distinction needs to be made between the selected sites themselves, for which the genealogy offers insufficient information, and linked neutral variation. This distinction seems to have been overlooked in previous papers, yet it has significant implications for the interpretation of data on DNA sequence variation. In particular, it predicts that under purifying selection, the frequency spectrum of neutral mutations will not reflect the skew toward rare polymorphisms at replacement sites even if there is no recombination between them. We caution, however, that the effect of weak selection on the genealogy is specific to the model; a (more realistic) model of multiple linked sites could lead to a more distorted genealogy than is observed for a single site.

New shuttle vectors for the introduction of cloned DNA in Desulfovibrio.

The pBG1 replicon from the cryptic plasmid of Desulfovibrio desulfuricans G100A was inserted into pTZ18U derivatives to generate a new family of shuttle vectors. These plasmids are stable both in Escherichia coli and in Desulfovibrio, they present a large number of unique restriction sites, and colonies of recombinant clones can be identified by blue/white screening in E. coli. The pBMC, pBMK, and pBMS series carry the cat, npt, or strAB genes as selectable markers, respectively. The pBMC6, pBMK6, and pBMS6 plasmids can be introduced both in D. desulfuricans and in Desulfovibrio fructosovorans by electrotransformation, and the pBMC7, pBMK7, and pBMS7 plasmids contain additional mobilization functions which makes them suitable for conjugation.