Search

Our ability to predict, control, or design biological function is fundamentally limited by poorly annotated gene function. This can be particularly challenging in non-model systems. Accordingly, there is motivation for new high-throughput methods for accurate functional annotation. Here, we used complementation of auxotrophs and DNA barcode sequencing (Coaux-Seq) to enable high-throughput characterization of protein function. Fragment libraries from eleven genetically diverse bacteria were tested in twenty different auxotrophic strains of Escherichia coli to identify genes that complement missing biochemical activity. We recovered 41% of expected hits, with effectiveness ranging per source genome, and observed success even with distant E. coli relatives like Bacillus subtilis and Bacteroides thetaiotaomicron. Coaux-Seq provided the first experimental validation for 53 proteins, of which 11 are less than 40% identical to an experimentally characterized protein. Among the unexpected function identified was a sulfate uptake transporter, an O-succinylhomoserine sulfhydrylase for methionine synthesis, and an aminotransferase. We also identified instances of cross-feeding wherein protein overexpression and nearby non-auxotrophic strains enabled growth. Altogether, Coaux-Seq's utility is demonstrated, with future applications in ecology, health, and engineering.

Studying plant–microbe–soil interactions is challenging due to their high complexity and variability in natural ecosystems. While fabricated ecosystems provide opportunities to recapitulate aspects of these systems in reduced complexity and controlled environments, inoculation can be a significant source of variation. To tackle this, we evaluated how different bacterial inoculation practices and plant harvesting time points affect the reproducibility of a synthetic microbial community (SynCom) in association with the model grass Brachypodium distachyon. We tested three microbial inoculation practices, seed inoculation, transplant inoculation, and seedling inoculation, and two harvesting points, early (14-day-old plants) and late (21 days postinoculation). We grew our plants and bacterial strains in sterile devices (EcoFABs) and characterized the microbial community from root, rhizosphere, and sand using 16S ribosomal RNA gene sequencing. The results showed that inoculation practices significantly affected the rhizosphere microbial community when harvesting at an early time point but not at the late time point. As the SynCom showed a persistent association with B. distachyon at 21 days postinoculation regardless of inoculation practices, we assessed the reproducibility of each inoculation method and found that transplant inoculation showed the highest reproducibility. Moreover, plant biomass was not adversely affected by transplant inoculation treatment. We concluded that bacterial inoculation while transplanting coupled with a later harvesting time point gives the most reproducible microbial community in the EcoFAB−B. distachyon−SynCom fabricated ecosystem. We recommend this method as a standardized protocol for use with fabricated ecosystem experimental systems.

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.

Exometabolomics is an approach to assess how microorganisms alter, or react to their environments through the depletion and production of metabolites. It allows the examination of how soil microbes transform the small molecule metabolites within their environment, which can be used to study resource competition and cross-feeding. This approach is most powerful when used with defined media that enable tracking of all metabolites. However, microbial growth media have traditionally been developed for the isolation and growth of microorganisms but not metabolite utilization profiling through Liquid Chromatography Tandem Mass Spectrometry (LC-MS/MS). Here, we describe the construction of a defined medium, the Northen Lab Defined Medium (NLDM), that not only supports the growth of diverse soil bacteria but also is defined and therefore suited for exometabolomic experiments. Metabolites included in NLDM were selected based on their presence in R2A medium and soil, elemental stoichiometry requirements, as well as knowledge of metabolite usage by different bacteria. We found that NLDM supported the growth of 108 of the 110 phylogenetically diverse (spanning 36 different families) soil bacterial isolates tested and all of its metabolites were trackable through LC-MS/MS analysis. These results demonstrate the viability and utility of the constructed NLDM medium for growing and characterizing diverse microbial isolates and communities.

Bacterial communities in water, soil, and humans play an essential role in environmental ecology and human health. PCR-based amplicon analysis, such as 16S rRNA sequencing, is a fundamental tool for quantifying and studying microbial composition, dynamics, and interactions. However, given the complexity of microbial communities, a substantial number of samples becomes necessary for analyses that parse the factors that determine microbial composition. A common bottleneck in performing these kinds of experiments is genomic DNA (gDNA) extraction, which is time-consuming, expensive, and often biased based on the types of species present. Direct PCR method is a potentially simpler and more accurate alternative to gDNA extraction methods that do not require the intervening purification step. In this study, we evaluated three variations of direct PCR methods using diverse heterogeneous bacterial cultures, including both Gram-positive and Gram-negative species, ZymoBIOMICS microbial community standards, and groundwater. By comparing direct PCR methods with DNeasy Blood and Tissue Kits for microbial isolates and DNeasy PowerSoil Kits for microbial communities, we found that a specific variant of the direct PCR method exhibits an overall efficiency comparable to that of the conventional DNeasy PowerSoil protocol in the circumstances we tested. We also found that the method showed higher efficiency for extracting gDNA from the Gram-negative strains compared to DNeasy Blood and Tissue protocol. This direct PCR method is 1,600 times less expensive ($0.34 for 96 samples) and 10 times simpler (15 min hands-on time for 96 samples) than the DNeasy PowerSoil protocol. The direct PCR method can also be fully automated and is compatible with small-volume samples, thereby permitting scaling of samples and replicates needed to support high-throughput large-scale bacterial community analysis. IMPORTANCE Understanding bacterial interactions and assembly in complex microbial communities using 16S rRNA sequencing normally requires a large experimental load. However, the current DNA extraction methods, including cell disruption and genomic DNA purification, are normally biased, costly, time-consuming, labor-intensive, and not amenable to miniaturization by droplets or 1,536-well plates due to the significant DNA loss during the purification step for tiny-volume and low-cell-density samples. A direct PCR method could potentially solve these problems. In this study, we developed a direct PCR method which exhibits similar efficiency as the widely used method, the DNeasy PowerSoil protocol, while being 1,600 times less expensive and 10 times faster to execute. This simple, cost-effective, and automation-friendly direct-PCR-based 16S rRNA sequencing method allows us to study the dynamics, microbial interaction, and assembly of various microbial communities in a high-throughput fashion.

