This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Blazejak, A.
Right arrow Articles by Dubilier, N.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Blazejak, A.
Right arrow Articles by Dubilier, N.
Agricola
Right arrow Articles by Blazejak, A.
Right arrow Articles by Dubilier, N.

 Previous Article  |  Next Article 

Applied and Environmental Microbiology, August 2006, p. 5527-5536, Vol. 72, No. 8
0099-2240/06/$08.00+0     doi:10.1128/AEM.02441-05
Copyright © 2006, American Society for Microbiology. All Rights Reserved.

Phylogeny of 16S rRNA, Ribulose 1,5-Bisphosphate Carboxylase/Oxygenase, and Adenosine 5'-Phosphosulfate Reductase Genes from Gamma- and Alphaproteobacterial Symbionts in Gutless Marine Worms (Oligochaeta) from Bermuda and the Bahamas

Anna Blazejak,1 Jan Kuever,1,{dagger} Christer Erséus,2 Rudolf Amann,1 and Nicole Dubilier1*

Max Planck Institute of Marine Microbiology, Celsiusstrasse 1, D-28359 Bremen, Germany,1 Department of Zoology, Göteborg University, Box 463, SE-405 30 Göteborg, Sweden2

Received 16 October 2005/ Accepted 17 April 2006


arrow
ABSTRACT
 
Gutless oligochaetes are small marine worms that live in obligate associations with bacterial endosymbionts. While symbionts from several host species belonging to the genus Olavius have been described, little is known of the symbionts from the host genus Inanidrilus. In this study, the diversity of bacterial endosymbionts in Inanidrilus leukodermatus from Bermuda and Inanidrilus makropetalos from the Bahamas was investigated using comparative sequence analysis of the 16S rRNA gene and fluorescence in situ hybridization. As in all other gutless oligochaetes examined to date, I. leukodermatus and I. makropetalos harbor large, oval bacteria identified as Gamma 1 symbionts. The presence of genes coding for ribulose-1,5-bisphosphate carboxylase/oxygenase form I (cbbL) and adenosine 5'-phosphosulfate reductase (aprA) supports earlier studies indicating that these symbionts are chemoautotrophic sulfur oxidizers. Alphaproteobacteria, previously identified only in the gutless oligochaete Olavius loisae from the southwest Pacific Ocean, coexist with the Gamma 1 symbionts in both I. leukodermatus and I. makropetalos, with the former harboring four and the latter two alphaproteobacterial phylotypes. The presence of these symbionts in hosts from such geographically distant oceans as the Atlantic and Pacific suggests that symbioses with alphaproteobacterial symbionts may be widespread in gutless oligochaetes. The high phylogenetic diversity of bacterial endosymbionts in two species of the genus Inanidrilus, previously known only from members of the genus Olavius, shows that the stable coexistence of multiple symbionts is a common feature in gutless oligochaetes.


arrow
INTRODUCTION
 
Gutless oligochaetes are small worms about 0.1 to 0.3 mm in diameter and 2 to 50 mm long that occur worldwide in marine sediments, with the highest diversity being found in tropical and subtropical coral reefs (11-15, 17, 20, 21). They form a monophyletic group that consists of only two genera, Inanidrilus and Olavius, within the subfamily Phallodrilinae (31). All worm species described so far, 25 Inanidrilus and 56 Olavius species (12), live in an obligate association with endosymbiotic bacteria. The lack of both a digestive and an excretory system led to the assumption that the symbiotic bacteria provide their hosts with a source of nutrition (19). Enzyme assays and uptake experiments with inorganic carbon indicated that at least some of the bacterial symbionts are thiotrophic; i.e., they use reduced sulfur compounds to fix CO2 into organic carbon compounds (16). It is assumed that the transfer of these organic compounds to the host is the main mode of energy transfer, although digestion of the bacteria may also supply nutrients.

The bacterial symbionts occur in a thick layer just below the cuticle of the worm in an extracellular space above the epidermal cells of the worm (23). A phylogenetically diverse assemblage of bacteria co-occurs within the symbiotic layer, with up to six rRNA phylotypes being identified in a single host species (3). All gutless oligochaetes harbor a large (3- by 8-µm) thiotrophic symbiont with numerous sulfur and poly-beta-hydroxyalkanoate inclusions (4). These symbionts are found in all host species and form a closely related cluster of 16S rRNA sequences within the Gammaproteobacteria (8). Coexisting with these primary symbionts are smaller bacteria (0.7 by 1.9 µm) without any conspicuous inclusions that can belong to the Gamma-, Alpha-, or Deltaproteobacteria (3, 7, 10). While the metabolism of the gamma- and alphaproteobacterial symbionts is not clear, the deltaproteobacterial symbionts have been identified as sulfate reducers and are assumed to be engaged in a syntrophic sulfur cycle with the thiotrophic symbiont (3, 10). In some host species, spirochetes have also been identified as members of the oligochaete symbiont community (3, 7, 22).

The phylogeny of symbionts has been described for three host species from the genus OlaviusOlavius loisae, from Australia (7), Olavius algarvensis, from the Mediterranean Sea (10), and Olavius crassitunicatus, from Peru (3)—but for only a single species from the genus Inanidrilus, Inanidrilus leukodermatus (9). I. leukodermatus occurs in high abundance in calcareous sands of coral reefs around the island Bermuda in the northwest Atlantic Ocean, and the morphology, physiology, and ecology of this species have been studied intensively (9, 16, 19, 20, 23, 25). The symbionts of this species were the first to be characterized using molecular methods, and only a single thiotrophic gammaproteobacterial phylotype was identified despite the presence of multiple bacterial morphotypes (9). More recent studies on other gutless oligochaete species using improved methods of 16S rRNA analysis and fluorescence in situ hybridization (FISH) have demonstrated that their multiple bacterial morphotypes correspond to multiple bacterial phylotypes (8). In the present study, we therefore reexamined the symbionts in I. leukodermatus and extended these analyses to include another Inanidrilus species from the northwest Atlantic Ocean, Inanidrilus makropetalos, that occurs in calcareous sands of coral reefs in the Bahamas (12).

