Genomic analysis reveals that Pseudomonas aeruginosa virulence is combinatorial

Background

Pseudomonas aeruginosa is a ubiquitous environmental bacterium and an important opportunistic human pathogen. Generally, the acquisition of genes in the form of pathogenicity islands distinguishes pathogenic isolates from nonpathogens. We therefore sequenced a highly virulent strain of P. aeruginosa, PA14, and compared it with a previously sequenced (and less pathogenic) strain, PAO1, to identify novel virulence genes.

Results

The PA14 and PAO1 genomes are remarkably similar, although PA14 has a slightly larger genome (6.5 megabses [Mb]) than does PAO1 (6.3 Mb). We identified 58 PA14 gene clusters that are absent in PAO1 to determine which of these genes, if any, contribute to its enhanced virulence in a Caenorhabditis elegans pathogenicity model. First, we tested 18 additional diverse strains in the C. elegans model and observed a wide range of pathogenic potential; however, genotyping these strains using a custom microarray showed that the presence of PA14 genes that are absent in PAO1 did not correlate with the virulence of these strains. Second, we utilized a full-genome nonredundant mutant library of PA14 to identify five genes (absent in PAO1) required for C. elegans killing. Surprisingly, although these five genes are present in many other P. aeruginosa strains, they do not correlate with virulence in C. elegans.

Conclusion

Genes required for pathogenicity in one strain of P. aeruginosa are neither required for nor predictive of virulence in other strains. We therefore propose that virulence in this organism is both multifactorial and combinatorial, the result of a pool of pathogenicity-related genes that interact in various combinations in different genetic backgrounds.     

Genomic Features and Niche-Adaptation of Enterococcus faecium Strains from Korean Soybean-Fermented Foods

Certain strains of Enterococcus faecium contribute beneficially to human health and food fermentation. However, other E. faecium strains are opportunistic pathogens due to the acquisition of virulence factors and antibiotic resistance determinants. To characterize E. faecium from soybean fermentation, we sequenced the genomes of 10 E. faecium strains from Korean soybean-fermented foods and analyzed their genomes by comparing them with 51 clinical and 52 non-clinical strains of different origins. Hierarchical clustering based on 13,820 orthologous genes from all E. faecium genomes showed that the 10 strains are distinguished from most of the clinical strains. Like non-clinical strains, their genomes are significantly smaller than clinical strains due to fewer accessory genes associated with antibiotic resistance, virulence, and mobile genetic elements. Moreover, we identified niche-associated gene gain and loss from the soybean strains. Thus, we conclude that soybean E. faecium strains might have evolved to have distinctive genomic features that may contribute to its ability to thrive during soybean fermentation.     

Comparative Genomics and Characterization of Hybrid Shigatoxigenic and Enterotoxigenic Escherichia coli (STEC/ETEC) Strains

Background

Shigatoxigenic Escherichia coli (STEC) and enterotoxigenic E. coli (ETEC) cause serious foodborne infections in humans. These two pathogroups are defined based on the pathogroup-associated virulence genes: stx encoding Shiga toxin (Stx) for STEC and elt encoding heat-labile and/or est encoding heat-stable enterotoxin (ST) for ETEC. The study investigated the genomics of STEC/ETEC hybrid strains to determine their phylogenetic position among E. coli and to define the virulence genes they harbor.

Methods

The whole genomes of three STEC/ETEC strains possessing both stx and est genes were sequenced using PacBio RS sequencer. Two of the strains were isolated from the patients, one with hemolytic uremic syndrome, and one with diarrhea. The third strain was of bovine origin. Core genome analysis of the shared chromosomal genes and comparison with E. coli and Shigella spp. reference genomes was performed to determine the phylogenetic position of the STEC/ETEC strains. In addition, a set of virulence genes and ETEC colonization factors were extracted from the genomes. The production of Stx and ST were studied.

Results

