Eyes underground: regression of visual protein networks in subterranean mammals

Unearthing the secrets of Australia's most enigmatic and cryptic mammal, the marsupial mole

The marsupial moles are arguably Australia's most enigmatic marsupials. Almost indistinguishable from placental (eutherian) moles, they provide a striking example of convergent evolution. Exploring the genome of the southern marsupial mole, we provide insights into its unusual biology. We show definitively by retrophylogenomic analysis that marsupial moles are most closely related to bandicoots and bilbies (order Peramelemorphia). We find evidence of a marked decline in marsupial mole effective population size, most likely preceding the arrival of humans in regions near its range, and potentially corresponding to periods of climatic change. Our analysis of loss of eye function-an adaptation to subterranean life-reveals a structured order of loss of gene function associated first with the lens, then cone, and finally rod cells. Last, we identify genetic changes suggestive of adaptation to an oxygen-poor environment and of its evolution of partially descended testes.

Reduced adult stem cell fate specification led to eye reduction in cave planarians

Eye loss occurs convergently in numerous animal phyla as an adaptation to dark environments. We investigate the cave planarian Girardia multidiverticulata (Gm), a representative species of the Spiralian clade, to study mechanisms of eye loss. We found that Gm, which was previously described as an eyeless species, retains rudimentary and functional eyes. Eyes are maintained in homeostasis and regenerated in adult planarians by stem cells, called neoblasts, through their fate specification to eye progenitors. The reduced number of eye cells in cave planarians is associated with a decreased rate of stem cell fate specification to eye progenitors during homeostasis and regeneration. Conversely, the homeostatic formation of new cells from stem cell-derived progenitors for other tissues, including for neurons, pharynx, and epidermis, is comparable between cave and surface species. These findings reveal a mode of evolutionary trait loss, with change in rate of fate specification in adult stem cells leading to tissue size reduction.

Parallel Spectral Tuning of a Cone Visual Pigment Provides Evidence for Ancient Deep-Sea Adaptations in Cetaceans

Dichromatic color vision is mediated by two cone visual pigments in many eutherian mammals. After reentry into the sea, early cetaceans lost their violet-sensitive visual pigment (short wavelength-sensitive 1) independently in the baleen and toothed whale ancestors and thus obtained only monochromatic cone vision. Subsequently, losses of the middle/long wavelength-sensitive (M/LWS) pigment have also been reported in multiple whale lineages, leading to rhodopsin (RH1)-mediated rod monochromatic vision. To further elucidate the phenotypic evolution of whale visual pigments, we assessed the spectral tuning of both M/LWS and RH1 from representative cetacean taxa. Interestingly, although the coding sequences for M/LWS are intact in both the pygmy right whale and the Baird's beaked whale, no spectral sensitivity was detected in vitro. Pseudogenization of other cone vision-related genes is observed in the pygmy right whale, suggesting a loss of cone-mediated vision. After ancestral sequence reconstructions, ancient M/LWS pigments from cetacean ancestors were resurrected and functionally measured. Spectral tuning of M/LWS from the baleen whale ancestor shows that it is green sensitive, with a 40-nm shift in sensitivity to a shorter wavelength. For the ancestor of sperm whales, although no spectral sensitivity could be recorded for its M/LWS pigment, a substantial sensitivity shift (20 to 30 nm) to a shorter wavelength may have also occurred before its functional inactivation. The parallel phenotypic evolution of M/LWS to shorter wavelength sensitivity might be visual adaptations in whales allowing more frequent deep-sea activities, although additional ecological differentiations may have led to their subsequent losses.

Dynamic Expansions and Retinal Expression of Spectrally Distinct Short-Wavelength Opsin Genes in Sea Snakes

The photopigment-encoding visual opsin genes that mediate color perception show great variation in copy number and adaptive function across vertebrates. An open question is how this variation has been shaped by the interaction of lineage-specific structural genomic architecture and ecological selection pressures. We contribute to this issue by investigating the expansion dynamics and expression of the duplicated Short-Wavelength-Sensitive-1 opsin (SWS1) in sea snakes (Elapidae). We generated one new genome, 45 resequencing datasets, 10 retinal transcriptomes, and 81 SWS1 exon sequences for sea snakes, and analyzed these alongside 16 existing genomes for sea snakes and their terrestrial relatives. Our analyses revealed multiple independent transitions in SWS1 copy number in the marine Hydrophis clade, with at least three lineages having multiple intact SWS1 genes: the previously studied Hydrophis cyanocinctus and at least two close relatives of this species; Hydrophis atriceps and Hydrophis fasciatus; and an individual Hydrophis curtus. In each lineage, gene copy divergence at a key spectral tuning site resulted in distinct UV and Violet/Blue-sensitive SWS1 subtypes. Both spectral variants were simultaneously expressed in the retinae of H. cyanocinctus and H. atriceps, providing the first evidence that these SWS1 expansions confer novel phenotypes. Finally, chromosome annotation for nine species revealed shared structural features in proximity to SWS1 regardless of copy number. If these features are associated with SWS1 duplication, expanded opsin complements could be more common in snakes than is currently recognized. Alternatively, selection pressures specific to aquatic environments could favor improved chromatic distinction in just some lineages.

Diversity and Molecular Evolution of Nonvisual Opsin Genes across Environmental, Developmental, and Morphological Adaptations in Frogs

Nonvisual opsins are transmembrane proteins expressed in the eyes and other tissues of many animals. When paired with a light-sensitive chromophore, nonvisual opsins form photopigments involved in various nonvisual, light-detection functions including circadian rhythm regulation, light-seeking behaviors, and seasonal responses. Here, we investigate the molecular evolution of nonvisual opsin genes in anuran amphibians (frogs and toads). We test several evolutionary hypotheses including the predicted loss of nonvisual opsins due to nocturnal ancestry and potential functional differences in nonvisual opsins resulting from environmental light variation across diverse anuran ecologies. Using whole-eye transcriptomes of 81 species, combined with genomes, multitissue transcriptomes, and independently annotated genes from an additional 21 species, we identify which nonvisual opsins are present in anuran genomes and those that are also expressed in the eyes, compare selective constraint among genes, and test for potential adaptive evolution by comparing selection between discrete ecological classes. At the genomic level, we recovered all 18 ancestral vertebrate nonvisual opsins, indicating that anurans demonstrate the lowest documented amount of opsin gene loss among ancestrally nocturnal tetrapods. We consistently found expression of 14 nonvisual opsins in anuran eyes and detected positive selection in a subset of these genes. We also found shifts in selective constraint acting on nonvisual opsins in frogs with differing activity periods, habitats, distributions, life histories, and pupil shapes, which may reflect functional adaptation. Although many nonvisual opsins remain poorly understood, these findings provide insight into the diversity and evolution of these genes across anurans, filling an important gap in our understanding of vertebrate opsins and setting the stage for future research on their functional evolution across taxa.

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.

Diversity and Evolution of Frog Visual Opsins: Spectral Tuning and Adaptation to Distinct Light Environments

