Interrogation of the Burkholderia pseudomallei Genome to Address Differential Virulence among Isolates

Infection by the Gram-negative pathogen Burkholderia pseudomallei results in the disease melioidosis, acquired from the environment in parts of southeast Asia and northern Australia. Clinical symptoms of melioidosis range from acute (fever, pneumonia, septicemia, and localized infection) to chronic (abscesses in various organs and tissues, most commonly occurring in the lungs, liver, spleen, kidney, prostate and skeletal muscle), and persistent infections in humans are difficult to cure. Understanding the basic biology and genomics of B. pseudomallei is imperative for the development of new vaccines and therapeutic interventions. This formidable task is becoming more tractable due to the increasing number of B. pseudomallei genomes that are being sequenced and compared.Here, we compared three B. pseudomallei genomes, from strains MSHR668, K96243 and 1106a, to identify features that might explain why MSHR668 is more virulent than K96243 and 1106a in a mouse model of B. pseudomallei infection. Our analyses focused on metabolic, virulence and regulatory genes that were present in MSHR668 but absent from both K96243 and 1106a. We also noted features present in K96243 and 1106a but absent from MSHR668, and identified genomic differences that may contribute to variations in virulence noted among the three B. pseudomallei isolates. While this work contributes to our understanding of B. pseudomallei genomics, more detailed experiments are necessary to characterize the relevance of specific genomic features to B. pseudomallei metabolism and virulence. Functional analyses of metabolic networks, virulence and regulation shows promise for examining the effects of B. pseudomallei on host cell metabolism and will lay a foundation for future prediction of the virulence of emerging strains. Continued emphasis in this area will be critical for protection against melioidosis, as a better understanding of what constitutes a fully virulent Burkholderia isolate may provide for better diagnostic and medical countermeasure strategies.     

Genomic diversity of Burkholderia pseudomallei clinical isolates: subtractive hybridization reveals a Burkholderia mallei-specific prophage in B. pseudomallei 1026b

Burkholderia pseudomallei is the etiologic agent of the disease melioidosis and is a category B biological threat agent. The genomic sequence of B. pseudomallei K96243 was recently determined, but little is known about the overall genetic diversity of this species. Suppression subtractive hybridization was employed to assess the genetic variability between two distinct clinical isolates of B. pseudomallei, 1026b and K96243. Numerous mobile genetic elements, including a temperate bacteriophage designated phi1026b, were identified among the 1026b-specific suppression subtractive hybridization products. Bacteriophage phi1026b was spontaneously produced by 1026b, and it had a restricted host range, infecting only Burkholderia mallei. It possessed a noncontractile tail, an isometric head, and a linear 54,865-bp genome. The mosaic nature of the phi1026b genome was revealed by comparison with bacteriophage phiE125, a B. mallei-specific bacteriophage produced by Burkholderia thailandensis. The phi1026b genes for DNA packaging, tail morphogenesis, host lysis, integration, and DNA replication were nearly identical to the corresponding genes in phiE125. On the other hand, phi1026b genes involved in head morphogenesis were similar to head morphogenesis genes encoded by Pseudomonas putida and Pseudomonas aeruginosa bacteriophages. Consistent with this observation, immunogold electron microscopy demonstrated that polyclonal antiserum against phiE125 reacted with the tail of phi1026b but not with the head. The results presented here suggest that B. pseudomallei strains are genetically heterogeneous and that bacteriophages are major contributors to the genomic diversity of this species. The bacteriophage characterized in this study may be a useful diagnostic tool for differentiating B. pseudomallei and B. mallei, two closely related biological threat agents.     

Genotypes Associated with Virulence in Environmental Isolates of Vibrio cholerae

