Patterns of morphological changes and hybridization between sympatric whitefish morphs (Coregonus spp.) in a Swiss lake: a role for eutrophication?

Whitefish, genus Coregonus, show exceptional levels of phenotypic diversity with sympatric morphs occurring in numerous postglacial lakes in the northern hemisphere. Here, we studied the effects of human-induced eutrophication on sympatric whitefish morphs in the Swiss lake, Lake Thun. In particular, we addressed the questions whether eutrophication (i) induced hybridization between two ecologically divergent summer-spawning morphs through a loss of environmental heterogeneity, and (ii) induced rapid adaptive morphological changes through changes in the food web structure. Genetic analysis based on 11 microsatellite loci of 282 spawners revealed that the pelagic and the benthic morph represent highly distinct gene pools occurring at different relative proportions on all seven known spawning sites. Gill raker counts, a highly heritable trait, showed nearly discrete distributions for the two morphs. Multilocus genotypes characteristic of the pelagic morph had more gill rakers than genotypes characteristic of benthic morph. Using Bayesian methods, we found indications of recent but limited introgressive hybridization. Comparisons with historical gill raker data yielded median evolutionary rates of 0.24 haldanes and median selection intensities of 0.27 for this trait in both morphs for 1948-2004 suggesting rapid evolution through directional selection at this trait. However, phenotypic plasticity as an alternative explanation for this phenotypic change cannot be discarded. We hypothesize that both the temporal shifts in mean gill raker counts and the recent hybridization reflect responses to changes in the trophic state of the lake induced by pollution in the 1960s, which created novel selection pressures with respect to feeding niches and spawning site preferences.

Identification of candidate genes and physiological pathways involved in gonad deformation in whitefish (Coregonus spp.) from Lake Thun, Switzerland

In 2000, fishermen reported the appearance of deformed reproductive organs in whitefish (Coregonus spp.) from Lake Thun, Switzerland. Despite intensive investigations, the causes of these abnormalities remain unknown. Using gene expression profiling, we sought to identify candidate genes and physiological processes possibly associated with the observed gonadal deformations, in order to gain insights into potential causes. Using in situ-synthesized oligonucleotide arrays, we compared the expression levels at 21,492 unique transcript probes in liver and head kidney tissue of male whitefish with deformed and normally developed gonads, respectively. The fish had been collected on spawning sites of two genetically distinct whitefish forms of Lake Thun. We contrasted the gene expression profiles of 56 individuals, i.e., 14 individuals of each phenotype and of each population. Gene-by-gene analysis revealed weak expression differences between normal and deformed fish, and only one gene, ictacalcin, was found to be up-regulated in head kidney tissue of deformed fish from both whitefish forms, However, this difference could not be confirmed with quantitative real-time qPCR. Enrichment analysis on the level of physiological processes revealed (i) the involvement of immune response genes in both tissues, particularly those linked to complement activation in the liver, (ii) proteolysis in the liver and (iii) GTPase activity and Ras protein signal transduction in the head kidney. In comparison with current literature, this gene expression pattern signals a chronic autoimmune disease in the testes. Based on the recent observations that gonad deformations are induced through feeding of zooplankton from Lake Thun we hypothesize that a xenobiotic accumulated in whitefish via the plankton triggering autoimmunity as the likely cause of gonad deformations. We propose several experimental strategies to verify or reject this hypothesis.

Divergence along a steep ecological gradient in lake whitefish (Coregonus sp.)

To understand mechanisms structuring diversity in young adaptive radiations, quantitative and unbiased information about genetic and phenotypic diversity is much needed. Here, we present the first in-depth investigation of whitefish diversity in a Swiss lake, with continuous spawning habitat sampling in both time and space. Our results show a clear cline like pattern in genetics and morphology of populations sampled along an ecological depth gradient in Lake Neuchâtel. Divergent natural selection appears to be involved in shaping this cline given that trait specific P(ST)-values are significantly higher than F(ST)-values when comparing populations caught at different depths. These differences also tend to increase with increasing differences in depth, indicating adaptive divergence along a depth gradient, which persists despite considerable gene flow between adjacent demes. It however remains unclear, whether the observed pattern is a result of currently stable selection-gene flow balance, incipient speciation, or reverse speciation due to anthropogenic habitat alteration causing two formerly divergent species to collapse into a single gene pool.

Unexplained gonad alterations in whitefish (Coregonus spp.) from Lake Thun, Switzerland: levels of persistent organic pollutants in different morphs