In addition to the phylogenetic characterization of the Inanidrilus symbionts, we investigated genes coding for ribulose 1,5-bisphosphate carboxylase/oxygenase (RubisCO), the key enzyme for autotrophic CO2 fixation in the Calvin-Benson-Bassham cycle, and adenosine 5'-phosphosulfate (APS) reductase, an enzyme involved in sulfur metabolism. Two distinct forms of RubisCO have been found in chemoautotrophic symbionts: form I, for which the cbbL gene is used as a functional marker, and form II, for which the cbbM gene is used (5, 25, 33, 38, 39). Form I RubisCO differs from form II structurally and is better adapted to CO2 fixation under aerobic conditions (40, 41).

APS reductase is used in both the reductive and oxidative modes of sulfur metabolism and is therefore found in both sulfate reducers and sulfur oxidizers (18, 24). In sulfate reducers, APS reductase catalyzes the two-electron reduction of APS to sulfite and AMP, and in sulfur oxidizers, it catalyzes the reverse reaction (35). APS reductase consists of an alpha and beta subunit, encoded by the genes aprA and aprB, respectively, that form a {alpha}1ß1 heterodimer. The aprA gene has been proposed as a useful phylogenetic marker for bacteria involved in oxidative and reductive sulfur metabolism (24). The phylogeny of the aprA gene has been well studied in sulfate reducers (18), but nothing is known about the phylogeny of aprA from sulfur-oxidizing bacteria. Only a few sequences from free-living sulfur oxidizers are currently available in the GenBank database, and prior to this study, none from symbiotic sulfur oxidizers were known.


arrow
MATERIALS AND METHODS
 
Specimen collection.
Inanidrilus leukodermatus specimens were collected in June 1998 in Harrington Sound, Bermuda (Fig. 1). Inanidrilus makropetalos specimens were collected in April 1999 at Lee Stocking Island in the Bahamas (Fig. 1). The worms were extracted from shallow-water (<3 m) sediments by decantation with seawater and identified under a microscope. Specimens were fixed in 95% ethanol for DNA analyses and for FISH as described previously (9) and stored at 4°C.


Figure 1
View larger version (22K):
[in this window]
[in a new window]
 
FIG. 1. Map of sampling sites in the northwest Atlantic Ocean. I. leukodermatus worms were collected in Bermuda, and I. makropetalos specimens were collected at Lee Stocking Island (the Bahamas).

DNA preparation.
Six I. leukodermatus and three I. makropetalos individuals were prepared individually for PCR. Specimens were rinsed three times in MilliQ water, and DNA was isolated as described by Schizas et al. (36) in a protocol using proteinase K for digestion and the reagent GeneReleaser (BioVentures, Murfreesboro, Tenn.) for DNA purification.

PCR amplification. (i) 16S rRNA gene.
Amplifications were performed with primers specific for the bacterial 16S rRNA gene (8F and 1492R [28]) using Taq DNA polymerase (Eppendorf, Hamburg, Germany). Template DNA (1 to 2 µl) was added after preheating the PCR mix (100 µl, total volume) to 80°C to avoid nonspecific annealing of the primers to nontarget DNA. The following thermocycling conditions were used: 1 cycle at 80°C for 5 min; 27 cycles at 95°C for 1 min, 40°C for 1 min, and 72°C for 3 min; and 1 cycle at 72°C for 10 min.

(ii) cbbL and aprA genes.
The genes cbbL, coding for RubisCO form I, cbbM, for RubisCO form II, and aprA, for APS reductase, were amplified using Taq DNA polymerase (Eppendorf, Hamburg, Germany). For specific amplification of the genes coding for RubisCO forms I and II, the following primers were designed from available cbbL and cbbM sequences in the GenBank database: cbbLF (5'-CACCTGGACCACVGTBTGG-3') and cbbLR (5'-CGGTGYATGTGCAGCAGCAT5CCG-3') for the cbbL gene and cbbMF (5'-ATCATCAARCCSAARCTSGGYCTGCG-3') and cbbM1R (5'-GAGGTSACSGCRCCRTGRCCRGCMCGRTG-3') or cbbM2R (5'-GAGGTSACSGCRCCRTGRCC-3') for the cbbM gene. For specific amplification of the aprA gene, the primers aps1F (5'-TGGCAGATCATGATYMAYGG-3') and aps4R (5'-GCGCCAACYGGRCCRTA-3') (J. Kuever, unpublished data) were used. Template DNA (0.5 µl) was added to the PCR mix (50 µl, total volume). The following thermocycling conditions were applied: 1 cycle at 95°C for 3 min; 27 cycles at 95°C for 1 min, 48°C (for cbbL and cbbM genes) or 54°C (for aprA genes) for 1 min, and 72°C for 3 min; and 1 cycle at 72°C for 5 min. The PCR products from each individual were purified using a QIAquick PCR purification kit (QIAGEN, Hilden, Germany) and directly sequenced in both directions (~696 bp for cbbL genes and ~396 bp for aprA genes).

Cloning and sequencing.
16S rRNA PCR products from each individual were cloned separately using a TA cloning kit (Invitrogen, Breda, The Netherlands) according to the manufacturer's protocol. For screening of 16S rRNA genes, clones were randomly picked and selected for the correct insert size (~1,500 bp) by PCR with the vector primers M13F and M13R. For template DNA, a small amount of cells from each clone colony was picked with a sterile toothpick and resuspended in 10 µl sterile water. After preheating of this suspension to 95°C for 5 min, template DNA was amplified by PCR as described above using a total volume of 30 µl. Clones with PCR products of the correct size were prepared using the MontagePlasmid Miniprep96 kit (Millipore, Bedford, Mass.) and screened by partial sequencing of 300 to 500 bp using the 16S rRNA primer GM1F. Sequencing reactions were run using ABI BigDye on an ABI Prism 3100 genetic analyzer (Applied Biosystems, Foster City, Calif.). Sequences were aligned and compared using the Bioedit program (www.mbio.ncsu.edu/BioEdit/bioedit.html). Sequences were grouped together in a clone family if they shared ≥ 99% sequence identity (percent identical nucleotides). For each host individual, a representative clone from each clone family was fully sequenced in both directions.

Phylogenetic analyses.
The 16S rRNA, cbbL, and aprA sequences were checked against sequences in GenBank using BLAST (1) for similarity searches. Chimeras were identified using CHIMERA_CHECK 2.7 from the Ribosomal Database Project (6) and by eye in sequence alignments and excluded from further analysis. The sequence data were analyzed using the ARB software package (www.arb-home.de).

