Vincent J. Denef

Burkholderia Xenovorans LB400 Harbors a Multi-Replicon, 9.73-Mbp Genome Shaped for Versatility

Burkholderia xenovorans LB400 (LB400), a well studied, effective polychlorinated biphenyl-degrader, has one of the two largest known bacterial genomes and is the first nonpathogenic Burkholderia isolate sequenced. From an evolutionary perspective, we find significant differences in functional specialization between the three replicons of LB400, as well as a more relaxed selective pressure for genes located on the two smaller vs. the largest replicon. High genomic plasticity, diversity, and specialization within the Burkholderia genus are exemplified by the conservation of only 44% of the genes between LB400 and Burkholderia cepacia complex strain 383. Even among four B. xenovorans strains, genome size varies from 7.4 to 9.73 Mbp. The latter is largely explained by our findings that >20% of the LB400 sequence was recently acquired by means of lateral gene transfer. Although a range of genetic factors associated with in vivo survival and intercellular interactions are present, these genetic factors are likely related to niche breadth rather than determinants of pathogenicity. The presence of at least eleven “central aromatic” and twenty “peripheral aromatic” pathways in LB400, among the highest in any sequenced bacterial genome, supports this hypothesis. Finally, in addition to the experimentally observed redundancy in benzoate degradation and formaldehyde oxidation pathways, the fact that 17.6% of proteins have a better LB400 paralog than an ortholog in a different genome highlights the importance of gene duplication and repeated acquirement, which, coupled with their divergence, raises questions regarding the role of paralogs and potential functional redundancies in large-genome microbes.

Strain-resolved community proteomics reveals recombining genomes of acidophilic bacteria

Microbes comprise the majority of extant organisms, yet much remains to be learned about the nature and driving forces of microbial diversification. Our understanding of how microorganisms adapt and evolve can be advanced by genome-wide documentation of the patterns of genetic exchange, particularly if analyses target coexisting members of natural communities. Here we use community genomic data sets to identify, with strain specificity, expressed proteins from the dominant member of a genomically uncharacterized, natural, acidophilic biofilm. Proteomics results reveal a genome shaped by recombination involving chromosomal regions of tens to hundreds of kilobases long that are derived from two closely related bacterial populations. Inter-population genetic exchange was confirmed by multilocus sequence typing of isolates and of uncultivated natural consortia. The findings suggest that exchange of large blocks of gene variants is crucial for the adaptation to specific ecological niches within the very acidic, metal-rich environment. Mass-spectrometry-based discrimination of expressed protein products that differ by as little as a single amino acid enables us to distinguish the behaviour of closely related coexisting organisms. This is important, given that microorganisms grouped together as a single species may have quite distinct roles in natural systems and their interactions might be key to ecosystem optimization. Because proteomic data simultaneously convey information about genome type and activity, strain-resolved community proteomics is an important complement to cultivation-independent genomic (metagenomic) analysis of microorganisms in the natural environment.

Systems Biology: Functional analysis of natural microbial consortia using community proteomics

We know very little about the metabolic functioning and evolutionary dynamics of microbial communities. Recent advances in comprehensive, sequencing-based methods, however, are laying a molecular foundation for new insights into how microbial communities shape the Earths biosphere. Here we explore the convergence of microbial ecology, genomics, biological mass spectrometry and informatics that form the new field of microbial community proteogenomics. We discuss the first applications of proteogenomics and its potential for studying the physiology, ecology and evolution of microbial populations and communities.

Functional analysis of natural microbial consortia using community proteomics

We know very little about the metabolic functioning and evolutionary dynamics of microbial communities. Recent advances in comprehensive, sequencing-based methods, however, are laying a molecular foundation for new insights into how microbial communities shape the Earths biosphere. Here we explore the convergence of microbial ecology, genomics, biological mass spectrometry and informatics that form the new field of microbial community proteogenomics. We discuss the first applications of proteogenomics and its potential for studying the physiology, ecology and evolution of microbial populations and communities.

Community Genomic and Proteomic Analyses of Chemoautotrophic Iron-Oxidizing “Leptospirillum rubarum” (Group II) and “Leptospirillum ferrodiazotrophum” (Group III) Bacteria in Acid Mine Drainage Biofilms

