High confidence prediction of essential genes in Burkholderia cenocepacia

Identification of discriminant features from stationary pattern of nucleotide bases and their application to essential gene classification

Introduction: Essential genes are essential for the survival of various species. These genes are a family linked to critical cellular activities for species survival. These genes are coded for proteins that regulate central metabolism, gene translation, deoxyribonucleic acid replication, and fundamental cellular structure and facilitate intracellular and extracellular transport. Essential genes preserve crucial genomics information that may hold the key to a detailed knowledge of life and evolution. Essential gene studies have long been regarded as a vital topic in computational biology due to their relevance. An essential gene is composed of adenine, guanine, cytosine, and thymine and its various combinations. Methods: This paper presents a novel method of extracting information on the stationary patterns of nucleotides such as adenine, guanine, cytosine, and thymine in each gene. For this purpose, some co-occurrence matrices are derived that provide the statistical distribution of stationary patterns of nucleotides in the genes, which is helpful in establishing the relationship between the nucleotides. For extracting discriminant features from each co-occurrence matrix, energy, entropy, homogeneity, contrast, and dissimilarity features are computed, which are extracted from all co-occurrence matrices and then concatenated to form a feature vector representing each essential gene. Finally, supervised machine learning algorithms are applied for essential gene classification based on the extracted fixed-dimensional feature vectors. Results: For comparison, some existing state-of-the-art feature representation techniques such as Shannon entropy (SE), Hurst exponent (HE), fractal dimension (FD), and their combinations have been utilized. Discussion: An extensive experiment has been performed for classifying the essential genes of five species that show the robustness and effectiveness of the proposed methodology.

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.

An overview and future prospects of recombinant protein production in Bacillus subtilis

Bacillus subtilis is a well-characterized Gram-positive bacterium and a valuable host for recombinant protein production because of its efficient secretion ability, high yield, and non-toxicity. Here, we comprehensively review the recent studies on recombinant protein production in B. subtilis to update and supplement other previous reviews. We have focused on several aspects, including optimization of B. subtilis strains, enhancement and regulation of expression, improvement of secretion level, surface display of proteins, and fermentation optimization. Among them, optimization of B. subtilis strains mainly involves undirected chemical/physical mutagenesis and selection and genetic manipulation; enhancement and regulation of expression comprises autonomous plasmid and integrated expression, promoter regulation and engineering, and fine-tuning gene expression based on proteases and molecular chaperones; improvement of secretion level predominantly involves secretion pathway and signal peptide screening and optimization; surface display of proteins includes surface display of proteins on spores or vegetative cells; and fermentation optimization incorporates medium optimization, process condition optimization, and feeding strategy optimization. Furthermore, we propose some novel methods and future challenges for recombinant protein production in B. subtilis.Key points• A comprehensive review on recombinant protein production in Bacillus subtilis.• Novel techniques facilitate recombinant protein expression and secretion.• Surface display of proteins has significant potential for different applications.

Harnessing model organism genomics to underpin the machine learning-based prediction of essential genes in eukaryotes - Biotechnological implications

The availability of high-quality genomes and advances in functional genomics have enabled large-scale studies of essential genes in model eukaryotes, including the 'elegant worm' (Caenorhabditis elegans; Nematoda) and the 'vinegar fly' (Drosophila melanogaster; Arthropoda). However, this is not the case for other, much less-studied organisms, such as socioeconomically important parasites, for which functional genomic platforms usually do not exist. Thus, there is a need to develop innovative techniques or approaches for the prediction, identification and investigation of essential genes. A key approach that could enable the prediction of such genes is machine learning (ML). Here, we undertake an historical review of experimental and computational approaches employed for the characterisation of essential genes in eukaryotes, with a particular focus on model ecdysozoans (C. elegans and D. melanogaster), and discuss the possible applicability of ML-approaches to organisms such as socioeconomically important parasites. We highlight some recent results showing that high-performance ML, combined with feature engineering, allows a reliable prediction of essential genes from extensive, publicly available 'omic data sets, with major potential to prioritise such genes (with statistical confidence) for subsequent functional genomic validation. These findings could 'open the door' to fundamental and applied research areas. Evidence of some commonality in the essential gene-complement between these two organisms indicates that an ML-engineering approach could find broader applicability to ecdysozoans such as parasitic nematodes or arthropods, provided that suitably large and informative data sets become/are available for proper feature engineering, and for the robust training and validation of algorithms. This area warrants detailed exploration to, for example, facilitate the identification and characterisation of essential molecules as novel targets for drugs and vaccines against parasitic diseases. This focus is particularly important, given the substantial impact that such diseases have worldwide, and the current challenges associated with their prevention and control and with drug resistance in parasite populations.

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.

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.

Identifying Bacterial Essential Genes Based on a Feature-Integrated Method

Essential genes are those genes of an organism that are considered to be crucial for its survival. Identification of essential genes is therefore of great significance to advance our understanding of the principles of cellular life. We have developed a novel computational method, which can effectively predict bacterial essential genes by extracting and integrating homologous features, protein domain feature, gene intrinsic features, and network topological features. By performing the principal component regression (PCR) analysis for Escherichia coli MG1655, we established a classification model with the average area under curve (AUC) value of 0.992 in ten times 5-fold cross-validation tests. Furthermore, when employing this new model to a distantly related organism-Streptococcus pneumoniae TIGR4, we still got a reliable AUC value of 0.788. These results indicate that our feature-integrated approach could have practical applications in accurately investigating essential genes from broad bacterial species, and also provide helpful guidelines for the minimal cell.

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.

Comprehensive review of the identification of essential genes using computational methods: focusing on feature implementation and assessment