Sequence similarity of the nucleotide sequences was calculated by distance analysis with the Jukes-Cantor correction. Phylogenetic trees for 16S rRNA sequences were estimated using parsimony, distance, and maximum-likelihood analyses with different filter sets, considering only sequences with a length of at least 1,255 bp.

Phylogenetic trees for the cbbL gene were generated from amino acid sequences using maximum-likelihood analyses with a 25% amino acid frequency filter. The leading and tailing amino acid stretches of full-length sequences were excluded from the analyses, eliminating those regions found in some but not all RubisCO sequences. The partial RubisCO form I sequences (cbbL, 211 amino acids) from this study were added to the maximum-likelihood tree using maximum parsimony. The phylogeny of the aprA gene was calculated from partial sequences (131 amino acids) using maximum-likelihood analyses with a 25% amino acid frequency filter.

FISH.
Six I. leukodermatus and three I. makropetalos specimens were prepared for FISH analyses of bacterial endosymbionts as described previously (3). Symbionts were detected by CARD (catalyzed reporter deposition) FISH with horseradish peroxidase (HRP)-labeled probes and tyramide signal amplification as described by Schönhuber et al. (37). Tissue sections were hybridized with the HRP-labeled probe for 3 h at 35°C. After washing for 15 min at 35°C in washing buffer (37), the sections were equilibrated for 15 min at room temperature in phosphate-buffered saline (PBS) buffer (pH 7.3). The moist tissue sections were incubated with amplification solution (1x phosphate-buffered saline, pH 7.3; 0.0015% [vol/vol] H2O2; and 1% Alexa Fluor 488, 546, or 633 dye [Molecular Probes, Leiden, The Netherlands]) for 10 min at 37°C in the dark and rinsed in PBS buffer for 15 min at room temperature. After air drying, tissue sections were embedded in the mounting fluid Vecta Shield (Vecta Laboratories, Burlingame, CA) and stored for microscopic evaluation at –20°C for <1 to 2 days. For dual and triple hybridizations, the CARD FISH protocol was repeated two or three times on the same sections using different probes and Alexa dyes. For this purpose, after the last washing step the tissue sections were covered with 0.01 M HCI for 10 min at room temperature to inactivate the HRP. After washing for 3 min in sterile water, tissue sections were hybridized with another probe as just described.

The oligonucleotide probes designed in this study target 16S rRNA sequences isolated from I. leukodermatus and I. makropetalos (Table 1). The probes were checked against sequences in GenBank using BLAST (1) and against small-subunit rRNA sequences in the Ribosomal Database Project using PROBE MATCH (32) and contained at least one mismatch relative to all entered sequences. The general bacterial probe EUB338, the general gammaproteobacterial probe GAM42a, and the general alphaproteobacterial probe ALF968 were used as positive controls, and the antisense probe NON338 was used as a negative control. All hybridizations were performed at formamide concentrations ensuring high specificity (Table 1).


View this table:
[in this window]
[in a new window]
 
TABLE 1. Oligonucleotide probes used in this study

Nucleotide sequence accession numbers.
The 16S rRNA sequences from Inanidrilus leukodermatus were submitted to GenBank under accession numbers AJ890100 (Gamma 1 symbiont), AJ890099 (Alpha 1a symbiont), AJ890098 (Alpha 1b symbiont), AJ890093 (Alpha 1c clone associated with I. leukodermatus), and AJ890097 (Alpha 3 symbiont), and those from Inanidrilus makropetalos were submitted under accession numbers AJ890094 (Gamma 1 symbiont), AJ890095 (Alpha 1b symbiont), and AJ890096 (Alpha 2 symbiont). Sequences of genes coding for RubisCO and APS reductase from Inanidrilus leukodermatus symbionts were submitted to GenBank under accession numbers AM228902 (aprA gene) and AM228899 (cbbL gene), and those from Inanidrilus makropetalos symbionts were submitted under accession numbers AM228901 (aprA gene) and AM228900 (cbbL gene).


arrow
RESULTS
 
Clone library analysis.
Bacterial 16S rRNA sequences from six I. leukodermatus specimens were grouped into five distinct clone families, with a total of 510 clones analyzed (Table 2). Within each clone family, sequence similarity was never less than 99.6%. Phylogenetic analyses (see below) revealed that in I. leukodermatus one clone family belongs to the Gammaproteobacteria (Gamma 1) and four to the Alphaproteobacteria (Alpha 1a, 1b, 1c, and 3). Sequences belonging to the Gamma 1 clone family were found in all individuals. For the Alpha 1a, Alpha 1b, and Alpha 3 clone families, at least three out of six I. leukodermatus individuals harbored sequences belonging to these groups. In contrast, sequences belonging to the Alpha 1c clone family were present in only a single individual.


View this table:
[in this window]
[in a new window]
 
TABLE 2. 16S rRNA clone libraries from six I. leukodermatus individuals and three I. makropetalos individuals

Bacterial 16S rRNA sequences from the three I. makropetalos specimens (318 clones analyzed) were grouped into three distinct clone families, with at least 99.6% sequence similarity within each clone family (Table 2). Phylogenetic analyses (see below) confirmed that in I. makropetalos one clone family belongs to the Gammaproteobacteria (Gamma 1) and two to the Alphaproteobacteria (Alpha 1b and 2). Sequences belonging to the Gamma 1 and Alpha 2 clone families were found in all three I. makropetalos individuals and from the Alpha 1b clone family in two out of the three specimens examined.

Phylogenetic analyses 16S rRNA gene.
Parsimony, distance, and maximum-likelihood analyses of the 16S rRNA sequences from the five I. leukodermatus and three I. makropetalos clone families confirmed that these belong to one bacterial group of the Gammaproteobacteria and five phylogenetically distinct bacterial groups of the Alphaproteobacteria. The five I. leukodermatus and three I. makropetalos sequences are unique to these hosts and differ from those of symbionts from other host species or free-living bacteria.

