Contrasting Gene Decay in Subterranean Vertebrates: Insights from Cavefishes and Fossorial Mammals

[Massive gene losses: the case of chemoreceptors in aquatic tetrapods]

The existence of multigene families is the result of gene duplication. In vertebrates, the genes coding for the chemoreceptors involved in olfaction often form large families, sometimes comprising several thousand genes. In this case, it is assumed that a large number of genes is essential to discriminate between a large number of odorant molecules, and that the ability to identify a large number of odors may be necessary for different purposes, such as finding food and sexual partners, and avoiding predators. However, it has been found that in some species the number of these genes is very small, resulting from the secondary loss of many genes. This massive loss of genes is not always clearly associated with a change in the biology of the species or its environment, but in some cases, it is associated with drastic changes, such as the return to aquatic life in tetrapods.

Evolution of the regulation of developmental gene expression in blind Mexican cavefish

Developmental evolution and diversification of morphology can arise through changes in the regulation of gene expression or protein-coding sequence. To unravel mechanisms underlying early developmental evolution in cavefish of the species Astyanax mexicanus, we compared transcriptomes of surface-dwelling and blind cave-adapted morphs at the end of gastrulation. Twenty percent of the transcriptome was differentially expressed. Allelic expression ratios in cave X surface hybrids showed that cis-regulatory changes are the quasi-exclusive contributors to inter-morph variations in gene expression. Among a list of 108 genes with change at the cis-regulatory level, we explored the control of expression of rx3, which is a master eye gene. We discovered that cellular rx3 levels are cis-regulated in a cell-autonomous manner, whereas rx3 domain size depends on non-autonomous Wnt and Bmp signalling. These results highlight how uncoupled mechanisms and regulatory modules control developmental gene expression and shape morphological changes. Finally, a transcriptome-wide search for fixed coding mutations and differential exon use suggested that variations in coding sequence have a minor contribution. Thus, during early embryogenesis, changes in gene expression regulation are the main drivers of cavefish developmental evolution.

Genome sequencing of Coryphaenoides yaquinae reveals convergent and lineage-specific molecular evolution in deep-sea adaptation

Abyssal (3501-6500 m) and hadal (>6500 m) fauna evolve under harsh abiotic stresses, characterized by high hydrostatic pressure, darkness and food shortage, providing unique opportunities to investigate mechanisms underlying environmental adaptation. Genomes of several hadal species have recently been reported. However, the genetic adaptation of deep sea species across a broad spectrum of ocean depths has yet to be thoroughly investigated, due to the challenges imposed by collecting the deep sea species. To elucidate the correlation between genetic innovation and vertical distribution, we generated a chromosome-level genome assembly of the macrourids Coryphaenoides yaquinae, which is widely distributed in the abyssal/hadal zone ranging from 3655 to 7259 m in depth. Genomic comparisons among shallow, abyssal and hadal-living species identified idiosyncratic and convergent genetic alterations underlying the extraordinary adaptations of deep-sea species including light perception, circadian regulation, hydrostatic pressure and hunger tolerance. The deep-sea fishes (Coryphaenoides Sp. and Pseudoliparis swirei) venturing into various ocean depths independently have undergone convergent amino acid substitutions in multiple proteins such as rhodopsin 1, pancreatic and duodenal homeobox 1 and melanocortin 4 receptor which are known or verified in zebrafish to be related with vision adaptation and energy expenditure. Convergent evolution events were also identified in heat shock protein 90 beta family member 1 and valosin-containing protein genes known to be related to hydrostatic pressure adaptation specifically in fishes found around the hadal range. The uncovering of the molecular convergence among the deep-sea species shed new light on the common genetic innovations required for deep-sea adaptation by the fishes.

Sporadic feeding regulates robust food entrainable circadian clocks in blind cavefish

The circadian clock represents a key timing system entrained by various periodic signals that ensure synchronization with the environment. Many investigations have pointed to the existence of two distinct circadian oscillators: one regulated by the light-dark cycle and the other set by feeding time. Blind cavefish have evolved under extreme conditions where they completely lack light exposure and experience food deprivation. Here, we have investigated feeding regulated clocks in two cavefish species, the Somalian cavefish Phreatichthys andruzzii and the Mexican cavefish Astyanax mexicanus, in comparison with the surface-dwelling zebrafish Danio rerio. Our results reveal that feeding represents an extremely strong synchronizer for circadian locomotor rhythmicity in subterranean cavefish. Indeed, we showed that consuming just one meal every 4 days is sufficient to entrain circadian rhythmicity in both cavefish species, but not in zebrafish. These profound adaptations to an extreme environment provide insight into the connections between feeding and circadian clocks.