The human STEC/ETEC strains clustered with strains representing ETEC, STEC, enteroaggregative E. coli, and commensal and laboratory-adapted E. coli. However, the bovine STEC/ETEC strain formed a remote cluster with two STECs of bovine origin. All three STEC/ETEC strains harbored several other virulence genes, apart from stx and est, and lacked ETEC colonization factors. Two STEC/ETEC strains produced both toxins and one strain Stx only.

Conclusions

This study shows that pathogroup-associated virulence genes of different E. coli can co-exist in strains originating from different phylogenetic lineages. The possibility of virulence genes to be associated with several E. coli pathogroups should be taken into account in strain typing and in epidemiological surveillance. Development of novel hybrid E. coli strains may cause a new public health risk, which challenges the traditional diagnostics of E. coli infections.     

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.     

Bioinformatics Analysis of the Complete Genome Sequence of the Mango Tree Pathogen Pseudomonas syringae pv. syringae UMAF0158 Reveals Traits Relevant to Virulence and Epiphytic Lifestyle

The genome sequence of more than 100 Pseudomonas syringae strains has been sequenced to date; however only few of them have been fully assembled, including P. syringae pv. syringae B728a. Different strains of pv. syringae cause different diseases and have different host specificities; so, UMAF0158 is a P. syringae pv. syringae strain related to B728a but instead of being a bean pathogen it causes apical necrosis of mango trees, and the two strains belong to different phylotypes of pv.syringae and clades of P. syringae. In this study we report the complete sequence and annotation of P. syringae pv. syringae UMAF0158 chromosome and plasmid pPSS158. A comparative analysis with the available sequenced genomes of other 25 P. syringae strains, both closed (the reference genomes DC3000, 1448A and B728a) and draft genomes was performed. The 5.8 Mb UMAF0158 chromosome has 59.3% GC content and comprises 5017 predicted protein-coding genes. Bioinformatics analysis revealed the presence of genes potentially implicated in the virulence and epiphytic fitness of this strain. We identified several genetic features, which are absent in B728a, that may explain the ability of UMAF0158 to colonize and infect mango trees: the mangotoxin biosynthetic operon mbo, a gene cluster for cellulose production, two different type III and two type VI secretion systems, and a particular T3SS effector repertoire. A mutant strain defective in the rhizobial-like T3SS Rhc showed no differences compared to wild-type during its interaction with host and non-host plants and worms. Here we report the first complete sequence of the chromosome of a pv. syringae strain pathogenic to a woody plant host. Our data also shed light on the genetic factors that possibly determine the pathogenic and epiphytic lifestyle of UMAF0158. This work provides the basis for further analysis on specific mechanisms that enable this strain to infect woody plants and for the functional analysis of host specificity in the P. syringae complex.     

Pathogenicity,Phylogenetic relationship and NGS based identification and assembly of tumorigenic Agrobacterium radiabacter plasmidic and chromosomic reads isolated from Prunus duclcis

Although many Agrobacterium radiobacter strains have already been identified, only a few genomes of strains belonging to genomovar G4 have been sequenced so far. In this study, we report the first virulent genome sequence of Agrobacterium radiobacter strain tun 183, which is highly virulent to almond specie. The genome size was estimated to be 5.53 Mb, with 57.9%GC content. In total, 6486 genes encoding proteins and 61 genes encoding RNAs were identified in this genome. Comparisons with the available sequenced genomes of genomovar G4 as well as with other A. sp. were conducted, revealing a hexapartite genome containing circular and linear chromosomes in addition to two accessory plasmids and a tumor inducing plasmid (pTi) in strain tun 183. The phylogenetic analysis of recA gene clearly showed the clustering of tun 183 strain within genomovar G4, supporting the monophyly within this genomovar.     

Comparative genomics of ocular Pseudomonas aeruginosa strains from keratitis patients with different clinical outcomes