Abstract Background In Drosophila, the Enhancer of split complex (E(spl)-C) comprises 11 bHLH and Bearded genes that function during Notch signaling to repress proneural identity in the developing peripheral nervous system. Comparison with other insects indicates that the basal state for Diptera is a single bHLH and Bearded homolog and that the expansion of the gene complex occurred in the lineage leading to Drosophila. However, comparative genomic data from other fly species that would elucidate the origin and sequence of gene duplication for the complex is lacking. Therefore, in order to examine the evolutionary history of the complex within Diptera, we reconstructed, using several fosmid clones, the entire E(spl)-complex in the stalk-eyed fly, Teleopsis dalmanni and collected additional homologs of E(spl)-C genes from searches of dipteran EST databases and the Glossina morsitans genome assembly. Results Comparison of the Teleopsis E(spl)-C gene organization with Drosophila indicates complete conservation in gene number and orientation between the species except that T. dalmanni contains a duplicated copy of E(spl)m5 that is not present in Drosophila. Phylogenetic analysis of E(spl)-complex bHLH and Bearded genes for several dipteran species clearly demonstrates that all members of the complex were present prior to the diversification of schizophoran flies. Comparison of upstream regulatory elements and 3' UTR domains between the species also reveals strong conservation for many of the genes and identifies several novel characteristics of E(spl)-C regulatory evolution including the discovery of a previously unidentified, highly conserved SPS+A domain between E(spl)mγ and E(spl)mβ. Conclusion Identifying the phylogenetic origin of E(spl)-C genes and their associated regulatory DNA is essential to understanding the functional significance of this well-studied gene complex. Results from this study provide numerous insights into the evolutionary history of the complex and will help refine the focus of studies examining the adaptive consequences of this gene expansion.

Abstract Background Despite considerable progress in our understanding of land plant phylogeny, several nodes in the green tree of life remain poorly resolved. Furthermore, the bulk of currently available data come from only a subset of major land plant clades. Here we examine early land plant evolution using complete plastome sequences including two previously unexamined and phylogenetically critical lineages. To better understand the evolution of land plants and their plastomes, we examined aligned nucleotide sequences, indels, gene and nucleotide composition, inversions, and gene order at the boundaries of the inverted repeats. Results We present the plastome sequences of Equisetum arvense, a horsetail, and of Isoetes flaccida, a heterosporous lycophyte. Phylogenetic analysis of aligned nucleotides from 49 plastome genes from 43 taxa supported monophyly for the following clades: embryophytes (land plants), lycophytes, monilophytes (leptosporangiate ferns + Angiopteris evecta + Psilotum nudum + Equisetum arvense), and seed plants. Resolution among the four monilophyte lineages remained moderate, although nucleotide analyses suggested that P. nudum and E. arvense form a clade sister to A. evecta + leptosporangiate ferns. Results from phylogenetic analyses of nucleotides were consistent with the distribution of plastome gene rearrangements and with analysis of sequence gaps resulting from insertions and deletions (indels). We found one new indel and an inversion of a block of genes that unites the monilophytes. Conclusions Monophyly of monilophytes has been disputed on the basis of morphological and fossil evidence. In the context of a broad sampling of land plant data we find several new pieces of evidence for monilophyte monophyly. Results from this study demonstrate resolution among the four monilophytes lineages, albeit with moderate support; we posit a clade consisting of Equisetaceae and Psilotaceae that is sister to the "true ferns," including Marattiaceae.

Abstract Background Tortula ruralis, a widely distributed species in the moss family Pottiaceae, is increasingly used as a model organism for the study of desiccation tolerance and mechanisms of cellular repair. In this paper, we present the chloroplast genome sequence of T. ruralis, only the second published chloroplast genome for a moss, and the first for a vegetatively desiccation-tolerant plant. Results The Tortula chloroplast genome is ~123,500 bp, and differs in a number of ways from that of Physcomitrella patens, the first published moss chloroplast genome. For example, Tortula lacks the ~71 kb inversion found in the large single copy region of the Physcomitrella genome and other members of the Funariales. Also, the Tortula chloroplast genome lacks petN, a gene found in all known land plant plastid genomes. In addition, an unusual case of nucleotide polymorphism was discovered. Conclusions Although the chloroplast genome of Tortula ruralis differs from that of the only other sequenced moss, Physcomitrella patens, we have yet to determine the biological significance of the differences. The polymorphisms we have uncovered in the sequencing of the genome offer a rare possibility (for mosses) of the generation of DNA markers for fine-level phylogenetic studies, or to investigate individual variation within populations.

