An integrated transcriptome-wide analysis of cave and surface dwelling Astyanax mexicanus

Metabolic shift toward ketosis in asocial cavefish increases social-like affinity

BACKGROUND

Social affinity and collective behavior are nearly ubiquitous in the animal kingdom, but many lineages feature evolutionarily asocial species. These solitary species may have evolved to conserve energy in food-sparse environments. However, the mechanism by which metabolic shifts regulate social affinity is not well investigated.

RESULTS

In this study, we used the Mexican tetra (Astyanax mexicanus), which features riverine sighted surface (surface fish) and cave-dwelling populations (cavefish), to address the impact of metabolic shifts on asociality and other cave-associated behaviors in cavefish, including repetitive turning, sleeplessness, swimming longer distances, and enhanced foraging behavior. After 1 month of ketosis-inducing ketogenic diet feeding, asocial cavefish exhibited significantly higher social affinity, whereas social affinity regressed in cavefish fed the standard diet. The ketogenic diet also reduced repetitive turning and swimming in cavefish. No major behavioral shifts were found regarding sleeplessness and foraging behavior, suggesting that other evolved behaviors are not largely regulated by ketosis. We further examined the effects of the ketogenic diet via supplementation with exogenous ketone bodies, revealing that ketone bodies are pivotal molecules positively associated with social affinity.

CONCLUSIONS

Our study indicated that fish that evolved to be asocial remain capable of exhibiting social affinity under ketosis, possibly linking the seasonal food availability and sociality.

Phylogeographic relationships and morphological evolution between cave and surface Astyanax mexicanus populations (De Filippi 1853) (Actinopterygii, Characidae)

The Astyanax mexicanus complex includes two different morphs, a surface- and a cave-adapted ecotype, found at three mountain ranges in Northeastern Mexico: Sierra de El Abra, Sierra de Guatemala and Sierra de la Colmena (Micos). Since their discovery, multiple studies have attempted to characterize the timing and the number of events that gave rise to the evolution of these cave-adapted ecotypes. Here, using RADseq and genome-wide sequencing, we assessed the phylogenetic relationships, genetic structure and gene flow events between the cave and surface Astyanax mexicanus populations, to estimate the tempo and mode of evolution of the cave-adapted ecotypes. We also evaluated the body shape evolution across different cave lineages using geometric morphometrics to examine the role of phylogenetic signal versus environmental pressures. We found strong evidence of parallel evolution of cave-adapted ecotypes derived from two separate lineages of surface fish and hypothesize that there may be up to four independent invasions of caves from surface fish. Moreover, a strong congruence between the genetic structure and geographic distribution was observed across the cave populations, with the Sierra de Guatemala the region exhibiting most genetic drift among the cave populations analysed. Interestingly, we found no evidence of phylogenetic signal in body shape evolution, but we found support for parallel evolution in body shape across independent cave lineages, with cavefish from the Sierra de El Abra reflecting the most divergent morphology relative to surface and other cavefish populations.

Discovery of putative long non-coding RNAs expressed in the eyes of Astyanax mexicanus (Actinopterygii: Characidae)

Astyanax mexicanus is a well-known model species, that has two morphotypes, cavefish, from subterranean rivers and surface fish, from surface rivers. They are morphologically distinct due to many troglomorphic traits in the cavefish, such as the absence of eyes. Most studies on A. mexicanus are focused on eye development and protein-coding genes involved in the process. However, lncRNAs did not get the same attention and very little is known about them. This study aimed to fill this knowledge gap, identifying, describing, classifying, and annotating lncRNAs expressed in the embryo's eye tissue of cavefish and surface fish. To do so, we constructed a concise workflow to assemble and evaluate transcriptomes, annotate protein-coding genes, ncRNAs families, predict the coding potential, identify putative lncRNAs, map them and predict interactions. This approach resulted in the identification of 33,069 and 19,493 putative lncRNAs respectively mapped in cavefish and surface fish. Thousands of these lncRNAs were annotated and identified as conserved in human and several species of fish. Hundreds of them were validated in silico, through ESTs. We identified lncRNAs associated with genes related to eye development. This is the case of a few lncRNAs associated with sox2, which we suggest being isomorphs of the SOX2-OT, a lncRNA that can regulate the expression of sox2. This work is one of the first studies to focus on the description of lncRNAs in A. mexicanus, highlighting several lncRNA targets and opening an important precedent for future studies focusing on lncRNAs expressed in A. mexicanus.

Novel genome sequence of Chinese cavefish (Triplophysa rosa) reveals pervasive relaxation of natural selection in cavefish genomes

All cavefishes, living exclusively in caves across the globe, exhibit similar phenotypic traits, including the characteristic loss of eyes. To understand whether such phenotypic convergence shares similar genomic bases, here we investigated genome-wide evolutionary signatures of cavefish phenotypes by comparing whole-genome sequences of three pairs of cavefishes and their surface fish relatives. Notably, we newly sequenced and generated a whole-genome assembly of the Chinese cavefish Triplophysa rosa. Our comparative analyses revealed several shared features of cavefish genome evolution. Cavefishes had lower mutation rates than their surface fish relatives. In contrast, the ratio of nonsynonymous to synonymous substitutions (ω) was significantly elevated in cavefishes compared to in surface fishes, consistent with the relaxation of purifying selection. In addition, cavefish genomes had an increased mutational load, including mutations that alter protein hydrophobicity profiles, which were considered harmful. Interestingly, however, we found no overlap in positively selected genes among different cavefish lineages, indicating that the phenotypic convergence in cavefishes was not caused by positive selection of the same sets of genes. Analyses of previously identified candidate genes associated with cave phenotypes supported this conclusion. Genes belonging to the lipid metabolism functional ontology were under relaxed purifying selection in all cavefish genomes, which may be associated with the nutrient-poor habitat of cavefishes. Our work reveals previously uncharacterized patterns of cavefish genome evolution and provides comparative insights into the evolution of cave-associated phenotypic traits.

Fish out of water: Genomic insights into persistence of rainbowfish populations in the desert

How populations of aquatic fauna persist in extreme desert environments is an enigma. Individuals often breed and disperse during favorable conditions. Theory predicts that adaptive capacity should be low in small populations, such as in desert fishes. We integrated satellite-derived surface water data and population genomic diversity from 20,294 single-nucleotide polymorphisms across 344 individuals to understand metapopulation persistence of the desert rainbowfish (Melanotaenia splendida tatei) in central Australia. Desert rainbowfish showed very small effective population sizes, especially at peripheral populations, and low connectivity between river catchments. Yet, there was no evidence of population-level inbreeding and a signal of possible adaptive divergence associated with aridity was detected. Candidate genes for local adaptation included functions related to environmental cues and stressful conditions. Eco-evolutionary modeling showed that positive selection in refugial subpopulations combined with connectivity during flood periods can enable retention of adaptive diversity. Our study suggests that adaptive variation can be maintained in small populations and integrate with neutral metapopulation processes to allow persistence in the desert.

Evolution in Sinocyclocheilus cavefish is marked by rate shifts, reversals, and origin of novel traits

Background

Natural model systems are indispensable for exploring adaptations in response to environmental pressures. Sinocyclocheilus of China, the most diverse cavefish clade in the world (75 species), provide unique opportunities to understand recurrent evolution of stereotypic traits (such as eye loss and sensory expansion) in the context of a deep and diverse phylogenetic group. However, they remain poorly understood in terms of their morphological evolution. Therefore, we explore key patterns of morphological evolution, habitat utilization and geographic distribution in these fishes.

