Phylogenomic analysis of a global collection of Escherichia coli ST38: evidence of interspecies and environmental transmission?

IMPORTANCE

Extraintestinal pathogenic Escherichia coli (ExPEC) sequence type (ST) 38 is one of the top 10 human pandemic lineages. Although a major cause of urinary tract and blood stream infections, ST38 has been poorly characterized from a global phylogenomic perspective. A comprehensive genome-scale analysis of 925 ST38 isolate genomes identified two broad ancestral clades and linkage of discrete ST38 clusters with specific bla CTX-M variants. In addition, the clades and clusters carry important virulence genes, with diverse but poorly characterized plasmids. Numerous putative interhost and environment transmission events were identified here by the presence of ST38 clones (defined as isolates with ≤35 SNPs) within humans, companion animals, food sources, urban birds, wildlife, and the environment. A small cluster of international ST38 clones from diverse sources, likely representing progenitors of a hospital outbreak that occurred in Brisbane, Australia, in 2017, was also identified. Our study emphasizes the importance of characterizing isolate genomes derived from nonhuman sources and geographical locations, without any selection bias.

Genome organization and genomics in Chlamydia: whole genome sequencing increases understanding of chlamydial virulence, evolution, and phylogeny

The genus Chlamydia contains important obligate intracellular bacterial pathogens to humans and animals, including C. trachomatis and C. pneumoniae. Since 1998, when the first Chlamydia genome was published, our understanding of how these microbes interact, evolved and adapted to different intracellular host environments has been transformed due to the expansion of chlamydial genomes. This review explores the current state of knowledge in Chlamydia genomics and how whole genome sequencing has revolutionised our understanding of Chlamydia virulence, evolution, and phylogeny over the past two and a half decades. This review will also highlight developments in multi-omics and other approaches that have complemented whole genome sequencing to advance knowledge of Chlamydia pathogenesis and future directions for chlamydial genomics.

Phylogenomics of two ST1 antibiotic-susceptible non-clinical Acinetobacter baumannii strains reveals multiple lineages and complex evolutionary history in global clone 1

Acinetobacter baumannii is an opportunistic pathogen that is difficult to treat due to its resistance to extreme conditions, including desiccation and antibiotics. Most strains causing outbreaks around the world belong to two main global lineages, namely global clones 1 and 2 (GC1 and GC2). Here, we used a combination of Illumina short read and MinION (Oxford Nanopore) long-read sequence data with a hybrid assembly approach to complete the genome sequence of two antibiotic-sensitive GC1 strains, Ex003 and Ax270, recovered in Lebanon from water and a rectal swab of a cat, respectively. Phylogenetic analysis of Ax270 and Ex003 with 186 publicly available GC1 genomes revealed two major clades, including five main lineages (L1–L5), and four single-isolate lineages outside of the two clades. Ax270 and Ex003, along with AB307-0294 and MRSN7213 (both predicted antibiotic-susceptible isolates) represent these individual lineages. Antibiotic resistance islands and transposons interrupting the comM gene remain important features in L1–L5, with L1 associated with the AbaR-type resistance islands, L2 with AbaR4, L3 strains containing either AbaR4 or its variants as well as Tn6022::ISAba42, and L4 and L5 associated with Tn6022 or its variants. Analysis of the capsule (KL) and outer core (OCL) polysaccharide loci further revealed a complex evolutionary history probably involving many recombination events. As more genomes become available, more GC1 lineages continue to emerge. However, genome sequence data from more diverse geographical regions are needed to draw a more accurate population structure of this globally distributed clone.

Two carbapenem-resistant ST1:ST231:KL1:OCL1 Acinetobacter baumannii strains recovered in Tehran, Iran, carry AbaR31 in the chromosome and AbaR4 and TnaphA6 in a RepAci6 plasmid

Objectives

To analyse the context of genes conferring antibiotic resistance in two carbapenem-resistant Acinetobacter baumannii isolates recovered in Tehran, Iran.

Methods

The antibiotic resistance phenotype for 28 antibiotics was determined using disc diffusion. The whole genome sequences of ABH008 and ABS200 were determined using the Illumina HiSeq X Ten platform. Resistance genes were identified using ResFinder and multilocus sequence types were determined using the Oxford and Institut Pasteur schemes.

Results

Isolates ABH008 and ABS200, recovered in 2012 and 2013, respectively, in two different Tehran hospitals, belong to the common global clone 1 lineage, ST1_IP and ST231_OX. They are resistant to sulfamethoxazole, tetracycline, gentamicin, amikacin, third-generation cephalosporins and carbapenems. Despite being isolated in different hospitals, phylogenetic analysis indicated they are closely related. Consistent with this, both isolates carry catA1, sul1, aacC1 and aadA1 in a novel variant of the AbaR3-type resistance island, named AbaR31. Both isolates are resistant to amikacin and carbapenems owing to aphA6 and oxa23, respectively. The oxa23 gene is located in the AbaR4 resistance island, and aphA6 in TnaphA6, and both mobile elements are in an ∼90 kbp plasmid encoding the putative RepAci6 replication initiation protein. Resistance to third-generation cephalosporins is due to the acquisition by homologous recombination of a 5 kb DNA segment that contains ISAba1-ampC from a ST623 strain.

