Rapid evolution of troglomorphic characters suggests selection rather than neutral mutation as a driver of eye reduction in cave crabs

Sebastian Klaus, José C. E. Mendoza, Jia Huan Liew, Martin Plath, Rudolf Meier, Darren C. J. Yeo

Abstract

This study asked whether reductive traits in cave organisms evolve at a slower pace (suggesting neutral evolution under relaxed selection) than constructive changes, which are likely to evolve under directional selection. We investigated 11 subterranean and seven surface populations of Sundathelphusa freshwater crabs on Bohol Island, Philippines, and examined constructive traits associated with improved food finding in darkness (increased leg and setae length) and reductive traits (reduced cornea size and eyestalk length). All changes occurred rapidly, given that the age of the most recent common ancestor was estimated to be 722–271 ka based on three mitochondrial markers. In order to quantify the speed of character change, we correlated the degree of morphological change with genetic distances between surface and subterranean individuals. The temporal pattern of character change following the transition to subterranean life was indistinguishable for constructive and reductive traits, characterized by an immediate onset and rapid evolutionary change. We propose that the evolution of these reductive traits—just like constructive traits—is most likely driven by strong directional selection.

1. Introduction

Cave animals are a promising model system for studying adaptation and ecological speciation in response to environmental gradients created by the presence or absence of light [14]. Many cave animals share two kinds of evolutionary modifications (‘troglomorphies’): characters that become functionless in darkness—most strikingly, the visual system and body pigmentation—are reduced (‘reductive traits’) [5], whereas ‘constructive traits’ are features that are related to improved non-visual orientation, navigation, food-finding and communication [6]. Typical constructive traits involve improved chemosensory capacities through elongated antennae, legs, fins, or setae [7,8].

It seems straightforward to assume that rapid evolutionary divergence in constructive characters is caused by strong, directional selection in the cave environment. Indeed, divergent selection may even promote ecological speciation [9]. Predictions regarding the speed of evolutionary reductive processes are less straightforward. In the classical literature on cave evolution, the view prevails that reductions are mainly due to the accumulation of neutral mutations under relaxed selection [1012]. If this was true, reductive traits should be fixed at a much slower pace than constructive evolutionary changes (see the electronic supplementary material, figure S1). Alternatively, reductive traits could also be under directional selection, if fitness advantages drive the reduction of expensive, dispensable traits [13,14] in the often resource-limited cave environment. For example, directional selection has been proposed for the amphipod Gammarus minus because eye reduction may allow improved processing of sensory input from the antennae in the brain [7]. Antagonistic pleiotropic effects have been suggested to affect eye reduction in the cave fish Astyanax mexicanus, where increasing the number of taste buds is developmentally correlated with eye reduction [6].

To date, no comprehensive study has compared the timing and degree of both constructive and reductive processes in closely related taxa that have repeatedly colonized caves at different points in time. A major obstacle (even for well-established model systems such as A. mexicanus) is the identification of the surface-living source population and, therefore, the timing of evolutionary change in different cave populations. Most studies can only provide broad time estimates for the most recent common ancestor of surface and subterranean forms [1517], and hence cannot correlate the degree of morphological change with the time spent in the cave environment.

Here, we investigate adaptations to the subterranean environment in gecarcinucid freshwater crabs on Bohol Island, Philippines. Crabs of the genus Sundathelphusa inhabit both surface streams and caves and show varying degrees of troglomorphism within a small area of less than 40 × 40 km (figure 1a), which led to the recognition of five species ([18]; figure 1c). Bohol Island emerged in the Plio- or Pleistocene (5.1–0.01 Mya) [19] indicating that its freshwater fauna is young. We sought to infer the onset and relative pace of constructive and reductive evolutionary change in cave forms of varying age and predicted an early onset in the evolution of constructive traits, whereas different evolutionary theories (see above) provide conflicting predictions for the timing of reductive processes.

Figure 1.

(a) Map of Bohol Island, Philippines with subterranean (black circles) and surface (white boxes) sampling sites. Geological formations are adopted from [19]. (b) Maximum clade credibility tree of the Bayesian analysis. Values on the branches are posterior probabilities, and bootstrap values of ML and MP analyses. Values at nodes represent 95% credibility intervals of divergence time estimates. Grey shading, subterranean taxa. (c) Different degrees of morphological adaptation to darkness. From top to bottom: Sundathelphusa boex; Sundathelphusa vedeniki; Sundathelphusa urichi; Sundathelphusa sottoae; Sundathelphusa cavernicola. Scale bars, 2 cm; frontal view not to scale.

2. Material and methods

Samples for the molecular genetic analyses (n = 38; 19 sites) are listed in the electronic supplementary material, table S1. Phylogenetic and age estimation was conducted with BEAST v. 1.6.1 [20] on a concatenated dataset of mitochondrial Cytb (408 bp), 12S rRNA (293 bp) and 16S rRNA (526 bp) sequences (EMBL accession no. HE999758HE999850), excluding site 19 (see the electronic supplementary material). For the 16S rRNA, we applied a mean (±s.d.) substitution rate of 2.03 ± 0.38% [21], and for the 12S rRNA and Cytb a mean of 2.0 per cent per million years with a broad s.d. of 2.0. For details on sample size, PCR conditions, phylogenetic analyses and divergence time estimation please refer to the electronic supplementary material.

