De Oliveira Etal 2019
De Oliveira Etal 2019
net/publication/329334032
CITATIONS READS
10 500
5 authors, including:
Virginia Abdala
National University of Tucuman
163 PUBLICATIONS 2,685 CITATIONS
SEE PROFILE
All content following this page was uploaded by Débora Lina Moreno Azócar on 07 May 2020.
Article
Abstract
The interaction between organisms and their environment is central in functional morphology.
Differences in habitat usage may imply divergent morphology of locomotor systems; thus,
detecting which morphological traits are conservative across lineages and which ones vary under
environmental pressure is important in evolutionary studies. We studied internal and external
morphology in 28 species of Neotropical anurans. Our aim was to determine if internal morphology
(muscle and tendons) shows lower phylogenetic signal than external morphology. In addition, we
wanted to know if morphology varies in relation to the habitat use and if there are different functional
groups. We found differences in the degree of phylogenetic signal on the groups of traits.
Interestingly, postaxial regions of the forelimb are evolutionarily more labile than the preaxial
regions. Phylomorphospace plots show that arboreal (jumpers and graspers) and swimmer frogs
cluster based on length of fingers and the lack of sesamoid, also reflected by the use of habitat.
These functional clusters are also related to phylogeny. Sesamoid and flexor plate dimensions to-
gether with digit tendons showed to be important to discriminate functional groups as well as use of
habitat classification. Our results allow us to identify a “grasping syndrome” in the hand of these
frogs, where palmar sesamoid and flexor plate are absent and a third metacarpal with a bony knob
are typical. Thus, a lighter skeleton, long fingers and a prensile hand may be key for arboreality.
Key words: comparative analyses, ecomorphology, neotropical frogs, function, habits, phylogeny
How an organism interacts with its environments has important for natural selection to act upon (Arnold 1983; Kingsolver and
implications for the selective forces shaping the phenotypes of the Huey 2003). In this sense, the relationship between morphology,
species. The potential association between morphology and func- locomotor performance and habitat use is one of the most studied
tional performance at different levels (individual, population and aspects of evolutionary phenotypic variation (Losos and Sinervo
species) may yield differences in fitness providing the raw material 1989; Losos 1990a, 1990b; Bonine and Garland 1999; Van Damme
and Vanhooydonck 2002; Goodman et al. 2008; Irschick et al. Fontanarrosa and Abdala 2014, 2016), we also included this struc-
2008). It was observed that differences in habitat occupation under- ture in our study. We test whether the external and muscle-
lie divergent evolution of the morphology of locomotor systems in tendinous morphologies of anuran forelimbs differ among taxa that
numerous taxa (Irschick et al. 2005; Calsbeek and Irschick 2007). exhibit different locomotor modes or habitat use. We also test
Among the studies on this issue, research on reptiles certainly out- whether such patterns can be explained solely by the phylogenetic
number those in other groups; for example, many studies have inves- relationships between the species examined or by another underly-
tigated the relationship between limb morphology and habitat use in ing factor. It should be noted that the muscle-tendinous structures
Anolis, liolaemids and tropidurine lizards (Losos 1990a, 1990b, have been seldom considered in ecomorphological studies of verte-
1990c; Kohlsdorf et al. 2001; Irschick et al. 2005; Calsbeek and brates (Abdala et al. 2008, 2014; Tulli et al. 2012b; Fabrezi et al.
Irschick 2007; Grizante et al. 2010; Tulli et al. 2009, 2011, 2012a, 2014; Carrizo et al. 2014; Fratani et al. 2018). We hypothesize that:
2012b). However, these associations are not always clear (Tulli a) internal morphology (muscle and tendons) will show a higher
et al. 2012b, 2016). In a broad evolutionary context, convergent phylogenetic signal pattern than external morphology as was shown
evolution of traits in species from different regions can lead to simi- in previous studies on other tetrapod taxa (Tulli et al. 2012a;
lar body shape among them, although they could be phylogenetical- Carrizo et al. 2014); b) morphological traits of terrestrial species
length, maximum width of muscles was also recorded to obtain an Mean values of all morphological variables (Supporting
estimate of the morphometric variation of each muscle as a whole to Information Supplementary Material) were log10 transformed for
allow testing if there exists correlation between aponeurosis and ten- further analyses. We then tested the data for phylogenetic signal
don dimensions that are associated with these muscles. (i.e., the consistency in trait values with the phylogeny); for this pur-
The different categories of habitat use and locomotor modes pose, we used the Pagel’s k value estimator for each variable (Pagel
were considered after Cei (1980), Wells (2007), Jorgensen and 1999) that varies between 0 and 1, where 0 means no phylogenetic
Reilly (2013) and personal observations. Therefore, we obtained 4 signal and 1 means that a variable is highly dependent of the phylo-
habitat use categories (fossorial, terrestrial, arboreal and aquatic), genetic structure. In addition, a randomization analysis was per-
and 6 locomotor mode categories (terrestrial walker, terrestrial formed for each variable to determine the probability of the
jumper, terrestrial burrower, arboreal jumper, arboreal grasper- estimator.
walker and swimmer; Figure 1). As the morphological variables studied here showed to be highly
Because species have a shared history (Figure 1), data of traits correlated with body size (r2 > 0.8), and to consider the phylogenet-
corresponding to these species cannot be considered as independent ic context of this study, we proceeded with phylogenetic size correc-
data-points (Harvey and Pagel 1991) and phylogenetically informed tion (Revell 2009). Through this, we obtained the residuals of each
statistical analyses are required. For this purpose, we used a compos- target morphological trait (least squares regression analysis trait vs.
ite tree of the phylogenetic relationship based on Pyron and Wiens body size - snout-vent length, SVL) while simultaneously controlling
(2011) tree and because branch lengths (BL) were not available, we for phylogenetic independence. Then, we used the residuals in subse-
used arbitrary branch lengths (e.g., all BL equal to one, BL trans- quent analyses as size-corrected morphological measurements.
formed with Grafen method and BL transformed with Pagel We studied 23 different morphological variables (Supplementary
method). We tested the adequacy of each one of these 3 arbitrary Material, Table S2). Since this amount of information is difficult to
branch lengths following Garland et al. (1992), by plotting the abso- analyze, we reduced the number of variables by running a phylogen-
lute value of each standardized independent contrast versus its etically based principal component analysis (phyl.pca). We collected
standard deviation. We used Mesquite v2.74 (Maddison and the scores of each principal component from the phylogenetically
Maddison 2015) and the PDAP PDTREE v1.15 modules (Midford informed PCA. We considered the first 4 principal components that
et al. 2009) and searched for significant differences. After plotting attained for 84% of the accumulated variance of the total analyses.
trees and BL for all traits, we found that the tree with BL equal to We also plotted the phylomorphospace and the contribution of all
one was the more adequate arrangement (BL ¼ 1 showed 4 out of morphological variables; for this descriptive purpose, we used the
26 significant plots, while Pagel BL transformation showed 8/26 and command phylomorphospace from the program Phytools (Revell
Grafen 11/26 significant relationships). 2012).
602 Current Zoology, 2019, Vol. 65, No. 5
We tested if internal or external variables showed variability in modes; these analyses were conducted through a phylogenetically
relation to the 6 functional groups considered here (fossorial informed ANOVA too. All these analyses were performed including
walkers, terrestrial walkers, terrestrial jumpers, arboreal jumpers, phylogenetic information through the following functions; phylosig,
arboreal walker or graspers and aquatic swimmers) as well as for phyl.resid, phyl.pca, phylomorphospace, phylANOVA of the
the type of habitats these frogs use (e.g., burrows, terrestrial, arbor- Phytools package version 6.0 (Revell 2012) in the open access envir-
eal and aquatic). In this case we used phylogenetic MANOVA by onment R (Version 3.4.0, R Core team 2017).
running Geiger program (Harmon et al. 2008), we concatenated the
internal or external variables in relation to the different factors
(functional groups or habitat types) as a formula to compute the Results
analyses. Finally, for testing if there were differences in each one of In all dissected specimens the m. flexor digitorum communis origi-
the different morphological traits (in total 23 external and internal nates from the humerus distal condyle, through a wide and short
traits) among species as a function of their habitat use (aquatic, ter- tendon from the aponeurosis covering the elbow. The m. flexor
restrial, arboreal or fossorial) a phylogenetically based analysis of carpi ulnaris originates on the medial epicondylus of the humerus,
variance was performed (phylANOVA of Phytools, Revell 2012) with some fibers connected to the m. flexor digitorum communis.
Figure 2. External characters. (A) body and arm external measurements, (B) manus external measurements (C) muscle and tendons showing sesamoid and flex-
or plate, (D) tendons and muscles where no flexor plate or sesamoid are present (drawings courtesy S. De Oliveira Lagôa).
De Oliveira-Lagôa et al. Muscle tendinous structure, locomotor modes and habitat use 603
Odontophrynidae, and most of Microhylidae present a flexor plate showed high values of phylogenetic signal estimators (0.999,
with an embedded palmar sesamoid, over which a palmar aponeur- Table 1).
osis is distinguishable (Figure 2C); conversely, hylids, Phylogenetic PCA showed that the following variables contrib-
Phyllomedusidae and some microhylids lack a palmar sesamoid uted most (based on higher absolute eigenvalues) to the principal
(Figure 2D). From the flexor plate originate the flexor tendons of components functions (Table 2): sesamoid length, sesamoid width,
digits III, IV, and V that insert onto the basal portion of the distal digit IV and V tendon length and width (respectively), flexor carpi
phalanx of each digit (Figure 2D). ulnaris tendon length and width (internal morphology variables;
In the case of bufonids, a tendinous sheet between the mm. flex- Table 2); similarly, forearm width and hand length show high load-
or digitorum longus and mm. opponens of digit II are distinguish- ings in the PCA (Table 2). In addition, phylomorphospace plots con-
able. In leptodactylids, telmatobiids, and odontophrynids, there is sidering external and internal morphology characters show that all
an aponeurosis between the flexor plate and the muscles of digit II. arboreal (jumpers and graspers) and one swimmer species cluster
In Leptodactylus fuscus, the flexor tendon of digit II is visible and according to the length of finger tendons and sesamoid (Figure 3A),
originates from the flexor plate inserting into the basal portion of which is reflected by the use of habitat where aquatic (partially) and
the distal phalanx. The tendon is covered by fascias and the mm. arboreal species clearly differ from terrestrial and burrower species
Locomotor mode
df Wilk’s F P
Habitat use
df Wilk’s F P
Discussion
Our data show an important pattern of phylogenetic signal: half of the
internal morphology traits show high phylogenetic signal estimator val-
ues (Table 1). The effect of the phylogeny on the internal morphology
was also verified in other tetrapod groups, such as lizards and anurans
(Tulli et al. 2012a; Abdala et al. 2014; Fratani et al. 2018). Interestingly,
digit tendons length was additionally an important trait that segregated
Figure 3. Phylomorphospace plots of (a) functional groups (terrestrial walker, arboreal and swimming species showing that this morphology–lifestyle
jumper, and burrower; swimmer and arboreal walker and jumper); (b) habitat relationship probably arose from features clustered in phylogenetic
use (terrestrial; fossorial, swimmer, and arboreal) and (c) family level groups with the same habitat use (Vanhooydonck and Van Damme
(Microhylidae, Hylidae, Leptodactylidae, Odontophynidae, Bufonidae, 1999). This niche conservatism (Wiens and Graham 2005) was also ap-
Telmatobiidae) among the 28 anuran species. Axes correspond to PC1 35.6%
parent among microhylids that exhibit a strong conservatism in ecology,
of variance (increase of finger tendon [III, IV, and V] length, decrease of ses-
morphology and performance (although, not in body size). This is espe-
amoid and flexor plate size) and PC2 28.0% of variance (increase of flexor
carpi ulnaris—FCU muscle and sesamoid length and decrease of FCU tendon cially remarkable, given that these species have been separated for more
size). than 65 million years (Moen et al. 2013). Taken together these data
De Oliveira-Lagôa et al. Muscle tendinous structure, locomotor modes and habitat use 605
Table 4. Phylogenetic ANOVA of external and internal variables as a function of locomotor mode (Tw ¼ terrestrial walker, Tj ¼ terrestrial
jumper, Bw ¼ fossorial walker, Aj ¼ arboreal jumper, Aw ¼ arboreal grasper or walker, and S ¼ aquatic swimmer) and habitat use (F ¼ fos-
sorial, T ¼ terrestrial, Ar ¼ arboreal, Aq ¼ aquatic)
F P posthoc F P posthoc
External Measures
Arm length 1.138 0.429 1.244 0.768
Forearm length 4.081 0.549 0.192 0.98
Palm length 4.084 0.545 0.974 0.808
Palm width 0.874 0.941 0.735 0.868
Hand length 1.631 0.848 2.851 0.549
Internal Measures
Sesamoid length 21.012 0.009 TwTjBwS AjAw 31.285 0.005 BT Aq Ar
Boldface denotes significant differences and post hoc (Holm method) shows the homogeneous groups formed between groups (when underlined together).
support our first hypothesis, since internal morphology traits tended to may play in locomotion (Wang et al. 2014). Previous studies have
have higher phylogenetic signal than external ones. shown that toads use their forelimbs almost exclusively to decelerate
Among external morphology traits, only arm length exhibited a and control the body at the moment of impact (Gillis et al. 2014),
phylogenetic pattern indicating that this trait is somehow more our results suggest that the postaxial region of the forelimbs could
unique to each lineage among the species studied here. Traits exhib- exert a differential role in landing deserving special consideration,
iting strong phylogenetic signal suggest gradual changes over time particularly in jumping species (all leptodactylid frogs and some
(e.g., Brownian Motion model of evolution), and traits with no hylids).
phylogenetic signal may be highly labile on the time scale of phyl- Our phylogenetic informed MANOVA and ANOVA show that
ogeny, or may not change at all (Revell 2008). some of the selected morphological variables are clearly related to
Other traits that showed low phylogenetic signal were the width the arboreal locomotion mode and may be characterized by the mor-
and length of the tendon of the flexor carpi ulnaris that also played phological pattern of the forelimb, such as lack of sesamoid bone
an important role in our phylomorphospace figures. These traits and flexor plate. Forearm structures in arboreal anurans show that
may segregate the jumper species from the rest (particularly among they present in their manus the same pervasive trend exhibited for
the terrestrial species), but we need more detailed information to all lineages of grasping tetrapods (Sustaita et al. 2013; Manzano
corroborate this assumption. The low phylogenetic signal in the et al. in press). A vertical position in terrestrial environments
flexor tendons width of digits IV and V, and flexor carpi ulnaris requires that the animals stay in balance by pulling their fore-feet to-
muscle suggests that the postaxial regions of the forelimb are evolu- ward the substrate, while their hind feet push on the substrate (Tulli
tionarily more prompted to vary than the preaxial regions. These et al. 2011). This ability allows organisms to balance over a vertical
results are interesting as the postaxial region corresponds to the pri- substrate and requires that the organisms maintain their center-of-
mary axis that passes through digit IV in early limb development mass aligned with their support. This can be achieved by developing
(Fabrezi et al. 2007) and it would be expected a more stable config- an adhesive or suction force between their body and the contact sur-
uration to this crucial limb zone. However, it should be considered face (Cartmill 1985) or grabbing a narrow stem of the surface, as
that those variables related to length of the postaxial (ulnar) struc- most arboreal anurans and other grasping tetrapods do (Sustaita
tures show a strong phylogenetic pattern, decoupled thus of their et al. 2013; Manzano et al. in press). These morphological descrip-
width, which exhibit a lower phylogenetic signal, contrasting with tors of the arboreal habitat use in the studied frogs allow us to up-
the observations of Herrel et al. (2014) for hind limb kinematics of hold our second and third predictions. Previous studies across many
Xenopus tropicalis. These authors found that hind limb length in organisms have shown a strong relationship between ecology (e.g.,
females were good descriptors of the peak resultant jumping force, microhabitat use and diet) and morphology (see reviews in
compared with longer iliums and the length of the longest toe in Wainwright 2007). Far fewer have considered performance (e.g.,
males. This probably indicates the different roles fore and hind limb Losos 1990a, 1990b, 1990c; Herrel et al. 2014; Tulli et al. 2011,
606 Current Zoology, 2019, Vol. 65, No. 5
2012a), which forms the link between ecology and morphology narrow branches or twigs (Manzano et al. in press). Aquatic frogs
(Wainwright 1991). In addition, Moen et al. (2013) found that frog such as Xenopus with a large palmar sesamoid suggest that their
species using the same microhabitat had both similar morphology presence is owed more to phylogeny than to ecological reasons.
and performance, regardless of geographical location. However, recent studies suggest that Xenopus, a genus lacking
A recent study on morphology and habitat use also found that grasping syndrome, exhibits also considerable manual abilities,
frogs that live in the bush and swim and jump, exhibit modifications including an intermediate or scissor grip (Anzeraey et al. 2017). This
in the postcranial skeleton such as broad proximal sacral diapophy- report defied the pervasive link between arboreality and manual
ses, broad vertebral bodies, and longer urostyles that relate their grasping and requires a new perspective to explain the genesis of
morphology and habits (Soliz and Ponssa 2016). These authors also grasping among anurans (Manzano et al. in press). In addition, it
suggest that a slender body would be better adjusted to face any has to be noted that the palmar sesamoid may be ambiguously pre-
challenging locomotor modes, which is typical of hylid frogs that sent or absent (Ponssa et al. 2010); consequently, our results should
possess slender hands too. In addition, arboreal and walker species be taken with caution because of a possible taxon sampling prob-
have smaller and shorter vertebral bodies in their column, weaker lem, probably as a consequence of choosing species that may coexist
coracoids and clavicles compared with jumper and swimmer species as they occur in Chaco and Monte biogeographic provinces
Burton TC, 1998. Are the distal extensor muscles of the fingers of anurans an Grizante MB, Navas CA, Garland T Jr, Kohlsdorf T, 2010. Morphological
adaptation to arboreality. J Herpetol 32:611–617. evolution in Tropidurinae squamates: an integrated view along a continuum
Burton TC, 2001. Variation in the foot muscles of frogs of the family of ecological settings. J Evol Biol 3:98–111.
Myobatrachidae. Australian J Zool 49:539–559. Haines RW, 1950. The flexor muscles of the forearm and hand in lizards and
Burton TC, 2004. Muscles of the pes of hylid frogs. J Morphol 260:209–233. mammals. J Anat 84:13–29.
Cabrera A, Willink A, 1980. Biogeografı́a de America Latina. Secretarı́a Harmon LJ, Weir JT, Brock CD, Glor RE, Challenger W, 2008. GEIGER:
General OEA, Serie Biológica, Monografı́a. investigating evolutionary radiations. Bioinformatics 24:129–131.
Calsbeek R, Irschick DJ, 2007. The quick and the dead: correlational selection Harvey PH, Pagel MD, 1991. The Comparative Method in Evolutionary
on morphology, performance, and habitat use in island lizards. Evolution Biology. Oxford: Oxford University Press.
61:2493–2503. Herrel A, Vasilopoulou-Kampitsi M, Bonneaud C, 2014. Jumping perform-
Carrizo LV, Tulli MJ, Dos Santos DA, Abdala V, 2014. Interplay between ance in the highly aquatic frog Xenopus tropicalis: sex-specific relationships
postcranial morphology and locomotor types in Neotropical sigmodontine between morphology and performance. PeerJ 2: e661.
rodents. J Anat 224:469–481. Hoyos JM, Mantilla C, Galindo D, Salgar L, 2014. Phylogenetic analysis with-
Cartmill M, 1985. Climbing. In: Hildebrand M, Bramble DM, Liem KF, in the Pristimantis unistrigatus (Anura, Craugastoridae) group based on
Wake DB, editors. Functional Vertebrate Morphology. Cambridge: Belknap morphological characters. Caldasia 36:107–124.
Pyron RA, Wiens JJ, 2011. A large-scale phylogeny of Amphibia including Tulli MJ, Abdala V, Cruz FB, 2011. Relationships among morphology, cling-
over 2800 species, and a revised classification of extant frogs, salamanders, ing performance and habitat use in Liolaemini lizards. J Evol Biol 24:
and caecilians. Mol Phyl Evol 61:543–583. 843–855.
R Core Team. 2017. R: A Language and Environment for Statistical Tulli MJ, Abdala V, Cruz FB, 2012a. Effects of different substrates on the
Computing. Vienna, Austria: R Foundation for Statistical Computing. sprint performance of lizards. J Exp Biol 215:774–784.
Available from: http://www.R-project.org/ Tulli MJ, Herrel A, Vanhooydonck B, Abdala V, 2012b. Is phylogeny driving
Rand AS, 1952. Jumping ability of certain anurans, with notes on endurance. tendon length in lizards? Acta Zool 93:319–329.
Copeia 1952:15–20. Tulli MJ, Cruz FB, Kohlsdorf T, Abdala V, 2016. When a general morphology
Reilly SM, Jorgensen ME, 2011. The evolution of jumping in frogs: morpho- allows many habitat uses. Integr Zool 11:473–489.
logical evidence for the basal anuran locomotor condition and radiation of Van Damme R, Vanhooydonck B, 2002. Speed versus maneuverability: associ-
locomotor system in crown group anurans. J Morphol 272:149–168. ation between vertebral number and hábitat structure in lacertid lizards.
Revell LJ, 2008. On the analysis of evolutionary change along single branches J Zool 258:327–334.
in a phylogeny. Am Nat 172:140–147. Vanhooydonck B, Van Damme R, 1999. Evolutionary relationships between
Revell LJ, 2009. Size-correction and principal components for interspecific body shape and habitat use in lacertid lizards. Evol Ecol Res 1:785–805.
comparative studies. Evolution 63:3258–3268. Vidal-Garcı́a M, Byrne PG, Roberts JD, Keogh JS, 2014. The role of phyl-