Skip to main content
  • Research article
  • Open access
  • Published:

Diet disparity among sympatric herbivorous cichlids in the same ecomorphs in Lake Tanganyika: amplicon pyrosequences on algal farms and stomach contents

Abstract

Background

Lake Tanganyika, an ancient lake in the Great Rift Valley, is famous for the adaptive radiation of cichlids. Five tribes of the Cichlidae family have acquired herbivory, with five ecomorphs: grazers, browsers, scrapers, biters and scoopers. Sixteen species of the herbivorous cichlids coexist on a rocky littoral slope in the lake. Seven of them individually defend feeding territories against intruding herbivores to establish algal farms. We collected epiphyton from these territories at various depths and also gathered fish specimens. Algal and cyanobacteria community structures were analysed using the amplicon-metagenomic method.

Results

Based on 454-pyrosequencing of SSU rRNA gene sequences, we identified 300 phototrophic taxa, including 197 cyanobacteria, 57 bacillariophytes, and 31 chlorophytes. Algal farms differed significantly in their composition among cichlid species, even in the same ecomorph, due in part to their habitat-depth segregation. The algal species composition of the stomach contents and algal farms of each species differed, suggesting that cichlids selectively harvest their farms. The stomach contents were highly diverse, even between species in the same tribe, in the same feeding ecomorph.

Conclusions

In this study, the amplicon-metagenomic approach revealed food niche separation based on habitat-depth segregation among coexisting herbivorous cichlids in the same ecomorphs in Lake Tanganyika.

Background

The cichlid species flock in Lake Tanganyika is a model system of adaptive radiation, defined as the rapid evolution of a multitude of species from a common ancestor as a consequence of their adaptation to various ecological niches. After the formation of the lake, 9 to 12 million years ago, several ancestral strains diversified into more than 200 species in 14 tribes [1]-[6]. Tanganyikan cichlids are unique because of the species richness of each of the ecomorphs [7],[8] that have evolved convergently and coexist in similar habitats, not only among different lakes but also within Lake Tanganyika [4]. Additionally, only six of the fourteen tribes present in Lake Tanganyika show sexual dichromatism; the others are monochromatic lineages [5],[9], in contrast to other African lakes in which the degree of cichlids’ adaptive radiation correlates highly with the presence of sexual dichromatism [5].

In Lake Tanganyika, five tribes of the family Cichlidae (Tropheini, Lamprologini, Ectodini, Eretmodini, and Tilapini) have acquired several ecomorphs that are closely related to feeding habits in herbivory such as grazing, browsing, scraping, biting, and scooping [10]-[13]. These tribes have no sexual dichromatism, and therefore, we can focus on the effect of ecological opportunity in the adaptive radiation of these herbivorous cichlids [14],[15]. Grazers, browsers, and scrapers are highly specialised and diversified, especially in these herbivorous habits [16]. Grazers comb unicellular algae from epilithic assemblages using multiple rows of similar-sized slender teeth with fork-like tricuspid tips [17],[18]. Browsers nip and nibble filamentous algae using their bicuspid teeth, which line the outermost of both jaws [11],[19]. Scrapers rub epiphyton from rock surfaces using several rows of chisel-like teeth [20]. The fishes of each ecomorph have distinct specialised trophic morphologies in their jaw structures and intestine lengths [21],[22], physiological abilities, such as secretion of digestive enzymes [21] and behaviours, such as cropping frequency [19].

In one rocky littoral area of the lake, 16 species of herbivorous cichlids coexist [23]. Seven of them individually defend their feeding territories against intruding herbivores to establish algal farms on which the territory-holder feeds [24]-[29]. Grazers and browsers utilise species-specific depth ranges, and species in the same depth range choose the substratum types on which they feed [10],[23],[30]. In particular, territorial cichlids of the same ecomorph exhibit separate habitat depths [29],[31],[32], and ecological character displacement in habitat depth and associated body shape were reported in sympatric populations of Tropheus moorii and Tropheus polli, compared to isolated populations of the two species [33].

Niche segregation leading to multi-species coexistence has at least three axes: space, food and time [34]. These herbivorous cichlids are diurnally active and do not have separate active times during the day [35]. In this study, we focused on food and habitat depth. How do feeding specialisation and habitat-depth selection cause food-source segregation, thereby enabling the coexistence of multiple species in the same ecomorph in this adaptive radiation system? To answer this question, the degree of diet disparity in algae and cyanobacteria between herbivorous cichlid species was analysed at the species-level using a metagenomic amplicon sequencing approach. Based on 454-pyrosequencing of SSU rRNA gene sequences, we described algal and cyanobacterial compositions in algal farms and the stomach contents of herbivorous cichlid species on a rocky slope of Lake Tanganyika. The composition of algal farms was compared among ecomorphs and among species within each ecomorph, and the effect of habitat depth on the composition was analysed. The stomach contents were compared among ecomorphs and among species within each ecomorph, and were compared with algal farm composition. Finally, the relationships between the phylogenetic distances between fish species and the similarity of algal farms and stomach contents were analysed.

Results

We sampled 16 species of herbivorous cichlids (Oreochromis tanganicae, Xenotilapia papilio, Eretmodus cyanostictus, Telmatochromis temporalis, Telmatochromis vittatus, Variabilichromis moorii, Interochromis loocki, Pseudosimochromis curvifrons, Petrochromis famula, Petrochromis fasciolatus, Petrochromis macrognathus, Petrochromis polyodon, Petrochromis horii, Petrochromis trewavasae, Simochromis diagramma, and Tropheus moorii) from Kasenga Point (8°43′ S, 31°08′ E) near Mpulungu, Zambia, on the southern tip of Lake Tanganyika in November 2010 (Table 1). Periphyton samples were simultaneously collected from five territories each of T. temporalis, V. moorii, P. macrognathus, P. polyodon, P. horii, P. trewavasae, T. moorii, dominant males of I. loocki and P. curvifrons (in these species females and inferior males do not have feeding territory, Table 1), and breeding pairs of E. cyanostictus and X. papilio. Another 13 periphyton samples were collected from outside the cichlid territories. Total DNA was extracted from each of the stomach content samples of the cichlids and periphyton samples, and pyrosequencing of SSU rRNA gene was conducted. The resulting 5,073 consensus sequences represented operational taxonomic units (OTUs; Additional file 1: Data S1). Of the 5,073 consensus reads, 2,663 reads were excluded as singletons because such sequences are putatively erroneous. Taxonomic assignment of the OTUs was performed by the Query-centric auto-k-nearest-neighbour (QCauto) method, using the software Claident v. 0.1.2012.05.21 [36],[37].

Table 1 Herbivorous cichlids in Lake Tanganyika and their ecomorphs, territoriality and the number of samples

In total, we identified 300 phototrophic OTUs (see Additional file 1: Data S1). The algal farms of cichlids, periphytons outside the territories and stomach contents comprised 31.6 ± 12.2 phototrophic OTUs (average ± SD), 30.1 ± 17.9 OTUs and 14.8 ± 15.6 OTUs, respectively. Of the 300 phototrophic OTUs observed, 197 were cyanobacteria, 57 were bacillariophytes, 31 were chlorophytes and the other 15 comprised euglenides, eustigmatophyceae, streptophytes, dinoflagellates, a rhodophyta and unknown viridiplantae (see Additional file 2: Table S1). Of the eukaryotic sequences, 47 were eukaryotic 18S sequences, 15 were chloroplast 16S sequences, 28 were mitochondrial 12S and 16S sequences, and 13 remained uncertain (see Additional file 1: Data S1). To reduce variance in alpha-diversity among samples that resulted from differences in sequencing effort (that is, variations in the number of sequencing reads among samples), each sample was rarefied to 40 reads using the rrarefy function in vegan v.2.0 (see Additional file 3: Data S2). We determined that the taxonomic composition within the sub-sampled 40 reads would provide a semi-quantitative measure of relative biomass within each sample.

