Evolution of the eye transcriptome under constant darkness in Sinocyclocheilus cavefish

Wall-following - Phylogenetic context of an enhanced behaviour in stygomorphic Sinocyclocheilus (Cypriniformes: Cyprinidae) cavefishes

With 75 known species, the freshwater fish genus Sinocyclocheilus is the largest cavefish radiation in the world and shows multiple adaptations for cave-dwelling (stygomorphic adaptations), which include a range of traits such as eye degeneration (normal-eyed, micro-eyed and eyeless), depigmentation of skin, and in some species, the presence of "horns". Their behavioural adaptations to subterranean environments, however, are poorly understood. Wall-following (WF) behaviour, where an organism remains in close contact with the boundary demarcating its habitat when in the dark, is a peculiar behaviour observed in a wide range of animals and is enhanced in cave dwellers. Hence, we hypothesise that wall-following is also present in Sinocyclocheilus, possibly enhanced in eyeless species compared to eye bearing (normal-/micro-eyed species). Using 13 species representative of Sinocyclocheilus radiation and eye morphs, we designed a series of assays, based on pre-existing methods for Astyanax mexicanus behavioural experiments, to examine wall-following behaviour under three conditions. Our results indicate that eyeless species exhibit significantly enhanced intensities of WF compared to normal-eyed species, with micro-eyed forms demonstrating intermediate intensities in the WF distance. Using a mtDNA based dated phylogeny (chronogram with four clades A-D), we traced the degree of WF of these forms to outline common patterns. We show that the intensity of WF behaviour is higher in the subterranean clades compared to clades dominated by normal-eyed free-living species. We also found that eyeless species are highly sensitive to vibrations, whereas normal-eyed species are the least sensitive. Since WF behaviour is presented to some degree in all Sinocyclocheilus species, and given that these fishes evolved in the late Miocene, we identify this behaviour as being ancestral with WF enhancement related to cave occupation. Results from this diversification-scale study of cavefish behaviour suggest that enhanced wall-following behaviour may be a convergent trait across all stygomorphic lineages.

Eyeless cave-dwelling Leptonetela spiders still rely on light

Subterranean animals living in perpetual darkness may maintain photoresponse. However, the evolutionary processes behind the conflict between eye loss and maintenance of the photoresponse remain largely unknown. We used Leptonetela spiders to investigate the driving forces behind the maintenance of the photoresponse in cave-dwelling spiders. Our behavioral experiments showed that all eyeless/reduced-eyed cave-dwelling species retained photophobic response and that they had substantially decreased survival at cave entrances due to weak drought resistance. The transcriptomic analysis demonstrated that nearly all phototransduction pathway genes were present and that all tested phototransduction pathway genes were subjected to strong functional constraints in cave-dwelling species. Our results suggest that cave-dwelling eyeless spiders still use light and that light detection likely plays a role in avoiding the cave entrance habitat. This study confirms that some eyeless subterranean animals have retained their photosensitivity due to natural selection and provides a case of mismatch between phenotype and genotype or physiological function in a long-term evolutionary process.

Trait Loss in Evolution: What Cavefish Have Taught Us about Mechanisms Underlying Eye Regression

Reduction or complete loss of traits is a common occurrence throughout evolutionary history. In spite of this, numerous questions remain about why and how trait loss has occurred. Cave animals are an excellent system in which these questions can be answered, as multiple traits, including eyes and pigmentation, have been repeatedly reduced or lost across populations of cave species. This review focuses on how the blind Mexican cavefish, Astyanax mexicanus, has been used as a model system for examining the developmental, genetic, and evolutionary mechanisms that underlie eye regression in cave animals. We focus on multiple aspects of how eye regression evolved in A. mexicanus, including the developmental and genetic pathways that contribute to eye regression, the effects of the evolution of eye regression on other traits that have also evolved in A. mexicanus, and the evolutionary forces contributing to eye regression. We also discuss what is known about the repeated evolution of eye regression, both across populations of A. mexicanus cavefish and across cave animals more generally. Finally, we offer perspectives on how cavefish can be used in the future to further elucidate mechanisms underlying trait loss using tools and resources that have recently become available.

Ecological Speciation Promoted by Divergent Regulation of Functional Genes Within African Cichlid Fishes

Rapid ecological speciation along depth gradients has taken place repeatedly in freshwater fishes, yet molecular mechanisms facilitating such diversification are typically unclear. In Lake Masoko, an African crater lake, the cichlid Astatotilapia calliptera has diverged into shallow-littoral and deep-benthic ecomorphs with strikingly different jaw structures within the last 1,000 years. Using genome-wide transcriptome data, we explore two major regulatory transcriptional mechanisms, expression and splicing-QTL variants, and examine their contributions to differential gene expression underpinning functional phenotypes. We identified 7,550 genes with significant differential expression between ecomorphs, of which 5.4% were regulated by cis-regulatory expression QTLs, and 9.2% were regulated by cis-regulatory splicing QTLs. We also found strong signals of divergent selection on differentially expressed genes associated with craniofacial development. These results suggest that large-scale transcriptome modification plays an important role during early-stage speciation. We conclude that regulatory variants are important targets of selection driving ecologically relevant divergence in gene expression during adaptive diversification.

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.

Quantitative lipidomics and spatial MS-Imaging uncovered neurological and systemic lipid metabolic pathways underlying troglomorphic adaptations in cave-dwelling fish

