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.     

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.     

Genomic Characterization of Burkholderia pseudomallei Isolates Selected for Medical Countermeasures Testing: Comparative Genomics Associated with Differential Virulence

Burkholderia pseudomallei is the causative agent of melioidosis and a potential bioterrorism agent. In the development of medical countermeasures against B. pseudomallei infection, the US Food and Drug Administration (FDA) animal Rule recommends using well-characterized strains in animal challenge studies. In this study, whole genome sequence data were generated for 6 B. pseudomallei isolates previously identified as candidates for animal challenge studies; an additional 5 isolates were sequenced that were associated with human inhalational melioidosis. A core genome single nucleotide polymorphism (SNP) phylogeny inferred from a concatenated SNP alignment from the 11 isolates sequenced in this study and a diverse global collection of isolates demonstrated the diversity of the proposed Animal Rule isolates. To understand the genomic composition of each isolate, a large-scale blast score ratio (LS-BSR) analysis was performed on the entire pan-genome; this demonstrated the variable composition of genes across the panel and also helped to identify genes unique to individual isolates. In addition, a set of ~550 genes associated with pathogenesis in B. pseudomallei were screened against the 11 sequenced genomes with LS-BSR. Differential gene distribution for 54 virulence-associated genes was observed between genomes and three of these genes were correlated with differential virulence observed in animal challenge studies using BALB/c mice. Differentially conserved genes and SNPs associated with disease severity were identified and could be the basis for future studies investigating the pathogenesis of B. pseudomallei. Overall, the genetic characterization of the 11 proposed Animal Rule isolates provides context for future studies involving B. pseudomallei pathogenesis, differential virulence, and efficacy to therapeutics.     

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.     

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 island identification in Vibrio vulnificus reveals significant genome plasticity in this human pathogen

Genomic islands (GIs) are large chromosomal regions present in a subset of bacterial strains that increase the fitness of the organism under specific conditions. We compared the complete genome sequences of two Vibrio vulnificus strains YJ016 and CMCP6 and identified 14 regions (ranging in size from 14 to 117 kb), which had the characteristics of GIs. Bioinformatic analysis of these 14 GI regions identified the presence of phage-like integrase genes, aberrant GC content and genome signature (dinucleotide frequency) within each GI compared with the core genome indicating that these regions were acquired from an anomalous source. We examined the distribution of the nine GIs from strain YJ016 among 27 V. vulnificus isolates and found that most GIs were absent from the majority of these isolates. The chromosomal insertion sites of three GIs were adjacent to tRNA sites, which contained novel horizontally acquired DNA in all six available sequenced Vibrionaceae genomes. Supplementary information: Supplementary data are available at Bioinformatics online.     

Comparative virulence of three different strains of Burkholderia pseudomallei in an aerosol non-human primate model

Melioidosis, caused by the Gram-negative bacterium Burkholderia pseudomallei, is a major cause of sepsis and mortality in endemic regions of Southeast Asia and Northern Australia. B. pseudomallei is a potential bioterrorism agent due to its high infectivity, especially via inhalation, and its inherent resistance to antimicrobials. There is currently no vaccine for melioidosis and antibiotic treatment can fail due to innate drug resistance, delayed diagnosis and treatment, or insufficient duration of treatment. A well-characterized animal model that mimics human melioidosis is needed for the development of new medical countermeasures. This study first characterized the disease progression of melioidosis in the African green monkey (AGM) and rhesus macaque (RM) for non-human primate model down-selection. All AGMs developed acute lethal disease similar to that described in human acute infection following exposure to aerosolized B. pseudomallei strain HBPUB10134a. Only 20% of RMs succumbed to acute disease. Disease progression, immune response and pathology of two other strains of B. pseudomallei, K96243 and MSHR5855, were also compared using AGMs. These three B. pseudomallei strains represent a highly virulent strain from Thailand (HBPUB101034a), a highly virulent strains from Australia (MSHR5855), and a commonly used laboratory strains originating from Thailand (K96243). Animals were observed for clinical signs of infection and blood samples were analyzed for cytokine responses, blood chemistry and leukocyte changes in order to characterize bacterial infection. AGMs experienced fever after exposure to aerosolized B. pseudomallei at the onset of acute disease. Inflammation, abscesses and/or pyogranulomas were observed in lung with all three strains of B. pseudomallei. Inflammation, abscesses and/or pyogranulomas were observed in lymph nodes, spleen, liver and/or kidney with B. pseudomallei, HBPUB10134a and K96243. Additionally, the Australian strain MSHR5855 induced brain lesions in one AGM similar to clinical cases of melioidosis seen in Australia. Elevated serum levels of IL-1β, IL-1 receptor antagonist, IL-6, MCP-1, G-CSF, HGF, IFNγ, MIG, I-TAC, and MIP-1β at terminal end points can be significantly correlated with non-survivors with B. pseudomallei infection in AGM. The AGM model represents an acute model of B. pseudomallei infection for all three strains from two geographical locations and will be useful for efficacy testing of vaccines and therapeutics against melioidosis. In summary, a dysregulated immune response leading to excessive persistent inflammation and inflammatory cell death is the key driver of acute melioidosis. Early intervention in these pathways will be necessary to counter B. pseudomallei and mitigate the pathological consequences of melioidosis.     