Results

We constructed phylogenies and categorized 49 species based on eye-related condition (Blind, Micro-eyed, and Normal-eyed), habitat types (Troglobitic—cave-restricted; Troglophilic—cave-associated; Surface—outside caves) and existence of horns. Geometric-morphometric analyses show Normal-eyed morphs with fusiform shapes segregating from Blind/Micro-eyed deeper bodied morphs along the first principal-component axis; second axis accounts for shape complexity related to horns. The body shapes showed a significant association with eye-related condition and horn, but not habitat types. Ancestral reconstructions suggest at least three independent origins of Blind morphs, each with different levels of modification in relation to their ancestral Normal-eyed morphs; Sinocyclocheilus are also pre-adapted for cave dwelling. Our geophylogeny shows an east-to-west diversification spanning Pliocene and Pleistocene, with early-diversifying Troglobitic species dominating subterranean habitats of karstic plains whereas predominantly Surface forms inhabit hills to the west. Evolutionary rates analyses suggest that lineages leading to Blind morphs were characterized by significant rate shifts, such as a slowdown in body size evolution and a 5–20 fold increase in rate of eye regression, possibly explained by limited resource availability. Body size and eye size have undergone reversals, but not horns, a trait entailing considerable time to form.

Conclusions

Sinocyclocheilus occupied cave habitats in response to drying associated with aridification of China during late Miocene and the Pliocene. The prominent cave-adaptations (eye-regression, horn-evolution) occur in clades associated with the extensive subterranean cave system in Guangxi and Guizhou provinces. Integration of morphology, phylogeny, rate analyses, molecular-dating and distribution show not only several remarkable patterns of evolution, but also interesting exceptions to these patterns signifying the diversification of Sinocyclocheilus as an invaluable model system to explore evolutionary novelty.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01776-y.

Phylogenetic analysis and expression differences of eye-related genes in cavefish genus Sinocyclocheilus

The adaptive evolution of visual systems has been observed in many cavefish. However, little is known about the molecular mechanisms underlying these adaptations, which include regressive changes such as eye degeneration. Here, we analyzed phylogenetic and expression patterns of 6 eye-related genes (crx, foxg1b, opn1sw2, otx2, rho and sox2) in 12 Sinocyclocheilus species from China, including 8 stygobionts and 4 stygophiles, and examined photoreceptor cell morphology of these species. Those eye-degenerated species of Sinocyclocheilus were polyphyletic and showed different degrees of photoreceptor defects in responses to cave environments. The eye loss and degeneration are the result of convergent evolution. Although S. anophthalmus grouped with the eye-normal species, it displayed not only a high degree of eye degeneration but also significant expression differences in eye-related genes compared with the eye-normal species. The gene foxg1b, which was determined to be under positive selection, might play an important role in the process of eye degeneration in S. anophthalmus based on differential expression. Eye-related gene expression and selection may have contributed to the polyphyly of the cave species. We examined gene expression and duplication in 6 eye-related genes and revealed that these genes displayed considerable diversity in relative expression in Sinocyclocheilus fishes. Otx2 and sox2 were significantly up-regulated in individual cave species, while the other 4 genes (crx, foxg1b, opn1sw2 and rho) were significantly down-regulated. These findings provide a valuable resource for elucidating molecular mechanisms associated with visual system evolution in cavefish.

Genomic Analysis of the Only Blind Cichlid Reveals Extensive Inactivation in Eye and Pigment Formation Genes

Trait loss represents an intriguing evolutionary problem, particularly when it occurs across independent lineages. Fishes in light-poor environments often evolve “troglomorphic” traits, including reduction or loss of both pigment and eyes. Here, we investigate the genomic basis of trait loss in a blind and depigmented African cichlid, Lamprologus lethops, and explore evolutionary forces (selection and drift) that may have contributed to these losses. This species, the only known blind cichlid, is endemic to the lower Congo River. Available evidence suggests that it inhabits deep, low-light habitats. Using genome sequencing, we show that genes related to eye formation and pigmentation, as well as other traits associated with troglomorphism, accumulated inactivating mutations rapidly after speciation. A number of the genes affected in L. lethops are also implicated in troglomorphic phenotypes in Mexican cavefish (Astyanax mexicanus) and other species. Analysis of heterozygosity patterns across the genome indicates that L. lethops underwent a significant population bottleneck roughly 1 Ma, after which effective population sizes remained low. Branch-length tests on a subset of genes with inactivating mutations show little evidence of directional selection; however, low overall heterozygosity may reduce statistical power to detect such signals. Overall, genome-wide patterns suggest that accelerated genetic drift from a severe bottleneck, perhaps aided by directional selection for the loss of physiologically expensive traits, caused inactivating mutations to fix rapidly in this species.

Dark-rearing uncovers novel gene expression patterns in an obligate cave-dwelling fish

Extreme environments often result in the evolution of dramatic adaptive features. The Mexican tetra, Astyanax mexicanus, includes 30 different populations of cave-dwelling forms that live in perpetual darkness. As a consequence, many populations have evolved eye loss, reduced pigmentation, and amplification of nonvisual sensory systems. Closely-related surface-dwelling morphs demonstrate typical vision, pigmentation, and sensation. Transcriptomic assessments in this system have revealed important developmental changes associated with the cave morph, however, they have not accounted for photic rearing conditions. Prior studies reared individuals under a 12:12 hr light/dark (LD) cycle. Here, we reared cavefish under constant darkness (DD) for 5+ years. From these experimental individuals, we performed mRNA sequencing and compared gene expression of surface fish reared under LD conditions to cavefish reared under DD conditions to identify photic-dependent gene expression differences. Gene Ontology enrichment analyses revealed a number of previously underappreciated cave-associated changes impacting blood physiology and olfaction. We further evaluated the position of differentially expressed genes relative to QTL positions from prior studies and found several candidate genes associated with these ecologically relevant lighting conditions. In sum, this work highlights photic conditions as a key environmental factor impacting gene expression patterns in blind cave-dwelling fish.

Unique transcriptional signatures of sleep loss across independently evolved cavefish populations

Animals respond to sleep loss with compensatory rebound sleep, and this is thought to be critical for the maintenance of physiological homeostasis. Sleep duration varies dramatically across animal species, but it is not known whether evolutionary differences in sleep duration are associated with differences in sleep homeostasis. The Mexican cavefish, Astyanax mexicanus, has emerged as a powerful model for studying the evolution of sleep. While eyed surface populations of A. mexicanus sleep approximately 8 hr each day, multiple blind cavefish populations have converged on sleep patterns that total as little as 2 hr each day, providing the opportunity to examine whether the evolution of sleep loss is accompanied by changes in sleep homeostasis. Here, we examine the behavioral and molecular response to sleep deprivation across four independent populations of A. mexicanus. Our behavioral analysis indicates that surface fish and all three cavefish populations display robust recovery sleep during the day following nighttime sleep deprivation, suggesting sleep homeostasis remains intact in cavefish. We profiled transcriptome-wide changes associated with sleep deprivation in surface fish and cavefish. While the total number of differentially expressed genes was not greater for the surface population, the surface population exhibited the highest number of uniquely differentially expressed genes than any other population. Strikingly, a majority of the differentially expressed genes are unique to individual cave populations, suggesting unique expression responses are exhibited across independently evolved cavefish populations. Together, these findings suggest sleep homeostasis is intact in cavefish despite a dramatic reduction in overall sleep duration.

