Genome-wide adaptive complexes to underground stresses in blind mole rats Spalax

Bank vole genomics links determinate and indeterminate growth of teeth

BACKGROUND

Continuously growing teeth are an important innovation in mammalian evolution, yet genetic regulation of continuous growth by stem cells remains incompletely understood. Dental stem cells responsible for tooth crown growth are lost at the onset of tooth root formation. Genetic signaling that initiates this loss is difficult to study with the ever-growing incisor and rooted molars of mice, the most common mammalian dental model species, because signals for root formation overlap with signals that pattern tooth size and shape (i.e., cusp patterns). Bank and prairie voles (Cricetidae, Rodentia, Glires) have evolved rooted and unrooted molars while retaining similar size and shape, providing alternative models for studying roots.

RESULTS

We assembled a de novo genome of Myodes glareolus, a vole with high-crowned, rooted molars, and performed genomic and transcriptomic analyses in a broad phylogenetic context of Glires (rodents and lagomorphs) to assess differential selection and evolution in tooth forming genes. Bulk transcriptomics comparisons of embryonic molar development between bank voles and mice demonstrated overall conservation of gene expression levels, with species-specific differences corresponding to the accelerated and more extensive patterning of the vole molar. We leverage convergent evolution of unrooted molars across the clade to examine changes that may underlie the evolution of unrooted molars. We identified 15 dental genes with changing synteny relationships and six dental genes undergoing positive selection across Glires, two of which were undergoing positive selection in species with unrooted molars, Dspp and Aqp1. Decreased expression of both genes in prairie voles with unrooted molars compared to bank voles supports the presence of positive selection and may underlie differences in root formation.

CONCLUSIONS

Our results support ongoing evolution of dental genes across Glires and identify candidate genes for mechanistic studies of root formation. Comparative research using the bank vole as a model species can reveal the complex evolutionary background of convergent evolution for ever-growing molars.

Splice Variant of Spalax Heparanase Skipping Exon 12

The subterranean blind mole rat, Spalax, has evolved significantly over 47 million years to thrive in its underground habitat. A key enzyme in this adaptation is heparanase, which degrades heparan sulfate (HS) in the extracellular matrix (ECM), facilitating angiogenesis and releasing growth factors for endothelial cells. Spalax heparanase has various splice variants influencing tumor growth and metastasis differently. We report a novel splice variant from a hypoxia-exposed kidney sample resulting from exon 12 skipping. This variant maintains the translation frame but lacks enzymatic activity, offering insights into Spalax's unique adaptations.

Thermal stress, p53 structures and learning from elephants

As species adapt to climatic changes, temperature-dependent functions of p53 in development, metabolism and cancer will adapt as well. Structural analyses of p53 epitopes interacting in response to environmental stressors, such as heat, may uncover physiologically relevant functions of p53 in cell regulation and genomic adaptations. Here we explore the multiple p53 elephant paradigm with an experimentally validated in silico model showing that under heat stress some p53 copies escape negative regulation by the MDM2 E3 ubiquitin ligase. Multiple p53 isoforms have evolved naturally in the elephant thus presenting a unique experimental system to study the scope of p53 functions and the contribution of environmental stressors to DNA damage. We assert that fundamental insights derived from studies of a historically heat-challenged mammal will provide important insights directly relevant to human biology in the light of climate change when 'heat' may introduce novel challenges to our bodies and health.

Vole genomics links determinate and indeterminate growth of teeth

Continuously growing teeth are an important innovation in mammalian evolution, yet genetic regulation of continuous growth by stem cells remains incompletely understood. Dental stem cells responsible for tooth crown growth are lost at the onset of tooth root formation. Genetic signaling that initiates this loss is difficult to study with the ever-growing incisor and rooted molars of mice, the most common mammalian dental model species, because signals for root formation overlap with signals that pattern tooth size and shape (i.e., cusp patterns). Different species of voles (Cricetidae, Rodentia, Glires) have evolved rooted and unrooted molars that have similar size and shape, providing alternative models for studying roots. We assembled a de novo genome of Myodes glareolus, a vole with high-crowned, rooted molars, and performed genomic and transcriptomic analyses in a broad phylogenetic context of Glires (rodents and lagomorphs) to assess differential selection and evolution in tooth forming genes. We identified 15 dental genes with changing synteny relationships and six dental genes undergoing positive selection across Glires, two of which were undergoing positive selection in species with unrooted molars, Dspp and Aqp1. Decreased expression of both genes in prairie voles with unrooted molars compared to bank voles supports the presence of positive selection and may underlie differences in root formation. Bulk transcriptomics analyses of embryonic molar development in bank voles also demonstrated conserved patterns of dental gene expression compared to mice, with species-specific variation likely related to developmental timing and morphological differences between mouse and vole molars. Our results support ongoing evolution of dental genes across Glires, revealing the complex evolutionary background of convergent evolution for ever-growing molars.

Relationship between Hypoxia and Hypercapnia Tolerance and Life Expectancy

The review discusses the potential relationship between hypoxia resistance and longevity, the influence of carbon dioxide on the mechanisms of aging of the mammalian organism, and intermittent hypercapnic-hypoxic effects on the signaling pathways of aging mechanisms. In the article, we focused on the potential mechanisms of the gero-protective efficacy of carbon dioxide when combined with hypoxia. The review summarizes the possible influence of intermittent hypoxia and hypercapnia on aging processes in the nervous system. We considered the perspective variants of the application of hypercapnic-hypoxic influences for achieving active longevity and the prospects for the possibilities of developing hypercapnic-hypoxic training methods.

