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.     

Genome-wide analysis reveals loci encoding anti-macrophage factors in the human pathogen Burkholderia pseudomallei K96243

Burkholderia pseudomallei is an important human pathogen whose infection biology is still poorly understood. The bacterium is endemic to tropical regions, including South East Asia and Northern Australia, where it causes melioidosis, a serious disease associated with both high mortality and antibiotic resistance. B. pseudomallei is a Gram-negative facultative intracellular pathogen that is able to replicate in macrophages. However despite the critical nature of its interaction with macrophages, few anti-macrophage factors have been characterized to date. Here we perform a genome-wide gain of function screen of B. pseudomallei strain K96243 to identify loci encoding factors with anti-macrophage activity. We identify a total of 113 such loci scattered across both chromosomes, with positive gene clusters encoding transporters and secretion systems, enzymes/toxins, secondary metabolite, biofilm, adhesion and signal response related factors. Further phenotypic analysis of four of these regions shows that the encoded factors cause striking cellular phenotypes relevant to infection biology, including apoptosis, formation of actin 'tails' and multi-nucleation within treated macrophages. The detailed analysis of the remaining host of loci will facilitate genetic dissection of the interaction of this important pathogen with host macrophages and thus further elucidate this critical part of its infection cycle.     

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.     

Genomic islands as a marker to differentiate between clinical and environmental Burkholderia pseudomallei

Burkholderia pseudomallei, as a saprophytic bacterium that can cause a severe sepsis disease named melioidosis, has preserved several extra genes in its genome for survival. The sequenced genome of the organism showed high diversity contributed mainly from genomic islands (GIs). Comparative genome hybridization (CGH) of 3 clinical and 2 environmental isolates, using whole genome microarrays based on B. pseudomallei K96243 genes, revealed a difference in the presence of genomic islands between clinical and environmental isolates. The largest GI, GI8, of B. pseudomallei was observed as a 2 sub-GI named GIs8.1 and 8.2 with distinguishable %GC content and unequal presence in the genome. GIs8.1, 8.2 and 15 were found to be more common in clinical isolates. A new GI, GI16c, was detected on chromosome 2. Presences of GIs8.1, 8.2, 15 and 16c were evaluated in 70 environmental and 64 clinical isolates using PCR assays. A combination of GIs8.1 and 16c (positivity of either GI) was detected in 70% of clinical isolates and 11.4% of environmental isolates (P<0.001). Using BALB/c mice model, no significant difference of time to mortality was observed between K96243 isolate and three isolates without GIs under evaluation (P>0.05). Some virulence genes located in the absent GIs and the difference of GIs seems to contribute less to bacterial virulence. The PCR detection of 2 GIs could be used as a cost effective and rapid tool to detect potentially virulent isolates that were contaminated in soil.     

Genomics of pyrrolnitrin biosynthetic loci: evidence for conservation and whole-operon mobility within gram-negative bacteria

Pyrrolnitrin (PRN) is a tryptophan-derived secondary metabolite produced by a narrow range of Gram-negative bacteria. The PRN biosynthesis by rhizobacteria presumably has a key role in their life strategies and in the biocontrol of plant diseases. The biosynthetic operon that encodes the pathway that converts tryptophan to PRN is composed of four genes, prnA through D , whose diversity, genomic context and spread over bacterial genomes are poorly understood. Therefore, we launched an endeavour aimed at retrieving, by in vitro and in silico means, diverse bacteria carrying the prnABCD biosynthetic loci in their genomes. Analysis of polymorphisms of the prnD gene sequences revealed a high level of conservation between Burkholderia , Pseudomonas and Serratia spp. derived sequences. Whole-operon- and prnD -based phylogeny resulted in tree topologies that are incongruent with the taxonomic status of the evaluated strains as predicted by 16S rRNA gene phylogeny. The genomic composition of c . 20 kb DNA fragments containg the PRN operon varied in different strains. Highly conserved and distinct transposase-encoding genes surrounding the PRN biosynthetic operons of Burkholderia pseudomallei strains were found. A prnABCD -deprived genomic region in B. pseudomallei strain K96243 contained the same gene composition as, and shared high homology with, the flanking regions of the PRN operon in B. pseudomallei strains 668, 1106a and 1710b. Our results strongly suggest that the PRN biosynthetic operon is mobile. The extent, frequency and promiscuity of this mobility remain to be understood.     