Gammaproteobacterial symbionts.
The 16S rRNA sequence of the I. leukodermatus Gamma 1 symbiont from this study differed from the I. leukodermatus Gamma sequence described in an earlier study (8) by three nucleotides (positions 82, 861, and 862 based on E. coli numbering). All three nucleotide discrepancies are in stem regions, and these led to mismatched stem structures in our earlier sequence (8), while the stems are conserved in the sequence from this study. Furthermore, these sites are conserved in all Gamma 1 symbionts and the nucleotides in this study are in consensus with all other Gamma 1 symbionts, suggesting that these three substitutions were caused by PCR or sequencing error in the earlier study (8).

In all phylogenetic analyses, the I. makropetalos and the I. leukodermatus Gamma 1 sequences consistently group together (98.5% sequence similarity) (Fig. 2a). The Inanidrilus Gamma 1 sequences group in the same cluster as the sequences of Gamma 1 symbionts from Olavius hosts and the ectosymbiont of the marine nematode Laxus sp. (≥95.7% sequence similarity) (Fig. 2a). These symbiotic sequences, together with a clone sequence isolated from the Kazan mud volcano in the eastern Mediterranean Sea, are most closely related to a clade of free-living, phototrophic, sulfur-oxidizing bacteria from the family Chromatiaceae.


Figure 2
View larger version (49K):
[in this window]
[in a new window]
 
FIG. 2. Phylogenetic placement of bacterial symbionts in I. leukodermatus and I. makropetalos based on 16S rRNA sequences. Maximum-likelihood trees of members of the (a) Gammaproteobacteria and (b) Alphaproteobacteria are shown. Symbionts of gutless oligochaetes are listed in bold type, with the I. leukodermatus symbionts boxed in light gray and I. makropetalos symbionts in dark gray. The bars represent 10% estimated sequence divergence.

Alphaproteobacterial symbionts.
The alphaproteobacterial sequences from I. leukodermatus and I. makropetalos belong to five phylogenetically distinct bacterial groups: Alpha 1a, 1b, and 1c, Alpha 2, and Alpha 3 lineages (Fig. 2b). The Alpha 1a to 1c sequences from both species are clustered (≥92% sequence similarity), whereas the I. makropetalos Alpha 2 and I. leukodermatus Alpha 3 sequences are phylogenetically separate from the Alpha 1 sequences (≥88.3% sequence similarity) and from each other.

Within the Alpha 1 group, I. leukodermatus contained sequences belonging to all three subgroups (Alpha 1a, 1b, and 1c), while in I. makropetalos only sequences belonging to the Alpha 1b group were found. The Alpha 1a sequence from I. leukodermatus is most closely related to the Alpha 1a endosymbionts of the gutless oligochaete Olavius loisae from the Australian Great Barrier Reef in all treeing methods (≥97.8% sequence similarity). The closest relatives of these sequences are the Alpha 1c sequence from I. leukodermatus (≥92% sequence similarity) and a cluster of clone sequences isolated from Atlantic Ocean and Arctic Ocean bacterioplankton communities and from Pacific Ocean deep-sea sediments (≥97.2% sequence similarity). The Alpha 1b I. leukodermatus and I. makropetalos sequences consistently grouped together in all three phylogenetic analyses (95.5% sequence similarity). These two sequences are distantly related to the Alpha 1a clade.

In addition to the Alpha 1 sequences, both host species contained another alphaproteobacterial sequence, called Alpha 2 in I. makropetalos and Alpha 3 in I. leukodermatus (Fig. 2b). The phylogeny of the Alpha 2 sequence from I. makropetalos was similar in maximum-likelihood and parsimony analyses, with a clone sequence from an Aegean Sea bacterioplankton community being the closest relative (91.1% sequence similarity). These two sequences are most closely related to the free-living halophilic bacteria Rhodovibrio salinarum and Rhodovibrio sodomensis (≥89.3% sequence similarity) and bacteria isolated from a dinoflagellate and coral mucus (≥90.2% sequence similarity). In distance analyses the I. makropetalos Alpha 2 sequence does not cluster with the Aegean Sea clone sequence, but otherwise the relationships to the next relatives are similar to those shown in Fig. 2b.

In I. leukodermatus the relationship of the Alpha 3 sequence to other 16S rRNA sequences varied with the phylogeny method used. In maximum-likelihood analyses the closest relatives were the nitrogen-fixing symbionts Sinorhizobium meliloti and Sinorhizobium fredii (≥95.7% sequence similarity), while parsimony and distance analyses placed the Alpha 3 sequence within a clade of bacteria belonging to the genera Mesorhizobium, Aminobacter, Hoeflea, and Agrobacterium and symbionts from the marine brittle star Ophiactis balli and the ant Tetraponera binghami.

cbbL and aprA genes.
The cbbL gene of RubisCO form I was amplified from all six I. leukodermatus and three I. makropetalos individuals. The sequences from each host species were identical between individuals. The cbbM gene of form II RubisCO was not detectable in either species despite multiple PCR assays under various conditions. Comparative phylogenetic analyses showed that the I. leukodermatus and I. makropetalos cbbL sequences are closely related to each other (99.5% amino acid sequence identity) (Fig. 3a). The closest relatives of these sequences are RubisCO form IA sequences from gammaproteobacterial sulfur-oxidizing endosymbionts of the snail Alvinoconcha hessleri, the clam Solemya velum, and the free-living sulfur-oxidizing bacterium Allochromatium vinosum (cbbL-2) (≥89% amino acid sequence identity). A. vinosum has two cbbL gene copies, of which only cbbL-2 is expressed (42).


Figure 3
View larger version (38K):
[in this window]
[in a new window]
 
FIG. 3. Phylogenetic placement of symbiotic RubisCO large-subunit (a) and APS reductase alpha-subunit (b) sequences from I. leukodermatus and I. makropetalos based on maximum-likelihood analyses. The bars represent 10% estimated sequence divergence, except for the estimated divergence between RubisCO form I and II sequences, which corresponds to 124%. Abbreviations: SRB, sulfate-reducing bacteria; SOX, sulfur-oxidizing bacteria; SRA, sulfate-reducing archaeon. Symbionts of chemosynthetic hosts are listed in bold type.

