Identification and Characterization of Putative Virulence Genes and Gene Clusters in Aeromonas hydrophila PPD134/91

Aeromonas hydrophila is a gram-negative opportunistic pathogen of animals and humans. The pathogenesis of A. hydrophila is multifactorial. Genomic subtraction and markers of genomic islands (GIs) were used to identify putative virulence genes in A. hydrophila PPD134/91. Two rounds of genomic subtraction led to the identification of 22 unique DNA fragments encoding 19 putative virulence factors and seven new open reading frames, which are commonly present in the eight virulence strains examined. In addition, four GIs were found, including O-antigen, capsule, phage-associated, and type III secretion system (TTSS) gene clusters. These putative virulence genes and gene clusters were positioned on a physical map of A. hydrophila PPD134/91 to determine their genetic organization in this bacterium. Further in vivo study of insertion and deletion mutants showed that the TTSS may be one of the important virulence factors in A. hydrophila pathogenesis. Furthermore, deletions of multiple virulence factors such as S-layer, serine protease, and metalloprotease also increased the 50% lethal dose to the same level as the TTSS mutation (about 1 log) in a blue gourami infection model. This observation sheds light on the multifactorial and concerted nature of pathogenicity in A. hydrophila. The large number of putative virulence genes identified in this study will form the basis for further investigation of this emerging pathogen and help to develop effective vaccines, diagnostics, and novel therapeutics.     

Use of proteomics to identify novel virulence determinants that are required for Edwardsiella tarda pathogenesis

Edwardsiella tarda is an important cause of haemorrhagic septicaemia in fish and also of gastro- and extraintestinal infections in humans. Using a combination of comparative proteomics and TnphoA mutagenesis, we have identified five proteins that may contribute to E. tarda PPD130/91 pathogenesis. Lowered protein secretion, impaired autoaggregation and the absence of six proteins were observed only in three highly attenuated mutants when cultured in Dulbecco's modified eagle medium (DMEM). Five out of six proteins could be identified by their mass spectra. Three proteins were identified as putative effector proteins (EseB, EseC and EseD) that are homologous to SseB, SseC and SseD of a type III secretion system (TTSS) in Salmonella species. The other two were EvpA and EvpC, homologous to Eip20 and Eip18 in Edwardsiella ictaluri. The complete sequencing and homology studies of evpA-H indicate that similar gene clusters are widely distributed in other pathogens such as Escherichia, Salmonella, Vibrio and Yersinia species with unknown functions. Insertional inactivation and deletion of evpB or evpC led to lower replication rates in gourami phagocytes, and reduced protein secretion and virulence in blue gourami. Complementation of these deletion mutants showed partial recovery in the above three phenotypes, indicating that these genes are vital for E. tarda pathogenesis. The transport of the EvpC protein may not use the TTSS in E. tarda. The expression of EvpA and EvpC as well as EseB, EseC and EseD was temperature dependent (suppressed at 37 degrees C), and disruption of esrB affected their expression. The present study identifies two possible secretion systems (TTSS and Evp) that are vital for E. tarda pathogenesis.     

Complete type III secretion system of a mesophilic Aeromonas hydrophila strain

We have investigated the existence and genetic organization of a functional type III secretion system (TTSS) in a mesophilic Aeromonas strain by initially using the Aeromonas hydrophila strain AH-3. We report for the first time the complete TTSS DNA sequence of an Aeromonas strain that comprises 35 genes organized in a similar disposition as that in Pseudomonas aeruginosa. Using several gene probes, we also determined the presence of a TTSS in clinical or environmental strains of different Aeromonas species: A. hydrophila, A. veronii, and A. caviae. By using one of the TTSS genes (ascV), we were able to obtain a defined insertion mutant in strain AH-3 (AH-3AscV), which showed reduced toxicity and virulence in comparison with the wild-type strain. Complementation of the mutant strain with a plasmid vector carrying ascV was fully able to restore the wild-type toxicity and virulence.     

Type III secretion system cluster 3 is required for maximal virulence of Burkholderia pseudomallei in a hamster infection model

