Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Are Diet Preferences Associated to Skulls Shape Diversification in Xenodontine Snakes?

  • Julia Klaczko ,

    jklaczko@gmail.com

    Affiliation Departamento de Biologia Animal, Instituto de Biologia, Universidade Estadual de Campinas, UNICAMP, Campinas, SP, Brazil

  • Emma Sherratt,

    Affiliation School of Environmental and Rural Science, University of New England, Armidale, New South Wales, Australia

  • Eleonore Z. F. Setz

    Affiliation Departamento de Biologia Animal, Instituto de Biologia, Universidade Estadual de Campinas, UNICAMP, Campinas, SP, Brazil

Abstract

Snakes are a highly successful group of vertebrates, within great diversity in habitat, diet, and morphology. The unique adaptations for the snake skull for ingesting large prey in more primitive macrostomatan snakes have been well documented. However, subsequent diversification in snake cranial shape in relation to dietary specializations has rarely been studied (e.g. piscivory in natricine snakes). Here we examine a large clade of snakes with a broad spectrum of diet preferences to test if diet preferences are correlated to shape variation in snake skulls. Specifically, we studied the Xenodontinae snakes, a speciose clade of South American snakes, which show a broad range of diets including invertebrates, amphibians, snakes, lizards, and small mammals. We characterized the skull morphology of 19 species of xenodontine snakes using geometric morphometric techniques, and used phylogenetic comparative methods to test the association between diet and skull morphology. Using phylogenetic partial least squares analysis (PPLS) we show that skull morphology is highly associated with diet preferences in xenodontine snakes.

Introduction

The origin and diversification of morphology are topics of great interest with the field of evolutionary biology, and the adaptation of organismal form to ecological conditions has been attributed as a primary driving force of morphological diversification [1]. Classic support for the hypothesis of adaptation by natural selection is evolutionary convergence. Evolutionary convergence occurs when similar phenotypes evolve in phylogenetically independent taxa as a response to similar ecological conditions [2, 3]. Due to the importance of the skull, and its direct link to an animal’s fitness, it is presumed that skull morphology is under considerably strong selection pressure [4]. Among the many ecological functions of the skull, feeding is one of the most essential; and presumably, diet can influence the skull [5, 6]. Many studies within vertebrates have corroborated this idea, showing strong correlations between diet and skull morphology [6, 7, 8]. Stayton [9] studied lizard skulls shape evolution across 17 families, using geometric morphometric tools. He showed morphological convergent evolution among lizards with similar diets. However, these patterns were secondary over the more general phylogenetic pattern of lizard skull diversity.

Snakes are a highly successful group of vertebrates, with large diversity in habits, environments, diet, and morphology [10]. They comprise more than 3000 species, widely distributed across temperate and tropical regions, with diverse dietary preferences [11]. Those historical shifts probably resulted in adaptive radiations that contributed to the high diversity of species observed today [12, 13].

Snakes, particularly macrostomatans, can feed on very large prey, much larger that the size of their heads. Snakes with macrostoman condition possess highly kinetic skulls, which allow for the transport of the entire prey through their oral cavities using ratcheting motions of the lower jaws [11, 14, 12, 15]. Genomic phylogenies suggest that this capacity have evolved twice, in tropidophiids, basal alethinophidians, and the families Bolyeriidae, Tropidophiidae, Boidae, and Pythonidae, plus the "advanced snakes" that include Acrochordidae and Colubroidea grouped as Caenophidia [16]. Several studies have suggested association between specific types of diet and the morphology of snakes [17, 18]. However, to our knowledge, no study has used geometric morphometric tools and phylogenetic comparative methods to explicitly test the association between diet and skull morphology in snakes.

