- Research article
- Open Access
Impact of transition to a subterranean lifestyle on morphological disparity and integration in talpid moles (Mammalia, Talpidae)
BMC Evolutionary Biology volume 19, Article number: 179 (2019)
Understanding the mechanisms promoting or constraining morphological diversification within clades is a central topic in evolutionary biology. Ecological transitions are of particular interest because of their influence upon the selective forces and factors involved in phenotypic evolution. Here we focused on the humerus and mandibles of talpid moles to test whether the transition to the subterranean lifestyle impacted morphological disparity and phenotypic traits covariation between these two structures.
Our results indicate non-subterranean species occupy a significantly larger portion of the talpid moles morphospace. However, there is no difference between subterranean and non-subterranean moles in terms of the strength and direction of phenotypic integration.
Our study shows that the transition to a subterranean lifestyle significantly reduced morphological variability in talpid moles. However, this reduced disparity was not accompanied by changes in the pattern of traits covariation between the humerus and the mandible, suggesting the presence of strong phylogenetic conservatism within this pattern.
Understanding why some clades achieve a large morphological, behavioral and ecological diversity, while others do not, represents a central aim in evolutionary biology. Studies addressing this question usually try to identify the factors allowing high phenotypic diversity, or constraining its realization . Key innovations and the transition to novel ecological niches are generally thought to promote morphological variation of clades, but might also reduce taxonomic diversification via niche specialization and stabilizing selection [1, 2]. These differential effects of trait acquisition or ecological transition on phenotypic diversity between clades are best exemplified by sister clades, where one clade shows higher phenotypic diversity than the others [3, 4]. In this context, it has been proposed that phenotypic traits covariation and modularity could promote differences in disparity among clades . Simulation studies have shown that traits covariation may drive morphological variability along different axes of variation, resulting in either less or more disparity depending on the relationship between selection pressures and the major axes of variation [1, 5]. Ecological transitions are of particular interest because of their large influence on phenotypic diversification [3, 4, 6]. Specifically, while species shifting into new niches are likely to evolve adaptations allowing them to exploit the new ecological settings, those remaining within the ancestral niche should retain plesiomorphic anatomical structures [7, 8]. Textbook examples of this are found among archosaurs , anuran and caecilian amphibians [4, 10], and primates , among others.
A notable ecological transition is represented by the colonization of the subterranean ecotope by different, unrelated taxa (rodents, insectivores and marsupials [12, 13]). Different studies on rodents have demonstrated how the colonization of the subterranean ecotope could dramatically influence species morphology and how digging specialization shapes the evolution of forelimbs and the cranio-dental complex [14,15,16]. One of the most spectacular examples of adaptation to subterranean life is represented by moles. The mammalian family Talpidae includes ambulatorial (Uropsilini), semi-aquatic (Desmanini and Codylurini), semi-fossorial (Urotrichini and Neurotrichini) and fully subterranean species (Scalopini and Talpini). The colonization of such a wide array of environments was realized by the combination of different behavioral, physiological and morphological adaptations [17,18,19,20,21] especially evident in the forelimbs of digging species [17, 18, 22, 23]. Such wide array of morphological adaptations makes talpids an ideal group for investigating the phenotypic effects of a major ecological transition, such as the colonization of the subterranean habitat.
In the present study, we provide a comprehensive morphological analysis of the humerus and mandible of all extant talpid genera. The humerus is widely known to be a good proxy for locomotor behavior in moles [19, 20, 23], while the mandible is obviously linked to feeding [24,25,26,27]. We used a 2D geometric morphometrics to quantitatively assess shape variation in both humeri and mandibles between subterranean and non-subterranean moles. Then, we investigated different aspects of these bones’ shape covariation under a phylogenetically-informed scenario. Specifically, we compared the strength and direction of phenotypic covariation between the two anatomical structures and between subterranean and non-subterranean moles. The need for coordination and integration between the feeding and locomotor apparatuses has been shown for different vertebrate taxa [28, 29]. In particular, highly specialized behaviors may require a more complex interaction between different parts of an organism, resulting in increased covariation . We hypothesized that the species adapted to the subterranean lifestyle might display stronger trait covariation between the humerus and the mandible due to the functional constraints imposed by highly demanding digging kinematics and by the limited range of food items available underground [12, 17]. We further determined if the potential differences in morphological disparity between the species falling in the two ecotopes could be ascribed to variation in levels of phenotypic covariation.
We examined left 365 mandibles and 463 left humeri belonging to adult individuals of 37 extant talpid species. Specimens are stored in the following institutions: ISEZ-PAN, Krakow, Poland, Tsukuba Natural History Museum, Tsukuba, Ibaraki, Japan; Museu de Historia Natural, Lisboa, Portugal; Natural History Museum, London, UK; BSPG, Munich, Germany; Wien Natural History Museum, Wien, Austria; LACM, Los Angeles, USA; UCMP, Berkeley, USA; “La Specola” Museo di Storia Naturale di Firenze, Italy and Museo di Anatomia Comparata G.B. Grassi, “Sapienza” Università di Roma, Rome, Italy. Details about the sample are summarized in Table 1. We separated the species under investigation into two groups: non-subterranean and subterranean. Following [12, 13, 31], subterranean species spend most of their life underground and come above the ground only incidentally, whereas non-subterranean species lack the extreme underground specializations and spend a considerable amount of time above the ground (e.g. foraging). Groupings at the species level are reported in Table 1 (Further details on the sampling effort can be found in Additional file 1). Clade assignation follows .
The phylogenetic history of Talpidae is highly debated. Despite a growing number of papers on the subject matter, there is still a lack of agreement between the different phylogenetic hypotheses [18, 21, 33,34,35,36,37,38,39,40]. In particular, morphological and molecular approaches conflict upon the position of Scalopini, the monophyly of Urotrichini and the position of Condylura [20, 32, 38,39,40,41]. Since a systematic revision of Talpidae phylogeny is beyond the scope of the present work, we decided to use two different phylogenetic hypotheses when using phylogenetic comparative methods. The first is based on molecular data, where Neurotrichini (Scaptonyx and Neurotrichus) are a polyphyletic group . However, because this phylogenetic tree does not include all the species under investigation in the present paper, we built a synthetic phylogeny (an informal supertree using the Mesquite software ) supplemented with additional data to resolve the relationships within the genus Talpa  and within the genus Euroscaptor . The second phylogenetic hypothesis is based on trees produced by maximum parsimony cladistic analysis of morphological characters based on a published character matrix we developed elsewhere . The character matrices were analysed using PAUP 4.0 a147 [43, 44] using a heuristic search and stepwise addition, with a random addition sequence of 1000 replicates. The phylogenetic comparative analyses (see below) were applied on a strict consensus tree computed on the three most parsimonious trees found. The character matrix and relative character list are presented in Additional file 2. The time calibration for branches has been derived from a thorough review of the palaeontological literature on the subject matter [20, 31, 39]. For the time calibration we considered: the ages of the first occurrence of extant species and the molecular clock estimate (when available). The time calibration has been performed using the Stratigraphic Tool in the Mesquite software . The two phylogenies are presented in Fig. 1, further details on the cladistics analysis are presented in Additional file 2, whereas detailed information about the time calibration are presented in Additional file 3.
Geometric morphometrics and shape analysis
Humeri and mandibles were photographed in caudal view at a fixed distance of 50 cm using a Nikon D100 camera with a Micro-Nikkor 105-mm lens. We digitized 22 landmarks and 14 semi-landmarks on the humerus and 12 landmarks and 26 semi-landmarks on the mandible (Fig. 2a and b) using the tpsDig2 software . The humeral and mandibular landmark configurations were derived from [40, 46,47,48] respectively. Semi-landmarks were used to capture the morphology of complex outlines where homologous anatomical points were missing. Semi-landmarks assume that curves or contours are homologous among specimens . Successively, a generalized Procrustes analysis (GPA) , implemented with the procSym() function in the R-package “Morpho” , was used to rotate, translate, and scale landmark configurations to unit centroid size (CS = the square root of the sum of squared distances of the landmarks from their centroid ). To visualize the multivariate ordination of the aligned specimens, we performed a between-group PCA (bgPCA), using the function groupPCA() included in the R-package “Morpho”, considering the species as groups. The bgPCA provides a projection of the data onto the principal components of the group means, resulting in an ordination of the shape variables between the group means. The new axes are orthogonal and can be computed even when the per-group data matrices are not of full rank. This method provides a good performance when the number of observations in each group is smaller than the number of variables [53, 54].
Because we have a different number of specimens for humeri and mandibles for each species, all of the following analyses were performed on per-species averaged data. The significance of the observed shape changes between subterranean and non-subterranean species was evaluated by performing a Procrustes ANOVA on aligned Procrustes coordinates using the function procD.lm() included in the R package “geomorph” . To measure shape disparity among non-subterranean and subterranean species we used Procrustes variance, which is the sum of the diagonal elements of the group covariance matrix divided by the number of observations in the group using the function morphol.disparity() from the R package “geomorph” [56, 57]. In order to visualize shape changes in ordination plots we chose to use the method described in . There it was suggested that a useful way to visualize local, infinitesimal variation within a deformation grid is to use the Jacobian (J) of the Thin Plate Spline interpolation function. J captures very local information as localized variation in the non-affine component of the deformation using derivatives of the used interpolation function (TPS in our case). In 2D J is a 2 × 2 matrix that can be evaluated at any point within a body. The logarithm of its determinant represents the change in the area in the region about the interpolation point. Values < 0 indicate that, with respect to the source (here the sample consensus), the target (here the PC’s extremes) experiences a reduction in the local area, while values > 0 indicate an enlargement.
Evolutionary allometry and size correction
Multivariate regressions between shape and size data were applied to determine the presence of evolutionary allometry in both the humerus and the mandible. To test for differences in slope among subterranean and non-subterranean species, we ran a permutational multivariate analysis of covariance, using Procrustes coordinates as dependent variables, centroid size (CS) as an independent variable and the two groups as factor [58, 59]. To test the effect of size on morphological disparity the CS was included in the model as a covariate. In order to test the potential effect of size on the pattern of morphological covariation we repeated each of the analyses described below after computing size-free Procrustes coordinates.
The measurement error associated with the digitization of landmarks was measured on three replicates of 100 specimens for each dataset (humeri plus mandibles). The mean Procrustes distances between all the combinations of pairs of specimens were computed for each replicated dataset using the TPSsmall software . We calculated the mean Procrustes distances for each triplet of the same subjects occurring in the three replicas. We then computed the averages of all the mean values of the minimum and maximum values of each triplet. The amount of digitization error, with respect to the total variation in the shape, can be expressed as a percentage. We calculated the ratio of the mean value for total digitization and the mean of the total dataset.
Trait covariation, strength and direction
It has been recently noted  that sliding semi-landmarks using the minimum bending energy (BEN) approach may result in increased covariation between modules. Because we used semi-landmarks in our dataset, we repeated all the following integration analyses using shape coordinates derived using both the minimum BEN and minimum Procrustes distances (PRD) approaches in order to evaluate any potential discrepancy in the results. We report here that we did not find any significant discrepancy when using either sliding methods, hence we present only the results obtained from the analyses performed on the shape coordinates derived after using the minimum BEN approach.
We assessed the covariation between the humerus and the mandible using partial least squares (PLS) analyses [52, 62]. PLS is suitable for the study of covariation between two sets of variables in several groups. We quantified the covariation for each pair of axes by means of a correlation coefficient, whose significance is addressed by means of permutation under the null hypothesis that the distribution of specimens on one axis has no bearing on the distribution on the other axis . Adams and Collyer  proposed a new strategy to compare the strength of PLS focusing on the first singular vector pair. The authors proposed a standardized test statistic (a z-score) for measuring the degree of morphological integration between sets of variables. The z-scores can be used to test for differences (via ANOVA) among groups. We used the compare.pls() function from the R package “geomorph” to compare the effect sizes, measured as standard deviates, z, and performs two-sample z-tests, using the pooled standard error from the sampling distributions of the PLS analyses. This tests for differences in the strength of covariation, whereas nothing is known about its direction. The orientation of integration patterns in space can be interpreted as the rate of shape changes in one module relative to the rate of shape changes in the other. This aspect is very important as it may reveal whether a common pattern of shape changes within modules exists between clades. In fact, groups may show similar integration coefficients but have different integration patterns [65, 66]. In order to investigate this issue, we performed separate major axis (MA) analyses on the different clades shapes on the space identified by the first pair of PLS axes [65, 66]. MA is particularly suitable here because of its ‘symmetry,’ i.e., residuals are computed orthogonally to the line of best fit and this is coherent with PLS aims. It does not require the classic assumption of dependence-independence relationship . MA slopes were then compared through pairwise ANOVA, using lifestyle categories as groups.
Recently, Bookstein  proposed a new method to evaluate morphological integration “intrinsically” to a structure. This method tests the null hypothesis of “self-similarity” (e.g., the absence of any interpretable change at any spatial scale) in a collection of shapes and it is based on a linear regression of log partial warps variance against their proper log bending energy (i.e., the log of eigenvalues of the bending energy matrix computed on the consensus). Here, a regression slope less than − 1 indicates “integration” whereas a slope greater than − 1 indicates “dis-integration”. If the regression slope is exactly − 1 data can be considered “self-similar” (for further details refer to  and to Additional file 2). Finally, we compared the resulting slopes between subterranean and non-subterranean species using the R package RRPP .
Phylogenetic non-independence and phylogenetic signal
The phylogenetic signal was calculated for the shape data using the Kmult statistic, a method that measures the similarity of trait values in relation to a Brownian motion model of evolution. It is specifically designed to address the challenge of working with high-dimensional landmark configurations .
The significance of the observed shape changes between subterranean and non-subterranean species was evaluated by performing a Procrustes ANOVA in a phylogenetic framework on aligned Procrustes coordinates using the function procD.pgls() included in the R package “geomorph” [49, 55]. We quantified the degree of phylogenetic morphological integration between the humerus and the mandible using partial least squares (PLS) analyses under the Brownian model of evolution using the function phylo.integration() from the geomorph R library .
The digitization errors in the humeral and mandibular datasets were as low as 0.8 and 0.6% of the total variation, respectively. Because the measurement error was smaller than 5% in both datasets it could be safely assumed its effect on the results was negligible .
Humerus. The bgPCA plot shows that subterranean and non-subterranean moles are well separated in the morphospace (Fig. 3). Procrustes ANOVA performed on the shape variables returned highly significant difference (p-value =0.001). Along the PC1 (87.4% of the total variance) there is a clear distinction between the subterranean (negative values) and non-subterranean (positive values) species. The humerus, at negative values shows a robust configuration with a relevant enlargement of the proximal region, in particular of the pectoral crest and the teres tubercle. At positive values the humerus shows a slender configuration with an overall reduced proximal region. Along the PC2 (3.7% of the total variance) it is possible to separate the Talpini (from consensus to negative values) from the Scalopini (positive values). At negative values the humerus shows a wide and elongated pectoral crest, while at negative values the pectoral crest is rounded and short.
Mandible. The bgPCA plot showed that subterranean and non-subterranean moles are well separated in the morphospace (Fig. 4). Procrustes ANOVA performed on the shape variables returned a highly significant result (p-value = 0.001). Along the PC1 (45.7% of the total variance) there is a clear distinction between the subterranean (positive values) and non-subterranean (negative values) species. At negative values the mandible shows a straight and robust horizontal ramus, while the coronoid process is large and elongated, the condylar process is short and the angular process is pointed and slender. At positive values the mandible shows a bent horizontal ramus, while the coronoid process is slender and pointed, the condylar process is elongated and the angular process is robust and rounded. Along the PC2 (16.1% of the total variance) it is possible to separate the taxa showing a semi-aquatic lifestyle from all of the other moles. In particular, Condylura occupies the region of the morphospace marked by positive values, distinguished by a very slender and bent horizontal ramus, the coronoid and angular processes are extremely reduced, while the condylar process is robust and elongated. Desmans occupy the region of the morphospace marked by negative values, where the mandible shows a robust and straight horizontal ramus and an overall robust condylar region.
The multivariate regressions revealed a significant interaction between shape and size for both the humerus and the mandible (F = 10.08, r2 = 0.22, p-value = 0.001; F = 14.25, r2 = 0.05, p-value = 0.001; respectively). Allometric trajectories were different between subterranean and non-subterranean taxa for both the humerus and the mandible (interaction p-value = 0.001, p-value = 0.003; respectively).
Results from pairwise comparison of Procrustes variance are summarised in Table 2. For each structure (humerus and mandible) non-subterranean species always showed a significantly higher morphological disparity. Finally, the humerus proved to have an overall greater disparity when compared to the mandible (Procrustes variances: 0.018 and 0.0046, respectively). The inclusion of size (CS) as covariate did not impact the results from the previous analyses for both the humerus (p-value = 0.004) and the mandible (p-value = 0.001).
Phenotypic trait covariation
Results from the PLS analyses are summarized in Table 3. We found that the humerus and the mandible are highly integrated with each other, and subterranean and non-subterranean taxa showed similarly high degree of covariation. When we repeated the PLS analyses on the size-corrected shape coordinates we did not find any difference in the significance level compared to the analyses conducted on the Procrustes coordinates (whole sample: r-pls = 0.90, p-value = 0.001; subterranean: r-pls = 0.87, p-value = 0.007; non-subterranean: r-pls = 0.92, p-value = 0.003).
Table 4 summarises the result of the PLSs effect sizes comparison and of the MA analyses. Overall, subterranean and non-subterranean taxa displayed non-significant differences in both strength and direction of covariation.
The results of PLS analyses are reported in Fig. 5. In the PLS plot of mandible and humerus the subterranean species placed at negative extremes of the two PLS axes, while non-subterranean taxa placed at positive extremes. At negative extremes the humerus showed the robust humeral configuration with the enlarged proximal region, while the mandible showed a thinner horizontal ramus, a shorter coronoid process and an enlarged angular process. At positive extremes the humerus showed the slender configuration with reduced proximal region, while the mandible showed a robust horizontal ramus, an elongated coronoid process and reduced angular process.
For both humerus and mandible the partial warp variance drops faster than the bending energy rises. The resulting slope of − 1.31 for the humerus suggests an integrated pattern. The same holds for the mandible where the regression slope of − 2.04 indicates a highly integrated pattern. We did not find any significant difference between subterranean and non-subterranean species in the degree of humeral integration (regression slopes: − 1.28; − 1.25, respectively; p-value = 0.573). Again, we did not find any significant difference between subterranean and non-subterranean species in the degree of mandibular integration (regression slopes: − 1.99; − 2.07, respectively; p-value = 0.642).
Phylogenetic non-independence and phylogenetic signal
Phylogenetic signal in the aligned Procrustes coordinates is high for both the humerus and the mandible. Similar results were obtained when using both morphological and molecular phylogenetic hypotheses (see Table 5). The correlation coefficient was non-significant for the phylogenetically-informed version of PLS when using both morphological and molecular phylogenetic hypotheses (r = 0.563, p-value = 0.268; r = 0.565, p-value = 0.216; respectively). However, we found the humerus and mandible to be significantly correlated in subterranean moles when using both morphological and molecular phylogenetic hypotheses (r = 0.895, p-value = 0.016; r = 0.893, p-value = 0.017; respectively), whereas correlation was not significant in non-subterranean moles (r = 0.726, p-value = 0.123; r = 0.718, p-value = 0.124; respectively). Similar results were obtained when we removed the size effect from the shape data when using both morphological and molecular phylogenetic hypotheses. The correlation was significant for subterranean moles (r = 0.815, p-value = 0.018; r = 0.895, p-value = 0.016; respectively), whereas correlation was not significant in non-subterranean moles (r = 0.706, p-value = 0.146; r = 0.701, p-value = 0.138; respectively).
Polly et al., 2013 noted that when running the phylogenetic version of the PLS analysis, shape changes associated with PLS axes could not be biologically interpretable, possibly reflecting the removal of the phylogenetic component associated with the morphological adaptive signal (see Additional file 2 for further details). Therefore, we will present and discuss only the PLS distribution, and associated shape changes, obtained prior to phylogenetic correction.
During their evolution, talpid moles diversified into a number of ecological niches and geographical areas [12, 18, 19, 71, 72]. The colonization of the subterranean environment is certainly the largest ecological transition ever experienced by the clade. It represents one of the most remarkable suites of adaptation showed by any mammalian group [17, 73]. Our study demonstrates that the transition to subterranean environments resulted in dramatically reduced shape disparity in both the humerus and the mandible of subterranean species (see Figs. 3 and 4). The subterranean ecotope is structurally simple, relatively stable and highly demanding in terms of locomotion [12, 13]. These features require a high degree of specialization, forcing species within a narrow ecological niche . Therefore, unrelated species evolving in this simple, but highly demanding environment are expected to display a high degree of phenotypic convergence .
Our results suggest that the two fully subterranean mole tribes (Talpini and Scalopini) have humeral [20, 39] and mandibular shape variation significantly reduced by functional constraints that are imposed by their highly specialized lifestyle.
The humerus of subterranean species is highly adapted for digging. Its evolution is characterized by slow evolutionary rates, convergent allometric trajectories among different talpid clades and, overall, a strong conservatism, as suggested by the presence of a strong phylogenetic signal [19,20,21, 39]. The same considerations hold true for the mandible, where the subterranean species display remarkable morphological similarity (see Fig. 4) despite the different degree of dental reduction in the two tribes . The diet of subterranean moles mainly includes earthworms (Annelida, up to 90% in some species; [75,76,77,78]), while non-subterranean moles have a broader range of food items including beetles, fish, crustaceans, plant material and seeds [45, 79,80,81]. However, the mandible showed a lower, yet significant, phylogenetic signal. This could be a consequence of the evolution of different adaptive strategies within a monophyletic group. Examples are represented by the evolution of hypsodont dentition in the urotrichine genus Urotrichus (which contributes to the ecological separation form Dymecodon ), and the independent evolution of the semi-aquatic lifestyle in desmans and Condylura, which separate them from the other talpids along PC2. Furthermore, the star-nosed mole (Condylura) displays a highly derived oral apparatus designed for underwater high-velocity food consumption . The derived condition in Condylura is somewhat typical for mammals adapted to a semi-aquatic lifestyle [3, 82,83,84].
We did not find any significant difference between subterranean and non-subterranean taxa in either strength or direction of trait covariation. We found a strong correlation between the humerus and the mandible in the whole sample and in the per-group separated analyses. Though, when accounting for phylogeny on the whole sample, we found that the correlation between the humerus and the mandible to be no longer significant. This could highlight the presence of a strong phylogenetic structure suggesting that the humerus and the mandible could have followed separated evolutionary pathways. However, we found the humerus and mandible to be significantly correlated only in subterranean moles, when accounting for shared ancestry. This result suggests, on the one hand, how the subterranean species’ shift into a different ecological niche was accompanied by concerted changes in both locomotion and dietary patterns. On the other hand, it reflects how the lack of covariation between the humerus and mandible of non-subterranean species might imply an independent evolution of the two structures, likely subject to lower environmental constraints or different functional trade-offs. In this case, greater evolutionary lability could have played a role in the adaptation of non-subterranean species to different environments [1, 4, 5, 85]. However, similar patterns of trait covariation between subterranean and non-subterranean taxa suggest the presence of a strong phylogenetic conservatism in talpids. Strong trait covariation is usually associated with reduced phenotypic variability, while modularity, by breaking down patterns of covariation between structures, might increase the number of possible axes of variation along which the phenotypes might diversify . However, module covariation can generate either more or less diversity according to the selective pressures acting on the principal axes of variation [1, 5, 86]. In the context of the present study, it is possible that high trait covariation might have constrained talpid evolution along lines of least evolutionary resistance, meaning that developmental processes might offer simple pathways to generate variation [87, 88]. It has been demonstrated that subterranean moles display an early autopodial chondrification as compared to non-subterranean taxa [89, 90]. This evidence suggests that subterranean moles might experience high levels of phenotypic covariation earlier during their ontogenetic development, resulting in the reduced morphological variability showed at the adult stage. A similar pattern was observed in marsupials where high integration could enhance the effect of the intense functional demand for continuous suckling earlier in their development, hence resulting in limited variability of the marsupial oral apparatus . Finally, we did not find any significant impact of size on the patterns of trait covariation in talpids, as well as on their morphological disparity. These results suggest that allometry could have played a major role in shaping talpids rates of morphological evolution , without significantly impacting the patterns of trait covariation. Nonetheless, changes in shape and proportions may represent another potential factor in restricting disparity [29, 30, 41, 91]. Sansalone et al.  demonstrated that the convergence of allometric trajectories between subterranean moles (Talpini and Scalopini) constrained the humeral shape to a restricted region of the morphospace. In this case, the response to a strong selective pressure may have resulted in the evolution of high covariation and evolutionary allometry. A consequence of these mechanisms is that subterranean mole morphology is extremely constrained [19, 20, 39].
Our study showed that the transition to the subterranean ecotope resulted in a significant loss of disparity in the subterranean clades compared to the non-subterranean species, probably triggered by high morphological covariation. We showed that a strong phylogenetic conservatism in the covariation patterns (strength and direction) might have played a fundamental role in constraining the axes of variation along which subterranean moles were able to attain the high degree of phenotypic specialization necessary to colonize the subterranean environment.
Availability of data and materials
The datasets supporting the conclusions of the present study are available at https://doi.org/10.25952/5bbbddb2d3522
Univariate analysis of variance.
between-group Principal Component Analysis
Generalized Procrustes Analysis
- Institutional abbreviations:
Bayerische Staatssamlung Für Paläontologie und Geologie
Instytut Systematyki i Ewolucji Zwierząt-Polskiej Akademii Nauk
Los Angeles County Museum
Partial Least Squares
University of California Museum of Paleontology
Goswami A, Randau M, Polly PD, Weisbecker V, Bennett CV, Hautier L, Sánchez-Villagra MR. Do developmental constraints and high integration limit the evolution of the marsupial oral apparatus? Int Comp Biol. 2016;56:404–15.
Raia P, Fortelius M. Cope’s law of the unspecialized, Cope’s rule, and weak directionality in evolution. Evol Ecol Res. 2013;15:747–56.
Jones KE, Smaers JB, Goswami A. Impact of the terrestrial-aquatic transition on disparity and rates of evolution in the carnivoran skull. BMC Evol Biol. 2015;15:8.
Vidal-García M, Keogh JS. Phylogenetic conservatism in skulls and evolutionary lability in limbs–morphological evolution across an ancient frog radiation is shaped by diet, locomotion and burrowing. BMC Evol Biol. 2017;17:165.
Goswami A, Smaers JB, Soligo C, Polly PD. The macroevolutionary consequences of phenotypic integration: from development to deep time. Phil Trans R Soc B. 2014;369:20130254.
Dial KP, Shubin N. Brainerd EL. Great transformations in vertebrate evolution: University of Chicago Press; 2015.
Losos JB. Adaptive radiation, ecological opportunity, and evolutionary determinism. Am Nat. 2010;175:623–39.
Losos JB. Convergence, adaptation, and constraint. Evolution. 2011;65:1827–40.
Brusatte SL, Nesbitt SJ, Irmis RB, Butler RJ, Benton MJ, Norell MA. The origin and early radiation of dinosaurs. Earth Sci Rev. 2010;101:68–100.
Sherratt E, Gower DJ, Klingenberg CP, Wilkinson M. Evolution of cranial shape in caecilians (Amphibia: Gymnophiona). Evol Biol. 2014;41:528–45.
Fleagle JG. Primate adaptation and evolution: Academic Press; 2013.
Nevo E. Adaptive convergence and divergence of subterranean mammals. Ann Rev Ecol Syst. 1979;10:269–308.
Nevo E. Mammalian evolution underground. The ecological-genetic-phenetic interfaces. Acta Theriol. 1995;40:9–31.
Steiner-Souza F, De Freitas TR, Cordeiro-Estrela P. Inferring adaptation within shape diversity of the humerus of subterranean rodent Ctenomys. Biol J Linn Soc. 2010;100:353–67.
Marcy AE, Hadly EA, Sherratt E, Garland K, Weisbecker V. Getting a head in hard soils: convergent skull evolution and divergent allometric patterns explain shape variation in a highly diverse genus of pocket gophers (Thomomys). BMC Evol Biol. 2016;16:207.
McIntosh AF, Cox PG. The impact of digging on the evolution of the rodent mandible. J Morphol. 2019;280:176–83.
Yates TL, Moore DW. Speciation and evolution in the family Talpidae (Mammalia: Insectivora). Prog Clin Biol Res. 1990;335:1.
Sánchez-Villagra MR, Horovitz I, Motokawa M. A comprehensive morphological analysis of talpid moles (Mammalia) phylogenetic relationships. Cladistics. 2006;22:59–88.
Piras P, Sansalone G, Teresi L, Kotsakis T, Colangelo P, Loy A. Testing convergent and parallel adaptations in talpids humeral mechanical performance by means of geometric morphometrics and finite element analysis. J Morphol. 2012;273:696–711.
Piras P, Sansalone G, Teresi L, Moscato M, Profico A, Eng R, Cox TC, Loy A, Colangelo P, Kotsakis T. Digging adaptation in insectivorous subterranean eutherians. The enigma of Mesoscalops montanensis unveiled by geometric morphometrics and finite element analysis. J Morphol. 2015;276:1157–71.
Schwermann AH, Thompson RS. Extraordinarily preserved talpids (Mammalia, Lipotyphla) and the evolution of fossoriality. J Vert Paleontol. 2015;4:35.
Barnosky AD. A skeleton of Mesoscalops (Mammalia, Insectivora) from the Miocene Deep River formation, Montana, and a review of the proscalopid moles: evolutionary, functional, and stratigraphic relationships. J Vert Paleontol. 1981;1:285–339.
Gambaryan P, Gasc JP, Renous S. Cinefluorographical study of the burrowing movements in the common mole, Talpa europaea (Lipotyphla, Talpidae). Russ J Theriol. 2003;1:91–109.
Meloro C, Raia P, Piras P, Barbera C, O'Higgins PA. The shape of the mandibular corpus in large fissiped carnivores: allometry, function and phylogeny. Zool J Linnean Soc. 2008;154:832–45.
Meloro C. Raia P (2010) cats and dogs down the tree: the tempo and mode of evolution in the lower carnassial of fossil and living Carnivora. Evol Biol. 2010;37:177–86.
Meloro C, Cáceres NC, Carotenuto F, Sponchiado J, Melo GL, Passaro F, Raia P. Chewing on the trees: constraints and adaptation in the evolution of the primate mandible. Evolution. 2015;69:1690–700.
Piras P, Maiorino L, Teresi L, Meloro C, Lucci F, Kotsakis T, Raia P. Bite of the cats: relationships between functional integration and mechanical performance as revealed by mandible geometry. Syst Biol. 2013;62:878–900.
Montuelle SJ, Herrel A, Libourel PA, Daillie S, Bels VL. Prey capture in lizards: differences in jaw–neck–forelimb coordination. Biol J Linn Soc. 2012;105:607–22.
Brown JG. Jaw function in Smilodon fatalis: a reevaluation of the canine shear-bite and a proposal for a new forelimb-powered class 1 lever model. PLoS One. 2014;9:e107456.
Kane EA, Higham TE. Complex systems are more than the sum of their parts: using integration to understand performance, biomechanics, and diversity. Int Comp Biol. 2015;55:146–65.
Nevo E. Mosaic evolution of subterranean mammals: regression, progression and global convergence.1999 Oxford Univ press. Oxford. .
Sansalone G, Kotsakis T, Schwermann AH, Van den Hoek Ostende LW, Piras P. When moles became diggers: Tegulariscaptor gen. Nov., from the early Oligocene of South Germany, and the evolution of talpid fossoriality. J Syst Palaeontol. 2018;16:645–57.
Motokawa M. Phylogenetic relationships within the family Talpidae (Mammalia: Insectivora). J Zool. 2004;263:147–57.
Shinohara A, Suzuki H, Tsuchiya K, Zhang YP, Luo J, Jiang XL, Wang YX, Campbell KL. Evolution and biogeography of talpid moles from continental East Asia and the Japanese Islands inferred from mitochondrial and nuclear gene sequences. Zool Sci. 2004;21:1177–85.
Cabria MT, Rubines J, Gómez-Moliner B, Zardoya R. On the phylogenetic position of a rare Iberian endemic mammal, the Pyrenean desman (Galemys pyrenaicus). Gene. 2006;375:1–3.
Colangelo P, Bannikova AA, Kryštufek B, Lebedev VS, Annesi F, Capanna E, Loy A. Molecular systematics and evolutionary biogeography of the genus Talpa (Soricomorpha: Talpidae). Mol Phyl Evol. 2010;55:72–80.
Hooker JJ. Skeletal adaptations and phylogeny of the oldest mole Eotalpa (Talpidae, Lipotyphla, Mammalia) from the UK Eocene: the beginning of fossoriality in moles. Palaeontology. 2016;59:195–216.
Bannikova AA, Zemlemerova ED, Colangelo P, Sözen M, Sevindik M, Kidov AA, Dzuev RI, Kryštufek B, Lebedev VS. An underground burst of diversity–a new look at the phylogeny and taxonomy of the genus Talpa Linnaeus, 1758 (Mammalia: Talpidae) as revealed by nuclear and mitochondrial genes. Zool J Linnean Soc. 2015;175:30–48.
Sansalone G, Colangelo P, Kotsakis T, Loy A, Castiglia R, Bannikova AA, Zemlemerova ED, Piras P. Influence of evolutionary allometry on rates of morphological evolution and disparity in strictly subterranean moles (Talpinae, Talpidae, Lipotyphla, Mammalia). J Mammal Evol. 2018;25:1–14.
He K, Shinohara A, Helgen KM, Springer MS, Jiang XL, Campbell KL. Talpid mole phylogeny unites shrew moles and illuminates overlooked cryptic species diversity. Mol Biol Evol. 2016;6:2–21.
He K, Shinohara A, Jiang XL, Campbell KL. Multilocus phylogeny of talpine moles (Talpini, Talpidae, Eulipotyphla) and its implications for systematics. Mol Phyl Evol. 2014;70:13–21.
Maddison WP, Maddison DR. Mesquite: a modular system for evolutionary analysis. Version. 2018:3.51 http://www.mesquiteproject.org.
Swofford DL. PAUP. Phylogenetic analysis using parsimony (and other methods). Version 4. Sinauer, Sunderland, Massachusetts, 2002 142 pp.
Piras P, Buscalioni A. Diplocynodon muelleri comb. nov., an Oligocene diplocynodontine alligatoroid from Catalonia (Ebro Basin, Lleida province, Spain). J Vert Paleont. 2006;26:608–20.
Rohlf FJ. The tps series of software. Hystrix. 2015;26.
Sansalone G, Kotsakis T, Piras P. New systematic insights about Plio-Pleistocene moles from Poland. Acta Palaeontol Pol. 2016;61:221–30.
Schwermann AH, He K, Peters BJ, Plogschties T, Sansalone G. Systematics and macroevolution of extant and fossil scalopine moles (Mammalia, Talpidae). Palaeontology. 2019. https://doi.org/10.1111/pala.12422.
Sansalone G. Evolution of hypsodonty reveals a long-standing ecological separation in the Japanese shrew-moles. J Zool. 2015;297:46–55.
Perez SI, Bernal V, Gonzalez PN. Differences between sliding semilandmark methods in geometric morphometrics, with an application to human cranio-facial and dental variation. J Anat. 2006;208:769–84.
Bookstein FL. Size and shape spaces for landmark data in two dimensions. Stat Sci. 1986;1:181–222.
Schlager S. Morpho: calculations and visualisations related to Geometric Morphometrics. 2014. R-package version 2.0.3–1.
Bookstein FL. Morphometric tools for landmark data. Cambridge: Cambridge University Press; 1991.
Boulesteix AL. A note on between-group PCA. Internatl J Pure Appl Math. 2005;19:359–66.
Mitteroecker P, Bookstein F. Linear discrimination, ordination, and the visualization of selection gradients in modern morphometrics. Evol Biol. 2011;38:100–14.
Adams DC, Otarola-Castillo E. Geomorph: an R package for the collection and analysis of geometric morphometric shape data. Methods Ecol Evol. 2013;4:393–9.
Zelditch ML, Swiderski DL, Sheets HD. Geometric Morphometrics for biologists: a primer. 2nd ed. Academic Press, San Diego: Elsevier; 2012.
Márquez EJ, Cabeen R, Woods RP, Houle D. The measurement of local variation in shape. Evol Biol. 2012;39:419–39.
Piras P, Teresi L, Buscalioni AD, Cubo J. The shadow of forgotten ancestors differently constrains the fate of Alligatoroidea and Crocodyloidea. Global Ecol Biogeography. 2009;18:30–40.
Maiorino L, Farke AA, Kotsakis T, Piras P. Is Torosaurus Triceratops? Geometric morphometric evidence of late Maastrichtian ceratopsid dinosaurs. PLoS One. 2013:1932–6203.
Rohlf FJ. TpsSmall ver. 1.29. 2014. Available at: life.bio.sunysb.edu/morph/soft-dataacq.html.
Cardini A. Integration and modularity in Procrustes shape data: is there a risk of spurious results? BioRxiv. 2018;1:371187.
Rohlf FJ, Corti M. Use of two-block partial least-squares to study covariation in shape. Syst Biol. 2000;49:40–53.
Adams DC. Evaluating modularity in morphometric data: challenges with the RV coefficient and a new test measure. Methods Ecol Evol. 2016;7:65–72.
Adams DC, Collyer ML. On the comparison of the strength of morphological integration across morphometric datasets. Evolution. 2016;70:23–31.
Piras P, Torromeo C, Evangelista A, Gabriele S, Esposito G, Nardinocchi P, Teresi L, Madeo A, Schiariti M, Varano V, Puddu PE. Homeostatic left heart integration and disintegration links atrio-ventricular covariation’s dyshomeostasis in hypertrophic cardiomyopathy. Sci Rep. 2017;7:6257.
Neaux D, Sansalone G, Ledogar JA, Ledogar SH, Luk TH, Wroe S. Basicranium and face: assessing the impact of morphological integration on primate evolution. J Hum Evol. 2018;118:43–55.
Bookstein FL. Integration, disintegration, and self-similarity: characterizing the scales of shape variation in landmark data. Evol Biol. 2015;42:395–426.
Collyer ML, Adams DC. RRPP: an R package for fitting linear models to high-dimensional data using residual randomization. Methods Ecol Evol. 2018;9:1772–9.
Adams DC. A generalized K statistic for estimating phylogenetic signal from shape and other high-dimensional multivariate data. Syst Biol. 2014;63:85–97.
Adams DC, Felice RN. Assessing trait covariation and morphological integration on phylogenies using evolutionary covariance matrices. PLoS One. 2014;9:e94335.
Ziegler R. Order Insectivora. In: Rössner GE, Heissig K, editors. The Miocene Land Mammals of Europe. München: Verlag Dr. Friedrich Pfeil; 1999. p. 53–74.
Gunnell GF, Bown TM, Hutchinson JH, Bloch JI. Lipotyphla. In: Janis CM, Gunnell GF, Uhen MD, editors. Evolution of tertiary mammals of North America, Vol 2: small mammals, Xenarthrans, and marine mammals. Cambridge University Press: Cambridge; 2008. p. 89–125.
Hildebrand M, Goslow GE, Hildebrand V. Analysis of vertebrate structure. New York: John Wiley & Sons; 1995.
Ziegler AC. Dental homologies and possible relationships of recent Talpidae. J Mammal. 1971;52:50–68.
Funmilayo O. Distribution and abundance of moles (Talpa europaea L.) in relation to physical habitat and food supply. Oecologia. 1977;30:277–83.
Funmilayo O. Food consumption, preferences and storage in the mole. Acta Theriol. 1979;24:79–89.
Takahashi M, Arai T, Ohdachi S. Stomach contents of Mogera wogura wogura (Temminck) collected in the Hanno City, Saitama prefecture. Japan Part I Bulletin of Kawagoe Senior High School. 1982;19:24–9.
Gorman ML, Stone RD. The natural history of the mole. London: Christopher Helm Ltd; 1990.
Dalquest WW, Orcutt DR. The biology of the least shrew-mole, Neurotrichus gibbsii minor. Am Midland Nat. 1942;27:387–401.
Catania KC, Remple FE. Asymptotic prey profitability drives star-nosed moles to the foraging speed limit. Nature. 2005;433:519.
Oparina OS, Filinova EI, Sonina EE, Malinina YA, Oparin ML. Current status of the Russian desman habitats in small rivers of the Don River basin in Saratov oblast and the abundance of this species. Biol Bull. 2013;40:54–61.
Dumont M, Wall CE, Botton-Divet L, Goswami A, Peigné S, Fabre AC. Do functional demands associated with locomotor habitat, diet, and activity pattern drive skull shape evolution in musteloid carnivorans? Biol J Linn Soc. 2016;117:858–78.
Arregoitia LD, Fisher DO, Schweizer M. Morphology captures diet and locomotor types in rodents. R Soc Open Sci. 2017;4:160957.
Botton-Divet L, Houssaye A, Herrel A, Fabre AC, Cornette R. Swimmers, diggers, climbers and more, a study of integration across the mustelids’ locomotor apparatus (Carnivora: Mustelidae). Evol Biol. 2018;45:182–95.
Hansen TF. Is modularity necessary for evolvability?: remarks on the relationship between pleiotropy and evolvability. Biosystems. 2003;69:83–94.
Goswami A, Polly PD. The influence of modularity on cranial morphological disparity in Carnivora and Primates (Mammalia). PLoS One. 2010;5:e9517.
Marroig G, Shirai LT, Porto A, de Oliveira FB, De Conto V. The evolution of modularity in the mammalian skull II: evolutionary consequences. Evol Biol. 2009;36:36–48.
Bright JA, Marugán-Lobón J, Cobb SN, Rayfield EJ. The shapes of bird beaks are highly controlled by nondietary factors. PNAS. 2016;201602683.
Bickelmann C, Mitgutsch C, Richardson MK, Jiménez R, de Bakker MA, Sánchez-Villagra MR. Transcriptional heterochrony in talpid mole autopods. EvoDevo. 2012;3:16.
Bickelmann C, Jiménez R, Richardson MK, Sánchez-Villagra MR. Humerus development in moles (Talpidae, Mammalia). Acta Zool. 2014;95:283–9.
Renaud S, Michaux J, Schmidt DN, Aguilar JP, Mein P, Auffray JC. Morphological evolution, ecological diversification and climate change in rodents. Proc R Soc B Biol Sci. 2005;272:609–17.
We are particularly grateful to Prof. Barbara Rzebik-Kowalska of ISEZ-PAN, Krakow, Poland, and Dr. Shin-Ichiro Kawada of the Tsukuba Natural History Museum, Tsukuba, Ibaraki, Japan, for their helpful suggestions during the manuscript preparation. We wish to thank Dr. Christiane Bastos-Silvera of the Museu de Historia Natural, Lisboa, Portugal; Dr. Paula Jenkins, Dr. Roberto Portella and Dr. Emma Bernard of NHM, London, UK; Germany; Dr. Gertrud Roessner of BSPG, Munich, Germany; Dr. Ursula Goelich of the Wien Natural History Museum, Wien, Austria; Dr. Jim Dines of LACM, Los Angeles, USA; Dr. Patricia Holroyd of UCMP, Berkeley, USA; and Dr. Paolo Agnelli of “La Specola” Museo di Storia Naturale di.
Firenze, Italy. They each allowed us to visit their collections and made our visits comfortable and pleasant. We are grateful to three anonymous reviewers for their comments, which greatly improved the quality of the manuscript.
G. Sansalone received support from the SYNTHESYS Project (http://www.synthesys.info), which is financed by the European Community Research Infrastructure Action under the FP7 “Capacities” Program (GB-TAF-2095 and AT-TAF-3415). The funding bodies had no role in study design, analysis and interpretation, or writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.