Since 2000, a surprisingly high number of macroscopical gonad alterations has been reported in whitefish (Coregonus spp.) from Lake Thun, Switzerland. This unique phenomenon is still unexplained and has received much public attention. As one possible trigger for these effects, the presence of persistent, bioaccumulative and toxic compounds acting as endocrine disruptors in the lake has been discussed. In this study, concentrations of selected persistent organic pollutants were examined in two morphs of whitefish from Lake Thun and their link to the observed abnormalities was investigated. Analyzed compound classes included polychlorinated biphenyls, polychlorinated dibenzo-p-dioxins and dibenzofurans, polychlorinated naphthalenes, polybrominated diphenyl ethers and hexabromocyclododecanes. The target substances were identified in all samples and concentrations of the analyzed compounds were highly correlated among each other. These correlations show that the analyzed substances have the same distribution pattern throughout the lake and that uptake, accumulation and elimination processes are similar. Significant differences in contaminant levels within the samples existed between the two analyzed morphs of whitefish, most likely due to different age, food patterns and growth rate. No difference in contaminant levels was observed between fish with abnormal gonads and fish with normal gonads, suggesting no causal link between the investigated lipophilic organohalogen compounds present in fish and the observed gonad abnormalities in whitefish from Lake Thun. A comparison to existing data shows that concentrations in Lake Thun whitefish are at the lower bound of contaminant levels in whitefish from Swiss lakes or from European waters.

Long-term estrogen exposure of whitefish Coregonus lavaretus induces intersex but not Lake Thun-typical gonad malformations

A high prevalence of gonad morphological variations has been observed in whitefish Coregonus lavaretus from Lake Thun (Switzerland). To clarify the role of endocrine disruption as a possible cause of the gonad alterations, whitefish were reared in a long-term laboratory experiment under exposure to 17 beta-estradiol (E2). Fish were fed from first-feeding until 3 yr of age at a daily rate of 0 (control), 0.5 or 50 microg E2 kg(-1) fish. E2 exposure resulted in a time- and concentration-dependent increase of prevalence and intensity of intersex gonads, i.e. gonads that macroscopically appeared as either testis or ovary but microscopically contained both male and female germ cells. Four types of intersex could be distinguished: Types 1 and 2 were composed of mainly male tissue, with Type 1 containing single oocytes and Type 2 displaying an ovary-like lamellar structure of the tissue. In Type 3, an increased percentage of the tissue was occupied by female germ cells, while in Type 4, the majority of the gonad tissue consisted of female germ cells. Chronic E2 exposure additionally resulted in a concentration-dependent shift of the sex ratio towards females, a reduced condition factor, retarded gonad growth together with delayed maturation of germ cells, and elevated levels of hepatic vitellogenin mRNA. However, Lake Thun-typical alterations of gonad morphology were not induced by chronic E2 exposure. The results provide evidence that estrogen-active compounds unlikely play a role in the etiology of gonad malformations in Lake Thun whitefish.

Adaptive plasticity and genetic divergence in feeding efficiency during parallel adaptive radiation of whitefish (Coregonus spp.)

Parallel phenotypic divergence in replicated adaptive radiations could either result from parallel genetic divergence in response to similar divergent selec- tion regimes or from equivalent phenotypically plastic response to the repeated occurrence of contrasting environments. In post-glacial fish, repli- cated divergence in phenotypes along the benthic-limnetic habitat axis is commonly observed. Here, we use two benthic-limnetic species pairs of whitefish from two Swiss lakes, raised in a common garden design, with reciprocal food treatments in one species pair, to experimentally measure whether feeding efficiency on benthic prey has a genetic basis or whether it underlies phenotypic plasticity (or both). To do so, we offered experimental fish mosquito larvae, partially burried in sand, and measured multiple feed- ing efficiency variables. Our results reveal both, genetic divergence as well as phenotypically plastic divergence in feeding efficiency, with the pheno- typically benthic species raised on benthic food being the most efficient forager on benthic prey. This indicates that both, divergent natural selection on genetically heritable traits and adaptive phenotypic plasticity, are likely important mechanisms driving phenotypic divergence in adaptive radiation.

Rapid parallel adaptive radiations from a single hybridogenic ancestral population

The Alpine lake whitefish (Coregonus lavaretus) species complex is a classic example of a recent radiation, associated with colonization of the Alpine lakes following the glacial retreat (less than 15 kyr BP). They have formed a unique array of endemic lake flocks, each with one to six described sympatric species differing in morphology, diet and reproductive ecology. Here, we present a genomic investigation of the relationships between and within the lake flocks. Comparing the signal between over 1000 AFLP loci and mitochondrial control region sequence data, we use phylogenetic tree-based and population genetic methods to reconstruct the phylogenetic history of the group and to delineate the principal centres of genetic diversity within the radiation. We find significant cytonuclear discordance showing that the genomically monophyletic Alpine whitefish clade arose from a hybrid swarm of at least two glacial refugial lineages. Within this radiation, we find seven extant genetic clusters centred on seven lake systems. Most interestingly, we find evidence of sympatric speciation within and parallel evolution of equivalent phenotypes among these lake systems. However, we also find the genetic signature of human-mediated gene flow and diversity loss within many lakes, highlighting the fragility of recent radiations.