Vibrio cholerae is an autochthonous inhabitant of riverine and estuarine environments and also is a facultative pathogen for humans. Genotyping can be useful in assessing the risk of contracting cholera, intestinal, or extraintestinal infections via drinking water and/or seafood. In this study, environmental isolates of V. cholerae were examined for the presence of ctxA, hlyA, ompU, stn/sto, tcpA, tcpI, toxR, and zot genes, using multiplex PCR. Based on tcpA and hlyA gene comparisons, the strains could be grouped into Classical and El Tor biotypes. The toxR, hlyA, and ompU genes were present in 100, 98.6, and 87.0% of the V. cholerae isolates, respectively. The CTX genetic element and toxin-coregulated pilus El Tor (tcpA ET) gene were present in all toxigenic V. cholerae O1 and V. cholerae O139 strains examined in this study. Three of four nontoxigenic V. cholerae O1 strains contained tcpA ET. Interestingly, among the isolates of V. cholerae non-O1/non-O139, two had tcpA Classical, nine contained tcpA El Tor, three showed homology with both biotype genes, and four carried the ctxA gene. The stn/sto genes were present in 28.2% of the non-O1/non-O139 strains, in 10.5% of the toxigenic V. cholerae O1, and in 14.3% of the O139 serogroups. Except for stn/sto genes, all of the other genes studied occurred with high frequency in toxigenic V. cholerae O1 and O139 strains. Based on results of this study, surveillance of non-O1/non-O139 V. cholerae in the aquatic environment, combined with genotype monitoring using ctxA, stn/sto, and tcpA ET genes, could be valuable in human health risk assessment.     

Comparative Genomics of Field Isolates of Mycobacterium bovis and M. caprae Provides Evidence for Possible Correlates with Bacterial Viability and Virulence

Mycobacteria of the Mycobacterium tuberculosis complex (MTBC) greatly affect humans and animals worldwide. The life cycle of mycobacteria is complex and the mechanisms resulting in pathogen infection and survival in host cells are not fully understood. Recently, comparative genomics analyses have provided new insights into the evolution and adaptation of the MTBC to survive inside the host. However, most of this information has been obtained using M. tuberculosis but not other members of the MTBC such as M. bovis and M. caprae. In this study, the genome of three M. bovis (MB1, MB3, MB4) and one M. caprae (MB2) field isolates with different lesion score, prevalence and host distribution phenotypes were sequenced. Genome sequence information was used for whole-genome and protein-targeted comparative genomics analysis with the aim of finding correlates with phenotypic variation with potential implications for tuberculosis (TB) disease risk assessment and control. At the whole-genome level the results of the first comparative genomics study of field isolates of M. bovis including M. caprae showed that as previously reported for M. tuberculosis, sequential chromosomal nucleotide substitutions were the main driver of the M. bovis genome evolution. The phylogenetic analysis provided a strong support for the M. bovis/M. caprae clade, but supported M. caprae as a separate species. The comparison of the MB1 and MB4 isolates revealed differences in genome sequence, including gene families that are important for bacterial infection and transmission, thus highlighting differences with functional implications between isolates otherwise classified with the same spoligotype. Strategic protein-targeted analysis using the ESX or type VII secretion system, proteins linking stress response with lipid metabolism, host T cell epitopes of mycobacteria, antigens and peptidoglycan assembly protein identified new genetic markers and candidate vaccine antigens that warrant further study to develop tools to evaluate risks for TB disease caused by M. bovis/M.caprae and for TB control in humans and animals.     

Molecular Characterization of Genetic Loci Required for Secretion of Exoproducts in Burkholderia pseudomallei

Previous studies have demonstrated that Burkholderia pseudomallei secretes protease, lipase, and phospholipase C (PLC) into the extracellular milieu, but their mechanisms of secretion and roles in pathogenesis have not been elucidated. In this study, we isolated and characterized 29 transposon mutants unable to secrete protease, lipase, and PLC.     

Comparative Genomic Characterization of Francisella tularensis Strains Belonging to Low and High Virulence Subspecies

