The Essential Genome of Burkholderia cenocepacia H111

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.

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.

The Small RNA NcS25 Regulates Biological Amine-Transporting Outer Membrane Porin BCAL3473 in Burkholderia cenocepacia

Regulation of porin expression in bacteria is complex and often involves small-RNA regulators. Several small-RNA regulators have been described for Burkholderia cenocepacia, and this study aimed to characterize the biological role of the conserved small RNA NcS25 and its cognate target, outer membrane protein BCAL3473. The B. cenocepacia genome carries a large number of genes encoding porins with yet-uncharacterized functions. Expression of the porin BCAL3473 is strongly repressed by NcS25 and activated by other factors, such as a LysR-type regulator and nitrogen-depleted growth conditions. The porin is involved in transport of arginine, tyrosine, tyramine, and putrescine across the outer membrane. Porin BCAL3473, with NcS25 as a major regulator, plays an important role in the nitrogen metabolism of B. cenocepacia. IMPORTANCE Burkholderia cenocepacia is a Gram-negative bacterium which causes infections in immunocompromised individuals and in people with cystic fibrosis. A low outer membrane permeability is one of the factors giving it a high level of innate resistance to antibiotics. Porins provide selective permeability for nutrients, and antibiotics can also traverse the outer membrane by this means. Knowing the properties and specificities of porin channels is therefore important for understanding resistance mechanisms and for developing new antibiotics and could help in overcoming permeability issues in antibiotic treatment.

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.

Bacterial genome reductions: Tools, applications, and challenges

Bacterial cells are widely used to produce value-added products due to their versatility, ease of manipulation, and the abundance of genome engineering tools. However, the efficiency of producing these desired biomolecules is often hindered by the cells’ own metabolism, genetic instability, and the toxicity of the product. To overcome these challenges, genome reductions have been performed, making strains with the potential of serving as chassis for downstream applications. Here we review the current technologies that enable the design and construction of such reduced-genome bacteria as well as the challenges that limit their assembly and applicability. While genomic reductions have shown improvement of many cellular characteristics, a major challenge still exists in constructing these cells efficiently and rapidly. Computational tools have been created in attempts at minimizing the time needed to design these organisms, but gaps still exist in modelling these reductions in silico. Genomic reductions are a promising avenue for improving the production of value-added products, constructing chassis cells, and for uncovering cellular function but are currently limited by their time-consuming construction methods. With improvements to and the creation of novel genome editing tools and in silico models, these approaches could be combined to expedite this process and create more streamlined and efficient cell factories.

Phenotype-Guided Comparative Genomics Identifies the Complete Transport Pathway of the Antimicrobial Lasso Peptide Ubonodin in Burkholderia

New antibiotics are needed as bacterial infections continue to be a leading cause of death, but efforts to develop compounds with promising antibacterial activity are hindered by a poor understanding of─and limited strategies for elucidating─their modes of action. We recently discovered a novel lasso peptide, ubonodin, that is active against opportunistic human lung pathogens from the Burkholderia cepacia complex (Bcc). Ubonodin inhibits RNA polymerase, but only select strains were susceptible, indicating that having a conserved cellular target does not guarantee activity. Given the cytoplasmic target, we hypothesized that cellular uptake of ubonodin determines susceptibility. Although Bcc strains harbor numerous nutrient uptake systems, these organisms lack close homologues of the single known lasso peptide membrane receptor, FhuA. Thus, a straightforward homology-driven approach failed to uncover the identity of the ubonodin transporter(s). Here, we used phenotype-guided comparative genomics to identify genes uniquely associated with ubonodin-susceptible Bcc strains, leading to the identification of PupB as the ubonodin outer membrane (OM) receptor in Burkholderia. The loss of PupB renders B. cepacia resistant to ubonodin, whereas expressing PupB sensitizes a resistant strain. We also examine how a conserved iron-regulated transcriptional pathway controls PupB to further tune ubonodin susceptibility. PupB is only the second lasso peptide OM receptor to be uncovered and the first outside of enterobacteria. Finally, we elucidate the full transport pathway for ubonodin by identifying its inner membrane receptor YddA in Burkholderia. Our work provides a complete picture of the mode of action of ubonodin and establishes a general framework for deciphering the transport pathways of other natural products with cytoplasmic targets.