Conclusions

The resistance gene complements of ABH008 and ABS200 were found in AbaR31 and a plasmid that encodes RepAci6. The close genetic relationship of ABH008 and ABS200, despite each being recovered from different hospitals, indicates transmission between the two hospitals.

Escherichia coli ST8196 is a novel, locally evolved, and extensively drug resistant pathogenic lineage within the ST131 clonal complex

The H30Rx subclade of Escherichia coli ST131 is a clinically important, globally dispersed pathogenic lineage that typically displays resistance to fluoroquinolones and extended spectrum β-lactams. Isolates EC233 and EC234, variants of ST131-H30Rx with a novel sequence type (ST) 8196, isolated from unrelated patients presenting with bacteraemia at a Sydney Hospital in 2014 are characterised here. EC233 and EC234 are phylogroup B2, serotype O25:H4A, and resistant to ampicillin, amoxicillin, cefoxitin, ceftazidime, ceftriaxone, ciprofloxacin, norfloxacin and gentamicin and are likely clonal. Both harbour an IncFII_2 plasmid (pSPRC_Ec234-FII) that carries most of the resistance genes on an IS26 associated translocatable unit, two small plasmids and a novel IncI1 plasmid (pSPRC_Ec234-I). SNP-based phylogenetic analysis of the core genome of representatives within the ST131 clonal complex places both isolates in a subclade with three clinical Australian ST131-H30Rx clade-C isolates. A MrBayes phylogeny analysis of EC233 and EC234 indicates ST8196 share a most recent common ancestor with ST131-H30Rx strain EC70 isolated from the same hospital in 2013. Our study identified genomic hallmarks that define the ST131-H30Rx subclade in the ST8196 isolates and highlights a need for unbiased genomic surveillance approaches to identify novel high-risk MDR E. coli pathogens that impact healthcare facilities.

Analysis of Complete Genome Sequence of Acinetobacter baumannii Strain ATCC 19606 Reveals Novel Mobile Genetic Elements and Novel Prophage

Acinetobacter baumannii isolate ATCC 19606 was recovered in the US prior to 1948. It has been used as a reference and model organism in many studies involving antibiotic resistance and pathogenesis of A. baumannii, while, until recently, a complete genome of this strain was not available. Here, we present an analysis of the complete 3.91-Mbp genome sequence, generated via a combination of short-read sequencing (Illumina) and long-read sequencing (MinION), and show it contains two small cryptic plasmids and a novel complete prophage of size 41.2 kb. We also characterised several regions of the ATCC 19606 genome, leading to the identification of a novel cadmium/mercury transposon, which was named Tn6551. ATCC 19606 is an antibiotic-sensitive strain, but a comparative analysis of all publicly available ST52 strains predicts a resistance to modern antibiotics by the accumulation of antibiotic-resistance genes via plasmids in recent isolates that belong to this sequence type.

Chromatin accessibility dynamics of Chlamydia-infected epithelial cells

Chlamydia are Gram-negative, obligate intracellular bacterial pathogens responsible for a broad spectrum of human and animal diseases. In humans, Chlamydia trachomatis is the most prevalent bacterial sexually transmitted infection worldwide and is the causative agent of trachoma (infectious blindness) in disadvantaged populations. Over the course of its developmental cycle, Chlamydia extensively remodels its intracellular niche and parasitises the host cell for nutrients, with substantial resulting changes to the host cell transcriptome and proteome. However, little information is available on the impact of chlamydial infection on the host cell epigenome and global gene regulation. Regions of open eukaryotic chromatin correspond to nucleosome-depleted regions, which in turn are associated with regulatory functions and transcription factor binding. We applied formaldehyde-assisted isolation of regulatory elements enrichment followed by sequencing (FAIRE-Seq) to generate temporal chromatin maps of C. trachomatis-infected human epithelial cells in vitro over the chlamydial developmental cycle. We detected both conserved and distinct temporal changes to genome-wide chromatin accessibility associated with C. trachomatis infection. The observed differentially accessible chromatin regions include temporally-enriched sets of transcription factors, which may help shape the host cell response to infection. These regions and motifs were linked to genomic features and genes associated with immune responses, re-direction of host cell nutrients, intracellular signalling, cell-cell adhesion, extracellular matrix, metabolism and apoptosis. This work provides another perspective to the complex response to chlamydial infection, and will inform further studies of transcriptional regulation and the epigenome in Chlamydia-infected human cells and tissues.

Genomic profiling of Escherichia coli isolates from bacteraemia patients: a 3-year cohort study of isolates collected at a Sydney teaching hospital