Tularemia is a geographically widespread, severely debilitating, and occasionally lethal disease in humans. It is caused by infection by a gram-negative bacterium, Francisella tularensis. In order to better understand its potency as an etiological agent as well as its potential as a biological weapon, we have completed draft assemblies and report the first complete genomic characterization of five strains belonging to the following different Francisella subspecies (subsp.): the F. tularensis subsp. tularensis FSC033, F. tularensis subsp. holarctica FSC257 and FSC022, and F. tularensis subsp. novicida GA99-3548 and GA99-3549 strains. Here, we report the sequencing of these strains and comparative genomic analysis with recently available public Francisella sequences, including the rare F. tularensis subsp. mediasiatica FSC147 strain isolate from the Central Asian Region. We report evidence for the occurrence of large-scale rearrangement events in strains of the holarctica subspecies, supporting previous proposals that further phylogenetic subdivisions of the Type B clade are likely. We also find a significant enrichment of disrupted or absent ORFs proximal to predicted breakpoints in the FSC022 strain, including a genetic component of the Type I restriction-modification defense system. Many of the pseudogenes identified are also disrupted in the closely related rarely human pathogenic F. tularensis subsp. mediasiatica FSC147 strain, including modulator of drug activity B (mdaB) (FTT0961), which encodes a known NADPH quinone reductase involved in oxidative stress resistance. We have also identified genes exhibiting sequence similarity to effectors of the Type III (T3SS) and components of the Type IV secretion systems (T4SS). One of the genes, msrA2 (FTT1797c), is disrupted in F. tularensis subsp. mediasiatica and has recently been shown to mediate bacterial pathogen survival in host organisms. Our findings suggest that in addition to the duplication of the Francisella Pathogenicity Island, and acquisition of individual loci, adaptation by gene loss in the more recently emerged tularensis, holarctica, and mediasiatica subspecies occurred and was distinct from evolutionary events that differentiated these subspecies, and the novicida subspecies, from a common ancestor. Our findings are applicable to future studies focused on variations in Francisella subspecies pathogenesis, and of broader interest to studies of genomic pathoadaptation in bacteria.     

Comparative Phenotypic, Molecular, and Virulence Characterization of Vibrio parahaemolyticus O3:K6 Isolates

Historically, Vibrio parahaemolyticus infections have been characterized by sporadic cases caused by multiple, diverse serotypes. However, since 1996, V. parahaemolyticus serotype O3:K6 strains have been associated with several large-scale outbreaks of illness, suggesting the emergence of a “new” group of organisms with enhanced virulence. We have applied three different molecular subtyping techniques to identify an appropriate method for differentiating O3:K6 isolates from other serotypes. Pulsed-field gel electrophoresis (PFGE) following NotI digestion differentiated seven closely related subtypes among O3:K6 and related strains, which were distinct from PFGE patterns for non-O3:K6 isolates. Ribotyping and tdh sequencing were less discriminatory than PFGE, but further confirmed close genetic relationships among recent O3:K6 isolates. In vitro adherence and cytotoxicity studies with human epithelial cells showed that O3:K6 isolates exhibited statistically higher levels of adherence and cytotoxicity to host cells than non-O3:K6 isolates. Epithelial cell cytotoxicity patterns were determined with a lactate dehydrogenase release assay. At 3 h postinfection, high relative cytotoxicities (>50% maximum lactate dehydrogenase activity) were found among a greater proportion of recently isolated O3:K6 and closely related strains (75%) than among the non-O3:K6 isolates (23%). A statistically significant relationship between adherence and cytotoxicity suggests that the pathogenic potential of some isolates may be associated with increased adherence to epithelial cells. Our findings suggest that enhanced adherence and cytotoxicity may contribute to the apparent unique pathogenic potential of V. parahaemolyticus O3:K6 strains.     

Identification of Gene Clusters Associated with Host Adaptation and Antibiotic Resistance in Chinese Staphylococcus aureus Isolates by Microarray-Based Comparative Genomics