Abstract Background Welwitschia mirabilis is the only extant member of the family Welwitschiaceae, one of three lineages of gnetophytes, an enigmatic group of gymnosperms variously allied with flowering plants or conifers. Limited sequence data and rapid divergence rates have precluded consensus on the evolutionary placement of gnetophytes based on molecular characters. Here we report on the first complete gnetophyte chloroplast genome sequence, from Welwitschia mirabilis, as well as analyses on divergence rates of protein-coding genes, comparisons of gene content and order, and phylogenetic implications. Results The chloroplast genome of Welwitschia mirabilis [GenBank: EU342371] is comprised of 119,726 base pairs and exhibits large and small single copy regions and two copies of the large inverted repeat (IR). Only 101 unique gene species are encoded. The Welwitschia plastome is the most compact photosynthetic land plant plastome sequenced to date; 66% of the sequence codes for product. The genome also exhibits a slightly expanded IR, a minimum of 9 inversions that modify gene order, and 19 genes that are lost or present as pseudogenes. Phylogenetic analyses, including one representative of each extant seed plant lineage and based on 57 concatenated protein-coding sequences, place Welwitschia at the base of all seed plants (distance, maximum parsimony) or as the sister to Pinus (the only conifer representative) in a monophyletic gymnosperm clade (maximum likelihood, bayesian). Relative rate tests on these gene sequences show the Welwitschia sequences to be evolving at faster rates than other seed plants. For these genes individually, a comparison of average pairwise distances indicates that relative divergence in Welwitschia ranges from amounts about equal to other seed plants to amounts almost three times greater than the average for non-gnetophyte seed plants. Conclusion Although the basic organization of the Welwitschia plastome is typical, its compactness, gene content and high nucleotide divergence rates are atypical. The current lack of additional conifer plastome sequences precludes any discrimination between the gnetifer and gnepine hypotheses of seed plant relationships. However, both phylogenetic analyses and shared genome features identified here are consistent with either of the hypotheses that link gnetophytes with conifers, but are inconsistent with the anthophyte hypothesis.

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.

Abstract Background The genus Cuscuta L. (Convolvulaceae), commonly known as dodders, are epiphytic vines that invade the stems of their host with haustorial feeding structures at the points of contact. Although they lack expanded leaves, some species are noticeably chlorophyllous, especially as seedlings and in maturing fruits. Some species are reported as crop pests of worldwide distribution, whereas others are extremely rare and have local distributions and apparent niche specificity. A strong phylogenetic framework for this large genus is essential to understand the interesting ecological, morphological and molecular phenomena that occur within these parasites in an evolutionary context. Results Here we present a well-supported phylogeny of Cuscuta using sequences of the nuclear ribosomal internal transcribed spacer and plastid rps2, rbcL and matK from representatives across most of the taxonomic diversity of the genus. We use the phylogeny to interpret morphological and plastid genome evolution within the genus. At least three currently recognized taxonomic sections are not monophyletic and subgenus Cuscuta is unequivocally paraphyletic. Plastid genes are extremely variable with regards to evolutionary constraint, with rbcL exhibiting even higher levels of purifying selection in Cuscuta than photosynthetic relatives. Nuclear genome size is highly variable within Cuscuta, particularly within subgenus Grammica, and in some cases may indicate the existence of cryptic species in this large clade of morphologically similar species. Conclusion Some morphological characters traditionally used to define major taxonomic splits within Cuscuta are homoplastic and are of limited use in defining true evolutionary groups. Chloroplast genome evolution seems to have evolved in a punctuated fashion, with episodes of loss involving suites of genes or tRNAs followed by stabilization of gene content in major clades. Nearly all species of Cuscuta retain some photosynthetic ability, most likely for nutrient apportionment to their seeds, while complete loss of photosynthesis and possible loss of the entire chloroplast genome is limited to a single small clade of outcrossing species found primarily in western South America.

Abstract Background Plastid genome content and protein sequence are highly conserved across land plants and their closest algal relatives. Parasitic plants, which obtain some or all of their nutrition through an attachment to a host plant, are often a striking exception. Heterotrophy can lead to relaxed constraint on some plastid genes or even total gene loss. We sequenced plastid genomes of two species in the parasitic genus Cuscuta along with a non-parasitic relative, Ipomoea purpurea, to investigate changes in the plastid genome that may result from transition to the parasitic lifestyle. Results Aside from loss of all ndh genes, Cuscuta exaltata retains photosynthetic and photorespiratory genes that evolve under strong selective constraint. Cuscuta obtusiflora has incurred substantially more change to its plastid genome, including loss of all genes for the plastid-encoded RNA polymerase. Despite extensive change in gene content and greatly increased rate of overall nucleotide substitution, C. obtusiflora also retains all photosynthetic and photorespiratory genes with only one minor exception. Conclusion Although Epifagus virginiana, the only other parasitic plant with its plastid genome sequenced to date, has lost a largely overlapping set of transfer-RNA and ribosomal genes as Cuscuta, it has lost all genes related to photosynthesis and maintains a set of genes which are among the most divergent in Cuscuta. Analyses demonstrate photosynthetic genes are under the highest constraint of any genes within the plastid genomes of Cuscuta, indicating a function involving RuBisCo and electron transport through photosystems is still the primary reason for retention of the plastid genome in these species.