Visual systems adapt to different light environments through several avenues including optical changes to the eye and neurological changes in how light signals are processed and interpreted. Spectral sensitivity can evolve via changes to visual pigments housed in the retinal photoreceptors through gene duplication and loss, differential and coexpression, and sequence evolution. Frogs provide an excellent, yet understudied, system for visual evolution research due to their diversity of ecologies (including biphasic aquatic-terrestrial life cycles) that we hypothesize imposed different selective pressures leading to adaptive evolution of the visual system, notably the opsins that encode the protein component of the visual pigments responsible for the first step in visual perception. Here, we analyze the diversity and evolution of visual opsin genes from 93 new eye transcriptomes plus published data for a combined dataset spanning 122 frog species and 34 families. We find that most species express the four visual opsins previously identified in frogs but show evidence for gene loss in two lineages. Further, we present evidence of positive selection in three opsins and shifts in selective pressures associated with differences in habitat and life history, but not activity pattern. We identify substantial novel variation in the visual opsins and, using microspectrophotometry, find highly variable spectral sensitivities, expanding known ranges for all frog visual pigments. Mutations at spectral-tuning sites only partially account for this variation, suggesting that frogs have used tuning pathways that are unique among vertebrates. These results support the hypothesis of adaptive evolution in photoreceptor physiology across the frog tree of life in response to varying environmental and ecological factors and further our growing understanding of vertebrate visual evolution.

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.

Three Blind Moles: Molecular Evolutionary Insights on the Tempo and Mode of Convergent Eye Degeneration in Notoryctes typhlops (Southern Marsupial Mole) and Two Chrysochlorids (Golden Moles)

Golden moles (Chrysochloridae) and marsupial moles (Notoryctidae) are textbook examples of convergent evolution. Both taxa are highly adapted to subterranean lifestyles and have powerful limbs for digging through the soil/sand, ears that are adapted for low-frequency hearing, vestigial eyes that are covered by skin and fur, and the absence of optic nerve connections between the eyes and the brain. The eyes of marsupial moles also lack a lens as well as retinal rods and cones. Two hypotheses have been proposed to account for the greater degeneracy of the eyes of marsupial moles than golden moles. First, marsupial moles may have had more time to adapt to their underground habitat than other moles. Second, the eyes of marsupial moles may have been rapidly and recently vestigialized to (1) reduce the injurious effects of sand getting into the eyes and (2) accommodate the enlargement of lacrimal glands that keep the nasal cavity moist and prevent the entry of sand into the nasal passages during burrowing. Here, we employ molecular evolutionary methods on DNA sequences for 38 eye genes, most of which are eye-specific, to investigate the timing of relaxed selection (=neutral evolution) for different groups of eye-specific genes that serve as proxies for distinct functional components of the eye (rod phototransduction, cone phototransduction, lens/cornea). Our taxon sampling included 12 afrothere species, of which two are golden moles (Amblysomus hottentotus, Chrysochloris asiatica), and 28 marsupial species including two individuals of the southern marsupial mole (Notoryctes typhlops). Most of the sequences were mined from databases, but we also provide new genome data for A. hottentotus and one of the two N. typhlops individuals. Even though the eyes of golden moles are less degenerate than the eyes of marsupial moles, there are more inactivating mutations (e.g., frameshift indels, premature stop codons) in their cone phototransduction and lens/cornea genes than in orthologous genes of the marsupial mole. We estimate that cone phototransduction recovery genes were inactivated first in each group, followed by lens/cornea genes and then cone phototransduction activation genes. All three groups of genes were inactivated earlier in golden moles than in marsupial moles. For the latter, we estimate that lens/cornea genes were inactivated ~17.8 million years ago (MYA) when stem notoryctids were burrowing in the soft soils of Australian rainforests. Selection on phototransduction activation genes was relaxed much later (5.38 MYA), during the early stages of Australia's aridification that produced coastal sand plains and eventually sand dunes. Unlike cone phototransduction activation genes, rod phototransduction activation genes are intact in both golden moles and one of the two individuals of N. typhlops. A second marsupial mole individual has just a single inactivating mutation in one of the rod phototransduction activation genes (PDE6B). One explanation for this result is that some rod phototransduction activation genes are pleiotropic and are expressed in extraocular tissues, possibly in conjunction with sperm thermotaxis.

Evolution of visual guanylyl cyclases and their activating proteins with respect to clade and species-specific visual system adaptation

Membrane guanylyl cyclase receptors are important regulators of local cGMP production, critically influencing cell growth and differentiation as well as ion transport, blood pressure and calcium feedback of vertebrate phototransduction. Currently, seven different subtypes of membrane guanylyl cyclase receptors have been characterized. These receptors have tissue specific expression and are activated either by small extracellular ligands, changing CO2 concentrations or, in the case of visual guanylyl cyclases, intracellularly interacting Ca2+-dependent activating proteins. In this report, we focus on the visual guanylyl cyclase receptors (GCs) GC-E (gucy2d/e) and GC-F (gucy2f) and their activating proteins (GCAP1/2/3; guca1a/b/c). While gucy2d/e has been detected in all analyzed vertebrates, GC-F receptors are missing in several clades (reptiles, birds, and marsupials) and/or individual species. Interestingly, the absence of GC-F in highly visual sauropsida species with up to 4 different cone-opsins is compensated by an increased number of guanylyl cyclase activating proteins, whereas in nocturnal or visually impaired species with reduced spectral sensitivity it is consolidated by the parallel inactivation of these activators. In mammals, the presence of GC-E and GC-F is accompanied by the expression of one to three GCAPs, whereas in lizards and birds, up to five different GCAPs are regulating the activity of the single GC-E visual membrane receptor. In several nearly blind species, a single GC-E enzyme is often accompanied by a single variant of GCAP, suggesting that one cyclase and one activating protein are both sufficient and required for conferring the basic detection of light.

Associative colour learning and discrimination in the South African Cape rock sengi Elephantulus edwardii (Macroscelidea, Afrotheria, Mammalia)

Beside insects, sengis also consume plant material such as leaves, fruits, seeds and floral nectar. It is known that they use olfaction for foraging, but little is known about their vision and visual learning capabilities. Colour vision has been tested in two species, showing that they are likely dichromats (green- and blue-sensitive retinal cone-photoreceptors, meaning red-green colour blind). Our aim was to examine the learning and colour discrimination abilities of another species, Elephantulus edwardii. Using training procedures and choice experiments, we tested the hypotheses that the animals can associate a reward with trained colours and that they can discriminate between different colour hues. The sengis preferred the trained colours over the others, indicating associative learning. They could discriminate between all tested colours (blue, red, green, yellow). The sengis’ colour choice behaviour indicates that the animals can use also colour features to find food plant material. Additionally, learning abilities most likely are essential for the sengis’ foraging activities, for instance by associating floral or fruit shape, colour or scent with nectar or ripe fruit, to increase the efficiency to locate food sources.

Accelerated evolution of dim-light vision-related arrestin in deep-diving amniotes