Xenodontines are a speciose clade of South American snakes that include 49 genera and approximately 330 species [19]. Molecular phylogenetic hypotheses corroborate the monophyly of the group [20], and are well known for their great diversity in morphology and ecological features (Fig 1, [21]). In particular, xenodontine snakes show a broad range of diets including invertebrates, amphibians, lizards, snakes and small mammals [22]; and so are a good representative of the diet preference diversity found among modern macrostomatan snakes. The high level of morphological and ecological diversity makes xenodontine snakes an ideal group to study mechanisms that promote lineage diversification and evolutionary radiations in macrostomatan snakes [21].

thumbnail
Fig 1. Example of morphological variation in Xenodontinae snake skulls.

Dorsal view (A) and lateral view (B) of Helicops angulatus UFMT 7818, an eater of fish, lizards, and anurans; Dorsal view (C) and lateral view (D) of Philodryas aestivus IBSP 6771, an eater of mammals and anurans; Dorsal view (E) and lateral view (F) of Phimophis guerini IBSP 66406, an eater of lizards and mammals; Dorsal view (G) and lateral view (H) of Phalotris mertensi ZUEC 486, an eater of elongate vertebrates (caecilians and amphisbaenians).

https://doi.org/10.1371/journal.pone.0148375.g001

Here we describe the variation in skull shape of 19 species of xenodontine snakes using geometric morphometric techniques [23; 24]. Using dietary preferences compiled from the literature with phylogenetic comparative methods we addressed the following question: is the variation of skull shape in xenodontine snakes influenced only by the speciation process, or has the diet fostered the diversification of skull shape?

Material and Methods

Data

We examined skull morphology in 19 species of South American xenodontine snakes (range 1–5 specimens per species, mean of 4), representing the diversity of diet preferences present in xenodontine snakes. The analyzed specimens were from the following museums: Museu Nacional, Rio de Janeiro, Brazil (MNRJ), Instituto Butantan, São Paulo, Brazil (IBSP), Universidade Federal do Mato Grosso, Cuiabá, Brazil (UFMT), Museu de Zoologia “Prof. Adão José Cardoso”, Unicamp, Campinas, São Paulo, Brazil (ZUEC). The complete list of the analyzed specimens can be found in the Supporting Information, Table A in S1 File. The skulls were dissected from museum specimens and skeletonized by hand. We obtained digital images of the skulls in dorsal and lateral views, using a Canon PowerShot S5 SI digital camera. Images were standardized for skull position, camera lens plane position, and distance between camera lens and skull.

One of us (JK) digitized a set of 20 landmarks in each view (dorsal and lateral) of the skull (Fig 2, for landmarks definitions see list on Table B in S1 File) using tpsDig 2.12 [25]. The landmarks were digitized only on one side of the skull, since our focus was not on the asymmetric component of shape [26]. The landmark data were aligned using Generalized Procrustes Analysis [27; 23]. Then the average shape for each species was calculated and used for subsequent analyses.

thumbnail
Fig 2. Cranial landmarks.

Cranial landmarks recorded from South American Xenodontinae snakes. Dorsal view (A), anatomical wire frame of dorsal view (B), lateral view (C), and anatomical wire frame of lateral view (D).

https://doi.org/10.1371/journal.pone.0148375.g002

Principal Component Analysis was performed on the average shape for all species Procrustes shape coordinates to visualize the variation among species. We obtained 18 principal axes (PCs) for each view that describes the shape variation. The differences in shape described by each principal axis (PC) were summarized using thin-plate spline deformation grids [28; 29]. The analyses were performed using Geomorph package v.2.1 [29; 30] for R software v.3.1.1 [31], and tpsRelw 1.46 software [25].

Diet preferences were obtained from literature, and are presented as the proportion of each item on the diet (Table 1).

thumbnail
Table 1. Proportion of prey items recorded for xenodontinae snakes; an = anuran amphibians; bi = birds; ca = caecilids or amphisbenids; eg = lizard eggs; ew = earthworms; fi = fishes; li = lizards; ma = mammals; mo = mollusks; sn = snakes.

https://doi.org/10.1371/journal.pone.0148375.t001

Phylogeny and Comparative Approach