Essential genes have attracted increasing attention in recent years due to the important functions of these genes in organisms. Among the methods used to identify the essential genes, accurate and efficient computational methods can make up for the deficiencies of expensive and time-consuming experimental technologies. In this review, we have collected researches on essential gene predictions in prokaryotes and eukaryotes and summarized the five predominant types of features used in these studies. The five types of features include evolutionary conservation, domain information, network topology, sequence component and expression level. We have described how to implement the useful forms of these features and evaluated their performance based on the data of Escherichia coli MG1655, Bacillus subtilis 168 and human. The prerequisite and applicable range of these features is described. In addition, we have investigated the techniques used to weight features in various models. To facilitate researchers in the field, two available online tools, which are accessible for free and can be directly used to predict gene essentiality in prokaryotes and humans, were referred. This article provides a simple guide for the identification of essential genes in prokaryotes and eukaryotes.

Variability of Bacterial Essential Genes Among Closely Related Bacteria: The Case of Escherichia coli

The definition of bacterial essential genes has been widely pursued using different approaches. Their study has impacted several fields of research such as synthetic biology, the construction of bacteria with minimal chromosomes, the search for new antibiotic targets, or the design of strains with biotechnological applications. Bacterial genomes are mosaics that only share a small subset of gene-sequences (core genome) even among members of the same species. It has been reported that the presence of essential genes is highly variable between closely related bacteria and even among members of the same species, due to the phenomenon known as “non-orthologous gene displacement” that refers to the coding for an essential function by genes with no sequence homology due to horizontal gene transfer (HGT). The existence of dormant forms among bacteria and the high incidence of HGT have been proposed to be driving forces of bacterial evolution, and they might have a role in the low level of conservation of essential genes among related bacteria by non-orthologous gene displacement, but this correlation has not been recognized. The aim of this mini-review is to give a brief overview of the approaches that have been taken to define and study essential genes, and the implications of non-orthologous gene displacement in bacterial evolution, focusing mainly in the case of Escherichia coli. To this end, we reviewed the available literature, and we searched for the presence of the essential genes defined by mutagenesis in the genomes of the 63 best-sequenced E. coli genomes that are available in NCBI database. We could not document specific cases of non-orthologous gene displacement among the E. coli strains analyzed, but we found that the quality of the genome-sequences in the database is not enough to make accurate predictions about the conservation of essential-genes among members of this bacterial species.

Genetic Determinants Associated With in Vivo Survival of Burkholderia cenocepacia in the Caenorhabditis elegans Model

A Burkholderia cenocepacia infection usually leads to reduced survival and fatal cepacia syndrome in cystic fibrosis patients. The identification of B. cenocepacia essential genes for in vivo survival is key to designing new anti-infectives therapies. We used the Transposon-Directed Insertion Sequencing (TraDIS) approach to identify genes required for B. cenocepacia survival in the model infection host, Caenorhabditis elegans. A B. cenocepacia J2315 transposon pool of ∼500,000 mutants was used to infect C. elegans. We identified 178 genes as crucial for B. cenocepacia survival in the infected nematode. The majority of these genes code for proteins of unknown function, many of which are encoded by the genomic island BcenGI13, while other gene products are involved in nutrient acquisition, general stress responses and LPS O-antigen biosynthesis. Deletion of the glycosyltransferase gene wbxB and a histone-like nucleoid structuring (H-NS) protein-encoding gene (BCAL0154) reduced bacterial accumulation and attenuated virulence in C. elegans. Further analysis using quantitative RT-PCR indicated that BCAL0154 modulates B. cenocepacia pathogenesis via transcriptional regulation of motility-associated genes including fliC, fliG, flhD, and cheB1. This screen has successfully identified genes required for B. cenocepacia survival within the host-associated environment, many of which are potential targets for developing new antimicrobials.

Targeting the Nonmevalonate Pathway in Burkholderia cenocepacia Increases Susceptibility to Certain β-Lactam Antibiotics

The nonmevalonate pathway is the sole pathway for isoprenoid biosynthesis in Burkholderia cenocepacia and is possibly a novel target for the development of antibacterial chemotherapy. The goals of the present study were to evaluate the essentiality of dxr, the second gene of the nonmevalonate pathway, in B. cenocepacia and to determine whether interfering with the nonmevalonate pathway increases susceptibility toward antibiotics. To this end, a rhamnose-inducible conditional dxr knockdown mutant of B. cenocepacia strain K56-2 (B. cenocepacia K56-2dxr) was constructed, using a plasmid which enables the delivery of a rhamnose-inducible promoter in the chromosome. Expression of dxr is essential for bacterial growth; the growth defect observed in the dxr mutant could be complemented by expressing dxr in trans under the control of a constitutive promoter, but not by providing 2-C-methyl-d-erythritol-4-phosphate, the reaction product of DXR (1-deoxy-d-xylulose 5-phosphate reductoisomerase). B. cenocepacia K56-2dxr showed markedly increased susceptibility to the β-lactam antibiotics aztreonam, ceftazidime, and cefotaxime, while susceptibility to other antibiotics was not (or was much less) affected; this increased susceptibility could also be complemented by in trans expression of dxr A similarly increased susceptibility was observed when antibiotics were combined with FR900098, a known DXR inhibitor. Our data confirm that the nonmevalonate pathway is essential in B. cenocepacia and suggest that combining potent DXR inhibitors with selected β-lactam antibiotics is a useful strategy to combat B. cenocepacia infections.

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.

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.

Identification and analysis of genomic islands in Burkholderia cenocepacia AU 1054 with emphasis on pathogenicity islands

BACKGROUND

Genomic islands (GIs) are genomic regions that reveal evidence of horizontal DNA transfer. They can code for many functions and may augment a bacterium's adaptation to its host or environment. GIs have been identified in strain J2315 of Burkholderia cenocepacia, whereas in strain AU 1054 there has been no published works on such regions according to our text mining and keyword search in Medline.

RESULTS

In this study, we identified 21 GIs in AU 1054 by combining two computational tools. Feature analyses suggested that the predictions are highly reliable and hence illustrated the advantage of joint predictions by two independent methods. Based on putative virulence factors, four GIs were further identified as pathogenicity islands (PAIs). Through experiments of gene deletion mutants in live bacteria, two putative PAIs were confirmed, and the virulence factors involved were identified as lipA and copR. The importance of the genes lipA (from PAI 1) and copR (from PAI 2) for bacterial invasion and replication indicates that they are required for the invasive properties of B. cenocepacia and may function as virulence determinants for bacterial pathogenesis and host infection.