A comparative genomic microarray comprising 2,457 genes from two whole genomes of S. aureus was employed for the comparative genome hybridization analysis of 50 strains of divergent clonal lineages, including methicillin-resistant S. aureus (MRSA), methicillin-susceptible S. aureus (MSSA), and swine strains in China. Large-scale validation was confirmed via polymerase chain reaction in 160 representative clinical strains. All of the 50 strains were clustered into seven different complexes by phylogenetic tree analysis. Thirteen gene clusters were specific to different S. aureus clones. Ten gene clusters, including seven known (vSa3, vSa4, vSaα, vSaβ, Tn5801, and phage ϕSa3) and three novel (C8, C9, and C10) gene clusters, were specific to human MRSA. Notably, two global regulators, sarH2 and sarH3, at cluster C9 were specific to human MRSA, and plasmid pUB110 at cluster C10 was specific to swine MRSA. Three clusters known to be part of SCCmec, vSa4 or Tn5801, and vSaα as well as one novel gene cluster C12 with homology with Tn554 of S. epidermidis were identified as MRSA-specific gene clusters. The replacement of ST239-spa t037 with ST239-spa t030 in Beijing may be a result of its acquisition of vSa4, phage ϕSa1, and ϕSa3. In summary, thirteen critical gene clusters were identified to be contributors to the evolution of host specificity and antibiotic resistance in Chinese S. aureus.     

Insights into the Dekkera bruxellensis Genomic Landscape: Comparative Genomics Reveals Variations in Ploidy and Nutrient Utilisation Potential amongst Wine Isolates

The yeast Dekkera bruxellensis is a major contaminant of industrial fermentations, such as those used for the production of biofuel and wine, where it outlasts and, under some conditions, outcompetes the major industrial yeast Saccharomyces cerevisiae. In order to investigate the level of inter-strain variation that is present within this economically important species, the genomes of four diverse D. bruxellensis isolates were compared. While each of the four strains was shown to contain a core diploid genome, which is clearly sufficient for survival, two of the four isolates have a third haploid complement of chromosomes. The sequences of these additional haploid genomes were both highly divergent from those comprising the diploid core and divergent between the two triploid strains. Similar to examples in the Saccharomyces spp. clade, where some allotriploids have arisen on the basis of enhanced ability to survive a range of environmental conditions, it is likely these strains are products of two independent hybridisation events that may have involved multiple species or distinct sub-species of Dekkera. Interestingly these triploid strains represent the vast majority (92%) of isolates from across the Australian wine industry, suggesting that the additional set of chromosomes may confer a selective advantage in winery environments that has resulted in these hybrid strains all-but replacing their diploid counterparts in Australian winery settings. In addition to the apparent inter-specific hybridisation events, chromosomal aberrations such as strain-specific insertions and deletions and loss-of-heterozygosity by gene conversion were also commonplace. While these events are likely to have affected many phenotypes across these strains, we have been able to link a specific deletion to the inability to utilise nitrate by some strains of D. bruxellensis, a phenotype that may have direct impacts in the ability for these strains to compete with S. cerevisiae.     

Occurrence of Virulence Genes Associated with Diarrheagenic Pathotypes in Escherichia coli Isolates from Surface Water

Escherichia coli isolates (n = 300) collected from six sites in subtropical Brisbane, Australia, prior to and after storm events were tested for the presence of 11 virulence genes (VGs) specific to diarrheagenic pathotypes. The presence of eaeA, stx₁, stx₂, and ehxA genes specific for the enterohemorrhagic E. coli (EHEC) pathotype was detected in 56%, 6%, 10%, and 13% of isolates, respectively. The VGs astA (69%) and aggR (29%), carried by enteroaggregative (EAEC) pathotypes, were frequently detected in E. coli isolates. The enteropathogenic E. coli (EPEC) gene bfp was detected in 24% of isolates. In addition, enteroinvasive E. coli (EIEC) VG ipaH was also detected in 14% of isolates. During dry periods, isolates belonging to the EAEC pathotype were most commonly detected (23%), followed by EHEC (11%) and EPEC (11%). Conversely, a more uniform prevalence of pathotypes, EPEC (14%), EAEC (12%), EIEC (10%), EHEC (7%), and ETEC (7%), was observed after the storm events. The results of this study highlight the widespread occurrence of potentially diarrheagenic pathotypes in the urban aquatic ecosystems. While the presence of VGs in E. coli isolates alone is insufficient to determine pathogenicity, the presence of diarrheagenic E. coli pathotypes in high frequency after the storm events could lead to increased health risks if untreated storm water were to be used for nonpotable purposes and recreational activities.     

