Stress, adaptation, and speciation in the evolution of the blind mole rat, Spalax, in Israel

Host diet shapes functionally differentiated gut microbiomes in sympatric speciation of blind mole rats in Upper Galilee, Israel

The gut microbiome is important for host nutrient metabolism and ecological adaptation. However, how the gut microbiome is affected by host phylogeny, ecology and diet during sympatric speciation remain unclear. Here, we compare and contrast the gut microbiome of two sympatric blind mole rat species and correlate them with their corresponding host phylogeny, ecology soil metagenomes, and diet to determine how these factors may influence their gut microbiome. Our results indicate that within the host microbiome there is no significant difference in community composition, but the functions between the two sympatric species populations vary significantly. No significant correlations were found between the gut microbiome differentiation and their corresponding ecological soil metagenomes and host phylogeny. Functional enrichment analysis suggests that the host diets may account for the functional divergence of the gut microbiome. Our results will help us understand how the gut microbiome changes with corresponding ecological dietary factors in sympatric speciation of blind subterranean mole rats.

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.

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.

Antioxidant response to severe hypoxia in Brandt's vole Lasiopodomys brandtii

The antioxidant defense system is essential for animals to cope with homeostasis disruption and overcome oxidative stress caused by adverse environmental conditions such as hypoxia. However, our understanding of how this system works in subterranean rodents remains limited. In this study, Brandt's vole Lasiopodomys brandtii was exposed to normoxia (21% O₂ ) or hypoxia (mild or severe hypoxia: 10% or 5% O₂ ) for 6 h. Changes in key enzymes of the classic enzymatic antioxidant system at both mRNA and enzyme activity levels, and tissue antioxidant levels of the low-molecular-weight antioxidant system were determined in brain, liver, and kidney. Transcript levels of the upstream regulator NF-E2-related factor 2 (Nrf2) were also measured. We found that the mRNA expression of Nrf2 and its downstream antioxidant enzyme genes in L. brandtii were relatively conserved in response to hypoxia in most tissues and genes tested, except in the liver. Hepatic Nrf2, Cu/Zn SOD, GPx1, and GPx3 levels were significantly upregulated in response to mild hypoxia, whereas Mn SOD level decreased significantly in severe hypoxia. Unmatched with changes at the RNA level, constitutively high and relatively stable antioxidant enzyme activities were maintained throughout. For the low-molecular-weight antioxidant system, an abrupt increase of cerebral ascorbic acid (AA) levels in hypoxia indicated a tissue-specific antioxidant response. Although hypoxia did not cause significant oxidative damage in most tissues tested, the significant decrease in antioxidant enzyme activities (GPX and GR) and increase in lipid peroxidation in the kidney suggest that prolonged hypoxia may pose a critical threat to this species.

Phylogeny and evolution of Lasiopodomys in subfamily Arvivolinae based on mitochondrial genomics

The species of Lasiopodomys Lataste 1887 with their related genera remains undetermined owing to inconsistent morphological characteristics and molecular phylogeny. To investigate the phylogenetic relationship and speciation among species of the genus Lasiopodomys, we sequenced and annotated the whole mitochondrial genomes of three individual species, namely Lasiopodomys brandtii Radde 1861, L. mandarinus Milne-Edwards 1871, and Neodon (Lasiopodomys) fuscus Büchner 1889. The nucleotide sequences of the circular mitogenomes were identical for each individual species of L. brandtii, L. mandarinus, and N. fuscus. Each species contained 13 protein-coding genes (PCGs), 22 transfer RNAs, and 2 ribosomal RNAs, with mitochondrial genome lengths of 16,557 bp, 16,562 bp, and 16,324 bp, respectively. The mitogenomes and PCGs showed positive AT skew and negative GC skew. Mitogenomic phylogenetic analyses suggested that L. brandtii, L. mandarinus, and L. gregalis Pallas 1779 belong to the genus Lasiopodomys, whereas N. fuscus belongs to the genus Neodon grouped with N. irene. Lasiopodomys showed the closest relationship with Microtus fortis Büchner 1889 and M. kikuchii Kuroda 1920, which are considered as the paraphyletic species of genera Microtus. T_MRCA and niche model analysis revealed that Lasiopodomys may have first appeared during the early Pleistocene epoch. Further, L. gregalis separated from others over 1.53 million years ago (Ma) and then diverged into L. brandtii and L. mandarinus 0.76 Ma. The relative contribution of climatic fluctuations to speciation and selection in this group requires further research.

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.

