Exposing the third chromosome of Burkholderia cepacia complex strains as a virulence plasmid

Transposon sequencing reveals the essential gene set and genes enabling gut symbiosis in the insect symbiont Caballeronia insecticola

Caballeronia insecticola is a bacterium belonging to the Burkholderia genus sensu lato, which is able to colonize multiple environments like soils and the gut of the bean bug Riptortus pedestris. We constructed a saturated Himar1 mariner transposon library and revealed by transposon-sequencing that 498 protein-coding genes constitute the essential genome of Caballeronia insecticola for growth in free-living conditions. By comparing essential gene sets of Caballeronia insecticola and seven related Burkholderia s.l. strains, only 120 common genes were identified, indicating that a large part of the essential genome is strain-specific. In order to reproduce specific nutritional conditions that are present in the gut of Riptortus pedestris, we grew the mutant library in minimal media supplemented with candidate gut nutrients and identified several condition-dependent fitness-defect genes by transposon-sequencing. To validate the robustness of the approach, insertion mutants in six fitness genes were constructed and their growth deficiency in media supplemented with the corresponding nutrient was confirmed. The mutants were further tested for their efficiency in Riptortus pedestris gut colonization, confirming that gluconeogenic carbon sources, taurine and inositol, are nutrients consumed by the symbiont in the gut. Thus, our study provides insights about specific contributions provided by the insect host to the bacterial symbiont.

Drosophila melanogaster as an organism model for studying cystic fibrosis and its major associated microbial infections

Cystic fibrosis (CF) is a human genetic disease caused by mutations in the cystic fibrosis transmembrane conductance regulator gene that encodes a chloride channel. The most severe clinical manifestation is associated with chronic pulmonary infections by pathogenic and opportunistic microbes. Drosophila melanogaster has become the invertebrate model of choice for modeling microbial infections and studying the induced innate immune response. Here, we review its contribution to the understanding of infections with six major pathogens associated with CF (Staphylococcus aureus, Pseudomonas aeruginosa, Burkholderia cepacia, Mycobacterium abscessus, Streptococcus pneumoniae, and Aspergillus fumigatus) together with the perspectives opened by the recent availability of two CF models in this model organism.

Independent origins and evolution of the secondary replicons of the class Gammaproteobacteria

Multipartite genomes, consisting of more than one replicon, have been found in approximately 10 % of bacteria, many of which belong to the phylum Proteobacteria. Many aspects of their origin and evolution, and the possible advantages related to this type of genome structure, remain to be elucidated. Here, we performed a systematic analysis of the presence and distribution of multipartite genomes in the class Gammaproteobacteria, which includes several genera with diverse lifestyles. Within this class, multipartite genomes are mainly found in the order Alteromonadales (mostly in the genus Pseudoalteromonas) and in the family Vibrionaceae. Our data suggest that the emergence of secondary replicons in Gammaproteobacteria is rare and that they derive from plasmids. Despite their multiple origins, we highlighted the presence of evolutionary trends such as the inverse proportionality of the genome to chromosome size ratio, which appears to be a general feature of bacteria with multipartite genomes irrespective of taxonomic group. We also highlighted some functional trends. The core gene set of the secondary replicons is extremely small, probably limited to essential genes or genes that favour their maintenance in the genome, while the other genes are less conserved. This hypothesis agrees with the idea that the primary advantage of secondary replicons could be to facilitate gene acquisition through horizontal gene transfer, resulting in replicons enriched in genes associated with adaptation to different ecological niches. Indeed, secondary replicons are enriched both in genes that could promote adaptation to harsh environments, such as those involved in antibiotic, biocide and metal resistance, and in functional categories related to the exploitation of environmental resources (e.g. carbohydrates), which can complement chromosomal functions.

Analysis of multipartite bacterial genomes using alignment free and alignment-based pipelines

Since the discovery of second chromosome in Rhodobacter sphaeroides 2.4.1 in 1989, multipartite genomes have been reported in over three hundred bacterial species under nine different phyla. This has shattered the unipartite (single chromosome) genome dogma in bacteria. Since then, many questions on various aspects of multipartite genomes in bacteria have been addressed. However, our understanding of how multipartite genomes emerge and evolve is still lacking. Importantly, the knowledge of genetic factors underlying the differences in multipartite and single-chromosome genomes is lacking. In this work, we have performed comparative evolutionary and functional genomics analyses to identify molecular factors that discriminate multipartite from unipartite bacteria, with the goal to decipher taxon-specific factors, and those that are prevalent across the taxa, underlying these traits. We assessed the roles of evolutionary mechanisms, specifically gene gain, in driving the divergence of bacteria with single and multiple chromosomes. In addition, we performed functional genomic analysis to garner support for our findings from comparative evolutionary analysis. We found genes such as those encoding conserved hypothetical proteins in Deinococcus radiodurans R1, and putative phage phi-C31 gp36 major capsid like and hypothetical proteins in Rhodobacter sphaeroides 2.4.1, which are located on accessory chromosomes in these bacteria but were not found in the inferred ancestral sequences, and on the primary chromosomes, as well as were not found in their closest relatives with single chromosome within the same clade. Our study shines a new light on the potential roles of the secondary chromosomes in helping bacteria with multipartite genomes to adapt to specialized environments or growth conditions.

Disruption of SMC-related genes promotes recombinant cholesterol esterase production in Burkholderia stabilis

Burkholderia stabilis strain FERMP-21014 secretes cholesterol esterase (BsChe), which is used in clinical settings to determine serum cholesterol levels. Previously, we constructed an expression plasmid with an endogenous constitutive promoter to enable the production of recombinant BsChe. In this study, we obtained one mutant strain with 13.1-fold higher BsChe activity than the wild type, using N-methyl-N'-nitro-N-nitrosoguanidine as a mutagen. DNA-sequencing analysis revealed that the strain had lost chromosome 3 (∆Chr3), suggesting that the genes hindering BsChe production may be encoded on Chr3. We also identified common mutations in the functionally unknown BSFP_068720/30 genes in the top 10 active strains generated during transposon mutagenesis. As BSFP_068720/30/40 comprised an operon on Chr3, we created the BSFP_068720/30/40 disruption mutant and confirmed that each disruption mutant containing the expression plasmid exhibited ~ 16.1-fold higher BsChe activity than the wild type. Quantitative PCR showed that each disruption mutant and ΔChr3 had a ~ 9.4-fold higher plasmid copy number than the wild type. Structural prediction models indicate that BSFP_068730/40 is structurally homologous to the structural maintenance of chromosomes (SMC) protein MukBE, which is responsible for chromosome segregation during cell division. Conversely, BSFP_068720/30/40 disruption did not lead to a Chr3 drop-out. These results imply that BSFP_068720/30/40 is not a SMC protein but is involved in destabilizing foreign plasmids to prevent the influx of genetic information from the environment. In conclusion, the disruption of BSFP_068720/30/40 improved plasmid stability and copy number, resulting in exceptionally high BsChe production. KEY POINTS: • Disruption of BSFP_068720/30/40 enabled mass production of Burkholderia Che/Lip. • BSFP_068730/40 is an SMC protein homolog not involved in chromosome retention. • BSFP_068720/30/40 is likely responsible for the exclusion of exogenous plasmids.

In silico discovery of biomarkers for the accurate and sensitive detection of Fusarium solani

Fusarium solani is worrisome because it severely threatens the agricultural productivity of certain crops such as tomatoes and peas, causing the general decline, wilting, and root necrosis. It has also been implicated in the infection of the human eye cornea. It is believed that early detection of the fungus could save these crops from the destructive activities of the fungus through early biocontrol measures. Therefore, the present work aimed to build a sensitive model of novel anti-Fusarium solani antimicrobial peptides (AMPs) against the fungal cutinase 1 (CUT1) protein for early, sensitive and accurate detection. Fusarium solani CUT1 receptor protein 2D secondary structure, model validation, and functional motifs were predicted. Subsequently, anti-Fusarium solani AMPs were retrieved, and the HMMER in silico algorithm was used to construct a model of the AMPs. After their structure predictions, the interaction analysis was analyzed for the Fusarium solani CUT1 protein and the generated AMPs. The putative anti-Fusarium solani AMPs bound the CUT1 protein very tightly, with OOB4 having the highest binding energy potential for HDock. The pyDockWeb generated high electrostatic, desolvation, and low van der Waals energies for all the AMPs against CUT1 protein, with OOB1 having the most significant interaction. The results suggested the utilization of AMPs for the timely intervention, control, and management of these crops, as mentioned earlier, to improve their agricultural productivity and reduce their economic loss and the use of HMMER for constructing models for disease detection.

The Sixth Element: a 102-kb RepABC Plasmid of Xenologous Origin Modulates Chromosomal Gene Expression in Dinoroseobacter shibae