The several virulence and drug-resistance mechanisms of Pseudomonas aeruginosa responsible for poor clinical outcomes in keratitis patients remain largely unknown. Here, we investigated the distribution of virulence factors and drug resistance by genes, mutations, efflux-pump systems of P. aeruginosa strains from keratitis patients with different clinical outcomes, included of whole-genome sequenced and annotated our five P. aeruginosa strains. Of the large number of virulence genes detected in all the genomes, MDR/XDR strains carry exoU and non-MDR strains carry exoS exotoxin of the type III secretion system, considered as main contributors of keratitis pathogenesis. However, several strain-specific virulence and resistance genes were detected in keratitis strains with poor outcome. Mainly, the flagellar genes fliC and fliD detected in the exoS carrying strains, reported to alter the host immune response, might impact the clinical outcome. This study may provide new insights into the genome of ocular strains and requires further functional studies.     

Complete Genome Sequence and Comparative Analysis of the Wild-type Commensal Escherichia coli Strain SE11 Isolated from a Healthy Adult

We sequenced and analyzed the genome of a commensal Escherichia coli (E. coli) strain SE11 (O152:H28) recently isolated from feces of a healthy adult and classified into E. coli phylogenetic group B1. SE11 harbored a 4.8 Mb chromosome encoding 4679 protein-coding genes and six plasmids encoding 323 protein-coding genes. None of the SE11 genes had sequence similarity to known genes encoding phage- and plasmid-borne virulence factors found in pathogenic E. coli strains. The comparative genome analysis with the laboratory strain K-12 MG1655 identified 62 poorly conserved genes between these two non-pathogenic strains and 1186 genes absent in MG1655. These genes in SE11 were mostly encoded in large insertion regions on the chromosome or in the plasmids, and were notably abundant in genes of fimbriae and autotransporters, which are cell surface appendages that largely contribute to the adherence ability of bacteria to host cells and bacterial conjugation. These data suggest that SE11 may have evolved to acquire and accumulate the functions advantageous for stable colonization of intestinal cells, and that the adhesion-associated functions are important for the commensality of E. coli in human gut habitat.Key words: Escherichia coli, commensal, human gut, genome sequencing     

Population Structure and Evolution of Non-O1/Non-O139 Vibrio cholerae by Multilocus Sequence Typing

Pathogenic non-O1/non-O139 Vibrio cholerae strains can cause sporadic outbreaks of cholera worldwide. In this study, multilocus sequence typing (MLST) of seven housekeeping genes was applied to 55 non-O1/non-O139 isolates from clinical and environmental sources. Data from five published O1 isolates and 17 genomes were also included, giving a total of 77 isolates available for analysis. There were 66 sequence types (STs), with the majority being unique, and only three clonal complexes. The V. cholerae strains can be divided into four subpopulations with evidence of recombination among the subpopulations. Subpopulations I and III contained predominantly clinical strains. PCR screening for virulence factors including Vibrio pathogenicity island (VPI), cholera toxin prophage (CTXΦ), type III secretion system (T3SS), and enterotoxin genes (rtxA and sto/stn) showed that combinations of these factors were present in the clinical isolates with 85.7% having rtxA, 51.4% T3SS, 31.4% VPI, 31.4% sto/stn (NAG-ST) and 11.4% CTXΦ. These factors were also present in environmental isolates but at a lower frequency. Five strains previously mis-identified as V. cholerae serogroups O114 to O117 were also analysed and formed a separate population with V. mimicus. The MLST scheme developed in this study provides a framework to identify sporadic cholera isolates by genetic identity.     

Genomic characterization of the Yersinia genus

Background

New DNA sequencing technologies have enabled detailed comparative genomic analyses of entire genera of bacterial pathogens. Prior to this study, three species of the enterobacterial genus Yersinia that cause invasive human diseases (Yersinia pestis, Yersinia pseudotuberculosis, and Yersinia enterocolitica) had been sequenced. However, there were no genomic data on the Yersinia species with more limited virulence potential, frequently found in soil and water environments.