The Sinocyclocheilus represents a rare, freshwater teleost genus endemic to China that comprises the river-dwelling surface fish and the cave-dwelling cavefish. Using a combinatorial approach of quantitative lipidomics and mass-spectrometry imaging (MSI), we demonstrated that neural compartmentalization of lipid distribution and lipid metabolism are associated with the evolution of troglomorphic traits in Sinocyclocheilus. Attenuated DHA biosynthesis via the Δ4 desaturase pathway led to reductions in docosahexaenoic acid (DHA)-phospholipids in cavefish cerebellum. Instead, cavefish accumulates arachidonic acid (ARA)-phospholipids that may disfavor retinotectal arbor growth. Importantly, MSI of sulfatides, coupled with immunostaining of myelin basic protein and transmission electron microscopy images of hindbrain axons revealed demyelination in cavefish raphe serotonergic neurons. Demyelination in cavefish parallels the loss of neuroplasticity governing social behavior such as aggressive dominance. Outside the brain, quantitative lipidomics and qRT-PCR revealed systemic reductions in membrane esterified DHAs in the liver, attributed to suppression of genes along the Sprecher pathway (elovl2, elovl5, acox1). Development of fatty livers was observed in cavefish, likely mediated by an impeded mobilization of storage lipids, as evident in the diminished expressions of pnpla2, lipea, lipeb, dagla and mgll; and suppressed β-oxidation of fatty acyls via both mitochondria and peroxisomes, reflected in the reduced expressions of cpt1ab, hadhaa, cpt2, decr1 and acox1. These neurological and systemic metabolic adaptations serve to reduce energy expenditure, forming the basis of recessive evolution that eliminates non-essential morphological and behavioral traits, giving cavefish a selective advantage to thrive in caves where proper resource allocation becomes a major determinant of survival.

Evolving in the darkness: Phylogenomics of Sinocyclocheilus cavefishes highlights recent diversification and cryptic diversity

Troglomorphism-any morphological adaptation enabling life to the constant darkness of caves, such as loss of pigment, reduced eyesight or blindness, over-developed tactile and olfactory organs-has long intrigued biologists. However, inferring the proximate and ultimate mechanisms driving the evolution of troglomorphism (stygomorphism) in freshwater fish requires a sound understanding of the evolutionary relationships between surface and stygomorphic lineages. We use Restriction Site Associated DNA Sequencing (RADseq) to better understand the evolution of the Sinocyclocheilus fishes of China. With a remarkable array of derived stygomorphic traits, they comprise the largest cavefish diversification in the world, emerging as a multi-species model system to study evolutionary novelty. We sequenced a total of 120 individuals throughout the Sinocyclocheilus distribution. The data comprised a total of 646,497 bp per individual, including 4378 loci and 67,983 SNPs shared across a minimum of 114 individuals at a given locus. Phylogenetic analyses using either the concatenated RAD loci (RAxML) or the SNPs under a coalescent model (SVDquartets, SNAPP) showed a high degree of congruence with similar topologies and high node support (>95 for most nodes in the phylogeny). The major clades recovered conform to a pattern previously established using Sanger-based mt-DNA sequences, with a few notable exceptions. We now recognize six major clades in this group, elevating the blind cavefish S. tianlinensis and the micro-eyed S. microphthalmus as two new distinct clades due to their deep divergence from other clades. PCA plots of the SNP data also support the recognition of six major clusters of species congruent with the identified clades in ordination space. A Bayes factor delimitation (BFD) analysis showed support for 21 species, recognizing 19 previously described species and two putative new cryptic ones. Two species whose identities were previously disputed, S. furcodorsalis and S. tianeensis, are supported here as distinct species. In addition, our multi-species calibrated tree in SNAPP suggests that the genus Sinocyclocheilus originated around 10.16 Mya, with most speciation events occurring in the last 2 Mya, likely favored by the uplift of the Qinghai-Tibetan Plateau and cave occupation induced by climate-driven aridification during this period. These results provide a firm basis for future comparative studies on the evolution of Sinocyclocheilus and its adaptations to cave life.

Unraveling the heritage of lost traits

We synthesize ontogenetic work spanning the past century that show evolutionarily lost structures are rarely entirely absent from earlier developmental stages. We discuss morphological and genetic insights from developmental studies reveal about the evolution of trait loss and regain.

Widespread retention of ohnologs in key developmental gene families following whole-genome duplication in arachnopulmonates

Whole-genome duplications (WGDs) have occurred multiple times during animal evolution, including in lineages leading to vertebrates, teleosts, horseshoe crabs, and arachnopulmonates. These dramatic events initially produce a wealth of new genetic material, generally followed by extensive gene loss. It appears, however, that developmental genes such as homeobox genes, signaling pathway components and microRNAs are frequently retained as duplicates (so-called ohnologs) following WGD. These not only provide the best evidence for WGD, but an opportunity to study its evolutionary consequences. Although these genes are well studied in the context of vertebrate WGD, similar comparisons across the extant arachnopulmonate orders are patchy. We sequenced embryonic transcriptomes from two spider species and two amblypygid species and surveyed three important gene families, Hox, Wnt, and frizzled, across these and 12 existing transcriptomic and genomic resources for chelicerates. We report extensive retention of putative ohnologs, further supporting the ancestral arachnopulmonate WGD. We also found evidence of consistent evolutionary trajectories in Hox and Wnt gene repertoires across three of the six arachnopulmonate orders, with interorder variation in the retention of specific paralogs. We identified variation between major clades in spiders and are better able to reconstruct the chronology of gene duplications and losses in spiders, amblypygids, and scorpions. These insights shed light on the evolution of the developmental toolkit in arachnopulmonates, highlight the importance of the comparative approach within lineages, and provide substantial new transcriptomic data for future study.

Comparative development and ocular histology between epigean and subterranean salamanders (Eurycea) from central Texas