Phototroph composition in algal farms varied significantly among cichlid species (Figure 1A; Adonis, P <0.05, Additional file 4: Table S2). Phototroph composition in stomach contents also varied significantly among cichlid species (Figure 1B; Adonis, P <0.01, Additional file 4: Table S2). These results indicate that variations within species are significantly smaller than variations among species. Among species of the same feeding ecomorph, the phototroph compositions in both their defended farms and stomach contents were widely scattered.

Figure 1
figure 1

Principal coordinate analysis among algal farms of herbivorous cichlids (A) and their stomach contents (B). Axes represent the first two principal coordinates maximising the variance in the data (PC 1 and PC 2). Dissimilarity was calculated using the Canberra distance index. The percentage of the total variance explained by each axis estimated under the broken stick model is shown in parentheses. Abbreviations are listed in Table 1. Orange, green, dark red, grey, and yellow coloured species indicate grazers, browsers, scrapers, scoopers and biters, respectively. Out = periphyton outside the cichlid territories.

Phototroph composition in algal farms and stomach contents mostly consisted of three phyla: Chlorophyta, Bacillariophyta and Cyanobacteria (Figure 2). Bacillariophytes were relatively more abundant in the stomachs of grazers compared to both the stomachs of browsers and the algal farms of grazers, indicating that grazers selectively ingested bacillariophytes. The stomach contents of browsers contained more cyanobacteria and chlorophytes than those of grazers. Cladophora sp. (OTU #3295) dominated the algal farms of all herbivorous cichlids, but most grazing herbivores seldom ingested this algal species. Only a territorial browser (T. temporalis) and non-territorial browser (T. vittatus) and a scraper (E. cyanostictus) ingested this alga. No cyanobacteria or algae other than the Cladophora sp. were dominant throughout the algal farms of herbivorous cichlids. As shown by principal coordinate analysis (PCoA) ordination, the phototroph compositions of algal farms and stomach contents differed considerably among species within the same feeding ecomorph.

Figure 2
figure 2

Phototroph compositions inside and outside the territories of herbivorous cichlids and their stomach contents. A) Algal farms of herbivorous cichlids and periphyton outside the cichlid territories; B) stomach contents of the cichlid fishes. The compositions were averaged over one to five samples. The number of samples is shown in Table 1 and Figure 4.

The habitat depths of cichlids differed significantly among species (Figure 3, Additional file 5: Table S3), and the depth significantly affected the phototrophic composition (Adonis, P <0.001, Additional file 4: Table S2). Two grazers (P. macrognathus and P. polyodon), two browsers (P. curvifrons and T. temporalis) and a scraper (E. cyanostictus) inhabited the shallowest zone, two grazers (I. loocki and P. trewavasae), and two browsers (V. moorii and T. moorii,) and a scooper (X. papilio) inhabited the intermediate depth, and a grazer (P. horii) occupied the deepest zone (Figure 3A). In addition to the depth segregation between cichlid species, phototrophic OTUs also showed significant variation in habitat depth (Figure 3B, Additional file 6: Table S4, see Additional file 7: Figure S1 for occurrence frequency of the dominant OTUs in various depths). The shallower inhabitants of the orders Nostocales (OTU #551) and Oscillatoriales (#4699), the subclass Oscillatoriophycideae (#3283), and two cyanobacteria (#553 and #661) occurred frequently in algal farms of the cichlids in the sallowest zone. The relatively deep inhabitants, Cladophora sp. (#3295), a Stigonematales (#3375), and a diatom (#3935) occurred in the territories of cichlids living in the intermediate depths. A cladophorales (#293) and a cyanobacteria (#485) occurred frequently in the farms of the deepest inhabitant, P. horii. Consequently, the standardised specialisation index (d’) was significantly higher in the algal farms of cichlid species, with the exceptions of I. loocki and T. temporalis, than in a random network provided by a null model.

Figure 3
figure 3

Habitat depths of cichlids (A) and phototrophic OTUs (B). Shared letters on box plots indicate no significant differences, and pairs that do not share any letters in common were significantly different by the Tukey post hoc test between cichlid species and between OTUs in (A) and (B), respectively. Species abbreviations are shown in Table 1. Ecomorphs of cichlids abbreviated in parentheses are as follows: G, grazer; B, browser; Sr, scraper; So, scooper. Colours of OTUs indicate their phyla, and patterns are the same as in Figure 2. Networks between algal farms of cichlid species and OTUs denote the frequency of occurrence with square-root transformation. Asterisks indicate significance in standardised specialisation based on d’ measures: *=5%, **=1%, ***=0.1%.

The bipartite graphs show the occurrence of each phototrophic OTU in the algal farms of herbivorous cichlid species (upper bipartite graphs), and the occurrence of those OTUs in the stomachs of the cichlid species (lower bipartite graphs) living in the shallower area (2.4 ± 1.6 m depth, Figure 4A) and deeper area (11.2 ± 5.4 m depth, Figure 4B). The network structures of the upper and lower graphs were quite different in both shallow and deep inhabitants (similarity matrices were compared between algal farms and stomach contents using the Mantel test, for both shallow and deep inhabitants, z =10.72 and 7.52, respectively, both P >0.05), which indicates that the cichlid species did not ingest their algal farms randomly, but they utilised phototrophs selectively. Further, for the upper graphs, the network level measure of specialisation H 2 values (0.20 and 0.27 in shallow and deep inhabitants, respectively) were significantly higher than those calculated for the null models (0.08 ± 0.01 and 0.10 ± 0.01 in shallow and deep inhabitants, respectively; both P <0.001), suggesting that a specific structure of phototrophs was established in the algal farms of each cichlid species within the same depth range. For the lower graphs, H 2 values (0.31 and 0.40 in shallow and deep inhabitants, respectively) were also significantly higher than those calculated for the null models (0.15 ± 0.01 and 0.17 ± 0.02 in shallow and deep inhabitants, respectively; both P <0.001), suggesting that the composition of stomach contents varied widely among cichlid species within the same depth range.

Figure 4
figure 4

Networks between herbivorous cichlids and phototrophs in their territories and between phototrophs inside stomachs and fish species. A) shallow inhabiting cichlids, B) deep inhabiting cichlids. Species abbreviations are shown in Table 1. Out = outside the territories. Link width indicates the frequency of occurrence of the OTU in algal farms of cichlids or their stomach contents. Algal and cyanobacteria OTUs that occurred less than 10 times in total are not shown to decrease complexity. Colours of OTUs indicate their phyla, and patterns are the same as in Figure 2. Asterisks indicate significance in standardised specialisation based on d’ measures: *=5%, **=1%, ***=0.1%. These networks were constructed using the composition averaged over one to five samples, and the numbers of samples are shown in the rectangles.

The lower bipartite graphs in Figure 4A,B show that different phototrophic species or different combinations of phototrophic species contributed to the stomach contents of different cichlid species even within the same depth range. All the cichlid species both living in shallow and deep zones had significant specific links with various phototrophic OTUs. Neither dissimilarity in the algal farm composition or stomach contents in species-pairs was related to the phylogenetic distance of the pairs (Figure 5, Mantel test: z = 1.63 for algal farm versus cichlid phylogeny; z = 3.19 for stomach contents versus cichlid phylogeny; both P >0.05), even between species in the Tropheini tribe, which comprises a single lineage of grazers and a few lineages of browsers (Mantel test: z = 0.35 for algal farm; z = 0.51 for stomach contents; both P >0.05). Dissimilarity in algal farms and stomach contents did not differ between same-tribe pairs and different-tribe pairs, or between same feeding-ecomorph pairs and different ecomorph pairs (generalized linear model (GLM), P >0.05 for both factors, Additional file 8: Table S5). On the other hand, dissimilarity in algal farm composition was positively correlated with differences in habitat depth (Figure 6; Mantel test, z = 10620.73, P <0.001; GLM, P <0.01, Additional file 8: Table S5). Dissimilarity in stomach contents of all the cichlids was not related to the difference in their habitat depths.

