The evolutionary forces driving the reduction of eyes and pigmentation in cave-adapted animals are unknown; Darwin famously questioned the role of natural selection in eye loss in cave fishes: "As it is difficult to imagine that eyes, although useless, could be in any way injurious to animals living in darkness, I attribute their loss wholly to disuse"[1]. We studied the genetics of eye and pigmentation regression in the Mexican cave tetra, Astyanax mexicanus, by mapping and quantitative trait loci (QTL) analysis. We also mapped QTL for the putatively constructive traits of jaw size, tooth number, and numbers of taste buds. The data suggest that eyes and pigmentation regressed through different mechanisms. Cave alleles at every eye or lens QTL we detected caused size reductions, consistent with evolution by natural selection but not with drift. QTL polarities for melanophore number were mixed, however, consistent with genetic drift. Arguments against a role for selection in the regression of cave-fish eyes cited the insignificant cost of their development [2, 3], but we argue that the energetic cost of their maintenance is sufficiently high for eyes to be detrimental in the cave environment. Regression can be caused either by selection or drift.

Understanding the molecular basis of repeatedly evolved phenotypes can yield key insights into the evolutionary process. Quantifying gene flow between populations is especially important in interpreting mechanisms of repeated phenotypic evolution, and genomic analyses have revealed that admixture occurs more frequently between diverging lineages than previously thought. In this study, we resequenced 47 whole genomes of the Mexican tetra from three cave populations, two surface populations, and outgroup samples. We confirmed that cave populations are polyphyletic and two Astyanax mexicanus lineages are present in our dataset. The two lineages likely diverged much more recently than previous mitochondrial estimates of 5–7mya. Divergence of cave populations from their phylogenetically closest surface population likely occurred between ~161k - 191k generations ago. The favored demographic model for most population pairs accounts for divergence with secondary contact and heterogeneous gene flow across the genome, and we rigorously identified gene flow among all lineages sampled. Therefore, the evolution of cave-related traits occurred more rapidly than previously thought, and trogolomorphic traits are maintained despite gene flow with surface populations. The recency of these estimated divergence events suggests that selection may drive the evolution of cave-derived traits, as opposed to disuse and drift. Finally, we show that a key trogolomorphic phenotype QTL is enriched for genomic regions with low divergence between caves, suggesting that regions important for cave phenotypes may be transferred between caves via gene flow. Our study shows that gene flow must be considered in studies of independent, repeated trait evolution.

The evolution of degenerate characteristics remains a poorly understood phenomenon. Only recently has the identification of mutations underlying regressive phenotypes become accessible through the use of genetic analyses. Focusing on the Mexican cave tetra Astyanax mexicanus, we describe, here, an analysis of the brown mutation, which was first described in the literature nearly 40 years ago. This phenotype causes reduced melanin content, decreased melanophore number, and brownish eyes in convergent cave forms of A. mexicanus. Crosses demonstrate non-complementation of the brown phenotype in F2 individuals derived from two independent cave populations: Pachón and the linked Yerbaniz and Japonés caves, indicating the same locus is responsible for reduced pigmentation in these fish. While the brown mutant phenotype arose prior to the fixation of albinism in Pachón cave individuals, it is unclear whether the brown mutation arose before or after the fixation of albinism in the linked Yerbaniz/Japonés caves. Using a QTL approach combined with sequence and functional analyses, we have discovered that two distinct genetic alterations in the coding sequence of the gene Mc1r cause reduced pigmentation associated with the brown mutant phenotype in these caves. Our analysis identifies a novel role for Mc1r in the evolution of degenerative phenotypes in blind Mexican cavefish. Further, the brown phenotype has arisen independently in geographically separate caves, mediated through different mutations of the same gene. This example of parallelism indicates that certain genes are frequent targets of mutation in the repeated evolution of regressive phenotypes in cave-adapted species.

Natural populations subjected to strong environmental selection pressures offer a window into the genetic underpinnings of evolutionary change. Cavefish populations, Astyanax mexicanus (Teleostei: Characiphysi), exhibit repeated, independent evolution for a variety of traits including eye degeneration, pigment loss, increased size and number of taste buds and mechanosensory organs, and shifts in many behavioural traits. Surface and cave forms are interfertile making this system amenable to genetic interrogation; however, lack of a reference genome has hampered efforts to identify genes responsible for changes in cave forms of A. mexicanus. Here we present the first de novo genome assembly for Astyanax mexicanus cavefish, contrast repeat elements to other teleost genomes, identify candidate genes underlying quantitative trait loci (QTL), and assay these candidate genes for potential functional and expression differences. We expect the cavefish genome to advance understanding of the evolutionary process, as well as, analogous human disease including retinal dysfunction.

When surface species colonize caves, a characteristic suite of traits eventually evolves over time, regardless of species. The genetic basis of the inevitable appearance of these very similar phenotypes was investigated through quantitative trait loci (QTL) mapping of 12 traits that differ significantly between the recently evolved (<1 Myr). Mexican cave tetra and its surface conspecific. The traits were a representative set, including eye size, pigment cell numbers, chemical sensitivity, body and skull morphology, standard length, and metabolism. We used both single- and multi-trait models for QTL mapping. QTL effects of these traits were significantly clustered in the genome. We mapped 13 regions in the genome with QTL effects on from three to nine traits. These clusters could be multigenic or could represent single locus with pleiotropic alleles. Given the relatively short time available to construct clusters from unlinked genes through genomic rearrangement, and the counterintuitive polarities of some of the substitution effects, we argue that at least some of the clusters must have a pleiotropic basis.