The Chlamydia psittaci genome: a comparative analysis of intracellular pathogens

Background

Chlamydiaceae are a family of obligate intracellular pathogens causing a wide range of diseases in animals and humans, and facing unique evolutionary constraints not encountered by free-living prokaryotes. To investigate genomic aspects of infection, virulence and host preference we have sequenced Chlamydia psittaci, the pathogenic agent of ornithosis.

Results

A comparison of the genome of the avian Chlamydia psittaci isolate 6BC with the genomes of other chlamydial species, C. trachomatis, C. muridarum, C. pneumoniae, C. abortus, C. felis and C. caviae, revealed a high level of sequence conservation and synteny across taxa, with the major exception of the human pathogen C. trachomatis. Important differences manifest in the polymorphic membrane protein family specific for the Chlamydiae and in the highly variable chlamydial plasticity zone. We identified a number of psittaci-specific polymorphic membrane proteins of the G family that may be related to differences in host-range and/or virulence as compared to closely related Chlamydiaceae. We calculated non-synonymous to synonymous substitution rate ratios for pairs of orthologous genes to identify putative targets of adaptive evolution and predicted type III secreted effector proteins.

Conclusions

This study is the first detailed analysis of the Chlamydia psittaci genome sequence. It provides insights in the genome architecture of C. psittaci and proposes a number of novel candidate genes mostly of yet unknown function that may be important for pathogen-host interactions.     

Genome sequences of Chlamydia trachomatis MoPn and Chlamydia pneumoniae AR39

The genome sequences of Chlamydia trachomatis mouse pneumonitis (MoPn) strain Nigg (1 069 412 nt) and Chlamydia pneumoniae strain AR39 (1 229 853 nt) were determined using a random shotgun strategy. The MoPn genome exhibited a general conservation of gene order and content with the previously sequenced C.trachomatis serovar D. Differences between C.trachomatis strains were focused on an ~50 kb ‘plasticity zone’ near the termination origins. In this region MoPn contained three copies of a novel gene encoding a >3000 amino acid toxin homologous to a predicted toxin from Escherichia coli 0157:H7 but had apparently lost the tryptophan biosyntheis genes found in serovar D in this region. The C.pneumoniae AR39 chromosome was >99.9% identical to the previously sequenced C.pneumoniae CWL029 genome, however, comparative analysis identified an invertible DNA segment upstream of the uridine kinase gene which was in different orientations in the two genomes. AR39 also contained a novel 4524 nt circular single-stranded (ss)DNA bacteriophage, the first time a virus has been reported infecting C.pneumoniae. Although the chlamydial genomes were highly conserved, there were intriguing differences in key nucleotide salvage pathways: C.pneumoniae has a uridine kinase gene for dUTP production, MoPn has a uracil phosphororibosyl transferase, while C.trachomatis serovar D contains neither gene. Chromosomal comparison revealed that there had been multiple large inversion events since the species divergence of C.trachomatis and C.pneumoniae, apparently oriented around the axis of the origin of replication and the termination region. The striking synteny of the Chlamydia genomes and prevalence of tandemly duplicated genes are evidence of minimal chromosome rearrangement and foreign gene uptake, presumably owing to the ecological isolation of the obligate intracellular parasites. In the absence of genetic analysis, comparative genomics will continue to provide insight into the virulence mechanisms of these important human pathogens.     

Intrastrain and interstrain genetic variation within a paralogous gene family in Chlamydia pneumoniae

Background

Chlamydia pneumoniae causes human respiratory diseases and has recently been associated with atherosclerosis. Analysis of the three recently published C. pneumoniae genomes has led to the identification of a new gene family (the Cpn 1054 family) that consists of 11 predicted genes and gene fragments. Each member encodes a polypeptide with a hydrophobic domain characteristic of proteins localized to the inclusion membrane.

Results

Comparative analysis of this gene family within the published genome sequences provided evidence that multiple levels of genetic variation are evident within this single collection of paralogous genes. Frameshift mutations are found that result in both truncated gene products and pseudogenes that vary among isolates. Several genes in this family contain polycytosine (polyC) tracts either upstream or within the terminal 5' end of the predicted coding sequence. The length of the polyC stretch varies between paralogous genes and within single genes in the three genomes. Sequence analysis of genomic DNA from a collection of 12 C. pneumoniae clinical isolates was used to determine the extent of the variation in the Cpn 1054 gene family.

Conclusions

These studies demonstrate that sequence variability is present both among strains and within strains at several of the loci. In particular, changes in the length of the polyC tract associated with the different Cpn 1054 gene family members are common within each tested C. pneumoniae isolate. The variability identified within this newly described gene family may modulate either phase or antigenic variation and subsequent physiologic diversity within a C. pneumoniae population.     