A Chromosome-Level Genome Assembly of the Anglerfish Lophius litulon

Anglerfishes are a highly diverse group of species with unique characteristics. Here, we report the first chromosome-level genome of a species in the order Lophiiformes, the yellow goosefish (Lophius litulon), obtained by whole genome shotgun sequencing and high-throughput chromatin conformation capture. Approximately 97.20% of the assembly spanning 709.23 Mb could be anchored to 23 chromosomes with a contig N50 of 164.91 kb. The BUSCO value was 95.4%, suggesting that the quality of the assembly was high. A comparative gene family analysis identified expanded and contracted gene families, and these may be associated with adaptation to the benthic environment and the lack of scales in the species. A majority of positively selected genes were related to metabolic processes, suggesting that digestive and metabolic system evolution expanded the diversity of yellow goosefish prey. Our study provides a valuable genetic resource for understanding the mechanisms underlying the unique features of the yellow goosefish and for investigating anglerfish evolution.

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.

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.

Altitudinal Effects on Innate Immune Response of a Subterranean Rodent

Immune defense is costly to maintain and deploy, and the optimal investment into immune defense depends on risk of infection. Altitude is a natural environmental factor that is predicted to affect parasite abundance, with lower parasite abundance predicted at higher altitudes due to stronger environmental stressors, which reduce parasite transmission. Using high and low altitude populations of the Turkish blind mole-rat (TBMR) Nannospalax xanthodon, we tested for effects of altitude on constitutive innate immune defense. Field studies were performed with 32 wild animals in 2017 and 2018 from two low- and one high-altitude localities in the Central Taurus Mountains, at respective altitudes of 1010 m, 1115 m, and 2900 m above sea level. We first compared innate standing immune defense as measured by the bacteria-killing ability of blood serum. We then measured corticosterone stress hormone levels, as stressful conditions may affect immune response. Finally, we compared prevalence and intensity of gastrointestinal parasites of field-captured TBMR. We found that the bacteria-killing ability of serum is greater in the mole-rat samples from high altitude. There was no significant difference in stress (corticosterone) levels between altitude categories. Coccidian prevalence and abundance were significantly higher in 2017 than 2018 samples, but there was no significant difference in prevalence, abundance, or intensity between altitudes, or between sexes. Small sample sizes may have reduced power to detect true differences; nevertheless, this study provides support that greater standing innate immunity in high altitude animals may reflect greater investment into constitutive defense.

Differential responses of Lasiopodomys mandarinus and Lasiopodomys brandtii to chronic hypoxia: a cross-species brain transcriptome analysis

BACKGROUND

Subterranean rodents have evolved many features to adapt to their hypoxic environment. The brain is an organ that is particularly vulnerable to damage caused by exposure to hypoxic conditions. To investigate the mechanisms of adaption to a hypoxic underground environment, we carried out a cross-species brain transcriptome analysis by RNA sequencing and identified genes that are differentially expressed between the subterranean vole Lasiopodomys mandarinus and the closely related above-ground species Lasiopodomys brandtii under chronic hypoxia [10.0% oxygen (O2)] and normoxia (20.9% O2).

RESULTS

A total of 355 million clean reads were obtained, including 69,611 unigenes in L. mandarinus and 69,360 in L. brandtii. A total of 235 and 92 differentially expressed genes (DEGs) were identified by comparing the hypoxic and control groups of L. mandarinus and L. brandtii, respectively. A Gene Ontology (GO) analysis showed that upregulated DEGs in both species had similar functions in response to hypoxia, whereas downregulated DEGs in L. mandarinus were enriched GO terms related to enzymes involved in aerobic reactions. In the Kyoto Encyclopedia of Genes and Genomes pathway analysis, upregulated DEGs in L. mandarinus were associated with angiogenesis and the increased O2 transport capacity of red blood cells, whereas downregulated DEGs were associated with immune responses. On the other hand, upregulated DEGs in L. brandtii were associated with cell survival, vascular endothelial cell proliferation, and neuroprotection, while downregulated genes were related to the synaptic transmission by neurons.

CONCLUSIONS

L. mandarinus actively adapts its physiological functions to hypoxic conditions, for instance by increasing O2 transport capacity and modulating O2 consumption. In contrast, L. brandtii reacts passively to hypoxia by decreasing overall activity in order to reduce O2 consumption. These results provide insight into hypoxia adaptation mechanisms in subterranean rodents that may be applicable to humans living at high altitudes or operating in other O2-poor environments.