Exploring the zoonotic potential of Mycobacterium avium subspecies paratuberculosis through comparative genomics

A comparative genomics approach was utilised to compare the genomes of Mycobacterium avium subspecies paratuberculosis (MAP) isolated from early onset paediatric Crohn's disease (CD) patients as well as Johne's diseased animals. Draft genome sequences were produced for MAP isolates derived from four CD patients, one ulcerative colitis (UC) patient, and two non-inflammatory bowel disease (IBD) control individuals using Illumina sequencing, complemented by comparative genome hybridisation (CGH). MAP isolates derived from two bovine and one ovine host were also subjected to whole genome sequencing and CGH. All seven human derived MAP isolates were highly genetically similar and clustered together with one bovine type isolate following phylogenetic analysis. Three other sequenced isolates (including the reference bovine derived isolate K10) were genetically distinct. The human isolates contained two large tandem duplications, the organisations of which were confirmed by PCR. Designated vGI-17 and vGI-18 these duplications spanned 63 and 109 open reading frames, respectively. PCR screening of over 30 additional MAP isolates (3 human derived, 27 animal derived and one environmental isolate) confirmed that vGI-17 and vGI-18 are common across many isolates. Quantitative real-time PCR of vGI-17 demonstrated that the proportion of cells containing the vGI-17 duplication varied between 0.01 to 15% amongst isolates with human isolates containing a higher proportion of vGI-17 compared to most animal isolates. These findings suggest these duplications are transient genomic rearrangements. We hypothesise that the over-representation of vGI-17 in human derived MAP strains may enhance their ability to infect or persist within a human host by increasing genome redundancy and conferring crude regulation of protein expression across biologically important regions.     

Genome reduction by deletion of paralogs in the marine cyanobacterium Prochlorococcus

Several isolates of the marine cyanobacterial genus Prochlorococcus have smaller genome sizes than those of the closely related genus Synechococcus. In order to test whether loss of protein-coding genes has contributed to genome size reduction in Prochlorococcus, we reconstructed events of gene family evolution over a strongly supported phylogeny of 12 Prochlorococcus genomes and 9 Synechococcus genomes. Significantly, more events both of loss of paralogs within gene families and of loss of entire gene families occurred in Prochlorococcus than in Synechococcus. The number of nonancestral gene families in genomes of both genera was positively correlated with the extent of genomic islands (GIs), consistent with the hypothesis that horizontal gene transfer (HGT) is associated with GIs. However, even when only isolates with comparable extents of GIs were compared, significantly more events of gene family loss and of paralog loss were seen in Prochlorococcus than in Synechococcus, implying that HGT is not the primary reason for the genome size difference between the two genera.     

Genomic variations define divergence of water/wildlife-associated Campylobacter jejuni niche specialists from common clonal complexes

Although the major food-borne pathogen Campylobacter jejuni has been isolated from diverse animal, human and environmental sources, our knowledge of genomic diversity in C. jejuni is based exclusively on human or human food-chain-associated isolates. Studies employing multilocus sequence typing have indicated that some clonal complexes are more commonly associated with particular sources. Using comparative genomic hybridization on a collection of 80 isolates representing diverse sources and clonal complexes, we identified a separate clade comprising a group of water/wildlife isolates of C. jejuni with multilocus sequence types uncharacteristic of human food-chain-associated isolates. By genome sequencing one representative of this diverse group (C. jejuni 1336), and a representative of the bank-vole niche specialist ST-3704 (C. jejuni 414), we identified deletions of genomic regions normally carried by human food-chain-associated C. jejuni. Several of the deleted regions included genes implicated in chicken colonization or in virulence. Novel genomic insertions contributing to the accessory genomes of strains 1336 and 414 were identified. Comparative analysis using PCR assays indicated that novel regions were common but not ubiquitous among the water/wildlife group of isolates, indicating further genomic diversity among this group, whereas all ST-3704 isolates carried the same novel accessory regions. While strain 1336 was able to colonize chicks, strain 414 was not, suggesting that regions specifically absent from the genome of strain 414 may play an important role in this common route of Campylobacter infection of humans. We suggest that the genomic divergence observed constitutes evidence of adaptation leading to niche specialization.     