Subterranean life: Behavior, metabolic, and some other adaptations of Astyanax cavefish

The ability of fishes to adapt to any aquatic environment seems limitless. It is enthralling how new species keep appearing at the deep sea or in subterranean environments. There are close to 230 known species of cavefishes, still today the best-known cavefish is Astyanax mexicanus, a Characid that has become a model organism, and has been studied and scrutinized since 1936. There are two morphotypes for A. mexicanus, a surface fish and a cavefish. The surface fish lives in central and northeastern Mexico and south of the United States, while the cavefish is endemic to the "Sierra del Abra-Tanchipa region" in northeast Mexico. The extensive genetic and genomic analysis depicts a complex origin for Astyanax cavefish, with multiple cave invasions and persistent gene flow among cave populations. The surface founder population prevails in the same region where the caves are. In this review, we focus on both morphotype's main morphological and physiological differences, but mainly in recent discoveries about behavioral and metabolic adaptations for subterranean life. These traits may not be as obvious as the troglomorphic characteristics, but are key to understand how Astyanax cavefish thrives in this environment of perpetual darkness.

Comparative Transcriptome Profiling of the Loaches Triplophysa bleekeri and Triplophysa rosa Reveals Potential Mechanisms of Eye Degeneration

Eye degeneration is one of the most obvious characteristics of organisms restricted to subterranean habitats. In cavefish, eye degeneration has evolved independently numerous times and each process is associated with different genetic mechanisms. To gain a better understanding of these mechanisms, we compared the eyes of adult individuals of the cave loach Triplophysa rosa and surface loach Triplophysa bleekeri. Compared with the normal eyes of the surface loach, those of the cave loach were found to possess a small abnormal lens and a defective retina containing photoreceptor cells that lack outer segments. Sequencing of the transcriptomes of both species to identify differentially expressed genes (DEGs) and genes under positive selection revealed 4,802 DEGs and 50 genes under positive selection (dN/dS > 1, FDR < 0.1). For cave loaches, we identified one Gene Ontology category related to vision that was significantly enriched in downregulated genes. Specifically, we found that many of the downregulated genes, including pitx3, lim2, crx, gnat2, rx1, rho, prph2, and β|γ-crystallin are associated with lens/retinal development and maintenance. However, compared with those in the surface loach, the lower dS rates but higher dN rates of the protein-coding sequences in T. rosa indicate that changes in amino acid sequences might be involved in the adaptation and visual degeneration of cave loaches. We also found that genes associated with light perception and light-stimulated vision have evolved at higher rates (some genes dN/dS > 1 but FDR > 0.1). Collectively, the findings of this study indicate that the degradation of cavefish vision is probably associated with both gene expression and amino acid changes and provide new insights into the mechanisms underlying the degeneration of cavefish eyes.

Developmental Transcriptomic Analysis of the Cave-Dwelling Crustacean, Asellus aquaticus

Cave animals are a fascinating group of species often demonstrating characteristics including reduced eyes and pigmentation, metabolic efficiency, and enhanced sensory systems. Asellus aquaticus, an isopod crustacean, is an emerging model for cave biology. Cave and surface forms of this species differ in many characteristics, including eye size, pigmentation, and antennal length. Existing resources for this species include a linkage map, mapped regions responsible for eye and pigmentation traits, sequenced adult transcriptomes, and comparative embryological descriptions of the surface and cave forms. Our ultimate goal is to identify genes and mutations responsible for the differences between the cave and surface forms. To advance this goal, we decided to use a transcriptomic approach. Because many of these changes first appear during embryonic development, we sequenced embryonic transcriptomes of cave, surface, and hybrid individuals at the stage when eyes and pigment become evident in the surface form. We generated a cave, a surface, a hybrid, and an integrated transcriptome to identify differentially expressed genes in the cave and surface forms. Additionally, we identified genes with allele-specific expression in hybrid individuals. These embryonic transcriptomes are an important resource to assist in our ultimate goal of determining the genetic underpinnings of the divergence between the cave and surface forms.

Expression of extraocular opsin genes and light-dependent basal activity of blind cavefish

Background

Animals living in well-lit environments utilize optical stimuli for detecting visual information, regulating the homeostatic pacemaker, and controlling patterns of body pigmentation. In contrast, many subterranean animal species without optical stimuli have evolved regressed binocular eyes and body pigmentation. Interestingly, some fossorial and cave-dwelling animals with regressed eyes still respond to light. These light-dependent responses may be simply evolutionary residuals or they may be adaptive, where negative phototaxis provides avoidance of predator-rich surface environments. However, the relationship between these non-ocular light responses and the underlying light-sensing Opsin proteins has not been fully elucidated.

Methods

To highlight the potential functions of opsins in a blind subterranean animal, we used the Mexican cave tetra to investigate opsin gene expression in the eyes and several brain regions of both surface and cave-dwelling adults. We performed database surveys, expression analyses by quantitative reverse transcription PCR (RT-qPCR), and light-dependent locomotor activity analysis using pinealectomized fish, one of the high-opsin expressing organs of cavefish.

Results

Based on conservative criteria, we identified 33 opsin genes in the cavefish genome. Surveys of available RNAseq data found 26 of these expressed in the surface fish eye as compared to 24 expressed in cavefish extraocular tissues, 20 of which were expressed in the brain. RT-qPCR of 26 opsins in surface and cavefish eye and brain tissues showed the highest opsin-expressing tissue in cavefish was the pineal organ, which expressed exo-rhodopsin at 72.7% of the expression levels in surface fish pineal. However, a pinealectomy resulted in no change to the light-dependent locomotor activity in juvenile cavefish and surface fish. Therefore, we conclude that, after 20,000 or more years of evolution in darkness, cavefish light-dependent basal activity is regulated by a non-pineal extraocular organ.

Expression analyses of cave mollies (Poecilia mexicana) reveal key genes involved in the early evolution of eye regression

Eye regression occurs across cave-dwelling populations of many species and is often coupled with a decrease or loss in eye function. Teleost fishes are among the few vertebrates to undergo widespread colonization of caves and often exhibit eye regression with blindness. Cave populations of the poeciliid fish Poecilia mexicana (cave molly) exhibit reduced-albeit functional-eyes, offering the opportunity to investigate partial eye regression. We sequenced eye transcriptomes of cave and surface populations of P. mexicana to identify differentially expressed genes that potentially underlie eye regression in cave mollies. We identified 28 significantly differentially expressed genes, 20 of which were directly related to light sensitivity, eye structure and visual signaling. Twenty-six of these genes were downregulated in cave compared to surface populations. Functional enrichment analysis revealed eye-related gene ontologies that were under-represented in cave mollies. In addition, a set of co-expressed genes related to vision and circadian rhythm was correlated with habitat type (cave versus surface). Our study suggests that differential gene expression plays a key role in the beginning evolutionary stages of eye regression in P. mexicana, shedding further light on regressive evolution in cavefish.

Stable transgenesis in Astyanax mexicanus using the Tol2 transposase system

BACKGROUND

Astyanax mexicanus is a well-established fish model system for evolutionary and developmental biology research. These fish exist as surface forms that inhabit rivers and 30 different populations of cavefish. Despite important progress in the deployment of new technologies, deep mechanistic insights into the genetic basis of evolution, development, and behavior have been limited by a lack of transgenic lines commonly used in genetic model systems.

RESULTS