The phylogenetic relationships between the xenodontine species were reconstructed using Dipsas indica as the outgroup. We used sequences of mitochondrial genes 12S and 16S rDNA, and nuclear gene Oocyte maturation factor-like (c-mos). All sequences were obtained from GenBank, and aligned using ClustalW [44] (List of used accession numbers–Table C in S1 File). jModelTest 0.1 [45] was employed to determine the most appropriate model of sequence evolution for each analyzed gene, estimated under the Bayesian Information Criterion (BIC). The best fit model for the genes 12S and 16S was TPM3uf + G, and for the gene c-mos was TPM2uf. The combined molecular data set was analyzed under the Bayesian Inference method, implementing the selected optimal sequence evolution model for each gene. The phylogenetic analysis was performed using MrBayes v3.1.2 [46]. Two runs were performed, with default heating for each of the four chains and sampling every 100 generations for 20,000,000 generations. The convergence was verified using Tracer v1.4 [47]; and the first 5000 topologies (25% of the sampled topologies) were discarded as burn-in.

The phylogenetic signal skull shape (the Procrustes shape coordinates) was evaluated using the multivariate K statistic (see [48]) in the R package Geomorph v.2.1. The K statistic can indicate little or no phylogenetic signal (K << 1), or phylogenetic signal (K ≈ 1), or greater than the expectation under a Brownian motion random-walk model of evolution (K >> 1). To test if the data contain a signal of phylogeny (i.e. K > 0), we randomly permuted the order of species on the tree 1000 times and recalculated K for each permutation. We then compared the observed K value to this null distribution to assess significance.

To visualize the evolutionary history of the analyzed of skull shape, the consensus phylogenetic tree was projected onto the shape space. We used the squared-change parsimony method implemented in MorphoJ [49], which is equivalent to maximum likelihood methods when branch lengths are present.

To characterize the evolutionary patterns of covariation between skull shape variation and diet preferences we used phylogenetic two-block partial least squares analysis (PPLS; [50]). Partial least squares is a statistical procedure that quantifies the degree of covariation between sets of variables, based on a singular value decomposition of the overall trait covariance matrix [51]. The phylogenetic PLS uses instead the evolutionary covariance matrix and assumes the expected lack of independence among samples as a result of phylogenetic relationship [51]. Since our consensus phylogenetic tree shows some polytomies, the PPLS analysis was conducted on 200 resolved trees randomly selected from MrBayes output. The significance of the model was assessed employing a permutation test. The PPLS was performed in the Geomorph package v.2.1.

Results

The principal components (PC) analysis of the skull shapes reveals that most of the shape variation is contained in few dimensions; in the dorsal view the first four PCs explained 91.5%, while in the lateral view the first four PCs explained 88.3% of the total shape variation.

In the dorsal view, PC1 explained 56.4% of the variation, and the species with negative values show an increase of skull length, especially of the parietal bone, and a relative decrease in orbit size; while positive values indicate a wider and shorter skull (Fig 3). On the other hand, in the lateral view the PC1 explained 54.4% of the variation; positive values of the PC1, on the right side of the graph, show a more robust and shortened skull, especially the maxilla bone, while species with negative values of PC1, show a more elongated jaw (Fig 4).

thumbnail
Fig 3. The phylomorphospace of dorsal cranial shape.

(A) Projection of the phylogenetic tree into the dorsal view PC morphospace. (B) Estimated changes in dorsal view shape are shown as deformations from the mean shape along the first and second principal components.

https://doi.org/10.1371/journal.pone.0148375.g003

thumbnail
Fig 4. The phylomorphospace of lateral cranial shape.

(A) Projection of the phylogenetic tree into the lateral view PC morphospace. (B) Estimated changes in lateral view shape are shown as deformations from the mean shape along the first and second principal components.

https://doi.org/10.1371/journal.pone.0148375.g004

