Sandeep J. Joseph

Genomic characterization of the Bacillus cereus sensu lato species: Backdrop to the evolution of Bacillus anthracis

The key genes required for Bacillus anthracis to cause anthrax have been acquired recently by horizontal gene transfer. To understand the genetic background for the evolution of B. anthracis virulence, we obtained high-redundancy genome sequences of 45 strains of the Bacillus cereus sensu lato (s.l.) species that were chosen for their genetic diversity within the species based on the existing multilocus sequence typing scheme. From the resulting data, we called more than 324,000 new genes representing more than 12,333 new gene families for this group. The core genome size for the B. cereus s.l. group was ∼1750 genes, with another 2150 genes found in almost every genome constituting the extended core. There was a paucity of genes specific and conserved in any clade. We found no evidence of recent large-scale gene loss in B. anthracis or for unusual accumulation of nonsynonymous DNA substitutions in the chromosome; however, several B. cereus genomes isolated from soil and not previously associated with human disease were degraded to various degrees. Although B. anthracis has undergone an ecological shift within the species, its chromosome does not appear to be exceptional on a macroscopic scale compared with close relatives.

Hypervirulent Chlamydia trachomatis Clinical Strain Is a Recombinant between Lymphogranuloma Venereum (L2) and D Lineages

ABSTRACT Chlamydia trachomatis is an obligate intracellular bacterium that causes a diversity of severe and debilitating diseases worldwide. Sporadic and ongoing outbreaks of lymphogranuloma venereum (LGV) strains among men who have sex with men (MSM) support the need for research on virulence factors associated with these organisms. Previous analyses have been limited to single genes or genomes of laboratory-adapted reference strain L2/434 and outbreak strain L2b/UCH-1/proctitis. We characterized an unusual LGV strain, termed L2c, isolated from an MSM with severe hemorrhagic proctitis. L2c developed nonfusing, grape-like inclusions and a cytotoxic phenotype in culture, unlike the LGV strains described to date. Deep genome sequencing revealed that L2c was a recombinant of L2 and D strains with conserved clustered regions of genetic exchange, including a 78-kb region and a partial, yet functional, toxin gene that was lost with prolonged culture. Indels (insertions/deletions) were discovered in an ftsK gene promoter and in the tarp and hctB genes, which encode key proteins involved in replication, inclusion formation, and histone H1-like protein activity, respectively. Analyses suggest that these indels affect gene and/or protein function, supporting the in vitro and disease phenotypes. While recombination has been known to occur for C. trachomatis based on gene sequence analyses, we provide the first whole-genome evidence for recombination between a virulent, invasive LGV strain and a noninvasive common urogenital strain. Given the lack of a genetic system for producing stable C. trachomatis mutants, identifying naturally occurring recombinants can clarify gene function and provide opportunities for discovering avenues for genomic manipulation. IMPORTANCE Lymphogranuloma venereum (LGV) is a prevalent and debilitating sexually transmitted disease in developing countries, although there are significant ongoing outbreaks in Australia, Europe, and the United States among men who have sex with men (MSM). Relatively little is known about LGV virulence factors, and only two LGV genomes have been sequenced to date. We isolated an LGV strain from an MSM with severe hemorrhagic proctitis that was morphologically unique in tissue culture compared with other LGV strains. Bioinformatic and statistical analyses identified the strain as a recombinant of L2 and D strains with highly conserved clustered regions of genetic exchange. The unique culture morphology and, more importantly, disease phenotype could be traced to the genes involved in recombination. The findings have implications for bacterial species evolution and, in the case of ongoing LGV outbreaks, suggest that recombination is a mechanism for strain emergence that results in significant disease pathology. Lymphogranuloma venereum (LGV) is a prevalent and debilitating sexually transmitted disease in developing countries, although there are significant ongoing outbreaks in Australia, Europe, and the United States among men who have sex with men (MSM). Relatively little is known about LGV virulence factors, and only two LGV genomes have been sequenced to date. We isolated an LGV strain from an MSM with severe hemorrhagic proctitis that was morphologically unique in tissue culture compared with other LGV strains. Bioinformatic and statistical analyses identified the strain as a recombinant of L2 and D strains with highly conserved clustered regions of genetic exchange. The unique culture morphology and, more importantly, disease phenotype could be traced to the genes involved in recombination. The findings have implications for bacterial species evolution and, in the case of ongoing LGV outbreaks, suggest that recombination is a mechanism for strain emergence that results in significant disease pathology.

Population Genomics of Chlamydia trachomatis: Insights on Drift, Selection, Recombination, and Population Structure