Results

We used high-throughput sequencing-by-synthesis instruments to obtain 25- to 42-fold average redundancy, whole-genome shotgun data from the type strains of eight species: Y. aldovae, Y. bercovieri, Y. frederiksenii, Y. kristensenii, Y. intermedia, Y. mollaretii, Y. rohdei, and Y. ruckeri. The deepest branching species in the genus, Y. ruckeri, causative agent of red mouth disease in fish, has the smallest genome (3.7 Mb), although it shares the same core set of approximately 2,500 genes as the other members of the species, whose genomes range in size from 4.3 to 4.8 Mb. Yersinia genomes had a similar global partition of protein functions, as measured by the distribution of Cluster of Orthologous Groups families. Genome to genome variation in islands with genes encoding functions such as ureases, hydrogeneases and B-12 cofactor metabolite reactions may reflect adaptations to colonizing specific host habitats.

Conclusions

Rapid high-quality draft sequencing was used successfully to compare pathogenic and non-pathogenic members of the Yersinia genus. This work underscores the importance of the acquisition of horizontally transferred genes in the evolution of Y. pestis and points to virulence determinants that have been gained and lost on multiple occasions in the history of the genus.     

Chromosome‐level genome assembly for the horned‐gall aphid provides insights into interactions between gall‐making insect and its host plant

The aphid Schlechtendalia chinensis is an economically important insect that can induce horned galls, which are valuable for the medicinal and chemical industries. Up to now, more than twenty aphid genomes have been reported. Most of the sequenced genomes are derived from free‐living aphids. Here, we generated a high‐quality genome assembly from a galling aphid. The final genome assembly is 271.52 Mb, representing one of the smallest sequenced genomes of aphids. The genome assembly is based on contig and scaffold N50 values of the genome sequence are 3.77 Mb and 20.41 Mb, respectively. Nine‐seven percent of the assembled sequences was anchored onto 13 chromosomes. Based on BUSCO analysis, the assembly involved 96.9% of conserved arthropod and 98.5% of the conserved Hemiptera single‐copy orthologous genes. A total of 14,089 protein‐coding genes were predicted. Phylogenetic analysis revealed that S. chinensis diverged from the common ancestor of Eriosoma lanigerum approximately 57 million years ago (MYA). In addition, 35 genes encoding salivary gland proteins showed differentially when S. chinensis forms a gall, suggesting they have potential roles in gall formation and plant defense suppression. Taken together, this high‐quality S. chinensis genome assembly and annotation provide a solid genetic foundation for future research to reveal the mechanism of gall formation and to explore the interaction between aphids and their host plants.     

Microarray-based genomic surveying of gene polymorphisms in <Emphasis Type="Italic">Chlamydia trachomatis</Emphasis>

By comparing two fully sequenced genomes of Chlamydia trachomatis using competitive hybridization on DNA microarrays, a logarithmic correlation was demonstrated between the signal ratio of the arrays and the 75-99% range of nucleotide identities of the genes. Variable genes within 14 uncharacterized strains of C. trachomatis were identified by array analysis and verified by DNA sequencing. These genes may be crucial for understanding chlamydial virulence and pathogenesis.     

Plasmid Characterization and Chromosome Analysis of Two netF+ Clostridium perfringens Isolates Associated with Foal and Canine Necrotizing Enteritis