The molecular phylogeny based on Bayesian Inference analysis was generally consistent with previous phylogenetic hypotheses, we found support for the monophyly of all the xenodontine tribes [19; 20], (Figure A in S1 File). Our phylogeny showed the clade Erythrolamprus and Lystrophis associated to the clade that contain Siphlophis, Oxyrhopus, Boiruna and Phimophis while Vidal et al. [19] showed Erythrolamprus and Lystrophis as a more inclusive clade sister group of the clade containing Helicops, Tomodon, Gomesophis, Siphlophis, Oxyrhopus, Boiruna and Phimophis. On the other hand, Pyron et al. [20] did not included Elapomorphus in their analysis. They showed Psomophis as a sister group of the clade formed by Uromacer, Erythrolamprus, and Lystrophis.

According to the phylogenetic signal metric K, there is significant phylogenetic signal in cranial shape. For the dorsal view, the mean K values from 200 trees is 0.95 (Pmean = 0.005) and similarly for lateral view, K mean is 0.98 (Pmean = 0.003) (see the histogram of K and P values in S1 Fig). Despite the strong phylogenetic signal results, projecting the consensus phylogenetic tree into the shape space for each view of the skull demonstrates a striking pattern (Figs 3 and 4). Some sister-taxa are widely separated in shape space, indicating they have very different morphologies. Species from Elapomorphini tribe are separated from all other species along PC1 and are more similar to each other in morphology, and share the same diet.

Xenodontine skull shape is highly correlated with diet. PPLS of dorsal skull shape on diet across the 200 trees resulted in a mean value of the degree of covariance accessed by permutation tests of 0.87 (Pmean = 0.003), and for lateral skull shape a mean value of the degree of covariance of 0.84 (Pmean = 0.01) (see the histogram of the values in S2 Fig).

In order to determine whether the species from Elapomorphini tribe (Apostolepis assimilis, Elapomorphus quinquelineatus, and Phalotris mertensi) were driving the association between diet preferences and skull morphology, we repeated the PPLS analyses of both skull views after excluding these taxa from the consensus tree. The dorsal skull shape resulted in a correlation of 0.80 (P = 0.007), and the lateral skull shape a correlation of 0.81 (P = 0.02), thus these results indicate that this clade has little influence on the overall results.

Discussion

Our results show that skull shape morphology in xenodontine snakes is strongly associated with diet [52]. Many studies have suggested association between snake skull morphology and diet in snakes. Maxillary teeth have been associated to different types of diet pronounced posterior ridges located on the posterior maxillary teeth are associated with slug predators [53] long, sharper teeth and elongated mandible bones are associated with piscivory [53], and enlarged anterior teeth and arched maxilla are thought to help Lycodon aulicus capucinus to ingest hard-bodied skinks [54]. Natricine piscivorous snakes were broadly studied, and the results show that fish-eating snakes have relatively longer skulls than frog-eating snakes, which tended to have broader skull components [55, 56, 57]. Studies analyzing the diversity of diet preferences usually focus on a broad scale analysis, for example, Vincent et al. [16] analyzed 12 monophyletic clades across macrostomatan snakes and found that head width is significantly related to the mean of the consumed prey mass, suggesting that skull in snakes is adapted to prey size (see also [14; 11]). Only few studies analyzed the evolution of skull morphology and diet at finer scales, such as the diversity of diet within a family or a genus [33].

Our study is the first to both describe snake skull morphology using geometric morphometric tools, and to correlate it to a broad spectrum of diet preferences using phylogenetic comparative methods in a highly diverse snake subfamily. However, our results alone are not able to state that diet drives the diversification of the skull morphology [52]. Our results highlighted a high correlation between diet and skull shape in xenodontine snakes; however we fail to define causality. To do so, it would be necessary to future studies of skull functionality and behavior directly associated to feeding.