Burkholderia pseudomallei, the etiological agent of melioidosis, is an animal pathogen capable of inducing a highly fatal septicemia. B. pseudomallei possesses three type III secretion system (TTSS) clusters, two of which (TTSS1 and TTSS2) are homologous to the TTSS of the plant pathogen Ralstonia solanacearum, and one (TTSS3) is homologous to the Salmonella SPI-1 mammalian pathogenicity island. We have demonstrated that TTSS3 is required for the full virulence of B. pseudomallei in a hamster model of infection. We have also examined the virulence of B. pseudomallei mutants deficient in several putative TTSS3 effector molecules, and found no significant attenuation of B. pseudomallei virulence in the hamster model.     

Detecting anomalous gene clusters and pathogenicity islands in diverse bacterial genomes.

A gene in a genome is defined as putative alien (pA) if its codon usage difference from the average gene exceeds a high threshold and codon usage differences from ribosomal protein genes, chaperone genes and protein-synthesis-processing factors are also high. pA gene clusters in bacterial genomes are relevant for detecting genomic islands (GIs), including pathogenicity islands (PAIs). Four other analyses appropriate to this task are G+C genome variation (the standard method); genomic signature divergences (dinucleotide bias); extremes of codon bias; and anomalies of amino acid usage. For example, the cagA domain of Helicobacter pylori is highly deviant in its genome signature and codon bias from the rest of the genome. Using these methods we can detect two potential PAIs in the Neisseria meningitidis genome, which contain hemagglutinin and/or hemolysin-related genes. Additionally, G+C variation and genome signature differences of the Mycobacterium tuberculosis genome indicate two pA gene clusters.     

Characterization of extracellular proteins produced by Aeromonas hydrophila AH-1

Aeromonas hydrophila is a ubiquitous Gram-negative bacterium which can cause motile aeromonad septicemia in both fish and humans. A. hydrophila secretes many extracellular proteins associated with pathogenicity and environmental adaptability. In this study, an extracellular proteome map of A. hydrophila AH-1 was constructed. The major extracellular virulence factors were characterized by comparing the proteomes of various deletion mutants with that of the wild type. The results suggested that serine protease was involved in the processing of a toxin and secreted enzymes such as hemolysin, glycerophospholipid-cholesterol acyltransferase and metalloprotease. We also showed that expressions of polar and lateral flagellins were under the control of temperature, FlhA, LafK, and RpoN. In addition, three novel proteins (potential effector proteins including one ExoT-like protein) were revealed to be secreted via the type III secretion system (TTSS) of A. hydrophila AH-1. Another novel finding was the demonstration of a crosstalk between the lateral flagellar system and the TTSS in A. hydrophila. These results showed that proteomics is a powerful tool for characterizing virulence factors. The construction of proteome maps will provide a valuable means of finding potential candidates for developing suitable diagnostics and therapeutics for this emerging pathogen.     

The Association of Virulence Factors with Genomic Islands

Background

It has been noted that many bacterial virulence factor genes are located within genomic islands (GIs; clusters of genes in a prokaryotic genome of probable horizontal origin). However, such studies have been limited to single genera or isolated observations. We have performed the first large-scale analysis of multiple diverse pathogens to examine this association. We additionally identified genes found predominantly in pathogens, but not non-pathogens, across multiple genera using 631 complete bacterial genomes, and we identified common trends in virulence for genes in GIs. Furthermore, we examined the relationship between GIs and clustered regularly interspaced palindromic repeats (CRISPRs) proposed to confer resistance to phage.

Methodology/Principal Findings

We show quantitatively that GIs disproportionately contain more virulence factors than the rest of a given genome (p<1E-40 using three GI datasets) and that CRISPRs are also over-represented in GIs. Virulence factors in GIs and pathogen-associated virulence factors are enriched for proteins having more “offensive” functions, e.g. active invasion of the host, and are disproportionately components of type III/IV secretion systems or toxins. Numerous hypothetical pathogen-associated genes were identified, meriting further study.

Conclusions/Significance