The aprA gene, coding for the alpha subunit of APS reductase, was found in all six I. leukodermatus and three I. makropetalos specimens and showed no sequence variation within each host species. The I. leukodermatus and I. makropetalos aprA sequences are closely related to each other (99.2% amino acid sequence identity) (Fig. 3b). In the absence of aprA sequences from other symbiotic sulfur-oxidizing bacteria, the closest relatives to the two Inanidrilus sequences are aprA sequences from free-living sulfur-oxidizing bacteria, such as the betaproteobacterium Thiobacillus denitrificans (86.6% amino acid sequence identity) and Chlorobium tepidum of the green sulfur bacteria phylum (60% amino acid sequence identity).

In situ identification.
FISH with oligonucleotide probes confirmed that four of the five 16S rRNA phylotypes isolated from I. leukodermatus and all three phylotypes isolated from I. makropetalos originated from bacteria in the symbiont-containing region between the cuticle and the epidermis of the worm (Fig. 4a and c). In I. leukodermatus, the gammaproteobacterial symbiont (Gamma 1) and three alphaproteobacterial symbionts (Alpha 1a, Alpha 1b, and Alpha 3) co-occurred in all six specimens. The Alpha 1c 16S rRNA sequence, found in a single I. leukodermatus specimen, was not observed in any of the six worms examined despite the design of six specific probes and multiple in situ hybridizations under different conditions. This suggests that the Alpha 1c sequence originated from a contaminant or from a symbiont that is very rare. In I. makropetalos, all three bacterial phylotypes, the gammaproteobacterial symbiont (Gamma 1) and the two alphaproteobacterial symbionts (Alpha 1b and Alpha 2), were observed to coexist in all examined host specimens.


Figure 4
View larger version (68K):
[in this window]
[in a new window]
 
FIG. 4. In situ identification of bacterial symbionts in I. leukodermatus (a and b) and I. makropetalos (c and d). Epifluorescence images show a cross section through the entire worm and the symbiont-containing region of the worm's body wall (all scale bars, 10 µm). (a) Dual hybridization with InaGAM and IleuAlpha1a-IleuAlpha1b-IleuAlpha3 probes, showing gammaproteobacterial symbionts in red and alphaproteobacterial symbionts in green. (Inset) Symbionts at a higher magnification. (b) Triple hybridization with Alpha 1a, Alpha 1b, and Alpha3 probes, showing the Alpha 1a symbionts in green, Alpha 1b symbionts in red, and Alpha 3 symbionts in blue. (c) Dual hybridization with InaGAM and ImakAlpha 1b-ImakAlpha2 probes, showing gammaproteobacterial symbionts in red and alphaproteobacterial symbionts in green. (inset) Symbionts at a higher magnification. (d) Dual hybridization with ImakAlpha1b and ImakAlpha2 probes, showing the Alpha 1b symbionts in green and Alpha 2 symbionts in red.

The hybridization patterns of the Gamma 1 symbionts in the symbiont-containing regions of I. leukodermatus and I. makropetalos were similar (Fig. 4a and c). In both species, the general probe for Gammaproteobacteria, GAM42a, and the specific probe InaGAM for the Gamma 1 sequences from these worms hybridized to large, oval bacteria throughout the symbiont-containing region (Fig. 4a and c), described as the large morphotype in I. leukodermatus (9, 23). This morphotype also hybridized with the specific probe IleuGAM for the Gamma 1 symbionts in I. leukodermatus (9, 23).

The hybridization signal of the general probe for Alphaproteobacteria, ALF968, was limited to small bacterial cells distributed throughout the entire symbiont region of I. leukodermatus and I. makropetalos. The hybridization patterns of the probes specific to the Alpha symbionts of I. leukodermatus, Alpha 1a, 1b, and 3 (Fig. 4b), and of I. makropetalos, Alpha 1b and 2 (Fig. 4d), corresponded to those from the general alphaproteobacterial probe ALF968, indicating that all alphaproteobacterial symbionts in the worms were identified with the specific probes.

While the distribution of the Gamma 1 symbionts in the symbiont-containing region was similar in the two host species, the distribution of the Alpha symbionts was different. The two Alpha symbionts of I. makropetalos, Alpha 1b and 2, were observed in all cross sections of the three individuals. In contrast, the distribution of three Alpha symbionts of I. leukodermatus, Alpha 1a, 1b, and 3, varied between individuals and within an individual, with some cross sections showing the coexistence of all three symbionts and some showing only one or two of the Alpha symbionts. These differences appeared to be random, with no obvious patterns of coexistence between the three Alpha symbionts.


arrow
DISCUSSION
 
Like all other gutless oligochaetes examined to date, I. leukodermatus and I. makropetalos harbor large, oval bacteria identified as Gamma 1 symbionts (Fig. 2a and 4). These symbionts are more closely related to each other than to the Gamma 1 symbionts from hosts belonging to the genus Olavius. The close phylogenetic relationship of the Inanidrilus Gamma 1 symbionts is congruent with the close relationship of host species belonging to the genus Inanidrilus (10) and could therefore have been caused by cospeciation. Alternatively, their biogeography may have played a role, as the two host species are geographically much closer to each other (both from the northwest Atlantic Ocean) than to Olavius species from the Mediterranean Sea and the Pacific Ocean. Additional studies on Gamma 1 symbionts from both Inanidrilus and Olavius host species to better understand the roles of biogeography and cospeciation in these associations are in progress.

The Gamma 1 symbionts of gutless oligochaetes have been identified as chemoautotrophic sulfur oxidizers based on their close phylogenetic relationship to free-living sulfur oxidizers, the presence of sulfur in their cells, and immunohistochemical studies showing the presence of form I RubisCO in these symbionts (10, 25). This study provides further evidence that these symbionts are thiotrophic, based on the presence of the cbbL gene of form I RubisCO and the aprA gene of APS reductase. The close relationship of the I. leukodermatus and I. makropetalos cbbL and aprA genes to those of free-living and symbiotic sulfur oxidizers indicates that the oligochaete sequences originated from their sulfur-oxidizing Gamma 1 symbionts. An alternative explanation, that some or all of the alphaproteobacterial symbionts of these hosts are sulfur-oxidizing chemoautotrophs and that they are the source of these genes, is unlikely. The oligochaete Alpha symbionts are not related to known alphaproteobacterial sulfur oxidizers, and it has been shown that these small bacteria in I. leukodermatus do not contain sulfur and are not labeled by a form I RubisCO antiserum (25).