The 'in's and out's' of descending pain modulation from the rostral ventromedial medulla

The descending-pain modulating circuit controls the experience of pain by modulating transmission of sensory signals through the dorsal horn. This circuit's key output node, the rostral ventromedial medulla (RVM), integrates 'top-down' and 'bottom-up' inputs that regulate functionally defined RVM cell types, 'OFF-cells' and 'ON-cells', which respectively suppress or facilitate pain-related sensory processing. While recent advances have sought molecular definition of RVM cell types, conflicting behavioral findings highlight challenges involved in aligning functional and molecularly defined types. This review summarizes current understanding, derived primarily from rodent studies but with corroborating evidence from human imaging, of the role of RVM populations in pain modulation and persistent pain states and explores recent advances outlining inputs to, and outputs from, RVM pain-modulating neurons.

Gene loss and relaxed selection of plaat1 in vertebrates adapted to low-light environments

Gene loss is an important mechanism for evolution in low-light or cave environments where visual adaptations often involve a reduction or loss of eyesight. The plaat gene family encodes phospholipases essential for the degradation of organelles in the lens of the eye. These phospholipases translocate to damaged organelle membranes, inducing them to rupture. This rupture is required for lens transparency and is essential for developing a functioning eye. Plaat3 is thought to be responsible for this role in mammals, while plaat1 is thought to be responsible in other vertebrates. We used a macroevolutionary approach and comparative genomics to examine the origin, loss, synteny and selection of plaat1 across bony fishes and tetrapods. We showed that plaat1 (probably ancestral to all bony fish + tetrapods) has been lost in squamates and is significantly degraded in lineages of low-visual-acuity and blind mammals and fishes. Our findings suggest that plaat1 is important for visual acuity across bony vertebrates, and that its loss through relaxed selection and pseudogenization may have played a role in the repeated evolution of visual systems in low-light environments. Our study sheds light on the importance of gene-loss in trait evolution and provides insights into the mechanisms underlying visual acuity in low-light environments.

A reference genome for the Andean cavefish Trichomycterus rosablanca (Siluriformes, Trichomycteridae): Building genomic resources to study evolution in cave environments

Animals living in caves are of broad relevance to evolutionary biologists interested in understanding the mechanisms underpinning convergent evolution. In the Eastern Andes of Colombia, populations from at least two distinct clades of Trichomycterus catfishes (Siluriformes) independently colonized cave environments and converged in phenotype by losing their eyes and pigmentation. We are pursuing several research questions using genomics to understand the evolutionary forces and molecular mechanisms responsible for repeated morphological changes in this system. As a foundation for such studies, here we describe a diploid, chromosome-scale, long-read reference genome for Trichomycterus rosablanca, a blind, depigmented species endemic to the karstic system of the department of Santander. The nuclear genome comprises 1 Gb in 27 chromosomes, with a 40.0× HiFi long-read genome coverage having an N50 scaffold of 40.4 Mb and N50 contig of 13.1 Mb, with 96.9% (Eukaryota) and 95.4% (Actinopterygii) universal single-copy orthologs (BUSCO). This assembly provides the first reference genome for the speciose genus Trichomycterus, serving as a key resource for research on the genomics of phenotypic evolution.

The nature and distribution of putative non-functional alleles suggest only two independent events at the origins of Astyanax mexicanus cavefish populations

BACKGROUND

Several studies suggested that cavefish populations of Astyanax mexicanus settled during the Late Pleistocene. This implies that the cavefish's most conspicuous phenotypic changes, blindness and depigmentation, and more cryptic characters important for cave life, evolved rapidly.

RESULTS

Using the published genomes of 47 Astyanax cavefish from la Cueva de El Pachón, El Sótano de la Tinaja, La Cueva Chica and El Sótano de Molino, we searched for putative loss-of-function mutations in previously defined sets of genes, i.e., vision, circadian clock and pigmentation genes. Putative non-functional alleles for four vision genes were identified. Then, we searched genome-wide for putative non-functional alleles in these four cave populations. Among 512 genes with segregating putative non-functional alleles in cavefish that are absent in surface fish, we found an enrichment in visual perception genes. Among cavefish populations, different levels of shared putative non-functional alleles were found. Using a subset of 12 genes for which putative loss-of-function mutations were found, we extend the analysis of shared pseudogenes to 11 cave populations. Using a subset of six genes for which putative loss-of-function mutations were found in the El Sótano del Toro population, where extensive hybridization with surface fish occurs, we found a correlation between the level of eye regression and the amount of putative non-functional alleles.

