Previous Article | Next Article 
Journal of Bacteriology, March 2007, p. 2128-2132, Vol. 189, No. 5
0021-9193/07/$08.00+0 doi:10.1128/JB.00378-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
Genomic Plasticity of the rrn-nqrF Intergenic Segment in the Chlamydiaceae
Zhi Liu,1
Roger Rank,2
Bernhard Kaltenboeck,3
Simone Magnino,4
Deborah Dean,5
Laurel Burall,1
Roger D. Plaut,1
Timothy D. Read,6
Garry Myers,7 and
Patrik M. Bavoil1*
Department of Biomedical Sciences, University of Maryland Dental School, Baltimore, Maryland,1
Department of Microbiology and Immunology, University of Arkansas, Little Rock, Arkansas,2
Department of Pathobiology, Auburn University, Auburn, Alabama,3
Istituto Zooprofilattico Sperimentale della Lombardia e dell'Emilia Romagna, Pavia, Italy,4
Children's Hospital Oakland Research Institute, Oakland, California,5
Biological Defense Research Directorate, Naval Research Center, Rockville, Maryland,6
The Institute for Genome Research, Rockville, Maryland7
Received 16 March 2006/
Accepted 28 November 2006

ABSTRACT
In
Chlamydiaceae, the nucleotide sequence between the 5S rRNA
gene and the gene for subunit F of the Na
+-translocating NADH-quinone
reductase (
nqrF or
dmpP) has varied lengths and gene contents.
We analyzed this site in 45
Chlamydiaceae strains having diverse
geographical and pathological origins and including members
of all nine species.