Arrestins are key molecules involved in the signaling of light-sensation initiated by visual pigments in retinal photoreceptor cells. Vertebrate photoreceptor cells have two types of arrestins, rod arrestin, which is encoded by SAG and is expressed in both rods and cones, and cone arrestin, encoded by ARR3 in cones. The arrestins can bind to visual pigments, and thus regulate either dim-light vision via interactions with rhodopsin or bright-light vision together with cone visual pigments. After adapting to terrestrial life, several amniote lineages independently went back to the sea and evolved deep-diving habits. Interestingly, the rhodopsins in these species exhibit specialized phenotypes responding to rapidly changing dim-light environments. However, little is known about whether their rod arrestin also experienced adaptive evolution associated with rhodopsin. Here, we collected SAG coding sequences from >250 amniote species, and examined changes in selective pressure experienced by the sequences from deep-diving taxa. Divergent patterns of evolution of SAG were observed in the penguin, pinniped and cetacean clades, suggesting possible co-adaptation with rhodopsin. After verifying pseudogenes, the same analyses were performed for cone arrestin (ARR3) in deep-diving species and only sequences from cetacean species, and not pinnipeds or penguins, have experienced changed selection pressure compared to other species. Taken together, this evidence for changes in selective pressures acting upon arrestin genes strengthens the suggestion that rapid dim-light adaptation for deep-diving amniotes require SAG, but not ARR3.

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.

Vision-related convergent gene losses reveal SERPINE3's unknown role in the eye

Despite decades of research, knowledge about the genes that are important for development and function of the mammalian eye and are involved in human eye disorders remains incomplete. During mammalian evolution, mammals that naturally exhibit poor vision or regressive eye phenotypes have independently lost many eye-related genes. This provides an opportunity to predict novel eye-related genes based on specific evolutionary gene loss signatures. Building on these observations, we performed a genome-wide screen across 49 mammals for functionally uncharacterized genes that are preferentially lost in species exhibiting lower visual acuity values. The screen uncovered several genes, including SERPINE3, a putative serine proteinase inhibitor. A detailed investigation of 381 additional mammals revealed that SERPINE3 is independently lost in 18 lineages that typically do not primarily rely on vision, predicting a vision-related function for this gene. To test this, we show that SERPINE3 has the highest expression in eyes of zebrafish and mouse. In the zebrafish retina, serpine3 is expressed in Müller glia cells, a cell type essential for survival and maintenance of the retina. A CRISPR-mediated knockout of serpine3 in zebrafish resulted in alterations in eye shape and defects in retinal layering. Furthermore, two human polymorphisms that are in linkage with SERPINE3 are associated with eye-related traits. Together, these results suggest that SERPINE3 has a role in vertebrate eyes. More generally, by integrating comparative genomics with experiments in model organisms, we show that screens for specific phenotype-associated gene signatures can predict functions of uncharacterized genes.

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.

Comparative Genomics Provides Insights into Adaptive Evolution in Tactile-Foraging Birds

Tactile-foraging birds have evolved an enlarged principal sensory nucleus (PrV) but smaller brain regions related to the visual system, which reflects the difference in sensory dependence. The “trade-off” may exist between different senses in tactile foragers, as well as between corresponding sensory-processing areas in the brain. We explored the mechanism underlying the adaptive evolution of sensory systems in three tactile foragers (kiwi, mallard, and crested ibis). The results showed that olfaction-related genes in kiwi and mallard and hearing-related genes in crested ibis were expanded, indicating they may also have sensitive olfaction or hearing, respectively. However, some genes required for visual development were positively selected or had convergent amino acid substitutions in all three tactile branches, and it seems to show the possibility of visual degradation. In addition, we may provide a new visual-degradation candidate gene PDLIM1 who suffered dense convergent amino acid substitutions within the ZM domain. At last, two genes responsible for regulating the proliferation and differentiation of neuronal progenitor cells may play roles in determining the relative sizes of sensory areas in brain. This exploration offers insight into the relationship between specialized tactile-forging behavior and the evolution of sensory abilities and brain structures.

Gene losses may contribute to subterranean adaptations in naked mole-rat and blind mole-rat

Background

Naked mole-rats (Heterocephalus glaber, NMRs) and blind mole-rats (Spalax galili, BMRs) are representative subterranean rodents that have evolved many extraordinary traits, including hypoxia tolerance, longevity, and cancer resistance. Although multiple candidate loci responsible for these traits have been uncovered by genomic studies, many of them are limited to functional changes to amino acid sequence and little is known about the contributions of other genetic events. To address this issue, we focused on gene losses (unitary pseudogenes) and systematically analyzed gene losses in NMRs and BMRs, aiming to elucidate the potential roles of pseudogenes in their adaptation to subterranean lifestyle.

Results

We obtained the pseudogene repertoires in NMRs and BMRs, as well as their respective aboveground relatives, guinea pigs and rats, on a genome-wide scale. As a result, 167, 139, 341, and 112 pseudogenes were identified in NMRs, BMRs, guinea pigs, and rats, respectively. Functional enrichment analysis identified 4 shared and 2 species-specific enriched functional groups (EFGs) in subterranean lineages. Notably, the pseudogenes in these EFGs might be associated with either regressive (e.g., visual system) or adaptive (e.g., altered DNA damage response) traits. In addition, several pseudogenes including TNNI3K and PDE5A might be associated with specific cardiac features observed in subterranean lineages. Interestingly, we observed 20 convergent gene losses in NMRs and BMRs. Given that the functional investigations of these genes are generally scarce, we provided functional evidence that independent loss of TRIM17 in NMRs and BMRs might be beneficial for neuronal survival under hypoxia, supporting the positive role of eliminating TRIM17 function in hypoxia adaptation. Our results also suggested that pseudogenes, together with positively selected genes, reinforced subterranean adaptations cooperatively.

Conclusions

Our study provides new insights into the molecular underpinnings of subterranean adaptations and highlights the importance of gene losses in mammalian evolution.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01243-0.

Hagfish to Illuminate the Developmental and Evolutionary Origins of the Vertebrate Retina

The vertebrate eye is a vital sensory organ that has long fascinated scientists, but the details of how this organ evolved are still unclear. The vertebrate eye is distinct from the simple photoreceptive organs of other non-vertebrate chordates and there are no clear transitional forms of the eye in the fossil record. To investigate the evolution of the eye we can examine the eyes of the most ancient extant vertebrates, the hagfish and lamprey. These jawless vertebrates are in an ideal phylogenetic position to study the origin of the vertebrate eye but data on eye/retina development in these organisms is limited. New genomic and gene expression data from hagfish and lamprey suggest they have many of the same genes for eye development and retinal neurogenesis as jawed vertebrates, but functional work to determine if these genes operate in retinogenesis similarly to other vertebrates is missing. In addition, hagfish express a marker of proliferative retinal cells (Pax6) near the margin of the retina, and adult retinal growth is apparent in some species. This finding of eye growth late into hagfish ontogeny is unexpected given the degenerate eye phenotype. Further studies dissecting retinal neurogenesis in jawless vertebrates would allow for comparison of the mechanisms of retinal development between cyclostome and gnathostome eyes and provide insight into the evolutionary origins of the vertebrate eye.