The large number of sexually transmitted diseases and ocular trachoma cases that are caused globally each year by Chlamydia trachomatis has made this organism a World Health Organization priority for vaccine development. However, there is no gene transfer system for Chlamydia to help identify potential vaccine targets. To accelerate discoveries toward this goal, here we analyzed the broadest diversity of C. trachomatis genomes to date, including 25 geographically dispersed clinical and seven reference strains representing 14 of the 19 known serotypes. Strikingly, all 32 genomes were found to have evidence of DNA acquisition by homologous recombination in their history. Four distinct clades were identified, which correspond to all C. trachomatis disease phenotypes: lymphogranuloma venereum (LGV; Clade 1); noninvasive urogenital infections (Clade 2); ocular trachoma (Clade 3); and protocolitis (Clade 4; also includes some noninvasive urogenital infections). Although the ancestral relationship between clades varied, most strains acted as donor and recipient of recombination with no evidence for barriers to genetic exchange. The niche-specific LGV and trachoma clades have undergone less recombination, although the opportunity for mixing with strains from other clades that infect the rectal and ocular mucosa, respectively, is evident. Furthermore, there are numerous occasions for gene conversion events through sequential infections at the same anatomic sites. The size of recombinant segments is relatively small (~357 bp) compared with in vitro experiments of various C. trachomatis strains but is consistent with in vitro estimates for other bacterial species including Escherichia coli and Helicobacter pylori. Selection has also played a crucial role during the diversification of the organism. Clade 2 had the lowest nonsynonymous to synonymous ratio (dN/dS) but the highest effect of recombination, which is consistent with the widespread occurrence of synonymous substitutions in recombined genomic segments. The trachoma Clade 3 had the highest dN/dS estimates, which may be caused by an increased effect of genetic drift from niche specialization and a reduced effective population size. The degree of drift, selection, and recombination in C. trachomatis suggests that the challenge will remain to identify genomic regions that are stable and cross protective for the development of an efficacious vaccine.

Characterization of Aeromonas hydrophila Wound Pathotypes by Comparative Genomic and Functional Analyses of Virulence Genes

ABSTRACT Aeromonas hydrophila has increasingly been implicated as a virulent and antibiotic-resistant etiologic agent in various human diseases. In a previously published case report, we described a subject with a polymicrobial wound infection that included a persistent and aggressive strain of A. hydrophila (E1), as well as a more antibiotic-resistant strain of A. hydrophila (E2). To better understand the differences between pathogenic and environmental strains of A. hydrophila, we conducted comparative genomic and functional analyses of virulence-associated genes of these two wound isolates (E1 and E2), the environmental type strain A. hydrophila ATCC 7966T, and four other isolates belonging to A. aquariorum, A. veronii, A. salmonicida, and A. caviae. Full-genome sequencing of strains E1 and E2 revealed extensive differences between the two and strain ATCC 7966T. The more persistent wound infection strain, E1, harbored coding sequences for a cytotoxic enterotoxin (Act), a type 3 secretion system (T3SS), flagella, hemolysins, and a homolog of exotoxin A found in Pseudomonas aeruginosa. Corresponding phenotypic analyses with A. hydrophila ATCC 7966T and SSU as reference strains demonstrated the functionality of these virulence genes, with strain E1 displaying enhanced swimming and swarming motility, lateral flagella on electron microscopy, the presence of T3SS effector AexU, and enhanced lethality in a mouse model of Aeromonas infection. By combining sequence-based analysis and functional assays, we characterized an A. hydrophila pathotype, exemplified by strain E1, that exhibited increased virulence in a mouse model of infection, likely because of encapsulation, enhanced motility, toxin secretion, and cellular toxicity. IMPORTANCE Aeromonas hydrophila is a common aquatic bacterium that has increasingly been implicated in serious human infections. While many determinants of virulence have been identified in Aeromonas, rapid identification of pathogenic versus nonpathogenic strains remains a challenge for this genus, as it is for other opportunistic pathogens. This paper demonstrates, by using whole-genome sequencing of clinical Aeromonas strains, followed by corresponding virulence assays, that comparative genomics can be used to identify a virulent subtype of A. hydrophila that is aggressive during human infection and more lethal in a mouse model of infection. This aggressive pathotype contained genes for toxin production, toxin secretion, and bacterial motility that likely enabled its pathogenicity. Our results highlight the potential of whole-genome sequencing to transform microbial diagnostics; with further advances in rapid sequencing and annotation, genomic analysis will be able to provide timely information on the identities and virulence potential of clinically isolated microorganisms. Aeromonas hydrophila is a common aquatic bacterium that has increasingly been implicated in serious human infections. While many determinants of virulence have been identified in Aeromonas, rapid identification of pathogenic versus nonpathogenic strains remains a challenge for this genus, as it is for other opportunistic pathogens. This paper demonstrates, by using whole-genome sequencing of clinical Aeromonas strains, followed by corresponding virulence assays, that comparative genomics can be used to identify a virulent subtype of A. hydrophila that is aggressive during human infection and more lethal in a mouse model of infection. This aggressive pathotype contained genes for toxin production, toxin secretion, and bacterial motility that likely enabled its pathogenicity. Our results highlight the potential of whole-genome sequencing to transform microbial diagnostics; with further advances in rapid sequencing and annotation, genomic analysis will be able to provide timely information on the identities and virulence potential of clinically isolated microorganisms.

