Evolution of an adaptive behavior and its sensory receptors promotes eye regression in blind cavefish
© Yoshizawa et al; licensee BioMed Central Ltd. 2012
Received: 3 October 2012
Accepted: 27 December 2012
Published: 27 December 2012
How and why animals lose eyesight during adaptation to the dark and food-limited cave environment has puzzled biologists since the time of Darwin. More recently, several different adaptive hypotheses have been proposed to explain eye degeneration based on studies in the teleost Astyanax mexicanus, which consists of blind cave-dwelling (cavefish) and sighted surface-dwelling (surface fish) forms. One of these hypotheses is that eye regression is the result of indirect selection for constructive characters that are negatively linked to eye development through the pleiotropic effects of Sonic Hedgehog (SHH) signaling. However, subsequent genetic analyses suggested that other mechanisms also contribute to eye regression in Astyanax cavefish. Here, we introduce a new approach to this problem by investigating the phenotypic and genetic relationships between a suite of non-visual constructive traits and eye regression.
Using quantitative genetic analysis of crosses between surface fish, the Pachón cavefish population and their hybrid progeny, we show that the adaptive vibration attraction behavior (VAB) and its sensory receptors, superficial neuromasts (SN) specifically found within the cavefish eye orbit (EO), are genetically correlated with reduced eye size. The quantitative trait loci (QTL) for these three traits form two clusters of congruent or overlapping QTL on Astyanax linkage groups (LG) 2 and 17, but not at the shh locus on LG 13. Ablation of EO SN in cavefish demonstrated a major role for these sensory receptors in VAB expression. Furthermore, experimental induction of eye regression in surface fish via shh overexpression showed that the absence of eyes was insufficient to promote the appearance of VAB or EO SN.
We conclude that natural selection for the enhancement of VAB and EO SN indirectly promotes eye regression in the Pachón cavefish population through an antagonistic relationship involving genetic linkage or pleiotropy among the genetic factors underlying these traits. This study demonstrates a trade-off between the evolution of a non-visual sensory system and eye regression during the adaptive evolution of Astyanax to the cave environment.
Keywordsanimal behavior regressive evolution constructive evolution neuromast hedgehog tradeoff quantitative trait locus eye QTL cluster adaptation
The dark and nutrient poor cave environment exerts substantial pressure upon cave-dwelling animals. Perhaps as a consequence of these limited resources, cave-adapted animals from most major phyla exhibit a remarkable convergence in morphological and physiological changes related to cave life, including features that are both constructive (the lengthening of legs, fins and antennae, the appearance of novel behaviors, and elaboration of non-visual sensory systems) and regressive (the reduction or loss of vision and pigmentation) [1, 2]. Although these changes have been studied in a diverse set of cave animals, the genetic and evolutionary mechanisms responsible for them remain poorly understood.
Although natural selection is probably involved in the evolution of constructive cave adapted phenotypes , the evolutionary forces driving regressive changes are less certain. Darwin suggested that eye degeneration - one of the most conspicuous traits found in cave animals - evolved due to disuse . This idea was refined by others to implicate neutral mutation and genetic drift as a consequence of relaxed selection for vision in the cave environment [5, 6]. The neutral hypothesis was favored until recently when the results of new genetic and developmental studies supported the adaptive evolution of eye regression in cave animals. Three competing hypotheses have been proposed for the adaptive evolution of eye regression: (1) direct natural selection against eyes to conserve energy in the resource poor cave environment ; (2) indirect selection against eyes to open sufficient space for the elaboration of constructive characters [2, 8–10]; and (3) indirect selection against eyes due to the enhancement of traits that are negatively linked to optic development by antagonistic pleiotropy [11, 12]. Distinguishing among these hypotheses has been difficult since the genetic basis of eye reduction is unknown in cave-dwelling animals.
The teleost Astyanax mexicanus is an excellent model organism for studying the evolution of traits associated with cave life, including eye regression [5, 13–16]. Within the past few million years, at least five independent colonizations by two different migrational waves of eyed surface fish have established 29 geographically isolated Astyanax cavefish populations in northeastern Mexico [17–20]. After subsequent radiation underground, the founder cavefish populations became isolated in separate caves and evolved eye regression, reduced pigmentation or albinism, enhanced sensory systems and behavioral changes associated with cave life [5, 11, 12, 21–32]. Despite this isolation, Astyanax surface fish and cavefish are interfertile in the laboratory, allowing the evolution of constructive and regressive traits to be studied by genetic analysis.
Previous studies reported that eye degeneration in Astyanax is triggered by lens apoptosis and dysfunction due to expanded sonic hedgehog (shh) gene expression along the embryonic midline [25, 26]. Additionally, shh hyper expression was shown to increase jaw width and taste bud number, and to mediate the expansion of the forebrain and hypothalamus [11, 22]. In a recent study, Elipot et al. found that shh modifies the hypothalamic serotonergic network and increases foraging efficiency in cavefish by shifting behavior from fighting to foraging . These experiments support the hypothesis that eye regression in Astyanax has evolved at least in part as a result of indirect selection against eyes in favor of increased feeding efficiency through pleiotropy of the shh genes. However, recent genetic studies have discovered 8 to 12 quantitative trait loci (QTL) involved in eye reduction in Astyanax cavefish [5, 7, 33], but none that are linked to either the shhA or shhB genes, suggesting that upstream modulators of the SHH signaling system and/or other genetic factor(s) may be important in eye regression (for review, see ref. ).
In the present study, we launch an alternative approach to address the evolutionary mechanisms involved in eye degeneration. Assuming a positive relationship between visual decay and the evolution of constructive changes in other sensory systems - as has been proposed by many previous investigators [2, 8, 10, 34, 35] - we investigated the genetic basis of a constructive trait, the vibration attraction behavior (VAB), to evaluate its possible relationship to eye degeneration. VAB is the swimming of cavefish toward an oscillating object, a behavior that has evolved repeatedly in different Astyanax cavefish populations [27, 36, 37] and which may be present in Amblyopsid cavefish as well [38, 39]. VAB is mediated by an increase in the number and size of cranial superficial neuromasts (SN) in cavefish . Although VAB is usually absent in surface fish, a small proportion of those raised in the laboratory can show a weak form of VAB. In wild populations of surface fish, VAB is presumably deleterious since it may be easily detectable by predators . In contrast, VAB is adaptive in cavefish since it increases foraging in an environment devoid of light with sparse food and no macroscopic predators [27, 40].
Here we report the results of genetic analyses of VAB, SN enhancement and eye size. Using this approach, we discovered that the QTL underlying VAB, a specific class of SN located within the cavefish eye orbit (EO), and reduced eyes form two distinct clusters of overlapping QTL that together explain a significant portion of the genetic variation underlying these traits. By ablation of EO SN, we discovered that these sensory receptors contribute to VAB. Moreover, by shh overexpression, we also showed that induction of eye degeneration in surface fish did not induce VAB or increase EO SN, suggesting that these traits are part of an antagonistic system impacting eye formation that is independent of SHH signaling, and that the extra space opened by eye regression is insufficient to promote the appearance of these constructive traits. Therefore, we propose that the adaptive evolution of VAB and EO SN enhancement has contributed to eye degeneration in Astyanax cavefish.
A small number of genetic factors control VAB and SN number
VAB and eye size are strongly correlated with EO SN but not SO-3 SN
Since SN are the sensory receptors that facilitate VAB , we expected the number of SN to be strongly correlated with the level of VAB in our hybrid families. In a previous study using cavefish, surface fish and their F1 progeny, we found a significant positive correlation between VAB and the number of SN at SO-3, which are located in the cranial region immediately ventral to the EO . However, in the present study using both F2 and F3 hybrids, we found that the correlation between these two traits was not significant (r = 0.13, P = 0.052; Figure 1G), although the correlation between VAB and SN diameter at SO-3 was significant (r = 0.14, P = 0.044; Figure 1H). This new result suggests that the number of SN at SO-3 may not be the most important parameter determining VAB. After further examination of SN within the cranial regions of surface fish and cavefish, we found that SN are located in the EO in cavefish (yellow dotted line in the bottom of Figure 1N) and in F2 and F3 hybrids with reduced eyes, but not in surface fish (yellow arrow in the bottom of Figure 1L, M: F3; Figure 1K: surface fish). In contrast to the results for SO-3 SN, we found that the number of EO SN was positively correlated with VAB (r = 0.29, P < 0.001; Figure 1F) and negatively correlated with eye size (r = -0.44, P < 0.001; Figure 1J), and that eye size was also negatively correlated with VAB (r = -0.26, P < 0.001; Figure 1I). These results suggest that VAB may be mediated by an increase in EO SN and that both traits are negatively correlated with eye size.
Genetic linkage mapping and QTL analysis
A summary of the location and effect size of significant QTLs for six traits
Effect size and direction
√NOA (0 to 7.68)
EO SN number
SN number (0 to 10)
Eye-size/SL (0.0558 to 0.1382)
Albino (0) - Pigmented (1)
Two QTL for cavefish VAB at 35 Hz but not other frequencies
Two QTL clusters control VAB, EO SN number and eye size
By performing stepwiseqtl, we detected nine QTL: two each for VAB and EO SN number on LG 2 and 17, and five for eye size on LGs 2, 4, 6, 16, and 17 (Figure 3, Table 1). The QTL for VAB and EO SN number were relatively consistent across all types of analyses, allowing for some variation since single-QTL models are not designed to detect QTL for complex traits. The eye size QTL on LG 2 and 17 were also consistent among analyses, but the eye QTL on LG 4, 6 and 16 were not evident in the single-QTL scans (Additional file 1B, left column). Although this result may be expected since eye size is genetically complex in Astyanax , these three QTL will need to be confirmed in future analyses. We did not detect significant epistasis among any of these QTL (Figure 3, Table 1). For each phenotype, the sum of detected QTL accounted for 19.8%, 18.5% and 40.1% of the total phenotypic variance in VAB, EO SN number and eye size, respectively (Table 1). Although the QTL for these three traits were found on five different linkage groups distributed throughout the Astyanax genome, the Bayesian credible intervals for six VAB, EO SN and eye size QTL formed congruent or overlapping clusters on LGs 2 and 17 (1 cM overlap at the 47 cM position in LG 2, and 5 cM overlap at the 3 to 8 cM position of LG 17; Figures 2 and 3B, C). Particularly noteworthy, the credible intervals for the two most significant eye size QTL overlap with the two QTL for VAB and EO SN number (Figure 3B, C). The observed number of QTL clusters for these three traits is significantly higher than expected by under a Poisson distribution (Χ2 = 98.2, df = 3, P = 3.8 × 10-21) . Since alleles located within short genomic regions are generally inherited together as a unit, these QTL clusters support the physical correlations we found for VAB, EO SN number and eye size (Figure 1F, I, J).
Further examination of these overlapping QTL supports the hypothesis that VAB is enhanced by expanding the number of SN within the EO at the cost of eyes. We found that the cavefish alleles at each QTL shift the VAB, SN and eye size phenotypes in the direction expected given the distribution of these phenotypes between the surface fish and cavefish P0; when VAB level was high there were more SN and smaller eyes than when VAB level was low or absent (Table 1; see insets of Figure 3B, C; and Additional file 3). None of the 28 candidate genes were associated with any of the 9 QTL. As reported previously [7, 33], the shhA locus was not located in or near any of the eye size QTL (see at LG13 of Figure 2). In addition, alphaA crsytallin (cryaa) and shroom2, which were located within eye QTL in a previous study , did not map to an eye size QTL in the present analysis. The significant clustering of VAB, EO SN and eye size QTL on LG 2 and 17 supports the conclusion that the genetic factors responsible for VAB and EO SN enhancement are also responsible for eye regression, either as a result of genetic hitchhiking or pleiotropy.
The EO SN are major receptors for VAB
shh-mediated eye degeneration is insufficient to facilitate VAB and EO SN
We have shown that eye regression in Astyanax cavefish is phenotypically and genetically correlated with the evolution of VAB, an adaptive trait controlling feeding efficiency , as well as the enhancement of its sensory receptors, SN specifically localized within the EO. Genetic mapping studies revealed clusters of congruent or overlapping QTL on LG 2 and 17 that control eye size, VAB and the number of EO SN. In addition, we have also demonstrated a significant role of the EO SN in VAB, and that the open space created by shh-induced eye degeneration is insufficient to promote VAB or SN enhancement. This experimental result further highlights the importance of the genetic linkage between eye loss and the VAB sensory system. Our results support the hypothesis that selection for an adaptive behavior and its underlying sensory receptors are enhanced in cavefish at the expense of eyes. To our knowledge, this is the first study to report phenotypic and genetic relationships between a regressive phenotype, an adaptive behavior and the sensory system responsible for this behavior.
We have discovered a significant genetic correlation between VAB and EO SN number but not between VAB and SO-3 SN number and size. This result was surprising in light of our previous studies, which found that removal of SO-3 SN reduced (but did not eliminate) VAB . A possible explanation for this discrepancy is that the physical SN ablation method used in our previous study did not discriminate effectively between EO SN and SO-3 SN, which are present in adjacent regions in the cavefish cranium, or that the actual contribution of SO-3 SN to VAB is minor compared to that of EO SN. Our current ablation study supports either explanation by showing that EO SN are the major VAB receptors, whereas SO-3 SN, despite their larger number, have only a minor role in VAB (Figures 1 and 4). Although our EO SN ablation experiment also ablated some SO-3 SN and canal neuromasts, we conclude that the latter receptors play a minor role in VAB for two reasons. First, in the SO-3 SN ablation experiment, more SN were ablated than EO SN. Yet, despite the loss of more SN in SO-3, the level of VAB was not attenuated significantly (Figure 4). In contrast, deleting a much smaller population of EO SN resulted in a dramatic and statistically significant reduction in VAB (Figure 4). Second, in several former studies, the involvement of canal neuromasts in VAB was excluded through three observations: (1) selective SN ablation was enough to decrease VAB even while maintaining canal neuromasts , (2) cavefish VAB was tuned to the frequency at the peak sensitivity of SN but not of canal neuromasts , and (3) the selective interference with canal neuromasts failed to attenuate the response . So what makes EO SN necessary for VAB? It is possible that there are structural or neural circuitry differences between the two populations of EO and SO-3 SN that may account for our results, or that the EO SN are better positioned above the eye socket to detect fluctuations in water due to the bulging shape of the adipose tissue plug that fills the eye orbits of cavefish  or bulging eye balls of F2 and F3 hybrids . However, since the level of VAB did not drop to zero in our EO SN ablation, it seems likely that other SN contribute to VAB as well, possibly including SN at SO-3 and on the dorsal trunk .
Our present investigation has also clarified the difference between the strong type of VAB observed in cavefish (10+ approaches towards a vibrating object) and the weaker form previously detected in a small proportion of surface fish raised in the laboratory (> 4 approaches) [27, 40]. Cavefish VAB peaks at 35 Hz, whereas surface fish VAB shows a continuous low activity from 5 to 35 Hz [27, 40]. In this study, VAB showed two genetic components detectable at 35 Hz, but none at 10 or 50 Hz, suggesting that the weak form of surface fish VAB appears without a major genomic change and could be accomplished through environmental effects and/or developmental plasticity. Considering that at least three cavefish populations (Pachón, Los Sabinos and Piedras cavefish) exhibit a strong form of VAB [27, 49], VAB may, therefore, represent a broadly distributed evolutionary trait in cavefish. The elaboration of SN at a new site, the EO, could have been a key factor in the evolution of VAB during adaptation to subterranean habitats [10, 40]. Determining the mechanisms underlying EO SN development and VAB among independently evolved cave populations [50, 51] could open a path to understanding the origins of neurological and behavioral novelty during evolution.
The major conclusions of this study are that EO SN number is strongly correlated with both VAB and eye size, and that all three traits are controlled by QTL on LGs 2 and 17 with congruent or overlapping Bayesian credible intervals. These results suggest that the gene loci responsible for VAB, EO SN number and eye-size may be subject to possible physical linkage and genetic hitchhiking or may in fact be the same genes with pleiotropic effects. The existence of QTL clusters could facilitate rapid adaptation to environmental changes by concurrently fixing multiple beneficial traits. Previous QTL analyses have also shown tight linkages and overlapping QTL for several other cavefish traits, including eye size, lens size, taste bud number, tooth number, melanophore number and chemical sense ability [7, 33]. This observation suggests that adaptation to cave life in Astyanax is the result of tradeoffs among many phenotypic traits, possibly through functional shifts within pluripotent neural crest cells which contribute to most of the above traits through developmental processes . However, one must note that widespread clustering of QTL does not necessarily entail a phenotypic tradeoff. In a previous scan of QTL in Astyanax, Protas et al. reported the overlap of six QTL for eye size and melanophore number (MelE), yet the phenotypic correlation between these traits is only r = -0.05 [7, 33].
The present investigation provides a clear example of a tradeoff between constructive traits (VAB and EO SN) and a regressive trait (eye size) in parallel to the SHH-signaling pathway. The short genomic distance between the QTL for eye regression, EO SN and VAB  could be a key factor in promoting the rapid adaptation of Astyanax to a novel environment. In other rapidly evolving animals, such as African cichlids and three-spined sticklebacks, QTL clusters have also been shown to control groups of potentially advantageous traits, including sex determination and body color, and body shape and plate number [53, 54].
We have discovered a significant genetic correlation among eye size, VAB and EO SN number. Our results, therefore, support the hypothesis that eye loss has evolved in Astyanax cavefish as a result of indirect selection against eyes due to a trade-off between two constructive traits, VAB and EO SN enhancement, based either on genetic hitchhiking or pleiotropy. There has been a long-standing debate on the neutrality or benefit of eye regression in cave animals [1, 2, 4, 34, 35]. Since cave-adapted arthropods exhibit reduced eyes and elongated sensory antennae [1, 2, 8], the loss of eyes as a result of selection for increased tactile sensitivity is an attractive general hypothesis for the convergent evolution of eye regression in cave fauna.
Biological materials and crosses
Astyanax mexicanus surface fish used in this study were laboratory raised descendants of original collections made in Balmorhea Springs State Park, Texas, and cavefish were laboratory raised descendants of original collections from Cueva de El Pachón (Pachón cavefish) in Tamaulipas, Mexico. We generated all hybrid progeny from an original mating between one pair of cavefish and surface fish [11, 41, 55, 56]. One pair of F1 hybrids from this cross was mated to generate 384 F2 individuals for phenotypic and genetic analysis, while another pair of F1 hybrids was crossed to generate a second group of about 200 F2 progeny. These F1 hybrids were randomly selected. We genotyped this second set of hybrids at markers 017a, 215d, 145a, 204d and 218e, which are near potential VAB (017a and 215d) and SN (145a, 204d and 218e) loci identified by a preliminary QTL scan using 384 F2 individuals. To overcome a reduced rate of spawning in this group, we chose three pairs of F2 progeny that were either homozygous or heterozygous for the cavefish alleles at the marker loci and crossed them to generate 91 F3 progeny. Although these hybrids exhibited a reduced rate of spawning, there was no noticeable difference in their survivability. Hybrids were fed living Artemia larvae and were maintained individually in 500 ml tanks or in groups of approximately 10 in 3 L tanks. We phenotyped each hybrid once for SN number and eye size and two or three times for VAB when they were one to two years old. Because of the limitations of tank size , hybrids of different ages approximate the same standard length (approximately 3 cm) and there is no positive relationship between age and SN number at either SO-3 or the EO (Additional file 4). Additionally, the results of our QTL scans did not differ when including age as a covariate (data not shown). The animal procedures used in this study were approved by the University of Maryland Animal Care and Use Committee and conform to NIH guidelines.
We isolated genomic DNA from fin-clips using the DNeasy Blood & Tissue Kit (Qiagen, Valencia, CA, USA) or the quick extraction protocol . We genotyped fish for 226 genome-wide polymorphic microsatellite loci described previously  and for polymorphisms within or near 20 candidate genes for SN number and behavioral variation. We successfully genotyped 93,208 (95.6%) of these microsatellite markers among 382 F2 and 11,375 (95.5%) of 125 selected microsatellite markers among 91 F3. Candidate gene homologues were identified in the fugu, medaka, stickleback, tetraodon and zebrafish genomes. Degenerate primers were designed to match the most conserved amino acid regions using CODEHOP  or SCPrimer . We performed degenerate PCR on genomic DNA extracted from the fin-clips of P0 surface fish and Pachón cavefish using Takara LA Taq DNA polymerase (Clontech, Mountain View, CA, USA) and the following touch-down program: 1 cycle at 94°C for 90 sec; 3 cycles each at 94°C for 20 sec, 63°C for 20 sec and 68°C for 2.5 minutes; 3 cycles each at 94°C for 20 sec, 61°C for 20 sec and 68°C for 2.5 minutes; 3 cycles each at 94°C for 20 sec, 59°C for 20 sec and 68°C for 2.5 minutes; 30 cycles at 94°C for 20 sec, 57°C for 20 sec and 68°C for 2.5 minutes; and finally 1 cycle at 72°C for 4 minutes, followed by a 4°C incubation. PCR fragments were characterized by either direct sequencing with the same degenerate primer used for PCR or by conventional sequencing after sub-cloning into the pCRII plasmid using the TOPO TA Cloning Kit (Invitrogen, Grand Island, NY, USA). The homology of each candidate gene fragment was confirmed by blastn. After sequencing each candidate gene region in surface fish and Pachón cavefish, the sequences were extended with primers designed from known sequences using the GenomeWalker™ Universal Kit (Clontech), and aligned using either Vector NTI (version 7, InforMax, Bethesda, MD, USA) or Sequencher™ (version 5, Gene Codes Corp., Ann Arbor, MI, USA) software to identify microsatellite and single nucleotide polymorphisms (SNPs). Microsatellite polymorphisms were genotyped as described by Protas et al.  after designing primers in Primer3Plus . SNPs were genotyped by either Taqman® (Roche Applied Science, Indianapolis, IN, USA) or HybProbe (Roche Applied Science) technology. Primers and fluorescent probes were designed using Primer Express (version 3.0, Applied Biosystems, Carlsbad, CA, USA) or LightCycler® Probe Design Software 2.0 (Roche Applied Science), respectively. The primer and probe sequences for these candidate genes are provided in Additional file 5.
Vibration attraction behavior
We assayed VAB as described previously [27, 49]. Briefly, four or five days before the beginning of an assay, individuals were acclimated in a cylindrical assay chamber (Pyrex 325 ml glass dish, 10 cm diameter × 5 cm high, Corning, Corning, NY, USA) filled with conditioned water (pH 6.8; conductivity approximately 600 μS). During the assays, vibration stimuli were created using a 7.5 mm-diameter glass rod vibrating at 35 Hz using a Leader LG1301 function generator (Leader Instruments Corp., Cypress, CA, USA) with an audio speaker (Pro Speakers, Apple, Cupertino, CA, USA). The number of approaches (NOA) to the vibrating rod was video recorded during a three-minute period under infrared illumination (880 nm wave length, BL41192-880 black light, Advanced Illumination, Rochester, VT, USA), and counted using ImageJ 1.42q software (NIH, Bethesda, MD, USA).
Eye size was determined from photographs of each fish by digitally measuring the diameter of each eyeball along its rostral-caudal axis using ImageJ. We standardized these measurements by dividing them by each fish's standard length (SL), the length of the body from the tip of the snout to the base of the caudal fin.
Neuromast vital staining
Neuromasts were vital stained as described previously [27, 49]. Briefly, fish were immersed in 25 μg/ml 2-(4-(dimethylamino)styryl)-N-ethylpyridinium iodide (DASPEI; Invitrogen, Eugene, OR, USA)  dissolved in conditioned water for one hour, followed by immersion in ice-cold 66.7 μg/ml Ethyl 3-aminobenzoate methanesulfonate salt (MS222, Sigma) in conditioned water. The specimens were viewed under a fluorescence microscope (Axioskop 2 equipped with 2.5× Plan-Neofluar lens with a numerical aperture of 0.075 and a FITC filter set; Zeiss, Göttingen, Germany); we then photographed the fish with a Zeiss Axiocam CCD camera. Neuromasts were quantified on images of DASPEI-stained fish using ImageJ software. SN were counted both in the epidermis over the cranial third suborbital (SO-3) bone  and within the orbit epidermis dorsal to the line of the suborbital canal neuromasts (indicated by red dashed line in Figure 1L-N). To determine SN size, the long diameter of the 10 largest SN in the same area was measured and averaged.
Linkage and QTL mapping
Following marker genotyping, we used MapManager QTXb20  software to construct an initial linkage map of the SNP and microsatellite markers using 384 F2 fish with complete genotype data and the Kosambi map function. We then used the program R/qtl  to identify statistically unlikely events, such as double crossovers in adjacent intervals or single crossovers in a small interval under maximum likelihood estimation , to refine marker order and to estimate inter-marker distances.
We performed QTL mapping using two strategies: (1) multiple interval mapping via the function stepwiseqtl in R/qtl ; and (2) single-QTL model mapping accounting for relatedness in the program QTLRel . All logarithms of the odds (LOD) significant thresholds were determined following 2,000 permutation tests after single-QTL model mapping with scanone implemented in R/qtl or scanOne in QTLRel. The 95% Bayesian credible intervals for QTL positions, which are analogous to 95% confidence intervals, were obtained using R/qtl. For multiple QTL mapping, we first calculated genotype probabilities for any missing genotype using the calc.genoprob function before performing the forward/backward model selection function implemented in stepwiseqtl. Following model selection, we chose the model with the maximum penalized LOD score generated from 1,200 permutation tests . For QTLRel mapping, we first generated pedigree charts for the F2 and F3 hybrid families and then calculated condensed identity coefficients and estimated variance components with the cic function. We used these variance components when scanning the genome for QTL using the scanOne function (see Additional file 2). The genome-wide significant thresholds (P < 0.05) for this latter analysis were calculated from 2,000 permutation tests.
SN were ablated by a modification of our previous method . Prior to ablation, we measured VAB and counted neuromasts stained with DASPEI. Only cavefish showing a high VAB level (more than 16 NOA) were selected for SN ablation. SN were ablated by applying Vetbond non-toxic tissue adhesive (3M, St. Paul, MN, USA) either to the EO or the SO-3 region using a micro nylon loop, a traditional embryological tool . After application of tissue adhesive to one side, cavefish were exposed to air for 10 seconds, and tissue adhesive was then applied to the same area on the opposite side, followed by a second cycle of air-drying. The treated fish were placed in a 10 cm-diameter cylindrical chamber containing conditioned water at room temperature. Within a day of active swimming the tissue adhesive usually peeled off the body, resulting in a void in the underlying field of SN. After four to six days of acclimation/recovery, SN ablated fish were subjected to the VAB assay followed immediately by staining for one hour with 25 μg/ml DASPEI.
Overexpression was carried out by injecting zebrafish shhA mRNA (20 to 80 pg) into two to four cell stage embryos as described previously . Capped shhA mRNA was synthesized from the pSP64T plasmid using the mMESSAGE mMACHINE® SP6 Kit (Life Technologies, Grand Island, NY, USA).
Correlation studies were conducted using IBM SPSS 20.0.0 software (IBM, Somers, NY, USA). Castle-Wright estimates were performed according to Lynch and Walsh . Significance of QTL clustering was performed using a goodness-of-fit . We divided the 1,500 cM Astyanax genome into 50 bins of 30 cM each and counted the number of bins with 0, 1, 2 and 3 QTL. We then compared this observed distribution to that expected from a Poisson distribution where the average number of QTL expected per bin is 0.18 (9 QTL divided by 50 bins of 30 cM each).
logarithms of the odds
melanophore number above the eye
number of approaches
quantitative trait locus
standard error of the mean
- shh :
standard length of fish
single nucleotide polymorphisms
suborbital bone 3
vibration attraction behavior.
We thank Z Sanford, R Fernandez, J Giarrusso, D Norman, S Ahmed, S Siddiqui, L Younger and J Park for fish maintenance; R Fernandez, D Norman, S Ahmed and J Park for genotyping and phenotyping; K Broman for helpful comments on QTL mapping; and A Parkhurst for technical assistance. A Japan Society for the Promotion of Science Postdoctoral Fellowship to MY and NSF (IBN-05384) and NIH (R01-EY014619) grants to WRJ supported this research.
- Culver DC: Cave Life, Evolution and Ecology. 1982, Cambridge: Harvard University Press, 189.View ArticleGoogle Scholar
- Culver DC, Pipan T: The Biology of Caves and other Subterranean Habitats. 2009, Oxford: Oxford University Press, 254.Google Scholar
- Protas M, Jeffery WR: Evolution and development in cave animals: from fish to crustaceans. Wiley Interdiscip Rev Dev Biol. 2012, 1: 823-845. 10.1002/wdev.61.PubMedPubMed CentralView ArticleGoogle Scholar
- Darwin C: On the Origin of Species Based on Natural Selection, or the Preservation of Favoured Races in the Struggle for Life. 1859, London: John Murray, 480-6Google Scholar
- Wilkens H: Evolution and genetics of epigean and cave Astyanax-fasciatus (Characidae, Pisces) - support for the neutral mutation theory. Evol Biol. 1988, 23: 271-367.Google Scholar
- Culver DC, Wilkens H: Critical review of the relevant theories of the evolution of subterranean animals. Ecosystems of the World. Edited by: Wilkens H, Culver DC, Humphreys WF. 2000, Amsterdam: Elsevier, 30: 381-398.Google Scholar
- Protas M, Conrad M, Gross JB, Tabin C, Borowsky R: Regressive evolution in the Mexican cave tetra, Astyanax mexicanus. Curr Biol. 2007, 17: 452-454. 10.1016/j.cub.2007.01.051.PubMedPubMed CentralView ArticleGoogle Scholar
- Jernigan RW, Culver DC, Fong DW: The dual role of selection and evolutionary history as reflected in genetic correlations. Evolution. 1994, 48: 587-596. 10.2307/2410471.View ArticleGoogle Scholar
- Schemmel C: Vergleichende Untersuchungen an den Hautsinnesorganen ober- und unterirdisch lebender Astyanax-Formen. Z Morph Tiere. 1967, 61: 255-316. 10.1007/BF00400988.View ArticleGoogle Scholar
- Franz-Odendaal TA, Hall BK: Modularity and sense organs in the blind cavefish, Astyanax mexicanus. Evol Dev. 2006, 8: 94-100. 10.1111/j.1525-142X.2006.05078.x.PubMedView ArticleGoogle Scholar
- Yamamoto Y, Byerly MS, Jackman WR, Jeffery WR: Pleiotropic functions of embryonic sonic hedgehog expression link jaw and taste bud amplification with eye loss during cavefish evolution. Dev Biol. 2009, 330: 200-211. 10.1016/j.ydbio.2009.03.003.PubMedPubMed CentralView ArticleGoogle Scholar
- Jeffery WR: Adaptive evolution of eye degeneration in the Mexican blind cavefish. J Hered. 2005, 96: 185-196. 10.1093/jhered/esi028.PubMedView ArticleGoogle Scholar
- Mitchell RW, Russell WH, Elliott WR: Mexican Eyeless Characin Fishes, Genus Astyanax: Environment, Distribution, and Evolution. 1977, Lubbock, Texas, USA: Texas Tech University Press, 89:Google Scholar
- Jeffery WR: Cavefish as a model system in evolutionary developmental biology. Dev Biol. 2001, 231: 1-12. 10.1006/dbio.2000.0121.PubMedView ArticleGoogle Scholar
- Jeffery WR: Evolution and development in the cavefish Astyanax. Curr Top Dev Biol. 2009, 86: 191-221.PubMedPubMed CentralView ArticleGoogle Scholar
- Jeffery WR: Emerging model systems in evo-devo: cavefish and microevolution of development. Evol Dev. 2008, 10: 265-272. 10.1111/j.1525-142X.2008.00235.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Strecker U, Hausdorf B, Wilkens H: Parallel speciation in Astyanax cave fish (Teleostei) in northern Mexico. Mol Phylogenet Evol. 2012, 62: 62-70. 10.1016/j.ympev.2011.09.005.PubMedView ArticleGoogle Scholar
- Ornelas-García CP, Domínguez-Domínguez O, Doadrio I: Evolutionary history of the fish genus Astyanax Baird & Girard (1854) (Actinopterygii, Characidae) in Mesoamerica reveals multiple morphological homoplasies. BMC Evol Biol. 2008, 8: 340-10.1186/1471-2148-8-340.PubMedPubMed CentralView ArticleGoogle Scholar
- Gross JB: The complex origin of Astyanax cavefish. BMC Evol Biol. 2012, 12: 105-10.1186/1471-2148-12-105.PubMedPubMed CentralView ArticleGoogle Scholar
- Bradic M, Beerli P, Garcia-de Leon FJ, Esquivel-Bobadilla S, Borowsky RL: Gene flow and population structure in the Mexican blind cavefish complex (Astyanax mexicanus). BMC Evol Biol. 2012, 12: 9-10.1186/1471-2148-12-9.PubMedPubMed CentralView ArticleGoogle Scholar
- Gross JB, Borowsky R, Tabin CJ: A novel role for Mc1r in the parallel evolution of depigmentation in independent populations of the cavefish Astyanax mexicanus. PLoS Genet. 2009, 5: e1000326-10.1371/journal.pgen.1000326.PubMedPubMed CentralView ArticleGoogle Scholar
- Menuet A, Alunni A, Joly JS, Jeffery WR, Rétaux S: Expanded expression of sonic hedgehog in Astyanax cavefish: multiple consequences on forebrain development and evolution. Development. 2007, 134: 845-855. 10.1242/dev.02780.PubMedView ArticleGoogle Scholar
- Protas ME, Hersey C, Kochanek D, Zhou Y, Wilkens H, Jeffery WR, Zon LI, Borowsky R, Tabin CJ: Genetic analysis of cavefish reveals molecular convergence in the evolution of albinism. Nat Genet. 2006, 38: 107-111. 10.1038/ng1700.PubMedView ArticleGoogle Scholar
- Varatharasan N, Croll RP, Franz-Odendaal T: Taste bud development and patterning in sighted and blind morphs of Astyanax mexicanus. Dev Dyn. 2009, 238: 3056-3064. 10.1002/dvdy.22144.PubMedView ArticleGoogle Scholar
- Yamamoto Y, Jeffery WR: Central role for the lens in cave fish eye degeneration. Science. 2000, 289: 631-633. 10.1126/science.289.5479.631.PubMedView ArticleGoogle Scholar
- Yamamoto Y, Stock DW, Jeffery WR: Hedgehog signalling controls eye degeneration in blind cavefish. Nature. 2004, 431: 844-847. 10.1038/nature02864.PubMedView ArticleGoogle Scholar
- Yoshizawa M, Gorički Š, Soares D, Jeffery WR: Evolution of a behavioral shift mediated by superficial neuromasts helps cavefish find food in darkness. Curr Biol. 2010, 20: 1631-1636. 10.1016/j.cub.2010.07.017.PubMedPubMed CentralView ArticleGoogle Scholar
- Duboué ER, Keene AC, Borowsky RL: Evolutionary convergence on sleep loss in cavefish populations. Curr Biol. 2011, 21: 671-676. 10.1016/j.cub.2011.03.020.PubMedView ArticleGoogle Scholar
- Coombs S, Patton P, Windsor S: Active wall following by Mexican blind cavefish (Astyanax mexicanus). J Comp Physiol A Neuroethol Sens Neural Behav Physiol. 2010, 196: 853-867. 10.1007/s00359-010-0567-8.PubMedView ArticleGoogle Scholar
- Sharma S, Coombs S, Patton P, De Perera TB: The function of wall-following behaviors in the Mexican blind cavefish and a sighted relative, the Mexican tetra (Astyanax). J Comp Physiol A Neuroethol Sens Neural Behav Physiol. 2009, 195: 225-240. 10.1007/s00359-008-0400-9.PubMedView ArticleGoogle Scholar
- Windsor SP, Norris SE, Cameron SM, Mallinson GD, Montgomery JC: The flow fields involved in hydrodynamic imaging by blind Mexican cave fish (Astyanax fasciatus). Part II: gliding parallel to a wall. J Exp Biol. 2010, 213: 3832-3842. 10.1242/jeb.040790.PubMedView ArticleGoogle Scholar
- Elipot Y, Hinaux H, Callebert J, Rétaux S: Evolutionary shift from fighting to foraging in blind cavefish through changes in the serotonin network. Curr Biol. 2012.Google Scholar
- Protas M, Tabansky I, Conrad M, Gross JB, Vidal O, Tabin CJ, Borowsky R: Multi-trait evolution in a cave fish, Astyanax mexicanus. Evol Dev. 2008, 10: 196-209. 10.1111/j.1525-142X.2008.00227.x.PubMedView ArticleGoogle Scholar
- Porter ML, Crandall KA: Lost along the way: the significance of evolution in reverse. Trends Ecol Evol. 2003, 18: 541-547. 10.1016/S0169-5347(03)00244-1.View ArticleGoogle Scholar
- Barr TC: Cave ecology and the evolution of troglobites. Evol Biol. 1968, 2: 35-102.Google Scholar
- Parzefall J: Field observation in epigean and cave populations of Mexican characid Astyanax mexicanus (Pisces, Characidae). Mém Biospéléol. 1983, 10: 171-176.Google Scholar
- Abdel-Latif H, Hassan ES, von Campenhausen C: Sensory performance of blind Mexican cave fish after destruction of the canal neuromasts. Naturwissenschaften. 1990, 77: 237-239. 10.1007/BF01138492.PubMedView ArticleGoogle Scholar
- Hill LG: Feeding and food habits of the spring cavefish, Chologaster agassizi. Am Midl Nat. 1969, 82: 110-116. 10.2307/2423821.View ArticleGoogle Scholar
- Eigenmann CH: Cave Vertebrates of America: a Study in Degenerative Evolution. 1909, Washington, DC: The Carnegie Institution of Washington Publications no.104View ArticleGoogle Scholar
- Yoshizawa M, Jeffery WR: Evolutionary tuning of an adaptive behavior requires enhancement of the neuromast sensory system. Commun Integr Biol. 2011, 4: 89-91.PubMedPubMed CentralView ArticleGoogle Scholar
- Yamamoto Y, Espinasa L, Stock DW, Jeffery WR: Development and evolution of craniofacial patterning is mediated by eye-dependent and -independent processes in the cavefish Astyanax. Evol Dev. 2003, 5: 435-446. 10.1046/j.1525-142X.2003.03050.x.PubMedView ArticleGoogle Scholar
- Gross JB, Protas M, Conrad M, Scheid PE, Vidal O, Jeffery WR, Borowsky R, Tabin CJ: Synteny and candidate gene prediction using an anchored linkage map of Astyanax mexicanus. Proc Natl Acad Sci USA. 2008, 105: 20106-20111. 10.1073/pnas.0806238105.PubMedPubMed CentralView ArticleGoogle Scholar
- Kirby RF, Thompson KW, Hubbs C: Karyotypic similarities between the Mexican and blind tetras. Copeia. 1977, 1977: 578-580. 10.2307/1443283.View ArticleGoogle Scholar
- Cheng R, Abney M, Palmer AA, Skol AD: QTLRel: an R Package for genome-wide association studies in which relatedness is a concern. BMC Genet. 2011, 12: 66.PubMedPubMed CentralView ArticleGoogle Scholar
- Broman KW, Wu H, Sen S, Churchill GA: R/qtl: QTL mapping in experimental crosses. Bioinformatics. 2003, 19: 889-890. 10.1093/bioinformatics/btg112.PubMedView ArticleGoogle Scholar
- Manichaikul A, Moon JY, Sen S, Yandell BS, Broman KW: A model selection approach for the identification of quantitative trait loci in experimental crosses, allowing epistasis. Genetics. 2009, 181: 1077-1086. 10.1534/genetics.108.094565.PubMedPubMed CentralView ArticleGoogle Scholar
- Albertson RC, Streelman JT, Kocher TD: Directional selection has shaped the oral jaws of Lake Malawi cichlid fishes. Proc Natl Acad Sci USA. 2003, 100: 5252-5257. 10.1073/pnas.0930235100.PubMedPubMed CentralView ArticleGoogle Scholar
- Rapo MA, Jiang HS, Grosenbaugh MA, Coombs S: Using computational fluid dynamics to calculate the stimulus to the lateral line of a fish in still water. J Exp Biol. 2009, 212: 1494-1505. 10.1242/jeb.026732.PubMedView ArticleGoogle Scholar
- Yoshizawa M, Ashida G, Jeffery WR: Parental genetic effects in a cavefish adaptive behavior explain disparity between nuclear and mitochondrial DNA. Evolution. 2012, 66: 2975-2982. 10.1111/j.1558-5646.2012.01651.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Wilkens H, Strecker U: Convergent evolution of the cavefish Astyanax (Characidae, Teleostei): genetic evidence from reduced eye-size and pigmentation. Biol J Linn Soc Lond. 2003, 80: 545-554. 10.1111/j.1095-8312.2003.00230.x.View ArticleGoogle Scholar
- Borowsky R: Restoring sight in blind cavefish. Curr Biol. 2008, 18: R23-R24. 10.1016/j.cub.2007.11.023.PubMedView ArticleGoogle Scholar
- Hall BK: Evolutionary Developmental Biology. 1999, New York: Springer, 2View ArticleGoogle Scholar
- Jones FC, Grabherr MG, Chan YF, Russell P, Mauceli E, Johnson J, Swofford R, Pirun M, Zody MC, White S, Birney E, Searle S, Schmutz J, Grimwood J, Dickson MC, Myers RM, Miller CT, Summers BR, Knecht AK, Brady SD, Zhang H, Pollen AA, Howes T, Amemiya C, Baldwin J, Bloom T, Jaffe DB, Nicol R, Wilkinson J, Lander ES, et al: The genomic basis of adaptive evolution in threespine sticklebacks. Nature. 2012, 484: 55-61. 10.1038/nature10944.PubMedPubMed CentralView ArticleGoogle Scholar
- Roberts RB, Ser JR, Kocher TD: Sexual conflict resolved by invasion of a novel sex determiner in Lake Malawi cichlid fishes. Science. 2009, 326: 998-1001. 10.1126/science.1174705.PubMedPubMed CentralView ArticleGoogle Scholar
- Yoshizawa M, Jeffery WR: Shadow response in the blind cavefish Astyanax reveals conservation of a functional pineal eye. J Exp Biol. 2008, 211: 292-299. 10.1242/jeb.012864.PubMedPubMed CentralView ArticleGoogle Scholar
- Jeffery WR, Yamamoto Y: The lens is a regulator of craniofacial development and evolution in the teleost Astyanax. Dev Biol. 2000, 222: 239.Google Scholar
- Gallo ND, Jeffery WR: Evolution of space dependent growth in the teleost Astyanax mexicanus. PLoS ONE. 2012, 7: e41443-10.1371/journal.pone.0041443.PubMedPubMed CentralView ArticleGoogle Scholar
- Nusslein-Volhard C, Dahm R: Zebrafish: A Practical Approach. 2002, New York: Oxford University Press, 303.Google Scholar
- Rose TM, Henikoff JG, Henikoff S: CODEHOP (COnsensus-DEgenerate Hybrid Oligonucleotide Primer) PCR primer design. Nucleic Acids Res. 2003, 31: 3763-3766. 10.1093/nar/gkg524.PubMedPubMed CentralView ArticleGoogle Scholar
- Jabado OJ, Palacios G, Kapoor V, Hui J, Renwick N, Zhai J, Briese T, Lipkin WI: Greene SCPrimer: a rapid comprehensive tool for designing degenerate primers from multiple sequence alignments. Nucleic Acids Res. 2006, 34: 6605-6611. 10.1093/nar/gkl966.PubMedPubMed CentralView ArticleGoogle Scholar
- Untergasser A, Nijveen H, Rao X, Bisseling T, Geurts R, Leunissen JAM: Primer3Plus, an enhanced web interface to Primer3. Nucleic Acids Res. 2007, 35: W71-W74. 10.1093/nar/gkm306.PubMedPubMed CentralView ArticleGoogle Scholar
- Jørgensen JM: Evolution of octavolateralis sensory cells. The Mechanosensory Lateral Line. Edited by: Coombs S, Görner P, Münz H. 1989, New York: Springer-Verlag, 115-145.View ArticleGoogle Scholar
- Manly KF, Cudmore RH, Meer JM: Map Manager QTX, cross-platform software for genetic mapping. Mamm Genome. 2001, 12: 930-932. 10.1007/s00335-001-1016-3.PubMedView ArticleGoogle Scholar
- Lincoln SE, Lander ES: Systematic detection of errors in genetic linkage data. Genomics. 1992, 14: 604-610. 10.1016/S0888-7543(05)80158-2.PubMedView ArticleGoogle Scholar
- Broman KW, Sen S: A Guide to QTL Mapping with R/qtl. 2009, New York: Springer-Verlag, 396-1View ArticleGoogle Scholar
- Lynch M, Walsh B: Genetics and Analysis of Quantitative Traits. 1998, Sunderland, Massachusetts, USA: Sinauer Associates, Inc.Google 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.