A natural history of vertebrate vision loss: Insight from mammalian vision for human visual function

Research on the origin of vision and vision loss in naturally "blind" animal species can reveal the tasks that vision fulfills and the brain's role in visual experience. Models that incorporate evolutionary history, natural variation in visual ability, and experimental manipulations can help disentangle visual ability at a superficial level from behaviors linked to vision but not solely reliant upon it, and could assist the translation of ophthalmological research in animal models to human treatments. To unravel the similarities between blind individuals and blind species, we review concepts of 'blindness' and its behavioral correlates across a range of species. We explore the ancestral emergence of vision in vertebrates, and the loss of vision in blind species with reference to an evolution-based classification scheme. We applied phylogenetic comparative methods to a mammalian tree to explore the evolution of visual acuity using ancestral state estimations. Future research into the natural history of vision loss could help elucidate the function of vision and inspire innovations in how to address vision loss in humans.

Genomic evidence for the parallel regression of melatonin synthesis and signaling pathways in placental mammals

Background: The study of regressive evolution has yielded a wealth of examples where the underlying genes bear molecular signatures of trait degradation, such as pseudogenization or deletion. Typically, it appears that such disrupted genes are limited to the function of the regressed trait, whereas pleiotropic genes tend to be maintained by natural selection to support their myriad purposes. One such set of pleiotropic genes is involved in the synthesis ( AANAT, ASMT) and signaling ( MTNR1A, MTNR1B) of melatonin, a hormone secreted by the vertebrate pineal gland. Melatonin provides a signal of environmental darkness, thereby influencing the circadian and circannual rhythmicity of numerous physiological traits. Therefore, the complete loss of a pineal gland and the underlying melatonin pathway genes seems likely to be maladaptive, unless compensated by extrapineal sources of melatonin. Methods: We examined AANAT, ASMT, MTNR1A and MTNR1B in 123 vertebrate species, including pineal-less placental mammals and crocodylians. We searched for inactivating mutations and modelled selective pressures (dN/dS) to test whether the genes remain functionally intact. Results: We report that crocodylians retain intact melatonin genes and express AANAT and ASMT in their eyes, whereas all four genes have been repeatedly inactivated in the pineal-less xenarthrans, pangolins, sirenians, and whales. Furthermore, colugos have lost these genes, and several lineages of subterranean mammals have partial melatonin pathway dysfunction. These results are supported by the presence of shared inactivating mutations across clades and analyses of selection pressure based on the ratio of non-synonymous to synonymous substitutions (dN/dS), suggesting extended periods of relaxed selection on these genes. Conclusions: The losses of melatonin synthesis and signaling date to tens of millions of years ago in several lineages of placental mammals, raising questions about the evolutionary resilience of pleiotropic genes, and the causes and consequences of losing melatonin pathways in these species.

Eye-Transcriptome and Genome-Wide Sequencing for Scolecophidia: Implications for Inferring the Visual System of the Ancestral Snake

Molecular genetic data have recently been incorporated in attempts to reconstruct the ecology of the ancestral snake, though this has been limited by a paucity of data for one of the two main extant snake taxa, the highly fossorial Scolecophidia. Here we present and analyze vision genes from the first eye-transcriptomic and genome-wide data for Scolecophidia, for Anilios bicolor, and A. bituberculatus, respectively. We also present immunohistochemistry data for retinal anatomy and visual opsin-gene expression in Anilios. Analyzed in the context of 19 lepidosaurian genomes and 12 eye transcriptomes, the new genome-wide and transcriptomic data provide evidence for a much more reduced visual system in Anilios than in non-scolecophidian (=alethinophidian) snakes and in lizards. In Anilios, there is no evidence of the presence of 7 of the 12 genes associated with alethinophidian photopic (cone) phototransduction. This indicates extensive gene loss and many of these candidate gene losses occur also in highly fossorial mammals with reduced vision. Although recent phylogenetic studies have found evidence for scolecophidian paraphyly, the loss in Anilios of visual genes that are present in alethinophidians implies that the ancestral snake had a better-developed visual system than is known for any extant scolecophidian.

Phylogenomic relationships and molecular convergences to subterranean life in rodent family Spalacidae

All extant species in the rodent family Spalacidae are subterranean and have evolved various traits for underground life. However, the phylogenomic relationships among its three subfamilies (Myospalacinae, Spalacinae, and Rhizomyinae) and the molecular basis underlying their adaptations to underground life remain poorly understood. Here, we inferred the phylogenomic relationships among these subfamilies based on de novo sequencing the genome of the hoary bamboo rat ( Rhizomys pruinosus). Analyses showed that ~50% of the identified 11 028 one-to-one orthologous protein-coding genes and the concatenated sequences of these orthologous genes strongly supported a sister relationship between Myospalacinae and Rhizomyinae. The three subfamilies diversified from each other within ~2 million years. Compared with the non-subterranean controls with similar divergence dates, the spalacids shared more convergent genes with the African subterranean mole-rats at the genomic scale due to more rapid protein sequence evolution. Furthermore, these convergent genes were enriched in the functional categories of carboxylic acid transport, vascular morphogenesis, and response to oxidative stress, which are closely associated with adaptations to the hypoxic-hypercapnic underground environment. Our study presents a well-supported phylogenomic relationship among the three subfamilies of Spalacidae and offers new insights into the molecular adaptations of spalacids living underground.

Selective Gene Loss of Visual and Olfactory Guanylyl Cyclase Genes Following the Two Rounds of Vertebrate-Specific Whole-Genome Duplications

Photoreceptors convey visual information and come in two flavors; dim-light and bright-light dedicated rod and cones. Both cell types feature highly specialized phototransduction cascades that convert photonic energy into intracellular signals. Although a substantial amount of phototransduction gene ohnologs are expressed either in rods or cones, visual guanylyl cyclases (GCs) involved in the calcium (Ca²⁺) dependent feedback regulation of phototransduction are neither rod nor cone specific. The co-existence of visual GCs in both photoreceptor types suggests that specialization of these ohnologs occurred despite their overlapping expression.

Here, we analyze gene retention and inactivation patterns of vertebrate visual and closely related olfactory GCs following two rounds (2R) of vertebrate-specific whole-genome duplication events (2R WGD). Although eutherians generally use two visual and one olfactory GC, independent inactivation occurred in some lineages. Sauropsids (birds, lizards, snakes, turtles, and crocodiles) generally have only one visual GC (GC-E). Additionally, turtles (testodes) also lost the olfactory GC (GC-D). Pseudogenization in mammals occurred in specific species/families likely according to functional needs (i.e., many species with reduced vision only have GC-E). Likewise, some species not relying on scent marks lack GC-D, the olfactory GC enzyme. Interestingly, in the case of fish, no species can be found with fewer than three (two visual and one olfactory) genes and the teleost-specific 3R WGD can increase this number to up to five. This suggests that vision in fish now requires at least two visual GCs.

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.

