Analysis of the differences in body and skull indices among populations
ANOVA revealed significant inter-population differences in multiple phenotypic traits (Table 1). Body weight showed no significant difference among regions (
P>0.05). Body length varied significantly (
P<0.01), with XX population the largest (significantly greater than WZ, LC, LF and ZQ) and ZQ/WZ the smallest. Tail length was greatest in QX and smallest in WT. Tail-to-body ratio differed highly significantly (
P<0.01), with LF, QX and ZQ significantly higher than WT, which was the lowest. Forelimb lengths differed significantly (
P<0.01), with WZ and KL significantly longer than SY and WT; hindlimb lengths also differed significantly (
P<0.01), with QX and WZ significantly longer than WT. Ear length differed significantly (
P<0.01), with KL and XX significantly longer than LC, WZ and ZQ (LC shortest). Ear width differed significantly (
P<0.01), with LF (left) and QX (right) largest and WT smallest.
Most cranial traits showed no significant regional differences (
P>0.05), except cranial height (
P<0.01), which was highest in XX and significantly higher than in LC, QX and WT. Other traits (cranial length, condylobasal length, occipitonasal length, zygomatic width, interorbital width, skull width, ear bulla length and upper/lower tooth row length) showed no significant differences (
P>0.05).
Analysis of phenotypic trait variation within and between populations
Variance component analysis for 22 traits (Table 2) revealed significant variance distribution differences among populations. Random error: 46.64%-91.01% (mean 69.91%); within-population: 4.71%-9.12% (mean 6.92%); among-population: 0-48.75% (mean 23.17%). Among-population variance generally exceeded within-population variance for most traits.
Body weight showed the lowest VST (6.96%), with random error variance at 90.39%. Body surface traits showed high differentiation: VST values for body length, tail length and tail-to-body ratio were 75.64%, 88.68% and 91.11%, respectively. Limb length VST ranged from 86.78% (left forelimb) to 91.37% (right hindlimb); ear-related traits all exceeded 86% (right ear width at 88.73%). For skull traits, skull width (SKW) and cranial height (CH) showed moderately high VST (77.90% and 79.01%), indicating stable inter-population differentiation. However, cranial length (33.02%), condylobasal length (18.50%) and occipitonasal length (30.56%) showed lower VST. Zygomatic width, interorbital width and upper/lower tooth row length had zero among-population variance and VST = 0, indicating no detectable inter-population differentiation.
CV for each trait varied significantly among populations (Table 3), with mean CV ranging from 2.17% to 19.08%. Body weight was highest (19.08%), followed by tail length (13.90%), tail-to-body ratio (13.77%) and body length (8.82%). Limb length mean CV: 8.14%-9.88%. Ear length: 10.56% (left) and 10.35% (right); ear width: 8.93% (left) and 9.31% (right). Skull traits showed low variation: cranial length (4.17%), condylobasal length (4.19%), occipitonasal length (3.84%), skull width (2.17%), cranial height (3.06%). Ear bulla length (6.54%) was relatively high among skull traits. Upper and lower tooth row lengths were lowest (2.63% and 2.54%).
Principal component analysis and cluster analysis of morphological indices
PCA revealed PC1, PC2 and PC3 explained 51.24%, 28.01% and 8.99% of variation (cumulative 88.24%, Fig 2). In PC1-PC2 space (Fig 2A), populations formed concentrated clusters with distinct shifts (especially PC1) and some overlap (PC2), suggesting phenotypic similarities. PC1-PC3 (Fig 2B) was consistent with PC1-PC2, with PC3 reflecting secondary trait differences. Hierarchical clustering (Fig 3) divided populations into Group I (QX, ZQ, LC, SY, WT) and Group II (WZ, LF, KL, XX); Group I (except SY) were mainly in the Taihang Mountains, Group II in the Lvliang Mountains.
Bergmann’s rule testing and analysis of C. longicaudatus morphology
Phenotypic indicators were grouped by function and morphology into body surface traits (BL, TL, TL/BL, BW, LIMBS, EAR_SIZE) and skull traits (further integrated into length, width, height and functional axes) to reduce single-indicator interference.
Spearman correlations between 22 traits and 6 environmental factors (altitude, humidity, latitude, longitude, temperature, wind speed) are shown in heatmap (Fig 4). Trait responses varied significantly. Wind speed and latitude showed relatively concentrated correlations with multiple traits (
P<0.05), while altitude, humidity, longitude and temperature showed no correlations with most traits. Latitude was significantly negatively correlated with cranial height (
P<0.01,
R = -0.46), skull width (
P<0.05,
R = 0.30), ear length (
P<0.01,
RLEL = -0.29,
RREL = -0.08), ear width (
P<0.01,
RLEW = -0.41,
RREW = -0.34) and forelimb length (
P<0.01,
RLFL = -0.30,
RRFL = -0.31). Wind speed was significantly positively correlated with skull width (
P<0.05,
R = -0.30), tail length (
P<0.01,
R = 0.27), tail-to-body ratio (
P<0.01,
R = 0.29), ear width (
P<0.01,
RLEW = 0.38,
RREW = 0.43) and limb length (
PLFL < 0.05,
RLFL = 0.26;
PRFL < 0.01,
RRFL = 0.27;
PLHL < 0.01,
RLHL = 0.44;
PRHL < 0.01,
RRHL = 0.43), but negatively correlated with body length (
P<0.05,
R = -0.23). Tail length was negatively correlated with altitude (
P < 0.01,
R = -0.29). Cranial height and ear width were negatively correlated with humidity (
P<0.05,
RCH = -0.32;
RLEW = -0.24,
RREW = -0.22). Ear bulla length was negatively correlated with longitude (
P<0.05,
R = -0.29), but positively correlated with temperature (
P<0.05,
R = 0.26), as was body length (
P<0.05,
R = 0.24). Other traits showed no significant correlations (
P>0.05).
Mantel tests (Fig 5) showed significant correlations mainly between specific environmental gradients and composite traits. Altitude correlated positively with LIMBS and EAR_SIZE; latitude correlated with body weight (strong), body length (weak) and tail/tail-body (weak but stable). For skull traits, latitude correlated positively with all four composite axes (length, width, height, functional), indicating synchronous skull changes along latitude. Longitude and humidity showed weaker/fewer correlations. Overall, trait-environment correlations concentrated between altitude and body surface traits and latitude and skull structure, reflecting differential trait responsiveness.
Bergmann’s rule (body size) and Allen’s rule (appendage length) (
Bergmann, 1848;
Allen, 1877) have species- and gradient-specific applicability (
Tabh and Nord, 2023;
Valladares-Gómez et al., 2024). Our latitude-body size correlations align with most rodents
(Ashton et al., 2000; Balčiauskas et al., 2025;
Tomlinson and Withers, 2009), supporting Bergmann’s rule and altitude-driven skull selection. However, Allen’s rule was contradicted - limb/ear sizes increased with altitude, not latitude, consistent with exceptions in Chilean rodents and
A. draco (
Valladares-Gómez et al., 2024;
Cui et al., 2020).
Differential growth rates drive morphological diversity, with plastic contributions nonlinear and energy-dependent (
Tabh and Nord, 2023). The appendage-altitude positive correlation in
C. longicaudatus may reflect allometric scaling with cooling, as in
A. agrarius,
Rattus norvegicus and
Niviventer confucianus, with climate effects mediated via body size (
Geng and Zhu, 2024). Bird analyses show similar variation across temperature gradients (
Frӧhlich et al., 2023), suggesting this correlation reflects locomotor maintenance rather than thermoregulation. Allometric ratios vary among rodent groups
(Rosalino et al., 2013), suggesting species-specific patterns.
Environmental gradients drive differentiation via selection and gene flow
(Price et al., 2003; Rodríguez-Cajarville et al., 2022;
Cao et al., 2026). High VST (>75%) in body/skull traits, with among-population variance exceeding within-population variance, mirrors patterns in
Toona ciliata (Wang et al., 2019), Tilia mandshurica (Mu et al., 2023) and rodent skulls
(Monteiro et al., 2005; Kang et al., 2024), indicating geographic isolation dominates. However, food/competition can outweigh temperature
(Cui et al., 2022), adjusting niche width
(Stephens et al., 2019) and foraging
(Manlick et al., 2021); thus ecological factors beyond climate matter. Body weight (VST = 6.96%) is season/reproduction-sensitive and not a sole indicator; food availability may drive size changes
(Alston et al., 2023), so multiple linear indicators are preferable.