Here, we expand the toolkit of transgenesis by characterizing two novel stable transgenic lines that were generated using the highly efficient Tol2 system, commonly used to generate transgenic zebrafish. A stable transgenic line consisting of the zebrafish ubiquitin promoter expresses enhanced green fluorescent protein ubiquitously throughout development in a surface population of Astyanax. To define specific cell-types, a Cntnap2-mCherry construct labels lateral line mechanosensory neurons in zebrafish. Strikingly, both constructs appear to label the predicted cell types, suggesting many genetic tools and defined promoter regions in zebrafish are directly transferrable to cavefish.

CONCLUSION

The lines provide proof-of-principle for the application of Tol2 transgenic technology in A. mexicanus. Expansion on these initial transgenic lines will provide a platform to address broadly important problems in the quest to bridge the genotype-phenotype gap.

The Kingdom of the Blind: Disentangling Fundamental Drivers in the Evolution of Eye Loss

Light is a fundamentally important biological cue used by almost every animal on earth, to maintain daily rhythms, navigate, forage, find mates, or avoid predators. But an enormous number of species live in darkness: in subterranean caves, deep oceans, underground burrows, and within parasitic host bodies, and the loss of eyes appears consistently across these ecosystems. However, the evolutionary mechanisms that lead to the reduction of the visual system remain the subject of great interest and debate more than 150 years after Darwin tackled the issue. Studies of model taxa have discovered significant roles for natural selection, neutral evolution, and pleiotropy, but the interplay between them remains unclear. To nail down unifying concepts surrounding the evolution of eye loss, we must embrace the enormous range of affected animals and habitats. The fine developmental details of model systems such as the Mexican cave tetra Astyanax mexicanus have transformed and enriched the field, but these should be complemented by wider studies to identify truly overarching patterns that apply throughout animals. Here, the major evolutionary drivers are placed within a conceptual cost-benefit framework that incorporates the fundamental constraints and forces that influence evolution in the dark. Major physiological, ecological, and environmental factors are considered within the context of this framework, which appears faithful to observed patterns in deep-sea and cavernicolous animals. To test evolutionary hypotheses, a comparative phylogenetic approach is recommended, with the goal of studying large groups exhibiting repeated reduction, and then comparing these across habitats, taxa, and lifestyles. Currently, developmental and physiological methods cannot feasibly be used on such large scales, but penetrative imaging techniques could provide detailed morphological data non-invasively and economically for large numbers of species. Comprehensive structural datasets can then be contextualized phylogenetically to examine recurrent trends and associations, and to reconstruct character histories through multiple independent transitions into darkness. By assessing these evolutionary trajectories within an energetic cost-benefit framework, the relationships between fundamental influences can be inferred and compared across different biological and physical parameters. However, substantial numbers of biological and environmental factors affect the evolutionary trajectory of loss, and it is critical that researchers make fair and reasonable comparisons between objectively similar groups.

Phototransduction Gene Expression and Evolution in Cave and Surface Crayfishes

In the absence of light in caves, animals have repeatedly evolved reduced eyes and visual systems. Whether the underlying genetic components remain intact in blind species remains unanswered across taxa. The freshwater crayfish have evolved to live in caves multiple times throughout their history; therefore, this system provides an opportunity to probe the genetic patterns and processes underlying repeated vision loss. Using transcriptomic data from the eyes of 14 species of cave and surface crayfishes, we identify the expression of 17 genes putatively related to visual phototransduction. We find a similarly complete repertoire of phototransduction gene families expressed in cave and surface species, but that the expression levels of those transcripts are consistently lower in cave species. We find statistical support for episodic positive selection, increased and decreased selection strength in caves, depending on the gene family. Analyses of gene expression evolution suggest convergent and possibly adaptive downregulation of these genes across eye-reduction events. Our results reveal a combination of evolutionary processes acting on the sequences and gene expression levels of vision-related genes underlying the loss of vision in caves.

Raising the Mexican Tetra Astyanax mexicanus for Analysis of Post-larval Phenotypes and Whole-mount Immunohistochemistry

River and cave-adapted populations of Astyanax mexicanus show differences in morphology, physiology, and behavior. Research focused on comparing adult forms has revealed the genetic basis of some of these differences. Less is known about how the populations differ at post-larval stages (at the onset of feeding). Such studies may provide insight into how cavefish survive through adulthood in their natural environment. Methods for comparing post-larval development in the laboratory require standardized aquaculture and feeding regimes. Here we describe how to raise fish on a diet of nutrient-rich rotifers in non-recirculating water for up to two-weeks post fertilization. We demonstrate how to collect post-larval fish from this nursery system and perform whole-mount immunostaining. Immunostaining is an attractive alternative to transgene expression analysis for investigating development and gene function in A. mexicanus. The nursery method can also be used as a standard protocol for establishing density-matched populations for growth into adults.

Microbiome differences between river-dwelling and cave-adapted populations of the fish Astyanax mexicanus (De Filippi, 1853)

Symbiotic relationships between host and microbiome can play a major role in local adaptation. Previous studies with freshwater organisms have shown that microbiome performs numerous important biochemical functions for the host, playing a key role in metabolism, physiology or health. Experimental studies in fish groups have found an effect of enzymatic activity of gut microbiota on a variety of metabolic processes. The goal of this study was to compare stomach microbiome from cave and surface Astyanax mexicanus, in order to evaluate the potential response of microbiota to contrasting environmental conditions and physiological adaptations of the host. Stomach microbiota was obtained from three different populations: Pachón cave, and two surface rivers (Rascón and Micos rivers). The stomach microbiome was analyzed using the Ion 16S Metagenomic kit considering seven variable regions: V2, V3, V4, V6-7, V8 and V9. A high diversity was observed across samples, including 16 phyla, 120 families and 178 genera. Gammaproteobacteria, Firmicutes, Bacteroidetes and Betaproteobacteria were the most abundant phyla across the samples. Although the relative abundance of the core OTUs at genus level were highly contrasting among populations, we did not recover differences in stomach microbiome between contrasting habitats (cave vs. surface rivers). Rather, we observed a consistent association between β-diversity and dissolved oxygen concentration in water. Therefore, and unexpectedly, the microbiota of A. mexicanus is not linked with the contrasting conditions of the habitat considered here but is related to water parameters.

Genetic analysis reveals candidate genes for activity QTL in the blind Mexican tetra, Astyanax mexicanus

Animal models provide useful tools for exploring the genetic basis of morphological, physiological and behavioral phenotypes. Cave-adapted species are particularly powerful models for a broad array of phenotypic changes with evolutionary, developmental and clinical relevance. Here, we explored the genetic underpinnings of previously characterized differences in locomotor activity patterns between the surface-dwelling and Pachón cave-dwelling populations of Astyanax mexicanus. We identified multiple novel QTL underlying patterns in overall levels of activity (velocity), as well as spatial tank use (time spent near the top or bottom of the tank). Further, we demonstrated that different regions of the genome mediate distinct patterns in velocity and tank usage. We interrogated eight genomic intervals underlying these activity QTL distributed across six linkage groups. In addition, we employed transcriptomic data and draft genomic resources to generate and evaluate a list of 36 potential candidate genes. Interestingly, our data support the candidacy of a number of genes, but do not suggest that differences in the patterns of behavior observed here are the result of alterations to certain candidate genes described in other species (e.g., teleost multiple tissue opsins, melanopsins or members of the core circadian clockwork). This study expands our knowledge of the genetic architecture underlying activity differences in surface and cavefish. Future studies will help define the role of specific genes in shaping complex behavioral phenotypes in Astyanax and other vertebrate taxa.