The model organism Dinoroseobacter shibae and many other marine Rhodobacterales (Roseobacteraceae, Alphaproteobacteria) are characterized by a multipartite genome organization. Here, we show that the original isolate (Dshi-6) contained six extrachromosomal replicons (ECRs), whereas the strain deposited at the DSMZ (Dshi-5) lacked a 102-kb plasmid. To determine the role of the sixth plasmid, we investigated the genomic and physiological differences between the two strains. Therefore, both genomes were (re)sequenced, and gene expression, growth, and substrate utilization were examined. For comparison, we included additional plasmid-cured strains in the analysis. In the Dshi-6 population, the conjugative 102-kb RepABC-9 plasmid was present in only about 50% of the cells, irrespective of its experimentally validated stability. In the presence of the sixth plasmid, copy number changes of other ECRs, in particular, a decrease of the 86-kb plasmid, were observed. The most conspicuous finding was the strong influence of plasmids on chromosomal gene expression, especially the repression of the CtrA regulon and the activation of the denitrification gene cluster. Expression is inversely controlled by either the presence of the 102-kb plasmid or the absence of the 86-kb plasmid. We identified regulatory genes on both plasmids, i.e., a sigma 70 factor and a quorum sensing synthase, that might be responsible for these major changes. The tremendous effects that were probably even underestimated challenge the current understanding of the relevance of volatile plasmids not only for the original host but also for new recipients after conjugation.

IMPORTANCE Plasmids are small DNA molecules that replicate independently of the bacterial chromosome. The common view of the role of plasmids is dominated by the accumulation of resistance genes, which is responsible for the antibiotic crisis in health care and livestock breeding. Beyond rapid adaptations to a changing environment, no general relevance for the host cell’s regulome was attributed to these volatile ECRs. The current study shows for the model organism D. shibae that its chromosomal gene expression is strongly influenced by two plasmids. We provide evidence that the gain or loss of plasmids not only results in minor alterations of the genetic repertoire but also can have tremendous effects on bacterial physiology. The central role of some plasmids in the regulatory network of the host could also explain their persistence despite fitness costs, which has been described as the “plasmid paradox.”

Identification of Key Factors for Anoxic Survival of B. cenocepacia H111

Burkholderia cenocepacia is an opportunistic pathogen that can lead to severe infections in patients suffering from cystic fibrosis (CF) and chronic granulomatous disease. Being an obligate aerobe, B. cenocepacia is unable to grow in the absence of oxygen. In this study, we show that the CF isolate B. cenocepacia H111 can survive in the absence of oxygen. Using a transposon sequencing (Tn-seq) approach, we identified 71 fitness determinants involved in anoxic survival, including a Crp-Fnr family transcriptional regulatory gene (anr₂), genes coding for the sensor kinase RoxS and its response regulator RoxR, the sigma factor for flagella biosynthesis (FliA) and subunits of a cytochrome bd oxidase (CydA, CydB and the potentially novel subunit CydP). Individual knockouts of these fitness determinants significantly reduced anoxic survival, and inactivation of both anr copies is shown to be lethal under anoxic conditions. We also show that the two-component system RoxS/RoxR and FliA are important for virulence and swarming/swimming, respectively.

Iron metabolism in bacterial cells: from physiological significance to a new class of antimicrobial agents

Infectious complications in the respiratory tract caused by microorganisms from the Burkholderia cepacia complex are the main cause of death among patients with cystic fibrosis. Natural and acquired resistance mechanisms allow Burkholderia cepacia complex pathogens to adapt to the conditions of regular antibiotic therapy, which necessitates the use of antibacterial drugs with an alternative mechanism of action. Studies on the importance of iron as an essential factor in the metabolism of bacteria and methods of its acquisition from the environment contributed to the development of a new antibiotic from a number of cephalosporins – cefiderocol. In the structure of cefiderocol, a fragment is formed that imitates siderophores – chelating molecules that ensure the transport of iron ions into the internal environment of the microorganism. A unique mechanism, described in the scientific literature as a “Trojan horse”, allows antibiotic molecules conjugated with siderophores to effectively penetrate into the bacterial cell, exerting a bactericidal effect. Thus, cefiderocol can be used to treat infectious complications in the lungs of patients with cystic fibrosis caused by bacteria from the Burkholderia cepacia complex, including multidrug-resistant strains. In addition, the spectrum of activity of cefiderocol allows the use of this antibiotic in the treatment of infections caused by nosocomial gram-negative bacteria such as Enterobacterales, Acinetobacter, Pseudomonas and Stenotrophomonas.

Pangenome inventory of Burkholderia sensu lato, Burkholderia sensu stricto, and the Burkholderia cepacia complex reveals the uniqueness of Burkholderia catarinensis

Here the pangenome analysis of Burkholderia sensu lato (s.l.) was performed for the first time, together with an updated analysis of the pangenome of Burkholderia sensu stricto, and Burkholderia cepacia complex (Bcc) focusing on the Bcc B. catarinensis specific features of its re-sequenced genome. The pangenome of Burkholderia s.l., Burkholderia s.s., and of the Bcc are open, composed of more than 96% of accessory genes, and more than 62% of unknown genes. Functional annotations showed that secondary metabolism genes belong to the variable portion of genomes, which might explain their production of several compounds with varied bioactivities. Taken together, this work shows the great variability and uniqueness of these genomes and reveals an underexplored unknown potential in poorly characterized genes. Regarding B. catarinensis 89^T, its genome harbors genes related to hydrolases production and plant growth promotion. This draft genome will be valuable for further investigation of its biotechnological potentials.

Methodological tools to study species of the genus Burkholderia

Bacteria belonging to the Burkholderia genus are extremely versatile and diverse. They can be environmental isolates, opportunistic pathogens in cystic fibrosis, immunocompromised or chronic granulomatous disease patients, or cause disease in healthy people (e.g., Burkholderia pseudomallei) or animals (as in the case of Burkholderia mallei). Since the genus was separated from the Pseudomonas one in the 1990s, the methodological tools to study and characterize these bacteria are evolving fast. Here we reviewed the techniques used in the last few years to update the taxonomy of the genus, to study gene functions and regulations, to deepen the knowledge on the drug resistance which characterizes these bacteria, and to elucidate their mechanisms to establish infections. The availability of these tools significantly impacts the quality of research on Burkholderia and the choice of the most appropriated is fundamental for a precise characterization of the species of interest.

Key points

• Updated techniques to study the genus Burkholderia were reviewed.

• Taxonomy, genomics, assays, and animal models were described.

• A comprehensive overview on recent advances in Burkholderia studies was made.

Precise measurement of the fitness effects of spontaneous mutations by droplet digital PCR in Burkholderia cenocepacia

Understanding how mutations affect survivability is a key component to knowing how organisms and complex traits evolve. However, most mutations have a minor effect on fitness and these effects are difficult to resolve using traditional molecular techniques. Therefore, there is a dire need for more accurate and precise fitness measurements methods. Here, we measured the fitness effects in Burkholderia cenocepacia HI2424 mutation accumulation (MA) lines using droplet-digital polymerase chain reaction (ddPCR). Overall, the fitness measurements from ddPCR-MA are correlated positively with fitness measurements derived from traditional phenotypic marker assays (r = 0.297, P = 0.05), but showed some differences. First, ddPCR had significantly lower measurement variance in fitness (F = 3.78, P < 2.6 × 10-13) in control experiments. Second, the mean fitness from ddPCR-MA measurements were significantly lower than phenotypic marker assays (-0.0041 vs -0.0071, P = 0.006). Consistent with phenotypic marker assays, ddPCR-MA measurements observed multiple (27/43) lineages that significantly deviated from mean fitness, suggesting that a majority of the mutations are neutral or slightly deleterious and intermixed with a few mutations that have extremely large effects. Of these mutations, we found a significant excess of mutations within DNA excinuclease and Lys R transcriptional regulators that have extreme deleterious and beneficial effects, indicating that modifications to transcription and replication may have a strong effect on organismal fitness. This study demonstrates the power of ddPCR as a ubiquitous method for high-throughput fitness measurements in both DNA- and RNA-based organisms regardless of cell type or physiology.

Conservation of Resistance-Nodulation-Cell Division Efflux Pump-Mediated Antibiotic Resistance in Burkholderia cepacia Complex and Burkholderia pseudomallei Complex Species

Burkholderia cepacia complex (Bcc) and Burkholderia pseudomallei complex (Bpc) species include pathogens that are typically multidrug resistant. Dominant intrinsic and acquired multidrug resistance mechanisms are efflux mediated by pumps of the resistance-nodulation-cell division (RND) family. From comparative bioinformatic and, in many instances, functional studies, we infer that RND pump-based resistance mechanisms are conserved in Burkholderia. We propose to use these findings as a foundation for adoption of a uniform RND efflux pump nomenclature.