CONCLUSIONS

This approach of in silico prediction of GIs and subsequent identification of potential virulence factors in the putative island regions with final validation using wet experiments could be used as an effective strategy to rapidly discover novel virulence factors in other bacterial species and strains.

The Escherichia coli rhaSR-PrhaBAD Inducible Promoter System Allows Tightly Controlled Gene Expression over a Wide Range in Pseudomonas aeruginosa

The araC-ParaBAD inducible promoter system is tightly controlled and allows gene expression to be modulated over a wide range in Escherichia coli, which has led to its widespread use in other bacteria. Although anecdotal evidence suggests that araC-ParaBAD is leaky in Pseudomonas aeruginosa, neither a thorough analysis of this inducible promoter system in P. aeruginosa nor a concerted effort to identify alternatives with improved functionality has been reported. Here, we evaluated the functionality of the araC-ParaBAD system in P. aeruginosa Using transcriptional fusions to a lacZ reporter gene, we determined that the noninduced expression from araC-ParaBAD is high and cannot be reduced by carbon catabolite repression as it can in E. coli Modulating translational initiation by altering ribosome-binding site strength reduced the noninduced activity but also decreased the maximal induced activity and narrowed the induction range. Integrating the inducible promoter system into the posttranscriptional regulatory network that controls catabolite repression in P. aeruginosa significantly decreased the noninduced activity and increased the induction range. In addition to these improvements in the functionality of the araC-ParaBAD system, we found that the lacI^q-Ptac and rhaSR-PrhaBAD inducible promoter systems had significantly lower noninduced expression and were inducible over a broader range than araC-ParaBAD We demonstrated that noninduced expression from the araC-ParaBAD system supported the function of genes involved in antibiotic resistance and tryptophan biosynthesis in P. aeruginosa, problems that were avoided with rhaSR-PrhaBAD. rhaSR-PrhaBAD is tightly controlled, allows gene expression over a wide range, and represents a significant improvement over araC-ParaBAD in P. aeruginosa IMPORTANCE: We report the shortcomings of the commonly used Escherichia coli araC-ParaBAD inducible promoter system in Pseudomonas aeruginosa, successfully reengineered it to improve its functionality, and show that the E. coli rhaSR-PrhaBAD system is tightly controlled and allows inducible gene expression over a wide range in P. aeruginosa.

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.

Integrative bacterial artificial chromosomes for DNA integration into the Bacillus subtilis chromosome

Bacillus subtilis is a well-characterized model bacterium frequently used for a number of biotechnology and synthetic biology applications. Novel strategies combining the advantages of B. subtilis with the DNA assembly and editing tools of Escherichia coli are crucial for B. subtilis engineering efforts. We combined Gibson Assembly and λ red recombineering in E. coli with RecA-mediated homologous recombination in B. subtilis for bacterial artificial chromosome-mediated DNA integration into the well-characterized amyE target locus of the B. subtilis chromosome. The engineered integrative bacterial artificial chromosome iBAC(cav) can accept any DNA fragment for integration into B. subtilis chromosome and allows rapid selection of transformants by B. subtilis-specific antibiotic resistance and the yellow fluorescent protein (mVenus) expression. We used the developed iBAC(cav)-mediated system to integrate 10kb DNA fragment from E. coli K12 MG1655 into B. subtilis chromosome. iBAC(cav)-mediated chromosomal integration approach will facilitate rational design of synthetic biology applications in B. subtilis.

Flagellar region 3b supports strong expression of integrated DNA and the highest chromosomal integration efficiency of the Escherichia coli flagellar regions

The Gram-negative bacterium Escherichia coli is routinely used as the chassis for a variety of biotechnology and synthetic biology applications. Identification and analysis of reliable chromosomal integration and expression target loci is crucial for E. coli engineering. Chromosomal loci differ significantly in their ability to support integration and expression of the integrated genetic circuits. In this study, we investigate E. coli K12 MG1655 flagellar regions 2 and 3b. Integration of the genetic circuit into seven and nine highly conserved genes of the flagellar regions 2 (motA, motB, flhD, flhE, cheW, cheY and cheZ) and 3b (fliE, F, G, J, K, L, M, P, R), respectively, showed significant variation in their ability to support chromosomal integration and expression of the integrated genetic circuit. While not reducing the growth of the engineered strains, the integrations into all 16 target sites led to the loss of motility. In addition to high expression, the flagellar region 3b supports the highest efficiency of integration of all E. coli K12 MG1655 flagellar regions and is therefore potentially the most suitable for the integration of synthetic genetic circuits.

Predicting Essential Genes and Proteins Based on Machine Learning and Network Topological Features: A Comprehensive Review

Essential proteins/genes are indispensable to the survival or reproduction of an organism, and the deletion of such essential proteins will result in lethality or infertility. The identification of essential genes is very important not only for understanding the minimal requirements for survival of an organism, but also for finding human disease genes and new drug targets. Experimental methods for identifying essential genes are costly, time-consuming, and laborious. With the accumulation of sequenced genomes data and high-throughput experimental data, many computational methods for identifying essential proteins are proposed, which are useful complements to experimental methods. In this review, we show the state-of-the-art methods for identifying essential genes and proteins based on machine learning and network topological features, point out the progress and limitations of current methods, and discuss the challenges and directions for further research.

High molecular weight DNA assembly in vivo for synthetic biology applications

DNA assembly is the key technology of the emerging interdisciplinary field of synthetic biology. While the assembly of smaller DNA fragments is usually performed in vitro, high molecular weight DNA molecules are assembled in vivo via homologous recombination in the host cell. Escherichia coli, Bacillus subtilis and Saccharomyces cerevisiae are the main hosts used for DNA assembly in vivo. Progress in DNA assembly over the last few years has paved the way for the construction of whole genomes. This review provides an update on recent synthetic biology advances with particular emphasis on high molecular weight DNA assembly in vivo in E. coli, B. subtilis and S. cerevisiae. Special attention is paid to the assembly of whole genomes, such as those of the first synthetic cell, synthetic yeast and minimal genomes.