Interplay of recombination and selection in the genomes of Chlamydia trachomatis.

BackgroundChlamydia trachomatis is an obligate intracellular bacterial parasite, which causes several severe and debilitating diseases in humans. This study uses comparative genomic analyses of 12 complete published C. trachomatis genomes to assess the contribution of recombination and selection in this pathogen and to understand the major evolutionary forces acting on the genome of this bacterium.ResultsThe conserved core genes of C. trachomatis are a large proportion of the pan-genome: we identified 836 core genes in C. trachomatis out of a range of 874-927 total genes in each genome. The ratio of recombination events compared to mutation (ρ/θ) was 0.07 based on ancestral reconstructions using the ClonalFrame tool, but recombination had a significant effect on genetic diversification (r/m = 0.71). The distance-dependent decay of linkage disequilibrium also indicated that C. trachomatis populations behaved intermediately between sexual and clonal extremes. Fifty-five genes were identified as having a history of recombination and 92 were under positive selection based on statistical tests. Twenty-three genes showed evidence of being under both positive selection and recombination, which included genes with a known role in virulence and pathogencity (e.g., ompA, pmps, tarp). Analysis of inter-clade recombination flux indicated non-uniform currents of recombination between clades, which suggests the possibility of spatial population structure in C. trachomatis infections.ConclusionsC. trachomatis is the archetype of a bacterial species where recombination is relatively frequent yet gene gains by horizontal gene transfer (HGT) and losses (by deletion) are rare. Gene conversion occurs at sites across the whole C. trachomatis genome but may be more often fixed in genes that are under diversifying selection. Furthermore, genome sequencing will reveal patterns of serotype specific gene exchange and selection that will generate important research questions for understanding C. trachomatis pathogenesis.ReviewersThis article was reviewed by Dr. Jeremy Selengut, Dr. Lee S. Katz (nominated by Dr. I. King Jordan) and Dr. Arcady Mushegian.

Comparative Analysis of Chlamydia psittaci Genomes Reveals the Recent Emergence of a Pathogenic Lineage with a Broad Host Range

ABSTRACT Chlamydia psittaci is an obligate intracellular bacterium. Interest in Chlamydia stems from its high degree of virulence as an intestinal and pulmonary pathogen across a broad range of animals, including humans. C. psittaci human pulmonary infections, referred to as psittacosis, can be life-threatening, which is why the organism was developed as a bioweapon in the 20th century and is listed as a CDC biothreat agent. One remarkable recent result from comparative genomics is the finding of frequent homologous recombination across the genome of the sexually transmitted and trachoma pathogen Chlamydia trachomatis. We sought to determine if similar evolutionary dynamics occurred in C. psittaci. We analyzed 20 C. psittaci genomes from diverse strains representing the nine known serotypes of the organism as well as infections in a range of birds and mammals, including humans. Genome annotation revealed a core genome in all strains of 911 genes. Our analyses showed that C. psittaci has a history of frequently switching hosts and undergoing recombination more often than C. trachomatis. Evolutionary history reconstructions showed genome-wide homologous recombination and evidence of whole-plasmid exchange. Tracking the origins of recombinant segments revealed that some strains have imported DNA from as-yet-unsampled or -unsequenced C. psittaci lineages or other Chlamydiaceae species. Three ancestral populations of C. psittaci were predicted, explaining the current population structure. Molecular clock analysis found that certain strains are part of a clonal epidemic expansion likely introduced into North America by South American bird traders, suggesting that psittacosis is a recently emerged disease originating in New World parrots. IMPORTANCE Chlamydia psittaci is classified as a CDC biothreat agent based on its association with life-threatening lung disease, termed psittacosis, in humans. Because of the recent remarkable findings of frequent recombination across the genome of the human sexually transmitted and ocular trachoma pathogen Chlamydia trachomatis, we sought to determine if similar evolutionary dynamics occur in C. psittaci. Twenty C. psittaci genomes were analyzed from diverse strains that may play a pathogenic role in human disease. Evolution of the strains revealed genome-wide recombination occurring at a higher rate than for C. trachomatis. Certain strains were discovered to be part of a recent epidemic clonal expansion originating in South America. These strains may have been introduced into the United States from South American bird traders, suggesting that psittacosis is a recently emerged disease originating in New World parrots. Our analyses indicate that C. psittaci strains have a history of frequently switching hosts and undergoing recombination. Chlamydia psittaci is classified as a CDC biothreat agent based on its association with life-threatening lung disease, termed psittacosis, in humans. Because of the recent remarkable findings of frequent recombination across the genome of the human sexually transmitted and ocular trachoma pathogen Chlamydia trachomatis, we sought to determine if similar evolutionary dynamics occur in C. psittaci. Twenty C. psittaci genomes were analyzed from diverse strains that may play a pathogenic role in human disease. Evolution of the strains revealed genome-wide recombination occurring at a higher rate than for C. trachomatis. Certain strains were discovered to be part of a recent epidemic clonal expansion originating in South America. These strains may have been introduced into the United States from South American bird traders, suggesting that psittacosis is a recently emerged disease originating in New World parrots. Our analyses indicate that C. psittaci strains have a history of frequently switching hosts and undergoing recombination.