Naked mole-rat hyaluronan

Naked mole rats (NMRs) are renowned for their exceptional longevity and remarkable maintenance of health throughout their lifetime. Their subterranean lifestyle has led to adaptations that have resulted in elevated levels of a very high molecular weight hyaluronan in their tissues. Hyaluronan, a glycosaminoglycan, is a key component of the extracellular matrix, which plays a critical role in maintaining tissue structure and regulating cell signaling pathways. This phenomenon in NMRs is attributed to a higher processing and production capacity by some of their hyaluronan synthases, along with lower degradation by certain hyaluronidases. Furthermore, this adaptation indirectly confers several advantages to NMRs, such as the preservation of skin elasticity and youthful appearance, accelerated wound healing, protection against oxidative stress, and resistance to conditions such as cancer and arthritis, largely attributable to CD44 signaling and other intricate mechanisms. Thus, the main objective of this study was to conduct a comprehensive study of the distinctive features of NMR hyaluronan, particularly emphasizing the currently known molecular mechanisms that contribute to its beneficial properties. Furthermore, this research delves into the potential applications of NMR hyaluronan in both cosmetic and therapeutic fields, as well as the challenges involved.

A Cryptic Subterranean Mammal Species, the Lesser Blind Mole Rat (Nannospalax leucodon syrmiensis)-Retreated but Not Extinct

Blind mole rats (genus Nannospalax) attract a great deal of attention because of their cancer resistance and longevity. Due to the high rate of chromosome rearrangements, 74 Nannospalax chromosomal forms have been discovered. The convergence of their external morphology complicates their taxonomy, and many cryptic species remain unrecognized. Thus, the European N. leucodon supersp. is listed in the IUCN Red List of Threatened Species with "Data Deficient" status. It is crucial for the conservation of biodiversity to clarify its taxonomy, to recognize each cryptic species, and assign to them the correct conservation status. Of the more than 20 chromosomal forms described within N. leucodon, five cryptic species occur in Serbia. The most threatened among them-N. l. syrmiensis, described and named 50 years ago in the regions of Srem, Belgrade and Mačva-has been declared extinct in the literature, which may have negative consequences for the conservation of wildlife genetic diversity. Through five years of fieldwork and comparison of 16SrRNA and MT-CYTB gene segments between old, archived teeth and recently collected material, we show that N. l. syrmiensis is not extinct. However, its habitat has been fragmented and reduced, owing primarily to anthropogenic impact. Therefore, detailed surveillance, population-structure studies, risk assessment, and appropriate conservation measures are needed.

Bat ASC2 suppresses inflammasomes and ameliorates inflammatory diseases

Bats are special in their ability to live long and host many emerging viruses. Our previous studies showed that bats have altered inflammasomes, which are central players in aging and infection. However, the role of inflammasome signaling in combating inflammatory diseases remains poorly understood. Here, we report bat ASC2 as a potent negative regulator of inflammasomes. Bat ASC2 is highly expressed at both the mRNA and protein levels and is highly potent in inhibiting human and mouse inflammasomes. Transgenic expression of bat ASC2 in mice reduced the severity of peritonitis induced by gout crystals and ASC particles. Bat ASC2 also dampened inflammation induced by multiple viruses and reduced mortality of influenza A virus infection. Importantly, it also suppressed SARS-CoV-2-immune-complex-induced inflammasome activation. Four key residues were identified for the gain of function of bat ASC2. Our results demonstrate that bat ASC2 is an important negative regulator of inflammasomes with therapeutic potential in inflammatory diseases.

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.

How voles adapt to subterranean lifestyle: Insights from RNA-seq

Life under the earth surface is highly challenging and associated with a number of morphological, physiological and behavioral modifications. Subterranean niche protects animals from predators, fluctuations in environmental parameters, but is characterized by high levels of carbon dioxide and low levels of oxygen and implies high energy requirements associated with burrowing. Moreover, it lacks most of the sensory inputs available above ground. The current study describes results from RNA-seq analysis of four subterranean voles from subfamily Arvicolinae: Prometheomys schaposchnikowi, Ellobius lutescens, Terricola subterraneus, and Lasiopodomys mandarinus. Original RNA-seq data were obtained for eight species, for nine species, SRA data were downloaded from the NCBI SRA database. Additionally assembled transcriptomes of Mynomes ochrogaster and Cricetulus griseus were included in the analysis. We searched for the selection signatures and parallel amino acid substitutions in a total of 19 species. Even within this limited data set, we found significant changes of dN/dS ratio by free-ratio model analysis for subterranean Arvicolinae. Parallel substitutions were detected in genes RAD23B and PYCR2. These genes are associated with DNA repair processes and response to oxidative stress. Similar substitutions were discovered in the RAD23 genes for highly specialized subterranean Heterocephalus glaber and Fukomys damarensis. The most pronounced signatures of adaptive evolution related to subterranean niche within species of Arvicolinae subfamily were detected for Ellobius lutescens. Our results suggest that genomic adaptations can occur very quickly so far as the amount of selection signatures was found to be compliant with the degree of specialization to the subterranean niche and independent from the evolutionary age of the taxon. We found that the number of genomic signatures of selection does not depend on the age of the taxon, but is positively correlated with the degree of specialization to the subterranean niche.

The Role of TP53 in Adaptation and Evolution