Resistance to DNA damage and enhanced DNA repair capacity in the hypoxia-tolerant blind mole rat Spalax carmeli

Blind mole rats of the genus Spalax are the only mammalian species to date for which spontaneous cancer has never been reported and resistance to carcinogen-induced cancers has been demonstrated. However, the underlying mechanisms are still poorly understood. The fact that Spalax spp. are also hypoxia-tolerant and long-lived species implies the presence of molecular adaptations to prevent genomic instability, which underlies both cancer and aging. We previously demonstrated the upregulation of transcripts related to DNA replication and repair pathways in Spalax Yet, to date, no direct experimental evidence for improved genomic maintenance has been demonstrated for this genus. Here, we show that compared with skin fibroblasts of the above-ground rat, Spalax carmeli skin fibroblasts in culture resist several types of genotoxic insult, accumulate fewer genotoxic lesions and exhibit an enhanced DNA repair capacity. Our results strongly support that this species has evolved efficient mechanisms to maintain DNA integrity as an adaptation to the stressful conditions in the subterranean habitat.

Adaptive Transcriptome Profiling of Subterranean Zokor, Myospalax baileyi, to High- Altitude Stresses in Tibet

Animals living at high altitudes have evolved distinct phenotypic and genotypic adaptations against stressful environments. We studied the adaptive patterns of altitudinal stresses on transcriptome turnover in subterranean plateau zokors (Myospalax baileyi) in the high-altitude Qinghai-Tibetan Plateau. Transcriptomes of zokors from three populations with distinct altitudes and ecologies (Low: 2846 m, Middle: 3282 m, High: 3,714 m) were sequenced and compared. Phylogenetic and principal component analyses classified them into three divergent altitudinal population clusters. Genetic polymorphisms showed that the population at H, approaching the uppermost species boundary, harbors the highest genetic polymorphism. Moreover, 1056 highly up-regulated UniGenes were identified from M to H. Gene ontologies reveal genes like EPAS1 and COX1 were overexpressed under hypoxia conditions. EPAS1, EGLN1, and COX1 were convergent in high-altitude adaptation against stresses in other species. The fixation indices (F_ST and G_ST)-based outlier analysis identified 191 and 211 genes, highly differentiated among L, M, and H. We observed adaptive transcriptome changes in Myospalax baileyi, across a few hundred meters, near the uppermost species boundary, regardless of their relatively stable underground burrows’ microclimate. The highly variant genes identified in Myospalax were involved in hypoxia tolerance, hypercapnia tolerance, ATP-pathway energetics, and temperature changes.

Exceptional Chromosomal Evolution and Cryptic Speciation of Blind Mole Rats Nannospalax leucodon (Spalacinae, Rodentia) from South-Eastern Europe

Mole rats are exclusively subterranean and highly specialized rodents. Their long lifespans, remarkable anti-cancer mechanisms, and various distinctive adaptive features make them a useful research model. Moreover, opposing convergence of morphological traits, they have developed extremely high karyotype variability. Thus, 74 chromosomal forms have been described so far and new ones are being revealed continuously. These evolved during the process of rapid radiation and occur in different biogeographical regions. During research into their reproductive biology we have already provided substantial evidence for species-level separation of these taxa. Here, we review diverse chromosomal forms of the lesser blind mole rat, Mediterranean Nannospalax leucodon, distributed in South-eastern Europe, their karyotype records, biogeography, origin, and phylogeny from our extensive research. In the light of new data from molecular genetic studies, we question some former valuations and propose a cryptospecies rank for seven reproductively isolated chromosomal forms with sympatric and parapatric distribution and clear ecogeographical discrepances in their habitats, as well as new experimental and theoretical methods for understanding the courses of speciation of these unique fossorial mammals.

Determinants of rodent longevity in the chaperone-protein degradation network

Proteostasis is an integral component of healthy aging, ensuring maintenance of protein structural and functional integrity with concomitant impact upon health span and longevity. In most metazoans, increasing age is accompanied by a decline in protein quality control resulting in the accrual of damaged, self-aggregating cytotoxic proteins. A notable exception to this trend is observed in the longest-lived rodent, the naked mole-rat (NMR, Heterocephalus glaber) which maintains proteostasis and proteasome-mediated degradation and autophagy during aging. We hypothesized that high levels of the proteolytic degradation may enable better maintenance of proteostasis during aging contributing to enhanced species maximum lifespan potential (MLSP). We test this by examining proteasome activity, proteasome-related HSPs, the heat-shock factor 1 (HSF1) transcription factor, and several markers of autophagy in the liver and quadriceps muscles of eight rodent species with divergent MLSP. All subterranean-dwelling species had higher levels of proteasome activity and autophagy, possibly linked to having to dig in soils rich in heavy metals and where underground atmospheres have reduced oxygen availability. Even after correcting for phylogenetic relatedness, a significant (p < 0.02) positive correlation between MLSP, HSP25, HSF1, proteasome activity, and autophagy-related protein 12 (ATG12) was observed, suggesting that the proteolytic degradation machinery and maintenance of protein quality play a pivotal role in species longevity among rodents.