The core and accessory genomes of Burkholderia pseudomallei: implications for human melioidosis

Natural isolates of Burkholderia pseudomallei (Bp), the causative agent of melioidosis, can exhibit significant ecological flexibility that is likely reflective of a dynamic genome. Using whole-genome Bp microarrays, we examined patterns of gene presence and absence across 94 South East Asian strains isolated from a variety of clinical, environmental, or animal sources. 86% of the Bp K96243 reference genome was common to all the strains representing the Bp "core genome", comprising genes largely involved in essential functions (eg amino acid metabolism, protein translation). In contrast, 14% of the K96243 genome was variably present across the isolates. This Bp accessory genome encompassed multiple genomic islands (GIs), paralogous genes, and insertions/deletions, including three distinct lipopolysaccharide (LPS)-related gene clusters. Strikingly, strains recovered from cases of human melioidosis clustered on a tree based on accessory gene content, and were significantly more likely to harbor certain GIs compared to animal and environmental isolates. Consistent with the inference that the GIs may contribute to pathogenesis, experimental mutation of BPSS2053, a GI gene, reduced microbial adherence to human epithelial cells. Our results suggest that the Bp accessory genome is likely to play an important role in microbial adaptation and virulence.     

Proteomic Profiling of Burkholderia cenocepacia Clonal Isolates with Different Virulence Potential Retrieved from a Cystic Fibrosis Patient during Chronic Lung Infection

Respiratory infections with Burkholderia cepacia complex (Bcc) bacteria in cystic fibrosis (CF) are associated with a worse prognosis and increased risk of death. In this work, we assessed the virulence potential of three B. cenocepacia clonal isolates obtained from a CF patient between the onset of infection (isolate IST439) and before death with cepacia syndrome 3.5 years later (isolate IST4113 followed by IST4134), based on their ability to invade epithelial cells and compromise epithelial monolayer integrity. The two clonal isolates retrieved during late-stage disease were significantly more virulent than IST439. Proteomic profiling by 2-D DIGE of the last isolate recovered before the patient’s death, IST4134, and clonal isolate IST439, was performed and compared with a prior analysis of IST4113 vs. IST439. The cytoplasmic and membrane-associated enriched fractions were examined and 52 proteins were found to be similarly altered in the two last isolates compared with IST439. These proteins are involved in metabolic functions, nucleotide synthesis, translation and protein folding, cell envelope biogenesis and iron homeostasis. Results are suggestive of the important role played by metabolic reprogramming in the virulence potential and persistence of B. cenocepacia, in particular regarding bacterial adaptation to microaerophilic conditions. Also, the content of the virulence determinant AidA was higher in the last 2 isolates. Significant levels of siderophores were found to be secreted by the three clonal isolates in an iron-depleted environment, but the two late isolates were more tolerant to low iron concentrations than IST439, consistent with the relative abundance of proteins involved in iron uptake.     

Molecular Screening of Enterococcus Virulence Determinants and Potential for Genetic Exchange between Food and Medical Isolates