Abstract Background Teleost fish have seven paralogous clusters of Hox genes stemming from two complete genome duplications early in vertebrate evolution, and an additional genome duplication during the evolution of ray-finned fish, followed by the secondary loss of one cluster. Gene duplications on the one hand, and the evolution of regulatory sequences on the other, are thought to be among the most important mechanisms for the evolution of new gene functions. Cichlid fish, the largest family of vertebrates with about 2500 species, are famous examples of speciation and morphological diversity. Since this diversity could be based on regulatory changes, we chose to study the coding as well as putative regulatory regions of their Hox clusters within a comparative genomic framework. Results We sequenced and characterized all seven Hox clusters of Astatotilapia burtoni, a haplochromine cichlid fish. Comparative analyses with data from other teleost fish such as zebrafish, two species of pufferfish, stickleback and medaka were performed. We traced losses of genes and microRNAs of Hox clusters, the medaka lineage seems to have lost more microRNAs than the other fish lineages. We found that each teleost genome studied so far has a unique set of Hox genes. The hoxb7a gene was lost independently several times during teleost evolution, the most recent event being within the radiation of East African cichlid fish. The conserved non-coding sequences (CNS) encompass a surprisingly large part of the clusters, especially in the HoxAa, HoxCa, and HoxDa clusters. Across all clusters, we observe a trend towards an increased content of CNS towards the anterior end. Conclusion The gene content of Hox clusters in teleost fishes is more variable than expected, with each species studied so far having a different set. Although the highest loss rate of Hox genes occurred immediately after whole genome duplications, our analyses showed that gene loss continued and is still ongoing in all teleost lineages. Along with the gene content, the CNS content also varies across clusters. The excess of CNS at the anterior end of clusters could imply a stronger conservation of anterior expression patters than those towards more posterior areas of the embryo.

Abstract Background The magnoliids with four orders, 19 families, and 8,500 species represent one of the largest clades of early diverging angiosperms. Although several recent angiosperm phylogenetic analyses supported the monophyly of magnoliids and suggested relationships among the orders, the limited number of genes examined resulted in only weak support, and these issues remain controversial. Furthermore, considerable incongruence resulted in phylogenetic reconstructions supporting three different sets of relationships among magnoliids and the two large angiosperm clades, monocots and eudicots. We sequenced the plastid genomes of three magnoliids, Drimys (Canellales), Liriodendron (Magnoliales), and Piper (Piperales), and used these data in combination with 32 other angiosperm plastid genomes to assess phylogenetic relationships among magnoliids and to examine patterns of variation of GC content. Results The Drimys, Liriodendron, and Piper plastid genomes are very similar in size at 160,604, 159,886 bp, and 160,624 bp, respectively. Gene content and order are nearly identical to many other unrearranged angiosperm plastid genomes, including Calycanthus, the other published magnoliid genome. Overall GC content ranges from 34–39%, and coding regions have a substantially higher GC content than non-coding regions. Among protein-coding genes, GC content varies by codon position with 1st codon > 2nd codon > 3rd codon, and it varies by functional group with photosynthetic genes having the highest percentage and NADH genes the lowest. Phylogenetic analyses using parsimony and likelihood methods and sequences of 61 protein-coding genes provided strong support for the monophyly of magnoliids and two strongly supported groups were identified, the Canellales/Piperales and the Laurales/Magnoliales. Strong support is reported for monocots and eudicots as sister clades with magnoliids diverging before the monocot-eudicot split. The trees also provided moderate or strong support for the position of Amborella as sister to a clade including all other angiosperms. Conclusion Evolutionary comparisons of three new magnoliid plastid genome sequences, combined with other published angiosperm genomes, confirm that GC content is unevenly distributed across the genome by location, codon position, and functional group. Furthermore, phylogenetic analyses provide the strongest support so far for the hypothesis that the magnoliids are sister to a large clade that includes both monocots and eudicots.

ABSTRACTSulfate-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.

Respiratory and catabolic genes are differentially distributed across microbial genomes. Thus, specific carbon sources may favor different respiratory processes. We profiled the influence of 94 carbon sources on the end products of nitrate respiration in microbial enrichment cultures from diverse terrestrial environments. We found that some carbon sources consistently favor dissimilatory nitrate reduction to ammonium (DNRA/nitrate ammonification) while other carbon sources favor nitrite accumulation or denitrification. For an enrichment culture from aquatic sediment, we sequenced the genomes of the most abundant strains, matched these genomes to 16S rDNA exact sequence variants (ESVs), and used 16S rDNA amplicon sequencing to track the differential enrichment of functionally distinct ESVs on different carbon sources. We found that changes in the abundances of strains with different genetic potentials for nitrite accumulation, DNRA or denitrification were correlated with the nitrite or ammonium concentrations in the enrichment cultures recovered on different carbon sources. Specifically, we found that either L-sorbose or D-cellobiose enriched for a Klebsiella nitrite accumulator, other sugars enriched for an Escherichia nitrate ammonifier, and citrate or formate enriched for a Pseudomonas denitrifier and a Sulfurospirillum nitrate ammonifier. Our results add important nuance to the current paradigm that higher concentrations of carbon will always favor DNRA over denitrification or nitrite accumulation, and we propose that, in some cases, carbon composition can be as important as carbon concentration in determining nitrate respiratory end products. Furthermore, our approach can be extended to other environments and metabolisms to characterize how selective parameters influence microbial community composition, gene content, and function.

[This corrects the article DOI: 10.1371/journal.pgen.1007147.].