This study sought to assess the genetic variability of Escherichia coli isolated from bloodstream infections (BSIs) presenting at Concord Hospital, Sydney during 2013-2016. Whole-genome sequencing was used to characterize 81 E. coli isolates sourced from community-onset (CO) and hospital-onset (HO) BSIs. The cohort comprised 64 CO and 17 HO isolates, including 35 multidrug-resistant (MDR) isolates exhibiting phenotypic resistance to three or more antibiotic classes. Phylogenetic analysis identified two major ancestral clades. One was genetically diverse with 25 isolates distributed in 16 different sequence types (STs) representing phylogroups A, B1, B2, C and F, while the other comprised phylogroup B2 isolates in subclades representing the ST131, ST73 and ST95 lineages. Forty-seven isolates contained a class 1 integron, of which 14 carried blaCTX -M-gene. Isolates with a class 1 integron carried more antibiotic resistance genes than isolates without an integron and, in most instances, resistance genes were localized within complex resistance loci (CRL). Resistance to fluoroquinolones could be attributed to point mutations in chromosomal parC and gyrB genes and, in addition, two isolates carried a plasmid-associated qnrB4 gene. Co-resistance to fluoroquinolone and broad-spectrum beta-lactam antibiotics was associated with ST131 (HO and CO), ST38 (HO), ST393 (CO), ST2003 (CO) and ST8196 (CO and HO), a novel ST identified in this study. Notably, 10/81 (12.3 %) isolates with ST95 (5 isolates), ST131 (2 isolates), ST88 (2 isolates) and a ST540 likely carry IncFII-IncFIB plasmid replicons with a full spectrum of virulence genes consistent with the carriage of ColV-like plasmids. Our data indicate that IncF plasmids play an important role in shaping virulence and resistance gene carriage in BSI E. coli in Australia.

Early Transcriptional Landscapes of Chlamydia trachomatis-Infected Epithelial Cells at Single Cell Resolution

Chlamydia are Gram-negative obligate intracellular bacterial pathogens responsible for a variety of disease in humans and animals worldwide. Chlamydia trachomatis causes trachoma in disadvantaged populations, and is the most common bacterial sexually transmitted infection in humans, causing reproductive tract disease. Antibiotic therapy successfully treats diagnosed chlamydial infections, however asymptomatic infections are common. High-throughput transcriptomic approaches have explored chlamydial gene expression and infected host cell gene expression. However, these were performed on large cell populations, averaging gene expression profiles across all cells sampled and potentially obscuring biologically relevant subsets of cells. We generated a pilot dataset, applying single cell RNA-Seq (scRNA-Seq) to C. trachomatis infected and mock-infected epithelial cells to assess the utility, pitfalls and challenges of single cell approaches applied to chlamydial biology, and to potentially identify early host cell biomarkers of chlamydial infection. Two hundred sixty-four time-matched C. trachomatis-infected and mock-infected HEp-2 cells were collected and subjected to scRNA-Seq. After quality control, 200 cells were retained for analysis. Two distinct clusters distinguished 3-h cells from 6- and 12-h. Pseudotime analysis identified a possible infection-specific cellular trajectory for Chlamydia-infected cells, while differential expression analyses found temporal expression of metallothioneins and genes involved with cell cycle regulation, innate immune responses, cytoskeletal components, lipid biosynthesis and cellular stress. We find that changes to the host cell transcriptome at early times of C. trachomatis infection are readily discernible by scRNA-Seq, supporting the utility of single cell approaches to identify host cell biomarkers of chlamydial infection, and to further deconvolute the complex host response to infection.

Bioinformatic analysis of bacteria and host cell dual RNA-sequencing experiments

Bacterial pathogens subvert host cells by manipulating cellular pathways for survival and replication; in turn, host cells respond to the invading pathogen through cascading changes in gene expression. Deciphering these complex temporal and spatial dynamics to identify novel bacterial virulence factors or host response pathways is crucial for improved diagnostics and therapeutics. Dual RNA sequencing (dRNA-Seq) has recently been developed to simultaneously capture host and bacterial transcriptomes from an infected cell. This approach builds on the high sensitivity and resolution of RNA sequencing technology and is applicable to any bacteria that interact with eukaryotic cells, encompassing parasitic, commensal or mutualistic lifestyles. Several laboratory protocols have been presented that outline the collection, extraction and sequencing of total RNA for dRNA-Seq experiments, but there is relatively little guidance available for the detailed bioinformatic analyses required. This protocol outlines a typical dRNA-Seq experiment, based on a Chlamydia trachomatis-infected host cell, with a detailed description of the necessary bioinformatic analyses with currently available software tools.

Dual RNA-Seq of Chlamydia and Host Cells