ABSTRACT We analyzed near-complete population (composite) genomic sequences for coexisting acidophilic iron-oxidizing Leptospirillum group II and III bacteria (phylum Nitrospirae) and an extrachromosomal plasmid from a Richmond Mine, Iron Mountain, CA, acid mine drainage biofilm. Community proteomic analysis of the genomically characterized sample and two other biofilms identified 64.6% and 44.9% of the predicted proteins of Leptospirillum groups II and III, respectively, and 20% of the predicted plasmid proteins. The bacteria share 92% 16S rRNA gene sequence identity and >60% of their genes, including integrated plasmid-like regions. The extrachromosomal plasmid carries conjugation genes with detectable sequence similarity to genes in the integrated conjugative plasmid, but only those on the extrachromosomal element were identified by proteomics. Both bacterial groups have genes for community-essential functions, including carbon fixation and biosynthesis of vitamins, fatty acids, and biopolymers (including cellulose); proteomic analyses reveal these activities. Both Leptospirillum types have multiple pathways for osmotic protection. Although both are motile, signal transduction and methyl-accepting chemotaxis proteins are more abundant in Leptospirillum group III, consistent with its distribution in gradients within biofilms. Interestingly, Leptospirillum group II uses a methyl-dependent and Leptospirillum group III a methyl-independent response pathway. Although only Leptospirillum group III can fix nitrogen, these proteins were not identified by proteomics. The abundances of core proteins are similar in all communities, but the abundance levels of unique and shared proteins of unknown function vary. Some proteins unique to one organism were highly expressed and may be key to the functional and ecological differentiation of Leptospirillum groups II and III.

Population Genomic Analysis of Strain Variation in Leptospirillum Group II Bacteria Involved in Acid Mine Drainage Formation

Deeply sampled community genomic (metagenomic) datasets enable comprehensive analysis of heterogeneity in natural microbial populations. In this study, we used sequence data obtained from the dominant member of a low-diversity natural chemoautotrophic microbial community to determine how coexisting closely related individuals differ from each other in terms of gene sequence and gene content, and to uncover evidence of evolutionary processes that occur over short timescales. DNA sequence obtained from an acid mine drainage biofilm was reconstructed, taking into account the effects of strain variation, to generate a nearly complete genome tiling path for a Leptospirillum group II species closely related to L. ferriphilum (sampling depth ∼20×). The population is dominated by one sequence type, yet we detected evidence for relatively abundant variants (>99.5% sequence identity to the dominant type) at multiple loci, and a few rare variants. Blocks of other Leptospirillum group II types (∼94% sequence identity) have recombined into one or more variants. Variant blocks of both types are more numerous near the origin of replication. Heterogeneity in genetic potential within the population arises from localized variation in gene content, typically focused in integrated plasmid/phage-like regions. Some laterally transferred gene blocks encode physiologically important genes, including quorum-sensing genes of the LuxIR system. Overall, results suggest inter- and intrapopulation genetic exchange involving distinct parental genome types and implicate gain and loss of phage and plasmid genes in recent evolution of this Leptospirillum group II population. Population genetic analyses of single nucleotide polymorphisms indicate variation between closely related strains is not maintained by positive selection, suggesting that these regions do not represent adaptive differences between strains. Thus, the most likely explanation for the observed patterns of polymorphism is divergence of ancestral strains due to geographic isolation, followed by mixing and subsequent recombination.

Persisting Viral Sequences Shape Microbial CRISPR-based Immunity