The salamander clade Eurycea from the karst regions of central Texas provides an ideal platform for comparing divergent nervous and sensory systems since some species exhibit extreme phenotypes thought to be associated with inhabiting a subterranean environment, including highly reduced eyes, while others retain an ancestral ocular phenotype appropriate for life above ground. We describe ocular morphology, comparing three salamander species representing two phenotypes-the surface-dwelling Barton Springs salamander (E. sosorum) and San Marcos salamander (E. nana) and the obligate subterranean Texas blind salamander (E. rathbuni) - in terms of structure and size of their eyes. Eyes were examined using confocal microscopy and measurements were made using ImageJ. Statistical analysis of data was carried out using R. We also provide a developmental series and track eye development and immunolocalization of Pax6 in E. sosorum and E. rathbuni. Adult histology of the surface-dwelling San Marcos salamander (E. nana) shows similarities to E. sosorum. The eyes of adults of the epigean species E. nana and E. sosorum appear fully developed with all the histological features of a fully functional eye. In contrast, the eyes of E. rathbuni adults have fewer layers, lack lenses and other features associated with vision as has been reported previously. However, in early developmental stages eye morphology did not differ significantly between E. rathbuni and E. sosorum. Parallel development is observed between the two phenotypes in terms of morphology; however, Pax6 labeling seems to decrease in the latter stages of development in E.rathbuni. We test for immunolabeling of the visual pigment proteins opsin and rhodopsin and observe immunolocalization around photoreceptor disks in E. nana and E. sosorum, but not in the subterranean E. rathbuni. Our results from examining developing salamanders suggest a combination of underdevelopment and degeneration contribute to the reduced eyes of adult E. rathbuni.

Chromosome-level assembly of southern catfish (silurus meridionalis) provides insights into visual adaptation to nocturnal and benthic lifestyles

The Southern catfish (Silurus meridionalis) is a nocturnal and benthic freshwater fish endemic to the Yangtze River and its tributaries. In this study, we constructed a chromosome-level draft genome of S. meridionalis using 69.7-Gb Nanopore long reads and 49.5-Gb Illumina short reads. The genome assembly was 741.2 Mb in size with a contig N50 of 13.19 Mb. An additional 116.4 Gb of Bionano and 77.4 Gb of Hi-C data were applied to assemble contigs into scaffolds and further into 29 chromosomes, resulting in a 738.9-Mb genome with a scaffold N50 of 28.04 Mb. A total of 22,965 protein-coding genes were predicted from the genome with 22,519 (98.06%) genes functionally annotated. Comparative genomic and transcriptomic analyses revealed a rod-dominated visual system which was responsible for scotopic vision. The absence of cone opsins SWS1 and SWS2 resulted in the lack of ultraviolet and blue violet sensitivity. Mutations at key amino acid sites of RH1.1, RH1.2 and RH2 resulted in spectral tuning good for dim light vision and narrow colour vision. A higher expression level of rod phototransduction genes than that of cone genes and higher rod-to-cone ratio led to higher optical sensitivity under dim light conditions. In addition, analysis of the genes involved in eye morphogenesis and development revealed the loss of some conserved noncoding elements, which might be associated with the small eyes in catfish. Together, our study provides important clues for the adaptation of the catfish visual system to the nocturnal and benthic lifestyles. The draft genome of S. meridionalis represents a valuable resource for studies of the molecular mechanisms of ecological adaptation.

Reproduction, development, asymmetry and late eye regression in the Brazilian cave catfish Ituglanis passensis (Siluriformes, Trichomycteridae): evidence contributing to the neutral mutation theory

The troglobitic (exclusively subterranean source population) catfish Ituglanis passensis (Siluriformes, Trichomycteridae) is endemic to the Passa Três Cave, São Domingos karst area, Rio Tocantins basin, Central Brazil. This unique population presents variably reduced eyes and melanic pigmentation. We describe reproduction and early development in this species based on a spontaneous (non-induced) reproductive-event that occurred in the laboratory in January–February, 2009, while simultaneously comparing with data from the cave-habitat and a previous reproductive event. Egg laying was parceled. Egg-size and number were within variations observed in epigean congeners. Larvae behavior and growth is described. A single surviving specimen was monitored over two years. Eye-regression started late, one year after birth, and followed a pattern of stasis phases intercalated with slow growth and fluctuating asymmetric rates. Late eye regression, associated with asymmetry in eye development and intra-population variability of troglomorphic traits, as shown by several Brazilian subterranean fishes, provide support for the Neutral Mutation Theory.

Comparative transcriptomics reveals the molecular genetic basis of pigmentation loss in Sinocyclocheilus cavefishes

Cave-dwelling animals evolve distinct troglomorphic traits, such as loss of eyes, skin pigmentation, and augmentation of senses following long-term adaptation to perpetual darkness. However, the molecular genetic mechanisms underlying these phenotypic variations remain unclear. In this study, we conducted comparative histology and comparative transcriptomics study of the skin of eight Sinocyclocheilus species (Cypriniformes: Cyprinidae) that included surface- and cave-dwelling species. We analyzed four surface and four cavefish species by using next-generation sequencing, and a total of 802,798,907 clean reads were generated and assembled into 505,495,009 transcripts, which contributed to 1,037,334 unigenes. Bioinformatic comparisons revealed 10,629 and 6,442 significantly differentially expressed unigenes between four different surface-cave fish groups. Further, tens of differentially expressed genes (DEGs) potentially related to skin pigmentation were identified. Most of these DEGs (including GNAQ, PKA, NRAS, and p38) are downregulated in cavefish species. They are involved in key signaling pathways of pigment synthesis, such as the melanogenesis, Wnt, and MAPK pathways. This trend of downregulation was confirmed through qPCR experiments. This study will deepen our understanding of the formation of troglomorphic traits in cavefishes.

