Chromosome translocation and its consequence in the genome of Burkholderia cenocepacia AU-1054

Chromids Aid Genome Expansion and Functional Diversification in the Family Burkholderiaceae

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

Genome rearrangements and selection in multi-chromosome bacteria Burkholderia spp

BACKGROUND

The genus Burkholderia consists of species that occupy remarkably diverse ecological niches. Its best known members are important pathogens, B. mallei and B. pseudomallei, which cause glanders and melioidosis, respectively. Burkholderia genomes are unusual due to their multichromosomal organization, generally comprised of 2-3 chromosomes.

RESULTS

We performed integrated genomic analysis of 127 Burkholderia strains. The pan-genome is open with the saturation to be reached between 86,000 and 88,000 genes. The reconstructed rearrangements indicate a strong avoidance of intra-replichore inversions that is likely caused by selection against the transfer of large groups of genes between the leading and the lagging strands. Translocated genes also tend to retain their position in the leading or the lagging strand, and this selection is stronger for large syntenies. Integrated reconstruction of chromosome rearrangements in the context of strains phylogeny reveals parallel rearrangements that may indicate inversion-based phase variation and integration of new genomic islands. In particular, we detected parallel inversions in the second chromosomes of B. pseudomallei with breakpoints formed by genes encoding membrane components of multidrug resistance complex, that may be linked to a phase variation mechanism. Two genomic islands, spreading horizontally between chromosomes, were detected in the B. cepacia group.

CONCLUSIONS

This study demonstrates the power of integrated analysis of pan-genomes, chromosome rearrangements, and selection regimes. Non-random inversion patterns indicate selective pressure, inversions are particularly frequent in a recent pathogen B. mallei, and, together with periods of positive selection at other branches, may indicate adaptation to new niches. One such adaptation could be a possible phase variation mechanism in B. pseudomallei.

A Comprehensive Overview of Online Resources to Identify and Predict Bacterial Essential Genes

Genes critical for the survival or reproduction of an organism in certain circumstances are classified as essential genes. Essential genes play a significant role in deciphering the survival mechanism of life. They may be greatly applied to pharmaceutics and synthetic biology. The continuous progress of experimental method for essential gene identification has accelerated the accumulation of gene essentiality data which facilitates the study of essential genes in silico. In this article, we present some available online resources related to gene essentiality, including bioinformatic software tools for transposon sequencing (Tn-seq) analysis, essential gene databases and online services to predict bacterial essential genes. We review several computational approaches that have been used to predict essential genes, and summarize the features used for gene essentiality prediction. In addition, we evaluate the available online bacterial essential gene prediction servers based on the experimentally validated essential gene sets of 30 bacteria from DEG. This article is intended to be a quick reference guide for the microbiologists interested in the essential genes.

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.

Three computational tools for predicting bacterial essential genes

Essential genes are those genes indispensable for the survival of any living cell. Bacterial essential genes constitute the cornerstones of synthetic biology and are often attractive targets in the development of antibiotics and vaccines. Because identification of essential genes with wet-lab ways often means expensive economic costs and tremendous labor, scientists changed to seek for alternative way of computational prediction. Aiming to help to solve this issue, our research group (CEFG: group of Computational, Comparative, Evolutionary and Functional Genomics, http://cefg.uestc.edu.cn) has constructed three online services to predict essential genes in bacterial genomes. These freely available tools are applicable for single gene sequences without annotated functions, single genes with definite names, and complete genomes of bacterial strains. To ensure reliable predictions, the investigated species should belong to the same family (for EGP) or phylum (for CEG_Match and Geptop) with one of the reference species, respectively. As the pilot software for the issue, predicting accuracies of them have been assessed and compared with existing algorithms, and note that all of other published algorithms have not any formed online services. We hope these services at CEFG will help scientists and researchers in the field of essential genes.

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.

Evolutionary effects of translocations in bacterial genomes

It has become clear that different genome regions need not evolve uniformly. This variation is particularly evident in bacterial genomes with multiple chromosomes, in which smaller, secondary chromosomes evolve more rapidly. We previously demonstrated that substitution rates and gene dispensability were greater on secondary chromosomes in many bacterial genomes. In Vibrio, the secondary chromosome is replicated later during the cell cycle, which reduces the effective dosage of these genes and hence their expression. More rapid evolution of secondary chromosomes may therefore reflect weaker purifying selection on less expressed genes. Here, we test this hypothesis by relating substitution rates of orthologs shared by multiple Burkholderia genomes, each with three chromosomes, to a study of gene expression in genomes differing by a major reciprocal translocation. This model predicts that expression should be greatest on chromosome 1 (the largest) and least on chromosome 3 (the smallest) and that expression should tend to decline within chromosomes from replication origin to terminus. Moreover, gene movement to the primary chromosome should associate with increased expression, and movement to secondary chromosomes should result in reduced expression. Our analysis supports each of these predictions, as translocated genes tended to shift expression toward their new chromosome neighbors despite inevitable cis-acting regulation of expression. This study sheds light on the early dynamics of genomes following rearrangement and illustrates how secondary chromosomes in bacteria may become evolutionary test beds.

The tRNAarg gene and engA are essential genes on the 1.7-Mb pSymB megaplasmid of Sinorhizobium meliloti and were translocated together from the chromosome in an ancestral strain

Bacterial genomes with two (or more) chromosome-like replicons are known, and these appear to be particularly frequent in alphaproteobacteria. The genome of the N(2)-fixing alfalfa symbiont Sinorhizobium meliloti 1021 contains a 3.7-Mb chromosome and 1.4-Mb (pSymA) and 1.7-Mb (pSymB) megaplasmids. In this study, the tRNA(arg) and engA genes, located on the pSymB megaplasmid, are shown to be essential for growth. These genes could be deleted from pSymB when copies were previously integrated into the chromosome. However, in the closely related strain Sinorhizobium fredii NGR234, the tRNA(arg) and engA genes are located on the chromosome, in a 69-kb region designated the engA-tRNA(arg)-rmlC region. This region includes bacA, a gene that is important for intracellular survival during host-bacterium interactions for S. meliloti and the related alphaproteobacterium Brucella abortus. The engA-tRNA(arg)-rmlC region lies between the kdgK and dppF2 (NGR_c24410) genes on the S. fredii chromosome. Synteny analysis showed that kdgK and dppF2 orthologues are adjacent to each other on the chromosomes of 15 sequenced strains of S. meliloti and Sinorhizobium medicae, whereas the 69-kb engA-tRNA(arg)-rmlC region is present on the pSymB-equivalent megaplasmids. This and other evidence strongly suggests that the engA-tRNA(arg)-rmlC region translocated from the chromosome to the progenitor of pSymB in an ancestor common to S. meliloti and S. medicae. To our knowledge, this work represents one of the first experimental demonstrations that essential genes are present on a megaplasmid.

[Current status of theoretical studies on essential genes in microbes]

Essential genes are indispensable for the survival of an organism in optimal conditions. Recently, study on essential gene is becoming a hot topic of microbiology, genomics, and bioinformatics. This paper described the experiments that determined essential genes in some microbes and the theoretical researches on essential genes were reviewed. The major content contained comparison of essential genes and non-essential genes based on information on evolutionary conservation and sequence composition, and in silico prediction of essential genes, and analysis of the chromosomal distributions of essential genes. Finally, related progresses were concluded and the open problems were pointed out.