Lizards as models to explore the ecological and neuroanatomical correlates of miniaturization

Extreme body size reductions bring about unorthodox anatomical arrangements and novel ways in which animals interact with the environment. Drawing from studies of vertebrates and invertebrates, we provide a theoretical framework for miniaturization to inform hypotheses using lizards as a study system. Through this approach, we demonstrate the repeated evolution of miniaturization across 11 families and a tendency for miniaturized species to occupy terrestrial microhabitats, possibly driven by physiological constraints. Differences in gross brain morphology between two gecko species demonstrate a proportionally larger telencephalon and smaller olfactory bulbs in the miniaturized species, though more data are needed to generalize this trend. Our study brings into light the potential contributions of miniaturized lizards to explain patterns of body size evolution and its impact on ecology and neuroanatomy. In addition, our findings reveal the need to study the natural history of miniaturized species, particularly in relation to their sensory and physiological ecology.

Echolocation in soft-furred tree mice

Echolocation is the use of reflected sound to sense features of the environment. Here, we show that soft-furred tree mice (Typhlomys) echolocate based on multiple independent lines of evidence. Behavioral experiments show that these mice can locate and avoid obstacles in darkness using hearing and ultrasonic pulses. The proximal portion of their stylohyal bone fuses with the tympanic bone, a form previously only seen in laryngeally echolocating bats. Further, we found convergence of hearing-related genes across the genome and of the echolocation-related gene prestin between soft-furred tree mice and echolocating mammals. Together, our findings suggest that soft-furred tree mice are capable of echolocation, and thus are a new lineage of echolocating mammals.

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

Evolution sometimes proceeds by loss, especially when structures and genes become dispensable after an environmental shift relaxes functional constraints. Subterranean vertebrates are outstanding models to analyze this process, and gene decay can serve as a readout. We sought to understand some general principles on the extent and tempo of the decay of genes involved in vision, circadian clock, and pigmentation in cavefishes. The analysis of the genomes of two Cuban species belonging to the genus Lucifuga provided evidence for the largest loss of eye-specific genes and nonvisual opsin genes reported so far in cavefishes. Comparisons with a recently evolved cave population of Astyanax mexicanus and three species belonging to the Chinese tetraploid genus Sinocyclocheilus revealed the combined effects of the level of eye regression, time, and genome ploidy on eye-specific gene pseudogenization. The limited extent of gene decay in all these cavefishes and the very small number of loss-of-function mutations per pseudogene suggest that their eye degeneration may not be very ancient, ranging from early to late Pleistocene. This is in sharp contrast with the identification of several vision genes carrying many loss-of-function mutations in ancient fossorial mammals, further suggesting that blind fishes cannot thrive more than a few million years in cave ecosystems.

Recapitulating Evolutionary Divergence in a Single Cis-Regulatory Element Is Sufficient to Cause Expression Changes of the Lens Gene Tdrd7

Mutations in cis-regulatory elements play important roles for phenotypic changes during evolution. Eye degeneration in the blind mole rat (BMR; Nannospalax galili) and other subterranean mammals is significantly associated with widespread divergence of eye regulatory elements, but the effect of these regulatory mutations on eye development and function has not been explored. Here, we investigate the effect of mutations observed in the BMR sequence of a conserved noncoding element upstream of Tdrd7, a pleiotropic gene required for lens development and spermatogenesis. We first show that this conserved element is a transcriptional repressor in lens cells and that the BMR sequence partially lost repressor activity. Next, we recapitulated evolutionary changes in this element by precisely replacing the endogenous regulatory element in a mouse line by the orthologous BMR sequence with CRISPR-Cas9. Strikingly, this repressor replacement caused a more than 2-fold upregulation of Tdrd7 in the developing lens; however, increased mRNA level does not result in a corresponding increase in TDRD7 protein nor an obvious lens phenotype, possibly explained by buffering at the posttranscriptional level. Our results are consistent with eye degeneration in subterranean mammals having a polygenic basis where many small-effect mutations in different eye-regulatory elements collectively contribute to phenotypic differences.

Myoglobin primary structure reveals multiple convergent transitions to semi-aquatic life in the world's smallest mammalian divers

The speciose mammalian order Eulipotyphla (moles, shrews, hedgehogs, solenodons) combines an unusual diversity of semi-aquatic, semi-fossorial, and fossorial forms that arose from terrestrial forbearers. However, our understanding of the ecomorphological pathways leading to these lifestyles has been confounded by a fragmentary fossil record, unresolved phylogenetic relationships, and potential morphological convergence, calling for novel approaches. The net surface charge of the oxygen-storing muscle protein myoglobin (ZMb), which can be readily determined from its primary structure, provides an objective target to address this question due to mechanistic linkages with myoglobin concentration. Here, we generate a comprehensive 71 species molecular phylogeny that resolves previously intractable intra-family relationships and then ancestrally reconstruct ZMb evolution to identify ancient lifestyle transitions based on protein sequence alone. Our phylogenetically informed analyses confidently resolve fossorial habits having evolved twice in talpid moles and reveal five independent secondary aquatic transitions in the order housing the world's smallest endothermic divers.

Evolutionary Constraint on Visual and Nonvisual Mammalian Opsins

Animals have evolved light-sensitive G protein-coupled receptors, known as opsins, to detect coherent and ambient light for visual and nonvisual functions. These opsins have evolved to satisfy the particular lighting niches of the organisms that express them. While many unique patterns of evolution have been identified in mammals for rod and cone opsins, far less is known about the atypical mammalian opsins. Using genomic data from over 400 mammalian species from 22 orders, unique patterns of evolution for each mammalian opsins were identified, including photoisomerases, RGR-opsin (RGR) and peropsin (RRH), as well as atypical opsins, encephalopsin (OPN3), melanopsin (OPN4), and neuropsin (OPN5). The results demonstrate that OPN5 and rhodopsin show extreme conservation across all mammalian lineages. The cone opsins, SWS1 and LWS, and the nonvisual opsins, OPN3 and RRH, demonstrate a moderate degree of sequence conservation relative to other opsins, with some instances of lineage-specific gene loss. Finally, the photoisomerase, RGR, and the best-studied atypical opsin, OPN4, have high sequence diversity within mammals. These conservation patterns are maintained in human populations. Importantly, all mammalian opsins retain key amino acid residues important for conjugation to retinal-based chromophores, permitting light sensitivity. These patterns of evolution are discussed along with known functions of each atypical opsin, such as in circadian or metabolic physiology, to provide insight into the observed patterns of evolutionary constraint.

A chromosome-level genome of Astyanax mexicanus surface fish for comparing population-specific genetic differences contributing to trait evolution

