Suppression-subtractive hybridisation reveals variations in gene distribution amongst the Burkholderia cepacia complex, including the presence in some strains of a genomic island containing putative polysaccharide production genes

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.

Insertion sequence evolutionary patterns highlight convergent genetic inactivations and recent genomic island acquisitions among epidemic Burkholderia cenocepacia

The Burkholderia cenocepacia B&B clone was found previously to be responsible for an epidemic outbreak within an intensive care unit in France. This clone belongs to the ST32 clonal complex, which is one of the most prevalent among French cystic fibrosis patients and is known to be related to the highly virulent ET12 clonal complex. Genomic repartition biases of insertion sequences (ISs) were investigated to improve our understanding of the evolutionary events leading to B. cenocepacia diversification and the emergence of such epidemic lineages. IS were used for tracking convergent genetic inactivations and recent DNA acquisitions. B. cenocepacia IS families and subgroups were compared in terms of genetic diversity and genomic architecture using fully sequenced genomes, PCR screening and DNA blot analysis. These analyses revealed several features shared by the B&B and ET12 epidemic clones. IS elements showed a frequent localization on genomic islands (GI) and indicated convergent evolution towards inactivation of certain loci. The IS407 subgroup of the IS3 family was identified as a good indicator of recently acquired GIs in clone ET12. Several IS407 elements showed strain-specific or clonal complex-specific localizations. IS407 DNA probing of a DNA library built from the B. cenocepacia B&B epidemic clone led to the identification of a recently acquired IS407-tagged GI likely to be conjugative and integrative. The B&B clone showed significant differences in its IS architecture from that of ST32 strains isolated from Czech cystic fibrosis patients.

Genomic island variability facilitates Prochlorococcus-virus coexistence

Prochlorococcus cyanobacteria are extremely abundant in the oceans, as are the viruses that infect them. How hosts and viruses coexist in nature remains unclear, although the presence of both susceptible and resistant cells may allow this coexistence. Combined whole-genome sequencing and PCR screening technology now enables us to investigate the effect of resistance on genome evolution and the genomic mechanisms behind the long-term coexistence of Prochlorococcus and their viruses. Here we present a genome analysis of 77 substrains selected for resistance to ten viruses, revealing mutations primarily in non-conserved, horizontally transferred genes that localize to a single hypervariable genomic island. Mutations affected viral attachment to the cell surface and imposed a fitness cost to the host, manifested by significantly lower growth rates or a previously unknown mechanism of more rapid infection by other viruses. The mutant genes are generally uncommon in nature yet some carry polymorphisms matching those found experimentally. These data are empirical evidence indicating that viral-attachment genes are preferentially located in genomic islands and that viruses are a selective pressure enhancing the diversity of both island genes and island gene content. This diversity emerges as a genomic mechanism that reduces the effective host population size for infection by a given virus, thus facilitating long-term coexistence between viruses and their hosts in nature.

Exploring the metabolic network of the epidemic pathogen Burkholderia cenocepacia J2315 via genome-scale reconstruction

BACKGROUND

Burkholderia cenocepacia is a threatening nosocomial epidemic pathogen in patients with cystic fibrosis (CF) or a compromised immune system. Its high level of antibiotic resistance is an increasing concern in treatments against its infection. Strain B. cenocepacia J2315 is the most infectious isolate from CF patients. There is a strong demand to reconstruct a genome-scale metabolic network of B. cenocepacia J2315 to systematically analyze its metabolic capabilities and its virulence traits, and to search for potential clinical therapy targets.

RESULTS