Prior to this study, alphaproteobacterial symbionts had been found in only a single gutless oligochaete species: O. loisae, from the Australian Great Barrier Reef (7). The presence of these symbionts in two species of the host genus Inanidrilus that are geographically very distant from O. loisae indicates that symbioses with alphaproteobacterial symbionts may be more widespread in gutless oligochaetes than previously assumed. The Alpha 1 symbionts occur in all three host species with alphaproteobacterial symbionts and are relatively closely related to each other (≥92.6% sequence similarity). Within this clade, different phylotypes can co-occur in the same species, such as Alpha 1a and 1b symbionts in O. loisae and I. leukodermatus. In addition to the Alpha 1 lineage of symbionts, I. leukodermatus and I. makropetalos harbor additional alphaproteobacterial symbionts, the former Alpha 2 and the latter Alpha 3 symbionts, that are phylogenetically distinct from each other and the Alpha 1 symbionts.

The metabolism of the alphaproteobacterial symbionts is not currently known. The closest free-living relatives of the Alpha 1 symbionts from both host species and the Alpha 2 symbionts from I. makropetalos are the halophilic bacteria Rhodovibrio salinarum and R. sodomensis (≥89.2% sequence similarity). Both species are photoheterotrophic under anoxic conditions, while only R. salinarum can also grow in the dark under aerobic conditions as a chemoheterotroph (26, 30). Tests for photoautotrophic growth with reduced sulfur compounds were negative for both species (26, 30). The common feature of these bacteria is their use of fermentation products such as lactate, acetate, succinate, malate, and pyruvate as electron and carbon sources during anaerobic growth (26, 30). These metabolites are produced by marine invertebrates when oxygen concentrations become limiting and are excreted in animals without symbionts. It was suggested that in gutless oligochaetes with sulfate-reducing deltaproteobacterial symbionts, these might take up the anaerobic waste products of their hosts, thus recycling these valuable carbon compounds (10, 34). It is intriguing that in all gutless oligochaetes, either alpha- or deltaproteobacterial symbionts coexist with the Gamma 1 symbionts, and it is tempting to speculate that the alphaproteobacterial symbionts might play a role similar to that of the deltaproteobacterial symbionts by recycling the anaerobic waste products of the worms.

The Alpha 3 symbiont from I. leukodermatus is most closely related to nitrogen-fixing Sinorhizobium symbionts of leguminous plants (≥95.7% sequence similarity), suggesting that the Alpha 3 symbionts might also fix N2. However, numerous attempts to amplify the nifH gene, involved in bacterial N2 fixation, were unsuccessful (N. Dubilier and J. Zehr, unpublished data). Since this symbiont occurs only in I. leukodermatus, it does not appear to be essential for nitrogen uptake, as other oligochaete hosts are clearly able to acquire nitrogen without this Alpha 3 symbiont.

The association between the multiple symbiotic bacteria and the gutless oligochaetes I. leukodermatus and I. makropetalos is highly specific and stable within each host species. FISH analyses showed that the symbionts occur regularly in all examined individuals. In addition, the 16S rRNA sequences of the Gamma 1 symbiont from I. leukodermatus worms collected in 1992 (9) and in 1998 (this study) are identical (with the exception of three nucleotide substitutions caused by sequencing error), indicating a high specificity and evolutionary stability of this symbiont over this period of time. A high phylogenetic diversity of symbiotic bacteria was first observed in the gutless oligochaete O. crassitunicatus that harbors up to six bacterial phylotypes (3). This study shows that the stable coexistence of multiple endosymbionts is not limited to a single species but rather appears to be a common feature in oligochaete symbioses.


arrow
ACKNOWLEDGMENTS
 
We are grateful to Jens Krieger for collecting worms and to Christina Zaluski and Silke Wetzel for technical assistance. In Bermuda, Wolfgang Sterrer, Director of the Bermuda Aquarium, Museum, and Zoo, graciously provided logistical support for collecting worms, and the Bermuda Biological Station of Research provided laboratory space. We thank the Caribbean Marine Research Center in the Bahamas for support and lodging during worm collection.

The Max Planck Society, Munich, Germany, provided financial support.

This study is contribution no. 107 of the Bermuda Biodiversity Project, Bermuda Aquarium, Museum, and Zoo.


arrow
FOOTNOTES
 
* Corresponding author. Mailing address: Max Planck Institute of Marine Microbiology, Celsiusstrasse 1, D-28359 Bremen, Germany. Phone: 49-421 2028-932. Fax: 49-421 2028-580. E-mail: ndubilie{at}mpi-bremen.de. Back

{dagger} Present address: Bremen Institute for Materials Testing, Paul-Feller-Str. 1, D-28199 Bremen, Germany. Back