During the infection of a host cell by a bacterial pathogen, a cascading series of gene expression changes occurs as each organism manipulates or responds to the other via defense or survival strategies. Unraveling this complex interplay is key for our understanding of bacterial virulence and host response pathways for the development of novel therapeutics. Dual RNA sequencing (dual RNA-Seq) has recently been developed to simultaneously capture host and bacterial transcriptomes from an infected cell. Leveraging the sensitivity and resolution allowed by RNA-seq, dual RNA-Seq can be applied to any bacteria-eukaryotic host interaction. We pioneered dual RNA-Seq to simultaneously capture Chlamydia and host expression profiles during an in vitro infection as proof of principle. Here we provide a detailed laboratory protocol and bioinformatics analysis guidelines for dual RNA-seq experiments focusing on Chlamydia as the organism of interest.

Development of a Chlamydia suis-specific antibody enzyme-linked immunosorbent assay based on the use of a B-cell epitope of the polymorphic membrane protein C

Chlamydia suis infections lead to economic loss in the pork industry. Chlamydia suis infections could be successfully treated with tetracyclines until the appearance of a tetracycline resistant phenotype, which was acquired via horizontal gene transfer of the tet(C) gene. Given the importance of C. suis as a swine pathogen and as a recently emerged tetracycline resistant pathogen with zoonotic potential, our aim was to develop a sensitive C. suis-specific antibody ELISA based on the polymorphic membrane proteins (Pmps). Chlamydia Pmps are important virulence factors and candidate antigens for serodiagnosis. We identified nine Pmps (PmpA to I) in C. suis strain MD56 using a recently developed Hidden-Markov model. PmpC was the most promising candidate for the development of a C. suis-specific antibody ELISA as the protein was absent in C. abortus, C. pecorum and C. psittaci which also infect pigs and as the protein contained C. suis-specific amino acid regions, absent in C. trachomatis PmpC. We identified an immunodominant B-cell epitope in C. suis PmpC using experimental porcine sera. The sensitivity and specificity of the PmpC ELISA was compared to the complement fixation test (CFT) and to a recombinant MOMP ELISA using experimental sera. The PmpC ELISA detected all positive control sera and was in contrast to CFT and the rMOMP ELISA 100% C. suis specific as positive control sera against other Chlamydia species did not react in the PmpC ELISA. The test was successfully validated using slaughterhouse sera and sera from clinically affected pigs. The PmpC ELISA could assist in diminishing the spread of C. suis infections in the pork industry.

A Laboratory Methodology for Dual RNA-Sequencing of Bacteria and their Host Cells In Vitro

Dual RNA-Sequencing leverages established next-generation sequencing (NGS)-enabled RNA-Seq approaches to measure genome-wide transcriptional changes of both an infecting bacteria and host cells. By simultaneously investigating both organisms from the same biological sample, dual RNA-Seq can provide unique insight into bacterial infection processes and reciprocal host responses at once. However, the difficulties involved in handling both prokaryotic and eukaryotic material require distinct, optimized procedures. We previously developed and applied dual RNA-Seq to measure prokaryotic and eukaryotic expression profiles of human cells infected with bacteria, using in vitro Chlamydia-infected epithelial cells as proof of principle. Here we provide a detailed laboratory protocol for in vitro dual RNA-Seq that is readily adaptable to any host-bacteria system of interest.

Genetic diversity in the plasticity zone and the presence of the chlamydial plasmid differentiates Chlamydia pecorum strains from pigs, sheep, cattle, and koalas

Background

Chlamydia pecorum is a globally recognised pathogen of livestock and koalas. To date, comparative genomics of C. pecorum strains from sheep, cattle and koalas has revealed that only single nucleotide polymorphisms (SNPs) and a limited number of pseudogenes appear to contribute to the genetic diversity of this pathogen. No chlamydial plasmid has been detected in these strains despite its ubiquitous presence in almost all other chlamydial species. Genomic analyses have not previously included C. pecorum from porcine hosts. We sequenced the genome of three C. pecorum isolates from pigs with differing pathologies in order to re-evaluate the genetic differences and to update the phylogenetic relationships between C. pecorum from each of the hosts.

Methods

Whole genome sequences for the three porcine C. pecorum isolates (L1, L17 and L71) were acquired using C. pecorum-specific sequence capture probes with culture-independent methods, and assembled in CLC Genomics Workbench. The pairwise comparative genomic analyses of 16 pig, sheep, cattle and koala C. pecorum genomes were performed using several bioinformatics platforms, while the phylogenetic analyses of the core C. pecorum genomes were performed with predicted recombination regions removed. Following the detection of a C. pecorum plasmid, a newly developed C. pecorum-specific plasmid PCR screening assay was used to evaluate the plasmid distribution in 227 C. pecorum samples from pig, sheep, cattle and koala hosts.

Results