The TP53 gene is a major player in cancer formation, and it is considered the most important tumor suppressor gene. The p53 protein acts as a transcription factor, and it is involved in DNA repair, senescence, cell-cycle control, autophagy, and apoptosis. Beyond cancer, there is evidence that TP53 is associated with fertility, aging, and longevity. Additionally, more evidence exists that genetic variants in TP53 are associated with environmental adaptation. Special TP53 amino-acid residues or pathogenic TP53 mutations seem to be adaptive for animals living in hypoxic and cold environments or having been exposed to starvation, respectively. At the somatic level, it has recently been proven that multiple cancer genes, including TP53, are under positive selection in healthy human tissues. It is not clear why these driver mutations do not transform these tissues into cancerous ones. Other studies have shown that elephants have multiple TP53 copies, probably this being the reason for the very low cancer incidence in these large animals. This may explain the famous Peto's paradox. This review discusses in detail the multilevel role of TP53 in adaptation, according to the published evidence. This role is complicated, and it extends from cells to individuals and to populations.

Unusual occurrence of domestication syndrome amongst African mole-rats: Is the naked mole-rat a domestic animal?

The Naked mole-rat (NMR) is becoming a prominent model organism due to its peculiar traits, such as eusociality, extreme longevity, cancer resistance, and reduced pain sensitivity. It belongs to the African mole-rats (AMR), a family of subterranean rodents that includes solitary, cooperative breeding and eusocial species. We identified and quantified the domestication syndrome (DS) across AMR, a set of morphological and behavioural traits significantly more common and pronounced amongst domesticated animals than in their wild counterparts. Surprisingly, the NMR shows apparent DS traits when compared to the solitary AMR. Animals can self-domesticate when a reduction of the fear response is naturally selected, such as in islands with no predators, or to improve the group’s harmony in cooperative breeding species. The DS may be caused by alterations in the physiology of the neural crest cells (NCC), a transient population of cells that generate a full range of tissues during development. The NCC contribute to organs responsible for transmitting the fear response and various other tissues, including craniofacial bones. Therefore, mutations affecting the NCC can manifest as behavioural and morphological alterations in many structures across the body, as seen in neurocristopathies. We observed that all social AMRs are chisel-tooth diggers, an adaption to hard soils that requires the flattening of the skull. We hypothesise that chisel-tooth digging could impose a selective pressure on the NCC that triggered the DS’s appearance, possibly facilitating the evolution of sociality. Finally, we discuss how DS traits are neutral or beneficial for the subterranean niche, strategies to test this hypothesis and report well-studied mutations in the NMR that are associated with the NCC physiology or with the control of the fear response. In conclusion, we argue that many of the NMR’s unconventional traits are compatible with the DS and provide a hypothesis about its origins. Our model proposes a novel avenue to enhance the understanding of the extraordinary biology of the NMR.

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.

Evolutionary consequences of genomic deletions and insertions in the woolly mammoth genome

Woolly mammoths had a set of adaptations that enabled them to thrive in the Arctic environment. Many mammoth-specific single nucleotide polymorphisms (SNPs) responsible for unique mammoth traits have been previously identified from ancient genomes. However, a multitude of other genetic variants likely contributed to woolly mammoth evolution. In this study, we sequenced two woolly mammoth genomes and combined these with previously sequenced mammoth and elephant genomes to conduct a survey of mammoth-specific deletions and indels. We find that deletions are highly enriched in non-coding regions, suggesting selection against structural variants that affect protein sequences. Nonetheless, at least 87 woolly mammoth genes contain deletions or indels that modify the coding sequence, including genes involved in skeletal morphology and hair growth. These results suggest that deletions and indels contributed to the unique phenotypic adaptations of the woolly mammoth, and were potentially critical to surviving in its natural environment.

•

Two new high-quality woolly mammoth genomes have been generated

•

A new method was used to identify deletions and insertions in woolly mammoths

•

At least 87 genes have been affected by deletions or indels in the mammoth lineage

•

Genes involved in skeletal morphology and hair growth are affected by deletions

Response to Different Oxygen Partial Pressures and Evolution Analysis of Apoptosis-Related Genes in Plateau Zokor (Myospalax baileyi)

The plateau zokor (Myospalax baileyi) is a native species of the Qinghai–Tibet Plateau that spends its entire life underground in sealed burrows with hypoxic conditions. The present study aimed to assess the sequence characteristics of apoptosis-related genes and the response to different oxygen partial pressures (pO₂) in plateau zokor and Sprague-Dawley rats. The sequences of the p53-induced protein with a death domain (Pidd), p53-upregulated modulator of apoptosis (Puma), insulin-like growth factor binding protein 3 (Igfbp3), and apoptosis protease-activating factor 1 (Apaf1) were evaluated concerning homology and convergent evolution sites, and their mRNA levels were evaluated in different tissues under 14.13 (3,300 m) and 16.12 kPa (2,260 m) pO₂ conditions. Our results showed that, (1) the sequences of the apoptosis-related genes in plateau zokor were highly similar to those of Nannospalax galili, followed by Rattus norvegicus; (2). Pidd, Puma, Igfbp3, and Apaf1 of plateau zokor were found to have five, one, two, and five convergent sites in functional domains with N. galili, respectively. Lastly (3), under low pO₂, the expression of Pidd and Puma was downregulated in the lung of plateau zokors. In turn, Igfbp3 and Apaf1 were upregulated in the liver and lung, and Puma was upregulated in the skeletal muscle of plateau zokor under low pO₂. In Sprague-Dawley rats, low pO₂ downregulated Puma and Apaf1 expression in the liver and downregulated Igfbp3 and Puma in the lung and skeletal muscle separately. In contrast, low pO₂ upregulated Pidd expression in the liver and skeletal muscle of Sprague-Dawley rats. Overall, the expression patterns of Apaf1, Igfbp3, and Puma showed the opposite pattern in the liver, lung, and skeletal muscle, respectively, of plateau zokor as compared with Sprague-Dawley rats. In conclusion, for the long-time adaptation to hypoxic environments, Pidd, Puma, Igfbp3, and Apaf1 of plateau zokor underwent convergent evolution, which we believe may have led to upregulation of their levels under low oxygen partial pressures to induce apoptosis, so as to suppress tumorigenesis under hypoxic environments in plateau zokor.