CONCLUSIONS

We confirm that very few putative non-functional alleles are present in a large set of vision genes, in accordance with the recent origin of Astyanax mexicanus cavefish. Furthermore, the genome-wide analysis indicates an enrichment of putative loss-of-function alleles in genes with vision-related GO-terms, suggesting that visual perception may be the function chiefly impacted by gene losses related to the shift from a surface to a cave environment. The geographic distribution of putative loss-of-function alleles newly suggests that cave populations from Sierra de Guatemala and Sierra de El Abra share a common origin, albeit followed by independent evolution for a long period. It also supports that populations from the Micos area have an independent origin. In El Sótano del Toro, the troglomorphic phenotype is maintained despite massive introgression of the surface genome.

The genome of the blind bee louse fly reveals deep convergences with its social host and illuminates Drosophila origins

Social insects' nests harbor intruders known as inquilines,1 which are usually related to their hosts.2,3 However, distant non-social inquilines may also show convergences with their hosts,4,5 although the underlying genomic changes remain unclear. We analyzed the genome of the wingless and blind bee louse fly Braula coeca, an inquiline kleptoparasite of the western honey bee, Apis mellifera.6,7 Using large phylogenomic data, we confirmed recent accounts that the bee louse fly is a drosophilid8,9 and showed that it had likely evolved from a sap-breeder ancestor associated with honeydew and scale insects' wax. Unlike many parasites, the bee louse fly genome did not show significant erosion or strict reliance on an endosymbiont, likely due to a relatively recent age of inquilinism. However, we observed a horizontal transfer of a transposon and a striking parallel evolution in a set of gene families between the honey bee and the bee louse fly. Convergences included genes potentially involved in metabolism and immunity and the loss of nearly all bitter-tasting gustatory receptors, in agreement with life in a protective nest and a diet of honey, pollen, and beeswax. Vision and odorant receptor genes also exhibited rapid losses. Only genes whose orthologs in the closely related Drosophila melanogaster respond to honey bee pheromone components or floral aroma were retained, whereas the losses included orthologous receptors responsive to the anti-ovarian honey bee queen pheromones. Hence, deep genomic convergences can underlie major phenotypic transitions during the evolution of inquilinism between non-social parasites and their social hosts.

Convergent genomic signatures associated with vertebrate viviparity

BACKGROUND

Viviparity-live birth-is a complex and innovative mode of reproduction that has evolved repeatedly across the vertebrate Tree of Life. Viviparous species exhibit remarkable levels of reproductive diversity, both in the amount of care provided by the parent during gestation, and the ways in which that care is delivered. The genetic basis of viviparity has garnered increasing interest over recent years; however, such studies are often undertaken on small evolutionary timelines, and thus are not able to address changes occurring on a broader scale. Using whole genome data, we investigated the molecular basis of this innovation across the diversity of vertebrates to answer a long held question in evolutionary biology: is the evolution of convergent traits driven by convergent genomic changes?

RESULTS

We reveal convergent changes in protein family sizes, protein-coding regions, introns, and untranslated regions (UTRs) in a number of distantly related viviparous lineages. Specifically, we identify 15 protein families showing evidence of contraction or expansion associated with viviparity. We additionally identify elevated substitution rates in both coding and noncoding sequences in several viviparous lineages. However, we did not find any convergent changes-be it at the nucleotide or protein level-common to all viviparous lineages.

CONCLUSIONS

Our results highlight the value of macroevolutionary comparative genomics in determining the genomic basis of complex evolutionary transitions. While we identify a number of convergent genomic changes that may be associated with the evolution of viviparity in vertebrates, there does not appear to be a convergent molecular signature shared by all viviparous vertebrates. Ultimately, our findings indicate that a complex trait such as viviparity likely evolves with changes occurring in multiple different pathways.

A collaborative backbone resource for comparative studies of subterranean evolution: The World Asellidae database