Alanine racemase mutants of Burkholderia pseudomallei and Burkholderia mallei and use of alanine racemase as a non-antibiotic-based selectable marker

Burkholderia pseudomallei and Burkholderia mallei are category B select agents and must be studied under BSL3 containment in the United States. They are typically resistant to multiple antibiotics, and the antibiotics used to treat B. pseudomallei or B. mallei infections may not be used as selective agents with the corresponding Burkholderia species. Here, we investigated alanine racemase deficient mutants of B. pseudomallei and B. mallei for development of non-antibiotic-based genetic selection methods and for attenuation of virulence. The genome of B. pseudomallei K96243 has two annotated alanine racemase genes (bpsl2179 and bpss0711), and B. mallei ATCC 23344 has one (bma1575). Each of these genes encodes a functional enzyme that can complement the alanine racemase deficiency of Escherichia coli strain ALA1. Herein, we show that B. pseudomallei with in-frame deletions in both bpsl2179 and bpss0711, or B. mallei with an in-frame deletion in bma1575, requires exogenous D-alanine for growth. Introduction of bpsl2179 on a multicopy plasmid into alanine racemase deficient variants of either Burkholderia species eliminated the requirement for D-alanine. During log phase growth without D-alanine, the viable counts of alanine racemase deficient mutants of B. pseudomallei and B. mallei decreased within 2 hours by about 1000-fold and 10-fold, respectively, and no viable bacteria were present at 24 hours. We constructed several genetic tools with bpsl2179 as a selectable genetic marker, and we used them without any antibiotic selection to construct an in-frame ΔflgK mutant in the alanine racemase deficient variant of B. pseudomallei K96243. In murine peritoneal macrophages, wild type B. mallei ATCC 23344 was killed much more rapidly than wild type B. pseudomallei K96243. In addition, the alanine racemase deficient mutant of B. pseudomallei K96243 exhibited attenuation versus its isogenic parental strain with respect to growth and survival in murine peritoneal macrophages.     

prhKLM genes of Ralstonia solanacearum encode novel activators of hrp regulon and are required for pathogenesis in tomato

Burkholderia pseudomallei, the causative agent of melioidosis, exploits the Bsa type III secretion system (T3SS) to deliver effector proteins into host cells. These effectors manipulate host cell functions; thus, contributing to the ability of the bacteria to evade the immune response and cause disease. Only two Bsa-secreted effectors have been conclusively identified to date. Here, we report the identification of the third B. pseudomallei type III secreted effector protein, designated BopC. BopC is encoded by the bpss1516 gene abutting bpss1517, which encodes its putative chaperone. The genes are located in the close proximity to the bsa T3SS gene cluster of B. pseudomallei K96243 (Fig. 1). BopC was secreted into culture supernatant by the wild-type B. pseudomallei strain, but its secretion was abolished in the bsaZ T3SS mutant. Using pull down and co-purification assays, we confirmed that BopC interacts with its putative chaperone, BPSS1517, in vitro. Furthermore, the first 20 N-terminal amino acids of BopC were found to be sufficient to mediate the T3SS-dependent translocation of a reporter protein from a heterologous enteropathogenic Escherichia coli host into mammalian cells. Finally, bopC mutant was found to be less invasive than the wild-type strain in the epithelial cells.     

Exome sequencing and arrayCGH detection of gene sequence and copy number variation between ILS and ISS mouse strains