TEXT
The complete genomes of six
Chlamydiaceae species (
1,
3,
12,
19-
21,
23,
24) have revealed a high level of synteny, with large-scale
rearrangements and niche-specific genes restricted to a "plasticity
zone" located at the replication termination region. Analysis
of the
rrn-nqrF intergenic segments in these sequences reveals
significant variation in length and gene content that is inconsistent
with host range, tissue tropism, DNA-based phylogenies, and
disease spectrum in humans and animals. Therefore, this segment
may represent a new genomic plasticity zone in these bacteria.
We were especially interested in a frameshifted invasin/intimin-like
protein gene (
ilp) at this location in
Chlamydia caviae (
20).
Fragments of
ilp in several
Chlamydia suis strains were also
identified by Dugan et al. (
5). These observations led us to
a comprehensive analysis of the
rrn-nqrF intergenic sequences
across the
Chlamydiaceae. Strains are listed in Table
1. Data
on PCR and sequencing primers are available upon request. Our
findings are provided in a graphic summary (Fig.
1).
Species of the
Chlamydiaceae can be separated into four groups
based on the contents of their
rrn-nqrF intergenic segments.
Group 1 includes isolates from
Chlamydia pneumoniae,
Chlamydia psittaci,
Chlamydia abortus, and
Chlamydia felis, infectious
to primate (respiratory tract and eyes), avian (respiratory
and digestive tracts and eyes), ovine and bovine (digestive
tract and placenta), and feline (respiratory tract and eyes)
species, respectively. The
rrn-nqrF segment is fully conserved
in four
C. pneumoniae strains (
12,
19,
21) and is only 46 to
49% identical to segments in the three other species. The
rrn-nqrF segment of the Japanese
C. felis strain Fe/C-56 (
1) is 99% identical
to that of
C. psittaci strain parakeet but only 63 and 61% identical
to those of
C. psittaci I-10 or 6BC and another
C. felis strain,
FEPN Baker, respectively. The last strain, however, displays
an
rrn-nqrF intergenic sequence closely related to those of
three
C. abortus strains, S26 (
24), B577, and LW508. Nucleotide
substitutions, insertions, and deletions and hypervariable segments
including stretches of C's and T's (Fig.
2) are found in group
1 strains, causing various translational frames. A single open
reading frame (ORF) at the
rrn-nqrF site, identified as CF0129
in the
C. felis strain Fe/C-56 genome (
1), is also found in
C. psittaci strains 6BC, parakeet, and I-10 and in
C. abortus strain B577. However, mutations disrupt this ORF in
C. abortus strains S26 and LW508,
C. felis strain FEPN Baker, and four
C. pneumoniae strains.
Group 2 includes isolates that carry at least partial
ilp sequences
at the
rrn-nqrF site:
C. caviae strain GPIC,
Chlamydia muridarum strain MoPn, and eight
C. suis strains, which infect guinea
pigs (eyes and genital tract), mice (lungs), and pigs (digestive
tract), respectively. Of those, only
C. suis strain S45 (
13)
has intact
ilp; all others have either fragments of the gene
(
C. muridarum) or a gene that is interrupted by mutation (
C. caviae and Tet
r C. suis). The sequence similarity of
C. suis Ilp to
Yersinia pseudotuberculosis invasin (
9) and enteropathogenic
Escherichia coli intimin (
11) is restricted to the transmembrane
and immunoglobulin-like domains (D1). However, the general structure
of the predicted Ilp polypeptide is similar to that of the invasin/intimin
family and includes an extended rodlike structure consisting
of several modular repeat domains (D2 to D7) and a distal domain
(D8). Several conserved beta strands in D2 to D7 are consistent
with the crystal structures of intimin D0 to D2 and invasin
D1 to D3 (
8,
16). Secondary-structure predictions, however,
reveal no conserved salient features in D8 (
http://www.compbio.dundee.ac.uk/
www-jpred/;
domain definitions are per GenBank [accession number DQ076148]).
We reasoned that ilp of the C. caviae reference strain may have been inactivated through mutation during long-term in vitro culture. To address this question, the rrn-nqrF segment was amplified directly from early, archived specimens of C. caviae from the laboratory of Murray (18) (Table 1). Two specimens of the primary isolate that were maintained by passage through live animals, a 1962 specimen corresponding to the first passage and a 1973 specimen, were analyzed. In both cases, two frameshift mutations identical to those obtained from genomic analysis are present.
Group 3 includes only strains of the human pathogen Chlamydia trachomatis. PCR amplicons from six genitotropic strains, including clinical isolates of serovars D, G, H, I, and K and one serovar E reference strain, were compared. Sequences obtained from serovars D, G, and K are identical to the published sequence of serovar D strain UW-3/Cx (23). However, the serovar E sequence contains a 7-nucleotide insertion and 11 additional mutations, 7 of them located within the putative rrn transcription terminator (not shown). The rrn-nrqF intergenic sequence of lymphogranuloma venereum strain L2/434 (12) is identical to that of ocular strain A/HAR-13 (3) but less identical (96 to 99%) to that of urogenital strains (Table 1). The rrn-nrqF intergenic segments of all C. trachomatis strains and serovars contain no recognizable ORF.
Group 4 includes 10 isolates of Chlamydia pecorum obtained from mammalian species on two continents. Amplicons of the segments from these isolates range between 400 and 700 bp (Fig. 3A). Sequence analysis reveals multiple copies (between 11 and 44) of the simple sequence repeat (SSR) 5'-TGCTTTAG-3' (Fig. 3B). BLAST searches reveal no matches, suggesting that this SSR is unique to C. pecorum.
Taken together, our results indicate an unusual degree of genetic
plasticity at the
rrn-nqrF locus in the family
Chlamydiaceae that transcends phylogenetic boundaries. For instance, the
ilp sequence, albeit degenerate or mutated, is common to
C. suis,
C. muridarum, and
C. caviae, i.e., three species representing
three different hosts and two phylogenetic lineages of the
Chlamydiaceae (
7). The apparent genetic plasticity of this site stands in
sharp contrast with that of the neighboring rRNA operon, whose
functional constancy and genetic stability are a basis of current
phylogenetic analysis. It is possible that the plasticity of
the
rrn-nqrF segment stems from its being frequently exposed
during chlamydial development. There are only one or two copies
of
rrn in the
Chlamydia genomes. Elevated transcription activity
of the
rrn operon would require frequent unwinding of supercoiled
DNA, which might promote physical damage (
15), consistent with
the site's being a preferred target for recombination. The presence
of an 89-bp sequence nearly identical to a chlamydial 23S rRNA
domain 1 sequence, flanking
ilp in
C. caviae, may indeed be
a vestige of a previous horizontal
ilp transfer event (Fig.
1). The identification of four genomic islands within
ilp in
seven Tet
r C. suis strains further supports this hypothesis.
Although the Tet
r cassettes are all recombined at the same position
within
ilp (
5), strain-specific differences in the sequences
of
ompA or the Tet
r islands suggest that these recombinations
result from independent integration events (
2,
5). Similarly,
in group 1, although rRNA sequences of
C. felis strains Fe/C-56
and FEPN Baker are almost identical (>99.9%), their
rrn-nqrF intergenic segments are nearly 40% divergent and align with
those of
C. psittaci and
C. abortus, respectively. Our results
suggest that genetic exchange has occurred at this site between
C. felis,
C. psittaci, and
C. abortus. Apart from recombination
events, modifications of the
rrn-nqrF segment include frameshifts
in
ilp of
C. caviae; truncation of
ilp in
C. muridarum; ORF
fusion and gene decay in CF0129-related sequences in
C. pneumoniae,
C. psittaci,
C. abortus, and
C. felis; and unique SSR duplications
and deletions in
C. pecorum.
Sequence variation in pathogens often occurs in genes encoding surface components. Among ORFs of the rrn-nqrF locus, CF0129-related ORFs and ilp encode predicted membrane proteins. The predicted product of the CF0129-related ORF CAB852 of C. abortus S26 includes three transmembrane domains (24) that are conserved in all tested strains of C. abortus, C. felis, and C. psittaci. Slipped mispairing at poly(C) and poly(T) tracts (Fig. 2), however, results in various translational stops and frameshifts. Likewise, although undisrupted ilp isolated from C. suis strain S45 is transcribed in vitro (not shown), ilp sequence interruption by Tetr cassettes in seven C. suis strains (5), frameshift mutations in C. caviae, and a sequence truncation in C. muridarum (19) indicate that ilp is not required for infection. The roles of other members of the invasin/intimin family in pathogenesis also vary. The yersinial invasin confers an advantage to the bacterium early in infection, but loss of invasin function does not prevent infection (10, 17).
Sequence analysis reveals that the 5'-TGCTTTAG-3' octamer occurs four to six times more than a random sequence octamer in all chlamydial genomes and is oriented relative to ori. These properties are similar to previously reported properties of the genus-specific "Chi" recombination hot spot of Escherichia coli (14) and other bacteria (4, 6, 22). However, the function of multiple repeats of the octamer at the rrn-nqrF site of C. pecorum is unknown.
In summary, our analysis has revealed highly polymorphic sequences at the rrn-nqrF site across 45 chlamydial strains from nine species of diverse geographical and pathological origins. Although the underlying mechanism(s) is not clear, the data presented here support the notion that this site is a hot spot for gene recombination in the Chlamydiaceae.