Comparative Transcriptomics Reveals the Molecular Genetic Basis of Cave Adaptability in Sinocyclocheilus Fish Species

Cavefish evolved a series of distinct survival mechanisms for adaptation to cave habitat. Such mechanisms include loss of eyesight and pigmentation, sensitive sensory organs, unique dietary preferences, and predation behavior. Thus, it is of great interest to understand the mechanisms underlying these adaptability traits of troglobites. The teleost genus Sinocyclocheilus (Cypriniformes: Cyprinidae) is endemic to China and has more than 70 species reported (including over 30 cavefish species). High species diversity and diverse phenotypes make the Sinocyclocheilus as an outstanding model for studying speciation and adaptive evolution. In this study, we conducted a comparative transcriptomics study on the brain tissues of two Sinocyclocheilus species (surface-dwelling species – Sinocyclocheilus malacopterus and semi-cave-dwelling species – Sinocyclocheilus rhinocerous living in the same water body. A total of 425,188,768 clean reads were generated, which contributed to 102,839 Unigenes. Bioinformatic analysis revealed a total of 3,289 differentially expressed genes (DEGs) between two species Comparing to S. malacopterus, 2,598 and 691 DEGs were found to be respectively, down-regulated and up-regulated in S. rhinocerous. Furthermore, it is also found tens of DEGs related to cave adaptability such as insulin secretion regulation (MafA, MafB, MafK, BRSK, and CDK16) and troglomorphic traits formation (CEP290, nmnat1, coasy, and pqbp1) in the cave-dwelling S. rhinocerous. Interestingly, most of the DEGs were found to be down-regulated in cavefish species and this trend of DEGs expression was confirmed through qPCR experiments. This study would provide an appropriate genetic basis for future studies on the formation of troglomorphic traits and adaptability characters of troglobites, and improve our understanding of mechanisms of cave adaptation.

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.

A Metagenomic Study of Intestinal Microbial Diversity in Relation to Feeding Habits of Surface and Cave-Dwelling Sinocyclocheilus Species

Light is completely absent in cave habitats, causing a shortage or lack of autochthonous photosynthesis. Thus, understanding the mechanisms underlying the ability of organisms to adapt to the unique cave habitat is of great interest. We used high-throughput sequencing of the 16S ribosomal RNA gene of intestinal microorganisms from 11 Sinocyclocheilus (Cypriniformes: Cyprinidae) species, to explore the characteristics of intestinal microorganisms and the adaptive mechanisms of Sinocyclocheilus cavefish and surface fish. We found that the α-diversity and richness of the intestinal microbiome were much higher in cavefish than in surface fish. Principal coordinate analysis showed that cavefish and surface fish formed three clusters because of different dominant gut microorganisms which are generated by different habitats. Based on PICRUSt-predicted functions, harmful substance degradation pathways were much more common in cavefish intestinal microorganisms than in those from surface fish. The intestinal microbiota of surface fish group 1 had a higher capacity for carbohydrate metabolism, whereas protein and amino acid metabolism and digestive pathways were more abundant in microorganisms from the cavefish group and surface fish group 2. Combined analysis of the intestinal microbial composition and functional predictions further revealed the structures and functions of intestinal microbial communities in Sinocyclocheilus cave and surface species. Moreover, based on their habits and intestinal microbial composition and intestinal microbial functional predictions, we inferred that the three fish groups were all omnivorous; however, surface fish group 1 preferred feeding on plants, while surface fish group 2 and cavefish preferred meat. This study improves our understanding of mechanisms of adaptation in cave habitats and may contribute to the protection of these habitats from water pollution.

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 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.

Eye Degeneration and Loss of otx5b Expression in the Cavefish Sinocyclocheilus tileihornes

Cave animals possess remarkable phenotypes associated with existence in their dark environments. The Chinese cavefish Sinocyclocheilus tileihornes shows substantial eye degeneration, a trait shared by most cave species. The extent to which independent evolution of troglomorphic traits uses convergent molecular genetic mechanisms is as yet unknown. We performed transcriptome-wide gene expression profiling in S. tileihornes eyes and compared results with those from the closely related surface species S. angustiporus and an independently derived congeneric cavefish, S. anophthalmus. In total, 52.85 million 100 bp long paired-end clean reads were generated for S. tileihornes, and we identified differentially expressed genes between the three possible pairs of species. Functional analysis of genes differentially expressed between S. tileihornes and S. angustiporus revealed that phototransduction (KEGG id: dre04744) was the most significantly enriched pathway, indicating the obvious differences in response to captured photons between the cavefish S. tileihornes and the surface species S. angustiporus. Analysis of key genes regulating eye development showed complete absence of otx5b (orthodenticle homolog 5) expression in S. tileihornes eyes, probably related to degradation of rods, but normal expression of crx (cone-rod homeobox). The enriched pathways and Otx5 are involved in phototransduction, photoreceptor formation, and regulation of photoreceptor-related gene expression. Unlike the S. tileihornes reported here, S. anophthalmus has reduced crx and otx5 expression. These results show that different species of cavefish within the same genus that independently evolved troglodyte characteristics can have different genetic mechanisms of eye degeneration.

The Evolution of Gene Expression Underlying Vision Loss in Cave Animals

Dissecting the evolutionary genetic processes underlying eye reduction and vision loss in obligate cave-dwelling organisms has been a long-standing challenge in evolutionary biology. Independent vision loss events in related subterranean organisms can provide critical insight into these processes as well as into the nature of convergent loss of complex traits. Advances in evolutionary developmental biology have illuminated the significant role of heritable gene expression variation in the evolution of new forms. Here, we analyze gene expression variation in adult eye tissue across the freshwater crayfish, representing four independent vision-loss events in caves. Species and individual expression patterns cluster by eye function rather than phylogeny, suggesting convergence in transcriptome evolution in independently blind animals. However, this clustering is not greater than what is observed in surface species with conserved eye function after accounting for phylogenetic expectations. Modeling expression evolution suggests that there is a common increase in evolutionary rates in the blind lineages, consistent with a relaxation of selective constraint maintaining optimal expression levels. This is evidence for a repeated loss of expression constraint in the transcriptomes of blind animals and that convergence occurs via a similar trajectory through genetic drift.

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.

Brain of the blind: transcriptomics of the golden-line cavefish brain

The genus Sinocyclocheilus (golden-line barbel) includes 25 species of cave-dwelling blind fish (cavefish) and more than 30 surface-dwelling species with normal vision. Cave environments are dark and generally nutrient-poor with few predators. Cavefish of several genera evolved convergent morphological adaptations in visual, pigmentation, brain, olfactory, and digestive systems. We compared brain morphology and gene expression patterns in a cavefish Sinocyclocheilus anophthalmus with those of a closely related surface-dwelling species S. angustiporus. Results showed that cavefish have a longer olfactory tract and a much smaller optic tectum than surface fish. Transcriptomics by RNA-seq revealed that many genes upregulated in cavefish are related to lysosomes and the degradation and metabolism of proteins, amino acids, and lipids. Genes downregulated in cavefish tended to involve "activation of gene expression in cholesterol biosynthesis" and cholesterol degradation in the brain. Genes encoding Srebfs (sterol regulatory element-binding transcription factors) and Srebf targets, including enzymes in cholesterol synthesis, were downregulated in cavefish brains compared with surface fish brains. The gene encoding Cyp46a1, which eliminates cholesterol from the brain, was also downregulated in cavefish brains, while the total level of cholesterol in the brain remained unchanged. Cavefish brains misexpressed several genes encoding proteins in the hypothalamus-pituitary axis, including Trh, Sst, Crh, Pomc, and Mc4r. These results suggest that the rate of lipid biosynthesis and breakdown may both be depressed in golden-line cavefish brains but that the lysosome recycling rate may be increased in cavefish; properties that might be related to differences in nutrient availability in caves.

Embryonic origin and genetic basis of cave associated phenotypes in the isopod crustacean Asellus aquaticus

Characteristics common to animals living in subterranean environments include the reduction or absence of eyes, lessened pigmentation and enhanced sensory systems. How these characteristics have evolved is poorly understood for the majority of cave dwelling species. In order to understand the evolution of these changes, this study uses an invertebrate model system, the freshwater isopod crustacean, Asellus aquaticus, to examine whether adult differences between cave and surface dwelling individuals first appear during embryonic development. We hypothesized that antennal elaboration, as well as eye reduction and pigment loss, would be apparent during embryonic development. We found that differences in pigmentation, eye formation, and number of segments of antenna II were all present by the end of embryonic development. In addition, we found that cave and surface hatchlings do not significantly differ in the relative size of antenna II and the duration of embryonic development. To investigate whether the regions responsible for eye and pigment differences could be genetically linked to differences in article number, we genotyped F2 hybrids for the four previously mapped genomic regions associated with eye and pigment differences and phenotyped these F2 hybrids for antenna II article number. We found that the region previously known to be responsible for both presence versus absence of pigment and eye size also was significantly associated with article number. Future experiments will address whether pleiotropy and/or genetic linkage play a role in the evolution of cave characteristics in Asellus aquaticus.

modPhEA: model organism Phenotype Enrichment Analysis of eukaryotic gene sets

Motivation

Genome-scale phenotypic data are available for many model organisms, yet existing tools to functionally interpret gene sets from these phenotypic data are largely based on mutagenesis-derived phenotypes observed in mouse or human.

Results

Data from both mutagenesis and knockdown experiments are incorporated into modPhEA to allow users to perform enrichment analyses based on phenotypes observed in budding yeast (Saccharomyces cerevisiae), roundworm (Caenorhabditis elegans), fruit fly (Drosophila melanogaster), zebrafish (Danio rerio), mouse (Mus musculus) and humans (Homo sapiens). The phenotypes analysed can be customized to investigate complex traits and gene sets from any fully sequenced animal or fungal genome are also supported by modPhEA.

Availability and implementation

Freely available on the web at http://evol.nhri.org.tw/modPhEA/.

Contact

liaoby@nhri.org.tw.

Supplementary information

Supplementary data are available at Bioinformatics online.

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.

The transcriptomes of cave and surface populations of Gammarus minus (Crustacea: Amphipoda) provide evidence for positive selection on cave downregulated transcripts

Gammarus minus, a freshwater amphipod living in the cave and surface streams in the eastern USA, is an excellent model for investigating evolutionary adaptation to the subterranean environment. RNA-Seq was conducted on one pair of morphologically distinct sister populations inhabiting surface and cave habitats to identify genes that were differentially expressed in the two populations, as well as to compare levels and patterns of genetic variation within and between populations. Of the 104,630 transcripts identified in the transcriptome assembly, 57% had higher average levels of expression in the cave population. After Benjamini-Hochberg correction for multiple tests, 1517 and 551 transcripts were significantly upregulated or downregulated, respectively, in the cave population, indicating an almost three-fold enrichment of cave-upregulated genes. The average level of nucleotide diversity across all transcripts was significantly lower in the cave population. Within the cave population, where the average nucleotide diversity of cave-downregulated transcripts was 75% that of the cave-upregulated transcripts, a highly significant difference, whereas within the spring population the nucleotide diversities of cave-downregulated and cave-upregulated transcripts was virtually identical. Three lines of evidence suggest that the reduced variation in cave downregulated transcripts is due to positive selection in the cave population: 1) the average neutrality index of cave-downregulated genes was < 1, consistent with positive selection, and significantly less than that of cave-upregulated genes; 2) Tajima’s D was positively correlated with the cave:surface expression ratio, and 3) cave-downregulated transcripts were significantly more likely to be highly diverged from their surface homologs than cave-upregulated transcripts. Five transcripts had fixed premature termination codons in the cave population. The expression patterns and sequence variation in one such transcript, encoding the DNA repair protein photolyase, were examined in more detail and provide the first evidence for the relaxation of functional constraint in this light-dependent protein in a subterranean population.

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.