We reconstructed the genome-scale metabolic network of B. cenocepacia J2315. An iterative reconstruction process led to the establishment of a robust model, iKF1028, which accounts for 1,028 genes, 859 internal reactions, and 834 metabolites. The model iKF1028 captures important metabolic capabilities of B. cenocepacia J2315 with a particular focus on the biosyntheses of key metabolic virulence factors to assist in understanding the mechanism of disease infection and identifying potential drug targets. The model was tested through BIOLOG assays. Based on the model, the genome annotation of B. cenocepacia J2315 was refined and 24 genes were properly re-annotated. Gene and enzyme essentiality were analyzed to provide further insights into the genome function and architecture. A total of 45 essential enzymes were identified as potential therapeutic targets.

CONCLUSIONS

As the first genome-scale metabolic network of B. cenocepacia J2315, iKF1028 allows a systematic study of the metabolic properties of B. cenocepacia and its key metabolic virulence factors affecting the CF community. The model can be used as a discovery tool to design novel drugs against diseases caused by this notorious pathogen.

A putative genomic island, PGI-1, in Ralstonia solanacearum biovar 2 revealed by subtractive hybridization

Ralstonia solanacearum biovar 2, a key bacterial pathogen of potato, has recently established in temperate climate waters. On the basis of isolates obtained from diseased (potato) plants, its genome has been assumed to be virtually clonal, but information on environmental isolates has been lacking. Based on differences in pulsed-field gel electrophoresis patterns, we compared the genomes of two biovar 2 strains with different life histories. Thus, genomic DNA of the novel environmental strain KZR-5 (The Netherlands) was compared to that of reference potato strain 715 (Bangladesh) by suppressive subtractive hybridization. Various strain-specific sequences were found, all being homologous to those found in the genome of reference potato strain 1609. Approximately 20% of these were related to genes involved in recombinational processes. We found a deletion of a 17.6-Kb region, denoted as a putative genomic island PGI-1, in environmental strain KZR-5. The deleted region was, at both extremes, flanked by a composite of two insertion sequence (IS) elements, identified as ISRso2 and ISRso3. The PGI-1 region contained open reading frames that putatively encoded a (p)ppGpp synthetase, a transporter protein, a transcriptional regulator, a cellobiohydrolase, a site-specific integrase/recombinase, a phage-related protein and seven hypothetical proteins. As yet, no phenotype could be assigned to the loss of PGI-1. The ecological behavior of strain KZR-5 was compared to that of reference strain 715. Strain KZR-5 showed enhanced tolerance to 4°C as compared to the reference strain, but was not affected in its virulence on tomato.

The genome of Burkholderia cenocepacia J2315, an epidemic pathogen of cystic fibrosis patients

Bacterial infections of the lungs of cystic fibrosis (CF) patients cause major complications in the treatment of this common genetic disease. Burkholderia cenocepacia infection is particularly problematic since this organism has high levels of antibiotic resistance, making it difficult to eradicate; the resulting chronic infections are associated with severe declines in lung function and increased mortality rates. B. cenocepacia strain J2315 was isolated from a CF patient and is a member of the epidemic ET12 lineage that originated in Canada or the United Kingdom and spread to Europe. The 8.06-Mb genome of this highly transmissible pathogen comprises three circular chromosomes and a plasmid and encodes a broad array of functions typical of this metabolically versatile genus, as well as numerous virulence and drug resistance functions. Although B. cenocepacia strains can be isolated from soil and can be pathogenic to both plants and man, J2315 is representative of a lineage of B. cenocepacia rarely isolated from the environment and which spreads between CF patients. Comparative analysis revealed that ca. 21% of the genome is unique in comparison to other strains of B. cenocepacia, highlighting the genomic plasticity of this species. Pseudogenes in virulence determinants suggest that the pathogenic response of J2315 may have been recently selected to promote persistence in the CF lung. The J2315 genome contains evidence that its unique and highly adapted genetic content has played a significant role in its success as an epidemic CF pathogen.

Plant host and sugar alcohol induced exopolysaccharide biosynthesis in the Burkholderia cepacia complex