The genome analysis of Candidatus Burkholderia crenata reveals that secondary metabolism may be a key function of the Ardisia crenata leaf nodule symbiosis

A majority of Ardisia species harbour Burkholderia sp. bacteria within specialized leaf nodules. The bacteria are transmitted hereditarily and have not yet been cultured outside of their host. Because the plants cannot develop beyond the seedling stage without their symbionts, the symbiosis is considered obligatory. We sequenced for the first time the genome of Candidatus Burkholderia crenata (Ca. B. crenata), the leaf nodule symbiont of Ardisia crenata. The genome of Ca. B. crenata is the smallest Burkholderia genome to date. It contains a large amount of insertion sequences and pseudogenes and displays features consistent with reductive genome evolution. The genome does not encode functions commonly associated with plant symbioses such as nitrogen fixation and plant hormone metabolism. However, we identified unique genes with a predicted role in secondary metabolism in the genome of Ca. B. crenata. Specifically, we provide evidence that the bacterial symbionts are responsible for the synthesis of compound FR900359, a cyclic depsipeptide with biomedical properties previously isolated from leaves of A. crenata.

Genome-wide transcription start site profiling in biofilm-grown Burkholderia cenocepacia J2315

BACKGROUND

Burkholderia cenocepacia is a soil-dwelling Gram-negative Betaproteobacterium with an important role as opportunistic pathogen in humans. Infections with B. cenocepacia are very difficult to treat due to their high intrinsic resistance to most antibiotics. Biofilm formation further adds to their antibiotic resistance. B. cenocepacia harbours a large, multi-replicon genome with a high GC-content, the reference genome of strain J2315 includes 7374 annotated genes. This study aims to annotate transcription start sites and identify novel transcripts on a whole genome scale.

METHODS

RNA extracted from B. cenocepacia J2315 biofilms was analysed by differential RNA-sequencing and the resulting dataset compared to data derived from conventional, global RNA-sequencing. Transcription start sites were annotated and further analysed according to their position relative to annotated genes.

RESULTS

Four thousand ten transcription start sites were mapped over the whole B. cenocepacia genome and the primary transcription start site of 2089 genes expressed in B. cenocepacia biofilms were defined. For 64 genes a start codon alternative to the annotated one was proposed. Substantial antisense transcription for 105 genes and two novel protein coding sequences were identified. The distribution of internal transcription start sites can be used to identify genomic islands in B. cenocepacia. A potassium pump strongly induced only under biofilm conditions was found and 15 non-coding small RNAs highly expressed in biofilms were discovered.

CONCLUSIONS

Mapping transcription start sites across the B. cenocepacia genome added relevant information to the J2315 annotation. Genes and novel regulatory RNAs putatively involved in B. cenocepacia biofilm formation were identified. These findings will help in understanding regulation of B. cenocepacia biofilm formation.

Pseudomonas aeruginosa essentials: an update on investigation of essential genes

Pseudomonas aeruginosa is the leading cause of nosocomial infections, particularly in immunocompromised, cancer, burn and cystic fibrosis patients. Development of novel antimicrobials against P. aeruginosa is therefore of the highest importance. Although the first reports on P. aeruginosa essential genes date back to the early 2000s, a number of more sensitive genomic approaches have been used recently to better define essential genes in this organism. These analyses highlight the evolution of the definition of an 'essential' gene from the traditional to the context-dependent. Essential genes, particularly those indispensable under the clinically relevant conditions, are considered to be promising targets of novel antibiotics against P. aeruginosa. This review provides an update on the investigation of P. aeruginosa essential genes. Special focus is on recently identified P. aeruginosa essential genes and their exploitation for the development of antimicrobials.

Draft Genome Sequences of Burkholderia contaminans, a Burkholderia cepacia Complex Species That Is Increasingly Recovered from Cystic Fibrosis Patients

Burkholderia contaminans belongs to the Burkholderia cepacia complex (BCC), a group of bacteria that are ubiquitous in the environment and capable of infecting the immunocompromised and people with cystic fibrosis. We report here draft genome sequences for the B. contaminans type strain LMG 23361 and an Argentinian cystic fibrosis sputum isolate.

Comparative analysis of essential genes in prokaryotic genomic islands

Essential genes are thought to encode proteins that carry out the basic functions to sustain a cellular life, and genomic islands (GIs) usually contain clusters of horizontally transferred genes. It has been assumed that essential genes are not likely to be located in GIs, but systematical analysis of essential genes in GIs has not been explored before. Here, we have analyzed the essential genes in 28 prokaryotes by statistical method and reached a conclusion that essential genes in GIs are significantly fewer than those outside GIs. The function of 362 essential genes found in GIs has been explored further by BLAST against the Virulence Factor Database (VFDB) and the phage/prophage sequence database of PHAge Search Tool (PHAST). Consequently, 64 and 60 eligible essential genes are found to share the sequence similarity with the virulence factors and phage/prophages-related genes, respectively. Meanwhile, we find several toxin-related proteins and repressors encoded by these essential genes in GIs. The comparative analysis of essential genes in genomic islands will not only shed new light on the development of the prediction algorithm of essential genes, but also give a clue to detect the functionality of essential genes in genomic islands.

Identification and validation of novel chromosomal integration and expression loci in Escherichia coli flagellar region 1

Escherichia coli is used as a chassis for a number of Synthetic Biology applications. The lack of suitable chromosomal integration and expression loci is among the main hurdles of the E. coli engineering efforts. We identified and validated chromosomal integration and expression target sites within E. coli K12 MG1655 flagellar region 1. We analyzed five open reading frames of the flagellar region 1, flgA, flgF, flgG, flgI, and flgJ, that are well-conserved among commonly-used E. coli strains, such as MG1655, W3110, DH10B and BL21-DE3. The efficiency of the integration into the E. coli chromosome and the expression of the introduced genetic circuit at the investigated loci varied significantly. The integrations did not have a negative impact on growth; however, they completely abolished motility. From the investigated E. coli K12 MG1655 flagellar region 1, flgA and flgG are the most suitable chromosomal integration and expression loci.

