Genetic bias in repeated evolution of pigment loss in cave populations of the Asellus aquaticus species complex

Similar phenotypes can evolve repeatedly under the same evolutionary pressures. A compelling example is the evolution of pigment loss and eye loss in cave-dwelling animals. While specific genomic regions or genes associated with these phenotypes have been identified in model species, it remains uncertain whether a bias towards particular genetic mechanisms exists. An isopod crustacean, Asellus aquaticus, is an ideal model organism to investigate this phenomenon. It inhabits surface freshwaters throughout Europe but has colonized groundwater on multiple independent occasions and evolved several cave populations with distinct ecomorphology. Previous studies have demonstrated that three different cave populations utilized common genetic regions, potentially the same genes, in the evolution of pigment and eye loss. Expanding on this, we conducted analysis on two additional cave populations, distinct either phylogenetically or biogeographically from those previously examined. We generated F2 hybrids from cave × surface crosses and tested phenotype-genotype associations, as well as conducted complementation tests by crossing individuals from different cave populations. Our findings revealed that pigment loss and orange eye pigment in additional cave populations were associated with the same genomic regions as observed in the three previously tested cave populations. Moreover, the lack of complementation across all cross combinations suggests that the same gene likely drives pigment loss. These results substantiate a genetic bias in the recurrent evolution of pigment loss in this model system. Future investigations should focus on the cause behind this bias, possibly arising from allele recruitment from ancestral surface populations' genetic variation or advantageous allele effects via pleiotropy.

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.

Transcriptomic analysis of cave, surface, and hybrid samples of the isopod Asellus aquaticus and identification of chromosomal location of candidate genes for cave phenotype evolution

BACKGROUND

Transcriptomic methods can be used to elucidate genes and pathways responsible for phenotypic differences between populations. Asellus aquaticus is a freshwater isopod crustacean with surface- and cave-dwelling ecomorphs that differ greatly in multiple phenotypes including pigmentation and eye size. Multiple genetic resources have been generated for this species, but the genes and pathways responsible for cave-specific characteristics have not yet been identified. Our goal was to generate transcriptomic resources in tandem with taking advantage of the species' ability to interbreed and generate hybrid individuals.

RESULTS

We generated transcriptomes of the Rakov Škocjan surface population and the Rak Channel of Planina Cave population that combined Illumina short-read assemblies and PacBio Iso-seq long-read sequences. We investigated differential expression at two different embryonic time points as well as allele-specific expression of F₁ hybrids between cave and surface individuals. RNAseq of F₂ hybrids, as well as genotyping of a backcross, allowed for positional information of multiple candidate genes from the differential expression and allele-specific analyses.

CONCLUSIONS

As expected, genes involved in phototransduction and ommochrome synthesis were under-expressed in the cave samples as compared to the surface samples. Allele-specific expression analysis of F₁ hybrids identified genes with cave-biased (cave allele has higher mRNA levels than the surface allele) and surface-biased expression (surface allele has higher mRNA levels than the cave allele). RNAseq of F₂ hybrids allowed for multiple genes to be placed to previously mapped genomic regions responsible for eye and pigmentation phenotypes. In the future, these transcriptomic resources will guide prioritization of candidates for functional analysis.

Standing genetic variation as a potential mechanism of novel cave phenotype evolution in the freshwater isopod, Asellus aquaticus

Novel phenotypes can come about through a variety of mechanisms including standing genetic variation from a founding population. Cave animals are an excellent system in which to study the evolution of novel phenotypes such as loss of pigmentation and eyes. Asellus aquaticus is a freshwater isopod crustacean found in Europe and has both a surface and a cave ecomorph which vary in multiple phenotypic traits. An orange eye phenotype was previously revealed by F₂ crosses and backcrosses to the cave parent within two examined Slovenian cave populations. Complete loss of pigmentation, both in eye and body, is epistatic to the orange eye phenotype and therefore the orange eye phenotype is hidden within the cave populations. Our goal was to investigate the origin of the orange eye alleles within the Slovenian cave populations by examining A. aquaticus individuals from Slovenian and Romanian surface populations and Asellus aquaticus infernus individuals from a Romanian cave population. We found orange eye individuals present in lab raised surface populations of A. aquaticus from both Slovenia and Romania. Using a mapping approach with crosses between individuals of two surface populations, we found that the region known to be responsible for the orange eye phenotype within the two previously examined Slovenian cave populations was also responsible within both the Slovenian and the Romanian surface populations. Complementation crosses between orange eye Slovenian and orange eye Romanian surface individuals suggest that the same gene is responsible for the orange eye phenotype in both surface populations. Additionally, we observed a low frequency phenotype of eye loss in crosses generated between the two surface populations and also in the Romanian surface population. Finally, in a cave population from Romania, A. aquaticus infernus, we found that the same region is also responsible for the orange eye phenotype as the Slovenian cave populations and the Slovenian and Romanian surface populations. Therefore, we present evidence that variation present in the cave populations could originate from standing variation present in the surface populations and/or transgressive hybridization of different surface phylogenetic lineages rather than de novo mutations.