The species that presently constitute the Burkholderia cepacia complex (Bcc) have multiple roles; they include soil and water saprophytes, bioremediators, and plant, animal and human pathogens. Since the first description of pathogenicity in the Bcc was based on sour skin rot of onion bulbs, this study returned to this plant host to investigate the onion-associated phenotype of the Bcc. Many Bcc isolates, which were previously considered to be non-mucoid, produced copious amounts of exopolysaccharide (EPS) when onion tissue was provided as the sole nutrient. EPS production was not species-specific, was observed in isolates from both clinical and environmental sources, and did not correlate with the ability to cause maceration of onion tissue. Chemical analysis suggested that the onion components responsible for EPS induction were primarily the carbohydrates sucrose, fructose and fructans. Additional sugars were investigated, and all alcohol sugars tested were able to induce EPS production, in particular mannitol and glucitol. To investigate the molecular basis for EPS biosynthesis, we focused on the highly conserved bce gene cluster thought to be involved in cepacian biosynthesis. We demonstrated induction of the bce gene cluster by mannitol, and found a clear correlation between the inability of representatives of the Burkholderia cenocepacia ET12 lineage to produce EPS and the presence of an 11 bp deletion within the bceB gene, which encodes a glycosyltransferase. Insertional inactivation of bceB in Burkholderia ambifaria AMMD results in loss of EPS production on sugar alcohol media. These novel and surprising insights into EPS biosynthesis highlight the metabolic potential of the Bcc and show that a potential virulence factor may not be detected by routine laboratory culture. Our results also highlight a potential hazard in the use of inhaled mannitol as an osmolyte to improve mucociliary clearance in individuals with cystic fibrosis.

Application of suppressive subtractive hybridization to uncover the metagenomic diversity of environmental samples

Metagenomics addresses the collective genetic structure and functional composition of a microbial environmental sample without the bias or necessity for culturing the microorganisms from the community in question. Metagenomic studies are now beginning to take advantage of the plethora of complete genome sequences and the associated tools, such as bacterial artificial chromosome (BAC) and fosmid vectors, to discover novel genes and survey the structure and function of microbial communities. Complementary and less expensive methods to compare genomes from individual microbes have been utilized in comparative genomic studies. Suppressive subtractive hybridization (SSH) is one such approach, which has been utilized to compare the genomic content of closely related species of bacteria. Recently, SSH has also been used as a comparative method to examine the microbial diversity (i.e., species composition) and functional differences (i.e., gene composition) in the genomic content of two different rumen environmental communities. Through a series of hybridizations and pblymerase chain reaction (PCR) amplifications, metagenomic differences between two environmental samples can be isolated by SSH. Subsequent DNA sequencing and bioinformatic analyses allow the putative identification of these differences.

Attenuation and protective efficacy of an O-antigen-deficient mutant of Francisella tularensis LVS

Francisella tularensis is a zoonotic, Gram-negative coccobacillus that causes tularemia in humans and animals. F. tularensis subspecies tularensis (type A) and F. tularensis subspecies holarctica (type B) are antigenically similar and more virulent than Francisella novicida in humans. The genetic locus that encodes the LPS O antigen was found to be substantially different between the type B live vaccine strain (LVS) and F. novicida. One LVS-specific gene with homology to a galactosyl transferase was selected for allelic replacement using a sacB-chloramphenicol expression suicide plasmid, and recombinants were screened for colony morphology on Congo red agar that matched that of F. novicida. Two mutants (WbtI(S187Y) and WbtI(G191V)) were isolated that contained substitutions in conserved motifs in the sugar transamine/perosamine synthetase (WbtI) of the O-antigen locus, and the latter mutant was extensively tested and characterized. WbtI(G191V) grew at the same rate as the parent strain in Chamberlain's defined medium, completely lacked O antigen, was serum-sensitive but could grow in a mouse macrophage cell line, had increased resistance to sodium deoxycholate, and was highly attenuated in mice. Complementation of WbtI(G191V) with the wild-type wbtI gene in trans restored normal LPS synthesis, phenotypic properties similar to the parent, and virulence in mice. Immunization with WbtI(G191V) protected mice against a relatively low-dose intraperitoneal challenge with LVS, but was less protective against a high-dose challenge. These results indicate that complete loss of O antigen alters the surface phenotype and abrogates virulence in F. tularensis, but also compromises the induction of full protective immunity against F. tularensis infection in mice.