Characterization and comparison of activity profiles exhibited by the cave and surface morphotypes of the blind Mexican tetra, Astyanax mexicanus

Departure from normal circadian rhythmicity and exposure to atypical lighting cues has been shown to adversely affect human health and wellness in a variety of ways. In contrast, adaptation to extreme environments has led many species to alter or even entirely abandon their reliance upon cyclic environmental inputs, principally daily cycles of light and darkness. The extreme darkness, stability and isolation of cave ecosystems has made cave-adapted species particularly attractive systems in which to study the consequences of life without light and the strategies that allow species to survive and even thrive in such environments. In order to further explore these questions, we have assessed the rhythmicity of locomotion in the blind Mexican tetra, Astyanax mexicanus, under controlled laboratory conditions. Using high-resolution video tracking assays, we characterized patterns in locomotor activity and spatial tank usage for members of the surface and Pachón cave populations. Here we demonstrate that cavefish have a higher overall level of activity and use the space within the trial tank differently than surface fish. Further, Pachón cavefish show circadian rhythmicity in both activity and spatial tank usage under a 12:12 light/dark cycle. We provide further evidence that these cavefish retain a weakly light-entrainable, endogenous circadian oscillator with limited capability to sustain rhythms in activity, but not spatial tank usage, in the absence of photic cues. Finally, we demonstrate a putative behavioral "masking effect" contributing to behavioral rhythms and provide evidence that exposure to constant darkness during development may alter behavioral patterns later in life.

The aquatic animals' transcriptome resource for comparative functional analysis

Background

Aquatic animals have great economic and ecological importance. Among them, non-model organisms have been studied regarding eco-toxicity, stress biology, and environmental adaptation. Due to recent advances in next-generation sequencing techniques, large amounts of RNA-seq data for aquatic animals are publicly available. However, currently there is no comprehensive resource exist for the analysis, unification, and integration of these datasets. This study utilizes computational approaches to build a new resource of transcriptomic maps for aquatic animals. This aquatic animal transcriptome map database dbATM provides de novo assembly of transcriptome, gene annotation and comparative analysis of more than twenty aquatic organisms without draft genome.

Results

To improve the assembly quality, three computational tools (Trinity, Oases and SOAPdenovo-Trans) were employed to enhance individual transcriptome assembly, and CAP3 and CD-HIT-EST software were then used to merge these three assembled transcriptomes. In addition, functional annotation analysis provides valuable clues to gene characteristics, including full-length transcript coding regions, conserved domains, gene ontology and KEGG pathways. Furthermore, all aquatic animal genes are essential for comparative genomics tasks such as constructing homologous gene groups and blast databases and phylogenetic analysis.

Conclusion

In conclusion, we establish a resource for non model organism aquatic animals, which is great economic and ecological importance and provide transcriptomic information including functional annotation and comparative transcriptome analysis. The database is now publically accessible through the URL http://dbATM.mbc.nctu.edu.tw/.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4463-x) contains supplementary material, which is available to authorized users.

Hypocretin underlies the evolution of sleep loss in the Mexican cavefish

The duration of sleep varies dramatically between species, yet little is known about the genetic basis or evolutionary factors driving this variation in behavior. The Mexican cavefish, Astyanax mexicanus, exists as surface populations that inhabit rivers, and multiple cave populations with convergent evolution on sleep loss. The number of Hypocretin/Orexin (HCRT)-positive hypothalamic neurons is increased significantly in cavefish, and HCRT is upregulated at both the transcript and protein levels. Pharmacological or genetic inhibition of HCRT signaling increases sleep in cavefish, suggesting enhanced HCRT signaling underlies the evolution of sleep loss. Ablation of the lateral line or starvation, manipulations that selectively promote sleep in cavefish, inhibit hcrt expression in cavefish while having little effect on surface fish. These findings provide the first evidence of genetic and neuronal changes that contribute to the evolution of sleep loss, and support a conserved role for HCRT in sleep regulation.

Complexities of gene expression patterns in natural populations of an extremophile fish (Poecilia mexicana, Poeciliidae)

Variation in gene expression can provide insights into organismal responses to environmental stress and physiological mechanisms mediating adaptation to habitats with contrasting environmental conditions. We performed an RNA-sequencing experiment to quantify gene expression patterns in fish adapted to habitats with different combinations of environmental stressors, including the presence of toxic hydrogen sulphide (H2 S) and the absence of light in caves. We specifically asked how gene expression varies among populations living in different habitats, whether population differences were consistent among organs, and whether there is evidence for shared expression responses in populations exposed to the same stressors. We analysed organ-specific transcriptome-wide data from four ecotypes of Poecilia mexicana (nonsulphidic surface, sulphidic surface, nonsulphidic cave and sulphidic cave). The majority of variation in gene expression was correlated with organ type, and the presence of specific environmental stressors elicited unique expression differences among organs. Shared patterns of gene expression between populations exposed to the same environmental stressors increased with levels of organismal organization (from transcript to gene to physiological pathway). In addition, shared patterns of gene expression were more common between populations from sulphidic than populations from cave habitats, potentially indicating that physiochemical stressors with clear biochemical consequences can constrain the diversity of adaptive solutions that mitigate their adverse effects. Overall, our analyses provided insights into transcriptional variation in a unique system, in which adaptation to H2 S and darkness coincide. Functional annotations of differentially expressed genes provide a springboard for investigating physiological mechanisms putatively underlying adaptation to extreme environments.

A Comparative Transcriptomic Analysis of Development in Two Astyanax Cavefish Populations

Organisms that are isolated into extreme environments often evolve extreme phenotypes. However, global patterns of dynamic gene expression changes that accompany dramatic environmental changes remain largely unknown. The blind Mexican cavefish, Astyanax mexicanus, has evolved a number of severe cave-associated phenotypes including loss of vision and pigmentation, craniofacial bone fusions, increased fat storage, reduced sleep, and amplified nonvisual sensory systems. Interestingly, surface-dwelling forms have repeatedly entered different caves throughout Mexico, providing a natural set of "replicate" instances of cave isolation. These surrogate "ancestral" surface-dwelling forms persist in nearby rivers, enabling direct comparisons to the "derived" cave-dwelling form. We evaluated changes associated with subterranean isolation by measuring differential gene expression in two geographically distinct cave-dwelling populations (Pachón and Tinaja). To understand the impact of these expression changes on development, we performed RNA-sequencing across four critical stages during which troglomorphic traits first appear in cavefish embryos. Gene ontology (GO) studies revealed similar functional profiles evolved in both independent cave lineages. However, enrichment studies indicated that similar GO profiles were occasionally mediated by different genes. Certain "master" regulators, such as Otx2 and Mitf, appear to be important loci for cave adaptation, as remarkably similar patterns of expression were identified in both independent cave lineages. This work reveals that adaptation to an extreme environment, in two distinct cavefish lineages, evolves through a combination of unique and shared gene expression patterns. Shared expression profiles reflect common environmental pressures, while unique expression likely reflects the fact that similar adaptive traits evolve through diverse genetic mechanisms.

Transcriptomic comparison of invasive bigheaded carps (Hypophthalmichthys nobilis and Hypophthalmichthys molitrix) and their hybrids

