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.     

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.     

Development of diagnostic and vaccine markers through cloning, expression, and regulation of putative virulence-protein-encoding genes of Aeromonas hydrophila

Aeromonas hydrophila is an opportunistic bacterial pathogen that is associated with a number of diseases in fish, amphibians, reptiles, and humans. In fish it causes several disease symptoms including tail and skin rot, and haemorrhagic septicemia; in human it causes soft-tissue wound infection and diarrhoea. The pathogenesis of A. hydrophila is multifactorial, but the mechanism is unknown so far. It is considered to be mediated by expression and secretion of extracellular proteins such as aerolysin, lipase, chitinase, amylase, gelatinase, hemolysins, and enterotoxins. A number of the putative virulence-protein-encoding genes that are present in the genome of A. hydrophila have been targeted by PCR for molecular diagnosis. These significant genes are also targeted for over-production of proteins by cloning and expression methods. In this review, we emphasize recent progress in the cloning, expression, and regulation of putative virulence-protein-encoding genes of A. hydrophila for a better understanding of the pathogenesis and also help to provide effective strategies for control of diseases.     

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.     

Genomic study of polyhydroxyalkanoates producing Aeromonas hydrophila 4AK4

The complete genome of Gram-negative Aeromonas hydrophila 4AK4 that has been used for industrial production of poly(3-hydroxybutyrate-co-3-hydroxyhexanoate) was sequenced and annotated. Its chromosome is 4,527,993 bp in size encoding 4,272 genes, including 28 rRNA genes and 104 tRNA genes. Comparative analysis indicated that genome of A. hydrophila 4AK4 was similar to that of the A. hydrophila ATCC 7966^T, an intensively studied aeromonad for its pathogenicity related to its genomic information. Genes possibly coming from other species or even other genus were identified in A. hydrophila 4AK4. A large number of putative virulent genes were predicted. However, a cytotonic enterotoxin (Ast) is absent in A. hydrophila 4AK4, allowing the industrial strain to be different from other A. hydrophila strains, indicating possible reduced virulence of strain 4AK4, which is very important for industrial fermentation. Genes involved in polyhydroxyalkanoate (PHA) metabolism were predicted and analyzed. The resulting genomic information is useful for improved production of PHA via metabolic engineering of A. hydrophila 4AK4.     

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.     

Implication of Lateral Genetic Transfer in the Emergence of Aeromonas hydrophila Isolates of Epidemic Outbreaks in Channel Catfish

To investigate the molecular basis of the emergence of Aeromonas hydrophila responsible for an epidemic outbreak of motile aeromonad septicemia of catfish in the Southeastern United States, we sequenced 11 A. hydrophila isolates that includes five reference and six recent epidemic isolates. Comparative genomics revealed that recent epidemic A. hydrophila isolates are highly clonal, whereas reference isolates are greatly diverse. We identified 55 epidemic-associated genetic regions with 313 predicted genes that are present in epidemic isolates but absent from reference isolates and 35% of these regions are located within genomic islands, suggesting their acquisition through lateral gene transfer. The epidemic-associated regions encode predicted prophage elements, pathogenicity islands, metabolic islands, fitness islands and genes of unknown functions, and 34 of the genes encoded in these regions were predicted as virulence factors. We found two pilus biogenesis gene clusters encoded within predicted pathogenicity islands. A functional metabolic island that encodes a complete pathway for myo-inositol catabolism was evident by the ability of epidemic A. hydrophila isolates to use myo-inositol as a sole carbon source. Testing of A. hydrophila field isolates found a consistent correlation between myo-inositol utilization as a sole carbon source and the presence of an epidemic-specific genetic marker. All epidemic isolates and one reference isolate shared a novel O-antigen cluster. Altogether we identified four different O-antigen biosynthesis gene clusters within the 11 sequenced A. hydrophila genomes. Our study reveals new insights into the evolutionary changes that have resulted in the emergence of recent epidemic A. hydrophila strains.     

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.     

Pathogenic and antimicrobial resistance genes in Streptococcus oralis strains revealed by comparative genome analysis

Streptococcus oralis is an early colonizer bacterium in dental plaques and is considered a potential pathogen of infective endocarditis (IE) disease. In this study, we built a complete genome map of Streptococcus oralis strain SOT, Streptococcus oralis strain SOD and Streptococcus infantis strain SO and performed comparative genomic analysis among these three strains. The results showed that there are five genomic islands (GIs) in strain SOT and one CRISPR in strain SOD. Each genome harbors various pathogenic genes related to diseases and drug resistance, while the antibiotic resistance genes in strains SOT and SOD were quite similar but different from those in strain SO. In addition, we identified 17 main virulence factors and capsule-related genes in three strains. These results suggest the pathogenic potential of Streptococcus strains, which lay a foundation for the prevention and treatment of a Streptococcus oralis infection.     

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.     

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.     

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.     