The phylogenetic partial least squares analyses showed high correspondence of skull shape and diet that can be interpreted as an indicative of evolutionary convergence. Morphological convergence is an excellent opportunity to study adaptation at the macroevolutionary scale, since it offers multiple independent tests of the morphological response to similar functional demands [58]. However, as Revell et al. [2] pointed, morphological similarities among taxa could be a result of several causal processes, including evolutionary convergence, and evolutionary parallelism. Evolutionary convergence occurs when lineages with different ancestral morphologies evolve to the same phenotype through different directions, while in evolutionary parallelism lineages with the same ancestral morphology evolve in the same direction toward the same phenotype. In other words, parallelism is a special case of convergent evolution where independent evolutionary lineages evolve the same trait using the same genes and/or developmental pathways, while convergence occurs when lineages evolve similar traits using different developmental pathways. Species that evolve under parallel evolution are likely to share similar patterns of genetic covariation, and as a result, convergent evolution of a trait may occur not because selection favored that trait, but because it favored the same correlated trait in each species. As a consequence, natural selection could be favoring a correlated trait, and not the focal one, and any conclusion should carefully draw from evolutionary convergence and natural selection [59]. A more comprehensive sampling of taxa, combined with a robust and well resolved phylogeny of xenodontine snakes, is required to distinguish between these alternative hypotheses.

Vincent et al. [14], based on the analysis of skull size variation across several snakes species, showed that the increase in head width in snakes is followed by changes in jaw length and lower jaw out-lever length, suggesting that the snakes skulls were highly integrated. Skull parts are integrated because each part tightly associated to the other. Morphological integration can be a result of developmental, functional, genetic, and/or evolutionary interaction and constraints [60]. On the other hand, the skull is also modular, since morphological integration is not uniform throughout the entire skull, and it can be divided in modules that are strongly integrated internally, but are relatively independent of other modules [61]. Evolutionary integration and modularity can be investigated by examining how evolutionary changes in multiple parts are coordinated across a set of related species, and using comparative methods to take into account the phylogenetic structure of the variation [60]. Future studies are needed to test the modularity and integration among snake skull parts to better picture the forces influencing skull shape evolution.

Supporting Information

S1 Fig. Phylogenetic Signal of skull shape across 200 phylogenetic trees.

(A) Lateral cranial shape Phylogenetic signal values and (B) P values. (C) Dorsal cranial shape Phylogenetic signal values and (D) P values.

https://doi.org/10.1371/journal.pone.0148375.s001

(TIF)

S2 Fig. Phylogenetic Partial Least Squares of skull shape on diet across 200 phylogenetic trees.

(A) Lateral cranial shape correlations values and (B) P values. (C) Dorsal cranial shape correlations values and (D) P values.

https://doi.org/10.1371/journal.pone.0148375.s002

(TIF)

S1 File. List of analyzed material (Table A). Cranial landmarks definitions recorded from South American Xenodontinae snakes (Table B). Genbank access numbers (Table C). South American Xenodontinae phylogeny.

Bayesian inference phylogeny (mtDNA 12S and 16S, and c-mos) of South American Xenodontinae. Numbers on the branches represent the posterior probabilities (Figure A).

https://doi.org/10.1371/journal.pone.0148375.s003

(DOCX)

Acknowledgments

We thank R. Fernandes (Museu Nacional, Rio de Janeiro, Brasil), F. Franco (Instituto Butantan, São Paulo, Brasil), M. Carvalho (Universidade Federal do Mato Grosso, Cuiabá, Brasil) and P. Manzani (Museu de Zoologia, Unicamp, Brasil) for granting access to the xenodontine specimens under their care.

Author Contributions

Conceived and designed the experiments: JK. Performed the experiments: JK. Analyzed the data: JK ES. Contributed reagents/materials/analysis tools: JK EZFS. Wrote the paper: JK ES EZFS.