A genome-scale metabolic model of Cupriavidus necator H16 integrated with TraDIS and transcriptomic data reveals metabolic insights for biotechnological applications

Exploiting biological processes to recycle renewable carbon into high value platform chemicals provides a sustainable and greener alternative to current reliance on petrochemicals. In this regard Cupriavidus necator H16 represents a particularly promising microbial chassis due to its ability to grow on a wide range of low-cost feedstocks, including the waste gas carbon dioxide, whilst also naturally producing large quantities of polyhydroxybutyrate (PHB) during nutrient-limited conditions. Understanding the complex metabolic behaviour of this bacterium is a prerequisite for the design of successful engineering strategies for optimising product yields. We present a genome-scale metabolic model (GSM) of C. necator H16 (denoted iCN1361), which is directly constructed from the BioCyc database to improve the readability and reusability of the model. After the initial automated construction, we have performed extensive curation and both theoretical and experimental validation. By carrying out a genome-wide essentiality screening using a Transposon-directed Insertion site Sequencing (TraDIS) approach, we showed that the model could predict gene knockout phenotypes with a high level of accuracy. Importantly, we indicate how experimental and computational predictions can be used to improve model structure and, thus, model accuracy as well as to evaluate potential false positives identified in the experiments. Finally, by integrating transcriptomics data with iCN1361 we create a condition-specific model, which, importantly, better reflects PHB production in C. necator H16. Observed changes in the omics data and in-silico-estimated alterations in fluxes were then used to predict the regulatory control of key cellular processes. The results presented demonstrate that iCN1361 is a valuable tool for unravelling the system-level metabolic behaviour of C. necator H16 and can provide useful insights for designing metabolic engineering strategies.

Genome-scale metabolic models (GSMs) provide a tool for unravelling the complex metabolic behaviour of bacteria and how they adapt to changing environments and genetic perturbations, and thus offer invaluable insights for biotechnology applications. For a GSM to be used efficiently for strain development purposes, however, the model must be easily readable and reusable by other researchers, whilst being able to predict metabolic behaviour with a high level of accuracy. In this work, we developed a GSM for Cupriavidus necator H16 that is linked to the BioCyc database, which provides an efficient way of application, model update, integration of experimental data and network visualisation for other researchers. Using our model, we demonstrate how integrating experimental observations, including Transposon-directed Insertion site Sequencing (TraDIS) and omics data, can be used to compensate for the lack of regulatory, kinetic and thermodynamic information in GSMs, and thus improve model accuracy. Importantly, we found that TraDIS in vivo screening and GSM analysis are complementary approaches, which can be used in combination to provide reliable gene essentiality predictions. Overall, our results offer an informed strategy for the deliberate manipulation of C. necator H16 metabolic capabilities, towards its industrial application to convert greenhouse gases into biochemicals and biofuels.

Gradients in gene essentiality reshape antibacterial research

Essential genes encode the processes that are necessary for life. Until recently, commonly applied binary classifications left no space between essential and non-essential genes. In this review, we frame bacterial gene essentiality in the context of genetic networks. We explore how the quantitative properties of gene essentiality are influenced by the nature of the encoded process, environmental conditions and genetic background, including a strain's distinct evolutionary history. The covered topics have important consequences for antibacterials, which inhibit essential processes. We argue that the quantitative properties of essentiality can thus be used to prioritize antibacterial cellular targets and desired spectrum of activity in specific infection settings. We summarize our points with a case study on the core essential genome of the cystic fibrosis pathobiome and highlight avenues for targeted antibacterial development.

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.

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.

Bacterial surface properties influence the activity of the TAT-RasGAP317-326 antimicrobial peptide