Exploring the molecular basis of adaptive evolution in hydrothermal vent crab Austinograea alayseae by transcriptome analysis

Elucidating the genetic mechanisms of adaptation to the hydrothermal vent in organisms at genomic level is significant for understanding the adaptive evolution process in the extreme environment. We performed RNA-seq on four different tissues of a vent crab species, Austinograea alayseae, producing 725,461 unigenes and 134,489 annotated genes. Genes related to sensory, circadian rhythm, hormone, hypoxia stress, metal detoxification and immunity were identified. It was noted that in the degenerated eyestalk, transcription of phototransduction related genes which are important for retinal function was greatly reduced; three crucial neuropeptide hormones, one molt-inhibiting and two crustacean hyperglycemic hormone precursors were characterized with conserved domains; hypoxia-inducible factor 1 and two novel isoforms of metallothioneins in the vent crabs were discovered. An analysis of 6,932 orthologs among three crabs A. alayseae, Portunus trituberculutus and Eriocheir sinensis revealed 19 positive selected genes (PSGs). Most of the PSGs were involved in immune responses, such as crustins and anti-lipopolysaccharide factor, suggesting their function in the adaptation to environment. The characterization of the first vent crab transcriptome provides abundant resources for genetic and evolutionary studies of this species, and paves the way for further investigation of vent adaptation process in crabs.