Figure 5
figure 5

Pairwise distance-contrast plots between the phylogenetic distance of fishes versus dissimilarity in phototroph compositions. A) The correlation between the phylogenetic distance of fish species versus dissimilarity in the phototroph compositions of algal farms defended by the herbivorous cichlids, and B) pairwise plots between the phylogenetic distance versus dissimilarity in phototroph composition of stomach contents of the herbivorous cichlids. Circles and squares indicate species-pairs of the same feeding ecomorph and of different ecomorphs, respectively. Closed and open symbols indicate species-pairs of the same tribe and of different tribes, respectively. Dissimilarity was calculated using the Canberra distance index. Species abbreviations are listed in Table 1.

Figure 6
figure 6

Pairwise distance-contrast plots between the difference in habitat depth versus dissimilarity in phototroph compositions. A) The correlation between the difference in habitat depth of fish species versus the dissimilarity in the phototroph compositions of algal farms defended by the herbivorous cichlids, and B) pairwise plots between the difference in habitat depth versus the dissimilarity in phototroph composition of stomach contents of the herbivorous cichlids. Circles and squares indicate species-pairs of the same feeding ecomorph and of different ecomorphs, respectively. Closed and open symbols indicate species-pairs of the same tribe and of different tribes, respectively. The solid line is the fitted line for GLM using the difference in depth as a fixed factor. Dissimilarity was calculated using the Canberra distance index. Species abbreviations are listed in Table 1. GLM, generalized linear model.

Discussion

Phototroph composition in algal farms defended by herbivorous cichlids significantly varied among cichlid species, even among species in the same ecomorphs. This was due in part to habitat depth segregation of cichlids, as the dissimilarity of algal farms increased significantly when the difference in habitat depth increased. Segregation of habitat range by depth among cichlid species has been described in several species pairs of the same ecomorphs [10],[23],[30],[33]. Additionally, the stomach contents of these cichlids were segregated by at least 0.75 in the Canberra dissimilarity index, even between species in the same feeding ecomorph inhabiting the same depth range. The stomach content compositions were quite different from those of the algal farms of the same species. These results indicate that these cichlid species are highly selective when ingesting phototrophs from their specific algal farms, resulting in specific utilisation of phototrophs by each cichlid species. Consequently, similarity in algal farms was relatively high in the same depth range, but the stomach contents differed among species living in the same depth range. Stomach contents of cichlids were so diverse that no significant effect of ecomorph, tribe, habitat depth, and their interactions on the similarities of the stomach contents was detected. These territorial cichlids ingest algae from specific points inside their territories [10],[17],[48], and, therefore, further study on phototrophic distribution within territories and detailed feeding sites will reveal the mechanism that enables diet segregation among cichlid species of each ecomorph living in the same depth range.

Stomachs of grazers had higher compositions of bacillariophytes, such as Gomphonema sp., Encyonema sp., and Bacillariophyceae OTU #3585, suggesting that grazers ingest these diatoms selectively. This result is in agreement with previous reports of high proportions of epiphytic unicellular algae in the stomachs of the grazers [10],[13],[48],[49], and with a study that found high activity of a digestive enzyme, laminarinase, which hydrolyses laminarine, the main polysaccharide of diatoms, in the grazer P. orthognathus[21]. On the other hand, browsers ingested more cyanobacteria, such as Chroococcales, Oscillatoriales and chlorophytes, such as Scenedesmaceae. This result is in agreement with previous reports of high proportions of filamentous algae and cyanobacteria in the stomachs of browsers based on visual observations of stomach contents [10],[13],[48],[49].

A filamentous green alga, Cladophora sp., dominated the farms of most cichlid species, but was not ingested often. Cladophora are known to be unpalatable due to their chemical defences and poor amino acid contents [50],[51]. Therefore, this alga seems to be utilised as a substratum harbouring epiphytic diatoms for grazers and is a less preferred food item for other herbivorous cichlids.

In Lake Tanganyika, multiple species of herbivorous cichlids coexist sympatrically, but the number of ecomorphs is limited; therefore, the lake has been thought to be a `species-saturated’ community where the number of species exceeds the number of available niches [4]. This study, however, shows that these cichlids segregate their habitats on a finer scale, just as multiple Anolis lizard species of the same ecomorph specialise in specific thermal microhabitats or specific prey sizes and coexist in the Caribbean islands [52]-[54].

Our results demonstrate that phylogenetically close cichlid species neither always defend similar or opposite-structured algal farms nor ingest similar or opposite-structured phototrophic compositions. This confirms the hypothesis that niche lability exceeds niche conservatism in a community in which species have a long evolutionary history of ecological interaction, as seen in the adaptive radiation of Caribbean Anolis lizards [55]. This also suggests that the flexibility to specialise in microhabitats may drive their divergence and coexistence.

Conclusions

To date, it has been difficult to classify unicellular and filamentous algae, especially diatoms and cyanobacteria, to the species level, especially in stomach contents, and evidence of dietary segregation in herbivorous cichlids in the same feeding ecomorph is scarce [4],[56]. This study revealed that the phototrophic composition of algal farms was quite different among species, even those in the same feeding ecomorph, because of habitat depth segregation among species. In addition, cichlid species selectively fed on phototrophs from their algal farms. As a result, the algal compositions of stomachs differed among species, even those in the same ecomorph inhabiting the same depth range. Therefore, a metagenomic approach revealed food niche separation based on habitat-depth segregation and food preference among coexisting herbivorous cichlids in the same ecomorphs on a rocky shore of Lake Tanganyika.

Methods

Sampling

We sampled 16 species of herbivorous cichlids (Oreochromis tanganicae, Xenotilapia papilio, Eretmodus cyanostictus, Telmatochromis temporalis, Telmatochromis vittatus, Variabilichromis moorii, Interochromis loocki, Pseudosimochromis curvifrons, Petrochromis famula, Petrochromis fasciolatus, Petrochromis macrognathus, Petrochromis polyodon, Petrochromis horii, Petrochromis trewavasae, Simochromis diagramma, and Tropheus moorii) from Kasenga Point (8°43′ S, 31°08′ E) near Mpulungu, Zambia, on the southern tip of Lake Tanganyika in November 2010, using a gill net (Table 1). Stomach contents were extracted from one to five individuals of each species and immediately preserved in 100% ethanol for molecular analysis. Periphyton samples were simultaneously collected by scuba divers from five territories each of T. temporalis, V. moorii, P. macrognathus, P. polyodon, P. horii, P. trewavasae, T. moorii, dominant males of I. loocki and P. curvifrons (in these species females and inferior males do not have feeding territory, Table 1) and breeding pairs of E. cyanostictus and X. papilio. Another 13 periphyton samples were collected from outside the cichlid territories using separate sampling bottles to avoid cross-contamination during collection. Each periphyton sample was collected from the territory of different cichlid individuals, and the depth of the sampling location was measured to the nearest 0.1 m using a diving computer. We defined the territory as the area where the territory holder fed on and defended against conspecific and heterospecific herbivores [57]. Whether a site was located within or outside of a cichlid fish territory was determined by 20 minutes of observation immediately prior to sampling. Periphyton samples were immediately preserved in 100% ethanol and stored at −25°C in the laboratory until molecular analysis.

Metagenomic amplicon sequencing

DNA extraction

From each algal farm or stomach content sample, a 2-mm3 subsample was pulverised by beating with 4-mm zirconium beads 20 times per second for 2 minutes using a TissueLyser II (Qiagen, Venlo, The Netherlands). Total DNA was extracted from the crushed tissue using the cetyltrimethylammonium bromide (CTAB) method [58].

Universal primers, PCR, and pyrosequencing

We sequenced SSU rRNA genes using tag-encoded massively parallel pyrosequencing analysis [59]-[61]. We aimed to focus on all phototrophs, including cyanobacteria and various algal taxa. The SSU rRNA gene is effective for DNA barcoding of organisms, especially bacteria, because of the massive public database [62]. The primer set, CYA_ALG_F (5′-AAA CTC AAA GRA ATT GAC GG-3′; Escherichia coli position 906 to 25) and CYA_ALG_R (5′-GYT ACC TTG TTA CGA CTT C-3′; E. coli position 1490 to 1508), was designed to amplify the highly variable V6 region of the 16S rRNA gene [63]. The coverage of the primers was evaluated based on in silico polymerase chain reaction (PCR) using the Probe Match function in the Ribosomal Database Project (RDP) release 11.2 [64] on 10 June 2014, and the TestPrime 1.0 function [65] within the SILVA release 115 database [62] on 22 May 2014. The coverage of the newly designed primer set was much higher than that of two other primer sets designed for `phyto-specific’ 16S rRNA gene amplification [66],[67] (see Additional file 9: Figure S2). Although the coverage for eukaryotic algal taxa was limited in the TestPrime results, these coverage values may be underestimated because our focus was on the 16S rRNA gene of plastids, while the eligible sequences in the SILVA database include the 18S rRNA genes of eukaryotes. The primer coverage was determined to be sufficient for demonstrating the differences between the algal/cyanobacterial communities in the algal farms and stomach contents of cichlid species.