Three porcine C. pecorum genomes were sequenced using C. pecorum-specific sequence capture probes with culture-independent methods. Comparative genomics of the newly sequenced porcine C. pecorum genomes revealed an increased average number of SNP differences (~11 500) between porcine and sheep, cattle, and koala C. pecorum strains, compared to previous C. pecorum genome analyses. We also identified a third copy of the chlamydial cytotoxin gene, found only in porcine C. pecorum isolates. Phylogenetic analyses clustered porcine isolates into a distinct clade, highlighting the polyphyletic origin of C. pecorum in livestock.

Most surprising, we also discovered a plasmid in the porcine C. pecorum genome. Using this novel C. pecorum plasmid (pCpec) sequence, a) we developed a pCpec screening assay to evaluate the plasmid distribution in C. pecorum from different hosts; and b) to characterise the pCpec sequences from available previously sequenced C. pecorum genome data. pCpec screening showed that the pCpec is common in all hosts of C. pecorum, however not all C. pecorum strains carry pCpec.

Conclusions

This study provides further insight into the complexity of C. pecorum epidemiology and novel genomic regions that may be linked to host specificity. C. pecorum plasmid characterisation may aid in improving our understanding of C. pecorum pathogenesis across the variety of host species this animal pathogen infects.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2053-8) contains supplementary material, which is available to authorized users.

Chlamydia caviae infection alters abundance but not composition of the guinea pig vaginal microbiota

In humans, the vaginal microbiota is thought to be the first line of defense again pathogens including Chlamydia trachomatis. The guinea pig has been extensively used as a model to study chlamydial infection because it shares anatomical and physiological similarities with humans, such as a squamous vaginal epithelium as well as some of the long-term outcomes caused by chlamydial infection. In this study, we aimed to evaluate the guinea pig-C. caviae model of genital infection as a surrogate for studying the role of the vaginal microbiota in the early steps of C. trachomatis infection in humans. We used culture-independent molecular methods to characterize the relative and absolute abundance of bacterial phylotypes in the guinea pig vaginal microbiota in animals non-infected, mock-infected or infected by C. caviae. We showed that the guinea pig and human vaginal microbiotas are of different bacterial composition and abundance. Chlamydia caviae infection had a profound effect on the absolute abundance of bacterial phylotypes but not on the composition of the guinea pig vaginal microbiota. Our findings compromise the validity of the guinea pig-C. caviae model to study the role of the vaginal microbiota during the early steps of sexually transmitted infection.

The vaginal microbiota of the guinea pig differs from that of humans and cannot prevent chlamydial infections efficiently.

Comparative genomics of koala, cattle and sheep strains of Chlamydia pecorum

BACKGROUND

Chlamydia pecorum is an important pathogen of domesticated livestock including sheep, cattle and pigs. This pathogen is also a key factor in the decline of the koala in Australia. We sequenced the genomes of three koala C. pecorum strains, isolated from the urogenital tracts and conjunctiva of diseased koalas. The genome of the C. pecorum VR629 (IPA) strain, isolated from a sheep with polyarthritis, was also sequenced.

RESULTS

Comparisons of the draft C. pecorum genomes against the complete genomes of livestock C. pecorum isolates revealed that these strains have a conserved gene content and order, sharing a nucleotide sequence similarity > 98%. Single nucleotide polymorphisms (SNPs) appear to be key factors in understanding the adaptive process. Two regions of the chromosome were found to be accumulating a large number of SNPs within the koala strains. These regions include the Chlamydia plasticity zone, which contains two cytotoxin genes (toxA and toxB), and a 77 kbp region that codes for putative type III effector proteins. In one koala strain (MC/MarsBar), the toxB gene was truncated by a premature stop codon but is full-length in IPTaLE and DBDeUG. Another five pseudogenes were also identified, two unique to the urogenital strains C. pecorum MC/MarsBar and C. pecorum DBDeUG, respectively, while three were unique to the koala C. pecorum conjunctival isolate IPTaLE. An examination of the distribution of these pseudogenes in C. pecorum strains from a variety of koala populations, alongside a number of sheep and cattle C. pecorum positive samples from Australian livestock, confirmed the presence of four predicted pseudogenes in koala C. pecorum clinical samples. Consistent with our genomics analyses, none of these pseudogenes were observed in the livestock C. pecorum samples examined. Interestingly, three SNPs resulting in pseudogenes identified in the IPTaLE isolate were not found in any other C. pecorum strain analysed, raising questions over the origin of these point mutations.

CONCLUSIONS

The genomic data revealed that variation between C. pecorum strains were mainly due to the accumulation of SNPs, some of which cause gene inactivation. The identification of these genetic differences will provide the basis for further studies to understand the biology and evolution of this important animal pathogen.

Genome sequencing and comparative analysis of three Chlamydia pecorum strains associated with different pathogenic outcomes

Background