The Burkholderia cepacia bceA gene encodes a protein with phosphomannose isomerase and GDP-D-mannose pyrophosphorylase activities

The bceA gene is part of the Burkholderia cepacia IST408 exopolysaccharide (EPS) biosynthetic cluster. It encodes a 55.3-kDa bifunctional protein (type II PMI family) with phosphomannose isomerase (PMI) and GDP-mannose pyrophosphorylase (GMP) activities. GMP activity is strongly dependent on the presence of Ca(2+) or Mn(2+), while PMI activity can use a broader variety of divalent cations (Ca(2+)>Mn(2+)>Mg(2+)>Co(2+)>Ni(2+)). The lack of a functional bceA gene does not affect EPS production yield in a non-polar insertion bceA mutant. The in silico search for putative bceA homologues revealed the presence of 2-5 bceA orthologues in the Burkholderia genomes available. This suggests that in B. cepacia IST408 putative bceA functional homologues may compensate the bceA mutation. However, the viscosity of aqueous solutions prepared with the EPS produced by the bceA mutant was significantly reduced compared with wild-type biopolymer and the mutant forms biofilms with a size reduced by 6-fold.

A suppression subtractive hybridization approach reveals niche-specific genes that may be involved in predator avoidance in marine Synechococcus isolates

Picocyanobacteria of the genus Synechococcus are important contributors to marine primary production and are ubiquitous in the world's oceans. This genus is genetically diverse, and at least 10 discrete lineages or clades have been identified phylogenetically. However, little if anything is known about the genetic attributes which characterize particular lineages or are unique to specific strains. Here, we used a suppression subtractive hybridization (SSH) approach to identify strain- and clade-specific genes in two well-characterized laboratory strains, Synechococcus sp. strain WH8103 (clade III) and Synechococcus sp. strain WH7803 (clade V). Among the genes that were identified as potentially unique to each strain were genes encoding proteins that may be involved in specific predator avoidance, including a glycosyltransferase in strain WH8103 and a permease component of an ABC-type polysaccharide/polyol phosphate export system in WH7803. During this work the genome of one of these strains, WH7803, became available. This allowed assessment of the number of false-positive sequences (i.e., sequences present in the tester genome) present among the SSH-enriched sequences. We found that approximately 9% of the WH8103 sequences were potential false-positive sequences, which demonstrated that caution should be used when this technology is used to assess genomic differences in genetically similar bacterial strains.

Spatial epidemiology and natural population structure of Campylobacter jejuni colonizing a farmland ecosystem

Recent progress in determining the population structure of Campylobacter jejuni, and discerning associations between genotypes and specific niches, has emphasized the shortfall in our understanding of the ecology and epidemiology of this bacterium. We examined the natural structure of the C. jejuni community associated with cattle farmland in the UK by structured spatiotemporal sampling of habitats, including livestock and wild animal faeces, environmental water and soil, over a 10-week period within a 100 km2 area. A total of 172 isolates were characterized using multilocus sequence typing into 65 sequence types (STs). Isolates from cattle faeces were significantly over-represented in the ST-61 complex, whereas isolates from wildlife faeces and water were more likely to belong to the ST-45 complex and a number of unusual STs, many of which were first encountered during this study. Sampling within a narrow spatiotemporal window permitted the application of novel statistical methods exploring the relationship between the genetic relatedness and spatial separation of isolates. This approach showed that isolates from the same sampling squares and squares separated by <1.0 km were genetically more similar than isolates separated by greater distances. Our study demonstrates the potential of multilocus sequence typing combined with spatial modelling in exploring natural transmission pathways for C. jejuni.