This is the first systematic analysis across diverse genera indicating that virulence factors are disproportionately associated with GIs. “Offensive” virulence factors, as opposed to host-interaction factors, may more often be a recently acquired trait (on an evolutionary time scale detected by GI analysis). Newly identified pathogen-associated genes warrant further study. We discuss the implications of these results, which cement the significant role of GIs in the evolution of many pathogens.     

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.     

The Erwinia chrysanthemi EC16 hrp/hrc gene cluster encodes an active Hrp type III secretion system that is flanked by virulence genes functionally unrelated to the Hrp system

Erwinia chrysanthemi is a host-promiscuous plant pathogen that possesses a type III secretion system (TTSS) similar to that of the host-specific pathogens E. amylovora and Pseudomonas syringae. The regions flanking the TTSS-encoding hrp/hrc gene clusters in the latter pathogens encode various TTSS-secreted proteins. DNA sequencing of the complete E. chrysanthemi hrp/hrc gene cluster and approximately 12 kb of the flanking regions (beyond the previously characterized hecA adhesin gene in the left flank) revealed that the E. chrysanthemi TTSS genes were syntenic and similar (>50% amino-acid identity) with their E. amylovora orthologs. However, the hrp/hrc cluster was interrupted by a cluster of four genes, only one of which, a homolog of lytic transglycosylases, is implicated in TTSS functions. Furthermore, the regions flanking the hrp/hrc cluster lacked genes that were likely to encode TTSS substrates. Instead, some of the genes in these regions predict ABC transporters and methyl-accepting chemotaxis proteins that could have alternative roles in virulence. Mutations affecting all of the genes in the regions flanking or interrupting the hrp/hrc cluster were constructed in E. chrysanthemi CUCPB5047, a mutant whose reduced pectolytic capacity can enhance the phenotype of minor virulence factors. Mutants were screened in witloof chicory leaves and then in potato tubers and Nicotiana clevelandii seedlings. Mu dII1734 insertion in one gene, designated virA, resulted in strongly reduced virulence in all three tests. virA is immediately downstream of hecA, has an unusually low G+C content of 38%, and predicts an unknown protein of 111 amino acids. The E. chrysanthemi TTSS was shown to be active by its ability to translocate AvrPto-Cya (a P. syringae TTSS effector fused to an adenylate cyclase reporter that is active in the presence of eukaryote calmodulin) into N. benthamiana leaf cells. However, VirA(1-61)-Cya was not translocated into plant cells, and virA expression was not affected by mutations in E. chrysanthemi Hrp regulator genes hrpL and hrpS. Thus, the 44-kb region of the E. chrysanthemi EC16 genome that is centered on the hrplhrc cluster encodes a potpourri of virulence factors, but none of these appear to be a TTSS effector.     

Comparative genomics analysis of statistically significant genomic islands of Helicobacter pylori strains for better understanding the disease prognosis

Bacterial virulence factors are often located in their genomic islands (GIs). Helicobacter pylori, a highly diverse organism is reported to be associated with several gastrointestinal diseases like, gastritis, gastric cancer (GC), peptic ulcer, duodenal ulcer (DU) etc. A novel similarity score (Sm)-based comparative analysis with GIs of 50 H. pylori strains revealed clear idea of the various factors which promote disease progression. Two putative pathogenic GIs in some of the H. pylori strains were identified. One GI, having a putative labile enterotoxin and other dynamin-like proteins (DLPs), is predicted to increase the release of toxin by membrane vesicular formation. Another island contains a virulence-associated protein D (vapD) which is a component of a type-II toxin–antitoxin system (TAs), leads to enhance the severity of the H. pylori infection. Besides the well-known virulence factors like Cytotoxin-associated gene A (CagA) and vacA, several GIs have been identified which showed to have direct or indirect impact on H. pylori clinical outcomes. One such GI, containing lipopolysaccharide (LPS) biosynthesis genes was revealed to be directly connected with disease development by inhibiting the immune response. Another collagenase-containing GI worsens ulcers by slowing down the healing process. GI consisted of fliD operon was found to be connected to flagellar assembly and biofilm production. By residing in biofilms, bacteria can avoid antibiotic therapy, resulting in chronic infection. Along with well-studied CagA and vacuolating toxin A (vacA) virulent genes, it is equally important to study these identified virulence factors for better understanding H. pylori-induced disease prognosis.     