Burkholderia in the genomic era: from taxonomy to the discovery of new antimicrobial secondary metabolites

Species of Burkholderia are highly versatile being found not only abundantly in soil, but also as plants and animals' commensals or pathogens. Their complex multireplicon genomes harbour an impressive number of polyketide synthase (PKS) and nonribosomal peptide-synthetase (NRPS) genes coding for the production of antimicrobial secondary metabolites (SMs), which have been successfully deciphered by genome-guided tools. Moreover, genome metrics supported the split of this genus into Burkholderia sensu stricto (s.s.) and five new other genera. Here, we show that the successful antimicrobial SMs producers belong to Burkholderia s.s. Additionally, we reviewed the occurrence, bioactivities, modes of action, structural, and biosynthetic information of thirty-eight Burkholderia antimicrobial SMs shedding light on their diversity, complexity, and uniqueness as well as the importance of genome-guided strategies to facilitate their discovery. Several Burkholderia NRPS and PKS display unusual features, which are reflected in their structural diversity, important bioactivities, and varied modes of action. Up to now, it is possible to observe a general tendency of Burkholderia SMs being more active against fungi. Although the modes of action and biosynthetic gene clusters of many SMs remain unknown, we highlight the potential of Burkholderia SMs as alternatives to fight against new diseases and antibiotic resistance.

Presence of the Hmq System and Production of 4-Hydroxy-3-Methyl-2-Alkylquinolines Are Heterogeneously Distributed between Burkholderia cepacia Complex Species and More Prevalent among Environmental than Clinical Isolates

The Burkholderia cepacia complex (Bcc) comprises several species of closely related, versatile bacteria. Some Bcc strains produce 4-hydroxy-3-methyl-2-alkylquinolines (HMAQs), analogous to the 4-hydroxy-2-alkylquinolines of Pseudomonas aeruginosa. Using in silico analyses, we previously estimated that the hmqABCDEFG operon, which encodes enzymes involved in the biosynthesis of HMAQs, is carried by about one-third of Bcc strains, with considerable inter- and intraspecies variability. In the present study, we investigated by PCR, using consensus primers, the distribution of hmqABCDEFG in a collection of 312 Bcc strains (222 of clinical and 90 of environmental origins) belonging to 18 Bcc species. We confirmed that this operon is not distributed evenly among Bcc species. Among the 30% of strains bearing the hmqABCDEFG operon, we found that 92% of environmental isolates and 82% of clinically isolated Bcc strains produce levels of HMAQs detectable by liquid chromatography-mass spectrometry in at least one of the tested culture conditions. Among the hmqABCDEFG-positive but HMAQ-negative strains, none expressed the hmqA gene under the specified culture conditions. Interestingly, the hmqABCDEFG operon is more prevalent among plant root environment species (e.g., Burkholderia ambifaria and Burkholderia cepacia) and absent in species commonly found in chronically colonized individuals with cystic fibrosis (e.g., Burkholderia cenocepacia and Burkholderia multivorans), suggesting a role for the Hmq system in niche adaptation. We investigated the impact of the Hmq system on plant growth promotion and found that Pisum sativum root development by B. ambifaria required a functional HMAQ system.

IMPORTANCE Environmental bacteria belonging to the various closely related species forming the Burkholderia cepacia complex (Bcc) can infect plants and animals, including humans. Their pathogenicity is regulated by intercellular communication, or quorum sensing, allowing them to collaborate instead of acting individually. Bcc organisms generally exploit interacting quorum sensing systems based on N-acyl-homoserine lactones as signaling molecules. Several Bcc strains also carry an hmqABCDEFG operon responsible for the biosynthesis of 4-hydroxy-3-methyl-2-alkylquinolines (HMAQs), molecules analogous to the Pseudomonas quinolone signal (PQS) system of P. aeruginosa. Our finding that the prevalences of the Hmq system and HMAQ production are very different between various Bcc species suggests a key role in niche adaptation or pathogenicity. This is supported by a significant reduction in plant growth promotion in the absence of HMAQ production for a beneficial Bcc strain.

An investigation of Burkholderia cepacia complex methylomes via SMRT sequencing and mutant analysis

Bacterial genomes can be methylated at particular motifs by methyltransferases (M). This DNA modification allows restriction endonucleases (R) to discriminate between self and foreign DNA. While the accepted primary function of such restriction modification (RM) systems is to degrade incoming foreign DNA, other roles of RM systems and lone R or M components have been found in genome protection, stability and the regulation of various phenotypes. The Burkholderia cepacia complex (Bcc) is a group of closely related opportunistic pathogens with biotechnological potential. Here, we constructed and analysed mutants lacking various RM components in the clinical Bcc isolate Burkholderia cenocepacia H111 and used SMRT sequencing of single mutants to assign the B. cenocepacia H111 Ms to their cognate motifs. DNA methylation is shown to affect biofilm formation, cell shape, motility, siderophore production and membrane vesicle production. Moreover, DNA methylation had a large effect on the maintenance of the Bcc virulence megaplasmid pC3. Our data also suggest that the gp51 M-encoding gene, which is essential in H111 and is located within a prophage, is required for maintaining the bacteriophage in a lysogenic state, thereby ensuring a constant, low level of phage production within the bacterial population.ImportanceWhile genome sequence determines an organism's proteins, methylation of the nucleotides themselves can confer additional properties. In bacteria, Ms modify specific nucleotide motifs to allow discrimination of 'self' from 'non-self' DNA, e.g. from bacteriophages. Restriction enzymes detect 'non-self' methylation patterns and cut foreign DNA. Furthermore, methylation of promoter regions can influence gene expression and hence affect various phenotypes. In this study, we determined the methylated motifs of four strains from the Burkholderia cepacia complex of opportunistic pathogens. We deleted all genes encoding the restriction and modification components in one of these strains, Burkholderia cenocepacia H111. It is shown that DNA methylation affects various phenotypic traits, the most noteworthy being lysogenicity of a bacteriophage and maintenance of a virulence megaplasmid.

Plasmid-chromosome cross-talks

Plasmids can be acquired by recipient bacteria at a significant cost while conferring them advantageous traits. To counterbalance the costs of plasmid carriage, both plasmids and host bacteria have developed a tight regulatory network that may involve a cross-talk between the chromosome and the plasmids. Although plasmid regulation by chromosomal regulators is generally well known, chromosome regulation by plasmid has been far less investigated. Yet, a growing number of studies have highlighted an impact of plasmids on their host bacteria. Here, we describe the plasmid-chromosome cross-talk from the plasmid point of view. We summarize data about the chromosomal adaptive mutations generated by plasmid carriage; the impact of the loss of a domesticated plasmid or the gain of a new plasmid. Then, we present the control of plasmid-encoded regulators on chromosomal gene expression. The involvement of regulators homologous to chromosome-encoded proteins is illustrated by the H-NS-like proteins, and by the Rap-Phr system. Finally, plasmid-specific regulators of chromosomal gene expression are presented, which highlight the involvement of transcription factors and sRNAs. A comprehensive analysis of the mechanisms that allow a given plasmid to impact the chromosome of bacterium will help to understand the tight cross-talk between plasmids and the chromosome.

Genomic analyses of Burkholderia cenocepacia reveal multiple species with differential host-adaptation to plants and humans

Background

Burkholderia cenocepacia is a human opportunistic pathogen causing devastating symptoms in patients suffering from immunodeficiency and cystic fibrosis. Out of the 303 B. cenocepacia strains with available genomes, the large majority were isolated from a clinical context. However, several isolates originate from other environmental sources ranging from aerosols to plant endosphere. Plants can represent reservoirs for human infections as some pathogens can survive and sometimes proliferate in the rhizosphere. We therefore investigated if B. cenocepacia had the same potential.

Results

We selected genome sequences from 31 different strains, representative of the diversity of ecological niches of B. cenocepacia, and conducted comparative genomic analyses in the aim of finding specific niche or host-related genetic determinants. Phylogenetic analyses and whole genome average nucleotide identity suggest that strains, registered as B. cenocepacia, belong to at least two different species. Core-genome analyses show that the clade enriched in environmental isolates lacks multiple key virulence factors, which are conserved in the sister clade where most clinical isolates fall, including the highly virulent ET12 lineage. Similarly, several plant associated genes display an opposite distribution between the two clades. Finally, we suggest that B. cenocepacia underwent a host jump from plants/environment to animals, as supported by the phylogenetic analysis. We eventually propose a name for the new species that lacks several genetic traits involved in human virulence.

Conclusion

