Effects of predation pressure and resource use on morphological divergence in omnivorous prey fish
© Scharnweber et al.; licensee BioMed Central Ltd. 2013
Received: 22 March 2013
Accepted: 19 June 2013
Published: 27 June 2013
Body shape is one of the most variable traits of organisms and responds to a broad array of local selective forces. In freshwater fish, divergent body shapes within single species have been repeatedly observed along the littoral-pelagic axes of lakes, where the structural complexity of near shore habitats provides a more diverse set of resources compared to the open-water zones. It remains poorly understood whether similar resource-driven polymorphism occurs among lakes that vary in structural complexity and predation pressure, and whether this variation is heritable. Here, we analyzed body shape in four populations of omnivorous roach (Rutilus rutilus) inhabiting shallow lakes. We tested the relationship between body shape, gradients of resources, predation pressure, and, in a subset of two lakes, diet composition. We used genome scans of 331 polymorphic AFLP markers to test whether there was a heritable component to the observed morphological diversification.
Body shape differed among lakes and was significantly correlated to differences in predation pressure. Roach from the lake with highest predation pressure were most divergent from the average body shape of all populations, characterized by a more streamlined body and caudally inserted dorsal fins; features that facilitate predator escape. Surprisingly, diet composition was not associated with morphology. AFLP analysis revealed weak genetic differentiation among lakes and no isolation by distance (IBD). Outlier analysis detected three loci under positive selection with differing frequencies in the four populations. General linear models did not support an association of lake-specific genotypes with morphological variation.
Body shape was divergent among lakes, suggesting that processes previously reported from within single lakes may also be operating at the scale of whole lakes. We found no evidence for body shape being heritable, although sample size was small in these natural populations. Rather than habitat structure and diet, we conclude that predation had a stronger effect on the prevalence of local morphotypes. A variable morphotype facilitating the efficient uptake of a variety of spatially and temporarily scattered resources seems to be favored in these small aquatic systems.
Populations can exhibit divergent body morphologies in response to environmental cues and these differences may have a significant influence on the performance of individuals . A common form of intraspecific divergence in body morphology is resource polymorphism, where morphological variation is associated with segregation in habitat and diet [2, 3]. In freshwater fish, divergence of morphs has been frequently shown along the pelagic-littoral habitat axis of lakes (i.e., from the near shore to open-water habitats) . Foraging in the open water is associated with a high search rate for widely distributed and conspicuous planktonic prey. A more streamlined body occurs in these morphs to facilitate a high attack speed. In contrast, foraging in the structurally complex littoral zones results in a lower search rate for more cryptic benthic prey. Thus, a deeper body supports higher manoeuvrability [4, 5]. Fish species with coexisting pelagic and benthic morphs include three-spined stickleback (Gasterosteus aculeatus) [6–8], Eurasian perch (Perca fluviatilis) [9–11] and Arctic charr (Salvelinus alpinus) [12–14].
Predation pressure is another potentially selective force affecting the morphology of organisms. “Inducible defenses”  have been intensively studied in cladocerans [16, 17] and mussels  and are known from a range of fish species. In the presence of predators, Crucian carp (Carassius carassius) develop a deeper body that increases the handling time by predators . Pumpkinseed sunfish (Lepomis gibbosus) increase body depth and dorsal spine length when stimulated by permanent predation cues of walleye (Sander vitreus) . Similar predator-induced responses have been described for perch .
Phenotypic variation can result from phenotypic plasticity [22, 23] or from heritable adaptations to local selective forces [24–27]. To distinguish between phenotypically plastic and inherited variation, population genomics are widely used to test for the genetic basis of adaptive divergence . Genome scans are an appropriate method to compare natural populations of non-model organisms that are nonetheless ecologically important, allowing simultaneous screening of many loci throughout the genome and the identification of loci that are putatively under selection . As an example, genome scans were used to identify loci associated with morphological variation in European minnows (Phoxinus phoxinus) living in different habitat types .
In the present study, we focused on a fish species for which we hypothesized substantial morphological variation due to a variety of locally specific selective forces. The omnivorous cyprinid roach (Rutilus rutilus) occurs in many lakes and rivers of the European temperate zone . Resource polymorphism in roach was recently demonstrated in two lakes : individuals inhabiting the pelagic zone were streamlined and fed on zooplankton, whereas those living in the littoral zone had a deeper body and were feeding on benthic invertebrates . The dominant feature of near-shore littoral habitats is their structural complexity, created by aquatic macrophytes or woody debris that accumulates in shallow water zones. Structural complexity enhances the biomass of prey organisms for omnivorous fish, as the physical structures provide a substrate for epiphytic algae, a primary food source for the invertebrate prey [33–36].
Much previous work has focused on structural complexity within lakes, but complexity also varies among lakes, with high complexity being a typical feature of shallow lakes in particular . Shallow lakes with low phytoplankton abundance (“clear” lakes) can contain extended zones of structural complexity because the clear water allows light penetration and facilitates the growth of submerged macrophytes. In contrast, shallow lakes with an abundant phytoplankton community (“turbid” lakes) can have low structural complexity because the lack of light prevents macrophyte establishment. As a result, the differential selective forces operating in littoral and pelagic habitats within single lakes and inducing intralacustrine (i.e. within lake) divergence may also operate in an interlacustrine (i.e. among lake) comparison of populations from differently structured shallow lakes.
Here we examined morphology (geometric morphometrics), diet (gut content and stable isotopes) and genetic variation (AFLP) of roach in four shallow lakes in order to assess ecomorphological divergence and local adaptation among lakes. Two lakes were clear and had extensive macrophyte cover while two were turbid and had almost no macrophyte growth. We predicted that roach in clear lakes would feed mainly on benthic macroinvertebrates and have a correspondingly deeper body shape, and that roach from turbid lakes would be more planktivorous and exhibit a more slender shape. We also expected that morphological divergence would be further enhanced by exposure to different predators. The four lakes exhibited a gradient of predation pressure, measured as the proportion of piscivorous fish biomass in total fish assemblages, thus providing a basis to test for an association of body morphology with relative predation pressure. Finally, we tested whether the roach populations differed genetically between the lakes, and whether genetic markers were associated with the interlacustrine morphological divergence.
Pairwise distances between the group means in shape space
Mahalanobis distance D
Genetic (AFLP) analysis
Genetic diversity of roach from four lakes of different structural complexity
% polymorphic loci
95% credibility interval of hs
Σ = 183
Total = 77.6
Hs = 0.271
Outlier analysis using BAYESCAN classified three loci (loci 31, 48 and 155; 0.91% of polymorphic loci examined) as outliers under divergent selection (probability P > 0.97 that the locus-specific component of FST, α, is different from zero, corresponding to Posterior Odds > 32, estimated false discovery rate = 0.0084). Probabilities for the remaining 328 loci ranged from 0.064 to 0.572. The three outlier loci were present in all four populations, but with differing frequencies. Locus 48 dominated in Gollinsee (97%) and Globsowsee (88%), but was rarer in Schulzensee (28%) and Döllnsee (15%). In contrast, locus 155 dominated in Döllnsee (95%) and Schulzensee (78%), but was rare in Globsowsee (19%) and Gollinsee (11%). Locus 31 was highly frequent in Döllnsee (100%), Gollinsee (97%) and Schulzensee (91%), but less frequent in Globsowsee (31%). Linkage disequilibrium (LD; ARLEQUIN) was estimated for all possible pairs of the three outlier loci, but among the 12 (4 lakes × 3 pairs) pair-wise tests, only one locus pair in Gollinsee showed a significant linkage (P = 0.009, all other P > 0.23).
Association of shape with predation pressure and genetic markers
Principal coordinate analysis (PCoA) of the 331 polymorphic AFLP loci aggregated 58.8% of variation along PCoA axis 1 (26.2%), PCoA2 (17.7%) and PCoA3 (14.9%). CV1 of shape analysis was significantly smaller in Globsowsee, Schulzensee and Gollinsee than in Döllnsee (GLM, F7,85 = 75.0, adjusted R2 = 0.85, P < 0.0001, coefficients of lake contrasts: t < −7.1, P < 0.0001), but aggregated genotypes (PCoA1) did not contribute to shape variations (coefficient for PCoA1, t = 0.36, P = 0.72). Interactions between PCoA1 and the lakes were likewise not significant (coefficients of interactions, t > 0.22, P > 0.83). Similarly, PCoA2 of genotypes and its interactions with lakes were not significant predictors of shape (CV1) (coefficient of PCoA2, t = −0.36, P = 0.34), but CV1 again significantly differed between Döllnsee and the other three lakes in this analysis (GLM, F7,85 = 78.5, adjusted R2 = 0.86, P < 0.0001; coefficients of lake contrasts, t < −7.1, P < 0.0001). General linear models on CV1 as the dependent dominant shape variable and lake and each of the three outlier loci as predictors confirmed that the genotype had no effect on shape of fish, because the effect of the outlier loci was not significant in any case (GLMs, F7,85 > 75.0, adj. R2 > 0.85, P < 0.0001; coefficients of the single outlier loci: locus 31, t = −0.27, P = 0.79; locus 48, t = 0.58, P = 0.56; locus 155, t = 0.43, P = 0.67).
Roach diet and association with morphological data
PLS analysis of shape as the first block and gut contents of roach from Gollinsee and Schulzensee as the second block revealed no correlation of shape with diet composition (RV = 0.0475; permutation test with 1000 randomization runs: P = 0.66). Similarly, PLS with shape as the first block and isotopic ratios as second block did not show a significant correlation (RV = 0.0205; permutation test with 1000 randomization runs: P = 0.70).
Body morphologies differed in all four studied roach populations in shallow lakes, in particular with respect to overall body form and anterior-posterior location of the fins. Our analyses suggest that these morphological differences reflect a gradient in exposure to predation risk among the lakes, as shape and predation pressure were significantly correlated. Contrary to our predictions, fish morphology was not correlated with resource use. Therefore, we suggest that the strong morphological divergence among lakes in roach was primarily driven by differences in predation risk. The four populations exhibited small genetic differences. However, we did not find evidence that genetic differentiation and morphological variation among lakes were associated; suggesting that shape variation in roach did not have a heritable component detectable in our analysis.
Roach body shapes were most distinct in Döllnsee, the lake characterized by highest predator abundance. There are several examples for induced morphological defenses in fish, including a plastic response to achieve a size refuge against predators limited by gape-size [19–21]. Here, we observed a more slender body and a more caudal position of the dorsal fin. In contrast to morphological changes that prevent handling of prey by predators [19–21] we consider this body form to be advantageous for flight. To escape a predator attack, most fish show a similar fast-start behavior of sudden, high-energy swimming bursts, called a “C-start” [38–40]. During a C-start, fish bend into a “C” shape (stage 1) and then produce a propulsive stroke of the caudal region in the opposite direction (stage 2). This is facilitated by a slender body form and a relatively shallow anterior body and head region, as these contribute minimally to drag [41–43]. Thrust can be increased if the dorsal fin is in a caudal position ; however, Domenici et al.  demonstrated higher speed and acceleration in Crucian carp of deeper body forms, attributable to a high percentage of muscle mass. Therefore, there seems to be a trade-off between a streamlined body that reduces drag and a deeper body of high muscle mass that increases burst.
A previous study compared morphology of one-year-old roach raised under experimental conditions in the presence and absence of pike (Esox) predators . The main differences were changes of dorsal fin displacement and morphological responses similar to our findings. The authors suggested that these morphological changes affect swimming speed and manoeuvrability in roach . Contrary to our study, they detected a wider anal fin in roach exposed to predators. Adaptations to faster swimming seem advantageous, as roach were found to increase activity and even show the tactic of jumping out of the water when threatened by a predator . Comparable adaptations to faster predator escape as a response to high predator abundances were also found in Western mosquitofish (Gambusia affinis)  and in Galaxias platei, where caudal regions were enlarged to enhance burst performance. A recent review  concluded that most of the observed predator-induced morphological defenses are purely a by-product of swimming activity. Prey swim less when predators are present, resulting in lower feeding rates or decreased metabolism, leading to reduced growth. We cannot exclude that variation in roach morphology between the lakes reflects differences in swimming activity (in response to differing predation risk); however, roach from Döllnsee, where predation risk was highest, had morphologies that appear better adapted to rapid and extended swimming, which contrasts the conclusion that predation risk always induces reduced swimming activity of prey.
Body morphology was not associated with structural complexity in the lakes created by macrophyte coverage, with little evidence for an association of morphology and diet composition in this omnivorous species. Gut content and stable isotopes both suggested that benthic invertebrates dominated the diet of roach in the two lakes analyzed, whereas shape data significantly differed between these populations. Similar studies of fishes mainly focused on intralacustrine (within-lake) morphological divergence and resource polymorphism, with specialization on benthic food usually associated with a preference for littoral habitats with high structural complexity. Intralacustrine resource polymorphism driven by intraspecific competition has been studied in roach of two Swedish lakes . Benthivorous individuals from the littoral zone had a deeper body compared to planktivorous individuals from the pelagic zone. Similarly, pumpkinseed sunfish produce molluscivorous and a planktivorous morphotypes to decrease intraspecific competition [49, 50]. Ours is the first study of interlacustrine (among-lake) divergence of body morphology in roach of which we are aware, but similar morphological divergence between lakes has been reported for populations of three-spined sticklebacks in Alaska, USA . Very shallow lakes, where most of the stickleback habitat is structurally complex, were inhabited by individuals of the benthic form, characterized by a deeper body that favors maneuverability, whereas the limnetic form with a narrower body that reduces drag during prolonged swimming was found in deeper lakes with large areas of open waters .
A possible explanation for the lack of correlation between resource use and body shape here could be the spatial dimension of the systems studied. The four lakes we studied were small (< 0.2 km2) and relatively shallow. Roach may be able to easily exploit the resources from all habitats because near shore and offshore sites are usually less than 100 m apart. Furthermore, seasonal and annual changes in resource abundance require flexibility in resource use, perhaps counteracting any selective advantage gained by specialization. A high diet flexibility of roach is supported by the fact that guts of roach in the lake with almost no macrophytes (Gollinsee) contained even more benthic prey than those from fish in the lake with intermediate macrophyte coverage (Schulzensee). A more balanced diet composition with resources from both littoral and pelagic origin was obtained by applying the mixing model based on stable isotope analysis. Therefore, in systems where resources are spatially dispersed and temporally variable, selection may favor morphological shapes that support omnivory rather than discrete morphotypes specialized on a particular prey and habitat type.
Many previous studies found morphological changes to be phenotypically plastic [19, 52], although a few studies have demonstrated a heritable component [43, 53, 54]. In our dataset, we did not find any statistical association of shape differences with aggregated or outlier-specific genotypes among the lakes. However, the three outlier loci occurred in all populations, although with differing frequencies, suggesting that the allele frequencies of these loci are not the result of neutral genetic drift. There is some inherent uncertainty in AFLP studies, such as false scoring of alleles due to differences in PCR efficiency (false positives) or peak overlap (false negatives) . Furthermore, reduced precision during sample preparation due to sampling mix-up, pipetting or contamination may cause genotyping errors . Nevertheless, we have no indication that false positives had a strong influence on our results, because the reproducibility was 98% within the randomly chosen repeated 17.4% of samples.
The proportion of outlier loci (0.91%) fits within the range of reported outliers of former studies (0.4 – 24.5%) . In a comparison of three available statistical programs, BAYESCAN (the method we used here) showed the lowest rates of false positives , and hence we consider the three loci found in our study to be reliable (see also the very low false discovery rate estimated by BAYESCAN for our data). AFLP outlier loci have been linked to morphological diversification in other studies of fish, including European minnows in lake and stream habitats  or in sympatric lake whitefish (Coregonus clupeaformis) ecotypes . A similarly strong association between morphotype and genotype could not be confirmed for roach in our study. It is possible that our AFLP markers simply failed to include genomic regions that are associated with the shape variations between the lakes. Alternatively, the lack of association between genotype and phenotype may suggest a plastic response to predation pressure. Compared to fixed genetic controls, plastic responses to environmental pressures are advantageous when predatory regimes are variable . Strong interannual variability of predation pressure is likely in the communities analyzed here because these small and shallow lakes are characterized by frequent fish kills during strong winters that can substantially reduce predator abundance .
Overall, the extent of genetic divergence between populations is determined through a mixture of adaptive and neutral processes . Local adaptation is achieved through environmental filtering of genotypes , whereas neutral processes include drift-migration equilibrium wherein landscape isolation plays a role in reducing gene flow (e.g., isolation by distance, IBD) . We did not find evidence for a dominance of neutral processes and IBD in our data. The geographical distance between the four lakes is very limited (<35 km), and the two lakes situated only 4 km apart (Gollinsee and Döllnsee) exhibited intermediate genetic divergence (θ = 0.0266). Thus, we conclude that interlacustrine morphological diversification of the four roach populations is mainly led through processes related to local environments of each lake (in particular predator abundance) rather than neutral process due to geographic isolation between the lakes.
We observed divergent morphologies in omnivorous fish populations among four shallow lakes. Our results support previous studies suggesting that predation risk may play a more important role than habitat and diet in shaping body morphology of this prey species. In the case of roach inhabiting small and shallow lakes, we consider these morphological changes advantageous for fast escapes and thus directly cued by predation. Furthermore, we did not find evidence that the observed morphological divergence had a heritable component.
Study areas and sampling
Characteristics of the four lakes included in the study
Mean depth (m)
No. fish species
Coverage of submerged macrophytes (%)
Relative abundance of piscivorous fish (%)
Fish were caught using Nordic multi-mesh gill nets (length 30 m, height 1.5 m; 12 mesh sizes from 5–55 mm; Lundgrens Fiskredskapsfrabrik AB) or electrofishing with 200–400 V DC current (4–8 electrofishing aggregates EFG 4000, EFGI 4000; Bretschneider Spezialelektronik, Breitenbrunn, Germany) with anodic handnets (4–4.5 m long, 40 cm diameter, 6 mm mesh size) during summer and early autumn 2010. After capture, fish were measured (mm total length, TL) and weighed (g wet mass, wm) and subsequently stored on ice. In total, 185 roach (69 from Gollinsee, 70 from Schulzensee, 25 from Globsowsee and 21 from Döllnsee) were processed. Due to methodological problems with some individuals, the number of fish per analysis was slightly lower.
Body morphology of individual roach was examined using a landmark-based geometric morphometric method . Fish were photographed from the left side in a standardized way using a Nikon DX40 and a 28 mm lens, which represents a normal lens when considering the sensor size of the camera. Optical distortions were minimized by using a normal lens, a moderate distance to the specimen, and an aperture close to the critical aperture. Individuals were placed in a bowl of fine white gravel to obtain a straight position and the fins were stretched out to make the fin base visible. A few specimens with injuries on the left side were photographed from the right side and then the photographs were digitally flipped horizontally prior to analysis. After taking the photo, individuals were cut open to determine the sex. The shape data of the two sexes did not differ significantly (Mahalanobis D = 1.0761; P = 0.53), therefore sexes were pooled for further analyses. All roach were at least two summers old, an age where roach in lakes at this latitude are sexually mature .
Digital photographs were transferred to TPSdig2 (all TPS- programs http://life.bio.sunysb.edu/morph) and 14 landmarks were digitized. To diminish measurement error, digitizing was always performed by the same person. Haphazardly chosen individuals were repeated and results were compared to determine reliability of the results. We checked for outliers with the “Find outliers” function of the software MorphoJ. TPSSmall was used to determine whether the amount of variation in shape in our data set was small enough to permit statistical analyses in the linear tangent space approximated to the non-linear Kendall's shape space. For all further analyses, MorphoJ was used. Landmark configurations were aligned by Procrustes Superimposition  to minimize effects of translation, rotation and scaling. Shape data were size-corrected using a regression of shape (i.e. Procrustes coordinates) on size (i.e. log centroid size) for each lake separately. Residuals obtained from this regression were used for all further analyses. A Discriminant Function analysis (DFA) and a Canonical Variate analysis (CVA) were used to assess significance of shape differences between groups. A pairwise comparison between all four lakes was conducted (Globsowsee vs. Döllnses vs. Gollinsee vs. Schulzensee).
Genetic (AFLP) analysis
Primer combinations of selective amplification
Number of loci
MseI + CAG
EcoRI + ACT + 6FAM™
MseI + CAG
EcoRI + ACA + VIC®
MseI + CAG
EcoRI + AAG + NED™
MseI + CTA
EcoRI + AAG + PET®
Signal processing and binning was conducted using Genemapper™ V.4.1 (Applied Biosystems). Presence (1) or absence (0) of fragments were scored between 50 and 450 bp using a peak intensity threshold > 100 relative florescent units. A bin width of 1 bp and a max peak of 1.5 bp yielded optimization and fragments were scored between 50 to 450 bp. A complete repetition of a random 17.4% of the samples yielded 97.8% reproducibility.
The input data files for statistical analyses were prepared by AFLPDAT . AFLP-SURV  was used to calculate the frequency of polymorphic fragments. To obtain information on genetic structure within and among all four populations, genetic differentiation was calculated by using a Bayesian approach (θ) as implemented in HICKORY , using 20,000 burn-in iterations and 100,000 final iterations. AFLP loci potentially under selection (outlier loci) were identified using the BAYESCAN software  by assuming that loci influenced by directional selection show larger genetic differentiation than neutral loci, and loci that have been subjected to balancing selection show a lower genetic differentiation. This software uses a Bayesian likelihood method by reversible jump Monte Carlo Markov Chains. From FST calculations, posterior probabilities of two models are estimated: one model including selection on a certain locus (the locus-specific component of FST, α, is different from zero), the other one not. The Posterior Odds (ratio between these probabilities) provides a detection level for a locus to be under selection. We retained loci with P > 0.97 (Posterior Odds > 32) as outliers presumably under directional selection. These calculations were run only with the 331 polymorphic AFLP loci. We further used ARLEQUIN 3.5 to test for pairwise linkage disequilibrium (LD) of the outlier loci, using 1000 steps in the Markov chain and a dememorization of 1000 steps.
To approximate the predation pressure by piscivores on roach for each lake, lake-specific abundances of piscivorous fish were multiplied with the individual, size-dependent predation risk of roach from that lake. For Globsowsee, Schulzensee and Gollinsee, relative abundances of piscivorous fish were obtained from standardized gill net catches (6–8 benthic multimesh gill nets set over night in autumn). We included perch (>15 cm TL) (see  for a definition of piscivory in perch) and pike as piscivorous fish species. In the fourth lake (Döllnsee), destructive gill-net fishing could not be conducted because the lake is equipped with a 3D-telemetry system for in-situ studies of fish behavior. Therefore, we used information on fish assemblage composition from gill-net samplings conducted in previous years [73, 74]. The proportion of piscivores in Döllnsee calculated from earlier samplings is conservative because the lake has been heavily stocked with pike during the last years, and hence roach presumably experience an even greater predation pressure. The lake-specific relative abundances of piscivorous fish were ranked from highest (predator abundance = 4) to lowest (predator abundance = 1) (Figure 2). Predation risk of prey fish is size-dependent due to gape limitation of piscivores . Therefore, predation risk of individual roach was ranked according to total length, with TL ≥ 25 cm: predation risk = 0; TL ≥ 20-25 cm: predation risk = 1; TL ≥ 15-20 cm: predation risk = 2; TL ≥ 10-15 cm: predation risk = 3; TL < 10 cm: predation risk = 4 (Figure 2). The individual predation risk was ultimately multiplied with the rank among lakes with respect to relative abundances of piscivores. Therefore, the smallest roach in Döllnsee had the highest relative risk (4 × 4 = 16), whereas all roach > 25 cm TL had a risk of zero, independent of lake of origin. For all subsequent analyses, the obtained product for predation risk ranging from 0 to 16 was log2(x + 1)-transformed to achieve a near-linear distribution.
Gut content and stable isotope analysis
Detailed data on roach diet composition was obtained from roach in Gollinsee (N = 42) and Schulzensee (N = 42). In the laboratory, roach guts were removed and stored in 5% formaldehyde for subsequent GC. Individuals with empty stomachs were removed and the anterior third of the gut was examined under a stereo microscope and the volume proportion (equivalent to area proportion at uniform width) of each prey category observed in the sample was estimated to the nearest 10%, adapted from Windell . The observed prey items were grouped into two categories, namely benthic prey (trichoptera, gastropods, isopods, chironomid larvae, bryozoans, algae and detritus) and pelagic prey (cladocerans, copepods, water mites, ostracods, rotifers, and chaoborid larvae). Eleven fish with empty guts were removed prior to analysis. A non-parametric Mann–Whitney U test was conducted to compare the percentage of benthic and pelagic prey between the two lakes.
For SIA of carbon (13C/ 12C) and nitrogen (15 N/ 14 N) from roach (N = 46 from Schulzensee and N = 47 from Gollinsee, including fish with empty guts), a small sample of dorsal muscle tissue was excised from each fish. To obtain baseline values of potential diet groups for the omnivorous roach, macroinvertebrate samples for SIA were collected in April and June from various sites in the littoral zone using a sweep net. Zooplankton samples for SIA were taken by hauling zooplankton nets (55 μm and 100 μm) vertically through the entire water column several times during summer. Samples were taken to the laboratory and transferred to clean tap water to allow the animals to void their guts overnight. Macroinvertebrates were then sorted into taxon groups and zooplankton samples were filtered (100 μm filters). All samples were dried at 60°C for 48 hours, then ground to a fine powder using a mortar and pestle and material of about 0.5 mg dry mass was loaded into tin cups. Carbon (%C) and nitrogen (%N) content and stable isotope ratios of C (δ13C) and N (δ15N) were analyzed on a FlashEA 1112 elemental analyzer coupled to a Thermo Finnigan DELTA Plus Advantage mass spectrometer (Thermo Electron Corporation, Waltham, MA, USA) at the University of Jyväskylä, Finland. Stable isotope data are expressed in the conventional delta notation as the relative difference between ratios of samples and international standards (PeeDee Belemnite for δ13C, atmospheric N for δ15N). Analytical precision (mean SD from in-house standard) for each run was always better than 0.3‰ for δ13C and δ15N.
To calculate relative contribution of benthic and pelagic food sources to the diet of roach in both lakes, the Bayesian isotope mixing model package SIAR (Stable Isotope analysis in R)  was used in R version 2.12.0 . Fractionation factors were set to 0.4 ± 1.3‰ for δ13C and 3.4 ± 1.0‰ for δ15N, as suggested in a previous study . The mean isotopic value of the zooplankton samples of each lake was taken as an end point for pelagic δ13C and δ15N values. Littoral δ13C and δ15N end points were calculated from the mean value of the macroinvertebrates, containing trichoptera, isopoda, chironomidae and gastropoda. These taxa represent the most abundant littoral macroinvertebrates consumed by roach.
Association of shape data with genetic data and ecological variables
To reduce the dimensionality of the genotype information from 331 polymorphic AFLP loci, we performed Principal Coordinate analysis (PCoA) by applying the standardized covariance method. Genetic distance between two individuals was calculated as the sum of loci with different character states, a method that automatically discards monomorphic loci. These analyses were conducted in GenAlEx 6.5 [79, 80]. To examine association of shape with genotype, general linear models (GLM) between the dominant canonical variate (CV1) of shape data as dependent variable and lake as group variable were calculated. We used contrasts between Döllnsee (the lake with the most deviating shape) and the other three lakes. The first two axes of PCoA on AFLP-data and their potential interactions with lake origin were then used as continuous predictors. Using these GLMs, we tested whether the shape of roach differed among the four lakes, and whether the aggregated lake-specific genotype contributed to this differentiation. In a second step, we used only the three outlier loci and their interactions with lake origin as binary predictors of shape. With these GLMs, we tested whether the genotype outliers were directly associated with the differences in shape between the lakes. GLM analyses were conducted in R version 2.15.1. .
A two-block partial least square analysis (PLS), as implemented in MorphoJ, was used to study the association of shape and ecological variables at the level of individual roach. In all calculations, only those individuals for which all the respective variables were available to facilitate paired comparisons were included. All variables used in the second blocks were adjusted to standard deviate ; where si is the standard deviation of row or column i. To study correlation with predation pressure, shape was used as a first block and relative predation pressure as second block. To correlate shape with diet composition, PLS were calculated using percentage of ingested diet items obtained from GC or isotope ratios. A permutation test against the null hypothesis of independence (1000 randomization runs) was used for all PLS, as implemented in the software MorphoJ.
Two anonymous reviewers provided insightful comments on an earlier draft of this manuscript. We thank A Türck, C Helms, M Emmrich and S Oksanen for their help in the field and A Vogt, K Preuss, M Gamboa, K Pohlmann and S Trotzer for their help in the laboratory. D Thompson greatly helped with photographic procedures and figure generation. We thank in particular CP Klingenberg for organizing the Geometric Morphometrics Workshop 2011 in Telč, Czech Republic and helping with the first analyses. R. Michels, H. Mauersberger (Biosphärenreservat Schorfheide-Chorin), R Mauersberger (Förderverein Feldberg-Uckermärkische Seen e.V.) and R Tischbier (Stiftung Pro Artenvielfalt) kindly provided background information and access to the lakes. This study was financed by the TERRALAC-project (http://terralac.igb-berlin.de/) of the Wissenschaftsgemeinschaft Leibniz (WGL, Pakt für Forschung und Innovation). JS was supported by the IGB Fellowship program in Freshwater Science. KW was funded by the European Commission (FP-7 Marie Curie IIF, Co 237026).
- Koehl MAR: When does morphology matter?. Annual Review of Ecology and Systematics. 1996, 27: 501-542. 10.1146/annurev.ecolsys.27.1.501.View ArticleGoogle Scholar
- Robinson BW, Wilson DS: Character release and displacement in fishes - a neglected literature. Am Nat. 1994, 144 (4): 596-627. 10.1086/285696.View ArticleGoogle Scholar
- Skulason S, Smith TB: Resource polymorphism in vertebrates. Trends Ecol Evol. 1995, 10 (9): 366-370. 10.1016/S0169-5347(00)89135-1.PubMedView ArticleGoogle Scholar
- Anderson O: Optimal foraging by largemouth bass in structured environments. Ecology. 1984, 65 (3): 851-861. 10.2307/1938059.View ArticleGoogle Scholar
- Ehlinger TJ: Learning and individual variation in bluegill foraging - habitat-specific techniques. Anim Behav. 1989, 38: 643-658. 10.1016/S0003-3472(89)80010-7.View ArticleGoogle Scholar
- Schluter D, McPhail JD: Ecological character displacement and speciation in sticklebacks. Am Nat. 1992, 140 (1): 85-108. 10.1086/285404.PubMedView ArticleGoogle Scholar
- Cresko WA, Baker JA: Two morphotypes of lacustrine threespine stickleback, Gasterosteus aculeatus, in Benka Lake, Alaska. Environ Biol Fishes. 1996, 45 (4): 343-350. 10.1007/BF00002526.View ArticleGoogle Scholar
- McPhail JD: Ecology and evolution of sympatric sticklebacks (Gasterosteus) - morphological and genetic evidence for a species pair in Enos Lake, British Columbia. Canadian Journal of Zoology-Revue Canadienne de Zoologie. 1984, 62 (7): 1402-1408. 10.1139/z84-201.View ArticleGoogle Scholar
- Olsson J, Eklov P: Habitat structure, feeding mode and morphological reversibility: factors influencing phenotypic plasticity in perch. Evolutionary Ecology Research. 2005, 7 (8): 1109-1123.Google Scholar
- Svanback R, Eklov P: Effects of habitat and food resources on morphology and ontogenetic growth trajectories in perch. Oecologia. 2002, 131 (1): 61-70. 10.1007/s00442-001-0861-9.View ArticleGoogle Scholar
- Svanback R, Eklov P: Morphology dependent foraging efficiency in perch: a trade-off for ecological specialization?. Oikos. 2003, 102 (2): 273-284. 10.1034/j.1600-0706.2003.12657.x.View ArticleGoogle Scholar
- Skulason S, Noakes DLG, Snorrason SS: Ontogeny of trophic morphology in 4 sympatric morphs of Arctic charr Salvenius alpinus in Thingvallavatn, Iceland. Biol J Linn Soc. 1989, 38 (3): 281-301. 10.1111/j.1095-8312.1989.tb01579.x.View ArticleGoogle Scholar
- Malmquist HJ, Snorrason SS, Skulason S, Jonsson B, Sandlund OT, Jonasson PM: Diet differentiation in polymorphic Arctic charr in Thingvallavatn, Iceland. J Anim Ecol. 1992, 61 (1): 21-35. 10.2307/5505.View ArticleGoogle Scholar
- Alekseyev SS, Samusenok VP, Matveev AN, Pichugin MY: Diversification, sympatric speciation, and trophic polymorphism of Arctic charr, Salvelinus alpinus complex, in Transbaikalia. Environ Biol Fishes. 2002, 64 (1–3): 97-114.Google Scholar
- Harvell CD: The ecology and evolution of inducible defenses. Q Rev Biol. 1990, 65 (3): 323-340. 10.1086/416841.PubMedView ArticleGoogle Scholar
- Dzialowski AR, Lennon JT, O'Brien WJ, Smith VH: Predator-induced phenotypic plasticity in the exotic cladoceran Daphnia lumholtzi. Freshwater Biology. 2003, 48 (9): 1593-1602. 10.1046/j.1365-2427.2003.01111.x.View ArticleGoogle Scholar
- Grant JWG, Bayly IAE: Predator induction of crests in morphs of the Daphnia carinata King complex. Limnol Oceanogr. 1981, 26 (2): 201-218. 10.4319/lo.1981.26.2.0201.View ArticleGoogle Scholar
- Cote IM: Effects of predatory crab effluent on byssus production in mussels. J Exp Mar Biol Ecol. 1995, 188 (2): 233-241. 10.1016/0022-0981(94)00197-L.View ArticleGoogle Scholar
- Bronmark C, Miner JG: Predator-induced phenotypical change in body morphology in Crucian carp. Science. 1992, 258 (5086): 1348-1350. 10.1126/science.258.5086.1348.PubMedView ArticleGoogle Scholar
- Januszkiewicz AJ, Robinson BW: Divergent walleye (Sander vitreus)-mediated inducible defenses in the centrarchid pumpkinseed sunfish (Lepomis gibbosus). Biol J Linn Soc. 2007, 90 (1): 25-36. 10.1111/j.1095-8312.2007.00708.x.View ArticleGoogle Scholar
- Persson L, Andersson J, Wahlstrom E, Eklov P: Size-specific interactions in lake systems: predator gape limitation and prey growth rate and mortality. Ecology. 1996, 77 (3): 900-911. 10.2307/2265510.View ArticleGoogle Scholar
- Day T, Pritchard J, Schluter D: A comparison of two sticklebacks. Evolution. 1994, 48 (5): 1723-1734. 10.2307/2410260.View ArticleGoogle Scholar
- Agrawal AA: Phenotypic plasticity in the interactions and evolution of species. Science. 2001, 294 (5541): 321-326. 10.1126/science.1060701.PubMedView ArticleGoogle Scholar
- Losos JB: Ecological character displacement and the study of adaption. Proc Natl Acad Sci U S A. 2000, 97 (11): 5693-5695. 10.1073/pnas.97.11.5693.PubMed CentralPubMedView ArticleGoogle Scholar
- Smith TB, Skulason S: Evolutionary significance of resource polymorphisms in fishes, amphibians, and birds. Annual Review of Ecology and Systematics. 1996, 27: 111-133. 10.1146/annurev.ecolsys.27.1.111.View ArticleGoogle Scholar
- Rundell RJ, Price TD: Adaptive radiation, nonadaptive radiation, ecological speciation and nonecological speciation. Trends Ecol Evol. 2009, 24 (7): 394-399. 10.1016/j.tree.2009.02.007.PubMedView ArticleGoogle Scholar
- Lundsgaard-Hansen B, Matthews B, Vonlanthen P, Taverna A, Seehausen O: Adaptive plasticity and genetic divergence in feeding efficiency during parallel adaptive radiation of whitefish (Coregonus spp.). J Evol Biol. 2013, 26: 483-498. 10.1111/jeb.12063.PubMedView ArticleGoogle Scholar
- Nosil P, Funk DJ, Ortiz-Barrientos D: Divergent selection and heterogeneous genomic divergence. Mol Ecol. 2009, 18 (3): 375-402. 10.1111/j.1365-294X.2008.03946.x.PubMedView ArticleGoogle Scholar
- Rice AM, Rudh A, Ellegren H, Qvarnstrom A: A guide to the genomics of ecological speciation in natural animal populations. Ecol Lett. 2011, 14 (1): 9-18. 10.1111/j.1461-0248.2010.01546.x.PubMedView ArticleGoogle Scholar
- Collin H, Fumagalli L: Evidence for morphological and adaptive genetic divergence between lake and stream habitats in European minnows (Phoxinus phoxinus, Cyprinidae). Mol Ecol. 2011, 20 (21): 4490-4502. 10.1111/j.1365-294X.2011.05284.x.PubMedView ArticleGoogle Scholar
- Kottelat M, Freyhof J: Handbook of European freshwater fishes. 2007, Cornol, Switzerland: Publications KottelatGoogle Scholar
- Svanback R, Eklov P, Fransson R, Holmgren K: Intraspecific competition drives multiple species resource polymorphism in fish communities. Oikos. 2008, 117 (1): 114-124. 10.1111/j.2007.0030-1299.16267.x.View ArticleGoogle Scholar
- Gresens SE: Grazer diversity, competition and the response of the periphyton community. Oikos. 1995, 73 (3): 336-346. 10.2307/3545957.View ArticleGoogle Scholar
- Hargeby A, Blindow I, Andersson G: Long-term patterns of shifts between clear and turbid states in Lake Krankesjon and Lake Takern. Ecosystems. 2007, 10 (1): 28-35.View ArticleGoogle Scholar
- Pardue WJ, Webb DH: A comparison of aquatic macroinvertebrates occurring in association with Eurasian watermilfoil (Myriophyllum spicatum L.) with those found in the open littoral-zone. Journal of Freshwater Ecology. 1985, 3 (1): 69-79. 10.1080/02705060.1985.9665093.View ArticleGoogle Scholar
- Beckett DC, Aartila TP, Miller AC: Contrasts in density of benthic invertebrates between macrophyte beds and open littoral patches in Eau-Galle Lake, Wisconsin. Am Midl Nat. 1992, 127 (1): 77-90. 10.2307/2426324.View ArticleGoogle Scholar
- Scheffer M, Hosper SH, Meijer ML, Moss B, Jeppesen E: Alternative equilibria in shallow lakes. Trends Ecol Evol. 1993, 8 (8): 275-279. 10.1016/0169-5347(93)90254-M.PubMedView ArticleGoogle Scholar
- Domenici P, Blake RW: The kinematics and performance of fish fast-start swimming. J Exp Biol. 1997, 200 (8): 1165-1178.PubMedGoogle Scholar
- Webb PW: Body form, locomotion and foraging in aquatic vertebrates. Am Zool. 1984, 24 (1): 107-120.View ArticleGoogle Scholar
- Weihs D: The mechanisms of rapid starting of slender fish. Biorheology. 1973, 10: 343-350.PubMedGoogle Scholar
- Walker JA: Ecological morphology of lacustrine threespine stickleback Gasterosteus aculeatus L (Gasterosteidae) body shape. Biol J Linn Soc. 1997, 61 (1): 3-50.Google Scholar
- Weihs D: Design-features and meachanics of axial locomotion in fish. Am Zool. 1989, 29 (1): 151-160.View ArticleGoogle Scholar
- Langerhans RB, Layman CA, Shokrollahi AM, DeWitt TJ: Predator-driven phenotypic diversification in Gambusia affinis. Evolution. 2004, 58 (10): 2305-2318.PubMedView ArticleGoogle Scholar
- Domenici P, Turesson H, Brodersen J, Bronmark C: Predator-induced morphology enhances escape locomotion in crucian carp. Proc R Soc B-Biol Sci. 2008, 275 (1631): 195-201. 10.1098/rspb.2007.1088.View ArticleGoogle Scholar
- Eklov P, Jonsson P: Pike predators induce morphological changes in young perch and roach. J Fish Biol. 2007, 70 (1): 155-164. 10.1111/j.1095-8649.2006.01283.x.View ArticleGoogle Scholar
- Christensen B, Persson L: Species-specific antipredatory behaviors - effects on prey choice in different habitats. Behav Ecol Sociobiol. 1993, 32 (1): 1-9.View ArticleGoogle Scholar
- Milano D, Cussac VE, Macchi PJ, Ruzzante DE, Alonso MF, Vigliano PH, Denegri MA: Predator associated morphology in Galaxias platei in Patagonian lakes. J Fish Biol. 2002, 61 (1): 138-156. 10.1111/j.1095-8649.2002.tb01742.x.View ArticleGoogle Scholar
- Bourdeau PE, Johansson F: Predator-induced morphological defences as by-products of prey behaviour: a review and prospectus. Oikos. 2012, 121 (8): 1175-1190. 10.1111/j.1600-0706.2012.20235.x.View ArticleGoogle Scholar
- Robinson BW, Wilson DS, Margosian AS, Lotito PT: Ecological and morphological differentiation of pumpinseed sunfish in lakes without bluegill sunfish. Evolutionary Ecology. 1993, 7 (5): 451-464. 10.1007/BF01237641.View ArticleGoogle Scholar
- Gerry SP, Vogelzang M, Ascher JM, Ellerby DJ: Variation in the diet and feeding morphology of polyphenic Lepomis macrochirus. J Fish Biol. 2013, 82 (1): 338-346. 10.1111/j.1095-8649.2012.03471.x.PubMedView ArticleGoogle Scholar
- Aguirre WE: Microgeographical diversification of threespine stickleback: body shape-habitat correlations in a small, ecologically diverse Alaskan drainage. Biol J Linn Soc. 2009, 98 (1): 139-151. 10.1111/j.1095-8312.2009.01267.x.View ArticleGoogle Scholar
- Frommen JG, Herder F, Engqvist L, Mehlis M, Bakker TCM, Schwarzer J, Thunken T: Costly plastic morphological responses to predator specific odour cues in three-spined sticklebacks (Gasterosteus aculeatus). Evolutionary Ecology. 2011, 25 (3): 641-656. 10.1007/s10682-010-9454-6.View ArticleGoogle Scholar
- O'Steen S, Cullum AJ, Bennett AF: Rapid evolution of escape ability in Trinidadian guppies (Poecilia reticulata). Evolution. 2002, 56 (4): 776-784.PubMedView ArticleGoogle Scholar
- Ghalambor CK, Reznick DN, Walker JA: Constraints on adaptive evolution: the functional trade-off between reproduction and fast-start swimming performance in the Trinidadian guppy (Poecilia reticulata). Am Nat. 2004, 164 (1): 38-50. 10.1086/421412.PubMedView ArticleGoogle Scholar
- Bonin A, Bellemain E, Eidesen PB, Pompanon F, Brochmann C, Taberlet P: How to track and assess genotyping errors in population genetics studies. Mol Ecol. 2004, 13 (11): 3261-3273. 10.1111/j.1365-294X.2004.02346.x.PubMedView ArticleGoogle Scholar
- Pompanon F, Bonin A, Bellemain E, Taberlet P: Genotyping errors: causes, consequences and solutions. Nat Rev Genet. 2005, 6 (11): 847-859. 10.1038/nrg1707.PubMedView ArticleGoogle Scholar
- Perez-Figueroa A, Garcia-Pereira MJ, Saura M, Rolan-Alvarez E, Caballero A: Comparing three different methods to detect selective loci using dominant markers. J Evol Biol. 2010, 23 (10): 2267-2276. 10.1111/j.1420-9101.2010.02093.x.PubMedView ArticleGoogle Scholar
- Campbell D, Bernatchez L: Generic scan using AFLP markers as a means to assess the role of directional selection in the divergence of sympatric whitefish ecotypes. Mol Biol Evol. 2004, 21 (5): 945-956. 10.1093/molbev/msh101.PubMedView ArticleGoogle Scholar
- Miner BG, Sultan SE, Morgan SG, Padilla DK, Relyea RA: Ecological consequences of phenotypic plasticity. Trends Ecol Evol. 2005, 20 (12): 685-692. 10.1016/j.tree.2005.08.002.PubMedView ArticleGoogle Scholar
- Ruuhijarvi J, Rask M, Vesala S, Westermark A, Olin M, Keskitalo J, Lehtovaara A: Recovery of the fish community and changes in the lower trophic levels in a eutrophic lake after a winter kill of fish. Hydrobiologia. 2010, 646 (1): 145-158. 10.1007/s10750-010-0186-y.View ArticleGoogle Scholar
- Storfer A, Murphy MA, Spear SF, Holderegger R, Waits LP: Landscape genetics: where are we now?. Mol Ecol. 2010, 19 (17): 3496-3514. 10.1111/j.1365-294X.2010.04691.x.PubMedView ArticleGoogle Scholar
- Manel S, Joost S, Epperson BK, Holderegger R, Storfer A, Rosenberg MS, Scribner KT, Bonin A, Fortin MJ: Perspectives on the use of landscape genetics to detect genetic adaptive variation in the field. Mol Ecol. 2010, 19 (17): 3760-3772. 10.1111/j.1365-294X.2010.04717.x.PubMedView ArticleGoogle Scholar
- Wright S: Isolation by distance. Genetics. 1943, 28 (2): 114-138.PubMed CentralPubMedGoogle Scholar
- Bookstein FL: Morphometric tools for landmark data. 1991, New York: Cambridge University PressGoogle Scholar
- Lappalainen J, Tarkan AS, Harrod C: A meta-analysis of latitudinal variations in life-history traits of roach, Rutilus rutilus, over its geographical range: linear or non-linear relationships?. Freshwater Biology. 2008, 53 (8): 1491-1501. 10.1111/j.1365-2427.2008.01977.x.View ArticleGoogle Scholar
- Rohlf FJ, Slice D: Extensions of the procrustes method for the optimal superimposition of landmarks. Syst Zool. 1990, 39 (1): 40-59. 10.2307/2992207.View ArticleGoogle Scholar
- Vos P, Hogers R, Bleeker M, Reijans M, Vandelee T, Hornes M, Frijters A, Pot J, Peleman J, Kuiper M, et al: AFLP - A new technique for DNA-fingerprinting. Nucleic Acids Res. 1995, 23 (21): 4407-4414. 10.1093/nar/23.21.4407.PubMed CentralPubMedView ArticleGoogle Scholar
- Ehrich D: AFLPDAT: a collection of R functions for convenient handling of AFLP data. Molecular Ecology Notes. 2006, 6 (3): 603-604. 10.1111/j.1471-8286.2006.01380.x.View ArticleGoogle Scholar
- Vekemans X: AFLP-SURV version 1.0: Distributed by the author. Laboratoire de Génétique et Ecologie Végétale. 2002, Belgium: Université Libre de BruxellesGoogle Scholar
- Holsinger KE, Lewis PO: Hickory: a package for analysis of population genetic data v 1.1. 2007, Storrs: Distributed by the author: Department of Ecology and Evolutionary Biology, University of ConneticutGoogle Scholar
- Foll M, Gaggiotti O: A genome-scan method to identify selected loci appropriate for both dominant and codominant Markers: a Bayesian perspective. Genetics. 2008, 180 (2): 977-993. 10.1534/genetics.108.092221.PubMed CentralPubMedView ArticleGoogle Scholar
- Persson L: Effects of reduced interspecific competition on resource utilization in perch (Perca fluviatilis). Ecology. 1986, 67 (2): 355-364. 10.2307/1938578.View ArticleGoogle Scholar
- Pagel T: Determination of individual reproductive success in natural pike (Esox lucius) population: a DNA-based parentage assignment approach. 2009, Master Thesis: Humboldt-University of Berlin, GermanyGoogle Scholar
- Radke RJ: Strukturbildende Prozesse in Fischartengemeinschaften mesotropher Seen des nordostdeutschen Tieflandes. 1998, Doctoral Thesis: University of Konstanz, GermanyGoogle Scholar
- Windell JT: Food analysis and rate of digestion. Methods for assessment of fish production in fresh water. Edited by: Ricker WE. 1968, Oxford: Blackwell Scientific Publications, 197-203.Google Scholar
- Parnell AC, Inger R, Bearhop S, Jackson AL: Source partitioning using stable isotopes: coping with too much variation. PLoS One. 2010, 5 (3): e9672-10.1371/journal.pone.0009672. 9610.1371/journal.pone.0009672PubMed CentralPubMedView ArticleGoogle Scholar
- R Development Core Team: R: A language and environment for statistical computing. 2012, Vienna, Austria: R foundation for statistical Computing, URL http://www.R-project.orgGoogle Scholar
- Post DM: Using stable isotopes to estimate trophic position: models, methods, and assumptions. Ecology. 2002, 83 (3): 703-718. 10.1890/0012-9658(2002)083[0703:USITET]2.0.CO;2.View ArticleGoogle Scholar
- Peakall R, Smouse PE: GENALEX 6: genetic analysis in excel. Population genetic software for teaching and research. Molecular Ecology Notes. 2006, 6 (1): 288-295. 10.1111/j.1471-8286.2005.01155.x.View ArticleGoogle Scholar
- Peakall R, Smouse PE: GenAlEx 6.5: genetic analysis in excel. Population genetic software for teaching and research-an update. Bioinformatics. 2012, 28 (19): 2537-2539. 10.1093/bioinformatics/bts460.PubMed CentralPubMedView ArticleGoogle 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.