Antibiotic resistance is an increasing threat for public health, underscoring the need for new antibacterial agents. Antimicrobial peptides (AMPs) represent an alternative to classical antibiotics. TAT-RasGAP_317-326 is a recently described AMP effective against a broad range of bacteria, but little is known about the conditions that may influence its activity. Using RNA-sequencing and screening of mutant libraries, we show that Escherichia coli and Pseudomonas aeruginosa respond to TAT-RasGAP_317-326 by regulating metabolic and stress response pathways, possibly implicating two-component systems. Our results also indicate that bacterial surface properties, in particular integrity of the lipopolysaccharide layer, influence peptide binding and entry. Finally, we found differences between bacterial species with respect to their rate of resistance emergence against this peptide. Our findings provide the basis for future investigation on the mode of action of TAT-RasGAP_317-326, which may help developing antimicrobial treatments based on this peptide.

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.

Identification of genes required for gold and silver tolerance in Burkholderia cenocepacia H111 by transposon sequencing

Members of the genus Burkholderia show remarkable abilities to adapt to a wide range of environmental conditions and is frequently isolated from soils contaminated with heavy metals. In this study, we used a transposon sequencing approach to identify 138 and 164 genes that provide a benefit for growth of the opportunistic pathogen Burkholderia cenocepacia H111 in the presence of silver and gold ions respectively. The data suggest that arginine metabolism and citrate biosynthesis are important for silver tolerance, while components of an ABC transporter (BCAL0307-BCAL0308) and de novo cysteine biosynthesis are required for tolerance to gold ions. We show that determinants that affect tolerance to both metal ions include the two-component systems BCAL0497/99 and BCAL2830/31 and genes that are involved in maintaining the integrity of the cell envelope, suggesting that membrane proteins represent important targets of silver and gold ions. Furthermore, we show that that the P-type ATPase CadA (BCAL0055), which confers tolerance to cadmium contributes to silver but not gold tolerance. Our results may be useful for improving the antibacterial effect of silver and gold ions to combat drug-resistant pathogens.

Identification of Genes Required for Resistance to Peptidomimetic Antibiotics by Transposon Sequencing

Pseudomonas aeruginosa is an opportunistic human pathogen and a leading cause of nosocomial infections. Due to its high intrinsic and adaptive resistance to antibiotics, infections caused by this organism are difficult to treat and new therapeutic options are urgently needed. Novel peptidomimetic antibiotics that target outer membrane (OM) proteins have shown great promise for the treatment of P. aeruginosa infections. Here, we have performed genome-wide mutant fitness profiling using transposon sequencing (Tn-Seq) to identify resistance determinants against the recently described peptidomimetics L27-11, compounds 3 and 4, as well as polymyxin B2 (PMB) and colistin (COL). We identified a set of 13 core genes that affected resistance to all tested antibiotics, many of which encode enzymes involved in the modification of the lipopolysaccharide (LPS) or control their expression. We also identified fitness determinants that are specific for antibiotics with similar structures that may indicate differences in their modes of action. These results provide new insights into resistance mechanisms against these peptide antibiotics, which will be important for future clinical development and efforts to further improve their potency.

Copper resistance genes of Burkholderia cenocepacia H111 identified by transposon sequencing

Copper is an essential element but in excess is highly toxic and therefore cytoplasmic levels must be tightly controlled. Member of the genus Burkholderia are highly resistant to various heavy metals and are often isolated from acidic soils where copper bioavailability is high. In this study, we employed transposon sequencing (Tn-Seq) to identify copper resistance genes in Burkholderia cenocepacia H111. We identified a copper efflux system that shares similarities with the plasmid-based copper detoxification systems found in Escherichia coli and Pseudomonas syringae. We also found that several of the identified resistance determinants are involved in maintaining the integrity of the cell envelope, suggesting that proteins located in the outer membrane and periplasmic space are particularly sensitive to copper stress. Given that several of the resistance genes are required for the repair and turnover of misfolded proteins, we suggest that copper toxicity is caused by protein damage rather than by oxidative stress.

Horizontal Gene Transfer to a Defensive Symbiont with a Reduced Genome in a Multipartite Beetle Microbiome