Genome-wide identification of genes probably relevant to the adaptation of schizothoracins (Teleostei: Cypriniformes) to the uplift of the Qinghai-Tibet Plateau

BACKGROUND

Molecular adaptation to the severe environments present during the uplift of the Qinghai-Tibet Plateau has attracted the attention of researchers. The divergence of the three specialization groups of schizothoracins (Primitive, Specialized and Highly Specialized) may correspond to the three phases of plateau uplift. Based on the transcripts of representative species of the three specialized groups and an outgroup, genes in schizothoracins that may have played important roles during the adaptation to new environments were investigated.

RESULTS

The contigs of Gymnodiptychus dybowskii and Schizothorax pseudaksaiensis were compared with those of Gymnocypris przewalskii ganzihonensis and the outgroup Sinocyclocheilus angustiporus, and 5,894 ortholog groups with an alignment length longer than 90 nt after deleting gaps were retained. Evolutionary analyses indicated that the average evolutionary rate of the branch leading to the Specialized group was faster than that of the branch leading to the Highly Specialized group. Moreover, the numbers of gene categories in which more than half of the genes evolved faster than the average values of the genome were 117 and 15 along the branches leading to the Specialized and Highly Specialized groups, respectively. A total of 40, 36, and 55 genes were likely subject to positive selection along the branches leading to the Primitive, Specialized and Highly Specialized groups, respectively, and many of these genes are likely relevant to adaptation to the cold temperatures, low oxygen concentrations, and strong ultraviolet radiation that result from elevation.

CONCLUSIONS

By selecting representative species of the three groups of schizothoracins and applying next-generation sequencing technology, several candidate genes corresponding to adaptation to the three phases of plateau uplift were identified. Some of the genes identified in this report that were likely subject to positive selection are good candidates for subsequent evolutionary and functional analyses of adaptation to high altitude.

Complete mitochondrial genome of a cavefish Sinocyclocheilus ronganensis (Cypriniformes: Cyprinidae)

In this study, we report the complete mitochondrial genome of Sinocyclocheilus ronganensis. The whole mitochondrial genome is 16,587 bp with an accession number KX778473, and consists of 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs genes, and a 936 bp control region. Phylogenetic analysis shows that S. ronganensis is close to cave-restricted S. anophthalmus and surface-dwelling S. grahami. The complete mitogenome of S. ronganensis may provide useful information for studying the genetic mechanism of cavefish, and enrich the fish mitochondrial genome resource for further research.

Cavefish and the basis for eye loss

Animals have colonized the entire world from rather moderate to the harshest environments, some of these so extreme that only few animals are able to survive. Cave environments present such a challenge and obligate cave animals have adapted to perpetual darkness by evolving a multitude of traits. The most common and most studied cave characteristics are the regression of eyes and the overall reduction in pigmentation. Studying these traits can provide important insights into how evolutionary forces drive convergent and regressive adaptation. The blind Mexican cavefish (Astyanax mexicanus) has emerged as a useful model to study cave evolution owing to the availability of genetic and genomic resources, and the amenability of embryonic development as the different populations remain fertile with each other. In this review, we give an overview of our current knowledge underlying the process of regressive and convergent evolution using eye degeneration in cavefish as an example.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Different contributions of local- and distant-regulatory changes to transcriptome divergence between stickleback ecotypes