Regardless of the method used, our studies resulted in a disunited perspective of the B. cenocepacia species. Strains currently affiliated to this taxon belong to at least two distinct species, one having lost several determining animal virulence factors.

Destabilization of the Tumor-Inducing Plasmid from an Octopine-Type Agrobacterium tumefaciens Lineage Drives a Large Deletion in the Co-resident At Megaplasmid

Bacteria with multi-replicon genome organizations, including members of the family Rhizobiaceae, often carry a variety of niche-associated functions on large plasmids. While evidence exists for cross-replicon interactions and co-evolution between replicons in many of these systems, remarkable strain-to-strain variation is also observed for extrachromosomal elements, suggesting increased genetic plasticity. Here, we show that curing of the tumor-inducing virulence plasmid (pTi) of an octopine-type Agrobacterium tumefaciens lineage leads to a large deletion in the co-resident At megaplasmid (pAt). The deletion event is mediated by a repetitive IS-element, IS66, and results in a variety of environment-dependent fitness consequences, including loss of independent conjugal transfer of the plasmid. Interestingly, a related and otherwise wild-type A. tumefaciens strain is missing exactly the same large pAt segment as the pAt deletion derivatives, suggesting a similar event over its natural history. Overall, the findings presented here uncover a novel genetic interaction between the two large plasmids of A. tumefaciens and provide evidence for cross-replicon integration and co-evolution of these plasmids.

Chromids Aid Genome Expansion and Functional Diversification in the Family Burkholderiaceae

Multipartite genomes, containing at least two large replicons, are found in diverse bacteria; however, the advantage of this genome structure remains incompletely understood. Here, we perform comparative genomics of hundreds of finished β-proteobacterial genomes to gain insights into the role and emergence of multipartite genomes. Almost all essential secondary replicons (chromids) of the β-proteobacteria are found in the family Burkholderiaceae. These replicons arose from just two plasmid acquisition events, and they were likely stabilized early in their evolution by the presence of core genes. On average, Burkholderiaceae genera with multipartite genomes had a larger total genome size, but smaller chromosome, than genera without secondary replicons. Pangenome-level functional enrichment analyses suggested that interreplicon functional biases are partially driven by the enrichment of secondary replicons in the accessory pangenome fraction. Nevertheless, the small overlap in orthologous groups present in each replicon's pangenome indicated a clear functional separation of the replicons. Chromids appeared biased to environmental adaptation, as the functional categories enriched on chromids were also overrepresented on the chromosomes of the environmental genera (Paraburkholderia and Cupriavidus) compared with the pathogenic genera (Burkholderia and Ralstonia). Using ancestral state reconstruction, it was predicted that the rate of accumulation of modern-day genes by chromids was more rapid than the rate of gene accumulation by the chromosomes. Overall, the data are consistent with a model where the primary advantage of secondary replicons is in facilitating increased rates of gene acquisition through horizontal gene transfer, consequently resulting in replicons enriched in genes associated with adaptation to novel environments.

Genome mining identifies cepacin as a plant-protective metabolite of the biopesticidal bacterium Burkholderia ambifaria

Beneficial microorganisms are widely used in agriculture for control of plant pathogens, but a lack of efficacy and safety information has limited the exploitation of multiple promising biopesticides. We applied phylogeny-led genome mining, metabolite analyses and biological control assays to define the efficacy of Burkholderia ambifaria, a naturally beneficial bacterium with proven biocontrol properties but potential pathogenic risk. A panel of 64 B. ambifaria strains demonstrated significant antimicrobial activity against priority plant pathogens. Genome sequencing, specialized metabolite biosynthetic gene cluster mining and metabolite analysis revealed an armoury of known and unknown pathways within B. ambifaria. The biosynthetic gene cluster responsible for the production of the metabolite cepacin was identified and directly shown to mediate protection of germinating crops against Pythium damping-off disease. B. ambifaria maintained biopesticidal protection and overall fitness in the soil after deletion of its third replicon, a non-essential plasmid associated with virulence in Burkholderia cepacia complex bacteria. Removal of the third replicon reduced B. ambifaria persistence in a murine respiratory infection model. Here, we show that by using interdisciplinary phylogenomic, metabolomic and functional approaches, the mode of action of natural biological control agents related to pathogens can be systematically established to facilitate their future exploitation.

Potential of the Burkholderia cepacia Complex to Produce 4-Hydroxy-3-Methyl-2-Alkyquinolines

A few Burkholderia species, especially Burkholderia pseudomallei, Burkholderia thailandensis, Burkholderia ambifaria, and Burkholderia cepacia, are known to produce and release various 4-hydroxy-3-methyl-2-alkylquinolines (HMAQs), a family of molecules analogous to the 4-hydroxy-2-alkylquinolines [aka 2-n-alkyl-4(1H)-quinolones] of Pseudomonas aeruginosa, which include the Pseudomonas quinolone signal (PQS). However, while these exoproducts play several roles in P. aeruginosa virulence and survival, the available literature is very limited on their distribution and function in Burkholderia. In this perspective article, we studied the distribution of the hmqABCDEFG operon, which encodes the enzymes involved in the biosynthesis of HMAQs, in the Burkholderia cepacia complex (Bcc) group. Based on the available sequence data, about one third of Bcc species carry a homolog of the hmqABCDEFG, and not all sequenced strains in a given species possess this operon. Looking at the synteny of genes surrounding the hmqABCDEFG operon, we found that for some species, the operon seems to have been deleted or replaced by other genes. Finally, we review the literature on the possible function of HMAQs. Understanding the Hmq system may provide clues concerning their functions in Bcc.

Genomic Features and Insights into the Taxonomy, Virulence, and Benevolence of Plant-Associated Burkholderia Species

The members of the Burkholderia genus are characterized by high versatility and adaptability to various ecological niches. With the availability of the genome sequences of numerous species of Burkholderia, many studies have been conducted to elucidate the unique features of this exceptional group of bacteria. Genomic and metabolic plasticity are common among Burkholderia species, as evidenced by their relatively large multi-replicon genomes that are rich in insertion sequences and genomic islands and contain a high proportion of coding regions. Such unique features could explain their adaptability to various habitats and their versatile lifestyles, which are reflected in a multiplicity of species including free-living rhizospheric bacteria, plant endosymbionts, legume nodulators, and plant pathogens. The phytopathogenic Burkholderia group encompasses several pathogens representing threats to important agriculture crops such as rice. Contrarily, plant-beneficial Burkholderia have also been reported, which have symbiotic and growth-promoting roles. In this review, the taxonomy of Burkholderia is discussed emphasizing the recent updates and the contributions of genomic studies to precise taxonomic positioning. Moreover, genomic and functional studies on Burkholderia are reviewed and insights are provided into the mechanisms underlying the virulence and benevolence of phytopathogenic and plant-beneficial Burkholderia, respectively, on the basis of cutting-edge knowledge.

The afc antifungal activity cluster, which is under tight regulatory control of ShvR, is essential for transition from intracellular persistence of Burkholderia cenocepacia to acute pro-inflammatory infection

The opportunistic pathogen Burkholderia cenocepacia is particularly life-threatening for cystic fibrosis (CF) patients. Chronic lung infections with these bacteria can rapidly develop into fatal pulmonary necrosis and septicaemia. We have recently shown that macrophages are a critical site for replication of B. cenocepacia K56-2 and the induction of fatal pro-inflammatory responses using a zebrafish infection model. Here, we show that ShvR, a LysR-type transcriptional regulator that is important for biofilm formation, rough colony morphotype and inflammation in a rat lung infection model, is also required for the induction of fatal pro-inflammatory responses in zebrafish larvae. ShvR was not essential, however, for bacterial survival and replication in macrophages. Temporal, rhamnose-induced restoration of shvR expression in the shvR mutant during intramacrophage stages unequivocally demonstrated a key role for ShvR in transition from intracellular persistence to acute fatal pro-inflammatory disease. ShvR has been previously shown to tightly control the expression of the adjacent afc gene cluster, which specifies the synthesis of a lipopeptide with antifungal activity. Mutation of afcE, encoding an acyl-CoA dehydrogenase, has been shown to give similar phenotypes as the shvR mutant. We found that, like shvR, afcE is also critical for the switch from intracellular persistence to fatal infection in zebrafish. The closely related B. cenocepacia H111 has been shown to be less virulent than K56-2 in several infection models, including Galleria mellonella and rats. Interestingly, constitutive expression of shvR in H111 increased virulence in zebrafish larvae to almost K56-2 levels in a manner that absolutely required afc. These data confirm a critical role for afc in acute virulence caused by B. cenocepacia that depends on strain-specific regulatory control by ShvR. We propose that ShvR and AFC are important virulence factors of the more virulent Bcc species, either through pro-inflammatory effects of the lipopeptide AFC, or through AFC-dependent membrane properties.