Symbiotic mutualisms of bacteria and animals are ubiquitous in nature, running a continuum from facultative to obligate from the perspectives of both partners. The loss of functions required for living independently but not within a host gives rise to reduced genomes in many symbionts. Although the phenomenon of genome reduction can be explained by existing evolutionary models, the initiation of the process is not well understood. Here, we describe the microbiome associated with the eggs of the beetle Lagria villosa, consisting of multiple bacterial symbionts related to Burkholderia gladioli, including a reduced-genome symbiont thought to be the exclusive producer of the defensive compound lagriamide. We show that the putative lagriamide-producing symbiont is the only member of the microbiome undergoing genome reduction and that it has already lost the majority of its primary metabolism and DNA repair pathways. The key step preceding genome reduction in the symbiont was likely the horizontal acquisition of the putative lagriamide lga biosynthetic gene cluster. Unexpectedly, we uncovered evidence of additional horizontal transfers to the symbiont's genome while genome reduction was occurring and despite a current lack of genes needed for homologous recombination. These gene gains may have given the genome-reduced symbiont a selective advantage in the microbiome, especially given the maintenance of the large lga gene cluster despite ongoing genome reduction.IMPORTANCE Associations between microorganisms and an animal, plant, or fungal host can result in increased dependence over time. This process is due partly to the bacterium not needing to produce nutrients that the host provides, leading to loss of genes that it would need to live independently and to a consequent reduction in genome size. It is often thought that genome reduction is aided by genetic isolation-bacteria that live in monocultures in special host organs, or inside host cells, have less access to other bacterial species from which they can obtain genes. Here, we describe exposure of a genome-reduced beetle symbiont to a community of related bacteria with nonreduced genomes. We show that the symbiont has acquired genes from other bacteria despite going through genome reduction, suggesting that isolation has not yet played a major role in this case of genome reduction, with horizontal gene gains still offering a potential route for adaptation.

In vitro exploration of the Xanthomonas hortorum pv. vitians genome using transposon insertion sequencing and comparative genomics to discriminate between core and contextual essential genes

The essential genome of a bacterium encompasses core genes associated with basic cellular processes and conditionally essential genes dependent upon environmental conditions or the genetic context. Comprehensive knowledge of those gene sets allows for a better understanding of fundamental bacterial biology and offers new perspectives for antimicrobial drug research against detrimental bacteria such as pathogens. We investigated the essential genome of Xanthomonas hortorum pv. vitians, a gammaproteobacterial plant pathogen of lettuce (Lactuca sativa L.) which belongs to the plant-pathogen reservoir genus Xanthomonas and is affiliated to the family Xanthomonadaceae. No practical means of disease control or prevention against this pathogen is currently available, and its molecular biology is virtually unknown. To reach a comprehensive overview of the essential genome of X. hortorum pv. vitians LM16734, we developed a mixed approach combining high-quality full genome sequencing, saturated transposon insertion sequencing (Tn-Seq) in optimal growth conditions, and coupled computational analyses such as comparative genomics, synteny assessment and phylogenomics. Among the 370 essential loci identified by Tn-Seq, a majority was bound to critical cell processes conserved across bacteria. The remaining genes were either related to specific ecological features of Xanthomonas or Xanthomonadaceae species, or acquired through horizontal gene transfer of mobile genetic elements and associated with ancestral parasitic gene behaviour and bacterial defence systems. Our study sheds new light on our usual concepts about gene essentiality and is pioneering in the molecular and genomic study of X. hortorum pv. vitians.

ePath: an online database towards comprehensive essential gene annotation for prokaryotes