References

  1. 1. Schluter D. The Ecology of Adaptive Radiation. Oxford University of Press, Oxford; 2000.
  2. 2. Revell LJ, Johnson MA, Schulte JA, Kolbe JJ, Losos JB. A Phylogenetic Test For Adaptive Convergence In Rock-Dwelling Lizards. Evolution, 2007; 61(12): 2898–2912. pmid:17894806
  3. 3. Serb JM, Alejandrino A, Otárola-Castillo E, Adams DC. Morphological convergence of shell shape in distantly related scallop species (Mollusca: Pectinidae). Zoological Journal of the Linnean Society, 2011; 163(2): 571–584.
  4. 4. Hampton PM. Comparison of cranial form and function in association with diet in natricine snakes. Journal of Morphology 2011; 272(12): 1435–1443. pmid:21780158
  5. 5. Herrel A, Vincent SE, Alfaro ME, van Wassenbergh S, Vanhooydonck B, Irschick DJ. Morphological convergence as a consequence of extreme functional demands: examples from the feeding system of natricine snakes. Journal of Evolutionary Biology 2008; 21: 1438–1448. pmid:18547353
  6. 6. Vincent SE, Brandley MC, Herrel A, Alfaro ME. Convergence in trophic morphology and feeding performance among piscivorous natricine snakes. Journal of Evolutionary Biology. 2009; 22(6): 1203–1211 pmid:19389153
  7. 7. Monteiro LR, Nogueira MR. Ecological variation in South American geophagine cichlids arose during an early burst of adaptive morphological and functional evolution. BMC Evolutionary Biology. 2011; 280(1763): 20130849–20130849.
  8. 8. Santana SE, Grosse IR, Dumont ER. Dietary hardness, loading behavior, and the evolution of skull form in bats. Evolution. 2012; 66(8): 2587–2598. pmid:22834755
  9. 9. Stayton T. Morphological Evolution of the Lizard Skull: A Geometric Morphometrics Survey. Journal of Morphology. 2005; 263:47–59 pmid:15536647
  10. 10. Savitzky AH. Coadapted character complexes among snakes: fossoriality, piscivory, and durophagy. American Zoologist. 1983; 23(2): 397–409.
  11. 11. Cundall D, Greene HW. Feeding in snakes. Feeding: form, function, and evolution in tetrapod vertebrates, 2000; 293–333.
  12. 12. Greene HW. Dietary correlates of the origin and radiation of snakes. American Zoologist 1983; 23(2): 431–441.
  13. 13. Colston TJ, Costa GC, Vitt LJ. Snake diets and the deep history hypothesis. Biological Journal of the Linnean Society. 2010; 101(2): 476–486.
  14. 14. Gans C. The feeding mechanism of snakes and its possible evolution. American Zoologist. 1961; 1(2): 217–227.
  15. 15. Vincent SE, Dang PD, Herrel A, Kley NJ. Morphological integration and adaptation in the snake feeding system: a comparative phylogenetic study. Journal of Evolutionary Biology 2006; 19(5): 1545–1554. pmid:16910984
  16. 16. Hsiang AY, Field DJ, Webster TH, Behlke AD, Davis MB, Racicot RA, et al. The origin of snakes: revealing the ecology, behavior, and evolutionary history of early snakes using genomics, phenomics, and the fossil record. BMC evolutionary biology, 2015; 15(1), 87.
  17. 17. Jackson K, Fritts TH. Dentitional specialisations for durophagy in the Common Wolf snake, Lycodon aulicus capucinus. Amphibia-Reptilia, 2004; 25(3): 247–254.
  18. 18. Knox A, Jackson K. Ecological and phylogenetic influences on maxillary dentition in snakes. Phyllomedusa: Journal of Herpetology, 2010; 9: 121–131.
  19. 19. Vidal N, Dewynter M, Gower DJ. Dissecting the major American snake radiation: a molecular phylogeny of the Dipsadidae Bonaparte (Serpentes, Caenophidia). Comptes Rendus Biologies, 2010; 333: 48–55. pmid:20176336
  20. 20. Pyron RA, Burbrink FT, Colli GR, de Oca ANM, Vitt LJ, Kuczynski CA, et al. The phylogeny of advanced snakes (Colubroidea), with discovery of a new subfamily and comparison of support methods for likelihood trees. Molecular Phylogenetics and Evolution, 2011; 58: 329–342. pmid:21074626
  21. 21. Cadle JE. Molecular Systematics of Neotropical Xenodontine snakes: I South American xenodontines. Herpetologica 1984; 40: 8–20
  22. 22. Marques OAV, Pereira DN, Barbo FE, Germano VJ, Sawaya RJ. Os Répteis do Município de São Paulo: diversidade e ecologia da fauna pretérita e atual. Biota Neotropica 2009; 9 (2): 139–150.
  23. 23. Bookstein FL. Morphometric Tools for Landmark Data: Geometry and Biology. Cambridge University Press, Cambridge; 1991.
  24. 24. Mitteroecker P, Gunz P. Advances in geometric morphometrics. Evolutionary Biology. 2009; 36: 235–247.
  25. 25. Rohlf FJ. Tpsdig, 2008; Version 2.12; tpsrelw, Version 1.46. Stony Brook, NY: State University of New York at Stony Brook, Available: http://life.bio.sunysb.edu/morph/
  26. 26. Larson PM. Chondrocranial morphology and ontogenetic allometry in larval Bufo americanus (Anura, Bufonidae). Zoomorphology 2004; 123: 95–106
  27. 27. Rohlf FJ, Slice D. Extensions of the Procrustes method for the optimal superimposition of landmarks. Systematic Zoology 1990; 39: 40–59.
  28. 28. Rohlf FJ. Relative-warp analysis and an example of its application to mosquito wings. In: Marcus L.F., Bello E. Garcia-Valdecasas A. (Eds.), Contributions to Morphometrics. Museo Nacional de Ciencias Naturales, Madrid, 1993; 131–159.
  29. 29. Adams DC, Collyer ML, Otarola-Castillo E, Sherratt E. Geomorph: Software for geometric morphometric analyses. R package version 2.1 2014; Available: http://cran.r-project.org/web/packages/geomorph/index.html.
  30. 30. Adams DC, Otárola-Castillo E. Geomorph: an r package for the collection and analysis of geometric morphometric shape data. Methods in Ecology and Evolution 2013; 4: 393–399.
  31. 31. R Development Core Team. R: a language and environment for statistical computing. Vienna, Austria. 2014.
  32. 32. França FGR, Araújo AFB. The conservation status of snakes in central Brazil. South American Journal of Herpetology, 2006; 1(1): 25–36.
  33. 33. Alencar LR, Gaiarsa MP, Martins M. The evolution of diet and microhabitat use in Pseudoboine snakes. South American Journal of Herpetology, 2013; 8(1): 60–66.
  34. 34. Travaglia-Cardoso SR, Puorto G, Marques OA. Elapomorphus quinquelineatus (five-lined burrowing snake): Feeding on squamate eggs. The Herpetological Bulletin, 2014; 129: 28–29.
  35. 35. Marques OAV, Puorto G. Dieta e comportamento alimentar de Erythrolamprus aesculapii, uma serpente ofiófaga. Rev. Bras. Biol. 1994; 54: 253–259.
  36. 36. Oliveira L, Borges M, Marques OAV. Gomesophis brasiliensis (NCN). Reproduction and diet. Herp. Rev. 2013; 34: 251–252.
  37. 37. Scartozzoni RR. Estratégias reprodutivas e ecologia alimentar de serpentes aquáticas da tribo Hydropsini (Dipsadidae, Xenodontinae). Programa de Pós-Graduação Interunidades em Biotecnologia USP/Instituto Butantan/IPT. Volume Doctorate. São Paulo. Universidade de São Paulo, Brasil 2009.
  38. 38. Strussmann C, Sazima I. The snake assemblage of the Pantanal at Poconé, western Brazil: faunal composition and ecological summary. Studies on Neotropical Fauna and Environment, 1993; 28(3): 157–168.
  39. 39. Oliveira RB, Di-Bernardo M, Funk Pontes GM. Dieta e comportamento alimentar da Cobra-nariguda Lystrophis dorbignyi (Dumeril, Bibrón & Duméril, 1854), no litoral norte do Rio Grande do Sul, Brasil. Cuad. Herpetol. 2001; 14 (2): 117–122.
  40. 40. Bizerra A, Marques OA, Sazima I. Reproduction and feeding of the colubrid snake Tomodon dorsatus from south-eastern Brazil. Amphibia-Reptilia, 2005; 26(1), 33–38.
  41. 41. De Oliveira FS. Ecologia alimentar e reprodutiva de duas espécies de Tropidodryas (Serpentes, Colubridae) da Mata Atlântica. Volume Doctorate. São Paulo. Universidade de São Paulo, Brasil 2008.
  42. 42. Henderson RW, Schwartz A, Noeske-Hallin TA. Food habits of three colubrid tree snakes (genus Uromacer) on Hispaniola. Herpetologica, 1987; 241–248.
  43. 43. Martins M, Oliveira EM. Natural history of snakes in forests of the Manaus region, Central Amazonia, Brazil. Herpetol. Nat. Hist. 1998; 6:78–150.
  44. 44. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics 2007; 23: 2947–2948. pmid:17846036
  45. 45. Posada D. jModelTest: phylogenetic model averaging. Molecular Biology and Evolution 2008; 25: 1253–1256. pmid:18397919
  46. 46. Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003; 19: 1572–1574. pmid:12912839
  47. 47. Rambaut A, Suchard MA, Xie D, Drummond AJ. Tracer v1.6, 2014Available: http://beast.bio.ed.ac.uk/Tracer
  48. 48. Adams DC. A generalized K statistic for estimating phylogenetic signal from shape and other high-dimensional multivariate data. Systematic Biology. 2014; 63: 685–697 pmid:24789073
  49. 49. Klingenberg CP. MorphoJ: an integrated software package for geometric morphometrics. Molecular Ecology Resources. 2011; 11: 353–357. pmid:21429143
  50. 50. Adams DC, Felice RN. Assessing trait covariation and morphological integration on phylogenies using evolutionary covariance matrices. PLOS ONE, 2014; 9(4): e94335. pmid:24728003
  51. 51. Rohlf FJ, Corti M. Use of two-block partial least squares to study covariation in shape. Syst. Biol. 2000; 49: 740–753. pmid:12116437
  52. 52. Schluter D, Grant PR. Determinants of morphological patterns in communities of Darwin's finches. American Naturalist, 1984; 175–196.
  53. 53. Britt EJ, Clark AJ, Bennett AF. Dental morphologies in gartersnakes (Thamnophis) and their connection to dietary preferences. Journal of Herpetology, 2009; 43(2): 252–259.
  54. 54. Jackson K, Fritts TH. Dentitional specialisations for durophagy in the Common Wolf snake, Lycodon aulicus capucinus Amphibia-Reptilia 2004; 25: 247–254
  55. 55. Vincent SE, Brandley MC, Herrel A, Alfaro ME. Convergence in trophic morphology and feeding performance among piscivorous natricine snakes. Journal of Evolutionary Biology. 2009; 22: 1203–1211. pmid:19389153
  56. 56. Hampton PM. Comparison of cranial form and function in association with diet in natricine snakes. Journal of Morphology 2011; 272:1435–1443. pmid:21780158
  57. 57. Brecko J, Vervust B, Herrel A, Van Damme R. Head Morphology and Diet in the Dice Snake (Natrix tessellata). Mertensiella 2011; 18: 20–29
  58. 58. Larson A, Losos JB. Phylogenetic systematics of adaptation. In Rose M.R. Lauder G.V., Eds., Adaptation. Pp. 187–220. Academic Press: San Diego 1996
  59. 59. Losos JB. Convergence, adaptation, and constraint. Evolution 2011; 65: 1827–1840. pmid:21729041
  60. 60. Klingenberg C. Cranial integration and modularity: insights into evolution and development from morphometric data. Hystrix, the Italian Journal of Mammalogy, 2013; 24: 43–58.
  61. 61. Klingenberg C. Morphological integration and developmental modularity. Annual review of ecology, evolution, and systematics 2008: 115–132.