Well-studied innate immune systems exist throughout bacteria and archaea, but a more recently discovered genomic locus may offer prokaryotes surprising immunological adaptability. Mediated by a cassette-like genomic locus termed Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR), the microbial adaptive immune system differs from its eukaryotic immune analogues by incorporating new immunities unidirectionally. CRISPR thus stores genomically recoverable timelines of virus-host coevolution in natural organisms refractory to laboratory cultivation. Here we combined a population genetic mathematical model of CRISPR-virus coevolution with six years of metagenomic sequencing to link the recoverable genomic dynamics of CRISPR loci to the unknown population dynamics of virus and host in natural communities. Metagenomic reconstructions in an acid-mine drainage system document CRISPR loci conserving ancestral immune elements to the base-pair across thousands of microbial generations. This ‘trailer-end conservation’ occurs despite rapid viral mutation and despite rapid prokaryotic genomic deletion. The trailer-ends of many reconstructed CRISPR loci are also largely identical across a population. ‘Trailer-end clonality’ occurs despite predictions of host immunological diversity due to negative frequency dependent selection (kill the winner dynamics). Statistical clustering and model simulations explain this lack of diversity by capturing rapid selective sweeps by highly immune CRISPR lineages. Potentially explaining ‘trailer-end conservation,’ we record the first example of a viral bloom overwhelming a CRISPR system. The polyclonal viruses bloom even though they share sequences previously targeted by host CRISPR loci. Simulations show how increasing random genomic deletions in CRISPR loci purges immunological controls on long-lived viral sequences, allowing polyclonal viruses to bloom and depressing host fitness. Our results thus link documented patterns of genomic conservation in CRISPR loci to an evolutionary advantage against persistent viruses. By maintaining old immunities, selection may be tuning CRISPR-mediated immunity against viruses reemerging from lysogeny or migration.

The dynamic genetic repertoire of microbial communities

Community genomic data have revealed multiple levels of variation between and within microbial consortia. This variation includes large-scale differences in gene content between ecosystems as well as within-population sequence heterogeneity. In the present review, we focus specifically on how fine-scale variation within microbial and viral populations is apparent from community genomic data. A major unresolved question is how much of the observed variation is due to neutral vs. adaptive processes. Limited experimental data hint that some of this fine-scale variation may be in part functionally relevant, whereas sequence-based and modeling analyses suggest that much of it may be neutral. While methods for interpreting population genomic data are still in their infancy, we discuss current interpretations of existing datasets in the light of evolutionary processes and models. Finally, we highlight the importance of virus–host dynamics in generating and shaping within-population diversity.

Ecological distribution and population physiology defined by proteomics in a natural microbial community

An important challenge in microbial ecology is developing methods that simultaneously examine the physiology of organisms at the molecular level and their ecosystem level interactions in complex natural systems. We integrated extensive proteomic, geochemical, and biological information from 28 microbial communities collected from an acid mine drainage environment and representing a range of biofilm development stages and geochemical conditions to evaluate how the physiologies of the dominant and less abundant organisms change along environmental gradients. The initial colonist dominates across all environments, but its proteome changes between two stable states as communities diversify, implying that interspecies interactions affect this organisms metabolism. Its overall physiology is robust to abiotic environmental factors, but strong correlations exist between these factors and certain subsets of proteins, possibly accounting for its wide environmental distribution. Lower abundance populations are patchier in their distribution, and proteomic data indicate that their environmental niches may be constrained by specific sets of abiotic environmental factors. This research establishes an effective strategy to investigate ecological relationships between microbial physiology and the environment for whole communities in situ.

Coping with Polychlorinated Biphenyl (PCB) Toxicity: Physiological and Genome-Wide Responses of Burkholderia xenovorans LB400 to PCB-Mediated Stress

ABSTRACT The biodegradation of polychlorinated biphenyls (PCBs) relies on the ability of aerobic microorganisms such as Burkholderia xenovorans sp. LB400 to tolerate two potential modes of toxicity presented by PCB degradation: passive toxicity, as hydrophobic PCBs potentially disrupt membrane and protein function, and degradation-dependent toxicity from intermediates of incomplete degradation. We monitored the physiological characteristics and genome-wide expression patterns of LB400 in response to the presence of Aroclor 1242 (500 ppm) under low expression of the structural biphenyl pathway (succinate and benzoate growth) and under induction by biphenyl. We found no inhibition of growth or change in fatty acid profile due to PCBs under nondegrading conditions. Moreover, we observed no differential gene expression due to PCBs themselves. However, PCBs did have a slight effect on the biosurface area of LB400 cells and caused slight membrane separation. Upon activation of the biphenyl pathway, we found growth inhibition from PCBs beginning after exponential-phase growth suggestive of the accumulation of toxic compounds. Genome-wide expression profiling revealed 47 differentially expressed genes (0.56% of all genes) under these conditions. The biphenyl and catechol pathways were induced as expected, but the quinoprotein methanol metabolic pathway and a putative chloroacetaldehyde dehydrogenase were also highly expressed. As the latter protein is essential to conversion of toxic metabolites in dichloroethane degradation, it may play a similar role in the degradation of chlorinated aliphatic compounds resulting from PCB degradation.