Group living in highland tuco-tucos (Ctenomys opimus) persists despite a catastrophic decline in population density.

Identifying the factors that favor group living is central to studies of animal social behavior. One demographic parameter that is expected to substantially shape spatial and social relationships is population density. Specifically, high population densities may favor group living by constraining opportunities to live alone. In contrast, low densities may allow individuals to spread out within the habitat, leading to a reduction in the prevalence or size of social groups. Abrupt changes in density following natural catastrophic events provide important opportunities to evaluate the effects of population density on patterns of spatial and social organization. As part of long-term studies of the behavioral ecology of a population of highland tuco-tucos (Ctenomys opimus) at Monumento Natural Laguna de los Pozuelos, Jujuy Province, Argentina, we monitored the demographic and behavioral consequences of a flood that inundated our study site during December 2012. Unlike most species of Ctenomys studied to date, highland tuco-tucos are group living, meaning that multiple adults share burrow systems and nest sites. Despite a post-flood reduction in population density of ~75%, animals present on the study site during the 2013 breeding season continued to live in multi-adult social units (groups). No differences between pre- and post-flood home range sizes were detected and although between-unit spatial overlap was reduced in 2013, overlap within social units did not differ from that in pre-flood years. Animals assigned to the same social unit in 2013 had not lived together during 2012, indicating that post-flood groups were not simply the remnants of those present prior to the flood. Collectively, these findings indicate that group living in highland tuco-tucos is not driven by the density of conspecifics in the habitat. In addition to enhancing understanding of the adaptive bases for group living in Ctenomys, our analyses underscore the power of catastrophic events to generate insights into fundamental aspects of social behavior.

Effect of estrogens on apoptosis in the pituitary of viscachas (Lagostomus maximus maximus).

Estrogens, acting through their receptors (ERα and ERß), regulate cell turnover in the pituitary gland, influencing cell proliferation and apoptosis across various species. However, their role in pituitary processes in seasonally reproducing animals remains poorly understood. This study aims to investigate the influence of estrogens, through the expression of their specific receptors, on the apoptosis of PD cells in relation to sexual maturity, the reproductive cycle, and pregnancy in a seasonal reproductive rodent (Lagostomus maximus maximus). ERα and caspase-3-cleaved (CASP3c) immunoreactive (-ir) cells were identified through immunohistochemistry. Apoptotic cells were detected using the TUNEL technique, with quantitative analysis facilitated by image analysis software, alongside measurement of serum estradiol levels using radioimmunoassay The immunostaining pattern for ERα included nuclear (ERαn) and cytoplasmic (ERαc) staining. In male viscachas, ERα expression significantly increases from immature to adult animals, correlating with the rise in serum estradiol levels and a decrease in the percentage of apoptotic cells. During the gonadal regression period in adult males, a decrease in the number of ER-ir cells and serum levels of estradiol corresponds with an increase in the number of apoptotic cells. In females, serum levels of estradiol peaked during mid-pregnancy, coinciding with a significant decrease in the number of apoptotic cells in the PD. Simultaneously, the percentage of ERαn-ir cells reaches its maximum value during late pregnancy, indicating the need to maintain the protective action of this gonadal hormone throughout the extensive pregnancy in these rodents. Regional ERα receptor expression and apoptotic cells appear to be associated with distinct PD cell populations and their hormonal responses. Finally, elevated estradiol levels coincide with diminished apoptotic cells in the male reproductive cycle and during pregnancy, suggesting an antiapoptotic role of estradiol in this species.

The Evolution of Rodent Tail Morphology.