In many environments, toxic compounds restrict which microorganisms persist. However, in complex mixtures of inhibitory compounds, it is challenging to determine which specific compounds cause changes in abundance and prevent some microorganisms from growing. We focused on a contaminated aquifer in Oak Ridge, Tennessee, USA that has large gradients of pH and widely varying concentrations of uranium, nitrate, and many other inorganic ions. In the most contaminated wells, the microbial community is enriched in the Rhodanobacter genus. Rhodanobacter abundance is positively correlated with low pH and high concentrations of uranium and 13 other ions and we sought to determine which of these ions are selective pressures that favor the growth of Rhodanobacter over other taxa. Of these ions, low pH and high UO22+, Mn2+, Al3+, Cd2+, Zn2+, Co2+, and Ni2+ are both (a) selectively inhibitory of a Pseudomonas isolate from an uncontaminated well vs. a Rhodanobacter isolate from a contaminated well, and (b) reach toxic concentrations (for the Pseudomonas isolate) in the Rhodanobacter-dominated wells. We used mixtures of ions to simulate the groundwater conditions in the most contaminated wells and verified that few isolates aside from Rhodanobacter can tolerate these eight ions. These results clarify which ions are likely causal factors that impact the microbial community at this field site and are not merely correlated with taxonomic shifts. Furthermore, our general high-throughput approach can be applied to other environments, isolates, and conditions to systematically help identify selective pressures on microbial communities.

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.

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.

Biosouring results from production of H2S by sulfate-reducing microorganisms (SRMs) in oil reservoirs. H2S is toxic, corrosive, and explosive, and as such, represents a significant threat to personnel, production facilities, and transportation pipelines. Since typical oil reservoir pressures can range from 10 to 50 MPa, understanding the role that pressure plays in SRM metabolism is important to improving souring containment strategies. To explore the impact of pressure, we grew an oil-field SRM isolate, Desulfovibrio alaskensis G20, under a range of pressures (0.1-14 MPa) at 30°C. The observed microbial growth rate was an inverse function of pressure with an associated slight reduction in sulfate and lactate consumption rate. Competitive fitness experiments with randomly bar-coded transposon mutant library sequencing (RB-TnSeq) identified several genes associated with flagellar biosynthesis and assembly that were important at high pressure. The fitness impact of specific genes was confirmed using individual transposon mutants. Confocal microscopy revealed that enhanced cell aggregation occurs at later stages of growth under pressure. We also assessed the effect of pressure on SRM inhibitor potency. Dose-response experiments showed a twofold decrease in the sensitivity of D. alaskensis to the antibiotic chloramphenicol at 14 MPa. Fortuitously, pressure had no significant influence on the inhibitory potency of the common souring controlling agent nitrate, or the emerging SRM inhibitors perchlorate, monofluorophosphate, or zinc pyrithione. Our findings improve the conceptual model of microbial sulfate reduction in high-pressure environments and the influence of pressure on souring inhibitor efficacy.

The prevalence of lipids devoid of phosphorus suggests that the availability of phosphorus limits microbial growth and activity in many anoxic, stratified environments. To better understand the response of anaerobic bacteria to phosphate limitation and starvation, this study combines microscopic and lipid analyses with the measurements of fitness of pooled barcoded transposon mutants of the model sulfate reducing bacterium Desulfovibrio alaskensis G20. Phosphate-limited G20 has lower growth rates and replaces more than 90% of its membrane phospholipids by a mixture of monoglycosyl diacylglycerol (MGDG), glycuronic acid diacylglycerol (GADG) and ornithine lipids, lacks polyphosphate granules, and synthesizes other cellular inclusions. Analyses of pooled and individual mutants reveal the importance of the high-affinity phosphate transport system (the Pst system), PhoR, and glycolipid and ornithine lipid synthases during phosphate limitation. The phosphate-dependent synthesis of MGDG in G20 and the widespread occurrence of the MGDG/GADG synthase among sulfate reducing ∂-Proteobacteria implicate these microbes in the production of abundant MGDG in anaerobic environments where the concentrations of phosphate are lower than 10 μM. Numerous predicted changes in the composition of the cell envelope and systems involved in transport, maintenance of cytoplasmic redox potential, central metabolism and regulatory pathways also suggest an impact of phosphate limitation on the susceptibility of sulfate reducing bacteria to other anthropogenic or environmental stresses.

Cupriavidus basilensis 4G11 was isolated from groundwater at the Oak Ridge Field Research Center (FRC) site. Here, we report the complete genome sequence and annotation of Cupriavidus basilensis 4G11. The genome contains 8,421,483 bp, 7,661 predicted protein-coding genes, and a total GC content of 64.4%.

We investigated perchlorate (ClO(4)(-)) and chlorate (ClO(3)(-)) (collectively (per)chlorate) in comparison with nitrate as potential inhibitors of sulfide (H(2)S) production by mesophilic sulfate-reducing microorganisms (SRMs). We demonstrate the specificity and potency of (per)chlorate as direct SRM inhibitors in both pure cultures and undefined sulfidogenic communities. We demonstrate that (per)chlorate and nitrate are antagonistic inhibitors and resistance is cross-inducible implying that these compounds share at least one common mechanism of resistance. Using tagged-transposon pools we identified genes responsible for sensitivity and resistance in Desulfovibrio alaskensis G20. We found that mutants in Dde_2702 (Rex), a repressor of the central sulfate-reduction pathway were resistant to both (per)chlorate and nitrate. In general, Rex derepresses its regulon in response to increasing intracellular NADH:NAD(+) ratios. In cells in which respiratory sulfate reduction is inhibited, NADH:NAD(+) ratios should increase leading to derepression of the sulfate-reduction pathway. In support of this, in (per)chlorate or nitrate-stressed wild-type G20 we observed higher NADH:NAD(+) ratios, increased transcripts and increased peptide counts for genes in the core Rex regulon. We conclude that one mode of (per)chlorate and nitrate toxicity is as direct inhibitors of the central sulfate-reduction pathway. Our results demonstrate that (per)chlorate are more potent inhibitors than nitrate in both pure cultures and communities, implying that they represent an attractive alternative for controlling sulfidogenesis in industrial ecosystems. Of these, perchlorate offers better application logistics because of its inhibitory potency, solubility, relative chemical stability, low affinity for mineral cations and high mobility in environmental systems.