Experimental techniques for identification of essential genes (EGs) in prokaryotes are usually expensive, time-consuming and sometimes unrealistic. Emerging in silico methods provide alternative methods for EG prediction, but often possess limitations including heavy computational requirements and lack of biological explanation. Here we propose a new computational algorithm for EG prediction in prokaryotes with an online database (ePath) for quick access to the EG prediction results of over 4,000 prokaryotes ( https://www.pubapps.vcu.edu/epath/ ). In ePath, gene essentiality is linked to biological functions annotated by KEGG Ortholog (KO). Two new scoring systems, namely, E_score and P_score, are proposed for each KO as the EG evaluation criteria. E_score represents appearance and essentiality of a given KO in existing experimental results of gene essentiality, while P_score denotes gene essentiality based on the principle that a gene is essential if it plays a role in genetic information processing, cell envelope maintenance or energy production. The new EG prediction algorithm shows prediction accuracy ranging from 75% to 91% based on validation from five new experimental studies on EG identification. Our overall goal with ePath is to provide a comprehensive and reliable reference for gene essentiality annotation, facilitating the study of those prokaryotes without experimentally derived gene essentiality information.

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.

The role of small proteins in Burkholderia cenocepacia J2315 biofilm formation, persistence and intracellular growth

Burkholderia cenocepacia infections are difficult to treat due to resistance, biofilm formation and persistence. B. cenocepacia strain J2315 has a large multi-replicon genome (8.06 Mb) and the function of a large fraction of (conserved) hypothetical genes remains elusive. The goal of the present study is to elucidate the role of small proteins in B. cenocepacia, focusing on genes smaller than 300 base pairs of which the function is unknown. Almost 10% (572) of the B. cenocepacia J2315 genes are smaller than 300 base pairs and more than half of these are annotated as coding for hypothetical proteins. For 234 of them no similarity could be found with non-hypothetical genes in other bacteria using BLAST. Using available RNA sequencing data obtained from biofilms, a list of 27 highly expressed B. cenocepacia J2315 genes coding for small proteins was compiled. For nine of them expression in biofilms was also confirmed using LC-MS based proteomics and/or expression was confirmed using eGFP translational fusions. Overexpression of two of these genes negatively impacted growth, whereas for four others overexpression led to an increase in biofilm biomass. Overexpression did not have an influence on the MIC for tobramycin, ciprofloxacin or meropenem but for five small protein encoding genes, overexpression had an effect on the number of persister cells in biofilms. While there were no significant differences in adherence to and invasion of A549 epithelial cells between the overexpression mutants and the WT, significant differences were observed in intracellular growth/survival. Finally, the small protein BCAM0271 was identified as an antitoxin belonging to a toxin-antitoxin module. The toxin was found to encode a tRNA acetylase that inhibits translation. In conclusion, our results confirm that small proteins are present in the genome of B. cenocepacia J2315 and indicate that they are involved in various biological processes, including biofilm formation, persistence and intracellular growth.

Comparative analysis of the Burkholderia cenocepacia K56-2 essential genome reveals cell envelope functions that are uniquely required for survival in species of the genus Burkholderia

Burkholderia cenocepacia K56-2 belongs to the Burkholderia cepacia complex, a group of Gram-negative opportunistic pathogens that have large and dynamic genomes. In this work, we identified the essential genome of B. cenocepacia K56-2 using high-density transposon mutagenesis and insertion site sequencing (Tn-seq circle). We constructed a library of one million transposon mutants and identified the transposon insertions at an average of one insertion per 27 bp. The probability of gene essentiality was determined by comparing of the insertion density per gene with the variance of neutral datasets generated by Monte Carlo simulations. Five hundred and eight genes were not significantly disrupted, suggesting that these genes are essential for survival in rich, undefined medium. Comparison of the B. cenocepacia K56-2 essential genome with that of the closely related B. cenocepacia J2315 revealed partial overlapping, suggesting that some essential genes are strain-specific. Furthermore, 158 essential genes were conserved in B. cenocepacia and two species belonging to the Burkholderia pseudomallei complex, B. pseudomallei K96243 and Burkholderia thailandensis E264. Porins, including OpcC, a lysophospholipid transporter, LplT, and a protein involved in the modification of lipid A with aminoarabinose were found to be essential in Burkholderia genomes but not in other bacterial essential genomes identified so far. Our results highlight the existence of cell envelope processes that are uniquely essential in species of the genus Burkholderia for which the essential genomes have been identified by Tn-seq.