Identifying the genetic factors that underlie complex traits is central to understanding the mechanistic underpinnings of evolution. Cave-dwelling Astyanax mexicanus populations are well adapted to subterranean life and many populations appear to have evolved troglomorphic traits independently, while the surface-dwelling populations can be used as a proxy for the ancestral form. Here we present a high-resolution, chromosome-level surface fish genome, enabling the first genome-wide comparison between surface fish and cavefish populations. Using this resource, we performed quantitative trait locus (QTL) mapping analyses and found new candidate genes for eye loss such as dusp26. We used CRISPR gene editing in A. mexicanus to confirm the essential role of a gene within an eye size QTL, rx3, in eye formation. We also generated the first genome-wide evaluation of deletion variability across cavefish populations to gain insight into this potential source of cave adaptation. The surface fish genome reference now provides a more complete resource for comparative, functional and genetic studies of drastic trait differences within a species.

Design principles of gene evolution for niche adaptation through changes in protein-protein interaction networks

In contrast to fossorial and above-ground organisms, subterranean species have adapted to the extreme stresses of living underground. We analyzed the predicted protein–protein interactions (PPIs) of all gene products, including those of stress-response genes, among nine subterranean, ten fossorial, and 13 aboveground species. We considered 10,314 unique orthologous protein families and constructed 5,879,879 PPIs in all organisms using ChiPPI. We found strong association between PPI network modulation and adaptation to specific habitats, noting that mutations in genes and changes in protein sequences were not linked directly with niche adaptation in the organisms sampled. Thus, orthologous hypoxia, heat-shock, and circadian clock proteins were found to cluster according to habitat, based on PPIs rather than on sequence similarities. Curiously, "ordered" domains were preserved in aboveground species, while "disordered" domains were conserved in subterranean organisms, and confirmed for proteins in DistProt database. Furthermore, proteins with disordered regions were found to adopt significantly less optimal codon usage in subterranean species than in fossorial and above-ground species. These findings reveal design principles of protein networks by means of alterations in protein domains, thus providing insight into deep mechanisms of evolutionary adaptation, generally, and particularly of species to underground living and other confined habitats.

A fully-automated method discovers loss of mouse-lethal and human-monogenic disease genes in 58 mammals

Gene losses provide an insightful route for studying the morphological and physiological adaptations of species, but their discovery is challenging. Existing genome annotation tools focus on annotating intact genes and do not attempt to distinguish nonfunctional genes from genes missing annotation due to sequencing and assembly artifacts. Previous attempts to annotate gene losses have required significant manual curation, which hampers their scalability for the ever-increasing deluge of newly sequenced genomes. Using extreme sequence erosion (amino acid deletions and substitutions) and sister species support as an unambiguous signature of loss, we developed an automated approach for detecting high-confidence gene loss events across a species tree. Our approach relies solely on gene annotation in a single reference genome, raw assemblies for the remaining species to analyze, and the associated phylogenetic tree for all organisms involved. Using human as reference, we discovered over 400 unique human ortholog erosion events across 58 mammals. This includes dozens of clade-specific losses of genes that result in early mouse lethality or are associated with severe human congenital diseases. Our discoveries yield intriguing potential for translational medical genetics and evolutionary biology, and our approach is readily applicable to large-scale genome sequencing efforts across the tree of life.

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.

Foraging shifts and visual preadaptation in ecologically diverse bats

Changes in behaviour may initiate shifts to new adaptive zones, with physical adaptations for novel environments evolving later. While new mutations are commonly considered engines of adaptive change, sensory evolution enabling access to new resources might also arise from standing genetic diversity, and even gene loss. We examine the relative contribution of molecular adaptations, measured by positive and relaxed selection, acting on eye-expressed genes associated with shifts to new adaptive zones in ecologically diverse bats from the superfamily Noctilionoidea. Collectively, noctilionoids display remarkable ecological breadth, from highly divergent echolocation to flight strategies linked to specialized insectivory, the parallel evolution of diverse plant-based diets (e.g., nectar, pollen and fruit) from ancestral insectivory, and-unusually for echolocating bats-often have large, well-developed eyes. We report contrasting levels of positive selection in genes associated with the development, maintenance and scope of visual function, tracing back to the origins of noctilionoids and Phyllostomidae (the bat family with most dietary diversity), instead of during shifts to novel diets. Generalized plant visiting was not associated with exceptional molecular adaptation, and exploration of these novel niches took place in an ancestral phyllostomid genetic background. In contrast, evidence for positive selection in vision genes was found at subsequent shifts to either nectarivory or frugivory. Thus, neotropical noctilionoids that use visual cues for identifying food and roosts, as well as for orientation, were effectively preadapted, with subsequent molecular adaptations in nectar-feeding lineages and the subfamily Stenodermatinae of fig-eating bats fine-tuning pre-existing visual adaptations for specialized purposes.

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.

Positive Selection and Inactivation in the Vision and Hearing Genes of Cetaceans

The transition to an aquatic lifestyle in cetaceans (whales and dolphins) resulted in a radical transformation in their sensory systems. Toothed whales acquired specialized high-frequency hearing tied to the evolution of echolocation, whereas baleen whales evolved low-frequency hearing. More generally, all cetaceans show adaptations for hearing and seeing underwater. To determine the extent to which these phenotypic changes have been driven by molecular adaptation, we performed large-scale targeted sequence capture of 179 sensory genes across the Cetacea, incorporating up to 54 cetacean species from all major clades as well as their closest relatives, the hippopotamuses. We screened for positive selection in 167 loci related to vision and hearing and found that the diversification of cetaceans has been accompanied by pervasive molecular adaptations in both sets of genes, including several loci implicated in nonsyndromic hearing loss. Despite these findings, however, we found no direct evidence of positive selection at the base of odontocetes coinciding with the origin of echolocation, as found in studies examining fewer taxa. By using contingency tables incorporating taxon- and gene-based controls, we show that, although numbers of positively selected hearing and nonsyndromic hearing loss genes are disproportionately high in cetaceans, counts of vision genes do not differ significantly from expected values. Alongside these adaptive changes, we find increased evidence of pseudogenization of genes involved in cone-mediated vision in mysticetes and deep-diving odontocetes.

Losses of human disease-associated genes in placental mammals

We systematically investigate whether losses of human disease-associated genes occurred in other mammals during evolution. We first show that genes lost in any of 62 non-human mammals generally have a lower degree of pleiotropy, and are highly depleted in essential and disease-associated genes. Despite this under-representation, we discovered multiple genes implicated in human disease that are truly lost in non-human mammals. In most cases, traits resembling human disease symptoms are present but not deleterious in gene-loss species, exemplified by losses of genes causing human eye or teeth disorders in poor-vision or enamel-less mammals. We also found widespread losses of PCSK9 and CETP genes, where loss-of-function mutations in humans protect from atherosclerosis. Unexpectedly, we discovered losses of disease genes (TYMP, TBX22, ABCG5, ABCG8, MEFV, CTSE) where deleterious phenotypes do not manifest in the respective species. A remarkable example is the uric acid-degrading enzyme UOX, which we found to be inactivated in elephants and manatees. While UOX loss in hominoids led to high serum uric acid levels and a predisposition for gout, elephants and manatees exhibit low uric acid levels, suggesting alternative ways of metabolizing uric acid. Together, our results highlight numerous mammals that are 'natural knockouts' of human disease genes.