Genome Reduction Promotes Increase in Protein Functional Complexity in Bacteria

Obligate pathogenic and endosymbiotic bacteria typically experience gene loss due to functional redundancy, asexuality, and genetic drift. We hypothesize that reduced genomes increase their functional complexity through protein multitasking, in which many genes adopt new roles to counteract gene loss. Comparisons of interaction networks among six bacteria that have varied genome sizes (Mycoplasma pneumoniae, Treponema pallidum, Helicobacter pylori, Campylobacter jejuni, Synechocystis sp., and Mycobacterium tuberculosis) reveal that proteins in small genomes interact with proteins from a wider range of functions than do their orthologs in larger genomes. This suggests that surviving proteins form increasingly complex functional relationships to compensate for genes that are lost.     

Genome sequence of the cat pathogen, Chlamydophila felis.

Chlamydophila felis (Chlamydia psittaci feline pneumonitis agent) is a worldwide spread pathogen for pneumonia and conjunctivitis in cats. Herein, we determined the entire genomic DNA sequence of the Japanese C. felis strain Fe/C-56 to understand the mechanism of diseases caused by this pathogen. The C. felis genome is composed of a circular 1,166,239 bp chromosome encoding 1005 protein-coding genes and a 7552 bp circular plasmid. Comparison of C. felis gene contents with other Chlamydia species shows that 795 genes are common in the family Chlamydiaceae species and 47 genes are specific to C. felis. Phylogenetic analysis of the common genes reveals that most of the orthologue sets exhibit a similar divergent pattern but 14 C. felis genes accumulate more mutations, implicating that these genes may be involved in the evolutional adaptation to the C. felis-specific niche. Gene distribution and orthologue analyses reveal that two distinctive regions, i.e. the plasticity zone and frequently gene-translocated regions (FGRs), may play important but different roles for chlamydial genome evolution. The genomic DNA sequence of C. felis provides information for comprehension of diseases and elucidation of the chlamydial evolution.     

Restriction of Chlamydia pneumoniae replication in human dendritic cell by activation of indoleamine 2,3-dioxygenase

Infection with Chlamydia pneumoniae, a human respiratory pathogen, has been associated with various chronic diseases such as asthma and atherosclerosis, possibly because the pathogen can exist in a persistent form. C. pneumoniae persistently infect DCs in a TNF-α dependent manner. The present study investigated whether C. pneumoniae infection can induce indoleamine 2,3-dioxygenase (IDO) activity in dendritic cells, and whether the restriction of chlamydial growth in the DCs by TNF-α is IDO dependent. Our data indicate that infection of DCs with C. pneumoniae resulted in the induction of IDO expression. Reporting on our use of anti-TNF-α antibody adalimumab and varying concentrations of TNF-α, we further demonstrate that IDO induction following infection of DCs with C. pneumoniae is TNF-α dependent. The anti-chlamydial activity induced by TNF-α and the expression of chlamydial 16S rRNA gene, euo, groEL1, ftsk and tal genes were correlated with induction of IDO. Addition of excess amounts of tryptophan to the DC cultures resulted in abrogation of the TNF-α-mediated chlamydial growth restriction. These findings suggest that infection of DCs by C. pneumoniae induces production of functional IDO, which subsequently causes depletion of tryptophan. This may represent a potential mechanism for DCs to restrict bacterial growth in chlamydial infections.     

Characterization of Chlamydia MurC-Ddl, a fusion protein exhibiting D-alanyl-D-alanine ligase activity involved in peptidoglycan synthesis and D-cycloserine sensitivity