Differential gene expression can play an important role in phenotypic evolution and divergent adaptation. Although differential gene expression can be caused by both local- and distant-regulatory changes, we know little about their relative contribution to transcriptome evolution in natural populations. Here, we conducted expression quantitative trait loci (eQTL) analysis to investigate the genetic architecture underlying transcriptome divergence between marine and stream ecotypes of threespine sticklebacks (Gasterosteus aculeatus). We identified both local and distant eQTLs, some of which constitute hotspots, regions with a disproportionate number of significant eQTLs relative to the genomic background. The majority of local eQTLs including those in the hotspots caused expression changes consistent with the direction of transcriptomic divergence between ecotypes. Genome scan analysis showed that many local eQTLs overlapped with genomic regions of high differentiation. In contrast, nearly half of the distant eQTLs including those in the hotspots caused opposite expression changes, and few overlapped with regions of high differentiation, indicating that distant eQTLs may act as a constraint of transcriptome evolution. Finally, a comparison between two salinity conditions revealed that nearly half of eQTL hotspots were environment specific, suggesting that analysis of genetic architecture in multiple conditions is essential for predicting response to selection.

Resolving Cypriniformes relationships using an anchored enrichment approach

BACKGROUND

Cypriniformes (minnows, carps, loaches, and suckers) is the largest group of freshwater fishes in the world (~4300 described species). Despite much attention, previous attempts to elucidate relationships using molecular and morphological characters have been incongruent. In this study we present the first phylogenomic analysis using anchored hybrid enrichment for 172 taxa to represent the order (plus three out-group taxa), which is the largest dataset for the order to date (219 loci, 315,288 bp, average locus length of 1011 bp).

RESULTS

Concatenation analysis establishes a robust tree with 97 % of nodes at 100 % bootstrap support. Species tree analysis was highly congruent with the concatenation analysis with only two major differences: monophyly of Cobitoidei and placement of Danionidae.

CONCLUSIONS

Most major clades obtained in prior molecular studies were validated as monophyletic, and we provide robust resolution for the relationships among these clades for the first time. These relationships can be used as a framework for addressing a variety of evolutionary questions (e.g. phylogeography, polyploidization, diversification, trait evolution, comparative genomics) for which Cypriniformes is ideally suited.

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.

"Wrecks of Ancient Life": Genetic Variants Vetted by Natural Selection

The Genetics Society of America’s George W. Beadle Award honors individuals who have made outstanding contributions to the community of genetics researchers and who exemplify the qualities of its namesake as a respected academic, administrator, and public servant. The 2015 recipient is John Postlethwait. He has made groundbreaking contributions in developing the zebrafish as a molecular genetic model and in understanding the evolution of new gene functions in vertebrates. He built the first zebrafish genetic map and showed that its genome, along with that of distantly related teleost fish, had been duplicated. Postlethwait played an integral role in the zebrafish genome-sequencing project and elucidated the genomic organization of several fish species. Postlethwait is also honored for his active involvement with the zebrafish community, advocacy for zebrafish as a model system, and commitment to driving the field forward.

The coiled-coil domain containing protein Ccdc136b antagonizes maternal Wnt/β-catenin activity during zebrafish dorsoventral axial patterning

The coiled-coil domain containing protein CCDC136 is a putative tumor suppressor and significantly down-regulated in gastric and colorectal cancer tissues. However, little is known about its biological functions during vertebrate embryo development. Zebrafish has two CCDC136 orthologs, ccdc136a and ccdc136b, but only ccdc136b is highly expressed during early embryonic development. In this study, we demonstrate that ccdc136b is required for dorsal-ventral axial patterning in zebrafish embryos. ccdc136b morphants display strongly dorsalized phenotypes. Loss- and gain-of-function experiments in zebrafish embryos and mammalian cells show that Ccdc136b is a crucial negative regulator of the Wnt/β-catenin signaling pathway, and plays a critical role in the establishment of the dorsal-ventral axis. We further find that Ccdc136b interacts with APC, promotes the binding affinity of APC with β-catenin and then facilitates the turnover of β-catenin. These results provide the first evidence that CCDC136 regulates zebrafish dorsal-ventral patterning by antagonizing Wnt/β-catenin signal transduction and suggest a potential mechanism underlying its suppressive activity in carcinogenesis.

De novo transcriptome assembly of the grapevine phylloxera allows identification of genes differentially expressed between leaf- and root-feeding forms

Background

Grapevine phylloxera, an insect related to true aphids, is a major historic pest of viticulture only controlled through the selection of resistant rootstocks or through quarantine regulations where grapevine is cultivated own-rooted. Transcriptomic data could help understand the bases of its original life-traits, including a striking case of polyphenism, with forms feeding on roots and forms feeding in leaf-galls. Comparisons with true aphids (for which complete genomes have been sequenced) should also allow to link differences in life-traits of the two groups with changes in gene repertoires or shifts in patterns of expression.

Results

We sequenced transcriptomes of the grapevine phylloxera (Illumina technology), choosing three life-stages (adults on roots or on leaf galls, and eggs) to cover a large catalogue of transcripts, and performed a de novo assembly. This resulted in 105,697 contigs, which were annotated: most contigs had a best blastx hit to the pea aphid (phylogenetically closest complete genome), while very few bacterial hits were recorded (except for Probionibacterium acnes). Coding sequences were predicted from this data set (17,372 sequences), revealing an extremely high AT-bias (at the third codon position). Differential expression (DE) analysis among root-feeding and gall-feeding showed that i) the root-feeding form displayed a much larger number of differentially expressed transcripts ii) root-feeding biased genes were enriched in some categories, for example cuticular proteins and genes associated with cell-cell signaling iii) leaf-galling-biased genes were enriched in genes associated with the nucleus and DNA-replication, suggesting a metabolism more oriented towards fast and active multiplication. We also identified a gene family with a very high expression level (copies totaling nearly 10 % of the reads) in the grapevine phylloxera (both in root and leaf galling forms), but usually expressed at very low levels in true aphids (except in sexual oviparous females). These transcripts thus appear to be associated with oviparity.

Conclusions