Alone, in the dark: The extraordinary neuroethology of the solitary blind mole rat

On the social scale, the blind mole rat (BMR; Spalax ehrenbergi) is an extreme. It is exceedingly solitary, territorial, and aggressive. BMRs reside underground, in self-excavated tunnels that they rarely leave. They possess specialized sensory systems for social communication and navigation, which allow them to cope with the harsh environmental conditions underground. This review aims to present the blind mole rat as an ideal, novel neuroethological model for studying aggressive and solitary behaviors. We discuss the BMR's unique behavioral phenotype, particularly in the context of 'anti-social' behaviors, and review the available literature regarding its specialized sensory adaptations to the social and physical habitat. To date, the neurobiology of the blind mole rat remains mostly unknown and holds a promising avenue for scientific discovery. Unraveling the neural basis of the BMR's behavior, in comparison to that of social rodents, can shed important light on the underlying mechanisms of psychiatric disorders in humans, in which similar behaviors are displayed.

Genomic insights into zokors' phylogeny and speciation in China

SignificancePhylogenies are the basis of many ecological and evolutionary studies. However, zokor phylogeny and speciation patterns are heavily debated. This study disentangled the phylogeny and speciation of zokors genomically. Six species of the Eospalax were separated into high-altitude E. baileyi and E. smithi and the rest four low-altitude species by recent Qinghai-Tibet Plateau uplift 3.6 million y ago. E. rothschildi and E. smithi speciated south of the Qinling-Huaihe Line, where refuges were supplied during glaciation. Introgression and incomplete lineage sorting led to the complex phylogeny of zokors. Genomic islands were formed due to ancient polymorphisms and divergence hitchhiking. This study concluded that climatic, geological, and tectonic events shaped the phylogeny and speciation of zokors in China.

Cryptic Diversity of the European Blind Mole Rat Nannospalax leucodon Species Complex: Implications for Conservation

We explored the cryptic speciation of the Nannospalax leucodon species complex, characterised by intense karyotype evolution and reduced phenotypic variability that has produced different lineages, out of which 25 are described as chromosomal forms (CFs), so many cryptic species remain unnoticed. Although some of them should be classified as threatened, they lack the official nomenclature necessary to be involved in conservation strategies. Reproductive isolation between seven CFs has previously been demonstrated. To investigate the amount and dynamics of genetic discrepancy that follows chromosomal changes, infer speciation levels, and obtain phylogenetic patterns, we analysed mitochondrial 16S rRNA and MT-CYTB nucleotide polymorphism among 17 CFs—the highest number studied so far. Phylogenetic trees delineated 11 CFs as separate clades. Evolutionary divergence values overlapped with acknowledged higher taxonomic categories, or sometimes exceeded them. The fact that CFs with higher 2n are evolutionary older corresponds to the fusion hypothesis of Nannospalax karyotype evolution. To participate in conservation strategies, N. leucodon classification should follow the biological species concept, and proposed cryptic species should be formally named, despite a lack of classical morphometric discrepancy. We draw attention towards the syrmiensis and montanosyrmiensis CFs, estimated to be endangered/critically endangered, and emphasise the need for detailed monitoring and population survey for other cryptic species.

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.

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.

Adaptation of mammals to hypoxia

Oxygen plays a pivotal role in the metabolism and activities of mammals. However, oxygen is restricted in some environments-subterranean burrow systems or habitats at high altitude or deep in the ocean-and this could exert hypoxic stresses such as oxidative damage on organisms living in these environments. In order to cope with these stresses, organisms have evolved specific strategies to adapt to hypoxia, including changes in physiology, gene expression regulation, and genetic mutations. Here, we review how mammals have adapted to the three high-altitude plateaus of the world, the limited oxygen dissolved in deep water habitats, and underground tunnels, with the aim of better understanding the adaptation of mammals to hypoxia.

The hypoxia adaptation of small mammals to plateau and underground burrow conditions

Oxygen is one of the important substances for the survival of most life systems on the earth, and plateau and underground burrow systems are two typical hypoxic environments. Small mammals living in hypoxic environments have evolved different adaptation strategies, which include increased oxygen delivery, metabolic regulation of physiological responses and other physiological responses that change tissue oxygen utilization. Multi-omics predictions have also shown that these animals have evolved different adaptations to extreme environments. In particular, vascular endothelial growth factor (VEGF) and erythropoietin (EPO), which have specific functions in the control of O2 delivery, have evolved adaptively in small mammals in hypoxic environments. Naked mole-rats and blind mole-rats are typical hypoxic model animals as they have some resistance to cancer. This review primarily summarizes the main living environment of hypoxia tolerant small mammals, as well as the changes of phenotype, physiochemical characteristics and gene expression mode of their long-term living in hypoxia environment.

Revelations About Aging and Disease from Unconventional Vertebrate Model Organisms