Creation and Characterization of a Genomically Hybrid Strain in the Nitrogen-Fixing Symbiotic Bacterium Sinorhizobium meliloti

Many bacteria, often associated with eukaryotic hosts and of relevance for biotechnological applications, harbor a multipartite genome composed of more than one replicon. Biotechnologically relevant phenotypes are often encoded by genes residing on the secondary replicons. A synthetic biology approach to developing enhanced strains for biotechnological purposes could therefore involve merging pieces or entire replicons from multiple strains into a single genome. Here we report the creation of a genomic hybrid strain in a model multipartite genome species, the plant-symbiotic bacterium Sinorhizobium meliloti. We term this strain as cis-hybrid, since it is produced by genomic material coming from the same species' pangenome. In particular, we moved the secondary replicon pSymA (accounting for nearly 20% of total genome content) from a donor S. meliloti strain to an acceptor strain. The cis-hybrid strain was screened for a panel of complex phenotypes (carbon/nitrogen utilization phenotypes, intra- and extracellular metabolomes, symbiosis, and various microbiological tests). Additionally, metabolic network reconstruction and constraint-based modeling were employed for in silico prediction of metabolic flux reorganization. Phenotypes of the cis-hybrid strain were in good agreement with those of both parental strains. Interestingly, the symbiotic phenotype showed a marked cultivar-specific improvement with the cis-hybrid strains compared to both parental strains. These results provide a proof-of-principle for the feasibility of genome-wide replicon-based remodelling of bacterial strains for improved biotechnological applications in precision agriculture.

Periodic Variation of Mutation Rates in Bacterial Genomes Associated with Replication Timing

The causes and consequences of spatiotemporal variation in mutation rates remain to be explored in nearly all organisms. Here we examine relationships between local mutation rates and replication timing in three bacterial species whose genomes have multiple chromosomes: Vibrio fischeri, Vibrio cholerae, and Burkholderia cenocepacia. Following five mutation accumulation experiments with these bacteria conducted in the near absence of natural selection, the genomes of clones from each lineage were sequenced and analyzed to identify variation in mutation rates and spectra. In lineages lacking mismatch repair, base substitution mutation rates vary in a mirrored wave-like pattern on opposing replichores of the large chromosomes of V. fischeri and V. cholerae, where concurrently replicated regions experience similar base substitution mutation rates. The base substitution mutation rates on the small chromosome are less variable in both species but occur at similar rates to those in the concurrently replicated regions of the large chromosome. Neither nucleotide composition nor frequency of nucleotide motifs differed among regions experiencing high and low base substitution rates, which along with the inferred ~800-kb wave period suggests that the source of the periodicity is not sequence specific but rather a systematic process related to the cell cycle. These results support the notion that base substitution mutation rates are likely to vary systematically across many bacterial genomes, which exposes certain genes to elevated deleterious mutational load.

That mutation rates vary within bacterial genomes is well known, but the detailed study of these biases has been made possible only recently with contemporary sequencing methods. We applied these methods to understand how bacterial genomes with multiple chromosomes, like those of Vibrio and Burkholderia, might experience heterogeneous mutation rates because of their unusual replication and the greater genetic diversity found on smaller chromosomes. This study captured thousands of mutations and revealed wave-like rate variation that is synchronized with replication timing and not explained by sequence context. The scale of this rate variation over hundreds of kilobases of DNA strongly suggests that a temporally regulated cellular process may generate wave-like variation in mutation risk. These findings add to our understanding of how mutation risk is distributed across bacterial and likely also eukaryotic genomes, owing to their highly conserved replication and repair machinery.

Replicate Once Per Cell Cycle: Replication Control of Secondary Chromosomes

Faithful vertical transmission of genetic information, especially of essential core genes, is a prerequisite for bacterial survival. Hence, replication of all the replicons is tightly controlled to ensure that all daughter cells get the same genome copy as their mother cell. Essential core genes are very often carried by the main chromosome. However they can occasionally be found on secondary chromosomes, recently renamed chromids. Chromids have evolved from non-essential megaplasmids, and further acquired essential core genes and a genomic signature closed to that of the main chromosome. All chromids carry a plasmidic replication origin, belonging so far to either the iterons or repABC type. Based on these differences, two categories of chromids have been distinguished. In this review, we focus on the replication initiation controls of these two types of chromids. We show that the sophisticated mechanisms controlling their replication evolved from their plasmid counterparts to allow a timely controlled replication, occurring once per cell cycle.

Metabolomics and Transcriptomics Identify Multiple Downstream Targets of Paraburkholderia phymatum σ54 During Symbiosis with Phaseolus vulgaris

RpoN (or σ⁵⁴) is the key sigma factor for the regulation of transcription of nitrogen fixation genes in diazotrophic bacteria, which include α- and β-rhizobia. Our previous studies showed that an rpoN mutant of the β-rhizobial strain Paraburkholderia phymatum STM815^T formed root nodules on Phaseolus vulgaris cv. Negro jamapa, which were unable to reduce atmospheric nitrogen into ammonia. In an effort to further characterize the RpoN regulon of P. phymatum, transcriptomics was combined with a powerful metabolomics approach. The metabolome of P. vulgaris root nodules infected by a P. phymatumrpoN Fix⁻ mutant revealed statistically significant metabolic changes compared to wild-type Fix⁺ nodules, including reduced amounts of chorismate and elevated levels of flavonoids. A transcriptome analysis on Fix⁻ and Fix⁺ nodules—combined with a search for RpoN binding sequences in promoter regions of regulated genes—confirmed the expected control of σ⁵⁴ on nitrogen fixation genes in nodules. The transcriptomic data also allowed us to identify additional target genes, whose differential expression was able to explain the observed metabolite changes in numerous cases. Moreover, the genes encoding the two-component regulatory system NtrBC were downregulated in root nodules induced by the rpoN mutant, and contained a putative RpoN binding motif in their promoter region, suggesting direct regulation. The construction and characterization of an ntrB mutant strain revealed impaired nitrogen assimilation in free-living conditions, as well as a noticeable symbiotic phenotype, as fewer but heavier nodules were formed on P. vulgaris roots.

Biosynthesis of fragin is controlled by a novel quorum sensing signal

Members of the diazeniumdiolate class of natural compounds show potential for drug development because of their antifungal, antibacterial, antiviral, and antitumor activities. Yet, their biosynthesis has remained elusive to date. Here, we identify a gene cluster directing the biosynthesis of the diazeniumdiolate compound fragin in Burkholderia cenocepacia H111. We provide evidence that fragin is a metallophore and that metal chelation is the molecular basis of its antifungal activity. A subset of the fragin biosynthetic genes is involved in the synthesis of a previously undescribed cell-to-cell signal molecule, valdiazen. RNA-Seq analyses reveal that valdiazen controls fragin biosynthesis and affects the expression of more than 100 genes. Homologs of the valdiazen biosynthesis genes are found in various bacteria, suggesting that valdiazen-like compounds may constitute a new class of signal molecules. We use structural information, in silico prediction of enzymatic functions and biochemical data to propose a biosynthesis route for fragin and valdiazen.

Fragin is a diazeniumdiolate metabolite with antifungal activity, produced by some bacteria. Here, Jenul et al. show that metal chelation is the molecular basis of fragin’s antifungal activity, and that a gene cluster directing fragin biosynthesis is also involved in the synthesis of a signal molecule.

A putative lateral flagella of the cystic fibrosis pathogen Burkholderia dolosa regulates swimming motility and host cytokine production

Burkholderia dolosa caused an outbreak in the cystic fibrosis clinic at Boston Children's Hospital and was associated with high mortality in these patients. This species is part of a larger complex of opportunistic pathogens known as the Burkholderia cepacia complex (Bcc). Compared to other species in the Bcc, B. dolosa is highly transmissible; thus understanding its virulence mechanisms is important for preventing future outbreaks. The genome of one of the outbreak strains, AU0158, revealed a homolog of the lafA gene encoding a putative lateral flagellin, which, in other non-Bcc species, is used for movement on solid surfaces, attachment to host cells, or movement inside host cells. Here, we analyzed the conservation of the lafA gene and protein sequences, which are distinct from those of the polar flagella, and found lafA homologs to be present in numerous β-proteobacteria but notably absent from most other Bcc species. A lafA deletion mutant in B. dolosa showed a greater swimming motility than wild-type due to an increase in the number of polar flagella, but did not appear to contribute to biofilm formation, host cell invasion, or murine lung colonization or persistence over time. However, the lafA gene was important for cytokine production in human peripheral blood mononuclear cells, suggesting it may have a role in recognition by the human immune response.

What matters in chronic Burkholderia cenocepacia infection in cystic fibrosis: Insights from comparative genomics