Bighead carp (Hypophthalmichthys nobilis) and silver carp (Hypophthalmichthys molitrix), collectively called bigheaded carps, are invasive species in the Mississippi River Basin (MRB). Interspecific hybridization between bigheaded carps has been considered rare within their native rivers in China; however, it is prevalent in the MRB. We conducted de novo transcriptome analysis of pure and hybrid bigheaded carps and obtained 40,759 to 51,706 transcripts for pure, F₁ hybrid, and backcross bigheaded carps. The search against protein databases resulted in 20,336–28,133 annotated transcripts (over 50% of the transcriptome) with over 13,000 transcripts mapped to 23 Gene Ontology biological processes and 127 KEGG metabolic pathways. More transcripts were detected in silver carp than in bighead carp; however, comparable numbers of transcripts were annotated. Transcriptomic variation detected between two F₁ hybrids may indicate a potential loss of fitness in hybrids. The neighbor‐joining distance tree constructed using over 2,500 one‐to‐one orthologous sequences suggests transcriptomes could be used to infer the history of introgression and hybridization. Moreover, we detected 24,792 candidate SNPs that can be used to identify different species. The transcriptomes, orthologous sequences, and candidate SNPs obtained in this study should provide further knowledge of interspecific hybridization and introgression.

Evolutionary Genetics of the Cavefish Astyanax mexicanus

Blind and depigmented fish belonging to the species Astyanax mexicanus are outstanding models for evolutionary genetics. During their evolution in the darkness of caves, they have undergone a number of changes at the morphological, physiological, and behavioral levels, but they can still breed with their river-dwelling conspecifics. The fertile hybrids between these two morphotypes allow forward genetic approaches, from the search of quantitative trait loci to the identification of the mutations underlying the evolution of troglomorphism. We review here the past 30years of evolutionary genetics on Astyanax: from the first crosses and the discovery of convergent evolution of different Astyanax cavefish populations to the most recent evolutionary transcriptomics and genomics studies that have provided researchers with potential candidate genes to be tested using functional genetic approaches. Although significant progress has been made and some genes have been identified, cavefish have not yet fully revealed the secret of their adaptation to the absence of light. In particular, the genetic determinism of their loss of eyes seems complex and still puzzles researchers. We also discuss future research directions, including searches for the origin of cave alleles and searches for selection genome-wide, as well as the necessary but missing information on the timing of cave colonization by surface fish.

Natural bone fragmentation in the blind cave-dwelling fish, Astyanax mexicanus: candidate gene identification through integrative comparative genomics

Animals that colonize dark and nutrient-poor subterranean environments evolve numerous extreme phenotypes. These include dramatic changes to the craniofacial complex, many of which are under genetic control. These phenotypes can demonstrate asymmetric genetic signals wherein a QTL is detected on one side of the face but not the other. The causative gene(s) underlying QTL are difficult to identify with limited genomic resources. We approached this task by searching for candidate genes mediating fragmentation of the third suborbital bone (SO3) directly inferior to the orbit of the eye. We integrated positional genomic information using emerging Astyanax resources, and linked these intervals to homologous (syntenic) regions of the Danio rerio genome. We identified a discrete, approximately 6 Mb, conserved region wherein the gene causing SO3 fragmentation likely resides. We interrogated this interval for genes demonstrating significant differential expression using mRNA-seq analysis of cave and surface morphs across life history. We then assessed genes with known roles in craniofacial evolution and development based on GO term annotation. Finally, we screened coding sequence alterations in this region, identifying two key genes: transforming growth factor β3 (tgfb3) and bone morphogenetic protein 4 (bmp4). Of these candidates, tgfb3 is most promising as it demonstrates significant differential expression across multiple stages of development, maps close (<1 Mb) to the fragmentation critical locus, and is implicated in a variety of other animal systems (including humans) in non-syndromic clefting and malformations of the cranial sutures. Both abnormalities are analogous to the failure-to-fuse phenotype that we observe in SO3 fragmentation. This integrative approach will enable discovery of the causative genetic lesions leading to complex craniofacial features analogous to human craniofacial disorders. This work underscores the value of cave-dwelling fish as a powerful evolutionary model of craniofacial disease, and demonstrates the power of integrative system-level studies for informing the genetic basis of craniofacial aberrations in nature.

Alterations in Mc1r gene expression are associated with regressive pigmentation in Astyanax cavefish

Diverse changes in coloration across distant taxa are mediated through alterations in certain highly conserved pigmentation genes. Among these genes, Mc1r is a frequent target for mutation, and many documented alterations involve coding sequence changes. We investigated whether regulatory mutations in Mc1r may also contribute to pigmentation loss in the blind Mexican cavefish, Astyanax mexicanus. This species comprises multiple independent cave populations that have evolved reduced (or absent) melanic pigmentation as a consequence of living in darkness for millions of generations. Among the most salient cave-associated traits, complete absence (albinism) or reduced levels of pigmentation (brown) have long been the focus of degenerative pigmentation research in Astyanax. These two Mendelian traits have been linked to specific coding mutations in Oca2 (albinism) and Mc1r (brown). However, four of the seven caves harboring the brown phenotype exhibit unaffected coding sequences compared to surface fish. Thus, diverse genetic changes involving the same genes likely impact reduced pigmentation among cavefish populations. Using both sequence and expression analyses, we show that certain cave-dwelling populations harboring the brown mutation have substantial alterations to the putative Mc1r cis-regulatory region. Several of these sequence mutations in the Mc1r 5' region were present across multiple, independent cave populations. This study suggests that pigmentation reduction in Astyanax cavefish evolves through a combination of both coding and cis-regulatory mutations. Moreover, this study represents one of the first attempts to identify regulatory alterations linked to regressive changes in cave-dwelling populations of A. mexicanus.

A Transcriptomic Analysis of Cave, Surface, and Hybrid Isopod Crustaceans of the Species Asellus aquaticus

Cave animals, compared to surface-dwelling relatives, tend to have reduced eyes and pigment, longer appendages, and enhanced mechanosensory structures. Pressing questions include how certain cave-related traits are gained and lost, and if they originate through the same or different genetic programs in independent lineages. An excellent system for exploring these questions is the isopod, Asellus aquaticus. This species includes multiple cave and surface populations that have numerous morphological differences between them. A key feature is that hybrids between cave and surface individuals are viable, which enables genetic crosses and linkage analyses. Here, we advance this system by analyzing single animal transcriptomes of Asellus aquaticus. We use high throughput sequencing of non-normalized cDNA derived from the head of a surface-dwelling male, the head of a cave-dwelling male, the head of a hybrid male (produced by crossing a surface individual with a cave individual), and a pooled sample of surface embryos and hatchlings. Assembling reads from surface and cave head RNA pools yielded an integrated transcriptome comprised of 23,984 contigs. Using this integrated assembly as a reference transcriptome, we aligned reads from surface-, cave- and hybrid- head tissue and pooled surface embryos and hatchlings. Our approach identified 742 SNPs and placed four new candidate genes to an existing linkage map for A. aquaticus. In addition, we examined SNPs for allele-specific expression differences in the hybrid individual. All of these resources will facilitate identification of genes and associated changes responsible for cave adaptation in A. aquaticus and, in concert with analyses of other species, will inform our understanding of the evolutionary processes accompanying adaptation to the subterranean environment.

The rise of Astyanax cavefish

