*Department of Microbiology and Evolutionary Biology, Faculty of Science,
and
CAOS/CAMM Center, University of Nijmegen, Nijmegen, the Netherlands;
and
Biologisches Institut, Universitat Stuttgart, Stuttgart, Germany
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Anaerobic heterotrichous ciliates are especially suited for analysis of the evolution of such endosymbiotic associations in the various environments: they are closely related to each other, but they thrive in the most divergent niches, such as marine and freshwater sediments and the intestinal tracts of animals (Hirt, Wilkinson, and Embley 1998
; van Hoek et al. 1998
, 1999
). All of the anaerobic heterotrichous ciliates studied to date possess hydrogenosomes that most likely evolved from the mitochondria of their aerobic ancestors (Embley, Horner, and Hirt 1997
; Hackstein et al. 1999
; van Hoek et al. 2000
). Therefore, one might assume that these ciliates hosted methanogenic endosymbionts from the very beginnings of their evolution because, as mentioned before, hydrogenosomes require a low intracellular hydrogen partial pressure for proper functioning.
We sequenced the small-subunit (SSU) rDNA genes of 18 anaerobic heterotrichous ciliates and their methanogenic endosymbionts in order to analyze the evolutionary relationships between anaerobic heterotrichous ciliates and their methanogenic endosymbionts. Here, we will show that coevolution between hosts and their endosymbionts can be demonstrated in only a few cases, notwithstanding a strictly vertical transmission of the methanogenic endosymbiont. We will also demonstrate that the endosymbionts from ciliates living in a particular ecological niche, e.g., freshwater sediments or intestinal tracts, respectively, are monophyletic and related to their free-living relatives. These observations strongly suggest that methanogenic endosymbionts were acquired prior to radiation but that endosymbiont replacements most likely accompanied the evolution of the anaerobic heterotrichous ciliates.
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Microscopical Procedures
Endosymbiotic methanogenic bacteria were studied with a Leitz epifluorescence microscope (Doddema and Vogels 1978
). The number of methanogenic endosymbionts was counted from color slides made from squash preparations of unfixed ciliates. The methanogenic endosymbionts were identified in situ by hybridization with a fluorescent archaea-specific probe (ArcR915, see below) according to Fokin et al. (1996)
.
Gas Chromatography
Approximately 2001,000 intestinal ciliates were isolated by electromigration and incubated in an anaerobic culture medium (Yamin 1978, modified). Methane production was measured with a Pye Unicam gas chromatograph (FID; Porapak Q column; 80/100 mesh). Ethane was used as an internal standard.
DNA Isolation, PCR Amplification, Sequencing, and Amplified Ribosomal DNA Restriction Analysis
DNA was isolated from single ciliates (van Hoek et al. 1998
). The complete 16S rDNA sequence of the methanogenic endosymbiont of Nyctotherus ovalis from Blaberus sp. var. Nijmegen was determined after PCR amplification with the primers ArcF7 (5'-GTTGATCCTGCCAGAGG-3') and ArcR1326 (5'-TGTGTGCAAGGAGCAGGGAC-3'). For all other methanogenic endosymbionts, partial 16S rDNA sequences (approximately 900 nt) were obtained with the primers ArcF26 (5'-TTGGGTTCGATTAAGCCAT-3') and ArcR915 (5'-GTGCTCCCCCGCCAATTCCT-3'). PCR amplification of the 18S rDNA genes of the ciliates was performed using the primers euk-forward and euk-1200R or euk-300F and euk-reverse (van Hoek et al. 1998
, 1999
).
PCR products of the 16S and 18S rRNA genes of the ciliates and their methanogenic endosymbionts were sequenced directly in both directions using an ABI Prism 310 Automated DNA sequencer (Perkin Elmer; van Hoek et al. 1998
). For the determination of the complete rDNA sequence of the above-mentioned methanogenic endosymbiont, the PCR product was cloned using the TA Cloning Kit (Invitrogen) and INV
F One Shot competent Escherichia coli. Two clones were sequenced using six different methanogenic-archaea-specific and two plasmid primers. The DNA sequences were deposited in the EMBL database (accession numbers for 16S rRNA: AJ132639AJ132655, AJ238131).
Amplified 16S rDNA of one to three individuals of all freshwater ciliates and two to nine individuals of the intestinal ciliates was subjected to amplified ribosomal DNA restriction analysis (ARDRA; van Hoek et al. 1998
) with the enzymes HaeIII, MspI, and NciI (GibcoBRL), respectively, in order to confirm the validity of DNA sequence analysis. Also, the homogeneity of the endosymbiont populations was tested by ARDRA. For each freshwater sampling place, individual representatives of the various ciliate genera were subjected to ARDRA analysis of both 18S rDNA (of the host) and 16S rDNA (of the methanogenic endosymbionts) to validate the endosymbiont-host relationships. PCR was also used to control contamination of the washing buffers by environmental methanogenic archaea. Contaminated samples were excluded from further analysis.
Phylogenetic Analysis
DNA sequences were aligned using the PileUp program of the Wisconsin package, version 8.1 (Genetics Computer Group [GCG], Madison, Wis.; see also van Hoek et al. 1998
). For distance analyses, alignments of 1,243 nt of the 18S rRNA genes and 770 nt of the 16S rRNA genes were used. In order to test the reliability of the different radiations of anaerobic ciliates, phylogenetic analysis was also performed with a data set that contained the concatenated partial 18S and complete ITS1, ITS2, and 5.8S rDNAs from at least one representative of the potential branches. EDNADIST (Rice et al. 1995
) was used to calculate the sequence similarity and evolutionary distances using the Jukes and Cantor (1969)
nucleotide substitution model. Phylogenetic trees were constructed using the neighbor-joining method (Saitou and Nei 1987
). The NUCML program from the MOLPHY package (version 2.2; Adachi and Hasegawa 1996
), PUZZLE (Strimmer and von Haeseler 1996
), and likelihood mapping (Strimmer and von Haeseler 1997
) were used to perform maximum-likelihood (ML) analysis. The TreeMap program (Page 1995
), obtained from the Internet site http://taxonomy.zoology.gla.ac.uk/rod/treemap, was used to test for cospeciation of the hosts and their symbionts.
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
|
Molecular Phylogenies of Anaerobic Heterotrichous Ciliates and Their Methanogenic Endosymbionts
Phylogenetic analysis of the 18S rRNA genes of 18 anaerobic heterotrichous ciliates with the neighbor-joining method revealed monophyly for the intestinal ciliates and the existence of two to five paraphyletic clusters of the anaerobic, heterotrichous ciliates from freshwater sediments and intestinal tracts (fig. 2 ). The Plagiopyla and Trimyema species form a separate branch. This conclusion is also supported by ML analysis using NUCML (star-decomposition and quick-add-OTU options), quartet-puzzling, and likelihood mapping (see Materials and Methods). Moreover, using concatenated partial 18S, 5.8S, ITS1, and ITS2 rDNA sequences from at least one representative for each branch of the phylogenetic tree, the paraphly of the Caenomorpha species, certain Caenomorpha-like species, Brachonella, and Nyctotherus species could be confirmed. Only the branching order of the Brachonella and Caenomorpha-like species could not be resolved completely when using ML analysis. Consistently, this analysis supported clustering of the Nyctotherus species from frogs and millipedes on the one side and the N. ovalis species from the cockroach guts on the other side.
|
DNA sequence analysis of the 16S rRNA genes consistently revealed that only one rDNA of methanogenic origin could be amplified from a single ciliate. Also, ARDRA analysis of the amplified DNA confirmed the homogeneity of the rDNA. As a rule, ciliates from one and the same isolate hosted identical endosymbionts. Also, ciliates that had identical 18S rRNA genes but were from different environments (i.e., from Blaberus sp. var. Amsterdam and Periplaneta americana var. Bayer) possessed identical methanogens. However, DNA sequence analysis of the 16S rRNA genes of the endosymbionts of the majority of intestinal and free-living ciliates studied here revealed a significant degree of divergence between the endosymbionts. Notably, phylogenetic analysis of the 16S rRNA genes of these endosymbionts, together with 163 sequences from environmental methanogens, revealed that the endosymbionts of the intestinal ciliates cluster among the Methanobacteriales, whereas those of the free-living ciliates cluster among the Methanomicrobiales (fig. 3 ). Also, three DNA sequences from endosymbionts from Plagiopyla nasuta, Metopus contortus, and Trimyema compressum cluster with Methanomicrobiales. Three more DNA sequences from Plagiopyla frontata, Metopus striatus, and Metopus palaeformis cluster with different methanogens from environmental sources.
|
|
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
The persistence of the endosymbionts through cell divisions and cyst (resting) stages permits a vertical transmission of the methanogenic endosymbionts. There is no evidence for an environmental transmission of the endosymbionts, since transfaunation experiments with ciliates possessing morphologically distinct endosymbionts (cf. fig. 1a, c, and d
) failed to reveal an uptake of endosymbiotic methanogens from the environment (van Hoek et al. 1999
). Also, the observation that closely related ciliate species from the same sampling place (Dekkerswald; fig. 4
) possess clearly divergent endosymbionts suggests a strictly vertical transmission of the endosymbionts and a high specificity of the association. Also, the fact that identical ciliates from different environments host identical endosymbionts (fig. 4
) argues for a vertical transmission. Moreover, the endosymbionts of ciliates living in the intestinal tracts of frogs (Rana ridibunda and Rana temporaria; figs. 3
and 4
) belong to the "intestinal" type of endosymbionts, notwithstanding that a wealth of free-living methanogens populates the particular habitat of the frogs and their larvae. Thus, all available evidence argues for a high specificity of the association and a long-lasting genetic isolation of the ciliates and their endosymbionts (fig. 4
).
Comparable symbiotic associations exist between prokaryotes and insects, e.g., between the eubacterium Buchnera and many aphids (Baumann et al. 1995
; Baumann, Moran, and Baumann 1997
), and between flavobacteria and their termite and cockroach hosts (Bandi et al. 1994
, 1997
). Phylogenetic analysis of these associations revealed a complete congruency between host and symbiont phylogenies. Therefore, it was very surprising that the phylogenetic analysis of the SSU rRNA genes of 18 anaerobic heterotrichous ciliates and their methanogenic endosymbionts did not provide any evidence for a significant match between the phylogenies of the ciliate hosts and their symbionts (fig. 4 ). Such an incongruence between host and symbiont phylogenies can have many reasons, for example, insufficient phylogenetic information (Distel 1998
; Moran and Telang 1998
). With our data, it is unlikely that the tree topology is biased by the use of partial SSU rDNA sequences: a phylogenetic analysis that includes the internal transcribed spacer (ITS-1 and ITS-2) DNA sequences and the 5.8S rDNA genes of the ciliates generates trees with an identical topology. We also have no indication that the topology of the phylogenetic trees is biased by the alignments or the tree-building algorithms (fig. 4
).
We cannot exclude the possibility that the molecular clocks of the eukaryotic hosts and their prokaryotic endosymbionts are different and responsible for the low degree of congruence of host and symbiont trees (Hafner et al. 1994
; Page and Hafner 1996
). However, the observed relationship between the endosymbionts and their free-living relatives strongly suggests that host switches or endosymbiont replacements occurred, albeit rarely, during the evolution of these symbiotic associations. The main argument for this hypothesis stems from the observation that the endosymbionts of ciliates from freshwater sediment and those from the intestinal tracts of animals are monophyletic but clearly different in the two ecosystems. Also, the endosymbionts of Plagiopyla and Trimyema are related to free-living methanogens from the particular environments. The observed vertical transmission, the degree of DNA sequence divergence of the methanogenic endosymbionts, and the clustering of the endosymbionts in the phylogenetic tree (fig. 4
) exclude a frequent uptake of methanogens from the environment. Notably, there is also no evidence for the presence of substantial numbers of identical or closely related free-living methanogens in the freshwater sediments studied. In the intestinal tract of the cockroach P. americana as well, there is no evidence for the presence of a Methanobrevibacter-like methanogen (W. Sprenger, personal communication; PCR and ribotype analysis, not shown). Such a distribution resembles the situation in certain bivalves that host chemoautotrophic eubacterial endosymbionts (Distel 1998
; Krueger and Cavanaugh 1997
). Also for these symbioses, a symbiont replacement has been postulated (Krueger and Cavanaugh 1997
). In the case of the anaerobic, heterotrichous ciliates and their endosymbiotic methanogens, physiology and phylogeny clearly support the assumption of endosymbiont replacements. Figure 5
shows that at least two (but potentially as many as four) different clades of anaerobic heterotrichous ciliates evolved in the freshwater sediments, where they acquired closely related methanogenic endosymbionts. These endosymbionts cluster among the free-living "environmental" Methanomicrobiales. The Nyctotherus species that inhabit the intestinal tracts of the various animals are monophyletic and host endosymbionts that are related to intestinal Methanobacteriales. The phylogeny of the hydrogenosomes of N. ovalis suggests an origin of these organelles from the mitochondria of their aerobic ancestors (van Hoek et al. 2000
).
|
![]() |
Sequence Availability |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
1 Abbreviations: ARDRA, amplified ribosomal DNA restriction analysis; ITS, internal transcribed spacer; LSM, laser scanning microscope.
2 Keywords: anaerobic ciliates
methanogenic archaea
symbiosis
endosymbiont replacement
hydrogenosomes
coevolution
3 Address for correspondence and reprints: J. Hackstein, Department of Microbiology and Evolutionary Biology, Toernooiveld 1, NL-6525ED Nijmegen, the Netherlands. E-mail: hack{at}sci.kun.nl
![]() |
literature cited |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Adachi, K., and M. Hasegawa. 1996. MOLPHY version 2.3: programs for molecular phylogenetics based on maximum likelihood. Comput. Sci. Monogr. 28:1150.
Akhmanova, A., F. Voncken, T. van Alen, A. van Hoek, B. Boxma, G. Vogels, M. Veenhuis, and J. H. P. Hackstein. 1998. A hydrogenosome with a genome. Nature 396:527528.
Bandi, C., G. Damiani, L. Magrassi, A. Grigolo, R. Fani, and L. Sacchi. 1994. Flavobacteria as intracellular symbionts in cockroaches. Proc. R. Soc. Lond. B Biol. Sci. 257:4248.
Bandi, C., M. Sironi, C. A. Nalepa, S. Corona, and L. Sacchi. 1997. Phylogenetic distant intracellular symbionts in termites. Parasitologia 39:7175.
Baumann, P., C.-Y. Lai, L. Baumann, D. Rouhbakhsh, N. A. Moran, and M. A. Clark. 1995. Mutualistic associations of aphids and prokaryotes: biology of the genus Buchnera. Appl. Environ. Microbiol. 61:17.
Baumann, P., N. A. Moran, and L. Baumann. 1997. The evolution and genetics of aphid endosymbionts. Bioscience 47:1220.
Biagini, G. A., B. J. Finlay, and D. Lloyd. 1997. Evolution of the hydrogenosome. FEMS Microbiol. Lett. 155:133140.[ISI][Medline]
Brynnel, E. U., C. G. Kurland, N. A. Moran, and S. G. E. Andersson. 1998. Evolutionary rates for tuf genes in endosymbionts of aphids. Mol. Biol. Evol. 15:574582.[Abstract]
Distel, D. L. 1998. Evolution of chemoautotrophic endosymbioses in bivalves. Bivalve-bacteria chemosymbioses are phylogenetically diverse but morphologically similar. BioScience 48:277286.
Doddema, H. J., and G. D. Vogels. 1978. Improved identification of methanogenic bacteria by fluorescence microscopy. Appl. Environ. Microbiol. 36:752754.[ISI][Medline]
Doolittle, W. F. 1998. You are what you eat: a gene transfer ratchet could account for bacterial genes in eukaryotic nuclear genomes. Trends Genet. 14:307311.[ISI][Medline]
Embley, T. M., B. J. Finlay, P. L. Dyal, R. P. Hirt, M. Wilkinson, and A. G. Williams. 1995. Multiple origins of anaerobic ciliates with hydrogenosomes within the radiation of aerobic ciliates. Proc. R. Soc. Lond. B Biol. Sci. 262:8793.[ISI][Medline]
Embley, T. M., D. A. Horner, and R. P. Hirt. 1997. Anaerobic eukaryote evolution: hydrogenosomes as biochemically modified mitochondria? Trends Ecol. Evol. 12:437441.
Embley, T. M., and W. Martin. 1998. A hydrogen-producing mitochondrion. Nature 396:517519.
Felsenstein, J. 1985. Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39:783791.
Fenchel, T., and B. J. Finlay. 1995. Ecology and evolution in anoxic worlds. Oxford University Press, Oxford, New York, and Tokyo.
Fokin, S. I., T. Brigge, J. Brenner, and H.-D. G;auortz. 1996. Holospora species infecting the nuclei of Paramecium appear to belong in two groups of bacteria. Eur. J. Protistol. 32:1924.
Hackstein, J. H. P., A. S. Akhmanova, B. Boxma, H. R. Harhangi, and F. Voncken. 1999. Hydrogenosomes: eukaryotic adaptations to anaerobic environments. Trends Microbiol. 7:441447.[ISI][Medline]
Hafner, M. S., P. D. Sudman, F. X. Villablanca, T. A. Spradling, J. W. Demastes, and S. A. Nadler. 1994. Disparate rates of molecular evolution in cospeciating hosts and parasites. Science 265:10871090.
Hirt, R. P., P. L. Dyal, M. Wilkinson, B. J. Finlay, D. M. Roberts, and T. M. Embley. 1995. Phylogenetic relationships among karyorelictids and heterotrichs inferred from small subunit rRNA sequences: Resolution at the base of the ciliate tree. Mol. Phylogenet. Evol. 4:7787.[Medline]
Hirt, R. P., M. Wilkinson, and T. M. Embley. 1998. Molecular and cellular evolution of ciliates: a phylogenetic perspective. Pp. 327340 in G. H. Coombs, K. Vickerman, M. A. Sleigh, and A. Warren, eds. Evolutionary relationships among protozoa. Chapman and Hall, London.
Hurst, L. D., and G. T. McVean. 1996. ... and scandalous symbionts. Nature 381:650651.
Jukes, T. H., and C. R. Cantor. 1969. Evolution of protein molecules. Pp. 21132 in H. N. Munro, ed. Mammalian protein metabolism. Academic Press, New York.
Krueger, D. M., and C. M. Cavanaugh. 1997. Phylogenetic diversity of bacterial symbionts of Solemya hosts based on comparative sequence analysis of 16S rRNA genes. Appl. Environ. Microbiol. 63:9198.[Abstract]
Martin, W., and M. M;auuller. 1998. The hydrogen hypothesis for the first eukaryote. Nature 392:3741.
Moran, N. A. 1996. Accelerated evolution and Muller's ratchet in endosymbiotic bacteria. Proc. Natl. Acad. Sci. USA 93:28732878.
Moran, N. A., and A. Telang. 1998. Bacteriocyte-associated symbionts of insects. Bioscience 48:295304.
Page, R. D. M. 1995. Paralell phylogenies: reconstructing the history of host-parasite assemblage. Cladistics 10:155173.
Page, R. D. M., and M. S. Hafner. 1996. Molecular phylogenies and host-parasite cospeciation: gophers and lice as a model system. Pp. 255270 in P. H. Harvey, A. J. Leigh-Brown, J. Maynard-Smith, and S. Nee, eds. New uses for new phylogenies. Oxford University Press, Oxford, New York, and Tokyo.
Rice, P., R. Lopez, R. Doelz, and J. A. M. Leunissen. 1995. EGCG 8.0. EMBnet.news 2(2):57.
Saitou, N., and M. Nei. 1987. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4:406425.[Abstract]
Strimmer, K., and A. von Haeseler. 1996. Quartet puzzling: a quartet maximum-likelihood method for reconstructing tree topologies. Mol. Biol. Evol. 13:964969.
. 1997. Likelihood-mapping: a simple method to visualize phylogenetic content of a sequence alignment. Proc. Natl. Acad. Sci. USA 94:68156819.
van Hoek, A. H. A. M., A. S. Akhmanova, M. A. Huynen, and J. H. P. Hackstein. 2000. A mitochondrial ancestry of the hydrogenosomes of Nyctotherus ovalis. Mol. Biol. Evol. 17:202206.
van Hoek, A. H. A. M., V. S. I. Sprakel, T. A. van Alen, A. P. R. Theuvenet, G. D. Vogels, and J. H. P. Hackstein. 1999. Voltage-dependent reversal of anodic galvanotaxis in Nyctotherus ovalis. J. Eukaryot. Microbiol. 46:427433.
van Hoek, A. H. A. M., T. A. van Alen, V. S. I. Sprakel, J. H. P. Hackstein, and G. D. Vogels. 1998. Evolution of anaerobic ciliates from the gastrointestinal tract: phylogenetic analysis of the ribosomal repeat from Nyctotherus ovalis and its relatives. Mol. Biol. Evol. 15:11951206.[Abstract]
Yamin, M. A. 1978. Axenic cultivation of the cellulolytic flagellate Trichomitopsis termopsidis (Cleveland) from the termite Zootermopsis. J. Protozool. 25:535538.