Previous Article | Next Article 
Journal of Clinical Microbiology, June 2001, p. 2126-2133, Vol. 39, No. 6
0095-1137/01/$04.00+0 DOI: 10.1128/JCM.39.6.2126-2133.2001
Copyright © 2001, American Society for Microbiology. All rights reserved.
Phylogeny of Pneumocystis carinii from
18 Primate Species Confirms Host Specificity and Suggests
Coevolution
Christine
Demanche,1
Madeleine
Berthelemy,1
Thierry
Petit,2
Bruno
Polack,1
Ann E.
Wakefield,3
Eduardo
Dei-Cas,4,5 and
Jacques
Guillot1,*
UMR 956 INRA-AFSSA-ENVA Biologie
Moléculaire et Immunologie Parasitaires et Fongiques, École
Nationale Vétérinaire d'Alfort,
Maisons-Alfort,1 Parc Zoologique de La
Palmyre, Le Mathes,2 and
Parasitologie-Mycologie, Faculté de Médecine et
CHRU de Lille,4 and Ecologie du
Parasitisme, Institut Pasteur de Lille,5 Lille,
France, and Molecular Infectious Diseases Group, Institute of
Molecular Medicine, University of Oxford, Oxford, United
Kingdom3
Received 18 January 2001/Returned for modification 8 March
2001/Accepted 8 April 2001
 |
ABSTRACT |
Primates are regularly infected by fungal organisms identified as
Pneumocystis carinii. They constitute a valuable population for the confirmation of P. carinii host specificity. In
this study, the presence of P. carinii was assessed by
direct examination and nested PCR at mitochondrial large subunit
(mtLSU) rRNA and dihydropteroate synthetase (DHPS) genes in 98 lung
tissue samples from captive or wild nonhuman primates. Fifty-nine air
samples corresponding to the environment of different primate species in zoological parks were also examined. Cystic forms of P. carinii were detected in smears from 7 lung tissue samples
corresponding to 5 New World primate species. Amplifications at the
mtLSU rRNA gene were positive for 29 lung tissue samples representing
18 different primate species or subspecies and 2 air samples
corresponding to the environment of two simian colonies. Amplifications
at the DHPS gene were positive for 8 lung tissue samples representing 6 different primate species. Direct sequencing of nested PCR products demonstrated that a specific mtLSU rRNA and DHPS sequence could be
attributed to each primate species or subspecies. No nonhuman primate
harbored the human type of P. carinii (P. carinii f. sp. hominis). Genetic divergence in
primate-derived P. carinii organisms varied in terms of the
phylogenetic divergence existing among the corresponding host species,
suggesting coevolution.
 |
INTRODUCTION |
Pneumocystis carinii
pneumonia (PCP) is still considered one of the most serious fungal
respiratory infections that can occur in immunocompromised patients,
especially human immunodeficiency virus-infected individuals. Molecular
comparisons of various gene sequences (19, 42) clearly
demonstrated that the single name P. carinii corresponds in
fact to a complex group of eukaryotic organisms which should be
assigned to the kingdom Fungi.
Because a continuous in vitro propagation system for P. carinii is still lacking, the basic biology of this fungal group
and subsequent epidemiology of PCP remain poorly understood (8, 11, 42). Airborne transmission appears likely. Hospital
outbreaks have been documented (27), and well-designed
animal studies have provided evidence of direct transmission from host
to host via airborne infective particles (9, 26, 41). It
has been hypothesized that the main source of P. carinii is
patients with PCP. Another source of transmission may be through
maternal-infant exposure (9). Domestic or wild animals are
not considered sources for humans as P. carinii organisms
seem to be characterized by strong host specificity. Several
Pneumocystis species may be distinguished, each of them
residing in a specific mammalian host (29). Isolates of
P. carinii infecting different mammals are divergent at the genetic level (specific karyotypic profiles and gene sequences) (29) and at the phenotypic level (antigenic differences,
ultrastructure, and isoenzymatic polymorphism) (36, 46).
The concept of close host specificity has been further supported by the
failure of most cross-infection experiments (1, 2, 3, 23, 24, 47). However, Sethi (39) described the possibility
for human-derived P. carinii to develop in SCID mice, and
one recent study suggested transient colonization of owl monkeys
(Aotus nancymai) by human-derived P. carinii
(6).
In this report, the genetic diversity of P. carinii from
primates is examined by analyzing mitochondrial large subunit (mtLSU) rRNA and dihydropteroate synthetase (DHPS) gene sequences. Each of the
18 nonhuman primate species or subspecies which were proved to harbor
P. carinii had its own type of organism with specific mtLSU
rRNA or DHPS sequences. No P. carinii f. sp.
hominis was found in the nonhuman primate lung tissue or air
samples examined in the present work.
 |