Transition to novel environments, such as groundwater colonization by surface organisms, provides an excellent research ground to study phenotypic evolution. However, interspecific comparative studies on evolution to groundwater life are few because of the challenge in assembling large ecological and molecular resources for species-rich taxa comprised of surface and subterranean species. Here, we make available to the scientific community an operational set of working tools and resources for the Asellidae, a family of freshwater isopods containing hundreds of surface and subterranean species. First, we release the World Asellidae database (WAD) and its web application, a sustainable and FAIR solution to producing and sharing data and biological material. WAD provides access to thousands of species occurrences, specimens, DNA extracts and DNA sequences with rich metadata ensuring full scientific traceability. Second, we perform a large-scale dated phylogenetic reconstruction of Asellidae to support phylogenetic comparative analyses. Of 424 terminal branches, we identify 34 pairs of surface and subterranean species representing independent replicates of the transition from surface water to groundwater. Third, we exemplify the usefulness of WAD for documenting phenotypic shifts associated with colonization of subterranean habitats. We provide the first phylogenetically controlled evidence that body size of males decreases relative to that of females upon groundwater colonization, suggesting competition for rare receptive females selects for smaller, more agile males in groundwater. By making these tools and resources widely accessible, we open up new opportunities for exploring how phenotypic traits evolve in response to changes in selective pressures and trade-offs during groundwater colonization.

Chromosome-level genome assembly of hadal snailfish reveals mechanisms of deep-sea adaptation in vertebrates

As the deepest vertebrate in the ocean, the hadal snailfish (Pseudoliparis swirei), which lives at a depth of 6,000-8,000 m, is a representative case for studying adaptation to extreme environments. Despite some preliminary studies on this species in recent years, including their loss of pigmentation, visual and skeletal calcification genes, and the role of trimethylamine N-oxide in adaptation to high-hydrostatic pressure, it is still unknown how they evolved and why they are among the few vertebrate species that have successfully adapted to the deep-sea environment. Using genomic data from different trenches, we found that the hadal snailfish may have entered and fully adapted to such extreme environments only in the last few million years. Meanwhile, phylogenetic relationships show that they spread into different trenches in the Pacific Ocean within a million years. Comparative genomic analysis has also revealed that the genes associated with perception, circadian rhythms, and metabolism have been extensively modified in the hadal snailfish to adapt to its unique environment. More importantly, the tandem duplication of a gene encoding ferritin significantly increased their tolerance to reactive oxygen species, which may be one of the important factors in their adaptation to high-hydrostatic pressure.

Gene loss and relaxed selection of plaat1 in vertebrates adapted to low-light environments

Gene loss is an important mechanism for evolution in low-light or cave environments where visual adaptations often involve a reduction or loss of eyesight. The plaat gene family are phospholipases essential for the degradation of organelles in the lens of the eye. They translocate to damaged organelle membranes, inducing them to rupture. This rupture is required for lens transparency and is essential for developing a functioning eye. Plaat3 is thought to be responsible for this role in mammals, while plaat1 is thought to be responsible in other vertebrates. We used a macroevolutionary approach and comparative genomics to examine the origin, loss, synteny, and selection of plaat1 across bony fishes and tetrapods. We show that plaat1 (likely ancestral to all bony fish + tetrapods) has been lost in squamates and is significantly degraded in lineages of low-visual acuity and blind mammals and fish. Our findings suggest that plaat1 is important for visual acuity across bony vertebrates, and that its loss through relaxed selection and pseudogenization may have played a role in the repeated evolution of visual systems in low-light-environments. Our study sheds light on the importance of gene-loss in trait evolution and provides insights into the mechanisms underlying visual acuity in low-light environments.

An Astyanax mexicanus mao knockout line uncovers the developmental roles of monoamine homeostasis in fish brain

Monoaminergic systems are conserved in vertebrates, yet they present variations in neuroanatomy, genetic components and functions across species. MonoAmine Oxidase, or MAO, is the enzyme responsible for monoamine degradation. While mammals possess two genes, MAO-A and MAO-B, fish possess one single mao gene. To study the function of MAO and monoamine homeostasis on fish brain development and physiology, here we have generated a mao knockout line in Astyanax mexicanus (surface fish), by CRISPR/Cas9 technology. Homozygote mao knockout larvae died at 13 days post-fertilization. Through a time-course analysis, we report that hypothalamic serotonergic neurons undergo fine and dynamic regulation of serotonin level upon loss of mao function, in contrast to those in the raphe, which showed continuously increased serotonin levels - as expected. Dopaminergic neurons were not affected by mao loss-of-function. At behavioral level, knockout fry showed a transient decrease in locomotion that followed the variations in the hypothalamus serotonin neuronal levels. Finally, we discovered a drastic effect of mao knockout on brain progenitors proliferation in the telencephalon and hypothalamus, including a reduction in the number of proliferative cells and an increase of the cell cycle length. Altogether, our results show that MAO has multiple and varied effects on Astyanax mexicanus brain development. Mostly, they bring novel support to the idea that serotonergic neurons in the hypothalamus and raphe of the fish brain are different in nature and identity, and they unravel a link between monoaminergic homeostasis and brain growth.