For each sample, SSU rRNA gene sequences were amplified using the forward primer CYA_ALG_F fused with the 454 pyrosequencing Adaptor A (5′-CCA TCT CAT CCC TGC GTG TCT CCG ACT CAG-3′) and the 8-mer molecular ID [61], and the reverse primer CYA_ALG_R, with the Ampdirect Plus (Shimadzu Corp., Kyoto, Japan) buffer system and BIOTAQ HS DNA Polymerase (Bioline, London, U.K.). All forward primers were tested using AmplifX_1.5.4, developed by Nicolas Jullien [68] to avoid potential primer-dimers. PCR was conducted using a temperature profile of: 95°C for 10 minutes; 30 cycles at 94°C for 20 seconds, 48°C for 30 seconds, and 72°C for 30 seconds; and a final extension at 72°C for 7 minutes. For each PCR amplicon, DNA concentrations were measured using a Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA, USA) and were diluted to 4 ng/μl with Milli-Q water. PCR amplicons were pooled and purified using ExoSAP-IT (GE Healthcare, Little Chalfont, Buckinghamshire, U.K.) and a QIAquick PCR Purification Kit (Qiagen). The SSU rRNA amplicons were subjected to pyrosequencing using a GS Junior sequencer (Roche, Basel, Switzerland). The sequencing was conducted according to the manufacturer’s instructions.

Assembling of pyrosequencing reads

In total, we obtained 92,559 SSU rRNA gene sequence reads from the pyrosequencing run. We trimmed low-quality 3′ tails, with a minimum quality value of 27. We discarded short reads that were less than 150 base pairs, excluding the forward primer and molecular ID positions. After the elimination of possible chimeras using UCHIME v4.2.40 [69] with a minimum score of 0.1 to report a chimera, a total of 87,898 reads remained [DDBJ: DRA002208]. Reads were assembled using Assams-assembler v0.1.2012.05.24 [70],[71], which is a highly parallelised extension of the Minimus assembly pipeline [72]. All reads were assembled with a minimum cut-off similarity of 97% to remove pyrosequencing errors. The resulting 5,073 consensus sequences represented OTUs (see Additional file 1: Data S1). Of the 5,073 consensus reads, 2,663 reads were excluded as singletons because such sequences are putatively erroneous. The rarefaction curves, which were drawn using the R package vegan_2.0 [73], showed that the number of OTUs did not reach saturation for many of the samples (see Additional file 10: Figure S3); however, the rare OTUs may have contributed little as a dietary source for herbivorous fishes.

Molecular identification of cyanobacteria and algae

Taxonomic assignment of the OTUs was performed by the QCauto method [37], which is known to return the most accurate taxonomic identification results among the existing automated DNA barcoding methods [37], using the software Claident v. 0.1.2012.05.21 [36],[37]. The `all_genus’ and `all_class’ reference sequence databases that were updated on 5 February 2014 and provided by Claident were used as references. The `all_genus’ and `all_class’ datasets contain all the sequences of the genus or class, respectively, and therefore, `all_genus’ was used preferentially, and `all_class’ was used complementarily. The benefit of using the QCauto method is that it enables accurate and fully automated taxonomic identification based on BLAST + searches [74] without setting any arbitrary thresholds of sequence identity percentages or E-values [37]. In the taxonomic identification process, the relaxed lowest common ancestor algorithm [60],[75] was used. We extracted all phototrophic taxa, as shown in Additional file 1: Data S1. The origins of eukaryotic OTUs were assigned using the Metaxa software to determine whether an OTU originated from nucleic 18S rRNA genes, chloroplast 16S rRNA genes or mitochondrial 12S/16S rRNA genes [76]. The numbers of phototrophic reads for each periphyton and stomach content sample were 159.4 ± 95.1 (average ± SD) and 49.8 ± 76.4, respectively.

To reduce variance in alpha-diversity among samples that resulted from differences in sequencing effort (that is, variations in the number of sequencing reads among samples), we rarefied each sample to 40 reads using the rrarefy function in vegan v.2.0 (see Additional file 3: Data S2). This rarefying method is common for comparisons among samples with various numbers of reads (for example, [77],[78]). Although DNA copy number will often differ between species, and primer efficiency will differ among DNA sequences, the number of DNA molecules detected from faeces or digestive tract contents by quantitative PCR can be used as a semi-quantitative measure of prey biomass [79]. Therefore, we determined that the taxonomic composition within the sub-sampled 40 reads would provide a semi-quantitative measure of relative biomass within each sample.

Network analysis and measure of standardised specialisation

The networks between herbivorous cichlid species and the phototrophs that occurred in their territories, and the networks between cichlid species and phototrophs found inside their stomachs were constructed among species in shallower areas and species in deeper areas respectively, using the R package, bipartite_1.18 [80]. To quantitatively evaluate the association specificity, the d’ index of the degree of specialisation in the network was calculated [81],[82]. The number of reads was summed among samples of each species within each depth zone to form a data matrix for this analysis because here we aimed to determine the species-level interactions between cichlids and phototrophs. Note that intraspecific variation was significantly smaller than interspecific variation, as described below. The observed d' index values were compared with those of randomised links in the `vaznull’ model [83] with 10,000 permutations. A d’ index higher than expected by chance indicates association specificity between a cichlid species and phototrophic OTUs, or specificity between a phototrophic OTU and a cichlid species. H 2 [81], a network-level measure of specialisation based on the deviation of a species, realised number of interactions from the species, expected total number of interactions, was also calculated in all bipartite graphs. Observed values were compared with those of the vaznull models with 10,000 permutations. The bipartite package was also used for calculation of indices and making null models [84].

Statistical analyses