Chlamydia pecorum is the causative agent of a number of acute diseases, but most often causes persistent, subclinical infection in ruminants, swine and birds. In this study, the genome sequences of three C. pecorum strains isolated from the faeces of a sheep with inapparent enteric infection (strain W73), from the synovial fluid of a sheep with polyarthritis (strain P787) and from a cervical swab taken from a cow with metritis (strain PV3056/3) were determined using Illumina/Solexa and Roche 454 genome sequencing.

Results

Gene order and synteny was almost identical between C. pecorum strains and C. psittaci. Differences between C. pecorum and other chlamydiae occurred at a number of loci, including the plasticity zone, which contained a MAC/perforin domain protein, two copies of a >3400 amino acid putative cytotoxin gene and four (PV3056/3) or five (P787 and W73) genes encoding phospholipase D. Chlamydia pecorum contains an almost intact tryptophan biosynthesis operon encoding trpABCDFR and has the ability to sequester kynurenine from its host, however it lacks the genes folA, folKP and folB required for folate metabolism found in other chlamydiae. A total of 15 polymorphic membrane proteins were identified, belonging to six pmp families. Strains possess an intact type III secretion system composed of 18 structural genes and accessory proteins, however a number of putative inc effector proteins widely distributed in chlamydiae are absent from C. pecorum. Two genes encoding the hypothetical protein ORF663 and IncA contain variable numbers of repeat sequences that could be associated with persistence of infection.

Conclusions

Genome sequencing of three C. pecorum strains, originating from animals with different disease manifestations, has identified differences in ORF663 and pseudogene content between strains and has identified genes and metabolic traits that may influence intracellular survival, pathogenicity and evasion of the host immune system.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-23) contains supplementary material, which is available to authorized users.

Simultaneous transcriptional profiling of bacteria and their host cells

We developed an RNA-Seq-based method to simultaneously capture prokaryotic and eukaryotic expression profiles of cells infected with intracellular bacteria. As proof of principle, this method was applied to Chlamydia trachomatis-infected epithelial cell monolayers in vitro, successfully obtaining transcriptomes of both C. trachomatis and the host cells at 1 and 24 hours post-infection. Chlamydiae are obligate intracellular bacterial pathogens that cause a range of mammalian diseases. In humans chlamydiae are responsible for the most common sexually transmitted bacterial infections and trachoma (infectious blindness). Disease arises by adverse host inflammatory reactions that induce tissue damage & scarring. However, little is known about the mechanisms underlying these outcomes. Chlamydia are genetically intractable as replication outside of the host cell is not yet possible and there are no practical tools for routine genetic manipulation, making genome-scale approaches critical. The early timeframe of infection is poorly understood and the host transcriptional response to chlamydial infection is not well defined. Our simultaneous RNA-Seq method was applied to a simplified in vitro model of chlamydial infection. We discovered a possible chlamydial strategy for early iron acquisition, putative immune dampening effects of chlamydial infection on the host cell, and present a hypothesis for Chlamydia-induced fibrotic scarring through runaway positive feedback loops. In general, simultaneous RNA-Seq helps to reveal the complex interplay between invading bacterial pathogens and their host mammalian cells and is immediately applicable to any bacteria/host cell interaction.

Genus-optimized strategy for the identification of chlamydial type III secretion substrates

Among chlamydial virulence factors are the type III secretion (T3S) system and its effectors. T3S effectors target host proteins to benefit the infecting chlamydiae. The assortment of effectors, each with a unique function, varies between species. This variation likely contributes to differences in host specificity and disease severity. A dozen effectors of Chlamydia trachomatis have been identified; however, estimates suggest that more exist. A T3S prediction algorithm, SVM-based Identification and Evaluation of Virulence Effectors (SIEVE), along with a Yersinia surrogate secretion system helped to identify a new T3S substrate, CT082, which rather than functioning as an effector associates with the chlamydial envelope after secretion. SIEVE was modified to improve/expand effector predictions to include all sequenced genomes. Additional adjustments were made to the existing surrogate system whereby the N terminus of putative effectors was fused to a known effector lacking its own N terminus and was tested for secretion. Expansion of effector predictions by cSIEVE and modification of the surrogate system have also assisted in identifying a new T3S substrate from C. psittaci. The expanded predictions along with modifications to improve the surrogate secretion system have enhanced our ability to identify novel species-specific effectors, which upon characterization should provide insight into the unique pathogenic properties of each species.

Unity in variety--the pan-genome of the Chlamydiae