Biogeography and eye size evolution of the ogre-faced spiders

Net-casting spiders (Deinopidae) comprise a charismatic family with an enigmatic evolutionary history. There are 67 described species of deinopids, placed among three genera, Deinopis, Menneus, and Asianopis, that are distributed globally throughout the tropics and subtropics. Deinopis and Asianopis, the ogre-faced spiders, are best known for their giant light-capturing posterior median eyes (PME), whereas Menneus does not have enlarged PMEs. Molecular phylogenetic studies have revealed discordance between morphology and molecular data. We employed a character-rich ultra-conserved element (UCE) dataset and a taxon-rich cytochrome-oxidase I (COI) dataset to reconstruct a genus-level phylogeny of Deinopidae, aiming to investigate the group's historical biogeography, and examine PME size evolution. Although the phylogenetic results support the monophyly of Menneus and the single reduction of PME size in deinopids, these data also show that Deinopis is not monophyletic. Consequently, we formally transfer 24 Deinopis species to Asianopis; the transfers comprise all of the African, Australian, South Pacific, and a subset of Central American and Mexican species. Following the divergence of Eastern and Western deinopids in the Cretaceous, Deinopis/Asianopis dispersed from Africa, through Asia and into Australia with its biogeographic history reflecting separation of Western Gondwana as well as long-distance dispersal events.

Color, activity period, and eye structure in four lineages of ants: Pale, nocturnal species have evolved larger eyes and larger facets than their dark, diurnal congeners

The eyes of insects display an incredible diversity of adaptations to enhance vision across the gamut of light levels that they experience. One commonly studied contrast is the difference in eye structure between nocturnal and diurnal species, with nocturnal species typically having features that enhance eye sensitivity such as larger eyes, larger eye facets, and larger ocelli. In this study, we compared eye structure between workers of closely related nocturnal and diurnal above ground foraging ant species (Hymenoptera: Formicidae) in four genera (Myrmecocystus, Aphaenogaster, Temnothorax, Veromessor). In all four genera, nocturnal species tend to have little cuticular pigment (pale), while diurnal species are heavily pigmented (dark), hence we could use cuticle coloration as a surrogate for activity pattern. Across three genera (Myrmecocystus, Aphaenogaster, Temnothorax), pale species, as expected for nocturnally active animals, had larger eyes, larger facet diameters, and larger visual spans compared to their dark, more day active congeners. This same pattern occurred for one pale species of Veromessor, but not the other. There were no consistent differences between nocturnal and diurnal species in interommatidial angles and eye parameters both within and among genera. Hence, the evolution of eye features that enhance sensitivity in low light levels do not appear to have consistent correlated effects on features related to visual acuity. A survey across several additional ant genera found numerous other pale species with enlarged eyes, suggesting these traits evolved multiple times within and across genera. We also compared the size of the anterior ocellus in workers of pale versus dark species of Myrmecocystus. In species with larger workers, the anterior ocellus was smaller in pale than in dark species, but this difference mostly disappeared for species with smaller workers. Presence of the anterior ocellus also was size-dependent in the two largest pale species.

Recent speciation and phenotypic plasticity within a parthenogenetic lineage of Levantine whip spiders (Chelicerata: Amblypygi: Charinidae)

Caves constitute ideal study systems for investigating adaptation and speciation, as the abiotic conditions shared by aphotic habitats exert a set of environmental filters on their communities. Arachnids constitute an important component of many cave ecosystems worldwide. We investigated the population genomics of two whip spider species: Sarax ioanniticus, a widely distributed parthenogenetic species found across the eastern Mediterranean; and S. israelensis, a recently described troglomorphic species that is endemic to caves in Israel. Here, we show that S. israelensis is completely genetically distinct from S. ioanniticus and most likely also constitutes a parthenogen. Counterintuitively, despite the lack of genetic variability within S. ioanniticus and S. israelensis, we discovered considerable variation in the degree of median eye reduction, particularly in the latter species. Natural history data from captive-bred specimens of S. israelensis validated the interpretation of parthenogenesis. Our results are most consistent with a scenario of a sexual ancestral species that underwent speciation, followed by independent transitions to apomictic parthenogenesis in each of the two daughter species. Moreover, the lack of genetic variability suggests that variation in eye morphology in S. israelensis is driven exclusively by epigenetic mechanisms.

Building on 150 Years of Knowledge: The Freshwater Isopod Asellus aquaticus as an Integrative Eco-Evolutionary Model System