MATERIALS AND METHODS |
Samples for analysis.
Postmortem lung tissues from nonhuman
primates were obtained in four French zoological parks (La Palmyre,
Jardin des Plantes de Paris, Parc Zoologique de Vincennes, and Parc
Zoologique de Mulhouse) and from the Primate Research Center of
Strasbourg. Additional lung tissues from wild monkeys were obtained
from the ONC (Office National de la Chasse) of French Guyana. The lung tissues were frozen after necropsic examination and stored at
20°C
prior to direct examination and DNA extraction. Air samples corresponding to the environment of different primate species in
zoological parks were obtained by using the CAP device (Arelco, Fontenay sous Bois, France), as described by Guillot et al.
(25). This device sampled airborne particles with a flow
rate of 10 liters/min. Particles were impacted on the surface of a
rotative cup. Air sampling was performed in the four zoological parks
and in the primatology research center of Strasbourg for 24 to 72 h in front of or inside each of the cages where simian colonies were maintained.
Direct examination of lung tissue samples.
Impression smears
of cross sections of lungs were stained with toluidine blue O for the
detection of cystic forms of P. carinii (10). A
small part of the lung tissue (from 300 to 500 mg) was then finely
minced, homogenized with crushing, and subjected to sequential membrane
filtration (9). The final filtrates were used for
direct examination and DNA extraction. Cystic forms of P. carinii were detected and counted in 2.5-µl air-dried smears of
the final filtrate, stained with toluidine blue O.
DNA extraction from lung tissue and air samples.
A volume of
100 µl of the final filtrate of lung extract was first frozen at
20°C and digested with proteinase K (Boehringer Mannheim) at a
final concentration of 0.34 mg/ml. A phenol-chloroform extraction was
then performed with a final precipitation in ethanol. For air samples,
the rotative cup (from the CAP apparatus) was washed with 600 µl of
extraction buffer (10 mM Tris, 0.5% sodium dodecyl sulfate, 25 mM
EDTA, 0.1 M NaCl). DNA was prepared by proteinase K digestion
(Boehringer Mannheim) at a final concentration of 0.28 mg/ml, followed
by phenol-chloroform extraction.
Primers and PCR amplifications.
The presence of P. carinii DNA in lung tissue and air samples was assessed by nested
PCR at the mtLSU rRNA gene and at the DHPS locus. For mtLSU rRNA gene
amplification, the primer sets pAZ102-H-pAZ102-E (5'-GAT GGG TGT
TTC CAA GCC CA-3' and 5'-GTG TAC GTT GCA AAG TAC TC-3')
and pAZ102-X/R1-pAZ102-Y/R1 (5'-GGG AAT TCG TGA AAT ACA AAT
CGG ACT AGG-3' and 5'-GGG AAT TCT CAC TTA ATA TTA ATT GGG
GAG C-3') were used (45). The thermocycling conditions for the first PCR round were as follows: each cycle consisted of denaturation for 30 s at 94°C, annealing for 1 min at
50°C, and extension for 2 min at 72°C for 30 cycles. The second round of PCR was performed with 5% (vol/vol) of the first-round mix.
The thermocycling conditions for the second PCR round were as follows:
each cycle consisted of denaturation for 30 s at 94°C, annealing for
1 min at 55°C, and extension for 2 min at 72°C for 30 cycles.
For the first round of PCR at the DHPS gene, the primer set
AHUM-BHUM (5'-GCG CCT ACA CAT ATT ATG GCC
ATT TTA AAT C-3' and 5'-CAT AAA CAT CAT GAA CCC G-3')
was used (31). The thermocycling conditions were as
follows: 10 first cycles consisted of denaturation for 30 s at
94°C, annealing for 1 min at 52°C, and extension for 1 min at
72°C; 25 additional cycles consisted of denaturation for 30 s at
94°C, annealing for 1 min at 42°C, and extension for 1 min at
72°C. The second round of PCR was performed with 5% (vol/vol) of the
first-round mix and the primer set CPRIM-DPRIM
(5'-CCC CCA CTT ATA TCA-3' and 5'-GGG GGT GTT CAT
TCA-3'). The thermocycling conditions for the second PCR round
were as follows: each cycle consisted of denaturation for 30 s at
94°C, annealing for 1 min at 50°C, and extension for 1 min at
72°C for 30 cycles. Negative controls were included in each
experiment, in both DNA extraction and PCR amplification, to monitor
for possible contamination.
Amplification products were purified in a 2% agarose gel
(Tris-borate-EDTA buffer) and extracted with the Geneclean II kit
(Ozyme, Montigny-le-Bretonneux, France) when nonspecific bands
were
detected. Amplification products were directly sequenced
from both ends
by using sets of internal primers on an automated
DNA sequencer
(GenomeExpress, Montreuil, France). The sequences
have been submitted
to
GenBank.
Molecular phylogenetic analysis.
Nucleotide (mtLSU and DHPS)
and derived amino acid (DHPS) sequences were initially aligned with the
computer program CLUSTAL X (version 1.63b, December 1997)
(44) and then by visual optimization. Only regions without
ambiguity were included in the phylogenetic analysis (alignments are
available upon request). The aligned sequences were converted to
distance matrix (% of differences). For mtLSU rRNA and DHPS sequences,
the distance trees (phenetic trees) were generated using the
neighbor-joining method (38). Both the branch-and-bound
and heuristic search options in the Phylogenetic Analysis Using
Parsimony Program (PAUP 4.0; distributed by the Illinois Natural
History Survey, Champaign, Ill.) were used for comparison of sequence
alignments and generation of parsimonious trees. The strength of the
internal branches from the resulting trees was statistically tested by
bootstrap analysis (20) from 500 bootstrap replications.
For mtLSU rRNA sequence analyses, P. carinii strains from
rat (P. carinii f. sp. carinii, GenBank accession
number U42914) and from mouse (P. carinii f. sp. muris, GenBank accession number U20169) were chosen as
outgroups. For DHPS sequence analyses, P. carinii from mouse
(P. carinii f. sp. muris, GenBank accession
number U66283) was chosen as the outgroup.
 |
RESULTS |
Direct examination and PCR amplification.
A total of
98 lung tissue samples from 33 different primate species were collected
(Table 1). Three zoological groups were represented, with 72 lung tissue samples from 15 New World primate species (Platyrrhini), 21 lung tissue samples from 14 Old World primate
species (Catarrhini), and 5 lung tissue samples from 4 lemurs
(Strepsirhini Lemuridae). Ninety-four lung tissue samples were obtained
from French zoological parks and four lung tissue samples were from
wild primates in French Guyana (two pale-headed sakis, one red-handed
tamarin, and one squirrel monkey). The most frequent causes of death
noted at necropsic examination included infectious diseases and trauma.
Some monkeys presented pulmonary lesions but no case of pneumocystosis
could be diagnosed. Cystic forms of P. carinii were detected
in smears from seven lung tissue samples corresponding to five New
World primate species (one squirrel monkey, two Geoffroy's marmosets,
one red-handed tamarin, one Goeldi's monkey, and two Weddell's
tamarins). The organisms (3.5 to 5.9 µm in diameter) were either
isolated or disposed in small clusters. Amplification at mtLSU rRNA and
DHPS genes was positive with the first round of PCR for all the lung
tissue samples in which cystic forms were detected. Additional positive
amplification was obtained with the second round of PCR for 22 lung
tissue samples (representing 14 different primate species) at the mtLSU
rRNA gene and for 1 sample (from a swamp guenon) at the DHPS gene
(Table 1).
A total of 59 air samples were collected from the environment of 31 nonhuman primate species. Amplification at the mtLSU rRNA
gene using
nested PCR was positive from only two samples: one
from the environment
of a colony comprising 10 common marmosets
(
Callithrix
jacchus) in La Palmyre zoological park and another
one from the
environment of a colony comprising 5 black lemurs
(
Maki
macaco) in Vincennes zoological park (Table
1). No amplification
at the DHPS gene was obtained from air
samples.
PCR products ranged from 313 to 324 bp (first round of PCR at the mtLSU
rRNA gene), from 234 to 282 bp (nested PCR at the
mtLSU rRNA gene), and
from 702 to 720 bp (first round of PCR at
the DHPS gene). Amplification
was positive with the nested PCR
at the DHPS gene for one lung tissue
sample (from a swamp guenon),
and the PCR product comprised 664
bp.
Mitochondrial LSU rRNA sequence comparison and genetic
groupings.
Direct sequencing of PCR products demonstrated that a
specific mtLSU rRNA sequence could be attributed to each primate
species or subspecies. In the case of the common marmoset (C. jacchus), the PCR amplification for mtLSU rRNA was positive for
several lung tissue samples and for one air sample. The corresponding sequences were shown to be identical. In the case of the red-handed tamarin (Saguinus midas midas), positive PCR amplifications
were obtained from both captive and wild animals and the corresponding mtLSU rRNA sequences were identical. Multiple sequence types were not
detected in a single host species except for the rhesus monkey (Macaca mulatta). On the other hand, the sequence of the
human-derived P. carinii was never obtained from lung tissue
or air samples examined in the present study.
Eighteen new
P. carinii mtLSU rRNA sequences were obtained
from nine New World monkey species, six Old World monkey species,
and
two lemurs. The sequences were identified by the GenBank accession
numbers
AF362454 (
P. carinii from
C. jacchus),
AF362456 (
P. carinii from
Callithrix geoffroyi),
AF362461 (
P. carinii from
Callimico goeldii),
AF362462 (
P. carinii from
Saguinus fuscicolis),
AF362455 (
P. carinii from
S. midas midas),
AF362453 (
P. carinii from
Saguinus oedipus
oedipus),
AF362465 (
P. carinii from
Saguinus
imperator),
AF362458 (
P. carinii from
Saimiri sciureus),
AF362470 (
P. carinii from
Pithecia
pithecia),
AF362464 (
P. carinii from
Allenopithecus nigroviridis),
AF362457 (
P. carinii from
Cercopithecus hamlyni),
AF362460 (
P. carinii from
Cercopithecus nictitans),
AF362469
(
P. carinii from
Macaca fascicularis),
AF362467
(
P. carinii from an
M. mulatta Chinese
subspecies),
AF362468 (
P. carinii from an
M. mulatta Indian
subspecies),
AF362466 (
P. carinii from
Macaca nemestrina),
AF362459 (
P. carinii from
Hapalemur griseus), and
AF362463 (
P. carinii from
M. macaco). In the case of the rhesus monkey,
two sequence
types were observed according to the subspecies (Chinese
or Indian
M. mulatta). In the region selected for the phylogenetic
analysis, only one base substitution distinguished the two sequence
types. A matrix analysis of sequence divergence between
P. carinii mtLSU rRNA sequences is depicted in Table
2. Comparison of the
mtLSU rRNA aligned
sequences was carried out on 287 positions,
including gaps, for a total
of 21 taxa: 18 original sequence types
and 3 sequences already
published,
P. carinii f. sp.
hominis (GenBank
S42926),
P. carinii f. sp.
muris (GenBank
U20169), and
P. carinii f. sp.
carinii (GenBank
U42914).
P. carinii organisms
from Old World primates
(including human) were characterized by
an insertion at position 36. For members of the genera
Allenopithecus and
Cercopithecus the insertion was very long (63 bp), whereas
for macaques and human-derived
P. carinii the insertion was
shorter
(from 13 to 48 bp). Lemur-derived
P. carinii
organisms were characterized
by a distinct insertion at position 132 (15 bp for the grey gentle
lemur and 22 bp for the black lemur).
View this table:
[in this window]
[in a new window]
|
TABLE 2.
Matrix of mtLSU rRNA sequence divergence for P. carinii from 17 primate species (including P. carinii
f. sp. hominis) and rodents
|
|
No insertion was detected in mtLSU rRNA sequences of New World
primate-derived
P. carinii. The two primate species of the
genus
Callithrix (Geoffroy's marmoset and common marmoset)
harbored
distinct
P. carinii with eight nucleotide
differences in the mtLSU
rRNA sequences but these differences were
located in a region
which was not included in the phylogenetic
analysis. The same
situation was observed for
P. carinii
from two primate species
belonging to the genus
Saguinus
(Emperor tamarin and red-handed
tamarin), with 10 nucleotide
differences in the amplified mtLSU
rRNA sequences. Within the group
composed of primate-derived
P. carinii, the maximum
divergence (28.2%) was observed between the
black lemur-derived and
the squirrel monkey-derived
P. carinii.
The extent of
divergence may not be greater between a primate-derived
and a
nonprimate-derived
P. carinii than between primate-derived
organisms themselves (Table
2).
The phylogenetic analysis inferred from mtLSU rRNA sequence comparison
demonstrated that genetic groupings within primate-derived
P. carinii were in accordance with those of the corresponding
hosts.
Phylogenetic analyses, neighbor joining, heuristic search,
and branch
and bound produced topologically similar trees, although
statistical
support among the branches varied (Fig.
1). All three
analyses indicated that
primate-derived
P. carinii formed a monophyletic
group with
three distinct clades. The two lemur-derived
P. carinii organisms were associated and occupied a basal position. All
P. carinii organisms from Old World nonhuman primates were included
in the same clade supported by a 95% bootstrap value. Within this
group,
P. carinii from macaques (
Macaca spp.) and
from guenons
(
Allenopithecus and
Cercopithecus
spp.) were clearly distinguished
with branches supported by 98 and
100% bootstrap values, respectively.
A branch associating
P. carinii f. sp.
hominis and Old World primate-derived
P. carinii was systematically produced (heuristic search and
branch
and bound) but with a low bootstrap value (58%). The last major
clade included New World monkey-derived
P. carinii with a
high
bootstrap value (80%). Within this group,
P. carinii
from tamarins
(genera
Saguinus and
Callimico)
clustered together.
P. carinii organisms from marmosets
(genus
Callithrix) were also associated.
P. carinii from pale-headed sakis (
P. pithecia) occupied a
paraphyletic
position in the New World monkeys group. The exact
position of
P. carinii from the New World species
S. sciureus (squirrel monkey)
could not be elucidated.

View larger version (33K):
[in this window]
[in a new window]
|
FIG. 1.
Phylogenetic relationships of P. carinii
organisms from primate species inferred from mtLSU rRNA sequences. The
phylogram presented resulted from bootstrapped data sets
(20) obtained by using parsimony analysis (heuristic
search option in PAUP 4.0). This tree was identical to the consensus of
18 most-parsimonious trees generated from the branch and bound
algorithm in PAUP 4.0. The percentages above the branches are the
frequencies with which a given branch appeared in 500 bootstrap
replications. Bootstrap values below 50% are not displayed. The
position of the squirrel monkey (S. sciureus)-derived
P. carinii was not firmly established (dotted line). Branch
lengths correspond to the total nucleotide changes assigned to each
branch by PAUP 4.0. P. carinii from rat (P. carinii f. sp. carinii, GenBank no. U42914) and from
mouse (P. carinii f. sp. muris, GenBank no.
U20169) were chosen as outgroups.
|
|
DHPS sequence comparison and genetic groupings.
Direct
sequencing of PCR products demonstrated that a specific DHPS sequence
could be attributed to each primate species. Six different DHPS
sequences were described and corresponded to P. carinii
derived from five New World monkey species and one Old World monkey
species. The sequences were identified by GenBank accession numbers
AF362758 (P. carinii from C. geoffroyi), AF362760
(P. carinii from C. goeldii), AF362761 (P. carinii from S. fuscicolis), AF362762 (P. carinii from S. midas midas), AF362759 (P. carinii from S. sciureus) and AF362757 (P. carinii from A. nigroviridis). The DHPS aligned
sequence comprised 217 positions, including gaps for a total of nine
taxa: six original sequence types and three sequences already
published, P. carinii f. sp. hominis (GenBank
U66282), P. carinii f. sp. muris (GenBank U66283), and the sequence from the owl monkey (A. nancymai)-derived P. carinii (6). The last
sequence was very short. The lowest divergence was observed between
Weddell's tamarin and Goeldi's monkey-derived P. carinii
(1.4%). The higher level of divergence was found between the squirrel
monkey-derived P. carinii and P. carinii f. sp.
muris (16.6%). Point mutations at positions 165 (A/G) and
171 (C/T) related to sulfa drug resistance in P. carinii f.
sp. hominis were not detected in any of the DHPS sequences obtained in the present study (31).
The deduced partial amino acid sequences of the DHPS gene from
primate-derived
P. carinii contain 72 residues. Very low
levels
of divergence were observed between
P. carinii
organisms from
New World monkeys, Weddell's tamarin, Goeldi's monkey,
owl monkey,
and Geoffroy's marmoset (1.4%). Higher levels of
divergence were
found between the squirrel monkey-derived
P. carinii and Old World
primate-derived
P. carinii
(9.9%), between the Geoffroy's marmoset-derived
P. carinii
and
P. carinii f. sp.
hominis (9.9%), and
between
P. carinii f. sp.
muris and the squirrel
monkey and Geoffroy's marmoset-derived
P. carinii (16.9%).
The phylogenetic analyses inferred from DHPS nucleotide sequences
yielded identical topologies and confirmed the results obtained
with
mtLSU rRNA sequence analyses (Fig.
2).
The monophyletic group
composed of primate-derived
P. carinii organisms was clearly divided
in two clades. A first clade
comprised the two
P. carinii organisms
from Old World
primates (human and swamp guenon-derived
P. carinii).
A
second clade included all
P. carinii organisms from New
World
monkeys. Within this group,
P. carinii organisms from
tamarins
(genus
Saguinus) clustered together. As already
observed with
the analyses of mtLSU rRNA sequences, the branch
corresponding
to
P. carinii organisms from New World monkeys
was supported by
a higher bootstrap value (98%) than that
corresponding to
P. carinii organisms from Old World
primates (52%).
P. carinii from the squirrel
monkey had a
divergent DHPS sequence. However, it was included
in the group of New
World monkeys and was systematically associated
with the owl
monkey-derived
P. carinii. When the short sequence
of
P. carinii from owl monkey (
6) was excluded,
the phylogenetic
analyses could be performed on much longer sequences
(629 bp)
but the tree topologies remained unchanged, with the squirrel
monkey-derived
P. carinii included in the New World monkeys
P. carinii group.

View larger version (19K):
[in this window]
[in a new window]
|
FIG. 2.
Phylogenetic relationships of P. carinii
organisms from primate species inferred from DHPS sequences. The
presented phylogram resulted from bootstrapped data sets
(20) obtained by using parsimony analysis (heuristic
search option in PAUP 4.0). This tree was identical to the consensus of
three most-parsimonious trees generated from the branch and bound
algorithm in PAUP 4.0. The percentages above the branches are the
frequencies with which a given branch appeared in 500 bootstrap
replications. Bootstrap values below 50% are not displayed. Branch
lengths correspond to the total nucleotide changes assigned to each
branch by PAUP 4.0. P. carinii from mouse (P. carinii f. sp. muris, GenBank no. U66283) was chosen as
the outgroup.
|
|
The phylogenetic analyses inferred from DHPS amino acid sequences were
partially in accordance with those inferred from nucleotide
sequences.
All primate-derived
P. carinii organisms were associated
and
five New World monkey (including the squirrel monkey)-derived
P. carinii organisms clustered together. Surprisingly, the red-handed
tamarin-derived
P. carinii clustered with the swamp
guenon-derived
P. carinii, and
P. carinii f. sp.
hominis was paraphyletic to
the nonhuman primate-derived
P. carinii.
 |
DISCUSSION |
P. carinii was originally considered to be a single
organism responsible for pulmonary colonization or infection in a very wide range of mammalian hosts. Frenkel was the first author to suspect
that the situation might be more complex and to suggest a distinction
between human and rodent-derived P. carinii
(21). We now know that the single name P. carinii consists in fact of a heterogeneous group of fungal
organisms (11, 42). Molecular techniques, including
karyotyping, DNA sequence analysis, or single-strand conformation
polymorphism, have played a significant role in demonstrating the
degree of heterogeneity between P. carinii organisms. To
date, the mtLSU rRNA sequence has been reported from P. carinii organisms infecting nine different mammalian species:
human, rat (two divergent sequence types), mouse, rabbit, pig, horse,
shrew, ferret (five divergent sequence types), and rhesus macaque
(17, 23, 37, 46). The genetic diversity of P. carinii was examined at other loci such as the mitochondrial small
subunit of rRNA (28), the nuclear ribosomal RNA operon
(32, 46),
- and
-tubulin (18, 43),
arom (5), thymidilate synthetase (30,
35), DHPS (31), and more recently
manganese-dependent superoxide dismutase (16). However,
only a small variety of mammalian hosts with P. carinii were
examined in former analyses. The present study is the first one to
include so many species and to propose a phylogenetic analysis for
primate-derived P. carinii. We have used the sequences from
two very different types of genes; one is a mitochondrial gene encoding
ribosomal RNA and the other is a nuclear gene encoding a protein coding
sequence. Both nucleotide sequence analyses provided the same
phylogenetic relationships. A few incongruities were observed with the
analyses of amino acid DHPS sequences but the number of residues was low.
In 1994, the Pneumocystis workshop (29)
proposed a tripartite nomenclature for P. carinii. In
accordance with this provisional system we suggest the creation of five
new variant names for P. carinii organisms from nonhuman
primates (Table 3). The
Pneumocystis workshop also proposed a genetic ranking system
based on sequence comparison at five loci: mtLSU rRNA, thymidylate
synthetase,
-tubulin, arom, and internal transcribed
spacers in the nuclear ribosomal RNA genes. Class III corresponded to
the highest level of divergence (from 15 to 50%) observed between
P. carinii organisms from different mammalian species. Class
II divergence ranged from 4 to 7% at the selected loci (except at the
much more variable internal transcribed spacer regions) and could be
observed between different P. carinii species in the same
mammalian host (for example, P. carinii f. sp.
carinii and P. carinii f. sp. ratti in
rat). The lowest level of divergence (lower than 1%), termed class I,
induced a few base substitutions in the sequences of isolates within
the same host species.
The results of the present work suggest that in primate-derived
P. carinii populations, genetic differences vary in terms of
the phylogenetic divergence existing among the host species. In other
words, the sequence divergence between two P. carinii subpopulations appears to be correlated with the phylogenetic relationship between the corresponding hosts. When host species belong
to different mammalian orders (Insectivora, Lagomorpha, Rodentia,
Primates, Carnivora), the class III divergence level is usually
observed. Within the order Primates, the divergence level was shown to
vary from less than 1% to more than 25%. Thus, a regular, progressive
increase in mtLSU rRNA sequence divergence was observed between
P. carinii organisms from the crab-eating macaque (M. fascicularis) and P. carinii from progressively distant species: 0.4% (P. carinii from M. nemestrina),
3.5% (P. carinii from M. mulatta), 9.5 to 10.3%
(P. carinii from Old World primates belonging to a genus
other than Macaca), 16.8 to 25.4% (P. carinii from New World monkeys), and 23.0 to 25.3% (P. carinii from
lemurs and rodents).
Class I divergence was revealed in P. carinii organisms from
rhesus monkeys. Two closely related mtLSU rRNA sequences were each
associated with the Chinese and Indian M. mulatta
subspecies. Similar sequences were reported previously by Furuta et al.
(23) and Durand-Joly et al. (17),
respectively, who noticed an insertion at the beginning of the
sequence. In the present work, this insertion appeared as typical to
Old World primate-derived P. carinii. Class II divergence
was not found among primate-derived P. carinii organisms in
this study. This could be due to the fact that in this first approach
direct sequencing was used.
Interestingly, all the variants of P. carinii found in lung
tissue samples from 16 primate species housed in French zoological parks for many years represented new P. carinii mtLSU rRNA
sequences except the two sequences isolated in M. mulatta
previously reported. Variants of P. carinii detected in air
samples from the environment of two primate species (common marmoset
and black lemur) also represented new sequences. These findings
suggested a persistent circulation of P. carinii in simian
ecosystems, including in captivity. P. carinii f. sp.
hominis, the unique species identified in Homo sapiens
sapiens (29), the sole primate species dwelling in
France at present, was not found in any of the nonhuman primate species examined. The capacity of human-derived P. carinii to infect
animals was not the focus of this study. In the past, attempts to
infect animals with parasites from humans were often performed using recipient hosts which were latently infected with P. carinii
(3, 6, 39), and molecular tools were not available in the
oldest studies (48). On the other hand, recent
cross-infection experiments, developed by using P. carinii-free recipient hosts and molecular methods to identify the
parasite isolates (15), have revealed a strong host
specificity in P. carinii strains from rat, mouse, rabbit,
ferret, or macaque (1, 2, 3, 23). However, narrow parasite
specificity might not be a universal feature of the genus
Pneumocystis. For instance, Dermatophyte genera like Trichophyton or Microsporum contain stenoxenic,
oligoxenic, or euryxenic species, i.e., with different degrees of
parasite host specificity (13). Furthermore, even if
human-derived P. carinii natural populations were found to
exhibit restricted host specificity, isolates of the parasite might
develop in unusual hosts under experimental conditions. For instance,
Plasmodium falciparum, a human Apicomplexa protist with
strong host specificity, can be propagated in New World primates under
experimental conditions (40).
Since culture methods for the routine isolation of P. carinii f. sp. hominis from patients do not yet exist,
the alternative of isolation of human-derived parasites using animal
models has to be explored. Although in the present work human-derived
P. carinii isolates showed relatively high divergence from
P. carinii isolates from lemurs or from some neotropical
primates (Table 2), on the whole P. carinii f. sp.
hominis was found to be phylogenetically closer to most
primate-derived P. carinii isolates, especially those from
Old World primates, than to rodent-derived P. carinii (Table
2; Fig. 1 and 2). It is noteworthy that current drug testing and
therapeutic protocols against PCP are performed using P. carinii strains of rodent origin (4, 7). For these
reasons, the host specificity of P. carinii f. sp.
hominis should be further explored using nonhuman primates
as experimental hosts.
Finally, the results from this study strengthen the view of P. carinii organisms as a new biological group with numerous species (15, 22) widely present in ecosystems (14).
Although for almost one century P. carinii was considered to
be a unique taxonomic entity, recent research has clearly shown that
P. carinii is in fact a group of heterogeneous populations
(34, 42, 47) genetically isolated from each other
(36) that have undergone a prolonged process of genetic
and functional adaptation to each mammal species. This process is
probably speciation, and the status of P. carinii natural
populations conforms therefore to the biological definition of species
(12, 33).
 |
ACKNOWLEDGMENTS |
We acknowledge P. Moisson (Parc Zoologique de Mulhouse), J. Rigoulet (Jardin des Plantes de Paris), A. Lécu and F. Ollivet (Parc Zoologique de Vincennes), A. Gessain (Institut Pasteur de Paris),
H. Contamin and M. Kazanji (Institut Pasteur de Cayenne), E. André and N. Herrenschmidt (Centre de Primatologie de
Strasbourg), and C. Gottini (Muséum National d'Histoire
Naturelle) for providing lung tissue samples from captive primates and
E. Hansen (Office Nationale de la Chasse de Guyane Française) for
providing lung tissue samples from wild primates. We also thank J. P. Hugot (Muséum National d'Histoire Naturelle) for his
assistance in sequence analyses and R. Chermette (Ecole Nationale
Vétérinaire d'Alfort) for helpful discussions.
This study was developed in the framework of the PRFMMIP project
(Programme de Recherche Fondamentale en Microbiologie et Maladies
Infectieuses et Parasitaires, French Ministry of Education, Research
and Technology) and of the European Key Action "Eurocarinii" (5th
PCRD, QLK2-CT2000-01369).
 |
FOOTNOTES |
*
Corresponding author. Mailing address: Service de
Parasitologie-Mycologie, École Nationale Vétérinaire
d'Alfort, 7, Avenue du Général de Gaulle, 94704 Maisons-Alfort, France. Phone: 331 43 96 71 57. Fax: 331 43 75 35 07. E-mail: j.guillot{at}vet-alfort.fr.
 |
REFERENCES |
| 1.
|
Aliouat, E. M.,
E. Mazars,
E. Dei-Cas,
J. Y. Cesbron, and D. Camus.
1993.
Intranasal inoculation of mouse, rat or rabbit-derived Pneumocystis in SCID mice.
J. Protozool. Res.
3:94-98.
|
| 2.
|
Aliouat, E. M.,
E. Mazars,
E. Dei-Cas,
P. Delcourt,
P. Billault, and D. Camus.
1994.
Pneumocystis cross infection experiments using SCID mice and nude rats as recipient host showed strong host-species specificity.
J. Eukaryot. Microbiol.
41:S71[CrossRef].
|
| 3.
|
Atzori, C.,
F. Agostoni,
E. Angeli,
A. Mainini,
V. Micheli, and A. Cargnel.
1999.
Pneumocystis carinii host specificity: attempt of cross infection with human derived strains in rats.
J. Eukaryot. Microbiol.
46:S112.
|
| 4.
|
Aviles, P.,
E. M. Aliouat,
A. Martinez,
E. Dei-Cas,
E. Herreros,
L. Dujardin, and D. Gargallo-Viola.
2000.
In vitro pharmacodynamic parameters of sordarin derivatives in comparison with those of marketed compounds against Pneumocystis carinii isolated from rats.
Antimicrob. Agents Chemother.
44:1284-1290[Abstract/Free Full Text].
|
| 5.
|
Banerji, S.,
A. E. Wakefield,
A. G. Alien,
D. J. Maskell,
S. E. Peters, and J. M. Hopkin.
1993.
The cloning and characterization of the arom gene of Pneumocystis carinii.
J. Gen. Microbiol.
139:2901-2914[Abstract/Free Full Text].
|
| 6.
|
Beard, B. C.,
V. M. Jennings,
W. G. Teague,
J. L. Carter,
J. Mabry,
H. Moura,
G. S. Visvesvara,
W. E. Collins, and T. R. Navin.
1999.
Experimental inoculation of immunosuppressed owl monkeys with Pneumocystis carinii f. sp. hominis.
J. Eukaryot. Microbiol.
46:S113-S115.
|
| 7.
|
Brun Pascaud, M.,
E. Herreros,
E. M. Aliouat, and E. Dei-Cas.
1998.
Evaluation of drug efficacy by using animal models or in vitro systems.
FEMS Immunol. Med. Microbiol.
22:173-179[CrossRef][Medline].
|
| 8.
|
Cailliez, J. C.,
N. Séguy,
C. M. Denis,
E. M. Aliouat,
E. Mazars,
L. Polonelli,
D. Camus, and E. Dei-Cas.
1996.
Pneumocystis carinii: an atypical fungal micro-organism.
J. Med. Vet. Mycol.
34:227-239[Medline].
|
| 9.
|
Ceré, N.,
B. Polack,
N. K. Chanteloup, and P. Coudert.
1997.
Natural transmission of Pneumocystis carinii in nonimmunosuppressed animals: early contagiousness of experimentally infected rabbits.
J. Clin. Microbiol.
35:2670-2672[Abstract].
|
| 10.
|
Chalvardjian, A. M., and L. A. Grave.
1973.
A new procedure for the identification of Pneumocystis carinii cysts in tissue sections and smears.
J. Clin. Pathol.
16:383-384.
|
| 11.
|
Cushion, M.
1998.
Chapter 34: Pneumocystis carinii, p. 645-683.
In
L. Ajello, and R. J. Hay (ed.), Topley and Wilson's microbiology and microbial infections, 9th, vol. 4. Mycology. Arnold, London, England.
|
| 12.
|
Dei-Cas, E.
2000.
Pneumocystis infections: the iceberg?
Med. Mycol.
38(Suppl. 1):23-32.
|
| 13.
|
Dei-Cas, E., and A. Vernes.
1986.
Parasitic adaptation of pathogenic fungi to the mammalian hosts.
Crit. Rev. Microbiol.
13:173-218[Medline].
|
| 14.
|
Dei-Cas, E., and J. C. Cailliez.
1998.
Editorial.
FEMS Immunol. Med. Microbiol.
22:1-4[CrossRef][Medline].
|
| 15.
|
Dei-Cas, E.,
E. Mazars,
E. M. Aliouat,
G. Nevez,
J. C. Cailliez, and D. Camus.
1998.
The host-specificity of Pneumocystis carinii.
J. Mycol. Méd.
8:1-6.
|
| 16.
|
Denis, C. M.,
E. Mazars,
K. Guyot,
C. Odberg-Ferragut,
E. Viscogliosi,
E. Dei-Cas, and A. E. Wakefield.
2000.
Genetic divergence at the SODA locus of six different formae speciales of Pneumocystis carinii.
Med. Mycol.
38:289-300[Medline].
|
| 17.
|
Durand-Joly, I.,
A. E. Wakefield,
R. J. Palmer,
C. M. Denis,
C. Creusy,
L. Fleurisse,
I. Ricard,
J. P. Gut, and E. Dei-Cas.
2000.
Ultrastructural and molecular characterization of Pneumocystis carinii isolated from a rhesus monkey (Macaca mulatta).
Med. Mycol.
38:61-72[Medline].
|
| 18.
|
Edlind, T. D.,
M. S. Bartlett,
G. A. Weinberg,
G. N. Prah, and J. W. Smith.
1992.
The beta-tubulin gene from rat and human isolates of Pneumocystis carinii.
Mol. Microbiol.
6:3365-3373[CrossRef][Medline].
|
| 19.
|
Edman, J. C.,
J. C. Kovacs,
H. Masur,
D. V. Santi,
H. J. Helwood, and M. L. Sogin.
1988.
Ribosomal RNA sequence shows Pneumocystis carinii to be a member of the Fungi.
Nature
334:519-522[CrossRef][Medline].
|
| 20.
|
Felsenstein, J.
1985.
Confidence limits on phylogenies: an approach using the bootstrap.
Evolution
39:783-791[CrossRef].
|
| 21.
|
Frenkel, J. K.
1976.
Pneumocystis jiroveci n. sp. from man: morphology, physiology and immunology in relation to pathology.
Natl. Cancer Inst. Monogr.
43:13-30.
|
| 22.
|
Frenkel, J. K.
1999.
Pneumocystis pneumonia, an immunodeficiency-dependent disease: a critical historical overview.
J. Eukaryot. Microbiol.
46:S89-S92.
|
| 23.
|
Furuta, T.,
M. Fujita,
R. Mukai,
I. Sakakibara,
T. Sata,
K. Miki,
M. Hayami,
S. Kojima, and Y. Yoshikawa.
1993.
Severe pulmonary pneumocystosis in simian acquired immunodeficiency syndrome induced by simian immunodeficiency virus: its characterization by the polymerase-chain-reaction method and failure of experimental transmission to immunodeficient animals.
Parasitol. Res.
79:624-628[CrossRef][Medline].
|
| 24.
|
Gigliotti, F.,
A. G. Harsen,
C. G. Haidaris, and P. J. Haidaris.
1993.
Pneumocystis carinii is not universally transmissible between mammalian species.
Infect. Immun.
61:2886-2890[Abstract/Free Full Text].
|
| 25.
|
Guillot, J.,
M. Berthelemy,
B. Polack,
V. Laine,
P. Lacube,
R. Chermette, and P. Roux.
1999.
Impaction versus filtration for the detection of Pneumocystis carinii DNA in air.
J. Eukaryot. Microbiol.
46:S100-S101.
|
| 26.
|
Hendley, J. O., and T. H. Weller.
1971.
Activation and transmission in rats of infection with Pneumocystis.
Proc. Soc. Exp. Biol. Med.
137:1401-1404[CrossRef][Medline].
|
| 27.
|
Hennequin, C.,
B. Page,
P. Roux,
C. Legendre, and H. Kreis.
1995.
Outbreak of Pneumocystis carinii pneumonia in a renal transplant unit.
Eur. J. Clin. Microbiol. Infect. Dis.
14:122-126[CrossRef][Medline].
|
| 28.
|
Hunter, J. A. C., and A. E. Wakefield.
1996.
Genetic divergence at the mitochondrial small subunit ribosomal RNA gene among isolates of Pneumocystis carinii from five mammalian host species.
J. Eukaryot. Microbiol.
43:S24-S25[CrossRef].
|
| 29.
|
Journal of Eukaryotic Microbiology.
1994.
The Pneumocystis Workshop.
J. Eukaryot. Microbiol.
41:S121-S122.
|
| 30.
|
Keely, S.,
H. J. Pai,
R. Baughman,
C. Sidman,
S. M. Sunkin,
J. R. Stringer, and S. L. Stringer.
1994.
Pneumocystis species inferred from analysis of multiple genes.
J. Eukaryot. Microbiol.
41:S94.
|
| 31.
|
Lane, B. R.,
J. C. Ast,
P. A. Flossler,
D. P. Mindell,
M. S. Bartlett,
J. W. Smith, and S. R. Meshnick.
1997.
Dihydropteroate synthetase polymorphism in Pneumocystis carinii.
J. Infect. Dis.
175:482-485[Medline].
|
| 32.
|
Liu, Y., and M. J. Leibowitz.
1993.
Variation and in vitro splicing of group I introns in rRNA genes of Pneumocystis carinii.
Nucleic Acids Res.
20:3763-3772[Abstract/Free Full Text].
|
| 33.
|
Mayr, E.
1963.
Animal species and evolution.
The Belknap Press of Harvard University Press, Cambridge, Mass.
|
| 34.
|
Mazars, E., and E. Dei-Cas.
1998.
Epidemiological and taxonomic impact of Pneumocystis biodiversity.
FEMS Immunol. Med. Microbiol.
22:75-80[CrossRef][Medline].
|
| 35.
|
Mazars, E.,
C. Odberg Ferragut,
E. Dei-Cas,
M. N. Fourmaux,
E. M. Aliouat,
M. Brun Pascaud,
G. Mougeot, and D. Camus.
1995.
Polymorphism of the thymidylate synthase gene of Pneumocystis carinii from different host species.
J. Eukaryot. Microbiol.
42:26-32[Medline].
|
| 36.
|
Mazars, E.,
K. Guyot,
I. Durand,
E. Dei-Cas,
S. Boucher,
S. Ben Abderrazak,
A. L. Banuls,
M. Tibayrenc, and D. Camus.
1997.
Isoenzyme diversity in Pneumocystis carinii from rats, mice and rabbits.
J. Infect. Dis.
175:655-660[Medline].
|
| 37.
|
Palmer, R. J.,
O. P. Settnes,
J. Lodal, and A. E. Wakefield.
2000.
Population structure of rat-derived Pneumocystis carinii in Danish wild rats.
In
Appl. Environ. Microbiol. 66:4954-4956.
|
| 38.
|
Saitou, N., and M. Nei.
1987.
The neighbor-joining method: a new method for reconstructing phylogenetic trees.
Mol. Biol. Evol.
4:406-425[Abstract].
|
| 39.
|
Sethi, K. K.
1992.
Multiplication of human-derived P. carinii in severe combined immunodeficient (SCID) mice.
Experimentia
48:63-67[CrossRef][Medline].
|
| 40.
|
Smyth, J. D.
1994.
Introduction to animal parasitology.
Cambridge University Press, Cambridge, England.
|
| 41.
|
Soulez, B.,
F. Palluault,
E. Dei-Cas,
E. M. Alliouat, and D. Camus.
1991.
Introduction of Pneumocystis carinii in a colony of SCID mice.
J. Protozool.
38:S123-S125.
|
| 42.
|
Stringer, J. R.
1996.
Pneumocystis carinii: what is it, exactly?
Clin. Microbiol. Rev.
9:489-498[Abstract].
|
| 43.
|
Stringer, J. R.,
S. L. Stringer,
J. Zhang,
R. Baughman,
A. G. Smulian, and M. T. Cushion.
1993.
Molecular genetic distinction of Pneumocystis carinii from rats and humans.
J. Eukaryot. Microbiol.
40:733-741[Medline].
|
| 44.
|
Thompson, J. D.,
T. S. Gibson,
F. Plewniak,
F. Jeanmougin, and D. G. Higgings.
1997.
The Clustal X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools.
Nucleic Acids Res.
25:4876-4882[Abstract/Free Full Text].
|
| 45.
|
Wakefield, A. E.
1996.
DNA sequences identical to Pneumocystis carinii f. sp. carinii and Pneumocystis carinii f. sp. hominis in samples of air spora.
J. Clin. Microbiol.
34:1754-1759[Abstract].
|
| 46.
|
Wakefield, A. E.
1998.
Genetic heterogeneity in Pneumocystis carinii: an introduction.
FEMS Immunol. Med. Microbiol.
22:5-13[CrossRef][Medline].
|
| 47.
|
Wakefield, A. E.,
J. R. Stringer,
E. Tamburrini, and E. Dei Cas.
1998.
Genetics, metabolism and host specificity of Pneumocystis carinii.
Med. Mycol.
36:S183-S193.
|
| 48.
|
Walzer, P. D.
1984.
Experimental models of Pneumocystis carinii infection, p. 37-43.
In
L. S. Young (ed.), Pneumocystis carinii pneumonia, 1st ed. M. Dekker, Inc., New York, N.Y.
|
Journal of Clinical Microbiology, June 2001, p. 2126-2133, Vol. 39, No. 6
0095-1137/01/$04.00+0 DOI: 10.1128/JCM.39.6.2126-2133.2001
Copyright © 2001, American Society for Microbiology. All rights reserved.
This article has been cited by other articles:
-
Demanche, C., Wanert, F., Barthelemy, M., Mathieu, J., Durand-Joly, I., Dei-Cas, E., Chermette, R., Guillot, J.
(2005). Molecular and serological evidence of Pneumocystis circulation in a social organization of healthy macaques (Macaca fascicularis). Microbiology
151: 3117-3125
[Abstract]
[Full Text]
-
Goldman, A. S., Goldman, L. R., Goldman, D. A.
(2005). What Caused the Epidemic of Pneumocystis Pneumonia in European Premature Infants in the Mid-20th Century?. Pediatrics
115: e725-e736
[Abstract]
[Full Text]
-
Miller, R, Huang, L
(2004). Pneumocystis jirovecii infection. Thorax
59: 731-733
[Full Text]
-
Keely, S. P., Fischer, J. M., Cushion, M. T., Stringer, J. R.
(2004). Phylogenetic identification of Pneumocystis murina sp. nov., a new species in laboratory mice. Microbiology
150: 1153-1165
[Abstract]
[Full Text]
-
Totet, A., Duwat, H., Magois, E., Jounieaux, V., Roux, P., Raccurt, C., Nevez, G.
(2004). Similar genotypes of Pneumocystis jirovecii in different forms of Pneumocystis infection. Microbiology
150: 1173-1178
[Abstract]
[Full Text]
-
Schaffzin, J. K., Stringer, J. R.
(2004). Expression of the Pneumocystis carinii major surface glycoprotein epitope is correlated with linkage of the cognate gene to the upstream conserved sequence locus. Microbiology
150: 677-686
[Abstract]
[Full Text]
-
Ambrose, H. E., Keely, S. P., Aliouat, E. M., Dei-Cas, E., Wakefield, A. E., Miller, R. F., Stringer, J. R.
(2004). Expression and complexity of the PRT1 multigene family of Pneumocystis carinii. Microbiology
150: 293-300
[Abstract]
[Full Text]
-
Norris, K. A., Wildschutte, H., Franko, J., Board, K. F.
(2003). Genetic Variation at the Mitochondrial Large-Subunit rRNA Locus of Pneumocystis Isolates from Simian Immunodeficiency Virus-Infected Rhesus Macaques. CVI
10: 1037-1042
[Abstract]
[Full Text]
-
Gigliotti, F., Harmsen, A. G., Wright, T. W.
(2003). Characterization of Transmission of Pneumocystis carinii f. sp. muris through Immunocompetent BALB/c Mice. Infect. Immun.
71: 3852-3856
[Abstract]
[Full Text]
-
Durand-Joly, I., Aliouat, E. M., Recourt, C., Guyot, K., Francois, N., Wauquier, M., Camus, D., Dei-Cas, E.
(2002). Pneumocystis carinii f. sp. hominis Is Not Infectious for SCID mice. J. Clin. Microbiol.
40: 1862-1865
[Abstract]
[Full Text]
-
Wakefield, A. E
(2002). Pneumocystis carinii: Role in childhood respiratory infections. Br Med Bull
61: 175-188
[Abstract]
[Full Text]