Despite the environmental and economic cost of microbial sulfidogenesis in industrial operations, few compounds are known as selective inhibitors of respiratory sulfate reducing microorganisms (SRM), and no study has systematically and quantitatively evaluated the selectivity and potency of SRM inhibitors. Using general, high-throughput assays to quantitatively evaluate inhibitor potency and selectivity in a model sulfate-reducing microbial ecosystem as well as inhibitor specificity for the sulfate reduction pathway in a model SRM, we screened a panel of inorganic oxyanions. We identified several SRM selective inhibitors including selenate, selenite, tellurate, tellurite, nitrate, nitrite, perchlorate, chlorate, monofluorophosphate, vanadate, molydate, and tungstate. Monofluorophosphate (MFP) was not known previously as a selective SRM inhibitor, but has promising characteristics including low toxicity to eukaryotic organisms, high stability at circumneutral pH, utility as an abiotic corrosion inhibitor, and low cost. MFP remains a potent inhibitor of SRM growing by fermentation, and MFP is tolerated by nitrate and perchlorate reducing microorganisms. For SRM inhibition, MFP is synergistic with nitrite and chlorite, and could enhance the efficacy of nitrate or perchlorate treatments. Finally, MFP inhibition is multifaceted. Both inhibition of the central sulfate reduction pathway and release of cytoplasmic fluoride ion are implicated in the mechanism of MFP toxicity.

The adaptation capability of Desulfovibrio to natural fluctuations in electron acceptor availability was evaluated by studying Desulfovibrio alaskensis strain G20 under varying respiratory, fermentative and methanogenic coculture conditions in chemostats. Transition from lactate to pyruvate in coculture resulted in a dramatic shift in the population structure and closer interspecies cell-to-cell interactions. Lower methane production rates in coculture than predicted from pyruvate input was attributed to redirection of electron flow to fumarate reduction. Without a methanogenic partner, accumulation of H₂and formate resulted in greater succinate production. Comparative transcript and gene fitness analysis in concert with physiological data of G20 wildtype and mutants demonstrated that pyruvate fermentation involves respiration of cytoplasmically formed fumarate using cytoplasmic and membrane-bound energy-conserving complexes, Rnf, Hdr-Flox-1 and Hmc. At the low H₂/formate levels maintained in coculture, Rnf likely functions as proton-pumping ferredoxin (Fd): type-I cytochrome c oxidoreductase, which transitions to a proton-pumping Fd(red):  nicotinamide adenine dinucleotide (NAD⁺) oxidoreductase at high H₂/formate levels during fermentation in monoculture. Hdr-Flox-1 is postulated to recycle Fd(red) via a flavin-based electron bifurcation involving NADH, Fdox and the thiol/disulphide-containing DsrC. In a menaquinone (MQ)-based electron confurcation reaction, the high-molecular-weight cytochrome-c₃complex, Hmc, is proposed to then couple DsrC(red) and periplasmic H₂/formate oxidation using the MQ pool to fuel a membrane-bound fumarate reductase.