Numerous animals have invaded subterranean caverns and evolved remarkably similar features. These features include loss of vision and pigmentation, and gains in nonvisual sensation. This broad convergence echoes smaller-scale convergence, in which members of the same species repeatedly evolve the same cave-associated phenotypes. The blind Mexican tetra of the Sierra de El Abra region of northeastern Mexico has a complex origin, having recurrently colonized subterranean environments through numerous invasions of surface-dwelling fish. These colonizations likely occurred ∼1-5 MYa. Despite evidence of historical and contemporary gene flow between cave and surface forms, the cave-associated phenotype appears to remain quite stable in nature. This model system has provided insight to the mechanisms of phenotypic regression, the genetic basis for constructive trait evolution, and the origin of behavioral novelties. Here, we document the rise of this model system from its discovery by a Mexican surveyor in 1936, to a powerful system for cave biology and contemporary genetic research. The recently sequenced genome provides exciting opportunities for future research, and will help resolve several long-standing biological problems. Developmental Dynamics 244:1031-1038, 2015. © 2015 Wiley Periodicals, Inc.

The developmental transcriptome of contrasting Arctic charr (Salvelinus alpinus) morphs

Species and populations with parallel evolution of specific traits can help illuminate how predictable adaptations and divergence are at the molecular and developmental level. Following the last glacial period, dwarfism and specialized bottom feeding morphology evolved rapidly in several landlocked Arctic charrSalvelinus alpinuspopulations in Iceland.  To study the genetic divergence between small benthic morphs and limnetic morphs, we conducted RNA-sequencing charr embryos at four stages in early development. We studied two stocks with contrasting morphologies: the small benthic (SB) charr from Lake Thingvallavatn and Holar aquaculture (AC) charr.The data reveal significant differences in expression of several biological pathways during charr development. There was also an expression difference between SB- and AC-charr in genes involved in energy metabolism and blood coagulation genes. We confirmed differing expression of five genes in whole embryos with qPCR, includinglysozymeandnatterin-likewhich was previously identified as a fish-toxin of a lectin family that may be a putative immunopeptide. We also verified differential expression of 7 genes in the developing head that associated consistently with benthic v.s.limnetic morphology (studied in 4 morphs). Comparison of single nucleotide polymorphism (SNP) frequencies reveals extensive genetic differentiation between the SB and AC-charr (~1300 with more than 50% frequency difference). Curiously, three derived alleles in the otherwise conserved 12s and 16s mitochondrial ribosomal RNA genes are found in benthic charr.The data implicate multiple genes and molecular pathways in divergence of small benthic charr and/or the response of aquaculture charr to domestication. Functional, genetic and population genetic studies on more freshwater and anadromous populations are needed to confirm the specific loci and mutations relating to specific ecological traits in Arctic charr.

Genome editing using TALENs in blind Mexican Cavefish, Astyanax mexicanus

Astyanax mexicanus, a teleost fish that exists in a river-dwelling surface form and multiple cave-dwelling forms, is an excellent system for studying the genetic basis of evolution. Cavefish populations, which independently evolved from surface fish ancestors multiple times, have evolved a number of morphological and behavioral traits. Quantitative trait loci (QTL) analyses have been performed to identify the genetic basis of many of these traits. These studies, combined with recent sequencing of the genome, provide a unique opportunity to identify candidate genes for these cave-specific traits. However, tools to test the requirement of these genes must be established to evaluate the role of candidate genes in generating cave-specific traits. To address this need, we designed transcription activator-like effector nucleases (TALENs) to target two genes that contain coding changes in cavefish relative to surface fish and map to the same location as QTL for pigmentation, oculocutaneous albinism 2 (oca2) and melanocortin 1 receptor (mc1r). We found that surface fish genes can be mutated using this method. TALEN-induced mutations in oca2 result in mosaic loss of melanin pigmentation visible as albino patches in F0 founder fish, suggesting biallelic gene mutations in F0s and allowing us to evaluate the role of this gene in pigmentation. The pigment cells in the albino patches can produce melanin upon treatment with L-DOPA, behaving similarly to pigment cells in albino cavefish and providing additional evidence that oca2 is the gene within the QTL responsible for albinism in cavefish. This technology has the potential to introduce a powerful tool for studying the role of candidate genes responsible for the evolution of cavefish traits.

The central nervous system transcriptome of the weakly electric brown ghost knifefish (Apteronotus leptorhynchus): de novo assembly, annotation, and proteomics validation

Background

The brown ghost knifefish (Apteronotus leptorhynchus) is a weakly electric teleost fish of particular interest as a versatile model system for a variety of research areas in neuroscience and biology. The comprehensive information available on the neurophysiology and neuroanatomy of this organism has enabled significant advances in such areas as the study of the neural basis of behavior, the development of adult-born neurons in the central nervous system and their involvement in the regeneration of nervous tissue, as well as brain aging and senescence. Despite substantial scientific interest in this species, no genomic resources are currently available.

Results

Here, we report the de novo assembly and annotation of the A. leptorhynchus transcriptome. After evaluating several trimming and transcript reconstruction strategies, de novo assembly using Trinity uncovered 42,459 unique contigs containing at least a partial protein-coding sequence based on alignment to a reference set of known Actinopterygii sequences. As many as 11,847 of these contigs contained full or near-full length protein sequences, providing broad coverage of the proteome. A variety of non-coding RNA sequences were also identified and annotated, including conserved long intergenic non-coding RNA and other long non-coding RNA observed previously to be expressed in adult zebrafish (Danio rerio) brain, as well as a variety of miRNA, snRNA, and snoRNA. Shotgun proteomics confirmed translation of open reading frames from over 2,000 transcripts, including alternative splice variants. Assignment of tandem mass spectra was greatly improved by use of the assembly compared to databases of sequences from closely related organisms. The assembly and raw reads have been deposited at DDBJ/EMBL/GenBank under the accession number GBKR00000000. Tandem mass spectrometry data is available via ProteomeXchange with identifier PXD001285.

Conclusions

Presented here is the first release of an annotated de novo transcriptome assembly from Apteronotus leptorhynchus, providing a broad overview of RNA expressed in central nervous system tissue. The assembly, which includes substantial coverage of a wide variety of both protein coding and non-coding transcripts, will allow the development of better tools to understand the mechanisms underlying unique characteristics of the knifefish model system, such as their tremendous regenerative capacity and negligible brain senescence.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1354-2) contains supplementary material, which is available to authorized users.

Behaviors of cavefish offer insight into developmental evolution

Many developmental processes have evolved through natural selection, yet in only a few cases do we understand if and how a change of developmental process produces a benefit. For example, many studies in evolutionary biology have investigated the developmental mechanisms that lead to novel structures in an animal, but only a few have addressed if these structures actually benefit the animal at the behavioral level of prey hunting and mating. As such, this review discusses an animal's behavior as the integrated functional output of its evolved morphological and physiological traits. Specifically, we focus on recent findings about the blind Mexican cavefish, Astyanax mexicanus, for which clear relationships exist between its physical traits and ecosystem. This species includes two morphotypes: an eyed surface dweller versus many conspecific types of blind cave dwellers, some of which evolved independently; all of the blind subtypes derived from eyed surface dwellers. The blind cavefish evolved under clear selection pressures: food is sparse and darkness is perpetual. Simulating the major aspects of a cave ecosystem in the laboratory is relatively easy, so we can use this species to begin resolving the relationships between evolved traits and selection pressures—relationships which are more complex for other animals models. This review discusses the recent advances in cavefish research that have helped us establish some key relationships between morphological evolution and environmental shifts. Mol. Reprod. Dev. 82: 268–280, 2015. © 2015 Wiley Periodicals, Inc.

A high-density linkage map for Astyanax mexicanus using genotyping-by-sequencing technology