Our study illustrated major intraspecific changes in transcriptome profiles, related with different life-styles (and the feeding on roots versus in leaf-galls). At a different scale, we could also illustrate one major shift in expression levels associated with changes in life-traits that occurred along evolution and that respectively characterize (strictly oviparous) grapevine phylloxera and (mostly viviparous) true aphids.

Electronic supplementary material

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

The Sinocyclocheilus cavefish genome provides insights into cave adaptation

BACKGROUND

An emerging cavefish model, the cyprinid genus Sinocyclocheilus, is endemic to the massive southwestern karst area adjacent to the Qinghai-Tibetan Plateau of China. In order to understand whether orogeny influenced the evolution of these species, and how genomes change under isolation, especially in subterranean habitats, we performed whole-genome sequencing and comparative analyses of three species in this genus, S. grahami, S. rhinocerous and S. anshuiensis. These species are surface-dwelling, semi-cave-dwelling and cave-restricted, respectively.

RESULTS

The assembled genome sizes of S. grahami, S. rhinocerous and S. anshuiensis are 1.75 Gb, 1.73 Gb and 1.68 Gb, respectively. Divergence time and population history analyses of these species reveal that their speciation and population dynamics are correlated with the different stages of uplifting of the Qinghai-Tibetan Plateau. We carried out comparative analyses of these genomes and found that many genetic changes, such as gene loss (e.g. opsin genes), pseudogenes (e.g. crystallin genes), mutations (e.g. melanogenesis-related genes), deletions (e.g. scale-related genes) and down-regulation (e.g. circadian rhythm pathway genes), are possibly associated with the regressive features (such as eye degeneration, albinism, rudimentary scales and lack of circadian rhythms), and that some gene expansion (e.g. taste-related transcription factor gene) may point to the constructive features (such as enhanced taste buds) which evolved in these cave fishes.

CONCLUSION

As the first report on cavefish genomes among distinct species in Sinocyclocheilus, our work provides not only insights into genetic mechanisms of cave adaptation, but also represents a fundamental resource for a better understanding of cavefish biology.

The complete mitochondrial genome sequence of a cavefish Sinocyclocheilus anshuiensis (Cypriniformes: Cyprinidae)

Sinocyclocheilus anshuiensis is a special cavefish that lives in the Southwestern China with many specific regressive features, such as rudimentary eyes and scales, and loss of pigmentation. In this study, we performed sequencing and assembly of its complete mitochondrial genome. We confirmed that total length of the mitochondrion is 16 618 bp with an AT ratio of 55.4%. The complete mitochondrial genome contains 13 protein-coding genes, 22 transfer RNAs, 2 ribosomal RNAs and a 963 bp control region. Our current data provide important resources for the research of cavefish mitochondrial evolution and energy metabolism.

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.

Transcriptome-Wide Differential Gene Expression in Bicyclus anynana Butterflies: Female Vision-Related Genes Are More Plastic

Vision is energetically costly to maintain. Consequently, over time many cave-adapted species downregulate the expression of vision genes or even lose their eyes and associated eye genes entirely. Alternatively, organisms that live in fluctuating environments, with different requirements for vision at different times, may evolve phenotypic plasticity for expression of vision genes. Here, we use a global transcriptomic and candidate gene approach to compare gene expression in the heads of a polyphenic butterfly. Bicyclus anynana have two seasonal forms that display sexual dimorphism and plasticity in eye morphology, and female-specific plasticity in opsin gene expression. Nonchoosy dry season females downregulate opsin expression, consistent with the high physiological cost of vision. To identify other genes associated with sexually dimorphic and seasonally plastic differences in vision, we analyzed RNA-sequencing data from whole head tissues. We identified two eye development genes (klarsicht and warts homologs) and an eye pigment biosynthesis gene (henna) differentially expressed between seasonal forms. By comparing sex-specific expression across seasonal forms, we found that klarsicht, warts, henna, and another eye development gene (domeless) were plastic in a female-specific manner. In a male-only analysis, white (w) was differentially expressed between seasonal forms. Reverse transcription polymerase chain reaction confirmed that warts and white are expressed in eyes only, whereas klarsicht, henna and domeless are expressed in both eyes and brain. We find that differential expression of eye development and eye pigment genes is associated with divergent eye phenotypes in B. anynana seasonal forms, and that there is a larger effect of season on female vision-related genes.

Transcriptome analysis revealed changes of multiple genes involved in immunity in Cynoglossus semilaevis during Vibrio anguillarum infection

Half-smooth tongue sole (Cynoglossus semilaevis) is one of the most valuable marine aquatic species in Northern China. Given to the rapid development of aquaculture industry, the C. semilaevis was subjected to disease-causing bacteria Vibrio anguillarum. It therefore is indispensable and urgent to understand the mechanism of C. semilaevis host defense against V. anguillarum infection. In the present study, the extensively analysis at the transcriptome level for V. Anguillarum disease in tongue sole was carried out. In total, 94,716 high quality contigs were generated from 75,884,572 clean reads in three libraries (HOSG, NOSG, and CG). 22,746 unigenes were identified when compared with SwissProt, an NR protein database and NT nucleotide database. 954 genes exhibiting the differentially expression at least one pair of comparison in all three libraries were identified. GO enrichment for these genes revealed gene response to biotic stimulus, immune system regulation, and immune response and cytokine production. Further, the pathways such as complement and coagulation cascades and Vibrio cholerae infection pathways were enriched in defensing of pathogen. Besides, 13,428 SSRs and 118,239 SNPs were detected in tongue sole, providing further support for genetic variation and marker-assisted selection in future. In summary, this study identifies several putative immune pathways and candidate genes deserving further investigation in the context of development of therapeutic regimens and lays the foundation for selecting resistant lines of C. semilaevis against V. anguillarum.