UnlabelledThe genomes of sulfate-reducing bacteria remain poorly characterized, largely due to a paucity of experimental data and genetic tools. To meet this challenge, we generated an archived library of 15,477 mapped transposon insertion mutants in the sulfate-reducing bacterium Desulfovibrio alaskensis G20. To demonstrate the utility of the individual mutants, we profiled gene expression in mutants of six regulatory genes and used these data, together with 1,313 high-confidence transcription start sites identified by tiling microarrays and transcriptome sequencing (5' RNA-Seq), to update the regulons of Fur and Rex and to confirm the predicted regulons of LysX, PhnF, PerR, and Dde_3000, a histidine kinase. In addition to enabling single mutant investigations, the D. alaskensis G20 transposon mutants also contain DNA bar codes, which enables the pooling and analysis of mutant fitness for thousands of strains simultaneously. Using two pools of mutants that represent insertions in 2,369 unique protein-coding genes, we demonstrate that the hypothetical gene Dde_3007 is required for methionine biosynthesis. Using comparative genomics, we propose that Dde_3007 performs a missing step in methionine biosynthesis by transferring a sulfur group to O-phosphohomoserine to form homocysteine. Additionally, we show that the entire choline utilization cluster is important for fitness in choline sulfate medium, which confirms that a functional microcompartment is required for choline oxidation. Finally, we demonstrate that Dde_3291, a MerR-like transcription factor, is a choline-dependent activator of the choline utilization cluster. Taken together, our data set and genetic resources provide a foundation for systems-level investigation of a poorly studied group of bacteria of environmental and industrial importance.ImportanceSulfate-reducing bacteria contribute to global nutrient cycles and are a nuisance for the petroleum industry. Despite their environmental and industrial significance, the genomes of sulfate-reducing bacteria remain poorly characterized. Here, we describe a genetic approach to fill gaps in our knowledge of sulfate-reducing bacteria. We generated a large collection of archived, transposon mutants in Desulfovibrio alaskensis G20 and used the phenotypes of these mutant strains to infer the function of genes involved in gene regulation, methionine biosynthesis, and choline utilization. Our findings and mutant resources will enable systematic investigations into gene function, energy generation, stress response, and metabolism for this important group of bacteria.

Sulfate-reducing bacteria play major roles in the global carbon and sulfur cycles, but it remains unclear how reducing sulfate yields energy. To determine the genetic basis of energy conservation, we measured the fitness of thousands of pooled mutants of Desulfovibrio alaskensis G20 during growth in 12 different combinations of electron donors and acceptors. We show that ion pumping by the ferredoxin:NADH oxidoreductase Rnf is required whenever substrate-level phosphorylation is not possible. The uncharacterized complex Hdr/flox-1 (Dde_1207:13) is sometimes important alongside Rnf and may perform an electron bifurcation to generate more reduced ferredoxin from NADH to allow further ion pumping. Similarly, during the oxidation of malate or fumarate, the electron-bifurcating transhydrogenase NfnAB-2 (Dde_1250:1) is important and may generate reduced ferredoxin to allow additional ion pumping by Rnf. During formate oxidation, the periplasmic [NiFeSe] hydrogenase HysAB is required, which suggests that hydrogen forms in the periplasm, diffuses to the cytoplasm, and is used to reduce ferredoxin, thus providing a substrate for Rnf. During hydrogen utilization, the transmembrane electron transport complex Tmc is important and may move electrons from the periplasm into the cytoplasmic sulfite reduction pathway. Finally, mutants of many other putative electron carriers have no clear phenotype, which suggests that they are not important under our growth conditions, although we cannot rule out genetic redundancy.

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.

The mineralization of organic matter in anoxic environments relies on the cooperative activities of hydrogen producers and consumers obligately linked by interspecies metabolite exchange in syntrophic consortia that may include sulfate reducing species such as Desulfovibrio. To evaluate the metabolic flexibility of syntrophic Desulfovibrio to adapt to naturally fluctuating methanogenic environments, we studied Desulfovibrio alaskensis strain G20 grown in chemostats under respiratory and syntrophic conditions with alternative methanogenic partners, Methanococcus maripaludis and Methanospirillum hungatei, at different growth rates. Comparative whole-genome transcriptional analyses, complemented by G20 mutant strain growth experiments and physiological data, revealed a significant influence of both energy source availability (as controlled by dilution rate) and methanogen on the electron transfer systems, ratios of interspecies electron carriers, energy generating systems, and interspecies physical associations. A total of 68 genes were commonly differentially expressed under syntrophic versus respiratory lifestyle. Under low-energy (low-growth-rate) conditions, strain G20 further had the capacity to adapt to the metabolism of its methanogenic partners, as shown by its differing gene expression of enzymes involved in the direct metabolic interactions (e.g., periplasmic hydrogenases) and the ratio shift in electron carriers used for interspecies metabolite exchange (hydrogen/formate). A putative monomeric [Fe-Fe] hydrogenase and Hmc (high-molecular-weight-cytochrome c3) complex-linked reverse menaquinone (MQ) redox loop become increasingly important for the reoxidation of the lactate-/pyruvate oxidation-derived redox pair, DsrC(red) and Fd(red), relative to the Qmo-MQ-Qrc (quinone-interacting membrane-bound oxidoreductase; quinone-reducing complex) loop. Together, these data underscore the high enzymatic and metabolic adaptive flexibility that likely sustains Desulfovibrio in naturally fluctuating methanogenic environments.

Gene regulation in bacteria is usually described as an adaptive response to an environmental change so that genes are expressed when they are required. We instead propose that most genes are under indirect control: their expression responds to signal(s) that are not directly related to the genes' function. Indirect control should perform poorly in artificial conditions, and we show that gene regulation is often maladaptive in the laboratory. In Shewanella oneidensis MR-1, 24% of genes are detrimental to fitness in some conditions, and detrimental genes tend to be highly expressed instead of being repressed when not needed. In diverse bacteria, there is little correlation between when genes are important for optimal growth or fitness and when those genes are upregulated. Two common types of indirect control are constitutive expression and regulation by growth rate; these occur for genes with diverse functions and often seem to be suboptimal. Because genes that have closely related functions can have dissimilar expression patterns, regulation may be suboptimal in the wild as well as in the laboratory.

Abstract Background Despite considerable progress in our understanding of land plant phylogeny, several nodes in the green tree of life remain poorly resolved. Furthermore, the bulk of currently available data come from only a subset of major land plant clades. Here we examine early land plant evolution using complete plastome sequences including two previously unexamined and phylogenetically critical lineages. To better understand the evolution of land plants and their plastomes, we examined aligned nucleotide sequences, indels, gene and nucleotide composition, inversions, and gene order at the boundaries of the inverted repeats. Results We present the plastome sequences of Equisetum arvense, a horsetail, and of Isoetes flaccida, a heterosporous lycophyte. Phylogenetic analysis of aligned nucleotides from 49 plastome genes from 43 taxa supported monophyly for the following clades: embryophytes (land plants), lycophytes, monilophytes (leptosporangiate ferns + Angiopteris evecta + Psilotum nudum + Equisetum arvense), and seed plants. Resolution among the four monilophyte lineages remained moderate, although nucleotide analyses suggested that P. nudum and E. arvense form a clade sister to A. evecta + leptosporangiate ferns. Results from phylogenetic analyses of nucleotides were consistent with the distribution of plastome gene rearrangements and with analysis of sequence gaps resulting from insertions and deletions (indels). We found one new indel and an inversion of a block of genes that unites the monilophytes. Conclusions Monophyly of monilophytes has been disputed on the basis of morphological and fossil evidence. In the context of a broad sampling of land plant data we find several new pieces of evidence for monilophyte monophyly. Results from this study demonstrate resolution among the four monilophytes lineages, albeit with moderate support; we posit a clade consisting of Equisetaceae and Psilotaceae that is sister to the "true ferns," including Marattiaceae.

Abstract Background Tortula ruralis, a widely distributed species in the moss family Pottiaceae, is increasingly used as a model organism for the study of desiccation tolerance and mechanisms of cellular repair. In this paper, we present the chloroplast genome sequence of T. ruralis, only the second published chloroplast genome for a moss, and the first for a vegetatively desiccation-tolerant plant. Results The Tortula chloroplast genome is ~123,500 bp, and differs in a number of ways from that of Physcomitrella patens, the first published moss chloroplast genome. For example, Tortula lacks the ~71 kb inversion found in the large single copy region of the Physcomitrella genome and other members of the Funariales. Also, the Tortula chloroplast genome lacks petN, a gene found in all known land plant plastid genomes. In addition, an unusual case of nucleotide polymorphism was discovered. Conclusions Although the chloroplast genome of Tortula ruralis differs from that of the only other sequenced moss, Physcomitrella patens, we have yet to determine the biological significance of the differences. The polymorphisms we have uncovered in the sequencing of the genome offer a rare possibility (for mosses) of the generation of DNA markers for fine-level phylogenetic studies, or to investigate individual variation within populations.

While there has been strong support for Amborella and Nymphaeales (water lilies) as branching from basal-most nodes in the angiosperm phylogeny, this hypothesis has recently been challenged by phylogenetic analyses of 61 protein-coding genes extracted from the chloroplast genome sequences of Amborella, Nymphaea and 12 other available land plant chloroplast genomes. These character-rich analyses placed the monocots, represented by three grasses (Poaceae), as sister to all other extant angiosperm lineages. We have extracted protein-coding regions from draft sequences for six additional chloroplast genomes to test whether this surprising result could be an artifact of long-branch attraction due to limited taxon sampling. The added taxa include three monocots (Acorus, Yucca and Typha), a water lily (Nuphar), a ranunculid (Ranunculus), and a gymnosperm (Ginkgo). Phylogenetic analyses of the expanded DNA and protein datasets together with microstructural characters (indels) provided unambiguous support for Amborella and the Nymphaeales as branching from the basal-most nodes in the angiosperm phylogeny. However, their relative positions proved to be dependent on method of analysis, with parsimony favoring Amborella as sister to all other angiosperms, and maximum likelihood and neighbor-joining methods favoring an Amborella + Nympheales clade as sister. The maximum likelihood phylogeny supported the later hypothesis, but the likelihood for the former hypothesis was not significantly different. Parametric bootstrap analysis, single gene phylogenies, estimated divergence dates and conflicting indel characters all help to illuminate the nature of the conflict in resolution of the most basal nodes in the angiosperm phylogeny. Molecular dating analyses provided median age estimates of 161mya for the most recent common ancestor of all extant angiosperms and 145 mya for the most recent common ancestor of monocots, magnoliids and eudicots. Whereas long sequences reduce variance in branch lengths and molecular dating estimates, the impact of improved taxon sampling on the rooting of the angiosperm phylogeny together with the results of parametric bootstrap analyses demonstrate how long-branch attraction can mislead genome-scale phylogenetic analyses.

We used a unique combination of techniques to sequence the first complete chloroplast genome of a lycophyte, Huperzia lucidula. This plant belongs to a significant clade hypothesized to represent the sister group to all other vascular plants. We used fluorescence-activated cell sorting (FACS) to isolate the organelles, rolling circle amplification (RCA) to amplify the genome, and shotgun sequencing to 8x depth coverage to obtain the complete chloroplast genome sequence. The genome is 154,373 bp, containing inverted repeats of 15,314 bp each, a large single-copy region of 104,088 bp, and a small single-copy region of 19,671 bp. Gene order is more similar to those of mosses, liverworts, and hornworts than to gene order for other vascular plants. For example, the Huperzia chloroplast genome possesses the bryophyte gene order for a previously characterized 30 kb inversion, thus supporting the hypothesis that lycophytes are sister to all other extant vascular plants. The lycophyte chloroplast genome data also enable a better reconstruction of the basal tracheophyte genome, which is useful for inferring relationships among bryophyte lineages. Several unique characters are observed in Huperzia, such as movement of the gene ndhF from the small single copy region into the inverted repeat. We present several analyses of evolutionary relationships among land plants by using nucleotide data, amino acid sequences, and by comparing gene arrangements from chloroplast genomes. The results, while still tentative pending the large number of chloroplast genomes from other key lineages that are soon to be sequenced, are intriguing in themselves, and contribute to a growing comparative database of genomic and morphological data across the green plants.