Comparative Whole-Genome Analysis of Clinical Isolates Reveals Characteristic Architecture of Mycobacterium tuberculosis Pangenome

The tubercle complex consists of closely related mycobacterium species which appear to be variants of a single species. Comparative genome analysis of different strains could provide useful clues and insights into the genetic diversity of the species. We integrated genome assemblies of 96 strains from Mycobacterium tuberculosis complex (MTBC), which included 8 Indian clinical isolates sequenced and assembled in this study, to understand its pangenome architecture. We predicted genes for all the 96 strains and clustered their respective CDSs into homologous gene clusters (HGCs) to reveal a hard-core, soft-core and accessory genome component of MTBC. The hard-core (HGCs shared amongst 100% of the strains) was comprised of 2,066 gene clusters whereas the soft-core (HGCs shared amongst at least 95% of the strains) comprised of 3,374 gene clusters. The change in the core and accessory genome components when observed as a function of their size revealed that MTBC has an open pangenome. We identified 74 HGCs that were absent from reference strains H37Rv and H37Ra but were present in most of clinical isolates. We report PCR validation on 9 candidate genes depicting 7 genes completely absent from H37Rv and H37Ra whereas 2 genes shared partial homology with them accounting to probable insertion and deletion events. The pangenome approach is a promising tool for studying strain specific genetic differences occurring within species. We also suggest that since selecting appropriate target genes for typing purposes requires the expected target gene be present in all isolates being typed, therefore estimating the core-component of the species becomes a subject of prime importance.     

Genomic Diversity of Escherichia Isolates from Diverse Habitats

Our understanding of the Escherichia genus is heavily biased toward pathogenic or commensal isolates from human or animal hosts. Recent studies have recovered Escherichia isolates that persist, and even grow, outside these hosts. Although the environmental isolates are typically phylogenetically distinct, they are highly related to and phenotypically indistinguishable from their human counterparts, including for the coliform test. To gain insights into the genomic diversity of Escherichia isolates from diverse habitats, including freshwater, soil, animal, and human sources, we carried out comparative DNA-DNA hybridizations using a multi-genome E. coli DNA microarray. The microarray was validated based on hybridizations with selected strains whose genome sequences were available and used to assess the frequency of microarray false positive and negative signals. Our results showed that human fecal isolates share two sets of genes (n>90) that are rarely found among environmental isolates, including genes presumably important for evading host immune mechanisms (e.g., a multi-drug transporter for acids and antimicrobials) and adhering to epithelial cells (e.g., hemolysin E and fimbrial-like adhesin protein). These results imply that environmental isolates are characterized by decreased ability to colonize host cells relative to human isolates. Our study also provides gene markers that can distinguish human isolates from those of warm-blooded animal and environmental origins, and thus can be used to more reliably assess fecal contamination in natural ecosystems.     

Acquired type III secretion system determines environmental fitness of epidemic Vibrio parahaemolyticus in the interaction with bacterivorous protists

Genome analyses of marine microbial communities have revealed the widespread occurrence of genomic islands (GIs), many of which encode for protein secretion machineries described in the context of bacteria-eukaryote interactions. Yet experimental support for the specific roles of such GIs in aquatic community interactions remains scarce. Here, we test for the contribution of type III secretion systems (T3SS) to the environmental fitness of epidemic Vibrio parahaemolyticus. Comparisons of V. parahaemolyticus wild types and T3SS-defective mutants demonstrate that the T3SS encoded on genome island VPaI-7 (T3SS-2) promotes survival of V. parahaemolyticus in the interaction with diverse protist taxa. Enhanced persistence was found to be due to T3SS-2 mediated cytotoxicity and facultative parasitism of V. parahaemolyticus on coexisting protists. Growth in the presence of bacterivorous protists and the T3SS-2 genotype showed a strong correlation across environmental and clinical isolates of V. parahaemolyticus. Short-term microcosm experiments provide evidence that protistan hosts facilitate the invasion of T3SS-2 positive V. parahaemolyticus into a coastal plankton community, and that water temperature and productivity further promote enhanced survival of T3SS-2 positive V. parahaemolyticus. This study is the first to describe the fitness advantage of GI-encoded functions in a microbial food web, which may provide a mechanistic explanation for the global spread and the seasonal dynamics of V. parahaemolyticus pathotypes, including the pandemic serotype cluster O3:K6, in aquatic environments.     

Identification of genes and genomic islands correlated with high pathogenicity in Streptococcus suis using whole genome tiling microarrays