Interactions between organisms and their environments are central to how biological diversity arises and how natural populations and ecosystems respond to environmental change. These interactions involve processes by which phenotypes are affected by or respond to external conditions (e.g., via phenotypic plasticity or natural selection) as well as processes by which organisms reciprocally interact with the environment (e.g., via eco-evolutionary feedbacks). Organism-environment interactions can be highly dynamic and operate on different hierarchical levels, from genes and phenotypes to populations, communities, and ecosystems. Therefore, the study of organism-environment interactions requires integrative approaches and model systems that are suitable for studies across different hierarchical levels. Here, we introduce the freshwater isopod Asellus aquaticus, a keystone species and an emerging invertebrate model system, as a prime candidate to address fundamental questions in ecology and evolution, and the interfaces therein. We review relevant fields of research that have used A. aquaticus and draft a set of specific scientific questions that can be answered using this species. Specifically, we propose that studies on A. aquaticus can help understanding (i) the influence of host-microbiome interactions on organismal and ecosystem function, (ii) the relevance of biotic interactions in ecosystem processes, and (iii) how ecological conditions and evolutionary forces facilitate phenotypic diversification.

Parallel morphological evolution and habitat-dependent sexual dimorphism in cave- vs. surface populations of the Asellus aquaticus (Crustacea: Isopoda: Asellidae) species complex

Studying parallel evolution (repeated, independent evolution of similar phenotypes in similar environments) is a powerful tool to understand environment-dependent selective forces. Surface-dwelling species that repeatedly and independently colonized caves provide unique models for such studies. The primarily surface-dwelling Asellus aquaticus species complex is a good candidate to carry out such research, because it colonized several caves in Europe. By comparing 17 functional morphological traits between six cave and nine surface populations of the A. aquaticus species complex, we investigated population divergence in morphology and sexual dimorphism. We found habitat-dependent population divergence in 10 out of 17 traits, likely reflecting habitat-driven changes in selection acting on sensory systems, feeding, grooming, and antipredator mechanisms. Sexual dimorphism was present in 15 traits, explained by sexual selection acting on male traits important in male-male agonistic behavior or mate guarding and fecundity selection acting on female traits affecting offspring number and nursing. In eight traits, the degree of sexual dimorphism was habitat dependent. We conclude that cave-related morphological changes are highly trait- and function-specific and that the strength of sexual/fecundity selection strongly differs between cave and surface habitats. The considerable population variation within habitat type warrants further studies to reveal cave-specific adaptations besides the parallel patterns.

Population divergence in aggregation and sheltering behaviour in surface- versus cave-adapted Asellus aquaticus (Crustacea: Isopoda)

Aggregation (gathering together) and sheltering (hiding in cover) are basic behaviours that can reduce the risk of predation. However, both behaviours have costs, such as increased competition over resources and high prevalence of contact-spread parasites (aggregation) or lost opportunities for foraging and mating (sheltering). Therefore, variation in these behaviours is expected between populations with varying levels of predation risk. We compared aggregation and sheltering in surface- (various predators) and cave-adapted (no predator) populations of the isopod Asellus aquaticus in a common garden experiment. Given that the cave environment is constantly dark, we also tested for population variation in light-induced behavioural plasticity. Variation in sheltering was explained by habitat type: cave individuals sheltered less than surface individuals. We found high between-population variation in aggregation with or without shelters and their light-induced plasticity, which was not explained by habitat type. Cave individuals decreased (habituation) whereas surface individuals increased sheltering with time (sensitization). We suggest that population variation in sheltering is driven by predation, whereas variation in aggregation must be driven by other, unaccounted environmental factors, in a similar manner to light-induced behavioural plasticity. Based on habituation/sensitization patterns, we suggest that predation-adapted populations are more sensitive to disturbance related to routine laboratory procedures.

Reinventing the wheel? Reassessing the roles of gene flow, sorting and convergence in repeated evolution

Biologists have long been intrigued by apparently predictable and repetitive evolutionary trajectories inferred across a variety of lineages and systems. In recent years, high-throughput sequencing analyses have started to transform our understanding of such repetitive shifts. While researchers have traditionally categorized such shifts as either "convergent" or "parallel," based on relatedness of the lineages involved, emerging genomic insights provide an opportunity to better describe the actual evolutionary mechanisms at play. A synthesis of recent genomic analyses confirms that convergence is the predominant driver of repetitive evolution among species, whereas repeated sorting of standing variation is the major driver of repeated shifts within species. However, emerging data reveal numerous notable exceptions to these expectations, with recent examples of de novo mutations underpinning convergent shifts among even very closely related lineages, while repetitive sorting processes have occurred among even deeply divergent taxa, sometimes via introgression. A number of very recent analyses have found evidence for both processes occurring on different scales within taxa. We suggest that the relative importance of convergent versus sorting processes depends on the interplay between gene flow among populations, and phylogenetic relatedness of the lineages involved.