For morphometric analyses of cave adaptation, we investigated n = 18 surface crabs from eight sites and n = 34 subterranean crabs from 12 caves (see the electronic supplementary material, table S2). We measured carapace height, width and length dextrally, and conducted a principal component analysis on these traits and extracted PC1 as a proxy of body size. Measurements of leg length, mean length of the setae on the dactylus of the second peraeopod, cornea area and eyestalk length were subjected to a general linear model with PC1 as covariate. The corresponding residuals were extracted as size corrected values.

To discriminate between slow, gradual trait reduction (i.e. relaxed selection; see the electronic supplementary material, figure S1), and rapid, immediate change of character states, we fitted a linear and a logarithmic curve onto trait residuals (y) of the subterranean populations as a function of the minimum uncorrected pairwise distance of the Cytb sequences to the closest surface relative (x) and compared both models with a likelihood ratio test. We tested for slope heterogeneity between traits using ANCOVA on log-transformed data (constructive characters with change of sign; see the electronic supplementary material). The results of this analysis could be biased if pre-adapted surface populations were overrepresented in the sample or if the closest surface relatives of the genetically most distant populations escaped sampling. We therefore tested the robustness of the analyses by analysing reduced datasets, excluding surface sites and/or those cave sites that were genetically most distant to the surface (greater than 4% pairwise distance).

3. Results

The age of the most recent common ancestor of the freshwater crabs of Bohol is compatible with the geological record and dates to the Pleistocene, with a 95% highest posterior density interval of 722–271 ka (figure 1b). The resulting substitution rate of Cytb was determined as 3.55–8.11% per Ma. Minimum uncorrected Cytb pairwise distances of cave crabs to their closest surface-dwelling relative ranged from 0 to 7.9 per cent (figure 2).

Figure 2.

Z-transformed residuals of morphological character states (corrected for body size) as a function of the minimum p-distances to the closest surface (black circles) or subterranean (white circles) relative (mean ± s.d. for each sampling site; horizontal line: mean value of surface sites).

Using Cytb pairwise distances as a proxy for the age of the respective cave populations revealed that the evolution of troglomorphic characters was fast for both reductive (cornea size, eyestalk length) and constructive (leg length, setae length) traits (figure 2). The onset of morphological changes appears to be immediate and there was no difference between constructive and reductive traits; indeed, character change was better explained by a logarithmic than by a linear function for all four traits (p < 0.001; ΔAIC 1.01–3.45; electronic supplementary material, figure S2).

The test for slope heterogeneity detected no significant difference in slopes between the four different traits (see the electronic supplementary material, table S3), suggesting that all troglomorphic traits evolved at a comparable speed. This result was not affected by the exclusion of surface sites or genetically distant cave populations.

4. Discussion

Using a comparative approach, we provide the first empirical test of whether eye reduction in cave organisms occurs at a slower pace (suggesting neutral evolution under relaxed selection) than ‘constructive’ changes that clearly evolve under directional selection [3]. The pace of character change after the shift from surface to cave-dwelling was virtually the same for constructive and reductive traits, with a rapid initial change. Also, the onset of reductive changes did not lag behind that of constructive changes. We propose that—for both kinds of modifications—this pattern is caused by strong, directional selection. This is in line with evo-devo studies on other cave organisms proposing that eye reduction is not simply caused by accumulation of selectively neutral mutations, but must be subject to selection [6,13,22].

Recent work on the underlying genetic and developmental mechanisms of eye reduction has been conducted in A. mexicanus [4], Asellus aquaticus [23] and G. minus [24]. These studies show that only few and homologous genes are involved in the eye reduction of vertebrates [4] and crustaceans [23,24]. More importantly, in A. mexicanus mutations only affect lens apoptosis, whereas all other genes involved in eye development seem to stay intact [4], which would be unlikely in the case of relaxed selection. Also, some genes involved in eye development (like sonic hedgehog) have an increased expression, and are not simply shut down [4]. In contrast to eye reduction, loss of pigmentation in Astyanax was proposed to be caused by neutral mutations, as the involved genes are both up- and downregulated, and thus not driven in one direction via selection [13]. This might also apply to cave crabs and should be addressed in the future based on fresh material.

Future studies should also investigate the proximate and ultimate factors responsible for the rapid evolutionary change reported here, considering both ecological factors and the developmental genetic basis of the observed morphological changes. Our study, however, is the first of its kind to show that eye reduction in a subterranean organism evolves at the same rapid pace as constructive changes, and this is most likely a direct effect of altered selective regimes in the subterranean environment.

Acknowledgements

We thank Peter K. L. Ng for his help and valuable discussion, and two anonymous reviewers for their constructive criticism. Financial support came from the National University of Singapore (grant no. R-154-000-519-720, R-154-000-465-133, R-154-000-476-112) and National Research Foundation and Economic Development Board (SPORE, COY-15-EWI-RCFSA/N197-1).

  • Received November 23, 2012.
  • Accepted January 4, 2013.

References

View Abstract