AlexandrusPS: A User-Friendly Pipeline for the Automated Detection of Orthologous Gene Clusters and Subsequent Positive Selection Analysis

The detection of adaptive selection in a system approach considering all protein-coding genes allows for the identification of mechanisms and pathways that enabled adaptation to different environments. Currently, available programs for the estimation of positive selection signals can be divided into two groups. They are either easy to apply but can analyze only one gene family at a time, restricting system analysis; or they can handle larger cohorts of gene families, but require considerable prerequisite data such as orthology associations, codon alignments, phylogenetic trees, and proper configuration files. All these steps require extensive computational expertise, restricting this endeavor to specialists. Here, we introduce AlexandrusPS, a high-throughput pipeline that overcomes technical challenges when conducting transcriptome-wide positive selection analyses on large sets of nucleotide and protein sequences. The pipeline streamlines 1) the execution of an accurate orthology prediction as a precondition for positive selection analysis, 2) preparing and organizing configuration files for CodeML, 3) performing positive selection analysis using CodeML, and 4) generating an output that is easy to interpret, including all maximum likelihood and log-likelihood test results. The only input needed from the user is the CDS and peptide FASTA files of proteins of interest. The pipeline is provided in a Docker image, requiring no program or module installation, enabling the application of the pipeline in any computing environment. AlexandrusPS and its documentation are available via GitHub (https://github.com/alejocn5/AlexandrusPS).

Acoustic signatures in Mexican cavefish populations inhabiting different caves

Complex patterns of acoustic communication exist throughout the animal kingdom, including underwater. The river-dwelling and the Pachón cave-adapted morphotypes of the fish Astyanax mexicanus are soniferous and share a repertoire of sounds. Their function and significance is mostly unknown. Here, we explored whether and how sounds produced by blind cavefishes inhabiting different Mexican caves may vary. We compared "Clicks" and "Serial Clicks" produced by cavefish in six different caves distributed in three mountain ranges in Mexico. We also sampled laboratory-bred cavefish lines originating from four of these caves. Sounds were extracted and analyzed using both a manual method and a machine learning-based automation tool developed in-house. Multi-parametric analyses suggest wild cave-specific acoustic signatures, or "accents". An acoustic code also existed in laboratory cavefish lines, suggesting a genetic basis for the evolution of this trait. The variations in acoustic parameters between caves of origin did not seem related to fish phenotypes, phylogeography or ecological conditions. We propose that the evolution of such acoustic signatures would progressively lead to the differentiation of local accents that may prevent interbreeding and thus contribute to speciation.

Genetic identification and reiterated captures suggest that the Astyanax mexicanus El Pachón cavefish population is closed and declining

The sizes of Astyanax mexicanus blind cavefish populations of North-East Mexico are demographic parameters of great importance for investigating a variety of ecological, evolutionary, and conservation issues. However, few estimates have been obtained. For these mobile animals living in an environment difficult to explore as a whole, methods based on capture-mark-recapture are appropriate, but their feasibility and interpretation of results depend on several assumptions that must be carefully examined. Here, we provide evidence that minimally invasive genetic identification from captures at different time intervals (three days and three years) can give insights into cavefish population size dynamics as well as other important demographic parameters of interest. We also provide tools to calibrate sampling and genotyping efforts necessary to reach a given level of precision. Our results suggest that the El Pachón cave population is currently very small, of an order of magnitude of a few hundreds of individuals, and is distributed in a relatively isolated area. The probable decline in population size in the El Pachón cave since the last census in 1971 raises serious conservation issues.

Mitochondrial phylogeography and molecular evolution of the rhodopsin visual pigment in troglobitic populations of Astyanax mexicanus (De Filippi, 1853)

Cave-adapted animals provide a unique opportunity to study the evolutionary mechanisms underlying phenotypic, metabolic, behavioral, and genetic evolution in response to cave environments. The Mexican tetra ( Astyanax mexicanus) is considered a unique model system as it shows both surface and cave-dwelling morphs. To date, at least 33 different cave populations have been identified, with phylogenetic studies suggesting an origin from at least two independent surface lineages, thereby providing a unique opportunity to study parallel evolution. In the present study, we carried out the most exhaustive phylogeographic study of A. mexicanus to date, including cave and surface localities, using two mitochondrial markers (cytochrome b (cyt b) and cytochrome c oxidase subunit I ( COI)) and nuclear rhodopsin visual pigment ( rho). Additionally, we inferred the molecular evolution of rho within the two contrasting environments (cave and surface) and across three geographic regions (Sierra de El Abra, Sierra de Guatemala, and Micos). In total, 267 individuals were sequenced for the two mitochondrial fragments and 268 individuals were sequenced for the rho visual pigment from 22 cave and 46 surface populations. Phylogeographic results based on the mitochondrial data supported the two-lineage hypothesis, except for the Pachón and Chica caves, whose introgression has been largely documented. The Sierra de El Abra region depicted the largest genetic diversity, followed by the Sierra de Guatemala region. Regarding the phylogeographic patterns of rho, we recovered exclusive haplogroups for the Sierra de El Abra (Haplogroup I) and Sierra de Guatemala regions (Haplogroup IV). Moreover, a 544 bp deletion in the rho gene was observed in the Escondido cave population from Sierra de Guatemala, reducing the protein from seven to three intramembrane domains. This change may produce a loss-of-function (LOF) but requires further investigation. Regarding nonsynonymous ( dN) and synonymous ( dS) substitution rates (omega values ω), our results revealed the prevailing influence of purifying selection upon the rho pigment for both cave and surface populations (ω<1), but relaxation at the El Abra region. Notably, in contrast to the other two regions, we observed an increase in the number of dN mutations for Sierra de El Abra. However, given that a LOF was exclusively identified in the Sierra de Guatemala region, we cannot dismiss the possibility of a pleiotropic effect on the Rho protein.

Parallel and convergent genomic changes underlie independent subterranean colonization across beetles

Adaptation to life in caves is often accompanied by dramatically convergent changes across distantly related taxa, epitomized by the loss or reduction of eyes and pigmentation. Nevertheless, the genomic underpinnings underlying cave-related phenotypes are largely unexplored from a macroevolutionary perspective. Here we investigate genome-wide gene evolutionary dynamics in three distantly related beetle tribes with at least six instances of independent colonization of subterranean habitats, inhabiting both aquatic and terrestrial underground systems. Our results indicate that remarkable gene repertoire changes mainly driven by gene family expansions occurred prior to underground colonization in the three tribes, suggesting that genomic exaptation may have facilitated a strict subterranean lifestyle parallelly across beetle lineages. The three tribes experienced both parallel and convergent changes in the evolutionary dynamics of their gene repertoires. These findings pave the way towards a deeper understanding of the evolution of the genomic toolkit in hypogean fauna.

Photoreceptor genes in a trechine beetle, Trechiama kuznetsovi, living in the upper hypogean zone

To address how organisms adapt to a new environment, subterranean organisms whose ancestors colonized subterranean habitats from surface habitats have been studied. Photoreception abilities have been shown to have degenerated in organisms living in caves and calcrete aquifers. Meanwhile, the organisms living in a shallow subterranean environment, which are inferred to reflect an intermediate stage in an evolutionary pathway to colonization of a deeper subterranean environment, have not been studied well. In the present study, we examined the photoreception ability in a trechine beetle, Trechiama kuznetsovi, which inhabits the upper hypogean zone and has a vestigial compound eye. By de novo assembly of genome and transcript sequences, we were able to identify photoreceptor genes and phototransduction genes. Specifically, we focused on opsin genes, where one long wavelength opsin gene and one ultraviolet opsin gene were identified. The encoded amino acid sequences had neither a premature stop codon nor a frameshift mutation, and appeared to be subject to purifying selection. Subsequently, we examined the internal structure of the compound eye and nerve tissue in the adult head, and found potential photoreceptor cells in the compound eye and nerve bundle connected to the brain. The present findings suggest that T. kuznetsovi has retained the ability of photoreception. This species represents a transitional stage of vision, in which the compound eye regresses, but it may retain the ability of photoreception using the vestigial eye.

Reduction of Paraoxonase Expression Followed by Inactivation across Independent Semiaquatic Mammals Suggests Stepwise Path to Pseudogenization

Convergent adaptation to the same environment by multiple lineages frequently involves rapid evolutionary change at the same genes, implicating these genes as important for environmental adaptation. Such adaptive molecular changes may yield either change or loss of protein function; loss of function can eliminate newly deleterious proteins or reduce energy necessary for protein production. We previously found a striking case of recurrent pseudogenization of the Paraoxonase 1 (Pon1) gene among aquatic mammal lineages-Pon1 became a pseudogene with genetic lesions, such as stop codons and frameshifts, at least four times independently in aquatic and semiaquatic mammals. Here, we assess the landscape and pace of pseudogenization by studying Pon1 sequences, expression levels, and enzymatic activity across four aquatic and semiaquatic mammal lineages: pinnipeds, cetaceans, otters, and beavers. We observe in beavers and pinnipeds an unexpected reduction in expression of Pon3, a paralog with similar expression patterns but different substrate preferences. Ultimately, in all lineages with aquatic/semiaquatic members, we find that preceding any coding-level pseudogenization events in Pon1, there is a drastic decrease in expression, followed by relaxed selection, thus allowing accumulation of disrupting mutations. The recurrent loss of Pon1 function in aquatic/semiaquatic lineages is consistent with a benefit to Pon1 functional loss in aquatic environments. Accordingly, we examine diving and dietary traits across pinniped species as potential driving forces of Pon1 functional loss. We find that loss is best associated with diving activity and likely results from changes in selective pressures associated with hypoxia and hypoxia-induced inflammation.

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.

Molecular evolution of vision-related genes may contribute to marsupial photic niche adaptations

Vision plays an essential role in the life of many animals. While most mammals are night-active (nocturnal), many have adapted to novel light environments. This includes diurnal (day-active) and crepuscular (twilight-active) species. Here, we used integrative approaches to investigate the molecular evolution of 112 vision-related genes across 19 genomes representing most marsupial orders. We found that four genes (GUCA1B, GUCY2F, RGR, and SWS2) involved in retinal phototransduction likely became functionally redundant in the ancestor of marsupials, a group of largely obligate nocturnal mammals. We also show evidence of rapid evolution and positive selection of bright-light vision genes in the common ancestor of Macropus (kangaroos, wallaroos, and wallabies). Macropus-specific amino acid substitutions in opsin genes (LWS and SWS1), in particular, may be an adaptation for crepuscular vision in this genus via opsin spectral sensitivity tuning. Our study set the stage for functional genetics studies and provides a stepping stone to future research efforts that fully capture the visual repertoire of marsupials.

[The toolbox of developmental evolution or how Mexican cave fishes lost their eyes]

The fish Astyanax mexicanus comes in two very different forms: a "normal" river morph, and a blind, depigmented cave morph, living in the total and permanent darkness of Mexican caves. This species is on the way to becoming a model of choice in evolutionary and comparative biology, both for the study of the evolution of behavior, physiology or morphology, and for molecular genetics or population genetics. Here, I present the advancement of knowledge in the field of the developmental evolution of the eye of the cave morph. By rewinding back in time its development from the eye of the larva to the retinal field at the end of gastrulation, the cave-dwelling Astyanax embryo reveals mechanisms and processes likely to contribute to evolutionary variations between species, but also to pathological variations in the morphogenesis of the optic region.

Parallel decay of vision genes in subterranean water beetles

In the framework of neutral theory of molecular evolution, genes specific to the development and function of eyes in subterranean animals living in permanent darkness are expected to evolve by relaxed selection, ultimately becoming pseudogenes. However, definitive empirical evidence for the role of neutral processes in the loss of vision over evolutionary time remains controversial. In previous studies, we characterized an assemblage of independently-evolved water beetle (Dytiscidae) species from a subterranean archipelago in Western Australia, where parallel vision and eye loss have occurred. Using a combination of transcriptomics and exon capture, we present evidence of parallel coding sequence decay, resulting from the accumulation of frameshift mutations and premature stop codons, in eight phototransduction genes (arrestins, opsins, ninaC and transient receptor potential channel genes) in 32 subterranean species in contrast to surface species, where these genes have open reading frames. Our results provide strong evidence to support neutral evolutionary processes as a major contributing factor to the loss of phototransduction genes in subterranean animals, with the ultimate fate being the irreversible loss of a light detection system.

BASE: A novel workflow to integrate nonubiquitous genes in comparative genomics analyses for selection

Inferring the selective forces that orthologous genes underwent across different lineages can help us understand the evolutionary processes that have shaped their extant diversity and the phenotypes they underlie. The most widespread metric to estimate the selection regimes of coding genes-across sites and phylogenies-is the ratio of nonsynonymous to synonymous substitutions (dN/dS, also known as ω). Nowadays, modern sequencing technologies and the large amount of already available sequence data allow the retrieval of thousands of orthologous genes across large numbers of species. Nonetheless, the tools available to explore selection regimes are not designed to automatically process all genes, and their practical usage is often restricted to the single-copy ones which are found across all species considered (i.e., ubiquitous genes). This approach limits the scale of the analysis to a fraction of single-copy genes, which can be as low as an order of magnitude in respect to those which are not consistently found in all species considered (i.e., nonubiquitous genes). Here, we present a workflow named BASE that-leveraging the CodeML framework-eases the inference and interpretation of gene selection regimes in the context of comparative genomics. Although a number of bioinformatics tools have already been developed to facilitate this kind of analyses, BASE is the first to be specifically designed to allow the integration of nonubiquitous genes in a straightforward and reproducible manner. The workflow-along with all relevant documentation-is available at github.com/for-giobbe/BASE.

Chromosome-Level Assembly of the Atlantic Silverside Genome Reveals Extreme Levels of Sequence Diversity and Structural Genetic Variation

The levels and distribution of standing genetic variation in a genome can provide a wealth of insights about the adaptive potential, demographic history, and genome structure of a population or species. As structural variants are increasingly associated with traits important for adaptation and speciation, investigating both sequence and structural variation is essential for wholly tapping this potential. Using a combination of shotgun sequencing, 10x Genomics linked reads and proximity-ligation data (Chicago and Hi-C), we produced and annotated a chromosome-level genome assembly for the Atlantic silverside (Menidia menidia)-an established ecological model for studying the phenotypic effects of natural and artificial selection-and examined patterns of genomic variation across two individuals sampled from different populations with divergent local adaptations. Levels of diversity varied substantially across each chromosome, consistently being highly elevated near the ends (presumably near telomeric regions) and dipping to near zero around putative centromeres. Overall, our estimate of the genome-wide average heterozygosity in the Atlantic silverside is among the highest reported for a fish, or any vertebrate (1.32-1.76% depending on inference method and sample). Furthermore, we also found extreme levels of structural variation, affecting ∼23% of the total genome sequence, including multiple large inversions (> 1 Mb and up to 12.6 Mb) associated with previously identified haploblocks showing strong differentiation between locally adapted populations. These extreme levels of standing genetic variation are likely associated with large effective population sizes and may help explain the remarkable adaptive divergence among populations of the Atlantic silverside.

Evolutionary dynamics of the OR gene repertoire in teleost fishes: evidence of an association with changes in olfactory epithelium shape

Teleost fishes perceive their environment through a range of sensory modalities, among which olfaction often plays an important role. Richness of the olfactory repertoire depends on the diversity of receptors coded by homologous genes classified into four families: OR, TAAR, VR1, and VR2. Herein, we focus on the OR gene repertoire. While independent large contractions of the OR gene repertoire associated with ecological transitions have been found in mammals, little is known about the diversity of the OR gene repertoire and its evolution in teleost fishes, a group that includes more than 34,000 living species. We analyzed genomes of 163 species representing diversity in this large group. We found a large range of variation in the number of functional OR genes, from 15 in the Broad-nose Pipefish Syngnathus typhle and the Ocean Sunfish Mola mola, to 429 in the Zig-zag Eel Mastacembelus armatus. The number of OR genes was higher in species when a multilamellar olfactory rosette was present. Moreover, the number of lamellae was correlated with the richness of the OR gene repertoire. While a slow and balanced birth-and-death process generally drives the evolution of the OR gene repertoire, we inferred several episodes of high rates of gene loss, sometimes followed by large gains in the number of OR genes. These gains coincide with morphological changes of the olfactory organ and suggest a strong functional association between changes in the morphology and the evolution of the OR gene repertoire.

Diversity of Olfactory Responses and Skills in Astyanax Mexicanus Cavefish Populations Inhabiting different Caves

Animals in many phyla are adapted to and thrive in the constant darkness of subterranean environments. To do so, cave animals have presumably evolved mechano- and chemosensory compensations to the loss of vision, as is the case for the blind characiform cavefish, Astyanax mexicanus. Here, we systematically assessed the olfactory capacities of cavefish and surface fish of this species in the lab as well as in the wild, in five different caves in northeastern Mexico, using an olfactory setup specially developed to test and record olfactory responses during fieldwork. Overall cavefish showed lower (i.e., better) olfactory detection thresholds than surface fish. However, wild adult cavefish from the Pachón, Sabinos, Tinaja, Chica and Subterráneo caves showed highly variable responses to the three different odorant molecules they were exposed to. Pachón and Subterráneo cavefish showed the highest olfactory capacities, and Chica cavefish showed no response to the odors presented. We discuss these data with regard to the environmental conditions in which these different cavefish populations live. Our experiments in natural settings document the diversity of cave environments inhabited by a single species of cavefish, A. mexicanus, and highlight the complexity of the plastic and genetic mechanisms that underlie cave adaptation.