The Mexican tetra, Astyanax mexicanus, is a unique model system consisting of cave-adapted and surface-dwelling morphotypes that diverged >1 million years (My) ago. This remarkable natural experiment has enabled powerful genetic analyses of cave adaptation. Here, we describe the application of next-generation sequencing technology to the creation of a high-density linkage map. Our map comprises more than 2200 markers populating 25 linkage groups constructed from genotypic data generated from a single genotyping-by-sequencing project. We leveraged emergent genomic and transcriptomic resources to anchor hundreds of anonymous Astyanax markers to the genome of the zebrafish (Danio rerio), the most closely related model organism to our study species. This facilitated the identification of 784 distinct connections between our linkage map and the Danio rerio genome, highlighting several regions of conserved genomic architecture between the two species despite ~150 My of divergence. Using a Mendelian cave-associated trait as a proof-of-principle, we successfully recovered the genomic position of the albinism locus near the gene Oca2. Further, our map successfully informed the positions of unplaced Astyanax genomic scaffolds within particular linkage groups. This ability to identify the relative location, orientation, and linear order of unaligned genomic scaffolds will facilitate ongoing efforts to improve on the current early draft and assemble future versions of the Astyanax physical genome. Moreover, this improved linkage map will enable higher-resolution genetic analyses and catalyze the discovery of the genetic basis for cave-associated phenotypes.

Transcriptomics of developing embryos and organs: A raising tool for evo-devo

Comparative transcriptomics has become an important tool for revisiting many evo-devo questions and exploring new ones, and its importance is likely to increase in the near future, partly because RNA-seq data open many new possibilities. The aim of this opinion piece is twofold. In the first section, we discuss the particularities of transcriptomic studies in evo-devo, focusing mainly on RNA-seq data. The preliminary processing steps (getting coding sequences as well as expression levels) are challenging, because many studied species do not have a sequenced genome. The next step (interpreting expression differences) is also challenging, due to several issues with interpreting expression levels in complex tissues, managing developmental stages and species heterochronies, and the problem of conceptualizing expression differences. In the second section, we discuss some past and possible future applications of transcriptomic approaches (using microarray or RNA-seq) to three major themes in evo-devo: the evolution of the developmental toolkit, the genetic and developmental basis for phenotypic changes, and the general rules of the evolution of development. We believe that conceptual and technical tools are necessary in order to fully exploit the richness of multispecies transcriptomic time-series data.

A new model army: Emerging fish models to study the genomics of vertebrate Evo-Devo

Many fields of biology--including vertebrate Evo-Devo research--are facing an explosion of genomic and transcriptomic sequence information and a multitude of fish species are now swimming in this "genomic tsunami." Here, we first give an overview of recent developments in sequencing fish genomes and transcriptomes that identify properties of fish genomes requiring particular attention and propose strategies to overcome common challenges in fish genomics. We suggest that the generation of chromosome-level genome assemblies--for which we introduce the term "chromonome"--should be a key component of genomic investigations in fish because they enable large-scale conserved synteny analyses that inform orthology detection, a process critical for connectivity of genomes. Orthology calls in vertebrates, especially in teleost fish, are complicated by divergent evolution of gene repertoires and functions following two rounds of genome duplication in the ancestor of vertebrates and a third round at the base of teleost fish. Second, using examples of spotted gar, basal teleosts, zebrafish-related cyprinids, cavefish, livebearers, icefish, and lobefin fish, we illustrate how next generation sequencing technologies liberate emerging fish systems from genomic ignorance and transform them into a new model army to answer longstanding questions on the genomic and developmental basis of their biodiversity. Finally, we discuss recent progress in the genetic toolbox for the major fish models for functional analysis, zebrafish, and medaka, that can be transferred to many other fish species to study in vivo the functional effect of evolutionary genomic change as Evo-Devo research enters the postgenomic era.

Next generation phylogeography of cave and surface Astyanax mexicanus

The loss of traits is a commonly observed evolutionary pattern in cave organisms, but due to extensive morphological convergence, inferring relationships between cave and surface populations can be difficult. For instance, Astyanax mexicanus (the blind Mexican cavefish) is thought to have repeatedly lost its eyes following colonization of cave environments, but the number of evolutionarily independent invasions of this species into caves remains unclear. Because of these repeated losses, it has become a model organism for studying the genetic basis of phenotypic trait loss. Here we reconstruct a high-resolution phylogeography for A. mexicanus inferred from both mitochondrial DNA and several thousand single nucleotide polymorphisms. We provide novel insight into the origin of cave populations from the Sabinos and Río Subterráneo caves and present evidence that the Sabinos cave population is part of a unique cave lineage unrelated to other A. mexicanus cave populations. Our results indicate A. mexicanus cave populations have at least four independent origins.

Complex craniofacial changes in blind cave-dwelling fish are mediated by genetically symmetric and asymmetric loci

The genetic regulators of regressive craniofacial morphologies are poorly understood. To shed light on this problem, we examined the freshwater fish Astyanax mexicanus, a species with surface-dwelling and multiple independent eyeless cave-dwelling forms. Changes affecting the skull in cavefish include morphological alterations to the intramembranous circumorbital bones encircling the eye. Many of these modifications, however, have evolved separately from eye loss, such as fragmentation of the third suborbital bone. To understand the genetic architecture of these eye-independent craniofacial alterations, we developed and scored 33 phenotypes in the context of an F2 hybrid mapping pedigree bred from Pachón cavefish and surface fish. We discovered several individuals exhibiting dramatic left-right differences in bone formation, such as extensive fragmentation on the right side only. This observation, along with well-known eye size asymmetry in natural cave-dwelling animals, led us to further evaluate left-right genetic differences for the craniofacial complex. We discovered three phenotypes, inclusive of bone fragmentation and fusion, which demonstrated a directional heritable basis only on one side. Interestingly, the overall areas of affected bones were genetically symmetric. Phenotypic effect plots of these novel craniofacial QTL revealed that cave alleles are associated with abnormal conditions such as bony fusion and fragmentation. Moreover, many linked loci overlapped with other cave-associated traits, suggesting regressive craniofacial changes may evolve through linkage or as antagonistic pleiotropic consequences of cave-associated adaptations. These novel findings illuminate significant craniofacial changes accompanying evolution in complete darkness and reveal complex changes to the skull differentially influenced by genetic changes affecting the left and right sides.

RNA-Seq technology and its application in fish transcriptomics

High-throughput sequencing technologies, also known as next-generation sequencing (NGS) technologies, have revolutionized the way that genomic research is advancing. In addition to the static genome, these state-of-art technologies have been recently exploited to analyze the dynamic transcriptome, and the resulting technology is termed RNA sequencing (RNA-seq). RNA-seq is free from many limitations of other transcriptomic approaches, such as microarray and tag-based sequencing method. Although RNA-seq has only been available for a short time, studies using this method have completely changed our perspective of the breadth and depth of eukaryotic transcriptomes. In terms of the transcriptomics of teleost fishes, both model and non-model species have benefited from the RNA-seq approach and have undergone tremendous advances in the past several years. RNA-seq has helped not only in mapping and annotating fish transcriptome but also in our understanding of many biological processes in fish, such as development, adaptive evolution, host immune response, and stress response. In this review, we first provide an overview of each step of RNA-seq from library construction to the bioinformatic analysis of the data. We then summarize and discuss the recent biological insights obtained from the RNA-seq studies in a variety of fish species.