The captured retroviral envelope syncytin-A and syncytin-B genes are conserved in the Spalacidae together with hemotrichorial placentation

Syncytins are fusogenic envelope (env) genes of retroviral origin that have been captured for a function in placentation. Multiple independent events of syncytin gene capture were found to have occurred in primates, rodents, lagomorphs, carnivores, and ruminants. In the mouse, two syncytin-A and -B genes are present, which trigger the formation of the two-layered placental syncytiotrophoblast at the maternal-fetal interface, a structure classified as hemotrichorial. Here, we identified syncytin-A and -B orthologous genes in the genome of all Muroidea species analyzed, thus dating their capture back to about at least 40 million years ago, with evidence that they evolved under strong purifying selection. We further show, in the divergent Spalacidae lineage (blind mole rats [Spalax]), that both syncytins have conserved placenta-specific expression, as revealed by RT-PCR analysis of a panel of Spalax galili tissues, and display fusogenic activity, using ex vivo cell-cell fusion assays. Refined analysis of the placental architecture and ultrastructure revealed that the Spalax placenta displays a hemotrichorial organization of the interhemal membranes, as similarly observed for other Muroidea species, yet with only one trophoblastic cell layer being clearly syncytialized. In situ hybridization experiments further localized syncytin transcripts at the level of these differentiated interhemal membranes. These findings argue for a role of syncytin gene capture in the establishment of the original hemotrichorial placenta of Muroidea, and more generally in the diversity of placental structures among mammals.

Local adaptation and the evolution of chromosome fusions

We use forward and coalescent models of population genetics to study chromosome fusions that reduce the recombination between two locally adapted loci. Under a continent-island model, a fusion spreads and reaches a polymorphic equilibrium when it causes recombination between locally adapted alleles to be less than their selective advantage. In contrast, fusions in a two-deme model always spread; whether it reaches a polymorphic equilibrium or becomes fixed depends on the relative recombination rates of fused homozygotes and heterozygotes. Neutral divergence around fusion polymorphisms is markedly increased, showing peaks at the point of fusion and at the locally adapted loci. Local adaptation could explain the evolution of many of chromosome fusions, which are some of the most common chromosome rearrangements in nature.

Possible incipient sympatric ecological speciation in blind mole rats (Spalax)

Sympatric speciation has been controversial since it was first proposed as a mode of speciation. Subterranean blind mole rats (Spalacidae) are considered to speciate allopatrically or peripatrically. Here, we report a possible incipient sympatric adaptive ecological speciation in Spalax galili (2n = 52). The study microsite (0.04 km(2)) is sharply subdivided geologically, edaphically, and ecologically into abutting barrier-free ecologies divergent in rock, soil, and vegetation types. The Pleistocene Alma basalt abuts the Cretaceous Senonian Kerem Ben Zimra chalk. Only 28% of 112 plant species were shared between the soils. We examined mitochondrial DNA in the control region and ATP6 in 28 mole rats from basalt and in 14 from chalk habitats. We also sequenced the complete mtDNA (16,423 bp) of four animals, two from each soil type. Remarkably, the frequency of all major haplotype clusters (HC) was highly soil-biased. HCI and HCII are chalk biased. HC-III was abundant in basalt (36%) but absent in chalk; HC-IV was prevalent in basalt (46.5%) but was low (20%) in chalk. Up to 40% of the mtDNA diversity was edaphically dependent, suggesting constrained gene flow. We identified a homologous recombinant mtDNA in the basalt/chalk studied area. Phenotypically significant divergences differentiate the two populations, inhabiting different soils, in adaptive oxygen consumption and in the amount of outside-nest activity. This identification of a possible incipient sympatric adaptive ecological speciation caused by natural selection indirectly refutes the allopatric alternative. Sympatric ecological speciation may be more prevalent in nature because of abundant and sharply abutting divergent ecologies.