Burkholderia cenocepacia causes severe pulmonary infections in cystic fibrosis (CF) patients. Since the bacterium is virtually untreatable by antibiotics, chronic infections persist for years and might develop into fatal septic pneumonia (cepacia syndrome, CS). To devise new strategies to combat chronic B. cenocepacia infections, it is essential to obtain comprehensive knowledge about their pathogenesis. We conducted a comparative genomic analysis of 32 Czech isolates of epidemic clone B. cenocepacia ST32 isolated from various stages of chronic infection in 8 CF patients. High numbers of large-scale deletions were found to occur during chronic infection, affecting preferentially genomic islands and nonessential replicons. Recombination between insertion sequences (IS) was inferred as the mechanism behind deletion formation; the most numerous IS group was specific for the ST32 clone and has undergone transposition burst since its divergence. Genes functionally related to transition metal metabolism were identified as hotspots for deletions and IS insertions. This functional category was also represented among genes where nonsynonymous point mutations and indels occurred parallelly among patients. Another category exhibiting parallel mutations was oxidative stress protection; mutations in catalase KatG resulted in impaired detoxification of hydrogen peroxide. Deep sequencing revealed substantial polymorphism in genes of both categories within the sputum B. cenocepacia ST32 populations, indicating extensive adaptive evolution. Neither oxidative stress response nor transition metal metabolism genes were previously reported to undergo parallel evolution during chronic CF infection. Mutations in katG and copper metabolism genes were overrepresented in patients where chronic infection developed into CS. Among professional phagocytes, macrophages use both hydrogen peroxide and copper for their bactericidal activity; our results thus tentatively point to macrophages as suspects in pathogenesis towards the fatal CS.

The large Burkholderia cenocepacia populations which persist in cystic fibrosis lungs during many years of chronic infections have an inherent potential for adaptive evolution. The results provided by comparative genomics are key in understanding the processes involved. Mutational events which have taken place allow us to deductively reconstruct the history of chronic infection and to identify driving forces acting upon the bacteria. Beyond the conventional point mutation analysis of next generation sequencing data, we observed interesting phenomena such as large deletions and transposable element movement which represent another facet of adaptive evolution of B. cenocepacia during chronic infection. We also found, unexpectedly, that adaptive evolution in B. cenocepacia strain ST32 affects a set of genes conspicuously different from related species B. dolosa; these appear to be linked to host immune response. Our study provides clues to the complex puzzle of chronic B. cenocepacia infection establishment, persistence and outcome in cystic fibrosis.

Subfunctionalization influences the expansion of bacterial multidrug antibiotic resistance

Background

Antibiotic resistance is a major problem for human health. Multidrug resistance efflux pumps, especially those of the Resistance-Nodulation-Cell Division (RND) family, are major contributors to high-level antibiotic resistance in Gram-negative bacteria. Most bacterial genomes contain several copies of the different classes of multidrug resistance efflux pumps. Gene duplication and gain of function by the duplicate copies of multidrug resistance efflux pump genes plays a key role in the expansion and diversification of drug-resistance mechanisms.

Results

We used two members of the Burkholderia RND superfamily as models to understand how duplication events affect the antibiotic resistance of these strains. First, we analyzed the conservation and distribution of these two RND systems and their regulators across the Burkholderia genus. Through genetic manipulations, we identified both the exact substrate range of these transporters and their eventual interchangeability. We also performed a directed evolution experiment, combined with next generation sequencing, to evaluate the role of antibiotics in the activation of the expression of these systems. Together, our results indicate that the first step to diversify the functions of these pumps arises from changes in their regulation (subfunctionalization) instead of functional mutations. Further, these pumps could rewire their regulation to respond to antibiotics, thus maintaining high genomic plasticity.

Conclusions

Studying the regulatory network that controls the expression of the RND pumps will help understand and eventually control the development and expansion of drug resistance.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4222-4) contains supplementary material, which is available to authorized users.

The Essential Genome of Burkholderia cenocepacia H111

The study of the minimum set of genes required to sustain life is a fundamental question in biological research. Recent studies on bacterial essential genes suggested that between 350 and 700 genes are essential to support autonomous bacterial cell growth. Essential genes are of interest as potential new antimicrobial drug targets; hence, our aim was to identify the essential genome of the cystic fibrosis (CF) isolate Burkholderia cenocepacia H111. Using a transposon sequencing (Tn-Seq) approach, we identified essential genes required for growth in rich medium under aerobic and microoxic conditions as well as in a defined minimal medium with citrate as a sole carbon source. Our analysis suggests that 398 genes are required for autonomous growth in rich medium, a number that represents only around 5% of the predicted genes of this bacterium. Five hundred twenty-six genes were required to support growth in minimal medium, and 434 genes were essential under microoxic conditions (0.5% O₂). A comparison of these data sets identified 339 genes that represent the minimal set of essential genes required for growth under all conditions tested and can be considered the core essential genome of B. cenocepacia H111. The majority of essential genes were found to be located on chromosome 1, and few such genes were located on chromosome 2, where most of them were clustered in one region. This gene cluster is fully conserved in all Burkholderia species but is present on chromosome 1 in members of the closely related genus Ralstonia, suggesting that the transfer of these essential genes to chromosome 2 in a common ancestor contributed toward the separation of the two genera.IMPORTANCE Transposon sequencing (Tn-Seq) is a powerful method used to identify genes that are essential for autonomous growth under various conditions. In this study, we have identified a set of "core essential genes" that are required for growth under multiple conditions, and these genes represent potential antimicrobial targets. We also identified genes specifically required for growth under low-oxygen and nutrient-limited environments. We generated conditional mutants to verify the results of our Tn-Seq analysis and demonstrate that one of the identified genes was not essential per se but was an artifact of the construction of the mutant library. We also present verified examples of genes that were not truly essential but, when inactivated, showed a growth defect. These examples have identified so-far-underestimated shortcomings of this powerful method.

The Divided Bacterial Genome: Structure, Function, and Evolution

Approximately 10% of bacterial genomes are split between two or more large DNA fragments, a genome architecture referred to as a multipartite genome. This multipartite organization is found in many important organisms, including plant symbionts, such as the nitrogen-fixing rhizobia, and plant, animal, and human pathogens, including the genera Brucella, Vibrio, and Burkholderia. The availability of many complete bacterial genome sequences means that we can now examine on a broad scale the characteristics of the different types of DNA molecules in a genome. Recent work has begun to shed light on the unique properties of each class of replicon, the unique functional role of chromosomal and nonchromosomal DNA molecules, and how the exploitation of novel niches may have driven the evolution of the multipartite genome. The aims of this review are to (i) outline the literature regarding bacterial genomes that are divided into multiple fragments, (ii) provide a meta-analysis of completed bacterial genomes from 1,708 species as a way of reviewing the abundant information present in these genome sequences, and (iii) provide an encompassing model to explain the evolution and function of the multipartite genome structure. This review covers, among other topics, salient genome terminology; mechanisms of multipartite genome formation; the phylogenetic distribution of multipartite genomes; how each part of a genome differs with respect to genomic signatures, genetic variability, and gene functional annotation; how each DNA molecule may interact; as well as the costs and benefits of this genome structure.

Phylogeographic, genomic, and meropenem susceptibility analysis of Burkholderia ubonensis

The bacterium Burkholderia ubonensis is commonly co-isolated from environmental specimens harbouring the melioidosis pathogen, Burkholderia pseudomallei. B. ubonensis has been reported in northern Australia and Thailand but not North America, suggesting similar geographic distribution to B. pseudomallei. Unlike most other Burkholderia cepacia complex (Bcc) species, B. ubonensis is considered non-pathogenic, although its virulence potential has not been tested. Antibiotic resistance in B. ubonensis, particularly towards drugs used to treat the most severe B. pseudomallei infections, has also been poorly characterised. This study examined the population biology of B. ubonensis, and includes the first reported isolates from the Caribbean. Phylogenomic analysis of 264 B. ubonensis genomes identified distinct clades that corresponded with geographic origin, similar to B. pseudomallei. A small proportion (4%) of strains lacked the 920kb chromosome III replicon, with discordance of presence/absence amongst genetically highly related strains, demonstrating that the third chromosome of B. ubonensis, like other Bcc species, probably encodes for a nonessential pC3 megaplasmid. Multilocus sequence typing using the B. pseudomallei scheme revealed that one-third of strains lack the "housekeeping" narK locus. In comparison, all strains could be genotyped using the Bcc scheme. Several strains possessed high-level meropenem resistance (≥32 μg/mL), a concern due to potential transmission of this phenotype to B. pseudomallei. In silico analysis uncovered a high degree of heterogeneity among the lipopolysaccharide O-antigen cluster loci, with at least 35 different variants identified. Finally, we show that Asian B. ubonensis isolate RF23-BP41 is avirulent in the BALB/c mouse model via a subcutaneous route of infection. Our results provide several new insights into the biology of this understudied species.