A functional enrichment test for molecular convergent evolution finds a clear protein-coding signal in echolocating bats and whales

Distantly related species entering similar biological niches often adapt by evolving similar morphological and physiological characters. How much genomic molecular convergence (particularly of highly constrained coding sequence) contributes to convergent phenotypic evolution, such as echolocation in bats and whales, is a long-standing fundamental question. Like others, we find that convergent amino acid substitutions are not more abundant in echolocating mammals compared to their outgroups. However, we also ask a more informative question about the genomic distribution of convergent substitutions by devising a test to determine which, if any, of more than 4,000 tissue-affecting gene sets is most statistically enriched with convergent substitutions. We find that the gene set most overrepresented (q-value = 2.2e-3) with convergent substitutions in echolocators, affecting 18 genes, regulates development of the cochlear ganglion, a structure with empirically supported relevance to echolocation. Conversely, when comparing to nonecholocating outgroups, no significant gene set enrichment exists. For aquatic and high-altitude mammals, our analysis highlights 15 and 16 genes from the gene sets most affected by molecular convergence which regulate skin and lung physiology, respectively. Importantly, our test requires that the most convergence-enriched set cannot also be enriched for divergent substitutions, such as in the pattern produced by inactivated vision genes in subterranean mammals. Showing a clear role for adaptive protein-coding molecular convergence, we discover nearly 2,600 convergent positions, highlight 77 of them in 3 organs, and provide code to investigate other clades across the tree of life.

Transcriptomic data support a nocturnal bottleneck in the ancestor of gecko lizards

Gecko lizards are a species-rich clade of primarily-nocturnal squamate reptiles. In geckos, adaptations to nocturnality have dramatically reshaped the eye. Perhaps the most notable change is the loss of rod cells in the retina and subsequent "transmutation" of cones into a rod-like morphology and physiology. While many studies have noted the absence of some rod-specific genes, such as the visual pigment Rhodopsin (RH1), these studies have focused on just a handful of species that are nested deep in the gecko phylogeny. Thus, it is not clear whether these changes arose through convergence, are homologous and ubiquitous across geckos, or restricted to a subset of species. Here, we used de novo eye transcriptomes from five gecko species, and genomes from two additional gecko species, representing the breadth of extant gecko diversity (i.e. 4 of the 7 gecko families, spanning the deepest divergence of crown Gekkota), to show that geckos lost expression of almost the entire suite of necessary rod-cell phototransduction genes in the eye, distinct from all other squamate reptiles. Geckos are the first vertebrate group to have lost their complete rod-cell expression pathway, not just the visual pigment. In addition, all sampled species have also lost expression of the cone-opsin SWS2 visual pigment. These results strongly suggest a single loss of rod cells and subsequent cone-to-rod transmutation that occurred prior to the diversification of extant geckos.

Complete Inactivation of Sebum-Producing Genes Parallels the Loss of Sebaceous Glands in Cetacea

Genomes are dynamic biological units, with processes of gene duplication and loss triggering evolutionary novelty. The mammalian skin provides a remarkable case study on the occurrence of adaptive morphological innovations. Skin sebaceous glands (SGs), for instance, emerged in the ancestor of mammals serving pivotal roles, such as lubrication, waterproofing, immunity, and thermoregulation, through the secretion of sebum, a complex mixture of various neutral lipids such as triacylglycerol, free fatty acids, wax esters, cholesterol, and squalene. Remarkably, SGs are absent in a few mammalian lineages, including the iconic Cetacea. We investigated the evolution of the key molecular components responsible for skin sebum production: Dgat2l6, Awat1, Awat2, Elovl3, Mogat3, and Fabp9. We show that all analyzed genes have been rendered nonfunctional in Cetacea species (toothed and baleen whales). Transcriptomic analysis, including a novel skin transcriptome from blue whale, supports gene inactivation. The conserved mutational pattern found in most analyzed genes, indicates that pseudogenization events took place prior to the diversification of modern Cetacea lineages. Genome and skin transcriptome analysis of the common hippopotamus highlighted the convergent loss of a subset of sebum-producing genes, notably Awat1 and Mogat3. Partial loss profiles were also detected in non-Cetacea aquatic mammals, such as the Florida manatee, and in terrestrial mammals displaying specialized skin phenotypes such as the African elephant, white rhinoceros and pig. Our findings reveal a unique landscape of “gene vestiges” in the Cetacea sebum-producing compartment, with limited gene loss observed in other mammalian lineages: suggestive of specific adaptations or specializations of skin lipids.

REforge Associates Transcription Factor Binding Site Divergence in Regulatory Elements with Phenotypic Differences between Species

Elucidating the genomic determinants of morphological differences between species is key to understanding how morphological diversity evolved. While differences in cis-regulatory elements are an important genetic source for morphological evolution, it remains challenging to identify regulatory elements involved in phenotypic differences. Here, we present Regulatory Element forward genomics (REforge), a computational approach that detects associations between transcription factor binding site divergence in putative regulatory elements and phenotypic differences between species. By simulating regulatory element evolution in silico, we show that this approach has substantial power to detect such associations. To validate REforge on real data, we used known binding motifs for eye-related transcription factors and identified significant binding site divergence in vision-impaired subterranean mammals in 1% of all conserved noncoding elements. We show that these genomic regions are significantly enriched in regulatory elements that are specifically active in mouse eye tissues, and that several of them are located near genes, which are required for eye development and photoreceptor function and are implicated in human eye disorders. Thus, our genome-wide screen detects widespread divergence of eye-regulatory elements and highlights regulatory regions that likely contributed to eye degeneration in subterranean mammals. REforge has broad applicability to detect regulatory elements that could be involved in many other phenotypes, which will help to reveal the genomic basis of morphological diversity.

To see or not to see: molecular evolution of the rhodopsin visual pigment in neotropical electric fishes

Functional variation in rhodopsin, the dim-light-specialized visual pigment, frequently occurs in species inhabiting light-limited environments. Variation in visual function can arise through two processes: relaxation of selection or adaptive evolution improving photon detection in a given environment. Here, we investigate the molecular evolution of rhodopsin in Gymnotiformes, an order of mostly nocturnal South American fishes that evolved sophisticated electrosensory capabilities. Our initial sequencing revealed a mutation associated with visual disease in humans. As these fishes are thought to have poor vision, this would be consistent with a possible sensory trade-off between the visual system and a novel electrosensory system. To investigate this, we surveyed rhodopsin from 147 gymnotiform species, spanning the order, and analysed patterns of molecular evolution. In contrast with our expectation, we detected strong selective constraint in gymnotiform rhodopsin, with rates of non-synonymous to synonymous substitutions lower in gymnotiforms than in other vertebrate lineages. In addition, we found evidence for positive selection on the branch leading to gymnotiforms and on a branch leading to a clade of deep-channel specialized gymnotiform species. We also found evidence that deleterious effects of a human disease-associated substitution are likely to be masked by epistatic substitutions at nearby sites. Our results suggest that rhodopsin remains an important component of the gymnotiform sensory system alongside electrolocation, and that photosensitivity of rhodopsin is well adapted for vision in dim-light environments.