Aging is a major risk factor for multiple diseases. Understanding the underlying mechanisms of aging would help to delay and prevent age-associated diseases. Short-lived model organisms have been extensively used to study the mechanisms of aging. However, these short-lived species may be missing the longevity mechanisms that are needed to extend the lifespan of an already long-lived species such as humans. Unconventional long-lived animal species are an excellent resource to uncover novel mechanisms of longevity and disease resistance. Here, we review mechanisms that evolved in nonmodel vertebrate species to counteract age-associated diseases. Some antiaging mechanisms are conserved across species; however, various nonmodel species also evolved unique mechanisms to delay aging and prevent disease. This variety of antiaging mechanisms has evolved due to the remarkably diverse habitats and behaviors of these species. We propose that exploring a wider range of unconventional vertebrates will provide important resources to study antiaging mechanisms that are potentially applicable to humans.

Lifespan Extension in Long-Lived Vertebrates Rooted in Ecological Adaptation

Contemporary studies on aging and longevity have largely overlooked the role that adaptation plays in lifespan variation across species. Emerging evidence indicates that the genetic signals of extended lifespan may be maintained by natural selection, suggesting that longevity could be a product of organismal adaptation. The mechanisms of adaptation in long-lived animals are believed to account for the modification of physiological function. Here, we first review recent progress in comparative biology of long-lived animals, together with the emergence of adaptive genetic factors that control longevity and disease resistance. We then propose that hitchhiking of adaptive genetic changes is the basis for lifespan changes and suggest ways to test this evolutionary model. As individual adaptive or adaptation-linked mutations/substitutions generate specific forms of longevity effects, the cumulative beneficial effect is largely nonrandom and is indirectly favored by natural selection. We consider this concept in light of other proposed theories of aging and integrate these disparate ideas into an adaptive evolutionary model, highlighting strategies in decoding genetic factors of lifespan control.

Genomic analysis unveils mechanisms of northward invasion and signatures of plateau adaptation in the Asian house rat

The Asian house rat (AHR), Rattus tanezumi, has recently invaded the northern half of China. The AHR is a highly adaptive rat species that has also successfully conquered the Qinghai-Tibet Plateau (QTP) and replaced the brown rat (BR), R. norvegicus, at the edge of the QTP. Here, we assembled a draft genome of the AHR and explored the mechanisms of its northward invasion and the genetic basis underlying plateau adaptation in this species. Population genomic analyses revealed that the northwardly invasive AHRs consisted of two independent and genetically distinct populations which might result from multiple independent primary invasion events. One invasive population exhibited reduced genetic diversity and distinct population structure compared with its source population, while the other displayed preserved genetic polymorphisms and little genetic differentiation from its source population. Genes involved in G-protein coupled receptors and carbohydrate metabolism may contribute to the local adaptation of northern AHRs. In particular, RTN4 was identified as a key gene for AHRs in the QTP that favours adaptation to high-altitude hypoxia. Coincidently, the physiological performance and transcriptome profiles of hypoxia-exposed rats both showed better hypoxia adaptation in AHRs than in BRs that failed to colonize the heart of the QTP, which may have facilitated the replacement of the BR population by the invading AHRs at the edge of the QTP. This study provides profound insights into the multiple origins of the northwardly invasive AHR and the great tolerance to hypoxia in this species.

Transposon-triggered innate immune response confers cancer resistance to the blind mole rat

Blind mole rats (BMRs) are small rodents, characterized by exceptionally long lifespan (> 21 years) and resistance to both spontaneous and induced tumorigenesis. Here we report that cancer resistance in the BMR is mediated by retrotransposable elements (RTEs). BMR cells and tissues express very low levels of DNA methyltransferase 1 (DNMT1). Upon cell hyperplasia, the BMR genome DNA loses methylation, resulting in activation of RTEs. Up-regulated RTEs form cytoplasmic RNA/DNA hybrids, which activate cGAS-STING pathway to induce cell death. Although this mechanism is enhanced in the BMR, we show that it functions in mice and human. We propose that RTEs were coopted to serve as tumor suppressors that monitor cell proliferation and are activated in premalignant cells to trigger cell death via activation of innate immune response. RTEs activation is a double-edged sword, serving as a tumor suppressor but in late life contributing to aging via induction of sterile inflammation.

Local adaptation of bitter taste and ecological speciation in a wild mammal

Sensory systems are attractive evolutionary models to address how organisms adapt to local environments that can cause ecological speciation. However, tests of these evolutionary models have focused on visual, auditory, and olfactory senses. Here, we show local adaptation of bitter taste receptor genes in two neighboring populations of a wild mammal—the blind mole rat Spalax galili—that show ecological speciation in divergent soil environments. We found that basalt-type bitter receptors showed higher response intensity and sensitivity compared with chalk-type ones using both genetic and cell-based functional analyses. Such functional changes could help animals adapted to basalt soil select plants with less bitterness from diverse local foods, whereas a weaker reception to bitter taste may allow consumption of a greater range of plants for animals inhabiting chalk soil with a scarcity of food supply. Our study shows divergent selection on food resources through local adaptation of bitter receptors, and suggests that taste plays an important yet underappreciated role in speciation.

Cloning of two splice variants of Spalax heparanase encoding for truncated proteins