The genomics of phenotypically differentiated Asellus aquaticus cave, surface stream and lake ecotypes

Organisms well suited for the study of ecotype formation have wide distribution ranges, where they adapt to multiple drastically different habitats repeatedly over space and time. Here we study such ecotypes in a Crustacean model, Asellus aquaticus, a commonly occurring isopod found in freshwater habitats as diverse as streams, caves and lakes. Previous studies focusing on cave vs. surface ecotypes have attributed depigmentation, eye loss and prolonged antennae to several south European cave systems. Likewise, surveys across multiple Swedish lakes have identified the presence of dark-pigmented "reed" and light-pigmented "stonewort" ecotypes, which can be found within the same lake. In this study, we sequenced the first draft genome of A. aquaticus, and subsequently use this to map reads and call variants in surface stream, cave and two lake ecotypes. In addition, the draft genome was combined with a RADseq approach to perform a quantitative trait locus (QTL) mapping study using a laboratory bred F₂ and F₄ cave × surface intercross. We identified genomic regions associated with body pigmentation, antennae length and body size. Furthermore, we compared genome-wide differentiation between natural populations and found several genes potentially associated with these habitats. The assessment of the cave QTL regions in the light-dark comparison of lake populations suggests that the regions associated with cave adaptation are also involved with genomic differentiation in the lake ecotypes. These demonstrate how troglomorphic adaptations can be used as a model for related ecotype formation.

Systemic paralogy and function of retinal determination network homologs in arachnids

BACKGROUND

Arachnids are important components of cave ecosystems and display many examples of troglomorphisms, such as blindness, depigmentation, and elongate appendages. Little is known about how the eyes of arachnids are specified genetically, let alone the mechanisms for eye reduction and loss in troglomorphic arachnids. Additionally, duplication of Retinal Determination Gene Network (RDGN) homologs in spiders has convoluted functional inferences extrapolated from single-copy homologs in pancrustacean models.

RESULTS

We investigated a sister species pair of Israeli cave whip spiders, Charinus ioanniticus and C. israelensis (Arachnopulmonata, Amblypygi), of which one species has reduced eyes. We generated embryonic transcriptomes for both Amblypygi species, and discovered that several RDGN homologs exhibit duplications. We show that duplication of RDGN homologs is systemic across arachnopulmonates (arachnid orders that bear book lungs), rather than being a spider-specific phenomenon. A differential gene expression (DGE) analysis comparing the expression of RDGN genes in field-collected embryos of both species identified candidate RDGN genes involved in the formation and reduction of eyes in whip spiders. To ground bioinformatic inference of expression patterns with functional experiments, we interrogated the function of three candidate RDGN genes identified from DGE using RNAi in the spider Parasteatoda tepidariorum. We provide functional evidence that one of these paralogs, sine oculis/Six1 A (soA), is necessary for the development of all arachnid eye types.

CONCLUSIONS

Our work establishes a foundation to investigate the genetics of troglomorphic adaptations in cave arachnids, and links differential gene expression to an arthropod eye phenotype for the first time outside of Pancrustacea. Our results support the conservation of at least one RDGN component across Arthropoda and provide a framework for identifying the role of gene duplications in generating arachnid eye diversity.

Hard to get, easy to lose: Evolution of mantle photoreceptor organs in bivalves (Bivalvia, Pteriomorphia)

Morphologically diverse eyes have evolved numerous times, yet little is known about how eye gain and loss is related to photic environment. The pteriomorphian bivalves (e.g., oysters, scallops, and ark clams), with a remarkable range of photoreceptor organs and ecologies, are a suitable system to investigate the association between eye evolution and ecological shifts. The present phylogenetic framework was based on amino acid sequences from transcriptome datasets and nucleotide sequences of five additional genes. In total, 197 species comprising 22 families from all five pteriomorphian orders were examined, representing the greatest taxonomic sampling to date. Morphological data were acquired for 162 species and lifestyles were compiled from the literature for all 197 species. Photoreceptor organs occur in 11 families and have arisen exclusively in epifaunal lineages, that is, living above the substrate, at least five times independently. Models for trait evolution consistently recovered higher rates of loss over gain. Transitions to crevice-dwelling habit appear associated with convergent gains of eyespots in epifaunal lineages. Once photoreceptor organs have arisen, multiple losses occurred in lineages that shift to burrowing lifestyles and deep-sea habitats. The observed patterns suggest that eye evolution in pteriomorphians might have evolved in association with light-guided behaviors, such as phototaxis, body posture, and alarm responses.

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.

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.