AbstractPopulation-level variation in rodent tail structures has been variously attributed to facilitating social communication, locomotion, thermoregulation, and predator avoidance. Little is known, however, about the applicability of these ecological and social correlates to explaining the tremendous interspecific diversity of this appendage. To investigate the potential drivers of rodent tail morphology at a macroevolutionary level, we first carefully reviewed the literature and constructed a list of major hypotheses regarding this variation. We then compiled a database of 11 different tail traits related to length, color, texture, and ecological characteristics for 2,101 species of rodents (order Rodentia) and examined their key evolutionary correlates. Using Bayesian phylogenetic mixed models across the entire order and additionally within the five rodent suborders, we found that tail length is correlated with both temperature (Allen's rule) and locomotory mode, that black tips are more common in brightly lit environments, that naked tails are often found in warmer climates, that fluffy-tipped tails are more common in smaller and/or arboreal species, that prehensility is predominant in arboreal species and/or species with longer tails, and that tail autotomy is more common in open environments. Most of our tested predictions, largely drawn from population-level studies, are not recapitulated across the entire order, potentially indicating a role of local ecological context in shaping tail morphology.

Masting promotes transformation from predation to mutualism in an oak-weevil-rodent system.

The significance of ecological non-monotonicity (a function whose first derivative changes signs) in shaping the structure and functions of the ecosystem has recently been recognized, but such studies involving high-order interactions are rare. Here, we have proposed a three-trophic conceptual diagram on interactions among trees, rodents, and insects in mast and non-mast years and tested the hypothesis that oak (Quercus wutaishanica) masting could result in increased mutualism and less predation in an oak-weevil-rodent system in a warm temperate forest of China. Our 14-year dataset revealed that mast years coincided with a relatively low rodent abundance but a high weevil abundance. Masting not only benefited seedling recruitment of oaks through increased dispersal by rodents but also a decrease in predation by rodents and weevils, as well as an increase in the overwintering survival of rodents. Masting appeared to have increased weevil survival by reducing predation of infested acorns by rodents. These results suggest that masting benefits all participants in the plant-insect-rodent system by increasing mutualism and reducing predation behavior (i.e., a non-monotonic function). Our study highlights the significance of masting in maintaining the diversity and function of the forest ecosystem by facilitating the transformation from predation to mutualism among trophic species.

Is membrane androgen and estrogen receptor signaling imperative in the governing function of the adrenal cortex in the Eurasian beaver (Castor fiber L.)?

There is a need to fully know the physiology of Eurasian beaver due to its essential role in environmental homeostasis. However, a "human factor" impacts this, including stress conditions and environmental pollution. Adrenal glands protect these all. The regulation of endocrine processes by nonclassical androgen and estrogen signaling, the first and fastest control, is still a matter of research. The specific analyses performed here in mature female and male beaver adrenals contained: anatomical and histological examinations, expression and localization of membrane androgen receptor (zinc transporter, Zinc- and Iron-like protein 9; ZIP9) and membrane estrogen receptor coupled with G protein (GPER), and measurement of zinc (Zn2+) and copper (Ca2+) ion levels and corticosterone levels. We revealed normal anatomical localization, size, and tissue histology in female and male beavers, respectively. Equally, ZIP9 and GPER were localized in the membrane of all adrenal cortex cells. The protein expression of these receptors was higher (p < 0.001) in male than female adrenal cortex cells. Similarly, Zn2+ and Ca2+ ion levels were higher (p < 0.05, p < 0.01) in male than female adrenal cortex. The increased corticosterone levels (p < 0.001) were detected in the adrenal cortex of females when compared to males. The present study is the first to report the presence of nonclassical androgen and estrogen signaling and its possible regulatory function in the adrenal cortex of Eurasian beavers. We assume that this first-activated and fast-transmitted regulation can be important in the context of the effect of environmental physical and chemical stressors especially on adrenal cortex cells. The beaver adrenals may constitute an additional supplementary model for searching for universal mechanisms of adrenal cortex physiology and diseases.

Short-term dynamics of beaver dam flow states.

Beavers (Castor canadensis and C. fiber) build dams that modify catchment and pond water balances, and it has been suggested that they can be a nature-based solution for reducing flood hydrographs, enhancing low flow hydrographs and restoring hydrological functioning of degraded streams. How water moves through a beaver dam is determined by its flow state (e.g., overflow, underflow). However, current conceptual models only consider flow state as changing over the beaver site occupation-abandonment cycle. To assess whether flow state changes at shorter timescales and identify possible triggers (e.g., rainfall, animals), we integrated camera trap imagery, machine learning, water level measurements, and hydrometeorological data at beaver dams in a montane peatland in the Canadian Rocky Mountains. Contrary to current models, we found that flow states changed frequently, changing a maximum 12 times during the 139-day study period, but that changes had limited synchronicity amongst the dams in the same stream. More than two-thirds of the changes coincided with rainfall events. We observed no changes in flow state in response to beaver activity or wildlife crossings perhaps due to the camera positioning. Our findings augment the long-term oriented framework, which links changes to the occupancy cycle of a beaver pond and frequent and hydrological-driven changes. To develop realistic predictions of hydrological impacts of beaver dams, ecohydrological models should update their representation of the influence of beaver dams to include short-term dynamism of flow states and potential triggers. Our study advances the understanding of the important, yet understudied, role of beaver dams in stream restoration and climate change initiatives.

Effects of acute maternal stress induced by predator cues on spatial learning and memory of offspring in the subterranean rodent Ctenomys talarum.

One of the main selection pressures to which animals are exposed in nature is predation, which affects a wide variety of biological traits. When the mother experiences this stressor during pregnancy and/or lactation, behavioral and physiological responses may be triggered in the offspring as well. Thus, in order to broaden and deepen knowledge on the transgenerational effects of predation stress, we evaluated how maternal stress experienced during pregnancy and/or lactation affects the spatial abilities of progeny at the onset of adulthood in the subterranean rodent Ctenomys talarum. The results showed that, contrary to what was observed in other rodent species, maternal exposure to predator cues during pregnancy and lactation did not negatively affect the spatial abilities of the offspring, even registering some minor positive effects. Concomitantly, no effects of predatory cues on physiological parameters associated with stress were observed in the progeny. This difference in results between the present study and previous works on maternal stress highlights the importance of considering the species to be evaluated (strain, age and origin-wild or captive-) and the type of stressor used (artificial or natural, intensity of exposure) in the evaluation of the possible transgenerational effects of maternal stress.

Embryonic-placental relationship in Lagostomus maximus as compared to other hystricognath rodents and eutherian mammals.

Reproductive specializations in caviomorphs (infraorder Hystricognathi), are very peculiar within the Order Rodentia. These include long gestations, the birth of offspring with an extreme degree of precociality, and short lactation periods. This study describes the embryo-placental relationship of viable implantation sites (IS) of the plains viscacha, Lagostomus maximus, after 46 post-coital days. The observations recorded in this study are comparatively discussed with those of other hystricognaths and eutherians. At this stage, the embryo resembles that of other eutherians. At this time of embryo development, the placenta exhibits a size, shape, and organization similar to that it will have in its mature state. Besides, the subplacenta is already highly folded. These characteristics are adequate to sustain the development of future precocial offspring. The mesoplacenta, a structure present in other hystricognaths and related to uterine regeneration is described for the first time in this species. This detailed description of the placental and embryonic structure contributes to the knowledge of the reproductive and developmental biology of the viscacha, as well as that of hystricognaths. These characteristics will allow testing other hypotheses related to the morphology and physiology of the placenta and subplacenta, and their relationship with the growth and development of precocial offspring in Hystricognathi.

A nonparametric test of group distributional differences for hierarchically clustered functional data.

Biological sex and gender are critical variables in biomedical research, but are complicated by the presence of sex-specific natural hormone cycles, such as the estrous cycle in female rodents, typically divided into phases. A common feature of these cycles are fluctuating hormone levels that induce sex differences in many behaviors controlled by the electrophysiology of neurons, such as neuronal membrane potential in response to electrical stimulus, typically summarized using a priori defined metrics. In this paper, we propose a method to test for differences in the electrophysiological properties across estrous cycle phase without first defining a metric of interest. We do this by modeling membrane potential data in the frequency domain as realizations of a bivariate process, also depending on the electrical stimulus, by adopting existing methods for longitudinal functional data. We are then able to extract the main features of the bivariate signals through a set of basis function coefficients. We use these coefficients for testing, adapting methods for multivariate data to account for an induced hierarchical structure that is a product of the experimental design. We illustrate the performance of the proposed approach in simulations and then apply the method to experimental data.

Effects of water restriction and dirt on grooming behavior in neotropical rodents (Trinomys setosus and T. yonenagae) (Echimyidae).

Grooming in rodents presents an evolutionarily conserved behavioral pattern that may cause water loss since saliva is deposited during washing onto large body surfaces. Trinomys yonenagae and Trinomys setosus are sister species of spiny rats occurring in Brazil, the former inhabiting a paleodesert of fixed dunes in the Caatinga, the latter being found in mesic environments of the Atlantic Forest. Consequently, it is expected that both species evolved under different selective pressures related to water balance, with T. yonenagae presenting mechanisms for dealing with water deprivation not found in T. setosus. Reduction of self-cleaning expression seems to offer a possible way to save water, as previously suggested by studies of the sand-dwelling spiny rat. Therefore, we propose to investigate grooming under four conditions: 'control' (C), a regimen of 'water restriction' (WR), of 'dirt' (D), and the combination of both conflicting stimuli (WR + D), in T. setosus, T. yonenagae, and Rattus norvegicus to compare the behavioral responses of these species. The main differences are observed in the forest dweller: T. setosus expresses a low relative duration of face washing under C, whose value is intermediate between the ones found in the two other species. WR treatment does not alter this pattern, however, the addition of dirt (D, WR + D) significantly increases the relative duration of washing in relation to C. Locomotor activity is decreased both in T. setosus and Wistar rats when they are under WR, a situation that could jeopardize antipredatory performance. T. yonenagae, the sand dweller, maintains a significantly lower expression of washing under C, as previously suggested, and under WR, D and WR + D. In addition, differently from the other two species the sand dweller maintains a normal activity level during all treatments. This study suggests differences in grooming as a strategy alluding to water balance by the two spiny rats inhabiting different ecosystems. A significantly clear pattern that saves water is observed in T. yonenagae, which probably has contributed to his evolution in one of the hottest semiarid areas of the world.

How beavers affect riverine aquatic macroinvertebrates: a review.

Background: As ecosystem engineers, the construction of dams by beavers alters stream habitat physically and biologically, making them a species of interest for habitat restoration. Beaver-created habitat changes affect a wide range of aquatic invertebrate species. However, despite numerous individual studies of how beavers affect aquatic macroinvertebrate assemblages, there has been no evaluation of the consensus of these effects across studies. Methodology: We collated and examined studies comparing beaver-created ponds to nearby lotic reaches to determine general trends in aquatic macroinvertebrate richness, density, biomass, and functional composition between habitats. From this evidence, we highight knowledge gaps in how beaver activity affects aquatic macroinvertebrates. Results: Overall, in the majority of studies, aquatic macroinvertebrate richness was higher in nearby lotic reaches compared to beaver-created ponds, but richness at coarser scales (gamma diversity) increased with the addition of beaver ponds due to increased habitat heterogeneity. Functional feeding group (FFG) patterns were highly context-dependent, though predator taxa were generally more abundant in beaver ponds than adjacent lotic reaches. Site-specific geomorphological changes, coupled with dam or riparian zone characteristics and resulting differences in basal food resources likely shape other FFG responses. Conclusions: We identify a lack of long-term studies at single or multiple sites and conclude that fine-scale approaches may improve our understanding of the dynamics of macroinvertebrates within the freshwater realm and beyond. Due to the context-dependent nature of each study, further systematic studies of beaver engineering effects across a wider variety of environmental conditions and wetland types will also help inform land and species management decisions, such as where to prioritize protection of beaver habitats in the face of a global freshwater biodiversity crisis, or where to restore beaver populations to deliver maximum benefit.

Rodent responses to volatile compounds provide insights into the function of floral scent in mammal-pollinated plants.

Flowers pollinated by mammals have evolved in many plant families. Several scent compounds that attract bats to flowers have been identified, but the chemical ecology of pollination mutualisms between plants and ground-dwelling mammals is poorly understood. Rodents are key pollinators in South Africa and rely heavily on olfaction to locate food. Our aim was to identify compounds that may function to attract rodents to flowers. Eighteen volatile compounds, including 14 that are prominent in the scent of rodent-pollinated flowers, were used in choice experiments involving wild-caught individuals of four native rodent species. Rodents were generally attracted to oxygenated aliphatic compounds, specifically ketones and esters, but not to some aromatic compounds common in floral scents of insect-pollinated species, nor to a sulfide compound that is attractive to bats. Associative conditioning using sugar solution as a reward had only weak effects on the attractiveness of compounds to rodents. The attractive effect of some compounds disappeared when they were blended with compounds that did not attract rodents. We conclude that aliphatic ketones and esters are likely to play a key role in attracting rodents to flowers. Deployment of these compounds may allow plants to exploit rodent sensory bias that evolved in other contexts such as intra-specific communication and searching for seeds. This article is part of the theme issue 'Natural processes influencing pollinator health: from chemistry to landscapes'.

Sex, not social behavior, predicts fecal glucocorticoid metabolite concentrations in a facultatively social rodent, the highland tuco-tuco (Ctenomys opimus).

Social relationships may influence circulating glucocorticoid levels, particularly in group-living species in which individuals regularly engage in interactions with conspecifics. The effects of such interactions appear to vary, with greater social contact being associated with increased glucocorticoid concentrations in some species but decreased concentrations in others. These distinct responses raise intriguing questions regarding relationships among social behavior, individual phenotypes, and glucocorticoid physiology. To explore such relationships in a free-living mammal with a dynamic social organization, we quantified variation in baseline glucocorticoids in a population of highland tuco-tucos (Ctenomys opimus) from Jujuy Province, Argentina. These subterranean rodents are facultatively social, with lone and group-living individuals regularly occurring within the same population. To assess potential endocrine correlates of this behavioral variability, we examined differences in baseline fecal glucocorticoid metabolite (fGCm) concentrations as a function of social group size and composition as well as several metrics of social behavior derived from social network analyses. Despite marked variability in social relationships among the 37 (12 male, 25 female) free-living tuco-tucos sampled, none of the measures of social behavior examined were significant predictors of variation in fGCm concentrations. In contrast, individual variation in glucocorticoid metabolites was best explained by sex, with males having higher fGCm concentrations than females. These analyses provide the first characterization of the glucocorticoid physiology of highland tuco-tucos and underscore the potential importance of intrinsic phenotypic factors (e.g., sex) in shaping glucocorticoid variation in free-living mammals.

Facilitative interaction promotes occupancy of a desert amphibian across a climate gradient.

Biotic interactions and environmental factors influence a species' occurrence. Facilitative interactions have the potential to expand species occupancy and mitigate abiotic stress, but are often not considered. The Great Basin clade of the Columbia spotted frog (Rana luteiventris) is an aquatic amphibian found in the arid and semi-arid Great Basin, USA. Aquatic resources in this system are limited and likely to change under future climate scenarios. American beavers (Castor canadensis) are ecosystem engineers that impact the surrounding landscape by increasing water availability. Consequently, beaver-created habitat can facilitate the presence of other species by providing the aquatic resources required by many aquatic and riparian animals. Our objective was to understand patterns of co-occurrence between Columbia spotted frogs and beavers across environmental gradients in the Great Basin, USA. We used environmental DNA detections to quantify the co-occurrence of Columbia spotted frog and beaver using a two-species occupancy framework. At the lowest annual precipitation (202 mm), Columbia spotted frogs were 3× (95% CI 2.35, 3.96) more likely to occur with beavers than without; however, they were less likely to occur with beavers once precipitation was > 380 mm. Thus, beaver activity may increase aquatic resources for Columbia spotted frogs in the Great Basin that may mitigate changing precipitation patterns under climate change. Facilitative interactions likely extend to other aquatic species in arid systems and highlight an important consideration for ecologists when evaluating a species' response to climate change, and may promote the formation of refugia for species with strict abiotic tolerances and dispersal limitations.

Autonomic control of pancreatic beta cells: What is known on the regulation of insulin secretion and beta-cell proliferation in rodents and humans.

Insulin secretion and pancreatic beta-cell proliferation are tightly regulated by several signals such as hormones, nutrients, and neurotransmitters. However, the autonomic control of beta cells is not fully understood. In this review, we describe mechanisms involved in insulin secretion as well as metabolic and mitogenic actions on its target tissues. Since pancreatic islets are physically connected to the brain by nerves, parasympathetic and sympathetic neurotransmitters can directly potentiate or repress insulin secretion and beta-cell proliferation. Finally, we highlight the role of the autonomic nervous system in metabolic diseases such as diabetes and obesity.

The impact of bio-logging on body weight change of the Eurasian beaver.

Bio-logging is a common method to collect ecological data on wild animals, but might also induce stress, reduce body condition, and alter behavior. Eurasian beavers (Castor fiber) are a semi-aquatic and nocturnal species that are challenging to observe in the wild. Bio-loggers are hence useful tools to study their behaviour and movements, but this raises concerns of potential negative impacts of tagging. To investigate the potential negative impacts of glue-on tags, we compared body weight change for tagged and untagged Eurasian beavers. We hypothesized that tagged beavers would gain less body weight compared to untagged beavers, and that weight change might be affected by tagging length, tag weight, water temperature and the season of tagging. Daily percentage body weight change in relation to initial body weight during the first capture was compared during 57 tagging periods (18±7 days) and 32 controls periods (64±47 days). Body weight change varied between the two groups, with untagged beavers on average gaining daily weight whilst tagged beavers on average lost weight daily, indicating a negative effect of tagging. The average reduction in percentage body weight change per day for tagged beavers was small (0.1 ± 0.3%), and with large individual variation. Neither tag weight, number of tagging days, nor season were important in explaining body weight change of tagged animals. In other words, we found that tagging reduced daily body weight during the tagging period but were unable to determine the mechanism(s) responsible for this decline. Detrimental effects of tagging have important implications for animal welfare and can introduce bias in data that are collected. This calls for careful consideration in the use of tags. We conclude that studies investigating the effects of tagging should consider individual variation in the effects of tagging and, where possible, compare tagged animals with a control group.

Seed traits and rodent community interact to determine seed fate: evidence from both enclosure and field experiments.

Animal-mediated seed dispersal is an important ecological process in which a strong mutualism between animals and plants can arise. However, few studies have examined how a community of potential seed dispersers interacts with sympatric seed trees. We employed a series of experiments in the Qinling Mountains in both semi-natural enclosure and the field to assess the interactions among 3 sympatric rodent species and 3 Fagaceae tree seeds. Seed traits all showed similar tannin levels but markedly different physical traits and nutritional contents. We found that seeds with heavy weight, thick coat, and high nutritional contents were less likely to be eaten in situ but more often to be eaten after dispersal or hoarded by rodents. These results support both the handling time hypothesis and the high nutrition hypothesis. Surprisingly, we also found that rodents, maybe, preferred to consume seeds with low levels of crude fiber in situ, and to harvest and hoard those with high levels of crude fiber for later consumption. The sympatric rodent species, Cansumys canus, the largest rodent in our study, harvested and hoarded more Quercus variabilis seeds with high physical and nutritional traits, while Apodemus draco, the smallest rodent, harvested more Q. serrata seeds with low physical and nutritional traits, and Niviventer confucianus harvested and hoarded more Q. aliena seeds with medium physical and nutritional traits. Our study demonstrates that different seed traits play different roles in influencing the seed fate and the shaping of mutualism and predation interactions within a community of rodent species.

Increased burrow oxygen levels trigger defensive burrow-sealing behavior by plateau zokors.

Defensive behaviors are a response to immediate and potential threats in the environment, including abiotic and biotic threats. Subterranean rodents exhibit morphological and physiological adaptions for life underground, and they will seal with mounds and additional plugs when their burrow opened. However, little is known about the factors driving this defensive behavior. In this study, we selected a subterranean rodent, plateau zokor (Myospalax fontanieri), as a species to investigate (both in the laboratory and in the field) the possible factors responsible for burrow-sealing behavior. Our results showed that: (1) In the laboratory, the burrow-sealing frequency of plateau zokor in response to five factors were as follows: oxygen (52.63%) > light (34.58%) > temperature (20.24%) > gas flow (6.48%) > sound/control (0%). Except for light, the burrow-sealing frequency in response to other factors was significantly lower than that in response to oxygen (P < 0.05). (2) Burrow-sealing behavior in response to each treatment did not differ significantly between males and females in the laboratory experiment. (3) In the field, during the animal's active periods in both the cold and warm season, the burrow-sealing frequency under the oxygen treatment was higher than that under the light and temperature treatments. Plateau zokors were found not to be sensitive to these treatments during their inactive periods during both the cold and warm season. (4) The latency to reseal the burrow showed no obvious differences between each treatment both in the laboratory and in the field. In conclusion, the main factor that influences the burrow-sealing behavior of plateau zokors is the variation in oxygen concentration, and this defensive behavior is related to their activity rhythm.

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

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