Multivariate analyses of phototroph community structures in cichlid territories were conducted by permutational multivariate analysis of variance using distance matrices calculated with the Canberra index and Bray-Curtis dissimilarities (Adonis; [85]) with the factors of habitat depth and cichlid species. The community structures of stomach contents were also analysed by Adonis, with a factor of cichlid species. The number of null permutations was set to 9,999. PCoA ordination was conducted using the Canberra index and Bray-Curtis dissimilarities among the OTU composition data because PCoA plotting using Canberra distance is recommended for finding clusters in microbial community samples [86]. To compare the habitat depths of each cichlid species, a GLM analysis was conducted with gamma distribution, using the habitat depth as a response variable and cichlid species as an explanatory variable. The Tukey post hoc test in R [87] was performed in the multcomp_1.3.3 package (function glht) for the comparisons among cichlid species. To compare the habitat depths of each phototrophic OTU, generalised linear mixed model (GLMM) analysis was conducted with gamma distribution using the depth as a response variable, OTU as a fixed-effect explanatory variable and the cichlid species from whose algal farms the OTU samples were collected as a random-effect explanatory variable. The Tukey post hoc test was performed for the comparisons among OTUs. To estimate the extent of divergence and convergence in phototroph utilisation within herbivorous cichlids, we compared the phylogenetic distance between each pair of species to their similarity in phototroph composition of algal farms/stomach contents using the Mantel test with 9,999 permutations. The phylogenetic distance was calculated by the cophenetic () function in the ape package for R based on the mitochondrial ND2 and nuclear ednrb1 and phpt1 sequences [4]. The Canberra distance in phototroph composition of algal farms and those of stomach contents were compared using the Mantel test with 9,999 permutations, and the dissimilarities of algal farms and stomach contents were compared with differences in habitat depth using the Mantel test with 9,999 permutations. To test the effect of the differences in tribe, ecomorph and habitat depth on the similarity of algal farm and stomach contents composition, a GLM analysis was conducted with gamma distribution, using the Canberra distance index as a response variable and the differences in tribe, ecomorph and habitat depth as explanatory variables. All analyses were conducted using R_3.1.0 [88]. The two indices of dissimilarity resulted in similar outcomes in all analyses; therefore, the results using Canberra are shown hereafter, and those using Bray-Curtis are provided in Additional file 11: Figure S4, Additional file 12: Figure S5, Additional file 13: Figure S6.

Authors’ contributions

HH, HT, MK, MH participated in the design of the study. MH, MK arranged the permission and conditions of field research. HH, HM conducted field research. HH, HT, AST, SY carried out the molecular genetic studies. HH, HT, AST analysed data. HH wrote the manuscript with input from all authors. All authors read and approved the final manuscript.

Additional files