It has been well documented that genetic factors can influence predisposition to develop alcoholism. While the underlying genomic changes may be of several types, two of the most common and disease associated are copy number variations (CNVs) and sequence alterations of protein coding regions. The goal of this study was to identify CNVs and single-nucleotide polymorphisms that occur in gene coding regions that may play a role in influencing the risk of an individual developing alcoholism. Toward this end, two mouse strains were used that have been selectively bred based on their differential sensitivity to alcohol: the Inbred long sleep (ILS) and Inbred short sleep (ISS) mouse strains. Differences in initial response to alcohol have been linked to risk for alcoholism, and the ILS/ISS strains are used to investigate the genetics of initial sensitivity to alcohol. Array comparative genomic hybridization (arrayCGH) and exome sequencing were conducted to identify CNVs and gene coding sequence differences, respectively, between ILS and ISS mice. Mouse arrayCGH was performed using catalog Agilent 1 × 244 k mouse arrays. Subsequently, exome sequencing was carried out using an Illumina HiSeq 2000 instrument. ArrayCGH detected 74 CNVs that were strain-specific (38 ILS/36 ISS), including several ISS-specific deletions that contained genes implicated in brain function and neurotransmitter release. Among several interesting coding variations detected by exome sequencing was the gain of a premature stop codon in the alpha-amylase 2B (AMY2B) gene specifically in the ILS strain. In total, exome sequencing detected 2,597 and 1,768 strain-specific exonic gene variants in the ILS and ISS mice, respectively. This study represents the most comprehensive and detailed genomic comparison of ILS and ISS mouse strains to date. The two complementary genome-wide approaches identified strain-specific CNVs and gene coding sequence variations that should provide strong candidates to contribute to the alcohol-related phenotypic differences associated with these strains.     

Strain-dependent regulation of neurotransmission and actin-remodelling proteins in the mouse hippocampus

Individual mouse strains differ significantly in terms of behaviour, cognitive function and long-term potentiation. Hippocampal gene expression profiling of eight different mouse strains points towards strain-specific regulation of genes involved in neuronal information storage. Protein expression with regard to strain- dependent expression of structures related to neuronal information storage has not been investigated yet. Herein, a proteomic approach based on two-dimensional gel electrophoresis coupled with mass spectrometry (MALDI-TOF/TOF) has been chosen to address this question by determining strain-dependent expression of proteins involved in neurotransmission and activity-induced actin remodelling in hippocampal tissue of five mouse strains. Of 31 spots representing 16 different gene products analysed and quantified, N-ethylmaleimide-sensitive fusion protein, N-ethylmaleimide-sensitive factor attachment protein-alpha, actin-like protein 3, profilin and cofilin were expressed in a strain-dependent manner. By treating protein expression as a phenotype, we have shown significant genetic variation in brain protein expression. Further experiments in this direction may provide an indication of the genetic elements that contribute to the phenotypic differences between the selected strains through the expressional level of the translated protein. In view of this, we propose that proteomic analysis enabling to concomitantly survey the expression of a large number of proteins could serve as a valuable tool for genetic and physiological studies of central nervous system function.     

A horizontal gene transfer event defines two distinct groups within Burkholderia pseudomallei that have dissimilar geographic distributions

Burkholderia pseudomallei is the etiologic agent of melioidosis. Many disease manifestations are associated with melioidosis, and the mechanisms causing this variation are unknown; genomic differences among strains offer one explanation. We compared the genome sequences of two strains of B. pseudomallei: the original reference strain K96243 from Thailand and strain MSHR305 from Australia. We identified a variable homologous region between the two strains. This region was previously identified in comparisons of the genome of B. pseudomallei strain K96243 with the genome of strain E264 from the closely related B. thailandensis. In that comparison, K96243 was shown to possess a horizontally acquired Yersinia-like fimbrial (YLF) gene cluster. Here, we show that the homologous genomic region in B. pseudomallei strain 305 is similar to that previously identified in B. thailandensis strain E264. We have named this region in B. pseudomallei strain 305 the B. thailandensis-like flagellum and chemotaxis (BTFC) gene cluster. We screened for these different genomic components across additional genome sequences and 571 B. pseudomallei DNA extracts obtained from regions of endemicity. These alternate genomic states define two distinct groups within B. pseudomallei: all strains contained either the BTFC gene cluster (group BTFC) or the YLF gene cluster (group YLF). These two groups have distinct geographic distributions: group BTFC is dominant in Australia, and group YLF is dominant in Thailand and elsewhere. In addition, clinical isolates are more likely to belong to group YLF, whereas environmental isolates are more likely to belong to group BTFC. These groups should be further characterized in an animal model.     

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.