The recent discovery of a novel beta-pore-forming toxin, NetF, which is strongly associated with canine and foal necrotizing enteritis should improve our understanding of the role of type A Clostridium perfringens associated disease in these animals. The current study presents the complete genome sequence of two netF-positive strains, JFP55 and JFP838, which were recovered from cases of foal necrotizing enteritis and canine hemorrhagic gastroenteritis, respectively. Genome sequencing was done using Single Molecule, Real-Time (SMRT) technology-PacBio and Illumina Hiseq2000. The JFP55 and JFP838 genomes include a single 3.34 Mb and 3.53 Mb chromosome, respectively, and both genomes include five circular plasmids. Plasmid annotation revealed that three plasmids were shared by the two newly sequenced genomes, including a NetF/NetE toxins-encoding tcp-conjugative plasmid, a CPE/CPB2 toxins-encoding tcp-conjugative plasmid and a putative bacteriocin-encoding plasmid. The putative beta-pore-forming toxin genes, netF, netE and netG, were located in unique pathogenicity loci on tcp-conjugative plasmids. The C. perfringens JFP55 chromosome carries 2,825 protein-coding genes whereas the chromosome of JFP838 contains 3,014 protein-encoding genes. Comparison of these two chromosomes with three available reference C. perfringens chromosome sequences identified 48 (~247 kb) and 81 (~430 kb) regions unique to JFP55 and JFP838, respectively. Some of these divergent genomic regions in both chromosomes are phage- and plasmid-related segments. Sixteen of these unique chromosomal regions (~69 kb) were shared between the two isolates. Five of these shared regions formed a mosaic of plasmid-integrated segments, suggesting that these elements were acquired early in a clonal lineage of netF-positive C. perfringens strains. These results provide significant insight into the basis of canine and foal necrotizing enteritis and are the first to demonstrate that netF resides on a large and unique plasmid-encoded locus.     

江苏水产养殖区拟态弧菌种群结构及遗传多样性

【背景】拟态弧菌(Vibrio mimicus)是一种常见的革兰氏阴性病原菌,广泛分布于水环境和水生动物体内,可导致多种水产动物和人类感染。多位点序列分型(multilocus sequence typing, MLST)已被应用于多种病原菌的分子分型,其通过分析不同菌株之间的遗传关系,监测细菌传播的时间和地理分布,确定感染和传播途径,但目前未见有关拟态弧菌MLST的报道。【目的】开发一种基于MLST的拟态弧菌分型方法,并用于江苏水产养殖区拟态弧菌的种群结构和遗传进化分析,为拟态弧菌感染所引起的疾病防治提供理论基础。【方法】选择拟态弧菌的7个管家基因dnaE、gyrB、mdh、recA、rpoD、pntA和pyrH作为靶点,对江苏水产养殖区分离的155株拟态弧菌进行PCR扩增和测序。将测序结果分配等位基因,制作等位基因谱,分配不同的序列类型(sequence type, ST),利用软件goeBURST-1.2.1和MEGA-X对分配的ST型进行克隆复合体和遗传进化树聚类分析;此外,利用Kirby-Bauer圆盘扩散法测试155株拟态弧菌的药敏特性。【结果】155株拟态弧菌被分为56个STs,其中ST11占比最高;在双位点变异(double locus variants, DLV)水平分析发现56个STs分为3个克隆复合体和3个单体;系统发育树显示,56个STs被分为3个集群(cluster I、cluster II、cluster III)。药敏结果显示,155株拟态弧菌对红霉素类抗生素的耐药性最高(88.39%, 137/155),对氯霉素类抗生素敏感性最高(91.61%, 142/155)。【结论】本研究建立的MLST方法具有良好的分辨率,可作为拟态弧菌系统发育和未来流行病学调查有用的分子分型工具。根据抗生素耐药谱结果,提示在养殖过程中可选用氟苯尼考等国家批准使用的专用抗菌药对拟态弧菌进行防治。     

The genome sequence of Salmonella enterica serovar Choleraesuis, a highly invasive and resistant zoonotic pathogen