References

  1. Salzburger W, Meyer A, Baric S, Verheyen E, Sturmbauer C: Phylogeny of the Lake Tanganyika cichlid species flock and its relationship to the Central and East African haplochromine cichlid fish faunas. Syst Biol. 2002, 51: 113-135. 10.1080/106351502753475907.

    Article  PubMed  Google Scholar 

  2. Salzburger W, Mack T, Verheyen E, Meyer A: Out of Tanganyika: genesis, explosive speciation, key-innovations and phylogeography of the haplochromine cichlid fishes. BMC Evol Biol. 2005, 5: 17-10.1186/1471-2148-5-17.

    Article  PubMed Central  PubMed  Google Scholar 

  3. Danley PD, Husemann M, Ding B, Dipietro LM, Beverly EJ, Peppe DJ: The impact of the geologic history and paleoclimate on the diversification of East African cichlids. Int J Evol Biol. 2012, 2012: 574851-10.1155/2012/574851.

    Article  PubMed Central  PubMed  Google Scholar 

  4. Muschick M, Indermaur A, Salzburger W: Convergent evolution within an adaptive radiation of cichlid fishes. Curr Biol. 2012, 22: 2362-2368. 10.1016/j.cub.2012.10.048.

    Article  CAS  PubMed  Google Scholar 

  5. Wagner CE, Harmon LJ, Seehausen O: Ecological opportunity and sexual selection together predict adaptive radiation. Nature. 2012, 487: 366-369. 10.1038/nature11144.

    Article  CAS  PubMed  Google Scholar 

  6. Takahashi T: Systematics of Tanganyikan cichlid fishes (Teleostei: Perciformes). Ichthyol Res. 2003, 50: 367-382. 10.1007/s10228-003-0181-7.

    Article  Google Scholar 

  7. Williams EE: The origin of faunas. Evolution of lizard congeners in a complex island fauna: a trial analysis. Evol Biol. 1972, 6: 47-89. 10.1007/978-1-4684-9063-3_3.

    Article  Google Scholar 

  8. Losos JB: Adaptive radiation, ecological opportunity, and evolutionary determinism. Am Nat. 2010, 175: 623-639. 10.1086/652433.

    Article  PubMed  Google Scholar 

  9. Tsuboi M, Gonzalez-Voyer A, Höglund J, Kolm N: Ecology and mating competition influence sexual dimorphism in Tanganyikan cichlids. Evol Ecol. 2012, 26: 171-185. 10.1007/s10682-011-9489-3.

    Article  Google Scholar 

  10. Takamura K: Interspecific relationships of aufwuchs-eating fishes in Lake Tanganyika. Environ Biol Fishes. 1984, 10: 225-241. 10.1007/BF00001476.

    Article  Google Scholar 

  11. Mbomba N: Comparative morphology of the feeding apparatus in cichlidian algal feeders of Lake Tanganyika. Afr Stud Monogr. 1983, 3: 1-23.

    Google Scholar 

  12. Mbomba N: Comparative feeding ecology of aufwuchs eating cichlid fishes in Lake Tanganyika with reference to their developmental changes. Physiol Ecol Japan. 1986, 23: 79-108.

    Google Scholar 

  13. Hori M, Gashagaza MM, Nshombo M, Kawanabe H: Littoral fish communities in Lake Tanganyika: irreplaceable diversity supported by intricate interactions among species. Conserv Biol. 1993, 7: 657-666. 10.1046/j.1523-1739.1993.07030657.x.

    Article  Google Scholar 

  14. Gante HF, Salzburger W: Evolution: cichlid models on the runaway to speciation. Curr Biol. 2012, 22: R956-R958. 10.1016/j.cub.2012.09.045.

    Article  CAS  PubMed  Google Scholar 

  15. Salzburger W: The interaction of sexually and naturally selected traits in the adaptive radiations of cichlid fishes. Mol Ecol. 2009, 18: 169-185. 10.1111/j.1365-294X.2008.03981.x.

    Article  PubMed  Google Scholar 

  16. Koblmüller S, Egger B, Sturmbauer C, Sefc KM: Rapid radiation, ancient incomplete lineage sorting and ancient hybridization in the endemic Lake Tanganyika cichlid tribe Tropheini. Mol Phylogenet Evol. 2010, 55: 318-334. 10.1016/j.ympev.2009.09.032.

    Article  PubMed  Google Scholar 

  17. Yamaoka K: Morphology and feeding behaviour of five species of genus Petrochromis (Teleostei, Cichlidae). Physiol Ecol Japan. 1982, 19: 57-75.

    Google Scholar 

  18. Yamaoka K: Comparative osteology of the jaw of algal-feeding cichlids (Pisces, Teleostei) from Lake Tanganyika. Reports of the Usa Marine Biological Institute, Kochi University. 1987, 9: 87-137.

    Google Scholar 

  19. Yamaoka K: Feeding behaviour and dental morphology of algae scraping cichlids (Pisces: Teleostei) in Lake Tanganyika. Afr Stud Monogr. 1983, 4: 77-89.

    Google Scholar 

  20. Yamaoka K, Hori M, Kuratani S: Ecomorphology of feeding in “goby-like” cichlid fishes in Lake Tanganyika. Physiol Ecol Japan. 1986, 23: 17-29.

    Google Scholar 

  21. Sturmbauer C, Mark W, Dallinger R: Ecophysiology of aufwuchs-eating cichlids in Lake Tanganyika: niche separation by trophic specialization. Environ Biol Fishes. 1992, 35: 283-290. 10.1007/BF00001895.

    Article  Google Scholar 

  22. Wagner CE, McIntyre PB, Buels KS, Gilbert DM, Michel E: Diet predicts intestine length in Lake Tanganyika’s cichlid fishes. Funct Ecol. 2009, 23: 1122-1131. 10.1111/j.1365-2435.2009.01589.x.

    Article  Google Scholar 

  23. Takeuchi Y, Ochi H, Kohda M, Sinyinza D, Hori M: A 20-year census of a rocky littoral fish community in Lake Tanganyika. Ecol Freshw Fish. 2010, 19: 239-248. 10.1111/j.1600-0633.2010.00408.x.

    Article  Google Scholar 

  24. Kuwamura T: Overlapping territories of Psedosimochromis curvifrons males and other herbivorous cichlid fishes in Lake Tanganyika. Ecol Res. 1992, 7: 43-53. 10.1007/BF02348596.

    Article  Google Scholar 

  25. Yanagisawa Y, Nishida M: The social and mating system of the maternal mouthbrooder Tropheus moorii (Cichlidae) in Lake Tanganyika. Jpn J Ichthyol. 1991, 38: 271-282.

    Google Scholar 

  26. Kohda M: Coexistence of permanently territorial cichlids of the genus Petrochromis through male-mating attack. Environ Biol Fishes. 1998, 52: 231-242. 10.1023/A:1007381217885.

    Article  Google Scholar 

  27. Ochi H, Takeyama T, Yanagisawa Y: Increased energy investment in testes following territory acquisition in a maternal mouthbrooding cichlid. Ichthyol Res. 2009, 56: 227-231. 10.1007/s10228-008-0088-4.

    Article  Google Scholar 

  28. Kohda M, Takemon Y: Group foraging by the herbivorous cichlid fish, Petrochromis fasciolatus, in Lake Tanganyika. Ichthyol Res. 1996, 43: 55-63. 10.1007/BF02589608.

    Article  Google Scholar 

  29. Sturmbauer C, Fuchs C, Harb G, Damm E, Duftner N, Maderbacher M, Koch M, Koblmueller S: Abundance, distribution, and territory areas of rock-dwelling Lake Tanganyika cichlid fish species. Hydrobiologia. 2008, 615: 57-68. 10.1007/s10750-008-9557-z.

    Article  Google Scholar 

  30. Yamaoka K: A revision of the cichlid fish genus Petrochromis from Lake Tanganyika, with description of a new species. Jpn J Ichthyol. 1983, 30: 129-141.

    Google Scholar 

  31. Kohda M: Intra- and interspecific social organization among three herbivorous cichlid fishes in Lake Tanganyika. Jpn J Ichthyol. 1991, 38: 147-163.

    Google Scholar 

  32. Kohda M, Yanagisawa Y: Vertical distributions of two herbivorous cichlid fishes of the genus Tropheus in Lake Tanganyika, Africa. Ecol Freshw Fish. 1992, 1: 99-103. 10.1111/j.1600-0633.1992.tb00078.x.

    Article  Google Scholar 

  33. Kerschbaumer M, Mitteroecker P, Sturmbauer C: Evolution of body shape in sympatric versus non-sympatric Tropheus populations of Lake Tanganyika. Heredity (Edinb). 2014, 112: 89-98. 10.1038/hdy.2013.78.

    Article  CAS  Google Scholar 

  34. Pianka ER: Evolutionary Ecology. 1999, Addison-Wesley, Boston

    Google Scholar 

  35. Rossiter A: The cichlid fish assemblages of lake tanganyika: ecology, behaviour and evolution of its species flocks. Adv Ecol Res. 1995, 26: 187-252. 10.1016/S0065-2504(08)60066-5.

    Article  Google Scholar 

  36. Tanabe AS: Claident v0.1.2012.03.14, a Software Distributed by Author. [], [http://www.fifthdimension.jp]

  37. Tanabe AS, Toju H: Two new computational methods for universal DNA barcoding: a benchmark using barcode sequences of bacteria, archaea, animals, fungi, and land plants. PLoS One. 2013, 8: e76910-10.1371/journal.pone.0076910.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  38. Trewavas E: Tilapiine Fishes of the Genera Sarotherodon, Oreochromis and Danakilia. 1983, British Museum (Natural History), London, U.K

    Google Scholar 

  39. Brichard P: Pierre Brichard’s Book of Cichlids and All the Other Fishes of Lake Tanganyika. 1989, TFH Publications, Neptune City, NJ

    Google Scholar 

  40. Yanagisawa Y: Parental care in a monogamous mouthbrooding cichlid Xenotilapia flavipinnis in Lake Tanganyika. Ichthyol Res. 1986, 33: 249-261.

    Google Scholar 

  41. Morley J, Balshine S: Faithful fish: territory and mate defence favour monogamy in an African cichlid fish. Behav Ecol Sociobiol. 2002, 52: 326-331. 10.1007/s00265-002-0520-0.

    Article  Google Scholar 

  42. Morley J, Balshine S: Reproductive biology of Eretmodus cyanostictus, a cichlid fish from Lake Tanganyika. Environ Biol Fishes. 2003, 66: 169-179. 10.1023/A:1023610905675.

    Article  Google Scholar 

  43. Poll M: Poissons Cichlidae, volume III (5B). 1956, Institut Royal des Sciences Naturelles de Belgique, Brussels

    Google Scholar 

  44. Ochi H, Yanagisawa Y: Commensalism between cichlid fishes through differential tolerance of guarding parents toward intruders. J Fish Biol. 1998, 52: 985-996. 10.1111/j.1095-8649.1998.tb00598.x.

    Article  Google Scholar 

  45. Karino K: Depth-related differences in territory size and defense in the herbivorous cichlid, Neolamprologus moorii, in Lake Tanganyika. Ichthyol Res. 1998, 45: 89-94. 10.1007/BF02678579.

    Article  Google Scholar 

  46. Yamaoka K, Hori M, Kuwamura T: Interochromis, a new genus of the Tanganyikan cichlid fish. S Afr J Sci. 1998, 94: 381-386.

    Google Scholar 

  47. Takahashi T, Koblmüller S: A new species of Petrochromis (Perciformes: Cichlidae) from Lake Tanganyika. Ichthyol Res. 2014, 61: 252-264. 10.1007/s10228-014-0396-9.

    Article  Google Scholar 

  48. Takamura K: Interspecific relationship between two aufwuchs eaters Petrochromis polyodon and Tropheus moorei (Pisces: Cichlidae) of Lake Tanganyika, with a discussion on the evolution and functions of a symbiotic relationship. Physiol Ecol Japan. 1983, 20: 59-69.

    Google Scholar 

  49. Hori M: Structure of littoral fish communities organized by their feeding activities. Fish Communities in Lake Tanganyika. Edited by: Kawanabe H, Hori M, Nagoshi M. 1997, Kyoto University Press, Kyoto, 275-298.

    Google Scholar 

  50. Erickson A, Paul V, Alstyne K, Kwiatkowski L: Palatability of macroalgae that use different types of chemical defenses. J Chem Ecol. 2006, 32: 1883-1895. 10.1007/s10886-006-9116-x.

    Article  CAS  PubMed  Google Scholar 

  51. Dodds WK, Gudder DA: The ecology of Cladophora. J Phycol. 1992, 28: 415-427. 10.1111/j.0022-3646.1992.00415.x.

    Article  Google Scholar 

  52. Losos JB: Lizards in an Evolutionary Tree: Ecology and Adaptive Radiation of Anoles. 2009, Univ of California Press, Barkeley, CA

    Google Scholar 

  53. Rabosky DL, Glor RE: Equilibrium speciation dynamics in a model adaptive radiation of island lizards. Proc Natl Acad Sci U S A. 2010, 107: 22178-22183. 10.1073/pnas.1007606107.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  54. Losos JB, Glor RE, Jason JK, Nicholson K: Adaptation, speciation, and convergence: a hierarchical analysis of adaptive radiation in Caribbean anolis lizards. Ann Missouri Bot Gard. 2006, 93: 24-33. 10.3417/0026-6493(2006)93[24:ASACAH]2.0.CO;2.

    Article  Google Scholar 

  55. Losos JB, Leal M, Glor RE, de Queiroz K, Hertz PE, Rodriguez Schettino L, Chamizo Lara A, Jackman TR, Larson A: Niche lability in the evolution of a Caribbean lizard community. Nature. 2003, 424: 542-545. 10.1038/nature01814.

    Article  CAS  PubMed  Google Scholar 

  56. Bouton N, Seehausen O, van Alphen JJ: Resource partitioning among rock-dwelling haplochromines (Pisces: Cichlidae) from Lake Victoria. Ecol Freshw Fish. 1997, 6: 225-240. 10.1111/j.1600-0633.1997.tb00165.x.

    Article  Google Scholar 

  57. Kohda M: Intra- and interspecific territoriality of a temperate damselfish Eupomacentrus altus, (Teleostei: Pomacentridae). Physiol Ecol Japan. 1984, 21: 35-52.

    Google Scholar 

  58. Sato H, Murakami N: Reproductive isolation among cryptic species in the ectomycorrhizal genus Strobilomyces: population-level CAPS marker-based genetic analysis. Mol Phylogenet Evol. 2008, 48: 326-334. 10.1016/j.ympev.2008.01.033.

    Article  CAS  PubMed  Google Scholar 

  59. Toju H, Sato H, Yamamoto S, Kadowaki K, Tanabe AS, Yazawa S, Nishimura O, Agata K: How are plant and fungal communities linked to each other in belowground ecosystems? A massively parallel pyrosequencing analysis of the association specificity of root-associated fungi and their host plants. Ecol Evol. 2013, 3: 3112-3124. 10.1002/ece3.706.

    Article  PubMed Central  PubMed  Google Scholar 

  60. Toju H, Yamamoto S, Sato H, Tanabe AS, Gilbert GS, Kadowaki K: Community composition of root-associated fungi in a Quercus-dominated temperate forest: “codominance” of mycorrhizal and root-endophytic fungi. Ecol Evol. 2013, 3: 1281-1293. 10.1002/ece3.546.

    Article  PubMed Central  PubMed  Google Scholar 

  61. Hamady M, Walker JJ, Harris JK, Gold NJ, Knight R: Error-correcting barcoded primers for pyrosequencing hundreds of samples in multiplex. Nat Methods. 2008, 5: 235-237. 10.1038/nmeth.1184.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  62. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glöckner FO: The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013, 41: D590-D596. 10.1093/nar/gks1219.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  63. Wang Y, Qian PY: Conservative fragments in bacterial 16S rRNA genes and primer design for 16S ribosomal DNA amplicons in metagenomic studies. PLoS One. 2009, 4: e7401-10.1371/journal.pone.0007401.

    Article  PubMed Central  PubMed  Google Scholar 

  64. Cole JR, Wang Q, Fish JA, Chai B, McGarrell DM, Sun Y, Brown CT, Porras-Alfaro A, Kuske CR, Tiedje JM: Ribosomal Database Project: data and tools for high throughput rRNA analysis. Nucleic Acids Res. 2014, 42: D633-D642. 10.1093/nar/gkt1244.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  65. Klindworth A, Pruesse E, Schweer T, Peplies J, Quast C, Horn M, Glöckner FO: Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2013, 41: e1-10.1093/nar/gks808.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  66. Stiller JW, McClanahan AN: Phyto-specific 16S rDNA PCR primers for recovering algal and plant sequences from mixed samples. Mol Ecol Notes. 2005, 5: 1-3. 10.1111/j.1471-8286.2004.00805.x.

    Article  CAS  Google Scholar 

  67. Betournay S, Marsh AC, Donello N, Stiller JW: Selective recovery of microalgae from diverse habitats using “phyto-specific” 16S rDNA primers. J Phycol. 2007, 43: 609-613. 10.1111/j.1529-8817.2007.00350.x.

    Article  CAS  Google Scholar 

  68. Jullien N: AmplifX 1.6.0, CNRS, Aix-Marseille Université.. [] [in French]., [http://crn2m.univ-mrs.fr/pub/amplifx-dist]

  69. Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R: UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011, 27: 2194-2200. 10.1093/bioinformatics/btr381.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  70. Tanabe AS: Assams v0.1.2012.03.14, a Software Distributed by the Author. [], [http://www.fifthdimension.jp]

  71. Toju H, Sato H, Tanabe AS: Diversity and spatial structure of belowground plant–fungal symbiosis in a mixed subtropical forest of ectomycorrhizal and arbuscular mycorrhizal plants. PLoS One. 2014, 9: e86566-10.1371/journal.pone.0086566.

    Article  PubMed Central  PubMed  Google Scholar 

  72. Sommer D, Delcher A, Salzberg S, Pop M: Minimus: a fast, lightweight genome assembler. BMC Bioinformatics. 2007, 8: 64-10.1186/1471-2105-8-64.

    Article  PubMed Central  PubMed  Google Scholar 

  73. Oksanen JF, Blanchet G, Kindt R, Legendre P, Minchin PR, O’Hara RB, Simpson GL, Solymos P, Stevens MHH, Wagner H: Vegan: community ecology package. R package version 2.0-9. [], [http://CRAN.R-project.org/package=vegan]

  74. Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL: BLAST+: architecture and applications. BMC Bioinformatics. 2009, 10: 421-10.1186/1471-2105-10-421.

    Article  PubMed Central  PubMed  Google Scholar 

  75. Huson DH, Auch AF, Qi J, Schuster SC: MEGAN analysis of metagenomic data. Genome Res. 2007, 17: 377-386. 10.1101/gr.5969107.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  76. Bengtsson J, Eriksson KM, Hartmann M, Wang Z, Shenoy B, Grelet GA, Abarenkov K, Petri A, Alm Rosenblad M, Nilsson RH: Metaxa: a software tool for automated detection and discrimination among ribosomal small subunit (12S/16S/18S) sequences of archaea, bacteria, eukaryotes, mitochondria, and chloroplasts in metagenomes and environmental sequencing datasets. Antonie Van Leeuwenhoek. 2011, 100: 471-475. 10.1007/s10482-011-9598-6.

    Article  PubMed  Google Scholar 

  77. Lundberg DS, Yourstone S, Mieczkowski P, Jones CD, Dangl JL: Practical innovations for high-throughput amplicon sequencing. Nat Methods. 2013, 10: 999-1002. 10.1038/nmeth.2634.

    Article  CAS  PubMed  Google Scholar 

  78. Lundberg DS, Lebeis SL, Paredes SH, Yourstone S, Gehring J, Malfatti S, Tremblay J, Engelbrektson A, Kunin V, del Rio TG, Edgar RC, Eickhorst T, Ley RE, Hugenholtz P, Tringe SG, Dangl JL: Defining the core Arabidopsis thaliana root microbiome. Nature. 2012, 488: 86-90. 10.1038/nature11237.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  79. King RA, Read DS, Traugott M, Symondson WO: Molecular analysis of predation: a review of best practice for DNA-based approaches. Mol Ecol. 2008, 17: 947-963. 10.1111/j.1365-294X.2007.03613.x.

    Article  CAS  PubMed  Google Scholar 

  80. Dormann CF, Fruend J, Bluethgen N, Gruber B: Introducing the bipartite package: analysing ecological networks. R news. 2008, 8: 8-11.

    Google Scholar 

  81. Blüthgen N, Menzel F, Blüthgen N: Measuring specialization in species interaction networks. BMC Ecol. 2006, 6: 9-10.1186/1472-6785-6-9.

    Article  PubMed Central  PubMed  Google Scholar 

  82. Blüthgen N, Menzel F, Hovestadt T, Fiala B, Blüthgen N: Specialization, constraints, and conflicting interests in mutualistic networks. Curr Biol. 2007, 17: 341-346. 10.1016/j.cub.2006.12.039.

    Article  PubMed  Google Scholar 

  83. Vázquez DP, Melián CJ, Williams NM, Blüthgen N, Krasnov BR, Poulin R: Species abundance and asymmetric interaction strength in ecological networks. Oikos. 2007, 116: 1120-1127. 10.1111/j.0030-1299.2007.15828.x.

    Article  Google Scholar 

  84. Dormann CF, Fründ J, Blüthgen N, Gruber B: Indices, graphs and null models: analyzing bipartite ecological networks. Open Ecol J. 2009, 2: 7-24. 10.2174/1874213000902010007.

    Article  Google Scholar 

  85. Anderson MJ: A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001, 26: 32-46.

    Google Scholar 

  86. Kuczynski J, Liu Z, Lozupone C, McDonald D, Fierer N, Knight R: Microbial community resemblance methods differ in their ability to detect biologically relevant patterns. Nat Methods. 2010, 7: 813-819. 10.1038/nmeth.1499.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  87. Hothorn T, Bretz F, Westfall P: Simultaneous inference in general parametric models. Biom J. 2008, 50: 346-363. 10.1002/bimj.200810425.

    Article  PubMed  Google Scholar 

  88. R: a Language and Environment for Statistical Computing. 2013, R Foundation for Statistical Computing, Vienna

Download references

Acknowledgements

We are grateful to members of the Maneno Team (Tanganyika Research Project Team) and staff of LTRU, Mpulungu, Zambia, for their support. We also thank two anonymous referees for critically reading the manuscript and providing helpful suggestions. This study is supported by Grant-in-Aid for Young Scientists (B, nos 22770024, 25840159) for HH, Grant-in-Aid for Scientific Research (B, no. 22405010) for MK, JSPS KAKENHI Grant (No. 26711026) and the Funding Program for Next Generation World-Leading Researchers of Cabinet Office, the Government of Japan (GS014) for HT. The research presented here was conducted under permission of fish research in Lake Tanganyika from the Zambian Ministry of Agriculture, Food, and Fisheries and complied with the current law in Zambia.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Hiroki Hata.

Additional information

Competing interests

The authors declare that they have no competing interests.

Electronic supplementary material

12915_2014_90_MOESM1_ESM.xlsx

Additional file 1: Data S1.: Summary of pyrosequenced SSU rRNA gene sequence reads that passed quality filtering. OTUs observed in the algal farms and stomachs of territorial herbivorous cichlid fishes in Lake Tanganyika. (XLSX 1 MB)

Additional file 2: Table S1.: Summary of phototrophic OTUs detected in our pyrosequencing analysis. (PDF 38 KB)

Additional file 3: Data S2.: Table of rarefied phototrophic OTUs used for analyses. (XLSX 90 KB)

12915_2014_90_MOESM4_ESM.pdf

Additional file 4: Table S2.: Result of Adonis on distance matrices of algal composition in algal farms defended by territorial cichlids, and the stomach contents of herbivorous cichlids. Distances were calculated by Canberra distance index and Bray-Curtis dissimilarity index. DF, degree of freedom; NS, not significant. (PDF 11 KB)

12915_2014_90_MOESM5_ESM.pdf

Additional file 5: Table S3.: Summary of GLM testing for the effect of cichlid species on their habitat depth. (PDF 5 KB)

12915_2014_90_MOESM6_ESM.pdf

Additional file 6: Table S4.: Summary of GLMM testing for the effects of phototrophic OTUs on their habitat depths. Cichlid individual, from whose territory each OTU was collected, was included as a random effect. (PDF 11 KB)

Additional file 7: Figure S1.: Occurrence frequencies of the dominant phototrophic OTUs in various depths. (PDF 129 KB)

12915_2014_90_MOESM8_ESM.pdf

Additional file 8: Table S5.: Summary of GLMs testing for the effect of difference in tribe, ecomorph, and habitat depth on the algal farm and stomach contents similarity. Ecomorph combination (pairs of the same ecomorph or those of different ecomorphs) and tribe combination (pairs of the same tribe or those of different tribes) are included as two fixed effects. A, Canberra index; B, Bray-Curtis dissimilarity. NS, not significant. (PDF 13 KB)

12915_2014_90_MOESM9_ESM.pdf

Additional file 9: Figure S2.: Coverage of primer pairs for algal and cyanobacteria SSU rDNA. Forward and reverse primers: CYA_ALG_F and CYA_ALG_R (A, B); PSf (5′-GGG ATT AGA TAC CCC WGT AGT CCT-3′) and Ur (5′-ACG GYT ACC TTG TTA CGA CTT-3′) from [66] (C, D); and PSf and a universal bacteria primer (5′-TAC GGY TAC CTT GTT ACG ACT T-3′) from [67] (E, F). The percentages of sequences amplified by Probe Match of RDP Release 11.2 are shown in A, C and E, and those amplified by TestPrime 1.0 of SILVA are shown in B, D and F. Zero, one, or two nucleotide mismatches between the target primer and database sequences were allowed in each analysis. Numbers in parentheses indicate the numbers of eligible sequences that attempted to match with the primer pair. (PDF 99 KB)

12915_2014_90_MOESM10_ESM.pdf

Additional file 10: Figure S3.: Rarefaction curves of OTUs based on the number of reads. A) algal farms of cichlid fishes, B) stomach contents and C) rarefaction curve of OTUs based on the number of samples for algal farms of cichlid fishes (blue curve) and stomach contents (red curve). Species abbreviations are listed in Table 1. The shaded area represents the standard deviation obtained from 100 shuffles of sample-ID order. (PDF 394 KB)

12915_2014_90_MOESM11_ESM.pdf

Additional file 11: Figure S4.: Principal coordinate analysis among algal farms of herbivorous cichlids (A) and their stomach contents (B). Dissimilarity was calculated using the Bray-Curtis dissimilarity index. Axes represent the first two principal coordinates maximising the variance in the data (PC 1 and PC 2). The percentage of the total variance explained by each axis estimated under the broken stick model is shown in parentheses. Species abbreviations are listed in Table 1. Orange-, green-, dark red-, grey-, and yellow-coloured species indicate grazers, browsers, scrapers, scoopers, and biters, respectively. Out = periphyton outside the cichlid territories. (PDF 49 KB)

12915_2014_90_MOESM12_ESM.pdf

Additional file 12: Figure S5.: Pairwise distance-contrast plots between phylogenetic distance of fishes versus dissimilarity in phototroph compositions. A) The correlation between the phylogenetic distance of fish species versus dissimilarity in the phototroph compositions of algal farms defended by the herbivorous cichlids, and B) pairwise plots between the phylogenetic distance versus dissimilarity in phototroph composition of stomach contents of the herbivorous cichlids. Circles and squares indicate species-pairs of the same feeding ecomorph and of different ecomorphs, respectively. Closed and open symbols indicate species-pairs of the same tribe and of different tribes, respectively. Dissimilarity was calculated using the Bray-Curtis dissimilarity index. Species abbreviations are listed in Table 1. (PDF 64 KB)

12915_2014_90_MOESM13_ESM.pdf

Additional file 13: Figure S6.: Pairwise distance-contrast plots between difference in habitat depth versus dissimilarity in phototroph compositions. A) The correlation between the difference in habitat depth of fish species versus the dissimilarity in the phototroph compositions of algal farms defended by the herbivorous cichlids, and B) pairwise plots between the difference in habitat depth versus the dissimilarity in phototroph composition of stomach contents of the herbivorous cichlids. Circles and squares indicate species-pairs of the same feeding ecomorph and of different ecomorphs, respectively. Closed and open symbols indicate species-pairs of the same tribe and of different tribes, respectively. Dissimilarity was calculated using the Bray-Curtis dissimilarity index. Species abbreviations are listed in Table 1. (PDF 59 KB)

Authors’ original submitted files for images

Rights and permissions

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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hata, H., Tanabe, A.S., Yamamoto, S. et al. Diet disparity among sympatric herbivorous cichlids in the same ecomorphs in Lake Tanganyika: amplicon pyrosequences on algal farms and stomach contents. BMC Biol 12, 90 (2014). https://0-doi-org.brum.beds.ac.uk/10.1186/s12915-014-0090-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12915-014-0090-4

Keywords