Bacterial population genomics and infectious disease diagnostics

New sequencing technologies have made the production of bacterial genome sequences increasingly easy, and it can be confidently forecasted that vast genomic databases will be generated in the next few years. Here, we detail how collections of bacterial genomes from a particular species (population genomics libraries) have already been used to improve the design of several diagnostic assays for bacterial pathogens. Genome sequencing itself is also becoming more commonly used for epidemiological, forensic and clinical investigations. There is an opportunity for the further development of bioinformatic tools to bring even further value to bacterial diagnostic genomics.

Population structure of Neisseria gonorrhoeae based on whole genome data and its relationship with antibiotic resistance

Neisseria gonorrhoeae is the causative agent of gonorrhea, a sexually transmitted infection (STI) of major importance. As a result of antibiotic resistance, there are now limited options for treating patients. We collected draft genome sequence data and associated metadata data on 76 N. gonorrhoeae strains from around the globe and searched for known determinants of antibiotics resistance within the strains. The population structure and evolutionary forces within the pathogen population were analyzed. Our results indicated a cosmopolitan gonoccocal population mainly made up of five subgroups. The estimated ratio of recombination to mutation (r/m = 2.2) from our data set indicates an appreciable level of recombination occurring in the population. Strains with resistance phenotypes to more recent antibiotics (azithromycin and cefixime) were mostly found in two of the five population subgroups.

Functional genomic characterization of virulence factors from necrotizing fasciitis-causing strains of Aeromonas hydrophila.

ABSTRACT The genomes of 10 Aeromonas isolates identified and designated Aeromonas hydrophila WI, Riv3, and NF1 to NF4; A. dhakensis SSU; A. jandaei Riv2; and A. caviae NM22 and NM33 were sequenced and annotated. Isolates NF1 to NF4 were from a patient with necrotizing fasciitis (NF). Two environmental isolates (Riv2 and -3) were from the river water from which the NF patient acquired the infection. While isolates NF2 to NF4 were clonal, NF1 was genetically distinct. Outside the conserved core genomes of these 10 isolates, several unique genomic features were identified. The most virulent strains possessed one of the following four virulence factors or a combination of them: cytotoxic enterotoxin, exotoxin A, and type 3 and 6 secretion system effectors AexU and Hcp. In a septicemic-mouse model, SSU, NF1, and Riv2 were the most virulent, while NF2 was moderately virulent. These data correlated with high motility and biofilm formation by the former three isolates. Conversely, in a mouse model of intramuscular infection, NF2 was much more virulent than NF1. Isolates NF2, SSU, and Riv2 disseminated in high numbers from the muscular tissue to the visceral organs of mice, while NF1 reached the liver and spleen in relatively lower numbers on the basis of colony counting and tracking of bioluminescent strains in real time by in vivo imaging. Histopathologically, degeneration of myofibers with significant infiltration of polymorphonuclear cells due to the highly virulent strains was noted. Functional genomic analysis provided data that allowed us to correlate the highly infectious nature of Aeromonas pathotypes belonging to several different species with virulence signatures and their potential ability to cause NF.

Genome-wide recombination in Chlamydia trachomatis

A new study reports comparative genomic analysis of 52 geographically diverse strains of Chlamydia trachomatis. The authors reconstruct a genome-wide phylogeny of the species and report extensive genome-wide recombination across multiple lineages of this intracellular bacterial pathogen.