ACKNOWLEDGMENTS
We thank Pier Giorgio Vigo, Iris Labalestra, and Kyle Burall
for technical assistance and an anonymous reviewer for his assistance
in the analysis.
This work was supported by a subcontract to P. M. Bavoil of National Institutes of Health (NIH) grant RO1 AI51472 to G. Myers (The Institute for Genomic Research). Z. Liu was supported by Chongqing Southwest Hospital, China, and NIH grant RO1 AI4310 to P. M. Bavoil. L. Burall was supported by training grant NIDCR T32 DE07309-08. B. Kaltenboeck was supported by NIH grant AI47202.

FOOTNOTES
* Corresponding author. Mailing address: Department of Biomedical Sciences, University of Maryland, Baltimore, Dental School, 650 W. Baltimore St., 7-South, Baltimore, MD 21201. Phone: (410) 706-6789. Fax: (410) 706-0865. E-mail:
pbavoil{at}umaryland.edu.

Published ahead of print on 8 December 2006. 

REFERENCES
1 - Azuma, Y., H. Hirakawa, A. Yamashita, Y. Cai, M. A. Rahman, H. Suzuki, S. Mitaku, H. Toh, S. Goto, T. Murakami, K. Sugi, H. Hayashi, H. Fukushi, M. Hattori, S. Kuhara, and M. Shirai. 2006. Genome sequence of the cat pathogen, Chlamydophila felis. DNA Res. 13:15-23.[Abstract/Free Full Text]
2 - Bush, R. M., and K. D. Everett. 2001. Molecular evolution of the Chlamydiaceae. Int. J. Syst. Evol. Microbiol. 51:203-220.[Abstract]
3 - Carlson, J. H., S. F. Porcella, G. McClarty, and H. D. Caldwell. 2005. Comparative genomic analysis of Chlamydia trachomatis oculotropic and genitotropic strains. Infect. Immun. 73:6407-6418.[Abstract/Free Full Text]
4 - Chedin, F., S. D. Ehrlich, and S. C. Kowalczykowski. 2000. The Bacillus subtilis AddAB helicase/nuclease is regulated by its cognate Chi sequence in vitro. J. Mol. Biol. 298:7-20.[CrossRef][Medline]
5 - Dugan, J., D. D. Rockey, L. Jones, and A. A. Andersen. 2004. Tetracycline resistance in Chlamydia suis mediated by genomic islands inserted into the chlamydial inv-like gene. Antimicrob. Agents Chemother. 48:3989-3995.[Abstract/Free Full Text]
6 - el Karoui, M., D. Ehrlich, and A. Gruss. 1998. Identification of the lactococcal exonuclease/recombinase and its modulation by the putative Chi sequence. Proc. Natl. Acad. Sci. USA 95:626-631.[Abstract/Free Full Text]
7 - Everett, K. D., R. M. Bush, and A. A. Andersen. 1999. Emended description of the order Chlamydiales, proposal of Parachlamydiaceae fam. nov. and Simkaniaceae fam. nov., each containing one monotypic genus, revised taxonomy of the family Chlamydiaceae, including a new genus and five new species, and standards for the identification of organisms. Int. J. Syst. Bacteriol. 49:415-440.[Abstract/Free Full Text]
8 - Hamburger, Z. A., M. S. Brown, R. R. Isberg, and P. J. Bjorkman. 1999. Crystal structure of invasin: a bacterial integrin-binding protein. Science 286:291-295.[Abstract/Free Full Text]
9 - Isberg, R. P., H. L. Voorhis, and S. Falkow. 1987. Identification of invasin: a protein that allows enteric bacteria to penetrate cultured mammalian cells. Cell 50:769-778.[CrossRef][Medline]
10 - Jepson, M. A., and M. A. Clark. 1998. Studying M cells and their role in infection. Trends Microbiol. 6:359-365.[CrossRef][Medline]
11 - Jerse, A. E., J. Yu, B. D. Tall, and J. B. Kaper. 1990. A genetic locus of enteropathogenic Escherichia coli necessary for the production of attaching and effacing lesions on tissue culture cells. Proc. Natl. Acad. Sci. USA 87:7839-7843.[Abstract/Free Full Text]
12 - Kalman, S., W. Mitchell, R. Marathe, C. Lammel, J. Fan, R. W. Hyman, L. Olinger, J. Grimwood, R. W. Davis, and R. S. Stephens. 1999. Comparative genomes of Chlamydia pneumoniae and C. trachomatis. Nat. Genet. 21:385-389.[CrossRef][Medline]
13 - Kaltenboeck, B., and J. Storz. 1992. Biological properties and genetic analysis of the ompA locus in chlamydiae isolated from swine. Am. J. Vet. Res. 53:1482-1487.[Medline]
14 - Kowalczykowski, S. C. 2000. Initiation of genetic recombination and recombination-dependent replication. Trends Biochem. Sci. 25:156-165.[CrossRef][Medline]
15 - Kuraoka, I., M. Endou, Y. Yamaguchi, T. Wada, H. Handa, and K. Tanaka. 2003. Effects of endogenous DNA base lesions on transcription elongation by mammalian RNA polymerase. II. Implications for transcription-coupled DNA repair and transcriptional mutagenesis. J. Biol. Chem. 278:7294-7299.[Abstract/Free Full Text]
16 - Luo, Y., E. A. Frey, R. A. Pfuetzner, A. L. Creagh, D. G. Knoechel, C. A. Haynes, B. B. Finlay, and N. C. Strynadka. 2000. Crystal structure of enteropathogenic Escherichia coli intimin-receptor complex. Nature 405:1073-1077.[CrossRef][Medline]
17 - Marra, A., and R. R. Isberg. 1997. Invasin-dependent and invasin-independent pathways for translocation of Yersinia pseudotuberculosis across the Peyer's patch intestinal epithelium. Infect. Immun. 65:3412-3421.[Abstract]
18 - Murray, E. S. 1964. Guinea pig inclusion conjunctivitis. I. Isolation and identification as a member of the Psittacosis-Lymphogranuloma-Trachoma group. J. Infect. Dis. 114:1-12.[Medline]
19 - Read, T. D., R. C. Brunham, C. Shen, S. R. Gill, J. F. Heidelberg, O. White, E. K. Hickey, J. Peterson, T. Utterback, K. Berry, S. Bass, K. Linher, J. Weidman, H. Khouri, B. Craven, C. Bowman, R. Dodson, M. Gwinn, W. Nelson, R. DeBoy, J. Kolonay, G. McClarty, S. Salzberg, J. Eisen, and C. M. Fraser. 2000. Genome sequences of Chlamydia trachomatis MoPn and Chlamydia pneumoniae AR39. Nucleic Acids Res. 28:1397-1406.[Abstract/Free Full Text]
20 - Read, T. D., G. S. Myers, R. C. Brunham, W. C. Nelson, I. T. Paulsen, J. Heidelberg, E. Holtzapple, H. Khouri, N. B. Federova, H. A. Carty, L. A. Umayam, D. H. Haft, J. Peterson, M. J. Beanan, O. White, S. L. Salzberg, R. C. Hsia, G. McClarty, R. G. Rank, P. M. Bavoil, and C. M. Fraser. 2003. Genome sequence of Chlamydophila caviae (Chlamydia psittaci GPIC): examining the role of niche-specific genes in the evolution of the Chlamydiaceae. Nucleic Acids Res. 31:2134-2147.[Abstract/Free Full Text]
21 - Shirai, M., H. Hirakawa, M. Kimoto, M. Tabuchi, F. Kishi, K. Ouchi, T. Shiba, K. Ishii, M. Hattori, S. Kuhara, and T. Nakazawa. 2000. Comparison of whole genome sequences of Chlamydia pneumoniae J138 from Japan and CWL029 from USA. Nucleic Acids Res. 28:2311-2314.[Abstract/Free Full Text]
22 - Sourice, S., V. Biaudet, M. El Karoui, S. D. Ehrlich, and A. Gruss. 1998. Identification of the Chi site of Haemophilus influenzae as several sequences related to the Escherichia coli Chi site. Mol. Microbiol. 27:1021-1029.[CrossRef][Medline]
23 - Stephens, R. S., S. Kalman, C. Lammel, J. Fan, R. Marathe, L. Aravind, W. Mitchell, L. Olinger, R. L. Tatusov, Q. Zhao, E. V. Koonin, and R. W. Davis. 1998. Genome sequence of an obligate intracellular pathogen of humans: Chlamydia trachomatis. Science 282:754-759.[Abstract/Free Full Text]
24 - Thomson, N. R., C. Yeats, K. Bell, M. T. Holden, S. D. Bentley, M. Livingstone, A. M. Cerdeno-Tarraga, B. Harris, J. Doggett, D. Ormond, K. Mungall, K. Clarke, T. Feltwell, Z. Hance, M. Sanders, M. A. Quail, C. Price, B. G. Barrell, J. Parkhill, and D. Longbottom. 2005. The Chlamydophila abortus genome sequence reveals an array of variable proteins that contribute to interspecies variation. Genome Res. 15:629-640.[Abstract/Free Full Text]
Journal of Bacteriology, March 2007, p. 2128-2132, Vol. 189, No. 5
0021-9193/07/$08.00+0 doi:10.1128/JB.00378-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
This article has been cited by other articles:
-
Burall, L. S., Rodolakis, A., Rekiki, A., Myers, G. S. A., Bavoil, P. M.
(2009). Genomic Analysis of an Attenuated Chlamydia abortus Live Vaccine Strain Reveals Defects in Central Metabolism and Surface Proteins. Infect. Immun.
77: 4161-4167
[Abstract]
[Full Text]