NtrC-dependent control of exopolysaccharide synthesis and motility in Burkholderia cenocepacia H111

Burkholderia cenocepacia is a versatile opportunistic pathogen that survives in a wide variety of environments, which can be limited in nutrients such as nitrogen. We have previously shown that the sigma factor σ54 is involved in the control of nitrogen assimilation and virulence in B. cenocepacia H111. In this work, we investigated the role of the σ54 enhancer binding protein NtrC in response to nitrogen limitation and in the pathogenicity of H111. Of 95 alternative nitrogen sources tested the ntrC showed defects in the utilisation of nitrate, urea, L-citrulline, acetamide, DL-lactamide, allantoin and parabanic acid. RNA-Seq and phenotypic analyses of an ntrC mutant strain showed that NtrC positively regulates two important phenotypic traits: exopolysaccharide (EPS) production and motility. However, the ntrC mutant was not attenuated in C. elegans virulence.

Use of Synthetic Hybrid Strains To Determine the Role of Replicon 3 in Virulence of the Burkholderia cepacia Complex

The Burkholderia cepacia complex (Bcc) displays a wealth of metabolic diversity with great biotechnological potential, but the utilization of these bacteria is limited by their opportunistic pathogenicity to humans. The third replicon of the Bcc, megaplasmid pC3 (0.5 to 1.4 Mb, previously chromosome 3), is important for various phenotypes, including virulence, antifungal, and proteolytic activities and the utilization of certain substrates. Approximately half of plasmid pC3 is well conserved throughout sequenced Bcc members, while the other half is not. To better locate the regions responsible for the key phenotypes, pC3 mutant derivatives of Burkholderia cenocepacia H111 carrying large deletions (up to 0.58 Mb) were constructed with the aid of the FLP-FRT (FRT, flippase recognition target) recombination system from Saccharomyces cerevisiae The conserved region was shown to confer near-full virulence in both Caenorhabditis elegans and Galleria mellonella infection models. Antifungal activity was unexpectedly independent of the part of pC3 bearing a previously identified antifungal gene cluster, while proteolytic activity was dependent on the nonconserved part of pC3, which encodes the ZmpA protease. To investigate to what degree pC3-encoded functions are dependent on chromosomally encoded functions, we transferred pC3 from Burkholderia cenocepacia K56-2 and Burkholderia lata 383 into other pC3-cured Bcc members. We found that although pC3 is highly important for virulence, it was the genetic background of the recipient that determined the pathogenicity level of the hybrid strain. Furthermore, we found that important phenotypes, such as antifungal activity, proteolytic activity, and some substrate utilization capabilities, can be transferred between Bcc members using pC3.IMPORTANCE The Burkholderia cepacia complex (Bcc) is a group of closely related bacteria with great biotechnological potential. Some strains produce potent antifungal compounds and can promote plant growth or degrade environmental pollutants. However, their agricultural potential is limited by their opportunistic pathogenicity, particularly for cystic fibrosis patients. Despite much study, their virulence remains poorly understood. The third replicon, pC3, which is present in all Bcc isolates and is important for pathogenicity, stress resistance, and the production of antifungal compounds, has recently been reclassified from a chromosome to a megaplasmid. In this study, we identified regions on pC3 important for virulence and antifungal activity and investigated the role of the chromosomal background for the function of pC3 by exchanging the megaplasmid between different Bcc members. Our results may open a new avenue for the construction of antifungal but nonpathogenic Burkholderia hybrids. Such strains may have great potential as biocontrol strains for protecting fungus-borne diseases of plant crops.

High intracellular c-di-GMP levels antagonize quorum sensing and virulence gene expression in Burkholderia cenocepacia H111

The opportunistic human pathogen Burkholderia cenocepacia H111 uses two chemically distinct signal molecules for controlling gene expression in a cell density-dependent manner: N-acyl-homoserine lactones (AHLs) and cis-2-dodecenoic acid (BDSF). Binding of BDSF to its cognate receptor RpfR lowers the intracellular c-di-GMP level, which in turn leads to differential expression of target genes. In this study we analysed the transcriptional profile of B. cenocepacia H111 upon artificially altering the cellular c-di-GMP level. One hundred and eleven genes were shown to be differentially expressed, 96 of which were downregulated at a high c-di-GMP concentration. Our analysis revealed that the BDSF, AHL and c-di-GMP regulons overlap for the regulation of 24 genes and that a high c-di-GMP level suppresses expression of AHL-regulated genes. Phenotypic analyses confirmed changes in the expression of virulence factors, the production of AHL signal molecules and the biosynthesis of different biofilm matrix components upon altered c-di-GMP levels. We also demonstrate that the intracellular c-di-GMP level determines the virulence of B. cenocepacia to Caenorhabditis elegans and Galleria mellonella.

Comparative genomics of Burkholderia multivorans, a ubiquitous pathogen with a highly conserved genomic structure

The natural environment serves as a reservoir of opportunistic pathogens. A well-established method for studying the epidemiology of such opportunists is multilocus sequence typing, which in many cases has defined strains predisposed to causing infection. Burkholderia multivorans is an important pathogen in people with cystic fibrosis (CF) and its epidemiology suggests that strains are acquired from non-human sources such as the natural environment. This raises the central question of whether the isolation source (CF or environment) or the multilocus sequence type (ST) of B. multivorans better predicts their genomic content and functionality. We identified four pairs of B. multivorans isolates, representing distinct STs and consisting of one CF and one environmental isolate each. All genomes were sequenced using the PacBio SMRT sequencing technology, which resulted in eight high-quality B. multivorans genome assemblies. The present study demonstrated that the genomic structure of the examined B. multivorans STs is highly conserved and that the B. multivorans genomic lineages are defined by their ST. Orthologous protein families were not uniformly distributed among chromosomes, with core orthologs being enriched on the primary chromosome and ST-specific orthologs being enriched on the second and third chromosome. The ST-specific orthologs were enriched in genes involved in defense mechanisms and secondary metabolism, corroborating the strain-specificity of these virulence characteristics. Finally, the same B. multivorans genomic lineages occur in both CF and environmental samples and on different continents, demonstrating their ubiquity and evolutionary persistence.

An In vitro Study of Bio-Control and Plant Growth Promotion Potential of Salicaceae Endophytes

Microbial communities in the endosphere of Salicaceae plants, poplar (Populus trichocarpa) and willow (Salix sitchensis), have been demonstrated to be important for plant growth promotion, protection from biotic and abiotic stresses, and degradation of toxic compounds. Our study aimed to investigate bio-control activities of Salicaceae endophytes against various soil borne plant pathogens including Rhizoctonia solani AG-8, Fusarium culmorum, Gaeumannomyces graminis var. tritici, and Pythium ultimum. Additionally, different plant growth promoting traits such as biological nitrogen fixation (BNF), indole-3-acetic acid (IAA) biosynthesis, phosphate solubilization, and siderophore production were assessed in all bio-control positive strains. Burkholderia, Rahnella, Pseudomonas, and Curtobacterium were major endophyte genera that showed bio-control activities in the in-vitro assays. The bio-control activities of Burkholderia strains were stronger across all tested plant pathogens as compared to other stains. Genomes of sequenced Burkholderia strains WP40 and WP42 were surveyed to identify the putative genes involved in the bio-control activities. The ocf and hcnABC gene clusters responsible for biosynthesis of the anti-fungal metabolites, occidiofungin and hydrogen cyanide, are present in the genomes of WP40 and WP42. Nearly all endophyte strains showing the bio-control activities produced IAA, solubilized tricalcium phosphate, and synthesized siderophores in the culture medium. Moreover, some strains reduced acetylene into ethylene in the acetylene reduction assay, a common assay used for BNF. Salicaceae endophytes could be useful for bio-control of various plant pathogens, and plant growth promotion possibly through the mechanisms of BNF, IAA production, and nutrient acquisition.

1H-NMR-Based Endometabolome Profiles of Burkholderia cenocepacia Clonal Variants Retrieved from a Cystic Fibrosis Patient during Chronic Infection