Heparanase is an endoglycosidase that degrades heparan sulfate side chains of heparan sulfate-proteoglycans. It liberates heparan sulfate-bound growth factors and thereby promotes blood vessel sprouting and angiogenesis. The subterranean blind mole rat, Spalax, is a wild mammal that lives most of its life in underground tunnels where it experiences sharp fluctuations in oxygen and carbon dioxide levels. We described two splice variants of heparanase from Spalax, Splice 7 and splice 36, both devoid of heparanase enzymatic activity. Splice 7 increases tumor growth, while splice 36 functions as a dominant negative to wild-type heparanase and decreases tumor growth and metastasis. Here, we describe two novel splice variants of Spalax heparanase, splice 67 and splice 612. These splice variants result in production of a shorter heparanase proteins that are similar to the wild-type native heparanase in their N-terminal but have unique C-terminals. Both splice 67 and 612 lack heparan sulfate degradation activity.

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.

Evolution of the Autism-Associated Neuroligin-4 Gene Reveals Broad Erosion of Pseudoautosomal Regions in Rodents

Variants in genes encoding synaptic adhesion proteins of the neuroligin family, most notably neuroligin-4, are a significant cause of autism spectrum disorders in humans. Although human neuroligin-4 is encoded by two genes, NLGN4X and NLGN4Y, that are localized on the X-specific and male-specific regions of the two sex chromosomes, the chromosomal localization and full genomic sequence of the mouse Nlgn4 gene remain elusive. Here, we analyzed the neuroligin-4 genes of numerous rodent species by direct sequencing and bioinformatics, generated complete drafts of multiple rodent neuroligin-4 genes, and examined their evolution. Surprisingly, we find that the murine Nlgn4 gene is localized to the pseudoautosomal region (PAR) of the sex chromosomes, different from its human orthologs. We show that the sequence differences between various neuroligin-4 proteins are restricted to hotspots in which rodent neuroligin-4 proteins contain short repetitive sequence insertions compared with neuroligin-4 proteins from other species, whereas all other protein sequences are highly conserved. Evolutionarily, these sequence insertions initiate in the clade eumuroidea of the infraorder myomorpha and are additionally associated with dramatic changes in noncoding sequences and gene size. Importantly, these changes are not exclusively restricted to neuroligin-4 genes but reflect major evolutionary changes that substantially altered or even deleted genes from the PARs of both sex chromosomes. Our results show that despite the fact that the PAR in rodents and the neuroligin-4 genes within the rodent PAR underwent massive evolutionary changes, neuroligin-4 proteins maintained a highly conserved core structure, consistent with a substantial evolutionary pressure preserving its physiological function.

Regulation of HIF-1α and p53 in stress responses in the subterranean rodents Lasiopodomys mandarinus and Lasiopodomys brandtii (Rodentia: Cricetidae)

The response mechanism and interaction patterns of HIF-1α and p53 in animals in an hypoxic environment are crucial for their hypoxic tolerance and adaptation. Many studies have shown that underground rodents have better hypoxic adaptation characteristics. However, the mechanism by which HIF-1α and p53 in underground rodents respond to hypoxic environments compared with in ground rodents remains unclear. Further, whether a synergy between HIF-1α and p53 enables animals tolerate extremely hypoxic environments is unclear. We studied HIF-1α and p53 expression in the brain tissue and cell apoptosis in the hippocampal CA1 region during 6 hours of acute hypoxia (5% oxygen) in Lasiopodomys mandarinus (Milne-Edwards, 1871) and Lasiopodomys brandtii (Radde, 1861), two closely related small rodents with different life characteristics (underground and aboveground, respectively), using a comparative biology method to determine the mechanisms underlying their adaptation to this environment. Our results indicate that HIF-1α and p53 expression is more rapid in L. mandarinus than in L. brandtii under acute hypoxic environments, resulting in a significant synergistic effect in L. mandarinus. Correlation analysis revealed that HIF-1α expression and the apoptotic index of the hippocampal CA1 regions of the brain tissues of L. mandarinus and L. brandtii, both under hypoxia, were significantly negatively and positively correlated, respectively. Long-term existence in underground burrow systems could enable better adaptation to hypoxia in L. mandarinus than in L. brandtii. We speculate that L. mandarinus can quickly eliminate resulting damage via the synergistic effect of p53 and HIF-1α in response to acute hypoxic environments, helping the organism quickly return to a normal state after the stress.

Distinct organization of adaptive immunity in the long-lived rodent Spalax galili

A balanced immune response is a cornerstone of healthy aging. Here, we uncover distinctive features of the long-lived blind mole-rat (Spalax spp.) adaptive immune system, relative to humans and mice. The T-cell repertoire remains diverse throughout the Spalax lifespan, suggesting a paucity of large long-lived clones of effector-memory T cells. Expression of master transcription factors of T-cell differentiation, as well as checkpoint and cytotoxicity genes, remains low as Spalax ages. The thymus shrinks as in mice and humans, while interleukin-7 and interleukin-7 receptor expression remains high, potentially reflecting the sustained homeostasis of naive T cells. With aging, immunoglobulin hypermutation level does not increase and the immunoglobulin-M repertoire remains diverse, suggesting shorter B-cell memory and sustained homeostasis of innate-like B cells. The Spalax adaptive immune system thus appears biased towards sustained functional and receptor diversity over specialized, long-lived effector-memory clones-a unique organizational strategy that potentially underlies this animal's extraordinary longevity and healthy aging.

From immune to olfactory expression: neofunctionalization of formyl peptide receptors