On the road to synthetic life: the minimal cell and genome-scale engineering

Synthetic biology employs rational engineering principles to build biological systems from the libraries of standard, well characterized biological parts. Biological systems designed and built by synthetic biologists fulfill a plethora of useful purposes, ranging from better healthcare and energy production to biomanufacturing. Recent advancements in the synthesis, assembly and "booting-up" of synthetic genomes and in low and high-throughput genome engineering have paved the way for engineering on the genome-wide scale. One of the key goals of genome engineering is the construction of minimal genomes consisting solely of essential genes (genes indispensable for survival of living organisms). Besides serving as a toolbox to understand the universal principles of life, the cell encoded by minimal genome could be used to build a stringently controlled "cell factory" with a desired phenotype. This review provides an update on recent advances in the genome-scale engineering with particular emphasis on the engineering of minimal genomes. Furthermore, it presents an ongoing discussion to the scientific community for better suitability of minimal or robust cells for industrial applications.

Bacillus subtilis and Escherichia coli essential genes and minimal cell factories after one decade of genome engineering

Investigation of essential genes, besides contributing to understanding the fundamental principles of life, has numerous practical applications. Essential genes can be exploited as building blocks of a tightly controlled cell ‘chassis’. Bacillus subtilis and Escherichia coli K-12 are both well-characterized model bacteria used as hosts for a plethora of biotechnological applications. Determination of the essential genes that constitute the B. subtilis and E. coli minimal genomes is therefore of the highest importance. Recent advances have led to the modification of the original B. subtilis and E. coli essential gene sets identified 10 years ago. Furthermore, significant progress has been made in the area of genome minimization of both model bacteria. This review provides an update, with particular emphasis on the current essential gene sets and their comparison with the original gene sets identified 10 years ago. Special attention is focused on the genome reduction analyses in B. subtilis and E. coli and the construction of minimal cell factories for industrial applications.

Escherichia coli flagellar genes as target sites for integration and expression of genetic circuits

E. coli is a model platform for engineering microbes, so genetic circuit design and analysis will be greatly facilitated by simple and effective approaches to introduce genetic constructs into the E. coli chromosome at well-characterised loci. We combined the Red recombinase system of bacteriophage λ and Isothermal Gibson Assembly for rapid integration of novel DNA constructs into the E. coli chromosome. We identified the flagellar region as a promising region for integration and expression of genetic circuits. We characterised integration and expression at four candidate loci, fliD, fliS, fliT, and fliY, of the E. coli flagellar region 3a. The integration efficiency and expression from the four integrations varied considerably. Integration into fliD and fliS significantly decreased motility, while integration into fliT and fliY had only a minor effect on the motility. None of the integrations had negative effects on the growth of the bacteria. Overall, we found that fliT was the most suitable integration site.

The complete genome of Burkholderia phenoliruptrix strain BR3459a, a symbiont of Mimosa flocculosa: highlighting the coexistence of symbiotic and pathogenic genes

BACKGROUND

Burkholderia species play an important ecological role related to xenobiosis, the promotion of plant growth, the biocontrol of agricultural diseases, and symbiotic and non-symbiotic biological nitrogen fixation. Here, we highlight our study as providing the first complete genome of a symbiotic strain of B. phenoliruptrix, BR3459a (=CLA1), which was originally isolated in Brazil from nodules of Mimosa flocculosa and is effective in fixing nitrogen in association with this leguminous species.

RESULTS

Genomic comparisons with other pathogenic and non-pathogenic Burkholderia strains grouped B. phenoliruptrix BR3459a with plant-associated beneficial and environmental species, although it shares a high percentage of its gene repertoire with species of the B. cepacia complex (Bcc) and "pseudomallei" group. The genomic analyses showed that the bce genes involved in exopolysaccharide production are clustered together in the same genomic region, constituting part of the Group III cluster of non-pathogenic bacteria. Regarding environmental stresses, we highlight genes that might be relevant in responses to osmotic, heat, cold and general stresses. Furthermore, a number of particularly interesting genes involved in the machinery of the T1SS, T2SS, T3SS, T4ASS and T6SS secretion systems were identified. The xenobiotic properties of strain BR3459a were also investigated, and some enzymes involved in the degradation of styrene, nitrotoluene, dioxin, chlorocyclohexane, chlorobenzene and caprolactam were identified. The genomic analyses also revealed a large number of antibiotic-related genes, the most important of which were correlated with streptomycin and novobiocin. The symbiotic plasmid showed high sequence identity with the symbiotic plasmid of B. phymatum. Additionally, comparative analysis of 545 housekeeping genes among pathogenic and non-pathogenic Burkholderia species strongly supports the definition of a new genus for the second branch, which would include BR3459a.

CONCLUSIONS

The analyses of B. phenoliruptrix BR3459a showed key property of fixing nitrogen that together with genes for high tolerance to environmental stresses might explain a successful strategy of symbiosis in the tropics. The strain also harbours interesting sets of genes with biotechnological potential. However, the resemblance of certain genes to those of pathogenic Burkholderia raise concerns about large-scale applications in agriculture or for bioremediation.

IFIM: a database of integrated fitness information for microbial genes

Knowledge of an organism’s fitness for survival is important for a complete understanding of microbial genetics and effective drug design. Current essential gene databases provide only binary essentiality data from genome-wide experiments. We therefore developed a new database that Integrates quantitative Fitness Information for Microbial genes (IFIM). The IFIM database currently contains data from 16 experiments and 2186 theoretical predictions. The highly significant correlation between the experiment-derived fitness data and our computational simulations demonstrated that the computer-generated predictions were often as reliable as the experimental data. The data in IFIM can be accessed easily, and the interface allows users to browse through the gene fitness information that it contains. IFIM is the first resource that allows easy access to fitness data of microbial genes. We believe this database will contribute to a better understanding of microbial genetics and will be useful in designing drugs to resist microbial pathogens, especially when experimental data are unavailable.

