Comparative metagenomics of Daphnia symbionts
© Qi et al; licensee BioMed Central Ltd. 2009
Received: 04 March 2008
Accepted: 21 April 2009
Published: 21 April 2009
Shotgun sequences of DNA extracts from whole organisms allow a comprehensive assessment of possible symbionts. The current project makes use of four shotgun datasets from three species of the planktonic freshwater crustaceans Daphnia: one dataset from clones of D. pulex and D. pulicaria and two datasets from one clone of D. magna. We analyzed these datasets with three aims: First, we search for bacterial symbionts, which are present in all three species. Second, we search for evidence for Cyanobacteria and plastids, which had been suggested to occur as symbionts in a related Daphnia species. Third, we compare the metacommunities revealed by two different 454 pyrosequencing methods (GS 20 and GS FLX).
In all datasets we found evidence for a large number of bacteria belonging to diverse taxa. The vast majority of these were Proteobacteria. Of those, most sequences were assigned to different genera of the Betaproteobacteria family Comamonadaceae. Other taxa represented in all datasets included the genera Flavobacterium, Rhodobacter, Chromobacterium, Methylibium, Bordetella, Burkholderia and Cupriavidus. A few taxa matched sequences only from the D. pulex and the D. pulicaria datasets: Aeromonas, Pseudomonas and Delftia. Taxa with many hits specific to a single dataset were rare. For most of the identified taxa earlier studies reported the finding of related taxa in aquatic environmental samples. We found no clear evidence for the presence of symbiotic Cyanobacteria or plastids. The apparent similarity of the symbiont communities of the three Daphnia species breaks down on a species and strain level. Communities have a similar composition at a higher taxonomic level, but the actual sequences found are divergent. The two Daphnia magna datasets obtained from two different pyrosequencing platforms revealed rather similar results.
Three clones from three species of the genus Daphnia were found to harbor a rich community of symbionts. These communities are similar at the genus and higher taxonomic level, but are composed of different species. The similarity of these three symbiont communities hints that some of these associations may be stable in the long-term.
Metagenomics is the field that infers the properties of a habitat through the analysis of genomic sequence information obtained from a sample usually collected from a single habitat. The sequences are usually compared to databases, with the aim to characterize the biological community of this habitat. Among the advantages of this explorative method are the free and uncomplicated sampling of the material, the possibility of obtaining sequences from unknown and unculturable organisms, the absence of any taxonomic restrictions and the relative ease of conducting such studies [1–4]. Metagenomics studies have been done in various habitats, including sea water , ice cores  and deep mine communities . Of particular recent interest has been the application of metagenomic approaches to study samples obtained from organisms, which harbor various symbionts, such as unknown and uncultuable bacteria, protozoa or viruses. For example, the symbiont communities of honey bees , the guts of mice  and humans , marine sponges , oligochaetes  and plant-rhizobacteria  revealed many new symbiont taxa. However, not only samples collected with the aim to find symbionts revealed previously unknown organisms, but also datasets from genome projects where one single genome was targeted may contain sequences of other species, presumably symbionts . Here we report on the bacterial communities associated with three clones each from one species of crustaceans of the genus Daphnia, which had been used in genome projects and revealed besides sequences to the targeted species, a rich body of sequences to other species. We use the term symbiont to include organisms that were found to be associated with the samples of these Daphnia, disregarding whether they are parasites, commensals or mutualists. We cannot rule out, that some of these organisms are independent of the Daphnia, e.g. free living bacteria in the water, parts of the ingested food or contaminants from handling the samples. For simplicity we use the term symbiont throughout this article.
Here we take advantage of shotgun sequences obtained from three laboratory clones (= iso-female lines) each from one Daphnia species to search for indications of bacterial and plastid symbionts. For this we compared the sequences against the NCBI-nt database on nucleotide sequences using BLASTN  and analyzed and ordered the results using the metagenomics software MEGAN . This software allows the exploration of the taxonomic content of a community sample based on the NCBI taxonomy. Community shotgun datasets represent sequences independently sampled from random regions of genomes randomly selected from a given community. These sequences can have very different levels of conservation. Without any assumptions about the functions of the sequences used, MEGAN associates each sequence to the lowest common ancestor of the set of taxa it hits. Thus, species specific sequences are assigned to low order taxa such as species or strains, while widely conserved sequences are assigned to high-order taxa. In other words, the taxonomical level of the assigned taxon reflects the level of conservation of the sequence. The strength of this statistical approach is that it makes use of all kind of sequences for taxon identification. Therefore, when using random sequences MEGAN, will usually show better taxonomic resolution than an analysis using only a small set of phylogenetic markers . This type of analysis is in particular useful when, as is the case here, datasets are analyzed, which were obtained by random shotgun sequencing, rather than targeted sequencing (see also ) and where the length of the sequence reads are short [20, 22].
Our choice to use the software MEGAN for the analysis of the datasets from the Daphnia projects is based on several aspects, which help to reduce known problems in comparative metagenomics. A known shortcoming of the assignment of sequences to taxonomic groups is its inability to deal with horizontally transferred genes and the inability of mapping sequences to internal nodes of the tree . However, these problems are mainly of concern when using "best-BLAST-hit" mapping. The software MEGAN was developed to avoid this problem (see previous paragraph). A further problem of assigning sequences to taxonomic groups is the well know bias in the taxon representation in our databases [24, 25]. This problem cannot be fully solved, but the ability of MEGAN to assign sequence to the lowest common ancestor, ameliorates the consequences of a database bias. Sequences will be assigned to the common ancestor of the true species in question and those being represented in the database. Novel sequences will not be assigned at all .
The aims of our analysis were first to compare the shotgun sequences of the prokaryote communities coming from three Daphnia species. Second to test if the shotgun sequences give evidence for a plastid symbiont in Daphnia as had been suggested . Third, to estimate the repeatability of a metagenomics approach using two different sequencing platforms, the pyrosequencers GS 20 and GS FLX  for one of the three Daphnia species.
Results and discussion
Number of sequences assigned and unassigned in the MEGAN analysis.
Assigned to cellular organisms
Assigned to Bacteria without Firmicutes1
Sequences without hits
D. magna GS 20
D. magna GS FLX
Summary of the four datasets included in this analysis.
D. magna GS 20
D. magna GS FLX
Original input data:
Possible bacterial scaffolds
Contigs and raw reads longer than 500 bps
Contigs longer than 100 bps
Contigs longer than 100 bps
No. of original sequences
Total length (bps)
Average length (mean ± stdev bps)
2,743 ± 7,205
987 ± 255
988 ± 2,830
919 ± 2,507
Median length (bps)
Minimum length (bps)
Maximum length (bps)
Sequence fragments subjected to BLASTN:
Total length (bps)
Average length (mean ± stdev bps)
521 ± 100
524 ± 195
459 ± 149
463 ± 131
Median length (bps)
Assignment of sequences to the bacteria, without Firmicutes and Cyanobacteria
Another genus of the Bacteroidetes, which was consistently found in all datasets is Cytophaga (Fig. 3) These are gliding bacteria found in freshwater and marine habitats, in soil and in decomposing organic matter. However, hits to this genus were never frequent (between 10 and 25 hits).
Taxa within the Proteobacteria, which attracted at least 1% of the sequences within at least one of the four datasets.
D. magna GS 20
D. magna GS FLX
A few other genera within the Betaproteobacteria attracted relatively high numbers of sequences across all or most of the datasets: Chromobacterium, Methylibium, Bordetella, Burkholderia and Cupriavidus (Table 3, Fig. 5). Of those Methylibium petroleiphilum was highly represented. However, a closer inspection of the sequence alignments indicates that the species in our datasets is not exactly this, but a related species.
Hits to species of the genus Aeromonas were found in large number in the D. pulicaria dataset, but hardly in the other sets (Table 3, Fig. 6). Hits were mainly to A. hydrophila and A. salmonicida, but similarities were below 100%. Both can live under aerobic or anaerobic conditions and are found in water. A. hydrophila is an opportunistic pathogen of humans, A. salmonicida causes the fish disease, furunculosis.
The single most often assigned genus in the entire analysis was Pseudomonas in the D. pulex dataset (10,994 assigned reads, 43.3%). These hits were mainly to the species P. fluorescens (7,067 reads), and in particularly to the strain PfO-1. Similar, but not as extreme was the presence of the same bacterium in the D. pulicaria sequences (Table 3, Fig. 6). The P. fluorescens PfO-1 genome project was run in the same genome center (The DOE Joint Genome Institute (JGI, http://www.jgi.doe.gov/) where the D. pulex and the D. pulicaria sequences were obtained and it seemed possible, that these hits reflect a contamination in the D. pulex scaffolds, rather than a symbiont of D. pulex. However, inspection of bit scores and sequence identity values in the BLASTN outputs indicated that the Daphnia symbiont is clearly not P. fluorescens PfO-1. The P. fluorescens group includes diverse bacteria that are found in soil, but also in aquatic environments.
A further contamination candidate is the Gammmaproteobacterium Serratia, to which we found 2,184 matched sequences in the D. pulex genome. However, it is hardly seen among the D. pulicaria sequences, and not seen at all among the D. magna sequences (Table 3, Fig. 6). The species to which most sequences were assigned is Serratia proteamaculans 568, whose genome was sequenced as well by the DOE Joint Genome Institute. Also here, the inspection of the BLASTN results indicated high similarity, but few perfect matches, excluding contamination at the JGI. Serratia are often associated with the human gut, but are not pathogenic.
Another genus with many hits to the D. pulex and the D. pulicaria sequences, but not to the D. magna sequences (Table 3), is the already mentioned Betaproteobacterium Delftia (Fig. 5). The DOE Joint Genome Institute sequenced Delftia acidovorans strain SPH-1, which is the strain most of the sequences were assigned to. However, inspection of the BLASTN results again showed that the Daphnia symbiont is clearly not D. acidovorans strain SPH-1.
Searching for Cyanobacteria and plastid sequences
The D. pulicaria dataset revealed 23 sequences assigned to plastids. One of them was a short sequence (100 bps) to the chloroplasts of the green algae Chlamydomonas, the other to the chloroplasts of flowering plants. Hits to the later came mostly from one scaffold and had high bit scores (> 500) and similarities of more than 90%. The D. pulex sequences revealed no hits to plastids, but this is not surprising, as the dataset had been sorted out to contain predominately prokaryote sequences. The D. magna GS 20 dataset did not reveal any hits to plastids. The D. magna GS FLX sequences contained a short sequence (104 bps) matched to a plastid, the chloroplast of the green algae Stigeoclonium helveticum.
The presence of plastid sequences in Daphnia shotgun datasets has however, to be looked at with care, as unicellular green algae are the main food of Daphnia, both in the field and in the laboratory [34, 35]. However, the few sequences assigned to plastids here seem not to correspond closely with the algae, which were used to feed the Daphnia in the cultures, before they were used for DNA extraction. The D. magna and the D. pulex clone had been kept on an exclusive diet of the green algae Scenedesmus sp. and the D. pulicaria clone on a diet of the green algae Ankistrodesmus falcatus.
All in all we consider this as rather weak evidence for plastid symbionts in these Daphnia samples. The original finding was done in D. obtusa , which was not included in our study. The authors had observed variation in the type and frequency of plastid occurrence in this species, so it may not be surprising that things are different in other species. Furthermore, the long maintenance of the Daphnia clones in laboratory cultures may have contributed to a loss of plastids. Therefore, the absence of evidence from our metagenomics analysis is certainly not evidence for the absence of possible plastid symbionts in Daphnia.
Searching for 16S rDNA sequences
16S rDNA sequences close to full length identified in the four datasets.
Best matched 16S
Description of the next three matches4
D. magna GS 20
Daphnia endosymbiotic bacterium2
uncultured Pasteuria sp., P. nishizawae, P. penetrans
aquatic bacterium R1-C1
uncultured Cytophagales bacterium, aquatic bacterium R1-C5, uncultured bacterium
D. magna GS FLX
Escherichia coli W3110, E. coli K12, E. coli
uncultured Cytophagales bacterium
uncultured bacterium, Flavobacterium sp. Nj-26, uncultured Flavobacteriales bacterium
Myxococcales str. NOSO-1, Chondromyces pediculatus, Polyangium thaxteri
Daphnia endosymbiotic bacterium2
uncultured Pasteuria sp., P. nishizawae, P. penetrans
Flavobacterium sp. MH45
Arctic sea ice bacterium ARK10164, uncultured bacterium, Flavobacterium succinicans
Comamonadaceae bacterium BP-1b,
uncultured Burkholderiales bacterium, Comamonadaceae bacterium BP-1b, uncultured proteobacterium
uncultured Burkholderiales bacterium
uncultured bacterium, Comamonadaceae bacterium BP-1b, Comamonadaceae bacterium BP-1
Flavobacterium sp. GOBB3-209
uncultured bacterium, uncultured Cytophagales bacterium, uncultured Sphingobacteriales bacterium
uncultured Burkholderiales bacterium
uncultured beta proteobacterium, uncultured organism, Rhodoferax ferrireducens T118
aquatic bacterium R1-B19
uncultured beta proteobacterium, aquatic bacterium R1-B6, uncultured Burkholderiales bacterium
Aeromonas sp. 'CDC 859-83', A. molluscorum, uncultured bacterium
uncultured bacterium, Modestobacter multiseptatus, Sporichthya polymorpha
Pseudomonas sp. R-25061
Pseudomonas sp. R-25209, uncultured bacterium, P. pseudoalcaligenes
uncultured gamma proteobacterium, uncultured Pseudomonas sp., Pseudomonas sp. G2
Pseudomonas sp. GD100
Pseudomonas sp. Pb1(2006), P. poae, P. lurida
uncultured Burkholderiales bacterium
uncultured bacterium, Variovorax paradoxus, uncultured bacterium SJA-62
uncultured Cytophagales bacterium
uncultured bacterium, uncultured Bacteroidetes bacterium, rhizosphere soil bacterium RSC-II-81
uncultured Cytophagales bacterium, uncultured Bacteroidetes bacterium, uncultured bacterium
Hydrogenophaga sp. AH-24, Hydrogenophaga sp. CL3, Hydrogenophaga sp. YED1-18
P. argentinensis, P. fluorescens PfO-1,
uncultured soil bacterium, uncultured Comamonadaceae bacterium, uncultured beta proteobacterium
Serratia proteamaculans 568
Serratia proteamaculans 568, uncultured bacterium, uncultured proteobacterium
uncultured bacterium, Chitinibacter tainanensis, uncultured proteobacterium
gamma proteobacterium GPTSA100-21
gamma proteobacterium GPTSA100-22, uncultured bacterium, gamma proteobacterium GPTSA100-26
Pseudomonas sp. Hsa.28
uncultured bacterium, uncultured Pseudomonas sp., P. anguilliseptica
aquatic bacterium R1-B19
uncultured beta proteobacterium, aquatic bacterium R1-B6, aquatic bacterium R1-B7
uncultured Pseudomonas sp.
gamma proteobacterium LC-G-2, Pseudomonas sp. 7-1, P. fluorescens
The 16S rDNA sequences identified only a small subset of the species/genus found in our main analysis based on comparison to NCBI-nt database. One likely explanation of this discrepancy is the low sequencing coverage within the 16S rDNA regions in the shotgun datasets. Another explanation could be that some of the earlier predictions were false positives. However, MEGAN associates a sequence to the lowest common ancestor of the set of taxa defined by all matches above defined thresholds. The amount of false predictions is predicted to be low since the algorithm makes higher amount of unspecific assignments to higher taxonomy levels . Certainly when taxa were inferred regardless if the matched sequence was a suitable phylogenetic marker or not, it could not be excluded that some of the predictions were results of horizontal gene transfer events. However, if this were the case, MEGAN would assign the hit to the least common ancestor of the species, which were involved in horizontal gene transfer, unless neither these species nor related species are in the NCBI database. It was predicted that computing taxonomic content based on sequence comparison to NCBI-nt database will show better resolution at all levels of the taxonomy than an analysis based on a small set of phylogenetic markers or on 16S rDNA sequences alone [20, 21]. Our results are consistent with this prediction.
Searching for identical and similar sequences common in four datasets
Although sequences in all datasets were assigned to similar bacterial taxa, it is not clear how similar the sequences are across datasets. To identify common sequences, we compared the D. magna GS 20 sequences with sequences from D. magna GS FLX, D. pulex, and D. pulicaria using BLASTN. Identical or nearly identical sequences were identified when a stretch longer than 80% of a query sequence can be aligned with over 98% nucleotide sequence identity to a hit sequence. With this criterion five D. magna GS 20 contigs (corresponding to six D. pulex scaffolds and 12 D. pulicaria reads) were identified. Hits identical to these sequences were all found in complete genome sequences of Escherichia coli W3110 (AP009048.1) and E. coli K12 MG1655 (U00096.2), which suggests that commensal E. coli strains carried by the three Daphnia species are highly similar.
With a less stringent criterion (a stretch longer than 50% of a query sequence can be aligned with over 90% nucleotide sequence identity to a hit sequence), similar sequences to about 80 GS 20 contig sequences were also identified across the datasets. These sequences mainly fall into taxa within the Proteobacteria, with a few sequences assigned to Flavobacterium.
The small number of similar sequences shared across the datatsets suggested the bacterial community carried by the three Daphnia clones from which our datasets originated might be diverse at species and strain level, despite very high homogeneousness observed at higher taxonomy nodes. It should be noted however, that our datasets do not originate directly from field samples, but from three clones, which had been kept in three different laboratories for several generations before the DNA was isolated. This may possibly influence our results in two ways. First, we cannot truly make statements about three Daphnia species, but only about three clones, each coming from a different Daphnia species. Including more clones, might reveal more bacterial symbionts. Second, while culturing these clones in the laboratory, the symbiont community may have changed both qualitatively and quantitatively. New bacterial species may have arrived with food or culture conditions, while other bacteria may have been lost due to the inappropriateness of the laboratory conditions for their culture. For the current analysis, no attempts have been undertaken to vary the culture conditions for any of the three clones and the bacteria associated with the food alga have not been analyzed.
Repeatability of the metagenomics approach
Using contigs instead of reads
For the D. pulicaria dataset, both contigs and singleton raw reads were included in our analysis. For the other three datasets, we used only sequences, which had previously been assembled to contigs or scaffolds. This reduced the number of sequences and thus the number of BLASTN searches considerably. Using large numbers of raw reads would have been beyond our computing power and the abilities of the MEGAN software within a reasonable time period. Using contigs and scaffolds influences the results in various ways. First, it strongly reduces redundancy in the dataset and therefore makes the analysis much quicker. Second, it compromises somewhat the usefulness of the number of assigned sequences as a measure for the abundance for the different taxa. The number of assigned sequences is still a relative measure for the frequency of a given taxa, but the larger the real number of hits would have been, the more strongly the value is reduced. Third, rare members of the symbiont community are likely to remain undetected, because the few reads sequenced for rare species, were unlikely to be assembled in contigs. Thus, our estimates of the number of taxa detected are likely to underestimate the true number of taxa in the community. This conclusion is also supported by the observation that the D. pulicaria dataset contained the highest number of taxa identified.
Our analysis of shotgun sequences of three clones, each from one Daphnia species revealed a rich bacterial community to be associated with these clones. The particular data structure of our analysis allows for certain conclusions to be drawn. First, the majority of the common bacterial taxa identified are found in all Daphnia datasets. While the D. pulex and D. pulicaria clone cultures from which DNA was isolated originated from laboratories in North America, the D. magna cultures originate from a laboratory in Switzerland. To the best of our knowledge, there was never a cross Atlantic exchange of cultures between laboratories by the time these samples had been taken. Thus, we speculate that the similarity of the symbiont communities in European and North American Daphnia samples, indicates a long lasting stability of these associations.
Second, the symbiont communities across the three Daphnia species are remarkable similar, yet, they are not identical. At sequence level, the similarity breaks down, indicating that each Daphnia species harbors different species or strains of bacterial symbionts.
Third, some bacterial taxa were found to be specific to the two datasets produced in the DOE Joint Genome Institute (JGI). Coincidentally, some of the published genomes in these taxa had been originally sequenced by JGI, leading to speculations of whether the JGI may have contaminated the Daphnia samples. Our analysis allows us clearly to reject this hypothesis. Whether the bacterial taxa found to be associated with specific Daphnia samples are contaminations of the laboratory where they were cultured previous to sequencing, or if they are natural symbionts of the Daphnia, cannot not be worked out here.
Fourth, there is no clear evidence for a stable cyanobacterial or plastid symbiont in the Daphnia species. The few scattered hits to some plastid and Cyanobacteria may have been a contamination with the algae food of the Daphnia. Plastid symbionts had been observed in D. obtusa . However, the long laboratory culture of the clones used in the genome study may have influenced the presence of such a photoactive symbiont.
The D. pulex dataset
The sequences of D. pulex are from the DGC whole genome sequencing project. The chosen D. pulex clone called The Chosen One was cultured at Indiana University, Bloomington, USA on a diet of the green algae Scenedesmus sp. The animals used to isolate the DNA for the genome project were treated with tetracycline (250 mg/L overnight) before DNA isolation to reduce their bacterial load. Sequencing was done at the DOE Joint Genome Institute (JGI) using the Sanger method. These sequences were obtained from the wFleaBase website http://wfleabase.org:7182/genome/Daphnia_pulex/current/genome-assembly-full-jazz_20060901/scaffolds/sequences/. Scaffolds included in this study were excluded scaffolds, prokaryotic scaffolds, and possible bacterial scaffolds in the current D. pulex genome assembly http://wfleabase.org:7182/genome/Daphnia_pulex/current/bacteria/dpulex_jgi060905_possible_bacterial.txt.
The D. pulicaria dataset
Daphnia pulicaria is closely related to D. pulex and forms with intermediate characters are frequently encountered, suggesting hybridization of these two species. Indeed, allozyme test for allelic variation at the lactate dehydrogenase loci show both fast and slow electromorphic alleles, indicating that the chosen D. pulicaria strain is a pulicaria/pulex hybrid. This chosen D. pulicaria clone was cultured at the Hubbard Center for Genome Studies at the University of New Hampshire, USA, on a diet of the green algae Ankistrodesmus falcatus. Previous to it's culturing at the University of New Hampshire it was maintained in a laboratory at Utah State University. The animals used to isolate the DNA for the genome project were treated with tetracycline (250 mg/L overnight) before DNA isolation to reduce their bacterial load. Sequencing of D. pulicaria was also done at the DOE Joint Genome Institute (JGI) using the Sanger method. A low coverage genome assembly of a D. pulicaria clone is available to DGC members, and others may request access to this data. As the DGC and JGI data agreements allow, this will be released for public access on the wfleabase database: http://wfleabase.org/genome/Daphnia_pulicaria/. For more information on the D. pulex and D. pulicaria genome data see http://wfleabase.org/.
The D. magna datasets
The sequences of D. magna originated from a shotgun sequencing project which aimed at sequencing the endoparasitic bacterium P. ramosa. During the analysis of the data large number of sequences clearly unrelated to the Firmicutes (the group to which P. ramosa belongs) showed up. Only these sequences are included in this paper. As these data are not yet published elsewhere, we describe here the DNA isolation, library construction and sequencing in detail.
Daphnia magna cultures were raised at the University of Fribourg, Switzerland on a diet of the green algae Scenedesmus sp. The Daphnia had been exposed to the gram-positive bacterium Pasteuria ramosa, an endo-parasite of Daphnia  when they were 3–5 days old. Most animals became infected and were shipped for further processing to the University of Florida, USA. One thousand P. ramosa infected D. magna were suspended in 5 ml of Buffer A (1.0 M NaCl, 50 mM Tris-HCl pH 8.0) and homogenized gently in a glass pestle and mortar. The homogenate was passed through a 50–100 micron metal mesh and 21 micron nylon mesh to remove Daphnia debris. About 5,000,000 P. ramosa cells were obtained and resuspended in 450 μl of Buffer A. These were added to an equal volume (450 μl) of 2% agarose for preparing a gel plug to embed the vegetative cells, and 10 gel plugs were produced. To disrupt cells gently, the gel plugs were transferred into Buffer B (0.2% sodium deoxycholate, 0.5% Brij 58, 0.5% sarcosine, 50 mM Tris-HCl pH 8.0, 100 mM EDTA pH 8.0, 0.40 M NaCl) and incubated at 37°C overnight. These were then transferred into 10 ml of Buffer C (100 mM NaCl, 50 mM Tris-HCl pH 8.0, 100 mM EDTA pH 8.0, 0.5% sarcosine, 0.2 mg/ml protease K) at room temperature. The gel plugs were transferred to 40 ml of Wash Buffer (10 mM Tris-HCl pH 8.0, 10 mM EDTA pH 8.0) and washed three times in a shaker at low speed for 1 hourrespectively to remove detergents. Gel plugs were transferred to 40 ml of PMSF Buffer (1.0 mM phenylmethylsulfonyl floride PMSF, 10 mM Tris-HCl pH 8.0, 10 mM EDTA pH 8.0) and incubated at room temperature for 1 hourwith gentle shaking; this process was repeated with fresh PMSF buffer. The plugs were then washed twice in 40 ml of Wash Buffer following incubation at 50°C for 20 minutes. The gel plugs were then transferred to 40 ml of 50 mM EDTA (pH 8.0) and stored at 4°C overnight. The DNA in the gel plugs was digested with 10 U of HindIII per plug at 37°C for 30 minutes.
The gel plugs with the partially digested DNA were cut into slurry. They were loaded onto a 1% agraose gel (Sigma, Type VII, low gelling temperature), and sealed on the top with agarose. Electrophoretic development occurred in 0.7 × TAE Buffer using a FIGE apparatus under Program 4 (BioRad, Hercules, CA 94547). Products ranging in size from 18 to 33 Kb were extracted from the gel (estimated 60 ng DNA total) following the protocol of GELase Agarose Gel-Digesting Preparation kit (Epicentre, Madison, WI 53713), and used to prepare the cosmid library.
The preparation of the cosmid library followed the procedures described by Bell et al. , with additional information described by Chow et al. . In brief to construct the cosmid library an estimated 60 ng of 18–33 Kb fragments recovered from gel were cloned into vector pCC1 which was digested with HindIII and then dephosphorylated with shrimp alkaline phosphatase followed the protocol (Roche, Indianapolis, IN 46250). The ligation products were packaged into bacteriophage particles using MaxPlax Lamda DNA packaging extracts (Epicentre, Madison, WI 53713) according to the protocol of the kit. Bacteriophage containing an estimated 5 × 103 particles in 50 μL were applied to infect 200 μl of EPI300 cells grown to exponential phase in LB liquid medium (Luria-Bertani medium) containing 10 mM MgSO4 and 0.2% maltose, which had been inoculated from the overnight culture grown in LB containing 10 mM MgSO4. After absorption following incubating at 37°C for 20 minutes, 1 ml of fresh LB medium was added and incubated for an additional 45 minutes. The infected cells were spread on LB 1% agar plates containing 12.5 μg/ml of chloramphenicol, 1 mM of IPTG and 40 μg/ml of X-gal for selection.
The cosmid library was used in two runs of 454 pyrosequencing . The first run was carried out on a GS 20 454 pyrosequencer, which gave read length around 90 basepairs (bps). The second run was done on a GS FLX 454 pyrosequencer, which gave reads length around 250 bps. Both pyrosequencing projects were done in the Interdisciplinary Center for Biotechnology Research at the University of Florida, Gainesville, USA. The reads obtained from the GS 20 and the GS FLX shotgun sequencing were separately assembled into contigs. These contigs were used in the analyses presented here.
Scanning electron microscopy
For scanning electron microscopic (SEM) D. magna was fixed in 3% glutaraldehyde in 0.1 M PB for 2 hours at 20°C. Sample was washed two times in distilled water for 5 to 10 seconds, dehydrated in graded ethanol series, and critical point dried (CPD) overnight (16 hours). The specimens were coated with gold (20 nm) and viewed using a Philips XL 30 ESEM under high volume conditions from 5 to 15 kv.
Sequences from the D. pulex, D. pulicaria and the two D. magna datasets included in this study are described in Table 2. Sequences were compared against the NCBI-nt database on nucleotide sequences using BLASTN  with the default settings in December 2007. Sequences longer than 1000 bps were divided into overlapping fragments around 500 bps. Sequences were homogenized to fragments of similar length so BLAST scores were comparable across different searches. Sequence comparison is computational challenging and was performed with an Opteron Linux high performance computer cluster established and maintained by the [BC]2 Basel Computational Biology Center at the Biozentrum University of Basel http://www.bc2.ch/center/index.htm. For the graphical presentation of the results we combined the two D. magna data sets.
For the analysis of the BLASTN results we used the metagenomics software MEGAN . This software allows exploring the taxonomic content of a sample based on the NCBI taxonomy. The blast files were imported into MEGAN using the import option BLASTN. The program then uses several thresholds to generate sequence-taxon matches. The "min-score" filter sets a bit-score cutoff value. The "top-percent" filter is used to retain hits whose scores lie within a given percentage of the highest bit score. The "min-support" filter is used to set a threshold for the minimum number of sequences that must be assigned to a taxon. We used all default parameter settings of the software (top-percent = 10, min-support = 2), except the minimal threshold for the bit score of hits, which were set at 100, following the recommendation of the authors . This reduces the number of reads assigned to a taxon, but avoids assignment based on weak homology. This analysis was done for all datasets between the 8. and the 11. January 2008.
While inspecting the data we ignored reads assigned to taxa other than plants and bacteria. Within the bacteria, we ignored the taxon Firmicutes (mostly gram-positive bacteria, many of which are endospore formers), because the two datasets of D. magna came from animals infected with the endospore forming pathogen, P. ramosa. The two other datasets (D. pulex and D. pulicaria), had only few sequences assigned to the Firmicutes (less than 0.2%). Thus, excluding the Firmicutes from the analysis did not influence the overall analysis.
In a separate analysis we manually inspected all four datasets for hits assigned to plant taxa (every taxon within and including the Viridiplantae), searching for hits to plastids (chloroplasts). For this analysis we set the MEGAN parameter minimum supported taxa to one.
Support for the preparation and characterization of cosmid DNA libraries for D. magna was provided by USDA/CSREES Project 50554, USDA/CSREES Multi-State Project NE1019, and the University of Florida IFAS Agricultural Experiment Station (CRIS Projects FLA-MCS-04353 and FLA-MCS-04080). The sequencing and portions of the analyses of the D. magna data were done at the Interdisciplinary Center for Biotechnology Research at the University of Florida, Gainesville, USA. We thank Li Liu for support and for the assembly of the contigs of the two D. magna datasets. The sequencing and portions of the analyses of the D. pulex and the D. pulicaria data were performed at the DOE Joint Genome Institute under the auspices of the U.S. Department of Energy's Office of Science, Biological and Environmental Research Program, and by the University of California, Lawrence Livermore National Laboratory under Contract No. W-7405-Eng-48, Lawrence Berkeley National Laboratory under Contract No. DE-AC02-05CH11231, Los Alamos National Laboratory under Contract No. W-7405-ENG-36 and in collaboration with the Daphnia Genomics Consortium (DGC) http://daphnia.cgb.indiana.edu. Additional analyses were performed by wFleaBase, developed at the Genome Informatics Lab of Indiana University with support to Don Gilbert from the National Science Foundation and the National Institutes of Health. Coordination infrastructure for the DGC is provided by The Center for Genomics and Bioinformatics at Indiana University, which is supported in part by the METACyt Initiative of Indiana University, funded in part through a major grant from the Lilly Endowment, Inc. We thank [BC]2 Basel Computational Biology Center at the Biozentrum University of Basel for hardware and software support. Our work benefits from, and contributes to the Daphnia Genomics Consortium. We are grateful to Daniel Mathys from the Zentrum für Mikroskopie Universität Basel for technical support with the SEM.
- Delwart EL: Viral metagenomics. Reviews in Medical Virology. 2007, 17 (2): 115-131. 10.1002/rmv.532.View ArticlePubMedGoogle Scholar
- Beardsley TM: Metagenomics reveals microbial diversity. Bioscience. 2006, 56 (3): 192-196. 10.1641/0006-3568(2006)056[0192:MRMD]2.0.CO;2.View ArticleGoogle Scholar
- Allen EE, Banfield JF: Community genomics in microbial ecology and evolution. Nature Reviews Microbiology. 2005, 3 (6): 489-498. 10.1038/nrmicro1157.View ArticlePubMedGoogle Scholar
- Streit WR, Schmitz RA: Metagenomics – the key to the uncultured microbes. Curr Opin Mircobiol. 2004, 7 (5): 492-498. 10.1016/j.mib.2004.08.002.View ArticleGoogle Scholar
- Venter JC, Remington K, Heidelberg JF, Halpern AL, Rusch D, Eisen JA, Wu DY, Paulsen I, Nelson KE, Nelson W, Fouts DE, Levy S, Knap AH, Lomas MW, Nealson K, White O, Peterson J, Hoffman J, Parsons R, Baden-Tillson H, Pfannkoch C, Rogers Y, Smith HO: Environmental genome shotgun sequencing of the Sargasso Sea. Science. 2004, 304 (5667): 66-74. 10.1126/science.1093857.View ArticlePubMedGoogle Scholar
- Bidle KD, Lee S, Marchant DR, Falkowski PG: Fossil genes and microbes in the oldest ice on Earth. Proc Natl Acad Sci USA. 2007, 104 (33): 13455-13460. 10.1073/pnas.0702196104.PubMed CentralView ArticlePubMedGoogle Scholar
- Edwards RA, Rodriguez-Brito B, Wegley L, Haynes M, Breitbart M, Peterson DM, Saar MO, Alexander S, Alexander EC, Rohwer F: Using pyrosequencing to shed light on deep mine microbial ecology. Bmc Genomics. 2006, 7: 57-10.1186/1471-2164-7-57.PubMed CentralView ArticlePubMedGoogle Scholar
- Cox-Foster DL, Conlan S, Holmes EC, Palacios G, Evans JD, Moran NA, Quan PL, Briese T, Hornig M, Geiser DM, Martinson V, vanEngelsdorp D, Kalkstein AL, Drysdale A, Hui J, Zhai J, Cui L, Hutchison SK, Simons JF, Egholm M, Pettis JS, Lipkin WI: A metagenomic survey of microbes in honey bee colony collapse disorder. Science. 2007, 318 (5848): 283-287. 10.1126/science.1146498.View ArticlePubMedGoogle Scholar
- Turnbaugh PJ, Baeckhed F, Fulton L, Gordon JI: Diet-induced obesity is linked to marked but reversible alterations in the mouse distal gut microbiome. Cell Host & Microbe. 2008, 3 (4): 213-223. 10.1016/j.chom.2008.02.015.View ArticleGoogle Scholar
- Booijink C, Zoetendal EG, Kleerebezem M, de Vos WM: Microbial communities in the human small intestine: coupling diversity to metagenomics. Future Microbiology. 2007, 2 (3): 285-295. 10.2217/174609126.96.36.1995.View ArticlePubMedGoogle Scholar
- Schmitt S, Wehrl M, Bayer K, Siegl A, Hentschel U: Marine sponges as models for commensal microbe-host interactions. Symbiosis. 2007, 44 (1–3): 43-50.Google Scholar
- Woyke T, Teeling H, Ivanova NN, Huntemann M, Richter M, Gloeckner FO, Boffelli D, Anderson IJ, Barry KW, Shapiro HJ, Szeto E, Kyrpides NC, Mussmann M, Amann R, Bergin C, Ruehland C, Rubin EM, Dubilier N: Symbiosis insights through metagenomic analysis of a microbial consortium. Nature. 2006, 443 (7114): 950-955. 10.1038/nature05192.View ArticlePubMedGoogle Scholar
- Leveau JHJ: The magic and menace of metagenomics: prospects for the study of plant growth-promoting rhizobacteria. European Journal of Plant Pathology. 2007, 119 (3): 279-300. 10.1007/s10658-007-9186-9.View ArticleGoogle Scholar
- Poinar HN, Schwarz C, Qi J, Shapiro B, MacPhee RDE, Buigues B, Tikhonov A, Huson DH, Tomsho LP, Auch A, Rampp M, Miller W, Schuster SC: Metagenomics to paleogenomics: Large-scale sequencing of mammoth DNA. Science. 2006, 311 (5759): 392-394. 10.1126/science.1123360.View ArticlePubMedGoogle Scholar
- Peters RH, Bernardi DR, eds: Daphnia. 1987, Verbania Pallanza: Consiglio Nazionale delle Ricerche Istituto Italiano di IdrobiologiaGoogle Scholar
- Green J: Parasites and epibionts of Cladocera. Trans Zool Soc Lond. 1974, 32: 417-515.View ArticleGoogle Scholar
- Ebert D: Ecology, epidemiology and evolution of parasitism in Daphnia. 2005, Bethesda (MD): National Library of Medicine (US), National Center for Biotechnology Information, [http://www.ncbi.nlm.nih.gov/books/bookres.fcgi/daph/screenA4.pdf]Google Scholar
- Chang N, Jenkins DG: Plastid endosymbionts in the freshwater crustacean Daphnia obtusa. J Crustac Biol. 2000, 20 (2): 231-238. 10.1651/0278-0372(2000)020[0231:PEITFC]2.0.CO;2.View ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.View ArticlePubMedGoogle Scholar
- Huson DH, Auch AF, Qi J, Schuster SC: MEGAN analysis of metagenomic data. Genome Research. 2007, 17 (3): 377-386. 10.1101/gr.5969107.PubMed CentralView ArticlePubMedGoogle Scholar
- Krause L, Diaz NN, Goesmann A, Kelley S, Nattkemper TW, Rohwer F, Edwards RA, Stoye J: Phylogenetic classification of short environmental DNA fragments. Nucleic Acids Research. 2008, 36 (7): 2230-2239. 10.1093/nar/gkn038.PubMed CentralView ArticlePubMedGoogle Scholar
- Pop M, Salzberg SL: Bioinformatics challenges of new sequencing technology. Trends Genet. 2008, 24 (3): 142-149.PubMed CentralView ArticlePubMedGoogle Scholar
- Raes J, Foerstner KU, Bork P: Get the most out of your metagenome: computational analysis of environmental sequence data. Curr Opin Mircobiol. 2007, 10 (5): 490-498. 10.1016/j.mib.2007.09.001.View ArticleGoogle Scholar
- McHardy A, Rigoutsos I: What's in the mix: phylogenetic classification of metagenome sequence samples. Curr Opin Mircobiol. 2007, 10 (5): 499-503. 10.1016/j.mib.2007.08.004.View ArticleGoogle Scholar
- Schloss PD, Handelsman J: A statistical toolbox for metagenomics: assessing functional diversity in microbial communities. Bmc Bioinformatics. 2008, 9:Google Scholar
- Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Berka J, Braverman MS, Chen YJ, Chen ZT, Dewell SB, Du L, Fierro JM, Gomes XV, Godwin BC, He W, Helgesen S, Ho CH, Irzyk GP, Jando SC, Alenquer ML, Jarvie TP, Jirage KB, Kim JB, Knight JR, Lanza JR, Leamon JH, Lefkowitz SM, Lei M, Li J, Lohman KL, Lu H, Makhijani VB, McDade KE, McKenna MP, Myers EW, Nickerson E, Nobile JR, Plant R, Puc BP, Ronan MT, Roth GT, Sarkis GJ, Simons JF, Simpson JW, Srinivasan M, Tartaro KR, Tomasz A, Vogt KA, Volkmer GA, Wang SH, Wang Y, Weiner MP, Yu P, Begley RF, Rothberg JM: Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005, 437 (7057): 376-380.PubMed CentralPubMedGoogle Scholar
- Fraune S, Bosch TCG: Long-term maintenance of species-specific bacterial microbiota in the basal metazoan Hydra. Proc Natl Acad Sci USA. 2007, 104: 13146-13151. 10.1073/pnas.0703375104.PubMed CentralView ArticlePubMedGoogle Scholar
- Bandi C, Damiani G, Magrassi L, Grigolo A, Fani R, Sacchi L: Flavobacteria as intracellular symbionts in cockroaches. Proc R Soc Lond B. 1994, 257: 43-48. 10.1098/rspb.1994.0092.View ArticleGoogle Scholar
- Hurst GDD, Hammarton TC, Bandi C, Majerus TMO, Bertrand D, Majerus MEN: The diversity of inherited parasites of insects: the male-killing agent of the ladybird beetle Coleomegilla maculata is a member of the Flavobacteria. Genet Res Camb. 1997, 70: 1-6. 10.1017/S0016672397002838.View ArticleGoogle Scholar
- Hurst GDD, Bandi C, Sacchi L, Cochrane AG, Bertrand D, Bernardet JF, Nakagawa Y, Holmes B, Karaca I, Majerus MEN: Adonia variegata (Coleoptera: Coccinellidae) bears maternally inherited Flavobacteria that kill males only. Parasitology. 1999, 118: 125-134. 10.1017/S0031182098003655.View ArticlePubMedGoogle Scholar
- Pinhassi J, Azam F, Hemphala J, Long R, Martinez J, Zweifel U, Hagström A: Coupling between bacterioplankton species composition, population dynamics, and organic matter degradation. Aquat Microb Ecol. 1999, 17: 13-26. 10.3354/ame017013.View ArticleGoogle Scholar
- Cottrell M, Kirchman D: Natural assemblages of marine proteobacteria and members of the Cytophaga-Flavobacter cluster consuming low- and high-molecular-weight dissolved organic matter. Appl Environ Microbiol. 2000, 66: 1692-1697. 10.1128/AEM.66.4.1692-1697.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Bernardet J, Segers P, Vancanneyt M, Berthe F, Kersters K, Vandamme P: Cutting a gordian knot: emended classification and description of the genus Flavobacterium, emended description of the family Flavobacteriaceae, and proposal of Flavobacterium hydatis nom. nov. (Basonym, Cytophaga aquatilis Strohl and Tait 1978). Int J Bacteriol. 1996, 46: 128-148.View ArticleGoogle Scholar
- Lampert W: Feeding and Nutrition in Daphnia. Mem Ist Ital Idrobiol. 1987, 45: 143-192.Google Scholar
- Wetzel RG: Limnology. 1975, Philadelphia, USA: Saunders College PublishingGoogle Scholar
- Cole J, Chai B, Farris R, Wang Q, Kulam-Syed-Mohideen A, McGarrell D, Bandela A, Cardenas E, Garrity G, Tiedje J: The ribosomal database project (RDP-II): introducing myRDP space and quality controlled public data. Nucleic Acids Res. 2007, 35: D169-D172. 10.1093/nar/gkl889.PubMed CentralView ArticlePubMedGoogle Scholar
- Chang HH, Shyu HF, Wang YM, Sun DS, Shyu RH, Tang SS, Huang YS: Facilitation of cell adhesion by immobilized dengue viral nonstructural protein 1 (NS1): Arginine-glycine-aspartic acid structural mimicry within the dengue viral NS1 antigen. J Infect Dis. 2002, 186 (6): 743-751. 10.1086/342600.View ArticlePubMedGoogle Scholar
- Bell KS, Avrova AO, Holeva MC, Cardle L, Morris W, DeJong W, Toth IK, Waugh R, Bryan GJ, Birch PRJ: Sample sequencing of a selected region of the genome of Erwinia carotovora subsp. atroseptica reveals candidate phytopathogenicity genes and allows comparison with Escherichia coli. Microbiology. 2002, 148: 1367-1378.View ArticlePubMedGoogle Scholar
- Chow V, Nong G, Preston JF: Structure, Function, and Regulation of the Aldouronate Utilization Gene Cluster from Paenibacillus sp. Strain JDR-2. J Bacteriol. 2007, 189: 8863-8870. 10.1128/JB.01141-07.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.