Chlamydiae are evolutionarily well-separated bacteria that live exclusively within eukaryotic host cells. They include important human pathogens such as Chlamydia trachomatis as well as symbionts of protozoa. As these bacteria are experimentally challenging and genetically intractable, our knowledge about them is still limited. In this study, we obtained the genome sequences of Simkania negevensis Z, Waddlia chondrophila 2032/99, and Parachlamydia acanthamoebae UV-7. This enabled us to perform the first comprehensive comparative and phylogenomic analysis of representative members of four major families of the Chlamydiae, including the Chlamydiaceae. We identified a surprisingly large core gene set present in all genomes and a high number of diverse accessory genes in those Chlamydiae that do not primarily infect humans or animals, including a chemosensory system in P. acanthamoebae and a type IV secretion system. In S. negevensis, the type IV secretion system is encoded on a large conjugative plasmid (pSn, 132 kb). Phylogenetic analyses suggested that a plasmid similar to the S. negevensis plasmid was originally acquired by the last common ancestor of all four families and that it was subsequently reduced, integrated into the chromosome, or lost during diversification, ultimately giving rise to the extant virulence-associated plasmid of pathogenic chlamydiae. Other virulence factors, including a type III secretion system, are conserved among the Chlamydiae to variable degrees and together with differences in the composition of the cell wall reflect adaptation to different host cells including convergent evolution among the four chlamydial families. Phylogenomic analysis focusing on chlamydial proteins with homology to plant proteins provided evidence for the acquisition of 53 chlamydial genes by a plant progenitor, lending further support for the hypothesis of an early interaction between a chlamydial ancestor and the primary photosynthetic eukaryote.

Comparison of koala LPCoLN and human strains of Chlamydia pneumoniae highlights extended genetic diversity in the species

Background

Chlamydia pneumoniae is a widespread pathogen causing upper and lower respiratory tract infections in addition to a range of other diseases in humans and animals. Previous whole genome analyses have focused on four essentially clonal (> 99% identity) C. pneumoniae human genomes (AR39, CWL029, J138 and TW183), providing relatively little insight into strain diversity and evolution of this species.

Results

We performed individual gene-by-gene comparisons of the recently sequenced C. pneumoniae koala genome and four C. pneumoniae human genomes to identify species-specific genes, and more importantly, to gain an insight into the genetic diversity and evolution of the species. We selected genes dispersed throughout the chromosome, representing genes that were specific to C. pneumoniae, genes with a demonstrated role in chlamydial biology and/or pathogenicity (n = 49), genes encoding nucleotide salvage or amino acid biosynthesis proteins (n = 6), and extrachromosomal elements (9 plasmid and 2 bacteriophage genes).

Conclusions

We have identified strain-specific differences and targets for detection of C. pneumoniae isolates from both human and animal origin. Such characterisation is necessary for an improved understanding of disease transmission and intervention.

Chlamydia pneumoniae is genetically diverse in animals and appears to have crossed the host barrier to humans on (at least) two occasions

Chlamydia pneumoniae is a common human and animal pathogen associated with a wide range of diseases. Since the first isolation of C. pneumoniae TWAR in 1965, all human isolates have been essentially clonal, providing little evolutionary insight. To address this gap, we investigated the genetic diversity of 30 isolates from diverse geographical locations, from both human and animal origin (amphibian, reptilian, equine and marsupial). Based on the level of variation that we observed at 23 discreet gene loci, it was clearly evident that the animal isolates were more diverse than the isolates of human origin. Furthermore, we show that C. pneumoniae isolates could be grouped into five major genotypes, A-E, with A, B, D and E genotypes linked by geographical location, whereas genotype C was found across multiple continents. Our evidence strongly supports two separate animal-to-human cross species transfer events in the evolutionary history of this pathogen. The C. pneumoniae human genotype identified in the USA, Canada, Taiwan, Iran, Japan, Korea and Australia (non-Indigenous) most likely originated from a single amphibian or reptilian lineage, which appears to have been previously geographically widespread. We identified a separate human lineage present in two Australian Indigenous isolates (independent geographical locations). This lineage is distinct and is present in Australian amphibians as well as a range of Australian marsupials.

Chlamydia pneumoniae is an intracellular bacterial pathogen with an extremely diverse host range (humans, amphibians, reptiles and marsupials). We selected 23 target genes in order to investigate genetic diversity: seven of these had been lost or gained by C. pneumoniae, a further six were conserved, four were polymorphic (defined by greater than 20 SNPs per 1 kbp; in this study), and six were truncated or length polymorphic in one strain or the other. Our research highlights that C. pneumoniae animal isolates are much more genetically diverse than C. pneumoniae human isolates, and have crossed the host barrier to humans on at least two occasions. Our study provides new insights into the evolution of this complex pathogen.

Genome sequence and identification of candidate vaccine antigens from the animal pathogen Dichelobacter nodosus

Dichelobacter nodosus causes ovine footrot, a disease that leads to severe economic losses in the wool and meat industries. We sequenced its 1.4-Mb genome, the smallest known genome of an anaerobe. It differs markedly from small genomes of intracellular bacteria, retaining greater biosynthetic capabilities and lacking any evidence of extensive ongoing genome reduction. Comparative genomic microarray studies and bioinformatic analysis suggested that, despite its small size, almost 20% of the genome is derived from lateral gene transfer. Most of these regions seem to be associated with virulence. Metabolic reconstruction indicated unsuspected capabilities, including carbohydrate utilization, electron transfer and several aerobic pathways. Global transcriptional profiling and bioinformatic analysis enabled the prediction of virulence factors and cell surface proteins. Screening of these proteins against ovine antisera identified eight immunogenic proteins that are candidate antigens for a cross-protective vaccine.