Evolution in Quantum Leaps: Multiple Combinatorial Transfers of HPI and Other Genetic Modules in Enterobacteriaceae

Horizontal gene transfer is a key step in the evolution of Enterobacteriaceae. By acquiring virulence determinants of foreign origin, commensals can evolve into pathogens. In Enterobacteriaceae, horizontal transfer of these virulence determinants is largely dependent on transfer by plasmids, phages, genomic islands (GIs) and genomic modules (GMs). The High Pathogenicity Island (HPI) is a GI encoding virulence genes that can be transferred between different Enterobacteriaceae. We investigated the HPI because it was present in an Enterobacter hormaechei outbreak strain (EHOS). Genome sequence analysis showed that the EHOS contained an integration site for mobile elements and harbored two GIs and three putative GMs, including a new variant of the HPI (HPI-ICEEh1). We demonstrate, for the first time, that combinatorial transfers of GIs and GMs between Enterobacter cloacae complex isolates must have occurred. Furthermore, the excision and circularization of several combinations of the GIs and GMs was demonstrated. Because of its flexibility, the multiple integration site of mobile DNA can be considered an integration hotspot (IHS) that increases the genomic plasticity of the bacterium. Multiple combinatorial transfers of diverse combinations of the HPI and other genomic elements among Enterobacteriaceae may accelerate the generation of new pathogenic strains.     

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.     

Correlation between the distribution pattern of virulence genes and virulence of Aeromonas hydrophila strains

Nine strains of Aeromonas hydrophila isolated from diseased fish or soft-shelled tortoise were tested for the presence of three virulence genes including the genes encoding aerolysin,hemolysin,and extracellular serine protease (i.e.,aerA,hlyA,and ahpA,respectively).These genes were investigated using polymerase chain reaction (PCR)with specific primers for each gene.And the pathogenicities to Carrassius auratus ibebio of these strains were also assayed.PCR results demonstrated that the distribution patterns of aerA,hlyA,and ahpA were different in these strains.6/9 of A.hydrophila strains were aer A positive,8/9 of strains hly A positive,7/9 of strains ahp A positive,respectively.However,the assay for pathogenesis showed that two strains (A.hydrophila XS91-4-1 and C2)were strong virulent,two strains (A.hydrophila ST78-3-3 and 58-20-9)avirulent and the rest middle virulent was to the fish.In conclusion,there are significant correlation between the distribution pattern of the three virulence genes and the pathogenicity to Carrassius auratus ibebio.All strong virulent A.hydrophila strains were aerA+hlyA+ahpA+genotype,and all aerA+hlyA+ahpA+strains were virulent.Strains with the genotype of aerA-hlyA-ahpA+have middle pathogenicity.In the present study,we found for the first time that all A.hydrophila isolated from the ahpA positive were virulent to Carrassius auratus ibebio.Additionally,there was a positive correlation between the virulence of A.hydrophila and the presence of aerA and ahpA.     

Genomic subtraction to identify and characterize sequences of Shiga toxin-producing Escherichia coli O91:H21

To identify Shiga toxin-producing Escherichia coli genes associated with severe human disease, a genomic subtraction technique was used with hemolytic-uremic syndrome-associated O91:H21 strain CH014 and O6:H10 bovine strains. The method was adapted to the Shiga toxin-producing E. coli genome: three rounds of subtraction were used to isolate DNA fragments specific to strain CH014. The fragments were characterized by genetic support analysis, sequencing, and hybridization to the genome of a collection of Shiga toxin-producing E. coli strains. A total of 42 fragments were found, 19 of which correspond to previously identified unique DNA sequences in the enterohemorrhagic E. coli EDL933 reference strain, including 7 fragments corresponding to prophage sequences and others encoding candidate virulence factors, such a SepA homolog protein and a fimbrial usher protein. In addition, the subtraction procedure yielded plasmid-related sequences from Shigella flexneri and enteropathogenic and Shiga toxin-producing E. coli virulence plasmids. We found that lateral gene transfer is extensive in strain CH014, and we discuss the role of genomic mobile elements, especially bacteriophages, in the evolution and possible transfer of virulence determinants.     

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.     