Independent erosion of conserved transcription factor binding sites points to shared hindlimb, vision and external testes loss in different mammals

Genetic variation in cis-regulatory elements is thought to be a major driving force in morphological and physiological changes. However, identifying transcription factor binding events that code for complex traits remains a challenge, motivating novel means of detecting putatively important binding events. Using a curated set of 1154 high-quality transcription factor motifs, we demonstrate that independently eroded binding sites are enriched for independently lost traits in three distinct pairs of placental mammals. We show that these independently eroded events pinpoint the loss of hindlimbs in dolphin and manatee, degradation of vision in naked mole-rat and star-nosed mole, and the loss of external testes in white rhinoceros and Weddell seal. We additionally show that our method may also be utilized with more than two species. Our study exhibits a novel methodology to detect cis-regulatory mutations which help explain a portion of the molecular mechanism underlying complex trait formation and loss.

The role of ecological factors in shaping bat cone opsin evolution

Bats represent one of the largest and most striking nocturnal mammalian radiations, exhibiting many visual system specializations for performance in light-limited environments. Despite representing the greatest ecological diversity and species richness in Chiroptera, Neotropical lineages have been undersampled in molecular studies, limiting the potential for identifying signatures of selection on visual genes associated with differences in bat ecology. Here, we investigated how diverse ecological pressures mediate long-term shifts in selection upon long-wavelength (Lws) and short-wavelength (Sws1) opsins, photosensitive cone pigments that form the basis of colour vision in most mammals, including bats. We used codon-based likelihood clade models to test whether ecological variables associated with reliance on visual information (e.g. echolocation ability and diet) or exposure to varying light environments (e.g. roosting behaviour and foraging habitat) mediated shifts in evolutionary rates in bat cone opsin genes. Using additional cone opsin sequences from newly sequenced eye transcriptomes of six Neotropical bat species, we found significant evidence for different ecological pressures influencing the evolution of the cone opsins. While Lws is evolving under significantly lower constraint in highly specialized high-duty cycle echolocating lineages, which have enhanced sonar ability to detect and track targets, variation in Sws1 constraint was significantly associated with foraging habitat, exhibiting elevated rates of evolution in species that forage among vegetation. This suggests that increased reliance on echolocation as well as the spectral environment experienced by foraging bats may differentially influence the evolution of different cone opsins. Our study demonstrates that different ecological variables may underlie contrasting evolutionary patterns in bat visual opsins, and highlights the suitability of clade models for testing ecological hypotheses of visual evolution.

Duplications in Corneous Beta Protein Genes and the Evolution of Gecko Adhesion

Corneous proteins are an important component of the tetrapod integument. Duplication and diversification of keratins and associated proteins are linked with the origin of most novel integumentary structures like mammalian hair, avian feathers, and scutes covering turtle shells. Accordingly, the loss of integumentary structures often coincides with the loss of genes encoding keratin and associated proteins. For example, many hair keratins in dolphins and whales have become pseudogenes. The adhesive setae of geckos and anoles are composed of both intermediate filament keratins (IF-keratins, formerly known as alpha-keratins) and corneous beta-proteins (CBPs, formerly known as beta-keratins) and recent whole genome assemblies of two gecko species and an anole uncovered duplications in seta-specific CBPs in each of these lineages. While anoles evolved adhesive toepads just once, there are two competing hypotheses about the origin(s) of digital adhesion in geckos involving either a single origin or multiple origins. Using data from three published gecko genomes, I examine CBP gene evolution in geckos and find support for a hypothesis where CBP gene duplications are associated with the repeated evolution of digital adhesion. Although these results are preliminary, I discuss how additional gecko genome assemblies, combined with phylogenies of keratin and associated protein genes and gene duplication models, can provide rigorous tests of several hypotheses related to gecko CBP evolution. This includes a taxon sampling strategy for sequencing and assembly of gecko genomes that could help resolve competing hypotheses surrounding the origin(s) of digital adhesion.

Shifts in Selective Pressures on Snake Phototransduction Genes Associated with Photoreceptor Transmutation and Dim-Light Ancestry

The visual systems of snakes are heavily modified relative to other squamates, a condition often thought to reflect their fossorial origins. Further modifications are seen in caenophidian snakes, where evolutionary transitions between rod and cone photoreceptors, termed photoreceptor transmutations, have occurred in many lineages. Little previous work, however, has focused on the molecular evolutionary underpinnings of these morphological changes. To address this, we sequenced seven snake eye transcriptomes and utilized new whole-genome and targeted capture sequencing data. We used these data to analyze gene loss and shifts in selection pressures in phototransduction genes that may be associated with snake evolutionary origins and photoreceptor transmutation. We identified the surprising loss of rhodopsin kinase (GRK1), despite a low degree of gene loss overall and a lack of relaxed selection early during snake evolution. These results provide some of the first evolutionary genomic corroboration for a dim-light ancestor that lacks strong fossorial adaptations. Our results also indicate that snakes with photoreceptor transmutation experienced significantly different selection pressures from other reptiles. Significant positive selection was found primarily in cone-specific genes, but not rod-specific genes, contrary to our expectations. These results reveal potential molecular adaptations associated with photoreceptor transmutation and also highlight unappreciated functional differences between rod- and cone-specific phototransduction proteins. This intriguing example of snake visual system evolution illustrates how the underlying molecular components of a complex system can be reshaped in response to changing selection pressures.

Regressed but Not Gone: Patterns of Vision Gene Loss and Retention in Subterranean Mammals

Regressive evolution involves the degradation of formerly useful traits as organisms invade novel ecological niches. In animals, committing to a strict subterranean habit can lead to regression of the eyes, likely due to a limited exposure to light. Several lineages of subterranean mammals show evidence of such degeneration, which can include decreased organization of the retina, malformation of the lens, and subcutaneous positioning of the eye. Advances in DNA sequencing have revealed that this regression co-occurs with a degradation of genomic loci encoding visual functions, including protein-coding genes. Other dim light-adapted vertebrates with normal ocular anatomy, such as nocturnal and aquatic species, also demonstrate evidence of visual gene loss, but the absence of comparative studies has led to the untested assumption that subterranean mammals are special in the degree of this genomic regression. Additionally, previous studies have shown that not all vision genes have been lost in subterranean mammals, but it is unclear whether they are under relaxed selection and will ultimately be lost, are maintained due to pleiotropy or if natural selection is favoring the retention of the eye and certain critical underlying loci. Here I report that vision gene loss in subterranean mammals tends to be more extensive in quantity and differs in distribution from other dim light-adapted mammals, although some committed subterranean mammals demonstrate significant overlap with nocturnal microphthalmic species. In addition, blind subterranean mammals retain functional orthologs of non-pleiotropic visual genes that are evolving at rates consistent with purifying selection. Together, these results suggest that although living underground tends to lead to major losses of visual functions, natural selection is maintaining genes that support the eye, perhaps as an organ for circadian and/or circannual entrainment.