Research Article |
|
Corresponding author: Tuba Yağcı Gurbanov ( tuba.yagci@bilecik.edu.tr ) Academic editor: Peter Trontelj
© 2025 Tuba Yağcı Gurbanov, Rafig Gurbanov, Uygar Kabaoğlu, Burcu Tekin, Nursel Aşan Baydemir.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Gurbanov TY, Gurbanov R, Kabaoğlu U, Tekin B, Baydemir NA (2025) Contributions of molecular and cranial specializations for sensory adaptation in the subterranean blind mole rats (Nannospalax xanthodon). Subterranean Biology 52: 201-226. https://doi.org/10.3897/subtbiol.52.155412
|
The blind mole rats which have adapted excellently to subterranean life, exhibit remarkable navigation abilities despite structural regression in their visual systems and decline in auditory functions. This study investigates the sensory capabilities enabling blind mole rat possible adaptation to dark environments through gene expression analyses, evolutionary sequence comparisons, and cranial morphometric examinations. Comparative analyses reveal a significant reduction in the expression of genes related to vision (Pax6, Impg) and hearing (Coch, Necap1) in subterranean blind mole rat Nannospalax xanthodon compared to surface dwelling Rattus rattus. In contrast, genes associated with sound perception (Foxp2, Bmp7, Kcnq4, Tmc1) show a significant increase in expression. These findings indicate that sensory losses are aligned with the classical model observed in species that compensate for sensory deficits through echolocation. Remarkably, the gene expression profiles of N. xanthodon show similarities to those of mammals with highly developed echolocation abilities, revealing common molecular patterns between seismic and ultrasonic echolocation. Phylogenetic analyses of echolocation related genes show that blind mole rats share specific amino acid motifs with certain echolocating bat and whale species. Notably, positive selection signal observed in the Kcnq4 gene suggest that this gene may play a critical role in seismic vibration detection. Cranial morphometry studies support the notion that N. xanthodon possess not only putative morphological adaptations specific to their digging lifestyle but also specialized anatomical structures for detecting seismic vibrations. This study provides contributions to understanding the sensory adaptation of subterranean mammals and evaluates blind mole rats in this context for the first time.
Dark life, echolocation, gene expression, morphometry, phylogenetic analysis, seismic vibration
It has been suggested that the visual and auditory adaptations acquired by nocturnal mammals to facilitate activity in the dark gave them a survival advantage over diurnal lineages that re-emerged during the Cenozoic era (
In blind mole rats (Nannospalax ehrenbergi) the central visual pathways, which serve form and motion perception in sighted mammals, are taken over by the auditory system, and the occipital cortex, while structurally similar, is functionally activated by auditory stimuli (
Echolocation is a specialized communication method based on sound production and the reception of returning echoes, allowing animals to gather information about their environment. Although the mechanism by which Nannospalax detects underground vibrations is debated, their extraordinary echo-location system, where the sender and receiver are the same organism, has been described as a form of echolocation (
The use of different organs by bats and dolphins to produce and emit sound has not prevented the development of similar sound signal detection and processing systems (
Behavioral evidence suggesting echolocation abilities in the common shrew (Sorex araneus), a nocturnal insectivore with poor vision, and the soft-furred tree mouse (Typhlomys), a nearly blind semi-arboreal rodent, is molecularly supported by the detection of convergent amino acid changes in hearing genes associated with echolocation (
This study comprehensively examines the putative adaptation strategies developed by blind mole rats (Nannospalax xanthodon) for acoustic communication and navigation, despite the loss of their visual abilities. The research adopts an integrative approach that includes echolocation mechanisms to explore the impact of ecological niches on morphogenetic processes. The study focused on specific genes, which are involved in sound perception and signal processing, and are known to play key roles in the evolutionary modifications observed in echolocating species. Similarly, other genes associated with visual and auditory functions were also included in the study to better understand the sensory characteristics of blind mole rats potentially adapted to subterranean life at the molecular level. The comparative analyses conducted reveal how the distinct molecular and cranial differences between surface dwelling rats (Rattus rattus) and blind mole rats reflect ecological niche differentiation. The findings offer new research perspectives for understanding the possible selection pressures exerted by specialized habitats. The study also emphasizes the need for a deeper investigation into the probable underground adaptations, particularly focusing on the neurogenetic basis and its connection to cranial morphology and behavioral ecology. These findings provide an important reference point for sensory evolution and ecological adaptation studies that may be associated with subterranean life.
RNA isolation was performed using brain tissue samples (20 mg) from a total of 16 animals (8 individuals per species), representing the rodent group (Nannospalax xanthodon and Rattus rattus), which were stored at -80 °C. The samples belong to wild individuals of N. xanthodon with 2n = 60 chromosomes and R. rattus, collected from open fields in the Western Anatolia Region (Bilecik province). The isolation was carried out following the manufacturer’s instructions using the innuPREP RNA Mini Kit 2.0 (Analytik Jena GmbH). The collection of these specimens was conducted under permits issued by the General Directorate of Nature Conservation and National Parks, an agency of the Turkish Ministry of Agriculture and Forestry (Permission No: 72784983-488.04), and with approval from the Kırıkkale University Local Ethics Committee for Animal Research (Approval No: 15/02-15-18).
cDNA was synthesized from RNAs using the SCRIPT cDNA Synthesis Kit (Jena Bioscience) according to the manufacturer’s instructions. Random primers (0.5 µL from 100 µM stock) were incubated for 5 min at 65 °C for RNA (10 µL from 10 ng/µL stock) binding and then the reaction media were taken on the ice. Then the prepared reaction mixture (Suppl. material
Gene expression analysis was conducted using the LightCycler 480 II system (Roche Diagnostics) with innuMIX qPCR MasterMix SyGreen Sensitive (Analytik Jena GmbH) and 96-well DNase/RNase-free plates. A total of 8 genes were analyzed in this study, including hearing-related genes (Coch [Cochlin] and Necap1 [Adaptin ear-binding coat-associated protein 1]), vision-related genes (Impg2 [Interphotoreceptor Matrix Proteoglycan 2] and Pax6 [Paired box 6]), and genes associated with auditory perception (Kcnq4 [Potassium Voltage-Gated Channel Subfamily Q Member 4], Foxp2 [Forkhead Box P2], Tmc1 [Transmembrane Channel Like 1], and Bmp7 [Bone Morphogenetic Protein 7]) (
Primer design was performed using mRNA sequences from six species: Nannospalax galili, Rattus norvegicus, Mus musculus, Myotis blythii, Myotis myotis, and Miniopterus schreibersii. Conserved regions were identified using the Clustal Omega Multiple Sequence Alignment tool (https://www.ebi.ac.uk/Tools/msa/clustalo/). Conserved regions were then identified by visually examining the alignments to detect stretches of identical nucleotide sequences shared across all species. Suitable regions for primer design were selected based on the following criteria: a length of 19–25 base pairs and identical nucleotide sequences across all species. A total of 23 primer regions were chosen for primer control and Primer-BLAST analysis. Subsequently, the NCBI Primer-BLAST tool (https://www.ncbi.nlm.nih.gov/tools/primer-blast/) was used to screen primers against species from the Rodentia and Chiroptera orders. Factors such as off-target binding, primer dimer formation, and temperature ranges were evaluated to optimize primer suitability before the study. The reagents and gene primers (designed specifically for this study) used in the protocol are provided in Suppl. material
Quantity (relative) values were calculated with the help of the standard graph of Cq (ΔRn) values that changed during gene expression using the LightCycler 480 software. For each species, 8 animals with three replicates, (n = 24) were used in the analysis. The differences in targeted cDNA amount were eliminated by proportioning the targeted gene quantity and reference gene (GAPDH) values in GraphPad Prism 8.0 software. The concentration (ng) of gene expression was normalized with GAPDH (reference gene) expression. The data were analyzed by applying Dunnett’s multiple comparisons test with two-way ANOVA. The degrees of significance were denoted as p ≤ 0.05 *, p ≤ 0.01 **, and p ≤ 0.001 ***. The trade-off between data sensitivity and specificity were analysed using ROC (Receiver Operating Characteristic) test in GraphPad Prism 8.0 software.
Coding sequences of associated with auditory perception genes identified in the literature as key genes for echolocation, were collected from GenBank (https://www.ncbi.nlm.nih.gov/) and Ensembl Genome Database (https://www.ncbi.nlm.nih.gov/). Echolocating (marked with red) and non-echolocating species (marked with black) were selected based on previous studies, with reference genome for blind mole rats N. galili (marked with blue) included as the focal species (
To detect signals of positive selection, the codeml module of the PAML software was used. The analysis was carried out using nucleotide sequence alignments and a corresponding phylogenetic tree that included N. galili as well as echolocating and non-echolocating species (
The study analyzed the skulls of 28 adult male blind mole rats (N. xanthodon) and 18 adult male rats (R. rattus). Selected external cranial parameters associated with sound perception (braincase length, braincase height, greatest skull length, greatest skull height, rostrum length, tympanic bulla length, tympanic bulla width, zygomatic arch length) were measured to 0.01 mm precision using digital calipers. For Principal Component Analysis (PCA), data from 7 cranial parameters were used, with the values being normalized to the greatest skull length to determine the relationship between the skulls of blind mole rats and rat populations. All cranial data were transferred to Unscrambler X 10.3 (Camo, NO) multivariate analysis (MVA) software. The PCA model was generated using the Singular Value Decomposition (SVD) algorithm, and the results were presented as scores and loadings plots. Hotelling’s T² statistics with a 0.5% confidence limit were applied to identify potential outliers in the PCA scores plot, ensuring that the main variations captured by the model represent the underlying data structure rather than anomalous observations. The loadings plot allowed the analysis of all cranial parameters, as they show the positions that contribute the most to the variation explained by the PC terms. The magnitude of the difference between animals is represented by the loadings or absolute values at a specific position.
The expression levels of genes related to vision, such as Pax6 and Impg, hearing, such as Coch and Necap1, and sensory perception of sound (echolocation), such as Kcnq4, Foxp2, Tmc1, and Bmp7, were comparatively analyzed in subterranean blind mole rats (Nannospalax xanthodon) and surface-dwelling rats (Rattus rattus) to identify changes in gene expression associated with subterranean living conditions. Four echolocation-related genes showed significantly higher expression levels in mole rats (N) compared to rats (R), while genes related to hearing and vision exhibited lower expression levels. In surface dwelling rats, the expression of echolocation-related genes was significantly lower than that of vision and hearing genes (Figs
To investigate potential convergent evolution related to echolocation, we analyzed six key echolocation genes (Kcnq4, Prestin, Bmp7, Pjvk, Otof, and Cdh23) in Nannospalax galili as well as in echolocating and non-echolocating mammals (
Phylogenetic relationships and amino acid substitutions in auditory sensory perception-related genes among selected mammalian species, including echolocating lineages and the blind mole rat. Species capable of echolocation are highlighted in magenta, while the blind mole rat lineage is shown in green. For each species, amino acid substitutions in four key auditory sensory perception-related proteins (KCNQ4, Prestin, Cdh23, and Otof) are indicated. Substitutions observed in the blind mole rat , and those that are shared between the blind mole rat and at least one echolocating species are highlighted with a green background.(The phylogenetic tree shown in the figure represents an illustrative summary of the ML tree based on Kcnq4 gene sequences (see Suppl. material
Sequence alignment of the Kcnq4 gene across 25 species revealed that N. galili shares amino acid substitutions in three different regions with echolocating species (Suppl. material
PAML site model analyses were conducted to examine the molecular evolutionary patterns of six hearing-related genes (Kcnq4, Prestin, Otof, Cdh23, Bmp7, and Pjvk) across mammalian species (Suppl. material
The morphometric data were normalized to greatest skull length and analyzed using PCA to statistically evaluate cranial differences and thereby possible auditory adaptations between the species. Morphometric analyses have revealed that the cranial features of N. xanthodon, representing a subterranean lifestyle, are significantly distinct from those of R. rattus based on characters functionally linked to auditory perception. The most prominent morphological differences between the species, according to PCA analysis, were found in Zygomatic arch length (0.66), Braincase width (0.50), and Rostrum length (0.42) (Fig.
Comparative analysis of external cranial measurements between Nannospalax xanthodon (N) and Rattus rattus (R). Measured parameters: BW: Braincase length, BH: Braincase height, RL: Rostrum length, TBL: Tympanic bulla length, TBW: Tympanic bulla width, ZAL: Zygomatic arch length, GSH: Greatest skull height.
Subterranean life is not confined solely to caves and other natural underground habitats, such as interstitial spaces, but has expanded to include a variety of underground environments with diverse morphologies and microenvironmental conditions. These include shallow subterranean habitats located near the surface (less than 10 meters deep). The defining feature common to all of these environments is the absence of light (
We compared the expression levels of certain genes implicated in the development of vision, hearing, and echolocation between blind mole rats (N. xanthodon) and rats (R. rattus). Our findings revealed that the expression levels of Pax6 and Impg genes, associated with vision, were significantly lower in the subterranean blind mole rats compared to surface-dwelling rats, confirming the regression of light-dependent retinal genes (Fig.
In blind mole rats, degenerations such as weak auditory sensitivity, high-frequency hearing loss, and the inability to localize short sounds are characteristics specific to subterranean mammals. It has been suggested that adaptations to subterranean living conditions such as the poor propagation of airborne sounds and the restricted directional movement within tunnels may lead to the degeneration of auditory abilities, similar to how the absence of light diminishes visual capabilities (
On the other hand, the significantly higher expression of Foxp2, Bmp7, Kcnq4, and Tmc1 genes in blind mole rats compared to rats is intriguing, as these genes play a crucial role in the development of echolocation in echolocating species (Fig.
In adult mice, Bmp7, which is associated with the cochlea’s ability to repair and regenerate lost sensory hair cells, has been examined as a marker gene for inner ear cells. This capability has been attributed to the lower expression levels of this gene postnatally (
In mouse models exposed to high and low-frequency acoustic noise, loss of Kcnq4 in outer hair cells (OHCs) has been observed, revealing that Kcnq4 protects OHCs from noise-induced Ca2+ overload. Mutations and reduced activity in Kcnq4 are believed to underlie the common molecular basis of deafness and noise-induced hearing loss. Therefore, the modulation of Kcnq4 activity is considered an effective strategy for the treatment and prevention of these conditions (
Tmc1, which is essential for sensory transduction in auditory and vestibular hair cells, has been shown to have therapeutic effects on a variety of auditory disorders, including the behavioral modulation of auditory function in mice through gene therapy applications (
In some cases, it is thought that selection may act on both coding sequences and gene regulation to provide enhanced hearing in echolocating mammals (
Common adaptive strategies in response to similar ecological niches are widespread among mammals and represent significant evolutionary traits that ensure their survival. Mammals, such as the bottlenose dolphin (Tursiops truncatus), common shrew (Sorex araneus), subterranean naked mole rat (Heterocephalus glaber), and some bat species, exhibit convergent genomic changes, including losses of antioxidant genes (e.g., GPX6), as adaptations to oxidative stress environments (
Blind mole rats (Nannospalax ehrenbergi) use self-generated seismic waves as an echolocation mechanism to determine the size and shape of obstacles. Other subterranean species that produce seismic signals by striking their heads, hind legs, or incisors include the African mole rat (Georychus capensis), Damaraland mole rat (Fukomys damarensis), giant mole rat (Fukomys mechowii), and demon African mole rat (Tachyoryctes daemon) (
The cranial features of N. xanthodon, adapted to subterranean life, significantly differ from R. rattus in terms of auditory-related parameters. The most notable differences were observed in the measurements of zygomatic arch length (0.66), braincase width (0.50), and rostrum length (0.42) (Fig.
This research suggests that gene expressions related to sensory communication in blind mole-rats may play a role in seismic vibration detection, a distinct variant of echolocation possibly adapted to subterranean life. While their visual functions have diminished, the decline in auditory functions may indicate their potential adaptations to the unique acoustic environment underground. However, significant gaps remain in our understanding of the mechanisms underlying seismic echolocation and its evolutionary trajectory. Future studies should focus on expanding the scope of the investigation to include a broader range of subterranean species, which will allow for comparative analyses and the identification of shared or divergent genetic adaptations. Additionally, tissue-specific investigations are crucial for pinpointing the precise roles of candidate genes in sensory perception and navigation. Beyond the genes associated with ultrasonic echolocation, there is a need to explore novel genetic factors that may support the development of alternative echolocation strategies tailored to subterranean environments. Such research not only holds the potential to uncover fundamental insights into the molecular basis of sensory evolution but also offers broader implications for understanding convergent evolution and ecological adaptation that may be associated with subterranean life. The structure of the skull may show potential adaptations that enhance the effectiveness of bone conduction, which is necessary for this type of acoustic transmission. However, it is important to note that the same anatomical features are also heavily influenced by digging behavior. Therefore, more comprehensive and complex analyses supported by a greater number of subterranean species are required to disentangle these two effects. By elucidating the interplay between genetics, environment, and behavior, future studies could provide a more comprehensive framework for examining the diversity of navigation strategies in subterranean habitats. Ultimately, these efforts will offer deep insights into how life adapts to even the most challenging and specialized niches on Earth.
This research was partially supported by the Scientific Research Project Fund of Bilecik Şeyh Edebali University (2018-01.BŞEÜ.04-04).
The authors have declared that no competing interests exists.
Supplementary figures nad tables
Data type: docx
Supplementary data
Data type: xlsx