Recent characterization of chlamydial genes encoding functional peptidoglycan (PG)-synthesis proteins suggests that the Chlamydiaceae possess the ability to synthesize PG yet biochemical evidence for the synthesis of PG has yet to be demonstrated. The presence of D-amino acids in PG is a hallmark of bacteria. Chlamydiaceae do not appear to encode amino acid racemases however, a D-alanyl-D-alanine (D-Ala-D-Ala) ligase homologue (Ddl) is encoded in the genome. Thus, we undertook a genetics-based approach to demonstrate and characterize the D-Ala-D-Ala ligase activity of chlamydial Ddl, a protein encoded as a fusion with MurC. The full-length murC-ddl fusion gene from Chlamydia trachomatis serovar L2 was cloned and placed under the control of the arabinose-inducible ara promoter and transformed into a D-Ala-D-Ala ligase auxotroph of Escherichia coli possessing deletions of both the ddlA and ddlB genes. Viability of the E. coliDeltaddlADeltaddlB mutant in the absence of exogenous D-Ala-D-Ala dipeptide became dependent on the expression of the chlamydial murC-ddl thus demonstrating functional ligase activity. Domain mapping of the full-length fusion protein and site-directed mutagenesis of the MurC domain revealed that the structure of the full fusion protein but not MurC enzymatic activity was required for ligase activity in vivo. Recombinant MurC-Ddl exhibited substrate specificity for D-Ala. Chlamydia growth is inhibited by D-cycloserine (DCS) and in vitro analysis provided evidence for the chlamydial MurC-Ddl as the target for DCS sensitivity. In vivo sensitivity to DCS could be reversed by addition of exogenous D-Ala and D-Ala-D-Ala. Together, these findings further support our hypothesis that PG is synthesized by members of the Chlamydiaceae family and suggest that D-amino acids, specifically D-Ala, are present in chlamydial PG.     

Limited microsynteny between the genomes of Pristionchus pacificus and Caenorhabditis elegans

Nematodes are an attractive group of organisms for studying the evolution of developmental processes. Pristionchus pacificus was established as a satellite organism for comparing vulva development and other processes to Caenorhabditis elegans. The generation of a genetic linkage map of P.pacificus has provided a first insight into the structure and organization of the genome of this species. Pristionchus pacificus and C.elegans are separated from one another by >100 000 000 years such that the structure of the genomes of these two nematodes might differ substantially. To evaluate the amount of synteny between the two genomes, we have obtained 126 kb of continuous genomic sequence of P.pacificus, flanking the developmental patterning gene pal-1. Of the 20 predicted open reading frames in this interval, 11 have C.elegans orthologs. Ten of these 11 orthologs are located on C.elegans chromosome III, indicating the existence of synteny. However, most of these genes are distributed over a 12 Mb interval of the C.elegans genome and only three pairs of genes show microsynteny. Thus, intrachromosomal rearrange ments occur frequently in nematodes, limiting the likelihood of identifying orthologous genes of P.pacificus and C.elegans based on positional information within the two genomes.     

The ac53, ac78, ac101, and ac103 Genes Are Newly Discovered Core Genes in the Family Baculoviridae

The family Baculoviridae is a large group of insect viruses containing circular double-stranded DNA genomes of 80 to 180 kbp, which have broad biotechnological applications. A key feature to understand and manipulate them is the recognition of orthology. However, the differences in gene contents and evolutionary distances among the known members of this family make it difficult to assign sequence orthology. In this study, the genome sequences of 58 baculoviruses were analyzed, with the aim to detect previously undescribed core genes because of their remote homology. A routine based on Multi PSI-Blast/tBlastN and Multi HaMStR allowed us to detect 31 of 33 accepted core genes and 4 orthologous sequences in the Baculoviridae which were not described previously. Our results show that the ac53, ac78, ac101 (p40), and ac103 (p48) genes have orthologs in all genomes and should be considered core genes. Accordingly, there are 37 orthologous genes in the family Baculoviridae.     

Interaction between components of the type III secretion system of Chlamydiaceae

Members of the family Chlamydiaceae possess at least 13 genes, distributed throughout the chromosome, that are homologous with genes of known type III secretion systems (TTS). The aim of this study was to use putative TTS proteins of Chlamydophila pneumoniae, whose equivalents in other bacterial TTS function as chaperones, to identify interactions between chlamydial proteins. Using the BacterioMatch Two-Hybrid Vector system (Stratagene, La Jolla, Calif.), lcrH-2 and sycE, positions 1021 and 0325, respectively, from C. pneumoniae CM-1 were used as "bait" to identify target genes (positions 0324, 0705, 0708, 0808 to 0810, 1016 to 1020, and 1022) in close proximity on the chromosome. Interaction between the products of the lcrH-2 (1021) and lcrE (copN) (0324) genes was detected and confirmed by pull-down experiments and enzyme immunoassays using recombinant LcrH-2 and LcrE. As further confirmation of this interaction, the homologous genes from Chlamydia trachomatis, serovar E, and Chlamydophila psittaci, Texas turkey, were also cloned in the two-hybrid system to determine if LcrH-2 and LcrE would interact with their orthologs in other species. Consistent with their genetic relatedness, LcrH-2 from C. pneumoniae interacted with LcrE produced from the three species of Chlamydiaceae; LcrH-2 from C. psittaci reacted with LcrE from C. pneumoniae but not from C. trachomatis; and C. trachomatis LcrH-2 did not react with LcrE from the other two species. Deletions from the N and C termini of LcrE from C. pneumoniae identified the 50 C-terminal amino acids as essential for the interaction with LcrH-2. Thus, it appears that in the Chlamydiaceae TTS, LcrH-2 interacts with LcrE, and therefore it may serve as a chaperone for this protein.     