Streptococcus suis is an important zoonotic pathogen that can cause meningitis and sepsis in both pigs and humans. Infections in humans have been sporadic worldwide but two severe outbreaks occurred in China in recent years, while infections in pigs are a major problem in the swine industry. Some S. suis strains are more pathogenic than others with 2 sequence types (ST), ST1 and ST7, being well recognized as highly pathogenic. We analyzed 31 isolates from 23 serotypes and 25 STs by NimbleGen tiling microarray using the genome of a high pathogenicity (HP) ST1 strain, GZ1, as reference and a new algorithm to detect gene content difference. The number of genes absent in a strain ranged from 49 to 225 with a total of 632 genes absent in at least one strain, while 1346 genes were found to be invariably present in all strains as the core genome of S. suis, accounting for 68% of the GZ1 genome. The majority of genes are located in chromosomal blocks with two or more contiguous genes. Sixty two blocks are absent in two or more strains and defined as regions of difference (RDs), among which 26 are putative genomic islands (GIs). Clustering and statistical analyses revealed that 8 RDs including 6 putative GIs and 21 genes within these RDs are significantly associated with HP. Three RDs encode known virulence related factors including the extracellular factor, the capsular polysaccharide and a SrtF pilus. The strains were divided into 5 groups based on population genetic analysis of multilocus sequence typing data and the distribution of the RDs among the groups revealed gain and loss of RDs in different groups. Our study elucidated the gene content diversity of S. suis and identified genes that potentially promote HP.     

What It Takes to Be a Pseudomonas aeruginosa? The Core Genome of the Opportunistic Pathogen Updated

Pseudomonas aeruginosa is an opportunistic bacterial pathogen able to thrive in highly diverse ecological niches and to infect compromised patients. Its genome exhibits a mosaic structure composed of a core genome into which accessory genes are inserted en bloc at specific sites. The size and the content of the core genome are open for debate as their estimation depends on the set of genomes considered and the pipeline of gene detection and clustering. Here, we redefined the size and the content of the core genome of P. aeruginosa from fully re-analyzed genomes of 17 reference strains. After the optimization of gene detection and clustering parameters, the core genome was defined at 5,233 orthologs, which represented ~ 88% of the average genome. Extrapolation indicated that our panel was suitable to estimate the core genome that will remain constant even if new genomes are added. The core genome contained resistance determinants to the major antibiotic families as well as most metabolic, respiratory, and virulence genes. Although some virulence genes were accessory, they often related to conserved biological functions. Long-standing prophage elements were subjected to a genetic drift to eventually display a G+C content as higher as that of the core genome. This contrasts with the low G+C content of highly conserved ribosomal genes. The conservation of metabolic and respiratory genes could guarantee the ability of the species to thrive on a variety of carbon sources for energy in aerobiosis and anaerobiosis. Virtually all the strains, of environmental or clinical origin, have the complete toolkit to become resistant to the major antipseudomonal compounds and possess basic pathogenic mechanisms to infect humans. The knowledge of the genes shared by the majority of the P. aeruginosa isolates is a prerequisite for designing effective therapeutics to combat the wide variety of human infections.     

Identification of Genes and Genomic Islands Correlated with High Pathogenicity in Streptococcus suis Using Whole Genome Tilling Microarrays

Streptococcus suis is an important zoonotic pathogen that can cause meningitis and sepsis in both pigs and humans. Infections in humans have been sporadic worldwide but two severe outbreaks occurred in China in recent years, while infections in pigs are a major problem in the swine industry. Some S. suis strains are more pathogenic than others with 2 sequence types (ST), ST1 and ST7, being well recognized as highly pathogenic. We analyzed 31 isolates from 23 serotypes and 25 STs by NimbleGen tiling microarray using the genome of a high pathogenicity (HP) ST1 strain, GZ1, as reference and a new algorithm to detect gene content difference. The number of genes absent in a strain ranged from 49 to 225 with a total of 632 genes absent in at least one strain, while 1346 genes were found to be invariably present in all strains as the core genome of S. suis, accounting for 68% of the GZ1 genome. The majority of genes are located in chromosomal blocks with two or more contiguous genes. Sixty two blocks are absent in two or more strains and defined as regions of difference (RDs), among which 26 are putative genomic islands (GIs). Clustering and statistical analyses revealed that 8 RDs including 6 putative GIs and 21 genes within these RDs are significantly associated with HP. Three RDs encode known virulence related factors including the extracellular factor, the capsular polysaccharide and a SrtF pilus. The strains were divided into 5 groups based on population genetic analysis of multilocus sequence typing data and the distribution of the RDs among the groups revealed gain and loss of RDs in different groups. Our study elucidated the gene content diversity of S. suis and identified genes that potentially promote HP.