Database URL:http://cefg.uestc.edu.cn/ifim/ or http://cefg.cn/ifim/

Identification of the flagellin glycosylation system in Burkholderia cenocepacia and the contribution of glycosylated flagellin to evasion of human innate immune responses

Burkholderia cenocepacia is an opportunistic pathogen threatening patients with cystic fibrosis. Flagella are required for biofilm formation, as well as adhesion to and invasion of epithelial cells. Recognition of flagellin via the Toll-like receptor 5 (TLR5) contributes to exacerbate B. cenocepacia-induced lung epithelial inflammatory responses. In this study, we report that B. cenocepacia flagellin is glycosylated on at least 10 different sites with a single sugar, 4,6-dideoxy-4-(3-hydroxybutanoylamino)-D-glucose. We have identified key genes that are required for flagellin glycosylation, including a predicted glycosyltransferase gene that is linked to the flagellin biosynthesis cluster and a putative acetyltransferase gene located within the O-antigen lipopolysaccharide cluster. Another O-antigen cluster gene, rmlB, which is required for flagellin glycan and O-antigen biosynthesis, was essential for bacterial viability, uncovering a novel target against Burkholderia infections. Using glycosylated and nonglycosylated purified flagellin and a cell reporter system to assess TLR5-mediated responses, we also show that the presence of glycan in flagellin significantly impairs the inflammatory response of epithelial cells. We therefore suggest that flagellin glycosylation reduces recognition of flagellin by host TLR5, providing an evasive strategy to infecting bacteria.

Construction of an inducible system for the analysis of essential genes in Yersinia pestis

Yersinia pestis, a Gram negative bacterium, causes bubonic and pneumonic plague. Emerging antibiotic resistance in clinical isolates is driving a need to develop novel antibiotics to treat infection by this transmissible and highly virulent pathogen. Proteins required for viability, so called essential genes, are attractive potential therapeutic targets, however, confirmation of essentiality is problematic. For the first time, we report the development of a system that allows the rapid determination of Y. pestis gene essentiality through mutagenesis and inducible expression of a plasmid borne copy of the target gene. Using this approach, we have confirmed the uridine monophosphate kinase PyrH as an essential protein in Y. pestis. This methodology and the tools we have developed will allow the confirmation of other putative essential genes in this dangerous pathogen, and facilitate the identification of novel targets for antimicrobial development.

Training set selection for the prediction of essential genes

Various computational models have been developed to transfer annotations of gene essentiality between organisms. However, despite the increasing number of microorganisms with well-characterized sets of essential genes, selection of appropriate training sets for predicting the essential genes of poorly-studied or newly sequenced organisms remains challenging. In this study, a machine learning approach was applied reciprocally to predict the essential genes in 21 microorganisms. Results showed that training set selection greatly influenced predictive accuracy. We determined four criteria for training set selection: (1) essential genes in the selected training set should be reliable; (2) the growth conditions in which essential genes are defined should be consistent in training and prediction sets; (3) species used as training set should be closely related to the target organism; and (4) organisms used as training and prediction sets should exhibit similar phenotypes or lifestyles. We then analyzed the performance of an incomplete training set and an integrated training set with multiple organisms. We found that the size of the training set should be at least 10% of the total genes to yield accurate predictions. Additionally, the integrated training sets exhibited remarkable increase in stability and accuracy compared with single sets. Finally, we compared the performance of the integrated training sets with the four criteria and with random selection. The results revealed that a rational selection of training sets based on our criteria yields better performance than random selection. Thus, our results provide empirical guidance on training set selection for the identification of essential genes on a genome-wide scale.

antibacTR: dynamic antibacterial-drug-target ranking integrating comparative genomics, structural analysis and experimental annotation

BACKGROUND

Development of novel antibacterial drugs is both an urgent healthcare necessity and a partially neglected field. The last decades have seen a substantial decrease in the discovery of novel antibiotics, which combined with the recent thrive of multi-drug-resistant pathogens have generated a scenario of general concern. The procedures involved in the discovery and development of novel antibiotics are economically challenging, time consuming and lack any warranty of success. Furthermore, the return-on-investment for an antibacterial drug is usually marginal when compared to other therapeutics, which in part explains the decrease of private investment.

RESULTS

In this work we present antibacTR, a computational pipeline designed to aid researchers in the selection of potential drug targets, one of the initial steps in antibacterial-drug discovery. The approach was designed and implemented as part of two publicly funded initiatives aimed at discovering novel antibacterial targets, mechanisms and drugs for a priority list of Gram-negative pathogens: Acinetobacter baumannii, Escherichia coli, Helicobacter pylori, Pseudomonas aeruginosa and Stenotrophomonas maltophilia. However, at present this list has been extended to cover a total of 74 fully sequenced Gram-negative pathogens. antibacTR is based on sequence comparisons and queries to multiple databases (e.g. gene essentiality, virulence factors) to rank proteins according to their potential as antibacterial targets. The dynamic ranking of potential drug targets can easily be executed, customized and accessed by the user through a web interface which also integrates computational analyses performed in-house and visualizable on-site. These include three-dimensional modeling of protein structures and prediction of active sites among other functionally relevant ligand-binding sites.

CONCLUSIONS

Given its versatility and ease-of-use at integrating both experimental annotation and computational analyses, antibacTR may effectively assist microbiologists, medicinal-chemists and other researchers working in the field of antibacterial drug-discovery. The public web-interface for antibacTR is available at 'http://bioinf.uab.cat/antibactr'.

The third replicon of members of the Burkholderia cepacia Complex, plasmid pC3, plays a role in stress tolerance