Chlamydia trachomatis Tarp Harbors Distinct G and F Actin Binding Domains That Bundle Actin Filaments

All species of Chlamydia undergo a unique developmental cycle that transitions between extracellular and intracellular environments and requires the capacity to invade new cells for dissemination. A chlamydial protein called Tarp has been shown to nucleate actin in vitro and is implicated in bacterial entry into human cells. Colocalization studies of ectopically expressed enhanced green fluorescent protein (EGFP)-Tarp indicate that actin filament recruitment is restricted to the C-terminal half of the effector protein. Actin filaments are presumably associated with Tarp via an actin binding alpha helix that is also required for actin nucleation in vitro, but this has not been investigated. Tarp orthologs from C. pneumoniae, C. muridarum, and C. caviae harbor between 1 and 4 actin binding domains located in the C-terminal half of the protein, but C. trachomatis serovar L2 has only one characterized domain. In this work, we examined the effects of domain-specific mutations on actin filament colocalization with EGFP-Tarp. We now demonstrate that actin filament colocalization with Tarp is dependent on two novel F-actin binding domains that endow the Tarp effector with actin-bundling activity. Furthermore, Tarp-mediated actin bundling did not require actin nucleation, as the ability to bundle actin filaments was observed in mutant Tarp proteins deficient in actin nucleation. These data shed molecular insight on the complex cytoskeletal rearrangements required for C. trachomatis entry into host cells.     

Complete Genome Sequences of Elephant Endotheliotropic Herpesviruses 1A and 1B Determined Directly from Fatal Cases

A highly lethal hemorrhagic disease associated with infection by elephant endotheliotropic herpesvirus (EEHV) poses a severe threat to Asian elephant husbandry. We have used high-throughput methods to sequence the genomes of the two genotypes that are involved in most fatalities, namely, EEHV1A and EEHV1B (species Elephantid herpesvirus 1, genus Proboscivirus, subfamily Betaherpesvirinae, family Herpesviridae). The sequences were determined from postmortem tissue samples, despite the data containing tiny proportions of viral reads among reads from a host for which the genome sequence was not available. The EEHV1A genome is 180,421 bp in size and consists of a unique sequence (174,601 bp) flanked by a terminal direct repeat (2,910 bp). The genome contains 116 predicted protein-coding genes, of which six are fragmented, and seven paralogous gene families are present. The EEHV1B genome is very similar to that of EEHV1A in structure, size, and gene layout. Half of the EEHV1A genes lack orthologs in other members of subfamily Betaherpesvirinae, such as human cytomegalovirus (genus Cytomegalovirus) and human herpesvirus 6A (genus Roseolovirus). Notable among these are 23 genes encoding type 3 membrane proteins containing seven transmembrane domains (the 7TM family) and seven genes encoding related type 2 membrane proteins (the EE50 family). The EE50 family appears to be under intense evolutionary selection, as it is highly diverged between the two genotypes, exhibits evidence of sequence duplications or deletions, and contains several fragmented genes. The availability of the genome sequences will facilitate future research on the epidemiology, pathogenesis, diagnosis, and treatment of EEHV-associated disease.     

Processing of Chlamydia abortus Polymorphic Membrane Protein 18D during the Chlamydial Developmental Cycle

Background

Chlamydia possess a unique family of autotransporter proteins known as the Polymorphic membrane proteins (Pmps). While the total number of pmp genes varies between Chlamydia species, all encode a single pmpD gene. In both Chlamydia trachomatis (C. trachomatis) and C. pneumoniae, the PmpD protein is proteolytically cleaved on the cell surface. The current study was carried out to determine the cleavage patterns of the PmpD protein in the animal pathogen C. abortus (termed Pmp18D).

Methodology/Principal Findings

Using antibodies directed against different regions of Pmp18D, proteomic techniques revealed that the mature protein was cleaved on the cell surface, resulting in a100 kDa N-terminal product and a 60 kDa carboxy-terminal protein. The N-terminal protein was further processed into 84, 76 and 73 kDa products. Clustering analysis resolved PmpD proteins into three distinct clades with C. abortus Pmp18D, being most similar to those originating from C. psittaci, C. felis and C. caviae.

Conclusions/Significance