Comparative Genome Analysis Provides Insights into the Pathogenicity of Flavobacterium psychrophilum

Flavobacterium psychrophilum is a fish pathogen in salmonid aquaculture worldwide that causes cold water disease (CWD) and rainbow trout fry syndrome (RTFS). Comparative genome analyses of 11 F. psychrophilum isolates representing temporally and geographically distant populations were used to describe the F. psychrophilum pan-genome and to examine virulence factors, prophages, CRISPR arrays, and genomic islands present in the genomes. Analysis of the genomic DNA sequences were complemented with selected phenotypic characteristics of the strains. The pan genome analysis showed that F. psychrophilum could hold at least 3373 genes, while the core genome contained 1743 genes. On average, 67 new genes were detected for every new genome added to the analysis, indicating that F. psychrophilum possesses an open pan genome. The putative virulence factors were equally distributed among isolates, independent of geographic location, year of isolation and source of isolates. Only one prophage-related sequence was found which corresponded to the previously described prophage 6H, and appeared in 5 out of 11 isolates. CRISPR array analysis revealed two different loci with dissimilar spacer content, which only matched one sequence in the database, the temperate bacteriophage 6H. Genomic Islands (GIs) were identified in F. psychrophilum isolates 950106-1/1 and CSF 259–93, associated with toxins and antibiotic resistance. Finally, phenotypic characterization revealed a high degree of similarity among the strains with respect to biofilm formation and secretion of extracellular enzymes. Global scale dispersion of virulence factors in the genomes and the abilities for biofilm formation, hemolytic activity and secretion of extracellular enzymes among the strains suggested that F. psychrophilum isolates have a similar mode of action on adhesion, colonization and destruction of fish tissues across large spatial and temporal scales of occurrence. Overall, the genomic characterization and phenotypic properties may provide new insights to the mechanisms of pathogenicity in F. psychrophilum.     

Integrative analysis of transcriptome and genome indicates two potential genomic islands are associated with pathogenesis of Mycobacterium tuberculosis

Mycobacterium tuberculosis (M.tb) is a successful human pathogen and widely prevalent throughout the world. Genomic islands (GIs) are thought to be related to pathogenicity. In this study, we predicted two potential genomic islands in M.tb genome, respectively named as GI-1 and GI-2. It is indicated that the genes belong to PE_PGRS family in GI-1 and genes involved in sulfolipid-1 (SL-1) synthesis in GI-2 are strongly associated with M.tb pathogenesis. Sequence analysis revealed that the five PGRS genes are more polymorphic than other PGRS members in full virulence M.tb complex strains at significance level 0.01 but not in attenuated strains. Expression analysis of microarrays collected from literatures displayed that GI-1 genes, especially Rv3508 might be correlated with the response to the inhibition of aerobic respiration. Microarray analysis also showed that SL-1 cluster genes are drastically down-expressed in attenuated strains relative to full virulence strains. We speculated that the effect of SL-1 on M.tb pathogenicity could be associated with long-term survival and persistence establishment during infection. Additionally, the gene Rv3508 in GI-1 was under positive selection. Rv3508 may involve the response of M.tb to the inhibition of aerobic respiration by low oxygen or drug PA-824, and it may be a common feature of genes in GI-1. These findings may provide some novel insights into M.tb physiology and pathogenesis.     

Structure of a sigma28-regulated nonflagellar virulence protein from Campylobacter jejuni

Campylobacter jejuni, a Gram-negative motile bacterium, is a leading cause of human gastrointestinal infections. Although the mechanism of C.jejuni-mediated enteritis appears to be multifactorial, flagella play complex roles in the virulence of this human pathogen. Cj0977 is a recently identified virulence factor in C. jejuni and is expressed by a σ²⁸ promoter that controls late genes in the flagellar regulon. A Cj0977 mutant strain is fully motile but significantly reduced in the invasion of intestinal epithelial cells in vitro. Here, we report the crystal structure of the major structural domain of Cj0977, which reveals a homodimeric “hot-dog” fold architecture. Of note, the characteristic hot-dog fold has been found in various coenzyme A (CoA) compound binding proteins with numerous oligomeric states. Structural comparison with other known hot-dog fold proteins locates a putative binding site for an acyl-CoA compound in the Cj0977 protein. Structure-based site-directed mutagenesis followed by invasion assays indicates that key residues in the putative binding site are indeed essential for the Cj0977 virulence function, suggesting a possible function of Cj0977 as an acyl-CoA binding regulatory protein.     