arrow
REFERENCES
 
    1
  1. Altschul, S. F., T. L. Madden, A. A. Schäffer, J. Zhang, Z. Zhang, W. Miller, and D. J. Lipman. 1997. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25:3389-3402.[Abstract/Free Full Text]
  2. 2
  3. Amann, R. I., B. J. Binder, R. J. Olson, S. W. Chisholm, R. Devereux, and D. A. Stahl. 1990. Combination of 16S rRNA-targeted oligonucleotide probes with flow cytometry for analyzing mixed microbial populations. Appl. Environ. Microbiol. 56:1919-1925.[Abstract/Free Full Text]
  4. 3
  5. Blazejak, A., C. Erséus, R. Amann, and N. Dubilier. 2005. Coexistence of bacterial sulfide-oxidizers, sulfate-reducers, and spirochetes in a gutless worm (Oligochaeta) from the Peru margin. Appl. Environ. Microbiol. 71:1553-1561.[Abstract/Free Full Text]
  6. 4
  7. Bright, M., and O. Giere. 2005. Microbial symbiosis in Annelida. Symbiosis 38:1-45.[Medline]
  8. 5
  9. Cavanaugh, C. M., and J. J. Robinson. 1996. CO2 fixation in chemoautotrophic-invertebrate symbioses: expression of form I and form II RubisCO. In M. E. Lidstrom and F. R. Tabita (ed.), Microbial. growth on C1 compounds: Proceedings of the 8th International Symposium on Microbial growth on C1 Compounds. Kluwer Academic Publishers, Dordrecht, The Netherlands.
  10. 6
  11. Cole, J. R., B. Chai, T. L. Marsh, R. J. Farris, Q. Wang, S. A. Kulam, S. Chandra, D. M. McGarrell, T. M. Schmidt, G. M. Garrity, and J. M. Tiedje. 2003. The Ribosomal Database Project (RDP-II): previewing a new autoaligner that allows regular updates and the new prokaryotic taxonomy. Nucleic Acids Res. 31:442-443.[Abstract/Free Full Text]
  12. 7
  13. Dubilier, N., R. Amann, C. Erséus, G. Muyzer, S. Y. Park, O. Giere, and C. M. Cavanaugh. 1999. Phylogenetic diversity of bacterial endosymbionts in the gutless marine oligochaete Olavius loisae (Annelida). Mar. Ecol. Prog. Ser. 178:271-280.[CrossRef]
  14. 8
  15. Dubilier, N., A. Blazejak, and C. Rühland. 2006. Symbioses between bacteria and gutless marine oligochaetes. In J. Overmann (ed.), Molecular basis of symbiosis. Springer-Verlag, New York, N.Y.
  16. 9
  17. Dubilier, N., O. Giere, D. L. Distel, and C. M. Cavanaugh. 1995. Characterization of chemoautotrophic bacterial symbionts in a gutless marine worm (Oligochaeta, Annelida) by phylogenetic 16S rRNA sequence analysis and in situ hybridization. Appl. Environ. Microbiol. 61:2346-2350.[Abstract]
  18. 10
  19. Dubilier, N., C. Mülders, T. Ferdelman, D. de Beer, A. Pernthaler, M. Klein, M. Wagner, C. Erséus, F. Thiermann, J. Krieger, O. Giere, and R. Amann. 2001. Endosymbiotic sulphate-reducing and sulphide-oxidizing bacteria in an oligochaete worm. Nature 411:298-302.[CrossRef][Medline]
  20. 11
  21. Erséus, C. 1992. A generic revision of the Phallodrilinae (Oligochaeta, Tubificidae). Zool. Scr. 21:5-48.
  22. 12
  23. Erséus, C. 2003. The gutless Tubificidae (Annelida, Oligochaeta) of the Bahamas. Meiofauna Mar. 12:59-84.
  24. 13
  25. Erséus, C. 1990. The marine Tubificidae (Oligochaeta) of the barrier reef ecosystems at Carrie Bow Cay, Belize, and other parts of the Caribbean Sea, with descriptions of twenty-seven new species and revision of Heterodrilus, Thalassodrilides and Smithsonidrilus. Zool. Scr. 19:243-303.[CrossRef]
  26. 14
  27. Erséus, C. 1984. Taxonomy and phylogeny of the gutless Phallodrilinae (Oligochaeta, Tubificidae), with descriptions of one new genus and twenty-two new species. Zool. Scr. 13:239-272.
  28. 15
  29. Erséus, C. 1991. Two new deep-water species of the gutless genus Olavius (Oligochaeta: Tubificidae) from both sides of North America. Proc. Biol. Soc. Wash. 104:627-630.
  30. 16
  31. Felbeck, H., G. Liebezeit, R. Dawson, and O. Giere. 1983. CO2 fixation in tissues of marine oligochaetes (Phallodrilus leukodermatus and P. planus) containing symbiotic, chemoautotrophic bacteria. Mar. Biol. 75:187-191.[CrossRef]
  32. 17
  33. Finogenova, N. 1986. Six new species of marine Tubificidae (Oligochaeta) from the continental shelf off Peru. Zool. Scr. 15:45-51.
  34. 18
  35. Friedrich, M. W. 2002. Phylogenetic analysis reveals multiple lateral transfers of adenosine-5'-phosphosulfate reductase genes among sulfate-reducing microorganisms. J. Bacteriol. 184:278-289.[Abstract/Free Full Text]
  36. 19
  37. Giere, O. 1981. The gutless marine oligochaete Phallodrilus leukodermatus. Structural studies on an aberrant tubificid associated with bacteria. Mar. Ecol. Prog. Ser. 5:353-357.[CrossRef]
  38. 20
  39. Giere, O. 1979. Studies on marine Oligochaeta from Bermuda, with emphasis on new Phallodrilus species (Tubificidae). Cah. Biol. Mar. 20:301-314.
  40. 21
  41. Giere, O., and C. Erséus. 2002. Taxonomy and new bacterial symbioses of gutless marine Tubificidae (Annelida, Oligochaeta) from the Island of Elba (Italy). Org. Divers. Evol. 2:289-297.[CrossRef]
  42. 22
  43. Giere, O., and J. Krieger. 2001. A triple bacterial endosymbiosis in a gutless oligochaete (Annelida): ultrastructural and immunocytochemical evidence. Invertebr. Biol. 120:41-49.
  44. 23
  45. Giere, O., C. Nieser, R. Windoffer, and C. Erséus. 1995. A comparative structural study on bacterial symbioses of Caribbean gutless Tubificidae (Annelida, Oligochaeta). Acta Zool. 76:281-290.
  46. 24
  47. Hipp, W. M., A. S. Pott, N. Thum-Schmitz, I. Faath, C. Dahl, and H. G. Trüper. 1997. Towards the phylogeny of APS reductases and sirohaem sulfite reductases in sulfate-reducing and sulfur-oxidizing prokaryotes. Microbiology 143:2891-2902.[Abstract/Free Full Text]
  48. 25
  49. Krieger, J., O. Giere, and N. Dubilier. 2000. Localization of RubisCO and sulfur in endosymbiotic bacteria of the gutless marine oligochaete Inanidrilus leukodermatus (Annelida). Mar. Biol. 137:239-244.[CrossRef]
  50. 26
  51. Mack, E. E., L. Mandelco, C. R. Woese, and M. T. Madigan. 1993. Rhodospirillum sodomense, sp. nov., a Dead Sea Rhodospirillum species. Arch. Microbiol. 160:363-371.
  52. 27
  53. Manz, W., R. Amann, W. Ludwig, M. Wagner, and K.-H. Schleifer. 1992. Phylogenetic oligodeoxynucleotide probes for the major subclasses of Proteobacteria: problems and solutions. Syst. Appl. Microbiol. 15:593-600.
  54. 28
  55. Muyzer, G., A. Teske, C. O. Wirsen, and H. W. Jannasch. 1995. Phylogenetic relationships of Thiomicrospira species and their identification in deep-sea hydrothermal vent samples by denaturing gradient gel electrophoresis of 16S rDNA fragments. Arch. Microbiol. 164:165-172.[CrossRef][Medline]
  56. 29
  57. Neef, A. 1997. Anwendung der in situ-Einzelzell-Identifizierung von Bakterien zur Populationsanalyse in komplexen mikrobiellen Biozönoses. Ph.D. dissertation. Technical University Munich, Munich, Germany.
  58. 30
  59. Nissen, H., and I. D. Dundas. 1984. Rhodospirillum salinarum sp. nov., a halophilic photosynthetic bacterium isolated from a Portuguese saltern. Arch. Microbiol. 138:251-256.[CrossRef]
  60. 31
  61. Nylander, J. A. A., C. Erséus, and M. Källersjö. 1999. A test of monophyly of the gutless Phallodrilinae (Oligochaeta, Tubificidae) and the use of a 573 bp region of the mitochondrial cytochrome oxidase I gene in analysis of annelid phylogeny. Zoo. Scr. 28:305-313.
  62. 32
  63. Olsen, G. J., R. Overbeck, N. Larson, T. L. Marsh, M. J. McCaughey, M. A. Maciukenas, W. M. Kuan, T. J. Macke, Y. Xing, and C. R. Woese. 1992. Ribosomal Database Project. Nucleic Acids Res. 20:2199-2200.
  64. 33
  65. Robinson, J. J., and C. M. Cavanaugh. 1995. Expression of form I and form II RubisCO in chemoautotrophic symbioses: implications for the interpretation of stable carbon isotope values. Limnol. Oceanogr. 40:1496-1502.
  66. 34
  67. Rühland, C., A. Blazejak, A. Loy, M. Wagner, C. Erséus, R. Amann, and N. Dubilier. Unpublished data.
  68. 35
  69. Sánchez, O., I. Ferrera, C. Dahl, and, M. J. 2001. In vivo role of APSreductase in the purple sulfur bacterium Allochromatium vinosum. Arch. Microbiol. 176:301-305.[CrossRef][Medline]
  70. 36
  71. Schizas, N. V., G. T. Street, B. C. Coull, G. T. Chandler, and J. M. Quattro. 1997. An efficient DNA extraction method for small metazoans. Mol. Mar. Biol. Biotechnol. 6:381-383.[Medline]
  72. 37
  73. Schönhuber, W., B. Zarda, S. Eix, R. Rippka, M. Herdman, W. Ludwig, and R. Amann. 1999. In situ identification of cyanobacteria with horseradish peroxidase-labeled, rRNA-targeted oligonucleotide probes. Appl. Environ. Microbiol. 65:1259-1267.[Abstract/Free Full Text]
  74. 38
  75. Schwedock, J., T. L. Harmer, K. M. Scott, H. J. Hektor, A. P. Seitz, M. C. Fontana, D. L. Distel, and C. M. Cavanaugh. 2004. Characterization and expression of genes from the RubisCO gene cluster of the chemoautotrophic symbiont of Solemya velum: cbbLSQO. Arch. Microbiol. 182:18-29.[CrossRef][Medline]
  76. 39
  77. Stein, J. L., M. Haygood, and H. Felbeck. 1990. Nucleotide sequence and expression of a deep-sea ribulose-1,5-bisphosphate carboxylase gene cloned from a chemoautotrophic bacterial endosymbiont. Proc. Natl. Acad. Sci. USA 87:8850-8854.[Abstract/Free Full Text]
  78. 40
  79. Tabita, F. R. 1995. The biochemistry and metabolic regulation of carbon metabolism and CO2 fixation in purple bacteria. In R. E. Blankenship, M. T. Madigan, and C. E. Bauer (ed.), Anoxygenic photosynthetic bacteria. Kluwer, Dordrecht, The Netherlands.
  80. 41
  81. Tabita, F. R. 1988. Molecular and cellular regulation of autotrophic carbon dioxide fixation in microorganisms. Microbiol. Rev. 52:155-189.[Free Full Text]
  82. 42
  83. Viale, A. M., H. Kobayashi, and T. Akazawa. 1989. Expressed genes for plant-type ribulose 1,5-bisphosphate carboxylase/oxygenase in the photosynthetic bacterium Chromatium vinosum, which possesses two complete sets of the genes. J. Bacteriol. 171:2391-2400.[Abstract/Free Full Text]
  84. 43
  85. Wallner, G., R. Amann, and W. Beisker. 1993. Optimizing fluorescent in situ-hybridization with rRNA-targeted oligonucleotide probes for flow cytometric identification of microorganisms. Cytometry 14:136-143.[CrossRef][Medline]