This study indicates that C. abortus Pmp18D is proteolytically processed at the cell surface similar to the proteins of C. trachomatis and C. pneumoniae. However, patterns of cleavage are species-specific, with low sequence conservation of PmpD across the genus. The absence of conserved domains indicates that the function of the PmpD molecule in chlamydia remains to be elucidated.     

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.     

Comparative genome analysis of Mycoplasma pneumoniae

Background

Mycoplasma pneumoniae is a common pathogen that causes upper and lower respiratory tract infections in people of all ages, responsible for up to 40 % of community-acquired pneumonias. It also causes a wide array of extrapulmonary infections and autoimmune phenomena. Phylogenetic studies of the organism have been generally restricted to specific genes or regions of the genome, because whole genome sequencing has been completed for only 4 strains. To better understand the physiology and pathogenicity of this important human pathogen, we performed comparative genomic analysis of 15 strains of M. pneumoniae that were isolated between the 1940s to 2009 from respiratory specimens and cerebrospinal fluid originating from the USA, China and England.

Results

Illumina MiSeq whole genome sequencing was performed on the 15 strains and all genome sequences were completed. Results from the comparative genomic analysis indicate that although about 1500 SNP and indel variants exist between type1 and type 2 strains, there is an overall high degree of sequence similarity among the strains (>99 % identical to each other). Within the two subtypes, conservation of most genes, including the CARDS toxin gene and arginine deiminase genes, was observed. The major variation occurs in the P1 and ORF6 genes associated with the adhesin complex. Multiple hsdS genes (encodes S subunit of type I restriction enzyme) with variable tandem repeat copy numbers were found in all 15 genomes.

Conclusions

These data indicate that despite conclusions drawn from 16S rRNA sequences suggesting rapid evolution, the M. pneumoniae genome is extraordinarily stable over time and geographic distance across the globe with a striking lack of evidence of horizontal gene transfer.

Electronic supplementary material

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

What It Takes to Be a Pseudomonas aeruginosa? The Core Genome of the Opportunistic Pathogen Updated

Pseudomonas aeruginosa is an opportunistic bacterial pathogen able to thrive in highly diverse ecological niches and to infect compromised patients. Its genome exhibits a mosaic structure composed of a core genome into which accessory genes are inserted en bloc at specific sites. The size and the content of the core genome are open for debate as their estimation depends on the set of genomes considered and the pipeline of gene detection and clustering. Here, we redefined the size and the content of the core genome of P. aeruginosa from fully re-analyzed genomes of 17 reference strains. After the optimization of gene detection and clustering parameters, the core genome was defined at 5,233 orthologs, which represented ~ 88% of the average genome. Extrapolation indicated that our panel was suitable to estimate the core genome that will remain constant even if new genomes are added. The core genome contained resistance determinants to the major antibiotic families as well as most metabolic, respiratory, and virulence genes. Although some virulence genes were accessory, they often related to conserved biological functions. Long-standing prophage elements were subjected to a genetic drift to eventually display a G+C content as higher as that of the core genome. This contrasts with the low G+C content of highly conserved ribosomal genes. The conservation of metabolic and respiratory genes could guarantee the ability of the species to thrive on a variety of carbon sources for energy in aerobiosis and anaerobiosis. Virtually all the strains, of environmental or clinical origin, have the complete toolkit to become resistant to the major antipseudomonal compounds and possess basic pathogenic mechanisms to infect humans. The knowledge of the genes shared by the majority of the P. aeruginosa isolates is a prerequisite for designing effective therapeutics to combat the wide variety of human infections.     

Comparison of whole genome sequences of Chlamydia pneumoniae J138 from Japan and CWL029 from USA

Chlamydia pneumoniae is a widespread pathogen of humans causing pneumonia and bronchitis. There are many reports of an association between C.pneumoniae infection and atherosclerosis. We determined the whole genome sequence of C.pneumoniae strain J138 isolated in Japan in 1994 and compared it with the sequence of strain CWL029 isolated in the USA before 1987. The J138 circular chromosome consists of 1 226 565 nt (40.7% G+C) with 1072 likely protein-coding genes that is 3665 nt shorter than the CWL029 genome. Plasmids, phage- or transposon-like sequences were not identified. The overall genomic organization, gene order and predicted proteomes of the two strains are very similar, suggesting a high level of structural and functional conservation between the two unrelated isolates. The most conspicuous differences in the J138 genome relative to the CWL029 genome are the absence of five DNA segments, ranging in size from 89 to 1649 nt, and the presence of three DNA segments, ranging from 27 to 84 nt. The complex organization of these ‘different zones’ may be attributable to a unique system of recombination.