Variations in gene expression patterns represent a powerful source of evolutionary innovation. In a rodent living about 70 million years ago, a genomic accident led an immune formyl peptide receptor (FPR) gene to hijack a vomeronasal receptor regulatory sequence. This gene shuffling event forced an immune pathogen sensor to transition into an olfactory chemoreceptor, which thus moved from sensing the internal world to probing the outside world. We here discuss the evolution of the FPR gene family, the events that led to their neofunctionalization in the vomeronasal organ and the functions of immune and vomeronasal FPRs.

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.

Microbiome and Longevity: High Abundance of Longevity-Linked Muribaculaceae in the Gut of the Long-Living Rodent Spalax leucodon

With a world population living longer as well as marked disparities in life expectancy, understanding the determinants of longevity is one of the priority research agendas in 21st century life sciences. To this end, the blind mole-rat (Spalax leucodon), a subterranean mammalian, has emerged as an exceptional model organism due to its astonishing features such as remarkable longevity, hypoxia and hypercapnia tolerance, and cancer resistance. The microbiome has been found to be a vital parameter for cellular physiology and it is safe to assume that it has an impact on life expectancy. Although the unique characteristics of Spalax make it an ideal experimental model for longevity research, there is limited knowledge of the bacterial composition of Spalax microbiome, which limits its in-depth utilization. In this study, using 16S rRNA amplicon sequencing, we report the gut and skin bacterial structure of Spalax for the first time. The diversity between fecal and skin samples was manifested in the distant clustering, as revealed by beta diversity analysis. Importantly, the longevity-linked Muribaculaceae bacterial family was found to be the dominating bacterial taxa in Spalax fecal samples. These new findings contribute toward further development of Spalax as a model for longevity research and potential linkages between microbiome composition and longevity.

Genome-wide adaptive evolution to underground stresses in subterranean mammals: Hypoxia adaption, immunity promotion, and sensory specialization

Life underground has provided remarkable examples of adaptive evolution in subterranean mammals; however, genome-wide adaptive evolution to underground stresses still needs further research. There are approximately 250 species of subterranean mammals across three suborders and six families. These species not only inhabit hypoxic and dark burrows but also exhibit evolved adaptation to hypoxia, cancer resistance, and specialized sensory systems, making them an excellent model of evolution. The adaptive evolution of subterranean mammals has attracted great attention and needs further study. In the present study, phylogenetic analysis of 5,853 single-copy orthologous gene families of five subterranean mammals (Nannospalax galili, Heterocephalus glaber, Fukomys damarensis, Condylura cristata, and Chrysochloris asiatica) showed that they formed fou distinct clusters. This result is consistent with the traditional systematics of these species. Furthermore, comparison of the high-quality genomes of these five subterranean mammalian species led to the identification of the genomic signatures of adaptive evolution. Our results show that the five subterranean mammalian did not share positively selected genes but had similar functional enrichment categories, including hypoxia tolerance, immunity promotion, and sensory specialization, which adapted to the environment of underground stresses. Moreover, variations in soil hardness, climate, and lifestyles have resulted in different molecular mechanisms of adaptation to the hypoxic environment and different degrees of visual degradation. These results provide insights into the genome-wide adaptive evolution to underground stresses in subterranean mammals, with special focus on the characteristics of hypoxia adaption, immunity promotion, and sensory specialization response to the life underground.

Comparative morphometric analysis of lungs of the semifossorial giant pouched rat (Cricetomys gambianus) and the subterranean Nigerian mole rat (Cryptomys foxi)

Lungs of the rodent species, the African giant pouched rat (Cricetomys gambianus) and the Nigerian mole rat (Cryptomys foxi) were investigated. Significant morphometric differences exist between the two species. The volume of the lung per unit body mass was 2.7 times larger; the respiratory surface area 3.4 times greater; the volume of the pulmonary capillary blood 2 times more; the harmonic mean thickness of the blood-gas (tissue) barrier (τht) ~29% thinner and; the total pulmonary morphometric diffusing capacity (DLo2) for O2 2.3 times more in C. foxi. C. gambianus occupies open burrows that are ventilated with air while C. foxi lives in closed burrows. The less morphometrically specialized lungs of C. gambianus may be attributed to its much larger body mass (~6 times more) and possibly lower metabolic rate and its semifossorial life whereas the 'superior' lungs of C. foxi may largely be ascribed to the subterranean hypoxic and hypercapnic environment it occupies. Compared to other rodents species that have been investigated hitherto, the τht was mostly smaller in the lungs of the subterranean species and C. foxi has the highest mass-specific DLo2. The fossorial- and the subterranean rodents have acquired various pulmonary structural specializations that relate to habitats occupied.

Proteomics of Long-Lived Mammals

Mammalian species differ up to 100-fold in their aging rates and maximum lifespans. Long-lived mammals appear to possess traits that extend lifespan and healthspan. Genomic analyses have not revealed a single pro-longevity function that would account for all longevity effects. In contrast, it appears that pro-longevity mechanisms may be complex traits afforded by connections between metabolism and protein functions that are impossible to predict by genomic approaches alone. Thus, metabolomics and proteomics studies will be required to understand the mechanisms of longevity. Several examples are reviewed that demonstrate the naked mole rat (NMR) shows unique proteomic signatures that contribute to longevity by overcoming several hallmarks of aging. SIRT6 is also discussed as an example of a protein that evolves enhanced enzymatic function in long-lived species. Finally, it is shown that several longevity-related proteins such as Cip1/p21, FOXO3, TOP2A, AKT1, RICTOR, INSR, and SIRT6 harbor posttranslational modification (PTM) sites that preferentially appear in either short- or long-lived species and provide examples of crosstalk between PTM sites. Prospects of enhancing lifespan and healthspan of humans by altering metabolism and proteoforms with drugs that mimic changes observed in long-lived species are discussed.

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.