Applied and Environmental Microbiology, August 2006, p. 5527-5536, Vol. 72, No. 8
0099-2240/06/$08.00+0     doi:10.1128/AEM.02441-05
Copyright © 2006, American Society for Microbiology. All Rights Reserved.




This article has been cited by other articles:

  • Behrens, S., Losekann, T., Pett-Ridge, J., Weber, P. K., Ng, W.-O., Stevenson, B. S., Hutcheon, I. D., Relman, D. A., Spormann, A. M. (2008). Linking Microbial Phylogeny to Metabolic Activity at the Single-Cell Level by Using Enhanced Element Labeling-Catalyzed Reporter Deposition Fluorescence In Situ Hybridization (EL-FISH) and NanoSIMS. Appl. Environ. Microbiol. 74: 3143-3150 [Abstract] [Full Text]  
  • Meyer, B., Kuever, J. (2007). Molecular analysis of the distribution and phylogeny of dissimilatory adenosine-5'-phosphosulfate reductase-encoding genes (aprBA) among sulfur-oxidizing prokaryotes. Microbiology 153: 3478-3498 [Abstract] [Full Text]  

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Blazejak, A.
Right arrow Articles by Dubilier, N.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Blazejak, A.
Right arrow Articles by Dubilier, N.
Agricola
Right arrow Articles by Blazejak, A.
Right arrow Articles by Dubilier, N.