The metabolically versatile Burkholderia cepacia complex (Bcc) occupies a variety of niches, including the plant rhizosphere and the cystic fibrosis lung (where it is often fatal to the patient). Bcc members have multipartite genomes, of which the third replicon, pC3 (previously chromosome 3), has been shown to be a nonessential megaplasmid which confers virulence and both antifungal and proteolytic activity on several strains. In this study, pC3 curing was extended to cover strains of 16 of the 17 members of the Bcc, and the phenotypes conferred by pC3 were determined. B. cenocepacia strains H111, MCO-3, and HI2424 were previously cured of pC3; however, this had not proved possible in the epidemic strain K56-2. Here, we investigated the mechanism of this unexpected stability and found that efficient toxin-antitoxin systems are responsible for maintaining pC3 of strain K56-2. Identification of these systems allowed neutralization of the toxins and the subsequent deletion of K56-2pC3. The cured strain was found to exhibit reduced antifungal activity and was attenuated in both the zebrafish and the Caenorhabditis elegans model of infection. We used a PCR screening method to examine the prevalence of pC3 within 110 Bcc isolates and found that this replicon was absent in only four cases, suggesting evolutionary fixation. It is shown that plasmid pC3 increases the resistance of B. cenocepacia H111 to various stresses (oxidative, osmotic, high-temperature, and chlorhexidine-induced stresses), explaining the prevalence of this replicon within the Bcc.

DEG 10, an update of the database of essential genes that includes both protein-coding genes and noncoding genomic elements

The combination of high-density transposon-mediated mutagenesis and high-throughput sequencing has led to significant advancements in research on essential genes, resulting in a dramatic increase in the number of identified prokaryotic essential genes under diverse conditions and a revised essential-gene concept that includes all essential genomic elements, rather than focusing on protein-coding genes only. DEG 10, a new release of the Database of Essential Genes (available at http://www.essentialgene.org), has been developed to accommodate these quantitative and qualitative advancements. In addition to increasing the number of bacterial and archaeal essential genes determined by genome-wide gene essentiality screens, DEG 10 also harbors essential noncoding RNAs, promoters, regulatory sequences and replication origins. These essential genomic elements are determined not only in vitro, but also in vivo, under diverse conditions including those for survival, pathogenesis and antibiotic resistance. We have developed customizable BLAST tools that allow users to perform species- and experiment-specific BLAST searches for a single gene, a list of genes, annotated or unannotated genomes. Therefore, DEG 10 includes essential genomic elements under different conditions in three domains of life, with customizable BLAST tools.

CEG: a database of essential gene clusters

Background

Essential genes are indispensable for the survival of living entities. They are the cornerstones of synthetic biology, and are potential candidate targets for antimicrobial and vaccine design.

Description

Here we describe the Cluster of Essential Genes (CEG) database, which contains clusters of orthologous essential genes. Based on the size of a cluster, users can easily decide whether an essential gene is conserved in multiple bacterial species or is species-specific. It contains the similarity value of every essential gene cluster against human proteins or genes. The CEG_Match tool is based on the CEG database, and was developed for prediction of essential genes according to function. The database is available at http://cefg.uestc.edu.cn/ceg.

Conclusions

Properties contained in the CEG database, such as cluster size, and the similarity of essential gene clusters against human proteins or genes, are very important for evolutionary research and drug design. An advantage of CEG is that it clusters essential genes based on function, and therefore decreases false positive results when predicting essential genes in comparison with using the similarity alignment method.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-769) contains supplementary material, which is available to authorized users.

Common duckweed (Lemna minor) is a versatile high-throughput infection model for the Burkholderia cepacia complex and other pathogenic bacteria

Members of the Burkholderia cepacia complex (Bcc) have emerged in recent decades as problematic pulmonary pathogens of cystic fibrosis (CF) patients, with severe infections progressing to acute necrotizing pneumonia and sepsis. This study presents evidence that Lemna minor (Common duckweed) is useful as a plant model for the Bcc infectious process, and has potential as a model system for bacterial pathogenesis in general. To investigate the relationship between Bcc virulence in duckweed and Galleria mellonella (Greater wax moth) larvae, a previously established Bcc infection model, a duckweed survival assay was developed and used to determine LD₅₀ values. A strong correlation (R² = 0.81) was found between the strains’ virulence ranks in the two infection models, suggesting conserved pathways in these vastly different hosts. To broaden the application of the duckweed model, enteropathogenic Escherichia coli (EPEC) and five isogenic mutants with previously established LD₅₀ values in the larval model were tested against duckweed, and a strong correlation (R² = 0.93) was found between their raw LD₅₀ values. Potential virulence factors in B. cenocepacia K56-2 were identified using a high-throughput screen against single duckweed plants. In addition to the previously characterized antifungal compound (AFC) cluster genes, several uncharacterized genes were discovered including a novel lysR regulator, a histidine biosynthesis gene hisG, and a gene located near the gene encoding the recently characterized virulence factor SuhB_Bc. Finally, to demonstrate the utility of this model in therapeutic applications, duckweed was rescued from Bcc infection by treating with bacteriophage at 6-h intervals. It was observed that phage application became ineffective at a timepoint that coincided with a sharp increase in bacterial invasion of plant tissue. These results indicate that common duckweed can serve as an effective infection model for the investigation of bacterial virulence factors and therapeutic strategies to combat them.

Comparative metabolic systems analysis of pathogenic Burkholderia

Burkholderia cenocepacia and Burkholderia multivorans are opportunistic drug-resistant pathogens that account for the majority of Burkholderia cepacia complex infections in cystic fibrosis patients and also infect other immunocompromised individuals. While they share similar genetic compositions, B. cenocepacia and B. multivorans exhibit important differences in pathogenesis. We have developed reconciled genome-scale metabolic network reconstructions of B. cenocepacia J2315 and B. multivorans ATCC 17616 in parallel (designated iPY1537 and iJB1411, respectively) to compare metabolic abilities and contextualize genetic differences between species. The reconstructions capture the metabolic functions of the two species and give insight into similarities and differences in their virulence and growth capabilities. The two reconstructions have 1,437 reactions in common, and iPY1537 and iJB1411 have 67 and 36 metabolic reactions unique to each, respectively. After curating the extensive reservoir of metabolic genes in Burkholderia, we identified 6 genes essential to growth that are unique to iPY1513 and 13 genes uniquely essential to iJB1411. The reconstructions were refined and validated by comparing in silico growth predictions to in vitro growth capabilities of B. cenocepacia J2315, B. cenocepacia K56-2, and B. multivorans ATCC 17616 on 104 carbon sources. Overall, we identified functional pathways that indicate B. cenocepacia can produce a wider array of virulence factors compared to B. multivorans, which supports the clinical observation that B. cenocepacia is more virulent than B. multivorans. The reconciled reconstructions provide a framework for generating and testing hypotheses on the metabolic and virulence capabilities of these two related emerging pathogens.