Complete genome sequence of the plant commensal Pseudomonas fluorescens Pf-5

Pseudomonas fluorescens Pf-5 is a plant commensal bacterium that inhabits the rhizosphere and produces secondary metabolites that suppress soilborne plant pathogens. The complete sequence of the 7.1-Mb Pf-5 genome was determined. We analyzed repeat sequences to identify genomic islands that, together with other approaches, suggested P. fluorescens Pf-5's recent lateral acquisitions include six secondary metabolite gene clusters, seven phage regions and a mobile genomic island. We identified various features that contribute to its commensal lifestyle on plants, including broad catabolic and transport capabilities for utilizing plant-derived compounds, the apparent ability to use a diversity of iron siderophores, detoxification systems to protect from oxidative stress, and the lack of a type III secretion system and toxins found in related pathogens. In addition to six known secondary metabolites produced by P. fluorescens Pf-5, three novel secondary metabolite biosynthesis gene clusters were also identified that may contribute to the biocontrol properties of P. fluorescens Pf-5.

Comparison of the genome of the oral pathogen Treponema denticola with other spirochete genomes

We present the complete 2,843,201-bp genome sequence of Treponema denticola (ATCC 35405) an oral spirochete associated with periodontal disease. Analysis of the T. denticola genome reveals factors mediating coaggregation, cell signaling, stress protection, and other competitive and cooperative measures, consistent with its pathogenic nature and lifestyle within the mixed-species environment of subgingival dental plaque. Comparisons with previously sequenced spirochete genomes revealed specific factors contributing to differences and similarities in spirochete physiology as well as pathogenic potential. The T. denticola genome is considerably larger in size than the genome of the related syphilis-causing spirochete Treponema pallidum. The differences in gene content appear to be attributable to a combination of three phenomena: genome reduction, lineage-specific expansions, and horizontal gene transfer. Genes lost due to reductive evolution appear to be largely involved in metabolism and transport, whereas some of the genes that have arisen due to lineage-specific expansions are implicated in various pathogenic interactions, and genes acquired via horizontal gene transfer are largely phage-related or of unknown function.

Genome sequence of Chlamydophila caviae (Chlamydia psittaci GPIC): examining the role of niche-specific genes in the evolution of the Chlamydiaceae

The genome of Chlamydophila caviae (formerly Chlamydia psittaci, GPIC isolate) (1 173 390 nt with a plasmid of 7966 nt) was determined, representing the fourth species with a complete genome sequence from the Chlamydiaceae family of obligate intracellular bacterial pathogens. Of 1009 annotated genes, 798 were conserved in all three other completed Chlamydiaceae genomes. The C.caviae genome contains 68 genes that lack orthologs in any other completed chlamydial genomes, including tryptophan and thiamine biosynthesis determinants and a ribose-phosphate pyrophosphokinase, the product of the prsA gene. Notable amongst these was a novel member of the virulence-associated invasin/intimin family (IIF) of Gram-negative bacteria. Intriguingly, two authentic frameshift mutations in the ORF indicate that this gene is not functional. Many of the unique genes are found in the replication termination region (RTR or plasticity zone), an area of frequent symmetrical inversion events around the replication terminus shown to be a hotspot for genome variation in previous genome sequencing studies. In C.caviae, the RTR includes several loci of particular interest including a large toxin gene and evidence of ancestral insertion(s) of a bacteriophage. This toxin gene, not present in Chlamydia pneumoniae, is a member of the YopT effector family of type III-secreted cysteine proteases. One gene cluster (guaBA-add) in the RTR is much more similar to orthologs in Chlamydia muridarum than those in the phylogenetically closest species C.pneumoniae, suggesting the possibility of horizontal transfer of genes between the rodent-associated Chlamydiae. With most genes observed in the other chlamydial genomes represented, C.caviae provides a good model for the Chlamydiaceae and a point of comparison against the human atherosclerosis-associated C.pneumoniae. This crucial addition to the set of completed Chlamydiaceae genome sequences is enabling dissection of the roles played by niche-specific genes in these important bacterial pathogens.

Role of mobile DNA in the evolution of vancomycin-resistant Enterococcus faecalis

The complete genome sequence of Enterococcus faecalis V583, a vancomycin-resistant clinical isolate, revealed that more than a quarter of the genome consists of probable mobile or foreign DNA. One of the predicted mobile elements is a previously unknown vanB vancomycin-resistance conjugative transposon. Three plasmids were identified, including two pheromone-sensing conjugative plasmids, one encoding a previously undescribed pheromone inhibitor. The apparent propensity for the incorporation of mobile elements probably contributed to the rapid acquisition and dissemination of drug resistance in the enterococci.