Genome sequencing and comparative analysis of three Chlamydia pecorum strains associated with different pathogenic outcomes
© Sait et al.; licensee BioMed Central Ltd. 2014
Received: 24 September 2013
Accepted: 6 January 2014
Published: 14 January 2014
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.
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.
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.
KeywordsChlamydia pecorum Genome sequence Polymorphic membrane proteins Plasticity zone Tryptophan metabolism Folate biosynthesis Clustered tandem repeats
Members of the genus Chlamydia are Gram-negative, obligate intracellular pathogens that share a biphasic developmental cycle. Chlamydia pecorum infects a broad host range, including small and large ruminants, swine, birds and marsupials. Seroprevalence and PCR-based studies suggest that infection or exposure to C. pecorum and/or C. abortus is almost ubiquitous in cattle and sheep [1–5]. In the majority of these cases, infection is subclinical, with C. pecorum being routinely detected in the intestine and genital tract. The incidence and severity of disease caused by C. pecorum appears to be heightened in koalas and is associated with clinical disease such as conjunctivitis, urinary- and reproductive tract disease, and infertility . Many chlamydial species, including C. pecorum can enter persistent states, characterised in vitro by enlarged, morphologically aberrant, non-fusogenic reticulate bodies (RBs). Persistence can be induced in vitro by antibiotic exposure , amino acid-  or iron-  deficiencies and exposure to IFN-γ  and it is likely that C. pecorum causes a persistent, subclinical infection in the host. Subclinical infections can have detrimental effects on the animal’s health. Animals with inapparent chlamydiae infections have higher body temperatures, lower body weights, reduced growth rates, reduced iron, haemoglobin, haematocrit and leukocyte levels and a higher incidence of follicular bronchiolitis [11–13]. C. pecorum can also cause clinical disease including encephalomyelitis, vaginitis, endometritis, mastitis, conjunctivitis, polyarthritis, pneumonia, enteritis, orchitis, pleuritis, infertility or pericarditis .
Genetic variation has been reported to occur between C. pecorum strains in ompA, the rrn-nqr F intergenic region, incA, rRNAs, a number of housekeeping genes and the hypothetical protein ORF663 [14–22]. These and other unidentified genomic differences may enable differentiation between strains isolated from asymptomatic or diseased animals. However, to date, only the genome sequence of a single C. pecorum strain (E58) has been published . The genetic factors responsible for the diverse host range, tissue tropism, disease outcomes and associated sequelae of C. pecorum infections are thus still poorly understood. In this study, we present the complete genome sequences of three C. pecorum strains isolated from animals exhibiting different disease manifestations and use comparative genomics to provide insights into the biology of C. pecorum and to identify both genus- and species-specific virulence factors.
Results and discussion
Genome features and comparative analysis
General features of C. pecorum PV3056/3, P787 and W73 compared with the type strain E58 [CP002608]
C. pecorum(E58) []
Date of isolation
Cow, cervical swab
Sheep, synovial fluid
Genome size (bp)
% GC of genome
Predicted no. of pseudogenes
No. of CDS with functional assignment
No. of pmp proteins
No. of tRNA genes
No. of rRNA operons
No. of sRNA molecules
Location of OriC region
1104229–147 nt (471 nt)
1106263–147 nt (296 nt)
1106389–147 nt (295 nt)
305941–306155 nt (215 nt)
Comparative genomics of chlamydial species has identified a number of genes coding for metabolic functions, such as tryptophan metabolism, biotin biosynthesis and folate biosynthesis, where subtle variations in gene content may contribute to growth of the organism in vivo and the ability to evade the host immune system [23, 25–29].
The 3 sequenced C. pecorum strains and E58 contain the biotin biosynthesis operon encoding bioBFDA. (CPE1_0687-CPE1_0690; CPE2_0688-CPE2_0691; CPE3_0688-CPE3_0691). This region shows significant variability between chlamydial species, being absent in C. caviae, C. trachomatis and C. muridarum but present in C. abortus, C. psittaci, C. felis and C. pneumoniae. The ability to synthesise biotin is hypothesised to assist in the colonization of biotin-limited niches and contribute to the tissue tropism differences observed in the chlamydiae . Upstream of bioBFDA, located between dapB and bioB, a series of genes encoding hypothetical proteins with unknown function and limited distribution across chlamydial species are present. Chlamydia abortus, C. psittaci and C. felis genomes contain four genes (in C. abortus, CAB681, CAB682, CAB683 and CAB684), C. pneumoniae contains 2 genes in this region that are homologues of CAB681 and CAB682, while C. pecorum contains one gene in this region that is homologous to CAB684 (Additional file 2: Figure S1).
Bacterial secretion systems
The Type III secretion system (T3SS) consists of 18 genes encoding the major structural components of the secretion apparatus, accessory proteins and chaperones and is arranged in 4 genetic loci (Additional file 1: Table S5). In sequenced chlamydial genomes, a number of putative T3SS effector proteins belonging to the Inc or transmembrane head (TMH) protein family are located in the region extending between pmpD and lpxB (Additional file 2: Figure S3). The distance between the 3′ ends of these genes in C. pecorum is ~2.8 kb (2 genes) compared to 18.1 kb C. abortus (11 genes), 17.7 kb in C. psittaci (11 genes), 16.4 kb in C. caviae (13 genes), 15.9 kb in C. felis (11 genes) and 1.7 kb in C. pneumoniae (1 gene). The two genes present in this region in C. pecorum (CPE1_0764 (pseudogene), CPE1_0765, CPE2_0765, CPE2_0766, CPE3_0765, CPE3_0766) possess an N-terminal signal sequence, a single N-terminal transmembrane domain and two domains of unknown function (DUF1539 and DUF1548). Members of this protein family are present in C. abortus, C. psittaci, C. caviae and C. felis (3 CDSs each) and C. pneumoniae (1 CDS).
Simple sequence repeats
A region of variability between C. pecorum and other chlamydial species is located immediately upstream of the 5S rRNA gene. This region, between the 3′ ends of the 5S rRNA and nqrF genes range in size from 261–269 bp in C. pecorum to 4464 bp in C. caviae. In C. caviae this region encodes a 1291aa residue pseudogene identified as a member of the virulence-associated invasion/intimin family of outer membrane proteins of Gram-negative bacteria. The genome of C. abortus contains two CDSs in place of the intimin family gene in this region, encoding a conserved membrane protein and a unique hypothetical protein. In C. psittaci, C. felis and C. muridarum these two proteins are fused to encode a single hypothetical protein. In C. pecorum there are no predicted CDSs in this region and the intergenic region between the 5S rRNA and nqrF genes comprises an 8 bp simple sequence repeat sequence AAAGCACT repeated 12 (W73, PV3056/3 and E58) or 13 times (P787) (Additional file 2: Figure S4).
Clustered tandem repeat (CTR) sequences observed in orthologs of hypothetical protein ORF663
CTR motif sequence
Clustered tandem repeat sequences (CTR) observed in IncA
CTR motif sequence
In chlamydial species, the plasticity zone is defined as the region between inosine-5′-monophosphate dehydrogenase (guaB) and acetyl-CoA carboxylase (accB) and is the region of the genome that is most variable in gene content and sequence. In C. pecorum, this region is 40.3-42.1 kb in size and contains 16 (PV3056/3) or 17 (W73 and P787) genes encoding GMP synthase, an adenosine deaminase superfamily-protein, a MAC/perforin domain-containing protein, 3 (PV3056/3) or 4 (W73 and P787) phospholipase D family proteins, 2 cytotoxins and 4 hypothetical proteins (Additional file 1: Table S6).
Flanking the cytotoxin genes in C. pecorum are 4 (PV3056/3) or 5 (P787, W73 and E58) phospholipase D (PLD) genes each containing the conserved HxKx4Dx6GSxN (HKD) motif essential for the initiation of phosphodiesterase activity and amino acid motifs that are responsible for catalytic activity. PLD genes identified in the plasticity zone of P787, W73 and E58 share 95-99% amino acid sequence identity (CPE2_0554, CPE3_0554, G5S_0938; CPE2_0553, CPE3_0553, G5S_0935; CPE2_0551, CPE3_0551, G5S_0931; CPE2_0550, CPE3_0550, G5S_0930) whereas orthologous PLD genes in PV3056/3 are more divergent (58-71% sequence identity) (CPE1_0553, CPE1_0551, CPE1_0550). The remaining PLD gene is almost identical in E58 and W73 (98% identity, CPE2_0556, G5S_0945) but divergent in the remaining strains (55-79% identity, CPE1_0555, CPE3_0556). The presence of poly(G) and poly(C) homopolymeric tracts ranging in size from 5–19 nucleotides within the PLD genes and the presence of intact variants in the sequence reads of pseudogenes could indicate that these proteins are subject to phase variation by slip-strand pairing . Whilst the function of PLD in C. pecorum is currently unknown, PLD can perform numerous functions ranging from DNA hydrolysis, to protein-protein interactions with host signalling pathways, to the more classic lipase function. In C. trachomatis, PLD genes located in the PZ have been associated with inclusion formation , whereas in other bacteria PLD has been identified as an important virulence determinant involved in dissemination, serum resistance and invasion of epithelial cells [48, 49].
The complete genome sequence of C. pecorum P787, W73 and PV3056/3 was determined by Illumina/Solexa and Roche 454 genome sequencing. Despite the differences in the clinical manifestations of infections caused by the strains, comparative analysis revealed a high level of sequence conservation, gene content and order between the genomes. Additional genomic analyses of strains originating from other non-ruminant host species, such as pig and koala, will determine if the high level of sequence similarity is common to all, or just ruminant strains of C. pecorum. In agreement with previous studies , differences in the number of clustered tandem repeat sequences in ORF663 were observed between strains isolated from diseased (PV3056/3 and P787) or asymptomatic (W73) animals however, no other genetic differences were observed that may account for the different disease manifestations. A number of metabolic traits were identified in C. pecorum that may contribute to its ability to evade the host immune system and enable persistent infections to be established in the host. Specifically, this study has particularly highlighted the absence of genes involved in folate biosynthesis and the presence of tryptophan and biotin biosynthesis pathways. The presence of clustered tandem repeats in surface expressed proteins, 15 polymorphic membrane proteins, two cytotoxin genes and multiple phospholipase D genes that are likely to be subject to phase variable expression may play a role in the invasion of host cells and trigger the switching between persistent and acute disease in the host.
C. pecorum strain information, propagation and preparation of gDNA
Three C. pecorum strains originating from different geographical regions and disease manifestations were selected for genome sequencing. Strain P787 was isolated in Scotland, in 1977, from the affected synovial fluid of a sheep with polyarthritis. Strain PV3056/3 was isolated in Italy, in 1991, from a cervical swab of a cow with purulent metritis and has subsequently been shown to induce a purulent metritis following inoculation into the uterine body and cervix of cattle . Strain W73 was isolated in Northern Ireland, in 1989, from the faeces of a sheep with an inapparent enteric infection and has subsequently been found to be non-invasive in a mouse model of infection .
Strains were propagated in Caco-2 cells grown in RPMI medium supplemented with 5% FBS and 1 μg/ml cyclohexamide. Genomic DNA from PV3056/3 and P787 was derived from the 7th tissue culture passage of original strains propagated in fertile hens’ eggs. W73 was derived from the 6th tissue culture passage of a strain propagated in fertile hens’ eggs, however the passage history prior to this is unknown. Flasks of infected cells were harvested using glass beads followed by centrifugation at 22,000 × g for 40 mins. Pellets were washed in ice-cold PBS and re-centrifuged as before. Pellets were resuspended in 20 mM Tris–HCl (pH 7.5)/150 mM KCl/1% sarkosyl and lightly homogenised using a ground glass homogeniser. Homogenised cells were layered onto cushions of 15% sucrose in 20 mM Tris–HCl (pH 7.5)/150 mM KCl/1% sarkosyl and centrifuged at 70,000 × g for 45 min at 4°C. Genomic DNA was extracted from pellets using the Wizard DNA extraction kit (Promega).
Genome sequencing was performed by The Gene Pool genomic facility in The University of Edinburgh using Roche 454 GS-FLX and Solexa/Illumina 35-bp paired-end sequencing on standard libraries constructed according to the manufacturers instructions. Reads were assembled using Newbler v2 (Roche) and Velvet v.0.7 , combined using minimus2 and mapped to the reference genome of C. pecorum E58  to generate 13, 10 and 9 contigs for P787, W73 and PV3056/3 respectively. In total, 12,926,259 (PV3056/3), 8,169,259 (W73) and 10,039,539 (P787) reads obtained from Solexa/Illumina sequencing and 95,683 (PV3056/3), 101,405 (W73) and 65,050 (P787) reads from Roche 454 GS-FLX sequencing were obtained. Following quality filtering, sequencing reads were mapped to the reference genome providing approximately 253× (PV3056/3), 136× (W73) and 59.4× (P787) sequencing coverage. Regions spanning the contig ends were PCR-amplified using Phusion High-fidelity DNA polymerase (NEB) and the sequence determined ensuring that each base was covered by sequence in each direction.
Sequence annotation and analysis
Protein-encoding genes were predicted using Prodigal  and open reading frames (ORFs) consisting of fewer than 30 codons or those overlapping larger open reading frames were eliminated. Frameshifts, point mutations and pseudogenes were corrected or confirmed by visual inspection of mapped reads using Tablet . The origin of replication was determined using Ori-finder  and the genomes were adjusted so that the first base was upstream of the hemB gene in the oriC region. Ribosomal RNA genes and tRNA genes were identified using RNAmmer and ARAGORN [56, 57]. Sequences of experimentally validated small non-coding RNAs (sRNA) from chlamydia were downloaded from BSRD  and identified in C. pecorum genomes using blastn. Functional assignments were made based on homology searches using blastp  against protein sequences present in the NCBI nr database and the identification of conserved domains using Pfam  and InterProScan protein databases . Signal sequences were predicted using the LipoP 1.0 . KEGG orthology assignments were performed using KAAS . Data collation and annotation was performed using Artemis .
Comparative analysis were performed using the following genomes: C. pecorum E58 [GenBank: CP002608] , C. abortus S26/3 [GenBank: CR848038] , C. caviae GPIC [GenBank: AE015925] , C. felis Fe/C-56 [GenBank: AP006861] , C. psittaci 6BC [GenBank: CP002586] , C trachomatis D/UW-3/CX [GenBank: AE001273] ), C. pneumoniae AR39 [GenBank: AE002161] and C. muridarum Nigg [GenBank: AE002160] . Global genomic comparisons were visualised using ACT  with input files generated by the tblastx function in DoubleAct http://www.hpa-bioinfotools.org.uk/pise/double_act.html# with a cutoff score of 0. Comparisons of regions flanking the PZ were performed using default blastn settings in EasyFig . Orthologous gene sets were identified by OrthoMCL-DB using reciprocal blastp with a cutoff of e-5 and 50% match . Genome maps were generated using the CGView Server .
Reference sequences were obtained from GenBank and aligned with relevant C. pecorum CDSs using MUSCLE . Phylogenetic alignments and tree files are available from the Dryad Digital repository http://doi.org/10.5061/dryad.np597. For ribosomal proteins, 48 individual alignments were concatenated into a single alignment for analysis. For the phylogenetic analysis of cytotoxin genes, GBlocks v 0.91  was used to eliminate regions that could not be unambiguously aligned resulting in 2845 positions being analysed (75% of the original 3766 positions). Phylogenetic analyses were performed using PhyML (for ribosomal proteins and polymorphic membrane proteins) or MrBayes (for cytotoxins) software  launched from the TOPALi v2.5package  generated using the JTT + G (ribosomal proteins), JTT + I + G (polymorphic membrane proteins) or WAG + I + G (cytotoxins) substitution model that was determined to be the model of best fit based on the BIC criterion. For MrBayes phylogeny, trees were generated using Markov Chain Monte Carlo (MCMC) settings of 2 runs of 625,000 generations with a burn-in of 125,000 generations with trees sampled every 100 runs. For PhyML phylogeny, bootstrap analysis was performed based on 100 replicate trees. Phylogenetic network analysis was performed using SplitsTree .
Nucleotide sequence accession number
Genome sequences of C. pecorum strains PV3056/3, W73 and P787 have been deposited in GenBank under the accession numbers CP004033, CP004034, and CP004035, respectively.
This work as well as MS, NW, LS and EMC was funded by grant no. BB/E018939/1 from the Biotechnology and Biological Sciences Research Council (BBSRC) and by the Scottish Government Rural and Environment Science and Analytical Services division (RESAS).
- Cavirani S, Cabassi CS, Donofrio G, De Iaco B, Taddei S, Flammini CF: Association between Chlamydia psittaci seropositivity and abortion in Italian dairy cows. Prev Vet Med. 2001, 50: 145-151. 10.1016/S0167-5877(01)00197-0.PubMedView ArticleGoogle Scholar
- DeGraves FJ, Gao D, Hehnen HR, Schlapp T, Kaltenboeck B: Quantitative detection of Chlamydia psittaci and C. pecorum by high-sensitivity real-time PCR reveals high prevalence of vaginal infection in cattle. J Clin Microbiol. 2003, 41: 1726-1729. 10.1128/JCM.41.4.1726-1729.2003.PubMed CentralPubMedView ArticleGoogle Scholar
- Jee J, Degraves FJ, Kim T, Kaltenboeck B: High prevalence of natural Chlamydophila species infection in calves. J Clin Microbiol. 2004, 42: 5664-5672. 10.1128/JCM.42.12.5664-5672.2004.PubMed CentralPubMedView ArticleGoogle Scholar
- Lenzko H, Moog U, Henning K, Lederbach R, Diller R, Menge C, Sachse K, Sprague LD: High frequency of chlamydial co-infections in clinically healthy sheep flocks. BMC Vet Res. 2011, 7: 29-10.1186/1746-6148-7-29.PubMed CentralPubMedView ArticleGoogle Scholar
- Wilson K, Sammin D, Harmeyer S, Nath M, Livingstone M, Longbottom D: Seroprevalence of chlamydial infection in cattle in Ireland. Vet J. 2012, 193: 583-585. 10.1016/j.tvjl.2011.12.018.PubMedView ArticleGoogle Scholar
- Yousef Mohamad K, Rodolakis A: Recent advances in the understanding of Chlamydophila pecorum infections, sixteen years after it was named as the fourth species of the Chlamydiaceae family. Vet Res. 2010, 41: 27-10.1051/vetres/2009075.PubMed CentralView ArticleGoogle Scholar
- Matsumoto A, Manire GP: Electron microscopic observations on the effects of penicillin on the morphology of Chlamydia psittaci. J Bacteriol. 1970, 101: 278-285.PubMed CentralPubMedGoogle Scholar
- Coles AM, Reynolds DJ, Harper A, Devitt A, Pearce JH: Low-nutrient induction of abnormal chlamydial development: a novel component of chlamydial pathogenesis?. FEMS Microbiol Lett. 1993, 106: 193-200. 10.1111/j.1574-6968.1993.tb05958.x.PubMedView ArticleGoogle Scholar
- Raulston JE: Response of Chlamydia trachomatis serovar E to iron restriction in vitro and evidence for iron-regulated chlamydial proteins. Infect Immun. 1997, 65: 4539-4547.PubMed CentralPubMedGoogle Scholar
- Pantoja LG, Miller RD, Ramirez JA, Molestina RE, Summersgill JT: Characterization of Chlamydia pneumoniae persistence in HEp-2 cells treated with gamma interferon. Infect Immun. 2001, 69: 7927-7932. 10.1128/IAI.69.12.7927-7932.2001.PubMed CentralPubMedView ArticleGoogle Scholar
- Reinhold P, Jaeger J, Liebler-Tenorio E, Berndt A, Bachmann R, Schubert E, Melzer F, Elschner M, Sachse K: Impact of latent infections with Chlamydophila species in young cattle. Vet J. 2008, 175: 202-211. 10.1016/j.tvjl.2007.01.004.PubMedView ArticleGoogle Scholar
- Jaeger J, Liebler-Tenorio E, Kirschvink N, Sachse K, Reinhold P: A clinically silent respiratory infection with Chlamydophila spp. in calves is associated with airway obstruction and pulmonary inflammation. Vet Res. 2007, 38: 711-728. 10.1051/vetres:2007027.PubMedView ArticleGoogle Scholar
- Poudel A, Elsasser TH, Rahman KS, Chowdhury EU, Kaltenboeck B: Asymptomatic endemic Chlamydia pecorum infections reduce growth rates in calves by up to 48 percent. PLoS One. 2012, 7: e44961-10.1371/journal.pone.0044961.PubMed CentralPubMedView ArticleGoogle Scholar
- Anderson IE, Baxter SI, Dunbar S, Rae AG, Philips HL, Clarkson MJ, Herring AJ: Analyses of the genomes of chlamydial isolates from ruminants and pigs support the adoption of the new species Chlamydia pecorum. Int J Syst Bacteriol. 1996, 46: 245-251. 10.1099/00207713-46-1-245.PubMedView ArticleGoogle Scholar
- Jackson M, Giffard P, Timms P: Outer membrane protein A gene sequencing demonstrates the polyphyletic nature of koala Chlamydia pecorum isolates. Syst Appl Microbiol. 1997, 20: 187-200. 10.1016/S0723-2020(97)80065-3.View ArticleGoogle Scholar
- Kaltenboeck B, Kousoulas KG, Storz J: Structures of and allelic diversity and relationships among the major outer membrane protein (ompA) genes of the four chlamydial species. J Bacteriol. 1993, 175: 487-502.PubMed CentralPubMedGoogle Scholar
- Fukushi H, Hirai K: Genetic diversity of avian and mammalian Chlamydia psittaci strains and relation to host origin. J Bacteriol. 1989, 171: 2850-2855.PubMed CentralPubMedGoogle Scholar
- Salinas J, Souriau A, De Sa C, Andersen AA, Rodolakis A: Serotype 2-specific antigens from ruminant strains of Chlamydia pecorum detected by monoclonal antibodies. Comp Immunol Microbiol Infect Dis. 1996, 19: 155-161. 10.1016/0147-9571(95)00029-1.PubMedView ArticleGoogle Scholar
- Liu Z, Rank R, Kaltenboeck B, Magnino S, Dean D, Burall L, Plaut RD, Read TD, Myers G, Bavoil PM: Genomic plasticity of the rrn-nqrF intergenic segment in the Chlamydiaceae. J Bacteriol. 2007, 189: 2128-2132. 10.1128/JB.00378-06.PubMed CentralPubMedView ArticleGoogle Scholar
- Yousef Mohamad K, Roche SM, Myers G, Bavoil PM, Laroucau K, Magnino S, Laurent S, Rasschaert D, Rodolakis A: Preliminary phylogenetic identification of virulent Chlamydophila pecorum strains. Infect Genet Evol. 2008, 8: 764-771. 10.1016/j.meegid.2008.06.009.PubMedView ArticleGoogle Scholar
- Yousef Mohamad K, Rekiki A, Myers G, Bavoil PM, Rodolakis A: Identification and characterisation of coding tandem repeat variants in incA gene of Chlamydophila pecorum. Vet Res. 2008, 39: 56-10.1051/vetres:2008032.PubMedView ArticleGoogle Scholar
- Jelocnik M, Frentiu FD, Timms P, Polkinghorne A: Multilocus sequence analysis provides insights into molecular epidemiology of Chlamydia pecorum infections in Australian sheep, cattle and koalas. J Clin Microbiol. 2013, 51: 2625-2632. 10.1128/JCM.00992-13.PubMed CentralPubMedView ArticleGoogle Scholar
- Mojica S, Huot Creasy H, Daugherty S, Read TD, Kim T, Kaltenboeck B, Bavoil P, Myers GS: Genome sequence of the obligate intracellular animal pathogen Chlamydia pecorum E58. J Bacteriol. 2011, 193: 3690-10.1128/JB.00454-11.PubMed CentralPubMedView ArticleGoogle Scholar
- Pannekoek Y, Dickx V, Beeckman DSA, Jolley KA, Keijzers WC, Vretou E, Maiden MCJ, Vanrompay D, van der Ende A: Multi locus sequence typing of Chlamydia reveals an association between Chlamydia psittaci genotypes and host species. Plos one. 2010, 5: e14179-10.1371/journal.pone.0014179.PubMed CentralPubMedView ArticleGoogle Scholar
- Thomson NR, Yeats C, Bell K, Holden MT, Bentley SD, Livingstone M, Cerdeño-Tárraga AM, Harris B, Doggett J, Ormond D, Mungall K, Clarke K, Feltwell T, Hance Z, Sanders M, Quail MA, Price C, Barrell BG, Parkhill J, Longbottom D: The Chlamydophila abortus genome sequence reveals an array of variable proteins that contribute to interspecies variation. Genome Res. 2005, 15: 629-640. 10.1101/gr.3684805.PubMed CentralPubMedView ArticleGoogle Scholar
- Read TD, Myers GSA, Brunham RC, Nelson WC, Paulsen IT, Heidelberg J, Holtzapple E, Khouri H, Federova NB, Carty HA, Umayam LA, Haft DH, Peterson J, Beanan MJ, White O, Salzberg SL, Hsia R-C, McClarty G, Rank RG, Bavoil PM, Fraser CM: Genome sequence of Chlamydophila caviae (Chlamydia psittaci GPIC): examining the role of niche-specific genes in the evolution of the Chlamydiaceae. Nucleic Acids Res. 2003, 31: 2134-2147. 10.1093/nar/gkg321.PubMed CentralPubMedView ArticleGoogle Scholar
- Azuma Y, Hirakawa H, Yamashita A, Cai Y, Rahman MA, Suzuki H, Mitaku S, Toh H, Goto S, Murakami T, Sugi K, Hayashi H, Fukushi H, Hattori M, Kuhara S, Shirai M: Genome sequence of the cat pathogen, Chlamydophila felis. DNA Res. 2006, 13: 15-23. 10.1093/dnares/dsi027.PubMedView ArticleGoogle Scholar
- Grinblat-Huse V, Drabek EF, Creasy HH, Daugherty SC, Jones KM, Santana-Cruz I, Tallon LJ, Read TD, Hatch TP, Bavoil P, Myers GS: Genome sequences of the zoonotic pathogens Chlamydia psittaci 6BC and Cal10. J Bacteriol. 2011, 193: 4039-4040. 10.1128/JB.05277-11.PubMed CentralPubMedView ArticleGoogle Scholar
- Stephens RS, Kalman S, Lammel C, Fan J, Marathe R, Aravind L, Mitchell W, Olinger L, Tatusov RL, Zhao Q, Koonin EV, Davis RW: Genome sequence of an obligate intracellular pathogen of humans: Chlamydia trachomatis. Science. 1998, 282: 754-759.PubMedView ArticleGoogle Scholar
- Taylor MW, Feng GS: Relationship between interferon-γ, indoleamine 2,3-dioxygenase, and tryptophan catabolism. FASEB J. 1991, 5: 2516-2522.PubMedGoogle Scholar
- Beatty WI, Belanger TA, Desai AA, Morrison RP, Byrne GI: Tryptophan depletion as a mechanism of gamma interferon-mediated chlamydial persistence. Infect Immun. 1994, 62: 3705-3711.PubMed CentralPubMedGoogle Scholar
- Eudes A, Erkens GB, Slotboom DJ, Rodionov DA, Naponelli V, Hanson AD: Identification of genes encoding the folate- and thiamine-binding membrane proteins in Firmicutes. J Bacteriol. 2008, 190: 7591-7594. 10.1128/JB.01070-08.PubMed CentralPubMedView ArticleGoogle Scholar
- Myllykallio H, Leduc D, Filee J, Liebl U: Life without dihydrofolate reductase FolA. Trends Microbiol. 2003, 11: 220-223. 10.1016/S0966-842X(03)00101-X.PubMedView ArticleGoogle Scholar
- Henderson IR, Lam AC: Polymorphic proteins of Chlamydia spp. – autotransporters beyond the Proteobacteria. Trends Microbiol. 2001, 9: 573-578. 10.1016/S0966-842X(01)02234-X.PubMedView ArticleGoogle Scholar
- Siboo IR, Chambers HF, Sullam PM: Role of SraP, a serine-rich surface protein of Staphylococcus aureus, in binding to human platelets. Infect Immun. 2005, 73: 2273-2280. 10.1128/IAI.73.4.2273-2280.2005.PubMed CentralPubMedView ArticleGoogle Scholar
- Seifert KN, Adderson EE, Whiting AA, Bohnsack JF, Crowley PJ, Brady LJ: A unique serine-rich repeat protein (Srr-2) and novel surface antigen (epsilon) associated with a virulent lineage of serotype III Stretococcus agalactiae. Microbiol. 2006, 152: 1029-1040. 10.1099/mic.0.28516-0.View ArticleGoogle Scholar
- Gravekamp C, Rosner B, Madoff LC: Deletion of repeats in the alpha C protein enhances the pathogenicity of group B streptococci in immune mice. Infect Immun. 1998, 66: 4347-4354.PubMed CentralPubMedGoogle Scholar
- Rockey DD, Grosenbach D, Hruby DE, Peacock MG, Helnzen RA, Hackstadt T: Chlamydia psittaci IncA is phosphorylated by the host cell and exposed on the cytoplasmic face of the developing inclusion. Mol Microbiol. 1997, 24: 217-228. 10.1046/j.1365-2958.1997.3371700.x.PubMedView ArticleGoogle Scholar
- Yousef Mohamad K, Rekiki A, Berri M, Rodolakis A: Recombinant 35-kDa inclusion membrane protein IncA as a candidate antigen for serodiagnosis of Chlamydophila pecorum. Vet Microbiol. 2010, 143: 424-428. 10.1016/j.vetmic.2009.11.017.View ArticleGoogle Scholar
- Hackstadt T, Scidmore-Carlson MA, Shaw EI, Fisher ER: The Chlamydia trachomatis IncA protein is required for homotypic vesicle fusion. Cell Microbiol. 1999, 1: 119-130. 10.1046/j.1462-5822.1999.00012.x.PubMedView ArticleGoogle Scholar
- Geisler WM, Suchland RJ, Rockey DD, Stamm WE: Epidemiology and clinical manifestations of unique Chlamydia trachomatis isolates that occupy nonfusogenic inclusions. J Infect Dis. 2001, 184: 879-884. 10.1086/323340.PubMedView ArticleGoogle Scholar
- Von Eichel-Streiber C, Boquet P, Sauerborn M, Thelestam M: Large clostridial cytotoxins-a family of glycosyltransferases modifying small GTP-binding proteins. Trends Microbiol. 1996, 4: 375-382. 10.1016/0966-842X(96)10061-5.PubMedView ArticleGoogle Scholar
- Klapproth JA, Scaletsky ICA, McNamara BP, Lai L, Malstrom C, James SP, Donnenberg MS: A large toxin from pathogenic Escherichia coli strains that inhibits lymphocyte activation. Infect Immun. 2000, 68: 2148-2155. 10.1128/IAI.68.4.2148-2155.2000.PubMed CentralPubMedView ArticleGoogle Scholar
- Entrican G, Brown J, Graham S: Cytokines and the protective host immune response to Chlamydia psittaci. Comp Immun Microbiol Infect Dis. 1998, 21: 15-26. 10.1016/S0147-9571(97)00020-9.View ArticleGoogle Scholar
- Shemer Y, Sarov I: Inhibition of growth of Chlamydia trachomatis by human gamma interferon. Infect Immun. 1985, 48: 592-596.PubMed CentralPubMedGoogle Scholar
- Viratyosin W, Campbell LA, Kuo CC, Rockey DD: Intrastrain and interstrain genetic variation within a paralogous gene family in Chlamydia pneumoniae. BMC Microbiol. 2002, 2: 38-10.1186/1471-2180-2-38.PubMed CentralPubMedView ArticleGoogle Scholar
- Nelson DE, Crane DD, Taylor LD, Dorward DW, Goheen MM, Caldwell HD: Inhibition of chlamydiae by primary alcohols correlates with the strain-specific complement of plasticity zone phospholipase D genes. Infect Immun. 2006, 74: 73-80. 10.1128/IAI.74.1.73-80.2006.PubMed CentralPubMedView ArticleGoogle Scholar
- Jacobs AC, Hood I, Boyd KL, Olson PD, Morrison JM, Carson S, Sayood K, Iwen PC, Skaar EP, Dunman PM: Inactivation of phospholipase D diminishes Acinetobacter baumannii pathogenesis. Infect Immun. 2010, 78: 1952-1962. 10.1128/IAI.00889-09.PubMed CentralPubMedView ArticleGoogle Scholar
- Edwards JL, Entz DD, Apicella MA: Gonococcal phospholipase D modulates the expression and function of complement receptor 3 in primary cervical epithelial cells. Infect Immun. 2003, 71: 6381-6391. 10.1128/IAI.71.11.6381-6391.2003.PubMed CentralPubMedView ArticleGoogle Scholar
- Jones GE, Machell DA, Biolatti B, Appino S: Experimental infections of the genital tract of cattle with Chlamydia psittaci and Chlamydia pecorum. Proceedings of the Ninth International Symposium on Human Chlamydial Infection. Edited by: Stevens RS, Byrne GI, Christianson G. 1998, San Francisco, 446-449.Google Scholar
- Denamur E, Sayada C, Souriau A, Orfila J, Rodolakis A, Elion J: Restriction pattern of the major outer-membrane protein gene provides evidence for a homogeneous invasive group among ruminant isolates of Chlamydia psittaci. J Gen Microbiol. 1991, 137: 2525-2530. 10.1099/00221287-137-11-2525.PubMedView ArticleGoogle Scholar
- Zerbino DR, Birney E: Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008, 18: 821-829. 10.1101/gr.074492.107.PubMed CentralPubMedView ArticleGoogle Scholar
- Hyatt D, Chen GL, Locascio PF, Land ML, Larimer FW, Hauser LJ: Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinforma. 2010, 11: 119-10.1186/1471-2105-11-119.View ArticleGoogle Scholar
- Milne I, Stephen G, Bayer M, Cock PJ, Pritchard L, Cardle L, Shaw P, Marshall D: Using tablet for visual exploration of second-generation sequencing data. Brief Bioinform. 2013, 14: 193-202. 10.1093/bib/bbs012.PubMedView ArticleGoogle Scholar
- Gao F, Zhang CT: Ori-Finder: a web-based system for finding oriCs in unannotated bacterial genomes. BMC Bioinforma. 2008, 9: 79-10.1186/1471-2105-9-79.View ArticleGoogle Scholar
- Lagesen K, Hallin P, Rødland EA, Staerfeldt HH, Rognes T, Ussery DW: RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007, 35: 3100-3108. 10.1093/nar/gkm160.PubMed CentralPubMedView ArticleGoogle Scholar
- Laslett D, Canback B: ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences. Nucleic Acids Res. 2004, 32: 11-16. 10.1093/nar/gkh152.PubMed CentralPubMedView ArticleGoogle Scholar
- Li L, Huang D, Cheung MK, Nong W, Huang Q, Kwan HS: BSRD: a repository for bacterial small regulatory RNA. Nucleic Acids Res. 2013, 41: D233-D238. 10.1093/nar/gks1264.PubMed CentralPubMedView ArticleGoogle Scholar
- Atschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.View ArticleGoogle Scholar
- Finn RD, Tate J, Mistry J, Coggill PC, Sammut JS, Hotz HR, Ceric G, Forslund K, Eddy SR, Sonnhammer EL, Bateman A: The Pfam protein families database. Nucleic Acids Res Database Issue. 2008, 36: D281-D288. 10.1093/nar/gkn226.View ArticleGoogle Scholar
- Zdobnov EM, Apweiler R: InterProScan - an integration platform for the signature-recognition methods in InterPro. Bioinformatics. 2001, 17: 847-848. 10.1093/bioinformatics/17.9.847.PubMedView ArticleGoogle Scholar
- Juncker AS, Willenbrock H, von Heijne G, Nielsen H, Brunak S, Krogh A: Prediction of lipoprotein signal peptides in Gram-negative bacteria. Protein Sci. 2003, 12: 1652-1662. 10.1110/ps.0303703.PubMed CentralPubMedView ArticleGoogle Scholar
- Moriya Y, Itoh M, Okuda S, Yoshizawa A, Kanehisa M: KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007, 35: W182-W185. 10.1093/nar/gkm321.PubMed CentralPubMedView ArticleGoogle Scholar
- Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream MA, Barrell B: Artemis: sequence visualization and annotation. Bioinformatics. 2000, 16: 944-945. 10.1093/bioinformatics/16.10.944.PubMedView ArticleGoogle Scholar
- Read TD, Brunham RC, Shen C, Gill SR, Heidelberg JF, White O, Hickey EK, Peterson J, Utterback T, Berry K, Bass S, Linher K, Weidman J, Khouri H, Craven B, Bowman C, Dodson R, Gwinn M, Nelson W, DeBoy R, Kolonay J, McClarty G, Salzberg SL, Eisen J, Fraser CM: Genome sequences of Chlamydia trachomatis MoPn and Chlamydia pneumoniae AR39. Nucleic Acids Res. 2000, 15: 1397-1406.View ArticleGoogle Scholar
- Carver TJ, Rutherford KM, Berriman M, Rajandream MA, Barrell BG, Parkhill J: ACT: the Artemis comparison tool. Bioinformatics. 2005, 21: 3422-3423. 10.1093/bioinformatics/bti553.PubMedView ArticleGoogle Scholar
- Sullivan MJ, Petty NK, Beatson SA: Easyfig: a genome comparison visualiser. Bioinformatics. 2011, 17: 1009-1010.View ArticleGoogle Scholar
- Chen F, Mackey AJ, Stoeckert CJ, Roos DS: OrthoMCL-DB: querying a comprehensive multi-species collection of ortholog groups. Nucleic Acids Res. 2006, 34: D363-D368. 10.1093/nar/gkj123.PubMed CentralPubMedView ArticleGoogle Scholar
- Grant JR, Stothard P: The CGView server: a comparative genomics tool for circular genomes. Nucleic Acids Res. 2008, 36: W181-W184. 10.1093/nar/gkn179.PubMed CentralPubMedView ArticleGoogle Scholar
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.PubMed CentralPubMedView ArticleGoogle Scholar
- Castresana J: Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000, 17: 540-552. 10.1093/oxfordjournals.molbev.a026334.PubMedView ArticleGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.PubMedView ArticleGoogle Scholar
- Milne I, Lindner D, Bayer M, Husmeier D, McGuire G, Marshall DF, Wright F: TOPALi v2: a rich graphical interface for evolutionary analyses of multiple alignments on HPC clusters and multi-core desktops. Bioinformatics. 2009, 25: 126-127. 10.1093/bioinformatics/btn575.PubMed CentralPubMedView ArticleGoogle Scholar
- Huson DH, Bryant D: Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006, 23: 254-267.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.