In silico comparative study of the genomic islands of Vibrio cholerae MJ1236 with those of Classical and El Tor N16961 strains of Vibrio cholerae

The evolution of microbial genomes is greatly influenced by horizontal gene transfer (HGT), where large blocks of horizontally acquired foreign sequences, often encoding virulence determinants, occur in chromosomes of pathogenic bacteria. A program design-island developed in our laboratory was used on three completely sequenced Vibrio cholerae genomes, V. cholerae Classical O395, El Tor N16961 and MJ1236, in order to identify the putative horizontally acquired regions. The putative genomic islands (GIs) were graphically represented and analyzed. The study identified distinct regions in the GIs of V. cholerae MJ1236 which were shared either with the Classical or the El Tor strain of V. cholerae. A cluster comprising of 38 ORFs was common to V. cholerae strains of MJ1236 and Classical O395 but absent in El Tor N16961. About 5% of the predicted GIs of V. cholerae MJ1236 were unique to itself. Among these unique ORFs, a region of mostly hypothetical genes was identified, where the ORFs were present in a large cluster. The results show that the HGT had played a significant role in the evolution and the differentiation of V. cholerae MJ1236.     

Virulence potential and genetic diversity of Aeromonas caviae, Aeromonas veronii, and Aeromonas hydrophila clinical isolates from Mexico and Spain: a comparative study

A comparative study of 109 Aeromonas clinical isolates belonging to the 3 species most frequently isolated from patients with diarrhea in Mexico and Spain was performed to investigate the distribution of 3 prominent toxin genes and the gene encoding flagellin of lateral flagella; 4 well-established virulence factors in the genus Aeromonas. The aerolysin-hemolysin toxin genes were the most prevalent, being present in 89% of the total isolates. The ast toxin gene was conspicuously absent from the Aeromonas caviae and Aeromonas veronii groups but was present in 91% of the Aeromonas hydrophila isolates. Both the alt toxin gene and the lafA flagellin gene also had a low incidence in A. caviae and A. veronii. Differences in the prevalence of alt and lafA were observed between isolates from Mexico and Spain, confirming genus heterogeneity according to geographic location. Carriage of multiple toxin genes was primarily restricted to A. hydrophila isolates, suggesting that A. caviae and A. veronii isolates circulating in Mexico and Spain possess a limited array of virulence genes. Enterobacterial repetitive intergenetic consensus - polymerase chain reaction showed that the Aeromonas populations sampled lack dominant clones and were genetically heterogeneous, with A. caviae being the most diverse species. Further surveys of virulence determinants in genetically heterogeneous populations of Aeromonas isolates circulating worldwide are required to enhance the understanding of their capacity to cause disease.     

On the origin and function of an insertion element VPaI-1 specific to post-1995 pandemic Vibrio parahaemolyticus strains

Since 1995, new virulent strains of Vibrio parahaemolyticus have emerged and spread throughout the world. These "pandemic" strains have four strain specific genomic islands (GIs), which are considered to be potential factors of the pandemicity. We investigated the origin and function of 24 genes in the so-called VPaI-1, one of the four GIs, by searching homologs in various species in Bacteria and Archaea. Of these 24 genes, two are found only in Vibrio vulnificus CMCP6 and Shewanella sp. MR-7. The genomic segment (- 8 kb) encompassing the two genes shows the synteny among the three species. Since many of the Shewanella species can grow at 4 degrees C, these two genes may be candidates of adaptation to temperature stress. Further, we found a candidate for a swarming gene, which is reported as the V. cholerae virulence gene. Based on these findings, we hypothesized the emergence of pandemicity and discuss the mechanism for how these strains spread throughout the world.