Convergent reduction of V1R genes in subterranean rodents

BACKGROUND

Vomeronasal type 1 receptor genes (V1Rs) are expected to detect intraspecific pheromones. It is believed that rodents rely heavily on pheromonal communication mediated by V1Rs, but pheromonal signals are thought to be confined in subterranean rodents that live in underground burrows. Thus, subterranean rodents may show a contrasting mode of V1R evolution compared with their superterranean relatives.

RESULTS

We examined the V1R evolution in subterranean rodents by analyzing currently available genomes of 24 rodents, including 19 superterranean and 5 subterranean species from three independent lineages. We identified a lower number of putatively functional V1R genes in each subterranean rodent (a range of 22-40) compared with superterranean species (a range of 63-221). After correcting phylogenetic inertia, the positive correlation remains significant between the small V1R repertoire size and the subterranean lifestyle. To test whether V1Rs have been relaxed from functional constraints in subterranean rodents, we sequenced 22 intact V1Rs in 29 individuals of one subterranean rodent (Spalax galili) from two soil populations, which have been proposed to undergo incipient speciation. We found 12 of the 22 V1Rs to show significant genetic differentiations between the two natural populations, indicative of diversifying selection.

CONCLUSION

Our study demonstrates convergent reduction of V1Rs in subterranean rodents from three independent lineages. Meanwhile, it is noteworthy that most V1Rs in the two Spalax populations are under diversifying selection rather than relaxed selection, suggesting that functional constraints on these genes may have retained in some subterranean species.

TFforge utilizes large-scale binding site divergence to identify transcriptional regulators involved in phenotypic differences

Changes in gene regulation are important for phenotypic and in particular morphological evolution. However, it remains challenging to identify the transcription factors (TFs) that contribute to differences in gene regulation and thus to phenotypic differences between species. Here, we present TFforge (Transcription Factor forward genomics), a computational method to identify TFs that are involved in the loss of phenotypic traits. TFforge screens an input set of regulatory genomic regions to detect TFs that exhibit a significant binding site divergence signature in species that lost a particular phenotypic trait. Using simulated data of modular and pleiotropic regulatory elements, we show that TFforge can identify the correct TFs for many different evolutionary scenarios. We applied TFforge to available eye regulatory elements to screen for TFs that exhibit a significant binding site decay signature in subterranean mammals. This screen identified interacting and co-binding eye-related TFs, and thus provides new insights into which TFs likely contribute to eye degeneration in these species. TFforge has broad applicability to identify the TFs that contribute to phenotypic changes between species, and thus can help to unravel the gene-regulatory differences that underlie phenotypic evolution.

Comparative genomics suggests loss of keratin K24 in three evolutionary lineages of mammals

Keratins are the main cytoskeletal proteins of epithelial cells and changes in the expression of keratins have contributed to the evolutionary adaptation of epithelia to different environments. Keratin K24 was proposed to be a differentiation marker of epidermal keratinocytes but the significance of K24 expression in the epidermis versus other tissues has remained elusive. Here, we show by RT-PCR, western blot, and immunofluorescence analyses that K24 is highly expressed in the epithelium of the cornea whereas its expression levels are significantly lower in other stratified epithelia including in the epidermis. To investigate the evolutionary history of K24, we screened the genome sequences of vertebrates for orthologs of the human KRT24 gene. The results of this comparative genomics study suggested that KRT24 originated in a common ancestor of amniotes and that it was lost independently in three clades of mammals, i.e. camels, cetaceans, and a subclade of pinnipeds comprising eared seals and the walrus. Together, the results of this study identify K24 as component of the cytoskeleton in the human corneal epithelium and reveal previously unknown differences of keratin gene content among mammalian species.

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.

What lies beneath? Molecular evolution during the radiation of caecilian amphibians

Background

Evolution leaves an imprint in species through genetic change. At the molecular level, evolutionary changes can be explored by studying ratios of nucleotide substitutions. The interplay among molecular evolution, derived phenotypes, and ecological ranges can provide insights into adaptive radiations. Caecilians (order Gymnophiona), probably the least known of the major lineages of vertebrates, are limbless tropical amphibians, with adults of most species burrowing in soils (fossoriality). This enigmatic order of amphibians are very distinct phenotypically from other extant amphibians and likely from the ancestor of Lissamphibia, but little to nothing is known about the molecular changes underpinning their radiation. We hypothesised that colonization of various depths of tropical soils and of freshwater habitats presented new ecological opportunities to caecilians.

Results

A total of 8540 candidate groups of orthologous genes from transcriptomic data of five species of caecilian amphibians and the genome of the frog Xenopus tropicalis were analysed in order to investigate the genetic machinery behind caecilian diversification. We found a total of 168 protein-coding genes with signatures of positive selection at different evolutionary times during the radiation of caecilians. The majority of these genes were related to functional elements of the cell membrane and extracellular matrix with expression in several different tissues. The first colonization of the tropical soils was connected to the largest number of protein-coding genes under positive selection in our analysis. From the results of our study, we highlighted molecular changes in genes involved in perception, reduction-oxidation processes, and aging that likely were involved in the adaptation to different soil strata.

Conclusions

The genes inferred to have been under positive selection provide valuable insights into caecilian evolution, potentially underpin adaptations of caecilians to their extreme environments, and contribute to a better understanding of fossorial adaptations and molecular evolution in vertebrates.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5694-1) contains supplementary material, which is available to authorized users.