Differentiation in Quinolone Resistance by Virulence Genotype in Pseudomonas aeruginosa

Pseudomonas aeruginosa is a leading pathogen that has become increasingly resistant to the fluoroquinolone antibiotics due to widespread prescribing. Adverse outcomes have been shown for patients infected with fluoroquinolone-resistant strains. The type III secretion system (TTSS) is a major virulence determinant during acute infections through the injection of effector toxins into host cells. Most strains exhibit a unique TTSS virulence genotype defined by the presence of either exoS or exoU gene encoding two of the effector toxins, ExoS and ExoU, respectively. Specific TTSS effector genotype has been shown previously to differentially impact virulence in pneumonia. In this study, we examined the relationship between TTSS effector genotype and fluoroquinolone resistance mechanisms in a collection of 270 respiratory isolates. We found that a higher proportion of exoU+ strains were fluoroquinolone-resistant compared to exoS+ strains (63% vs 49%, p = 0.03) despite its lower overall prevalence (38% exoU+ vs 56% exoS+). Results from sequencing the quinolone resistance determining regions (QRDRs) of the 4 target genes (gyrA, gyrB, parC, parE) indicated that strains containing the exoU gene were more likely to acquire ≥2 mutations than exoS+ strains at MICs ≤8 µg/ml (13% vs none) and twice as likely to have mutations in both gyrA and parC than exoS+ strains (48% vs 24% p = 0.0439). Our findings indicate that P. aeruginosa strains differentially develop resistance-conferring mutations that correlate with TTSS effector genotype and the more virulent exoU+ subpopulation. Differences in mutational processes by virulence genotype that were observed suggest co-evolution of resistance and virulence traits favoring a more virulent genotype in the quinolone-rich clinical environment.     

A comparative genomic analysis of Xanthomonas arboricola pv. juglandis strains reveal hallmarks of mobile genetic elements in the adaptation and accelerated evolution of virulence

Xanthomonas arboricola pv. juglandis (Xaj) is the most significant aboveground walnut bacterial pathogen. Disease management uses copper-based pesticides which induce pathogen resistance. We examined the genetic repertoire associated with adaptation and virulence evolution in Xaj. Comparative genomics of 32 Xaj strains reveal the possible acquisition and propagation of virulence factors via insertion sequences (IS). Fine-scale annotation revealed a Tn3 transposon (TnXaj417) encoding copper resistance genes acquired by horizontal gene transfer and associated with adaptation and tolerance to metal-based pesticides commonly used to manage pathogens in orchard ecosystems. Phylogenomic analysis reveals IS involvement in acquisition and diversification of type III effector proteins ranging from two to eight in non-pathogenic strains, 16 to 20 in pathogenic strains, besides six other putative effectors with a reduced identity degree found mostly among pathogenic strains. Yersiniabactin, xopK, xopAI, and antibiotic resistance genes are also located near ISs or inside genomic islands and structures resembling composite transposons.     

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.     

A novel strategy for the identification of genomic islands by comparative analysis of the contents and contexts of tRNA sites in closely related bacteria

We devised software tools to systematically investigate the contents and contexts of bacterial tRNA and tmRNA genes, which are known insertion hotspots for genomic islands (GIs). The strategy, based on MAUVE-facilitated multigenome comparisons, was used to examine 87 Escherichia coli MG1655 tRNA and tmRNA genes and their orthologues in E.coli EDL933, E.coli CFT073 and Shigella flexneri Sf301. Our approach identified 49 GIs occupying ~1.7 Mb that mapped to 18 tRNA genes, missing 2 but identifying a further 30 GIs as compared with Islander [Y. Mantri and K. P. Williams (2004), Nucleic Acids Res., 32, D55–D58]. All these GIs had many strain-specific CDS, anomalous GC contents and/or significant dinucleotide biases, consistent with foreign origins. Our analysis demonstrated marked conservation of sequences flanking both empty tRNA sites and tRNA-associated GIs across all four genomes. Remarkably, there were only 2 upstream and 5 downstream deletions adjacent to the 328 loci investigated. In silico PCR analysis based on conserved flanking regions was also used to interrogate hotspots in another eight completely or partially sequenced E.coli and Shigella genomes. The tools developed are ideal for the analysis of other bacterial species and will lead to in silico and experimental discovery of new genomic islands.