Use of suppression-subtractive hybridization to identify genes in the Burkholderia cepacia complex that are unique to Burkholderia cenocepacia

We have previously shown differences in virulence between species of the Burkholderia cepacia complex using the alfalfa infection model and the rat agar bead chronic infection model. Burkholderia cenocepacia strains were more virulent in these two infection models than Burkholderia multivorans and Burkholderia stabilis strains. In order to identify genes that may account for the increased virulence of B. cenocepacia, suppression-subtractive hybridization was performed between B. cenocepacia K56-2 and B. multivorans C5393 and between B. cenocepacia K56-2 and B. stabilis LMG14294. Genes identified included DNA modification/phage-related/insertion sequences and genes involved in cell membrane/surface structures, resistance, transport, metabolism, regulation, secretion systems, as well as genes of unknown function. Several of these genes were present in the ET12 lineage of B. cenocepacia but not in other members of the B. cepacia complex. Virulence studies in a chronic lung infection model determined that the hypothetical YfjI protein, which is unique to the ET12 clone, contributes to lung pathology. Other genes specific to B. cenocepacia and/or the ET12 lineage were shown to play a role in biofilm formation and swarming or swimming motility.

Identification and physical organization of the gene cluster involved in the biosynthesis of Burkholderia cepacia complex exopolysaccharide

Bacteria belonging to the Burkholderia cepacia complex (BCC) are important opportunistic pathogens in patients with cystic fibrosis (CF). Since approximately 80% of the CF isolates examined produce exopolysaccharide (EPS), it was hypothesized that this EPS may play a role in the colonization and persistence of these bacteria in the CF lung. The present study describes the identification and physical organization of the EPS biosynthetic gene cluster. This bce gene cluster was identified following the isolation of three EPS-defective mutants from the highly mucoid CF isolate IST408, belonging to BCC genomovar I, based on random plasposon insertion mutagenesis and comparison of the nucleotide sequence of the interrupted genes with the available genome of Burkholderia cenocepacia J2315. This 16.2 kb cluster includes 12 genes and is located on chromosome 2. Database searches for homologous proteins and secondary structure analysis for the deduced Bce amino acid sequences revealed genes predicted to encode enzymes required for the formation of nucleotide sugar precursors, glycosyltransferases involved in the repeat-unit assembly, and other proteins involved in polymerization and export of bacterial surface polysaccharides.

Molecular typing of, and distribution of genetic markers among, Burkholderia cepacia complex isolates from Brazil

PCR tests were used to assign genomovar status to 39 non-cystic fibrosis (non-CF) and 11 CF Burkholderia cepacia complex isolates from patients in hospitals in Recife, Brazil. Non-CF isolates were assigned to genomovar IIIA (71.8%), genomovar I (15.4%), B. vietnamiensis (7.7%), and B. multivorans (5.1%). CF isolates were assigned to genomovar IIIA (18.2%), B. vietnamiensis (18.2%), and genomovar I (9.1%). Six CF isolates sharing recA PCR-restriction fragment length polymorphism (RFLP) and randomly amplified polymorphic DNA (RAPD) patterns could not be assigned to a genomovar. 16S rDNA sequence obtained from these isolates indicated a closest relationship to B. anthina, but the recA sequence was equally divergent from several genomovars. PCR screening indicated the presence of cblA in only two isolates, whereas the B. cepacia epidemic strain marker was found in 22 of 28 genomovar IIIA isolates. A type III secretion gene was detected in all but genomovar I isolates. RAPD and PCR-RFLP assays, targeting both recA and fliC, indicated a large amount of genetic variability among the isolates, with many novel patterns being observed. Nine genomovar IIIA isolates from different non-CF patients and clinical sources had identical genotypes, indicating the presence of a common clone.