Enterococci are used as starter and probiotic cultures in foods, and they occur as natural food contaminants. The genus Enterococcus is of increased significance as a cause of nosocomial infections, and this trend is exacerbated by the development of antibiotic resistance. In this study, we investigated the incidence of known virulence determinants in starter, food, and medical strains of Enterococcus faecalis, E. faecium, and E. durans. PCR and gene probe strategies were used to screen enterococcal isolates from both food and medical sources. Different and distinct patterns of incidence of virulence determinants were found for the E. faecalis and E. faecium strains. Medical E. faecalis strains had more virulence determinants than did food strains, which, in turn, had more than did starter strains. All of the E. faecalis strains tested possessed multiple determinants (between 6 and 11). E. faecium strains were generally free of virulence determinants, with notable exceptions. Significantly, esp and gelE determinants were identified in E. faecium medical strains. These virulence determinants have not previously been identified in E. faecium strains and may result from regional differences or the evolution of pathogenic E. faecium. Phenotypic testing revealed the existence of apparently silent gelE and cyl genes. In E. faecalis, the trend in these silent genes mirrors that of the expressed determinants. The potential for starter strains to acquire virulence determinants by natural conjugation mechanisms was investigated. Transconjugation in which starter strains acquired additional virulence determinants from medical strains was demonstrated. In addition, multiple pheromone-encoding genes were identified in both food and starter strains, indicating their potential to acquire other sex pheromone plasmids. These results suggest that the use of Enterococcus spp. in foods requires careful safety evaluation.     

Diversity of Cultivated Endophytic Bacteria from Sugarcane: Genetic and Biochemical Characterization of Burkholderia cepacia Complex Isolates

Bacteria were isolated from the rhizosphere and from inside the roots and stems of sugarcane plants grown in the field in Brazil. Endophytic bacteria were found in both the roots and the stems of sugarcane plants, with a significantly higher density in the roots. Many of the cultivated endophytic bacteria were shown to produce the plant growth hormone indoleacetic acid, and this trait was more frequently found among bacteria from the stem. 16S rRNA gene sequence analysis revealed that the selected isolates of the endophytic bacterial community of sugarcane belong to the genera of Burkholderia, Pantoea, Pseudomonas, and Microbacterium. Bacterial isolates belonging to the genus Burkholderia were the most predominant among the endophytic bacteria. Many of the Burkholderia isolates produced the antifungal metabolite pyrrolnitrin, and all were able to grow at 37°C. Phylogenetic analyses of the 16S rRNA gene and recA gene sequences indicated that the endophytic Burkholderia isolates from sugarcane are closely related to clinical isolates of the Burkholderia cepacia complex and clustered with B. cenocepacia (gv. III) isolates from cystic fibrosis patients. These results suggest that isolates of the B. cepacia complex are an integral part of the endophytic bacterial community of sugarcane in Brazil and reinforce the hypothesis that plant-associated environments may act as a niche for putative opportunistic human pathogenic bacteria.     

Leveraging structure determination with fragment screening for infectious disease drug targets: MECP synthase from Burkholderia pseudomallei

As part of the Seattle Structural Genomics Center for Infectious Disease, we seek to enhance structural genomics with ligand-bound structure data which can serve as a blueprint for structure-based drug design. We have adapted fragment-based screening methods to our structural genomics pipeline to generate multiple ligand-bound structures of high priority drug targets from pathogenic organisms. In this study, we report fragment screening methods and structure determination results for 2C-methyl-D-erythritol-2,4-cyclo-diphosphate (MECP) synthase from Burkholderia pseudomallei, the gram-negative bacterium which causes melioidosis. Screening by nuclear magnetic resonance spectroscopy as well as crystal soaking followed by X-ray diffraction led to the identification of several small molecules which bind this enzyme in a critical metabolic pathway. A series of complex structures obtained with screening hits reveal distinct binding pockets and a range of small molecules which form complexes with the target. Additional soaks with these compounds further demonstrate a subset of fragments to only bind the protein when present in specific combinations. This ensemble of fragment-bound complexes illuminates several characteristics of MECP synthase, including a previously unknown binding surface external to the catalytic active site. These ligand-bound structures now serve to guide medicinal chemists and structural biologists in rational design of novel inhibitors for this enzyme.