Geptop: a gene essentiality prediction tool for sequenced bacterial genomes based on orthology and phylogeny

Integrative genomics predictors, which score highly in predicting bacterial essential genes, would be unfeasible in most species because the data sources are limited. We developed a universal approach and tool designated Geptop, based on orthology and phylogeny, to offer gene essentiality annotations. In a series of tests, our Geptop method yielded higher area under curve (AUC) scores in the receiver operating curves than the integrative approaches. In the ten-fold cross-validations among randomly upset samples, Geptop yielded an AUC of 0.918, and in the cross-organism predictions for 19 organisms Geptop yielded AUC scores between 0.569 and 0.959. A test applied to the very recently determined essential gene dataset from the Porphyromonas gingivalis, which belongs to a phylum different with all of the above 19 bacterial genomes, gave an AUC of 0.77. Therefore, Geptop can be applied to any bacterial species whose genome has been sequenced. Compared with the essential genes uniquely identified by the lethal screening, the essential genes predicted only by Gepop are associated with more protein-protein interactions, especially in the three bacteria with lower AUC scores (<0.7). This may further illustrate the reliability and feasibility of our method in some sense. The web server and standalone version of Geptop are available at http://cefg.uestc.edu.cn/geptop/ free of charge. The tool has been run on 968 bacterial genomes and the results are accessible at the website.

Burkholderia cenocepacia conditional growth mutant library created by random promoter replacement of essential genes

Identification of essential genes by construction of conditional knockouts with inducible promoters allows the identification of essential genes and creation of conditional growth (CG) mutants that are then available as genetic tools for further studies. We used large-scale transposon delivery of the rhamnose-inducible promoter, PrhaB followed by robotic screening of rhamnose-dependent growth to construct a genomic library of 106 Burkholderia cenocepacia CG mutants. Transposon insertions were found where PrhaB was in the same orientation of widely conserved, well-characterized essential genes as well as genes with no previous records of essentiality in other microorganisms. Using previously reported global gene-expression analyses, we demonstrate that PrhaB can achieve the wide dynamic range of expression levels required for essential genes when the promoter is delivered randomly and mutants with rhamnose-dependent growth are selected. We also show specific detection of the target of an antibiotic, novobiocin, by enhanced sensitivity of the corresponding CG mutant (PrhaB controlling gyrB expression) within the library. Modulation of gene expression to achieve 30-60% of wild-type growth created conditions for specific hypersensitivity demonstrating the value of the CG mutant library for chemogenomic experiments. In summary, CG mutants can be obtained on a large scale by random delivery of a tightly regulated inducible promoter into the bacterial chromosome followed by a simple screening for the CG phenotype, without previous information on gene essentiality.

Identification of quorum sensing-controlled genes in Burkholderia ambifaria

The Burkholderia cepacia complex (Bcc) comprises strains with a virulence potential toward immunocompromised patients as well as plant growth-promoting rhizobacteria (PGPR). Owing to the link between quorum sensing (QS) and virulence, most studies among Bcc species have been directed toward QS of pathogenic bacteria. We have investigated the QS of B. ambifaria, a PGPR only infrequently recovered from patients. The cepI gene, responsible for the synthesis of the main signaling molecule N-octanoylhomoserine lactone (C8 -HSL), was inactivated. Phenotypes of the B. ambifaria cepI mutant we observed, such as increased production of siderophores and decreased proteolytic and antifungal activities, are in agreement with those of other Bcc cepI mutants. The cepI mutant was then used as background strain for a whole-genome transposon-insertion mutagenesis strategy, allowing the identification of 20 QS-controlled genes, corresponding to 17 loci. The main functions identified are linked to antifungal and antimicrobial properties, as we have identified QS-controlled genes implicated in the production of pyrrolnitrin, burkholdines (occidiofungin-like molecules), and enacyloxins. This study provides insights in the QS-regulated functions of a PGPR, which could lead to beneficial potential biotechnological applications.

Combining functional and structural genomics to sample the essential Burkholderia structome

Background

The genus Burkholderia includes pathogenic gram-negative bacteria that cause melioidosis, glanders, and pulmonary infections of patients with cancer and cystic fibrosis. Drug resistance has made development of new antimicrobials critical. Many approaches to discovering new antimicrobials, such as structure-based drug design and whole cell phenotypic screens followed by lead refinement, require high-resolution structures of proteins essential to the parasite.

Methodology/Principal Findings

We experimentally identified 406 putative essential genes in B. thailandensis, a low-virulence species phylogenetically similar to B. pseudomallei, the causative agent of melioidosis, using saturation-level transposon mutagenesis and next-generation sequencing (Tn-seq). We selected 315 protein products of these genes based on structure-determination criteria, such as excluding very large and/or integral membrane proteins, and entered them into the Seattle Structural Genomics Center for Infection Disease (SSGCID) structure determination pipeline. To maximize structural coverage of these targets, we applied an “ortholog rescue” strategy for those producing insoluble or difficult to crystallize proteins, resulting in the addition of 387 orthologs (or paralogs) from seven other Burkholderia species into the SSGCID pipeline. This structural genomics approach yielded structures from 31 putative essential targets from B. thailandensis, and 25 orthologs from other Burkholderia species, yielding an overall structural coverage for 49 of the 406 essential gene families, with a total of 88 depositions into the Protein Data Bank. Of these, 25 proteins have properties of a potential antimicrobial drug target i.e., no close human homolog, part of an essential metabolic pathway, and a deep binding pocket. We describe the structures of several potential drug targets in detail.

Conclusions/Significance

This collection of structures, solubility and experimental essentiality data provides a resource for development of drugs against infections and diseases caused by Burkholderia. All expression clones and proteins created in this study are freely available by request.