During cystic fibrosis (CF) chronic lung infections, bacteria of the Burkholderia cepacia complex (Bcc) are exposed for several years to a stressful and changing environment. These environmental challenges results in genetic changes of the initial infecting strain with the consequent diversification of genotypes and phenotypes. The exploitation of functional and comparative genomic approaches has suggested that such diversification is associated with massive metabolic remodeling but these alterations are poorly understood. In the present work, we have explored a high resolution ¹H-NMR-based metabolomic approach coupled to multivariate analysis to compare the endometabolome of three B. cenocepacia clonal variants retrieved from a CF patient from the onset of infection (IST439) until death with cepacia syndrome after 3.5 years (IST4113 and IST4134), to complement former proteomic and transcriptomic analyses. A fourth clonal variant (IST4129) retrieved from the same CF patient when the clinical condition worsened during the last months of life, was also examined since it was found to lack the third replicon. The metabolomic profiles obtained, based on the complete ¹H-NMR spectra, highlight the separation of the four clonal variants examined, the most distinct profile corresponding to IST4129. Results indicate a variable content of several amino acids in the different isolates examined and suggest that glycolysis and the glyoxylate shunt are favored in late variants. Moreover, the concentration of two metabolites with demonstrated cellular protective functions against stress, glycine-betaine and trehalose, is different in the different isolates examined. However, no clear correlation could be established between their content and stress tolerance. For example, IST4113, previously found to be the most resistant variant to antimicrobials of different classes, exhibits low levels of trehalose and glycine-betaine but the highest resistance to heat and oxidative stress. Also, IST4129, with a high level of glycine-betaine but lacking the third replicon, previously associated with stress resistance and virulence, exhibits the highest susceptibility to all the stresses tested. Taken together, results from this study provide insights into the metabolic diversification of B. cenocepacia clonal variants during long-term infection of the CF airways.

Complete Genome Sequences for Three Chromosomes of the Burkholderia stabilis Type Strain (ATCC BAA-67)

We report here the complete annotated genome sequence of the Burkholderia stabilis type strain ATCC BAA-67. There were three circular chromosomes with a combined size of 8,527,947 bp and G+C composition of 66.4%. These characteristics closely resemble the genomes of other sequenced members of the Burkholderia cepacia complex.

The Fitness Effects of Spontaneous Mutations Nearly Unseen by Selection in a Bacterium with Multiple Chromosomes

Mutation accumulation (MA) experiments employ the strategy of minimizing the population size of evolving lineages to greatly reduce effects of selection on newly arising mutations. Thus, most mutations fix within MA lines independently of their fitness effects. This approach, more recently combined with genome sequencing, has detailed the rates, spectra, and biases of different mutational processes. However, a quantitative understanding of the fitness effects of mutations virtually unseen by selection has remained an untapped opportunity. Here, we analyzed the fitness of 43 sequenced MA lines of the multi-chromosome bacterium Burkholderia cenocepacia that had each undergone 5554 generations of MA and accumulated an average of 6.73 spontaneous mutations. Most lineages exhibited either neutral or deleterious fitness in three different environments in comparison with their common ancestor. The only mutational class that was significantly overrepresented in lineages with reduced fitness was the loss of the plasmid, though nonsense mutations, missense mutations, and coding insertion-deletions were also overrepresented in MA lineages whose fitness had significantly declined. Although the overall distribution of fitness effects was similar between the three environments, the magnitude and even the sign of the fitness of a number of lineages changed with the environment, demonstrating that the fitness of some genotypes was environmentally dependent. These results present an unprecedented picture of the fitness effects of spontaneous mutations in a bacterium with multiple chromosomes and provide greater quantitative support for the theory that the vast majority of spontaneous mutations are neutral or deleterious.

Variation of Burkholderia cenocepacia virulence potential during cystic fibrosis chronic lung infection

During long-term lung infection in cystic fibrosis (CF) patients, Burkholderia cenocepacia faces multiple selective pressures in this highly stressful and fluctuating environment. As a consequence, the initial infecting strain undergoes genetic changes that result in the diversification of genotypes and phenotypes. Whether this clonal expansion influences the pathogenic potential is unclear. The virulence potential of 39 sequential B. cenocepacia (recA lineage IIIA) isolates, corresponding to 3 different clones retrieved from 3 chronically infected CF patients was compared in this study using the non-mammalian infection hosts Galleria mellonella and Caenorhabditis elegans. The isolates used in this retrospective study were picked randomly from selective agar plates as part of a CF Center routine, from the onset of infection until patients' death after 3.5 and 7.5 y or the more recent isolation date after 12.5 y of chronic infection. The infection models proved useful to assess virulence potential diversification, but for some isolates the relative values diverged in C. elegans and G. mellonella. Results also reinforce the concept of the occurrence of clonal diversification and co-existence of multiple phenotypes within the CF lungs, also with respect to pathogenicity. No clear trend of decrease (or increase) of the virulence potential throughout long-term infection was found but there is an apparent tendency for a clone/patient-dependent decrease of virulence when the G. mellonella model was used. The sole avirulent variant in both infection hosts was found to lack the small third replicon previously associated to virulence. Although possible, the in vivo loss of this nonessential megaplasmid was found to be a rare event (1 among a total of 64 isolates examined).

Candidate Essential Genes in Burkholderia cenocepacia J2315 Identified by Genome-Wide TraDIS

Burkholderia cenocepacia infection often leads to fatal cepacia syndrome in cystic fibrosis patients. However, antibiotic therapy rarely results in complete eradication of the pathogen due to its intrinsic resistance to many clinically available antibiotics. Recent attention has turned to the identification of essential genes as the proteins encoded by these genes may serve as potential targets for development of novel antimicrobials. In this study, we utilized TraDIS (Transposon Directed Insertion-site Sequencing) as a genome-wide screening tool to facilitate the identification of B. cenocepacia genes essential for its growth and viability. A transposon mutant pool consisting of approximately 500,000 mutants was successfully constructed, with more than 400,000 unique transposon insertion sites identified by computational analysis of TraDIS datasets. The saturated library allowed for the identification of 383 genes that were predicted to be essential in B. cenocepacia. We extended the application of TraDIS to identify conditionally essential genes required for in vitro growth and revealed an additional repertoire of 439 genes to be crucial for B. cenocepacia growth under nutrient-depleted conditions. The library of B. cenocepacia mutants can subsequently be subjected to various biologically related conditions to facilitate the discovery of genes involved in niche adaptation as well as pathogenicity and virulence.

Orderly Replication and Segregation of the Four Replicons of Burkholderia cenocepacia J2315

Bacterial genomes typically consist of a single chromosome and, optionally, one or more plasmids. But whole-genome sequencing reveals about ten per-cent of them to be multipartite, with additional replicons which by size and indispensability are considered secondary chromosomes. This raises the questions of how their replication and partition is managed without compromising genome stability and of how such genomes arose. Vibrio cholerae, with a 1 Mb replicon in addition to its 3 Mb chromosome, is the only species for which maintenance of a multipartite genome has been investigated. In this study we have explored the more complex genome of Burkholderia cenocepacia (strain J2315). It comprises an extra replicon (c2) of 3.21 Mb, comparable in size to the3.87Mb main chromosome (c1), another extra replicon(c3) of 0.87 Mb and a plasmid of 0.09 Mb. The replication origin of c1 is typically chromosomal and those of c2 and c3 are plasmid-like; all are replicated bidirectionally. Fluorescence microscopy of tagged origins indicates that all initiate replication at mid-cell and segregate towards the cell quarter positions sequentially, c1-c2-p1/c3. c2 segregation is as well-phased with the cell cycle as c1, implying that this plasmid-like origin has become subject to regulation not typical of plasmids; in contrast, c3 segregates more randomly through the cycle. Disruption of individual Par systems by deletion of parAB or by addition of parS sites showed each Par system to govern the positioning of its own replicon only. Inactivation of c1, c2 and c3 Par systems not only reduced growth rate, generated anucleate cells and compromised viability but influenced processes beyond replicon partition, notably regulation of replication, chromosome condensation and cell size determination. In particular, the absence of the c1 ParA protein altered replication of all three chromosomes, suggesting that the partition system of the main chromosome is a major participant in the choreography of the cell cycle.

Unlike higher organisms, bacteria typically carry their genetic information on a single chromosome. But in a few bacterial families the genome includes one to three additional chromosome-like DNA molecules. Because these families are rich in pathogenic and environmentally versatile species, it is important to understand how their split genomes evolved and how their maintenance is managed without confusion. We find that mitotic segregation (partition) of all three chromosomes of the cystic fibrosis type strain, Burkholderia cenocepacia J2315, proceeds from mid-cell to cell quarter positions, but that it occurs in a sequential manner, from largest chromosome to smallest. Positioning of each chromosome is specified solely by its own partition proteins. Nevertheless, the partition system of the largest chromosome appears also to play a global role in the cell cycle, by modulating the timing of initiation of replication. In addition, disrupting the partition systems of all three chromosomes induced specific cell abnormalities. Hence, although such bacteria are governed mainly by the largest, housekeeping chromosome, all the Par systems have insinuated themselves into cell cycle regulation to become indispensable for normal growth. Exploration of the underlying mechanisms should allow us to understand their full importance to bacterial life.