Speciation leads to divergent methylmercury accumulation in sympatric whitefish

Central European lake whitefish (Coregonus spp.) colonized Swiss lakes following the last glacial retreat and have undergone rapid speciation and adaptive radiation. Up to six species have been shown to coexist in some lakes, and individual species occupy specific ecological niches and have distinct feeding and reproductive ecologies. We studied methylmercury (MeHg) accumulation in sympatric whitefish species from seven Swiss lakes to determine if ecological divergence has led to different rates of MeHg bioaccumulation. In four of seven lakes, sympatric species had distinctly different MeHg levels, which varied by up to a factor of two between species. Generally, species with greater MeHg levels were smaller in body size and planktivorous, and species with lower MeHg were larger and benthivorous. While modest disparities in trophic position between species might be expected a priori to explain the divergence in MeHg, δ15N of bulk tissue did not correlate with fish MeHg in five of seven lakes. Results of a nested ANCOVA analysis across all lakes indicated that only two factors (species, lake) explained substantial portions of the variance, with species accounting for more variance (52 %) than inter-lake differences (32 %). We suggest that differences in MeHg accumulation were likely caused by diverging metabolic traits between species, such as differences in energy partitioning between anabolism and catabolism, potentially interacting with species-specific prey resource utilization. These results indicate substantial variability in MeHg accumulation between closely related fish species, illustrating that ecological speciation in fish can lead to divergent MeHg accumulation patterns.

Intergenomic epistasis causes asynchronous hatch times in whitefish hybrids, but only when parental ecotypes differ

Support for the theory of ecological speciation requires evidence for ecological divergence between species which directly or indirectly causes reproductive isolation. This study investigates effects of ecological vs. genetic disparity of parental species on the presence of endogenous selection (deformation and mortality rates) and potential sources of exogenous selection (growth rates and hatch timing) on hybrids. Hybrid embryonic development is analysed in a common-garden full-sib cross of three species belonging to two different ecotypes within the Coregonus lavaretus species flock in the central Alpine region of Europe. Although hatch timing was similar across the three species, embryonic growth rates and egg sizes differed between ecotypes. This led to a mismatch between embryonic growth rate and egg size in hybrid crosses that reveals epistasis between the maternal and embryonic genomes and transgressive hatch times that were asynchronous with control crosses. A strong constraint of egg size to embryo size at late development was also evident. We argue that this demonstrates potential for coadaptation of a maternal trait (egg size) with offspring growth rate to be an important source of selection against hybridization between ecotypes with different egg sizes. Implications for the measurement and quantification of early life-history traits affected by this additive relationship, such as hatch day and larval size, are also discussed.

Genomic signatures of relaxed disruptive selection associated with speciation reversal in whitefish

Background: Speciation reversal: the erosion of species differentiation via an increase in introgressive hybridization due to the weakening of previously divergent selection regimes, is thought to be an important, yet poorly understood, driver of biodiversity loss. Our study system, the Alpine whitefish (Coregonus spp.) species complex is a classic example of a recent postglacial adaptive radiation: forming an array of endemic lake flocks, with the independent origination of similar ecotypes among flocks. However, many of the lakes of the Alpine radiation have been seriously impacted by anthropogenic nutrient enrichment, resulting in a collapse in neutral genetic and phenotypic differentiation within the most polluted lakes. Here we investigate the effects of eutrophication on the selective forces that have shaped this radiation, using population genomics. We studied eight sympatric species assemblages belonging to five independent parallel adaptive radiations, and one species pair in secondary contact. We used AFLP markers, and applied FST outlier (BAYESCAN, DFDIST) and logistic regression analyses (MATSAM), to identify candidate regions for disruptive selection in the genome and their associations with adaptive traits within each lake flock. The number of outlier and adaptive trait associated loci identified per lake were then regressed against two variables (historical phosphorus concentration and contemporary oxygen concentration) representing the strength of eutrophication. Results: Whilst we identify disruptive selection candidate regions in all lake flocks, we find similar trends, across analysis methods, towards fewer disruptive selection candidate regions and fewer adaptive trait/candidate loci associations in the more polluted lakes. Conclusions: Weakened disruptive selection and a concomitant breakdown in reproductive isolating mechanisms in more polluted lakes has lead to increased gene flow between coexisting Alpine whitefish species. We hypothesize that the resulting higher rates of interspecific recombination reduce either the number or extent of genomic islands of divergence surrounding loci evolving under disruptive natural selection. This produces the negative trend seen in the number of selection candidate loci recovered during genome scans of whitefish species flocks, with increasing levels of anthropogenic eutrophication: as the likelihood decreases that AFLP restriction sites will fall within regions of heightened genomic divergence and therefore be classified as FST outlier loci. This study explores for the first time the potential effects of human-mediated relaxation of disruptive selection on heterogeneous genomic divergence between coexisting species.