Salmonella enterica serovar Choleraesuis (S.Choleraesuis), a highly invasive serovar among non-typhoidal Salmonella, usually causes sepsis or extra-intestinal focal infections in humans. S.Choleraesuis infections have now become particularly difficult to treat because of the emergence of resistance to multiple antimicrobial agents. The 4.7 Mb genome sequence of a multidrug-resistant S.Choleraesuis strain SC-B67 was determined. Genome wide comparison of three sequenced Salmonella genomes revealed that more deletion events occurred in S.Choleraesuis SC-B67 and S.Typhi CT18 relative to S.Typhimurium LT2. S.Choleraesuis has 151 pseudogenes, which, among the three Salmonella genomes, include the highest percentage of pseudogenes arising from the genes involved in bacterial chemotaxis signal-transduction pathways. Mutations in these genes may increase smooth swimming of the bacteria, potentially allowing more effective interactions with and invasion of host cells to occur. A key regulatory gene of TetR/AcrR family, acrR, was inactivated through the introduction of an internal stop codon resulting in overexpression of AcrAB that appears to be associated with ciprofloxacin resistance. While lateral gene transfer providing basic functions to allow niche expansion in the host and environment is maintained during the evolution of different serovars of Salmonella, genes providing little overall selective benefit may be lost rapidly. Our findings suggest that the formation of pseudogenes may provide a simple evolutionary pathway that complements gene acquisition to enhance virulence and antimicrobial resistance in S.Choleraesuis.     

Comparative genomic analysis of two Burkholderia glumae strains from different geographic origins reveals a high degree of plasticity in genome structure associated with genomic islands

Burkholderia glumae is the major causal agent of bacterial panicle blight of rice, a growing disease problem in global rice production. To better understand its genome-scale characteristics, the genome of the highly virulent B. glumae strain 336gr-1 isolated from Louisiana, USA was sequenced using the Illumina Genome Analyser II system. De novo assembled 336gr-1 contigs were aligned and compared with the previously sequenced genome of B. glumae strain BGR1, which was isolated from an infected rice plant in South Korea. Comparative analysis of the whole genomes of B. glumae 336gr-1 and B. glumae BGR1 revealed numerous unique genomic regions present only in one of the two strains. These unique regions contained accessory genes including mobile elements and phage-related genes, and some of the unique regions in B. glumae BGR1 corresponded to predicted genomic islands. In contrast, little variation was observed in known and potential virulence genes between the two genomes. The considerable amount of plasticity largely based on accessory genes and genome islands observed from the comparison of the genomes of these two strains of B. glumae may explain the versatility of this bacterial species in various environmental conditions and geographic locations.     

The Conserved nhaAR Operon Is Drastically Divergent between B2 and Non-B2 Escherichia coli and Is Involved in Extra-Intestinal Virulence

The Escherichia coli species is divided in phylogenetic groups that differ in their virulence and commensal distribution. Strains belonging to the B2 group are involved in extra-intestinal pathologies but also appear to be more prevalent as commensals among human occidental populations. To investigate the genetic specificities of B2 sub-group, we used 128 sequenced genomes and identified genes of the core genome that showed marked difference between B2 and non-B2 genomes. We focused on the gene and its surrounding region with the strongest divergence between B2 and non-B2, the antiporter gene nhaA. This gene is part of the nhaAR operon, which is in the core genome but flanked by mobile regions, and is involved in growth at high pH and high sodium concentrations. Consistently, we found that a panel of non-B2 strains grew faster than B2 at high pH and high sodium concentrations. However, we could not identify differences in expression of the nhaAR operon using fluorescence reporter plasmids. Furthermore, the operon deletion had no differential impact between B2 and non-B2 strains, and did not result in a fitness modification in a murine model of gut colonization. Nevertheless, sequence analysis and experiments in a murine model of septicemia revealed that recombination in nhaA among B2 strains was observed in strains with low virulence. Finally, nhaA and nhaAR operon deletions drastically decreased virulence in one B2 strain. This effect of nhaAR deletion appeared to be stronger than deletion of all pathogenicity islands. Thus, a population genetic approach allowed us to identify an operon in the core genome without strong effect in commensalism but with an important role in extra-intestinal virulence, a landmark of the B2 strains.