Complex Patterns of Diversification in the Gray Zone of Speciation: Model-Based Approaches Applied to Patagonian Liolaemid Lizards (Squamata: Liolaemus kingii clade)

In this study we detangled the evolutionary history of the Patagonian lizard clade Liolaemus kingii, coupling dense geographic sampling and novel computational analytical approaches. We analyzed nuclear and mitochondrial data (restriction site-associated DNA sequencing and cytochrome b) to hypothesize and evaluate species limits, phylogenetic relationships, and demographic histories. We complemented these analyses with posterior predictive simulations to assess the fit of the genomic data to the multispecies coalescent model. We also employed a novel approach to time-calibrate a phylogenetic network. Our results show several instances of mito-nuclear discordance and consistent support for a reticulated history, supporting the view that the complex evolutionary history of the kingii clade is characterized by extensive gene flow and rapid diversification events. We discuss our findings in the contexts of the "gray zone" of speciation, phylogeographic patterns in the Patagonian region, and taxonomic outcomes. [Model adequacy; multispecies coalescent; multispecies network coalescent; phylogenomics; species delimitation.].

Native Lizards Living in Brazilian Cities: Effects of Developmental Environments on Thermal Sensitivity and Morpho-Functional Associations of Locomotion

Environmental conditions often affect developmental processes and consequently influence the range of phenotypic variation expressed at population level. Expansion of urban sites poses new challenges for native species, as urbanization usually affects the intensity of solar exposure and shade availability, determining the thermal regimes organisms are exposed to. In this study, we evaluate the effects of different developmental conditions in a Tropidurus lizard commonly found in Brazilian urban sites. After incubating embryos of Tropidurus catalanensis in two different thermal regimes (Developmental Environments [DE]: cold 24°C and warm 30°C), we measured morphological traits in the neonates and quantified locomotor performance in horizontal and vertical surfaces at three temperatures [Test Temperatures (TT) = 24°C, 30°C and 36°C]. Results indicate effects of developmental temperatures on morphological features, expressing functional implications that might be decisive for the viability of T. catalanensis in urbanized areas. Lizards ran similarly on horizontal and vertical surfaces, and isolated analyses did not identify significant effects of DE or TT on the sprint speeds measured. Absolute Vmax (i.e., the maximum sprint speed reached among all TTs) positively correlated with body size (SVL), and neonates from the warm DE (30°C) were larger than those from the cold DE (24°C). Morpho-functional associations of absolute Vmax also involved pelvic girdle width and forelimb, hindlimb, trunk, and tail lengths. Emerging discussions aim to understand how animals cope with abrupt environmental shifts, a likely common challenge in urbanized sites. Our findings add a new dimension to the topic, providing evidence that temperature, an environmental parameter often affected by urbanization, influences the thermal sensitivity of locomotion and the morphological profile of T. catalanensis neonates. Thermal sensitivity of specific developmental processes may influence the ability of these lizards to remain in habitats that change fast, as those suffering rapid urbanization due to city growth.

Environmental correlates of phenotypic evolution in ecologically diverse Liolaemus lizards

Evolutionary correlations between phenotypic and environmental traits characterize adaptive radiations. However, the lizard genus Liolaemus, one of the most ecologically diverse terrestrial vertebrate radiations on earth, has so far shown limited or mixed evidence of adaptive diversification in phenotype. Restricted use of comprehensive environmental data, incomplete taxonomic representation and not considering phylogenetic uncertainty may have led to contradictory evidence. We compiled a 26‐taxon dataset for the Liolaemus gracilis species group, representing much of the ecological diversity represented within Liolaemus and used environmental data to characterize how environments occupied by species' relate to phenotypic evolution. Our analyses, explicitly accounting for phylogenetic uncertainty, suggest diversification in phenotypic traits toward the present, with body shape evolution rapidly evolving in this group. Body shape evolution correlates with the occupation of different structural habitats indicated by vegetation axes suggesting species have adapted for maximal locomotory performance in these habitats. Our results also imply that the effects of phylogenetic uncertainty and model misspecification may be more extensive on univariate, relative to multivariate analyses of evolutionary correlations, which is an important consideration in analyzing data from rapidly radiating adaptive radiations.

A guide to incubate eggs of Tropidurus lizards under laboratory conditions

Studies in Evo-Devo benefit from the use of a variety of organisms, as comparative approaches provide a better understanding of Biodiversity and Evolution. Standardized protocols to incubate eggs and manipulate embryo development enable postulation of additional species as suitable biological systems for research in the field. In the past decades, vertebrate lineages such as Squamata (lizards, snakes, and amphisbaenians) emerged as crucial study systems for addressing topics as diverse as phenotypic evolution and climate change. However, protocols for maintaining gravid females and incubating eggs in the lab under experimental conditions are available to only a few squamate species. This resource article presents a simple incubation guide that standardizes conditions to maintain embryos of Tropidurus catalanensis (Squamata: Tropiduridae) under different experimental conditions, manipulating relevant environmental factors like temperature and humidity. We identified associated effects relating the egg incubation condition to developmental stage, incubation time, hatching success, and resulting morphotypes. Temperature and humidity play a key role in development and require attention when establishing the experimental design. Current literature comprises information for Tropidurus lizards that ponders how general in Squamata are the ecomorphs originally described for Anolis. Studies evaluating phenotypic effects of developmental environments suggest plasticity in some of the traits that characterize the ecomorphological associations described for this family. We expect that this incubation guide encourages future studies using Tropidurus lizards to address Evo-Devo questions.

Rapid radiation and rampant reticulation: Phylogenomics of South American Liolaemus lizards

Understanding the factors that cause heterogeneity among gene trees can increase the accuracy of species trees. Discordant signals across the genome are commonly produced by incomplete lineage sorting (ILS) and introgression, which in turn can result in reticulate evolution. Species tree inference using the multispecies coalescent is designed to deal with ILS and is robust to low levels of introgression, but extensive introgression violates the fundamental assumption that relationships are strictly bifurcating. In this study, we explore the phylogenomics of the iconic Liolaemus subgenus of South American lizards, a group of over 100 species mostly distributed in and around the Andes mountains. Using mitochondrial DNA (mtDNA) and genome-wide restriction-site associated DNA sequencing (RADseq; nDNA hereafter), we inferred a time-calibrated mtDNA gene tree, nDNA species trees, and phylogenetic networks. We found high levels of discordance between mtDNA and nDNA, which we attribute in part to extensive ILS resulting from rapid diversification. These data also reveal extensive and deep introgression, which combined with rapid diversification, explain the high level of phylogenetic discordance. We discuss these findings in the context of Andean orogeny and glacial cycles that fragmented, expanded, and contracted species distributions. Finally, we use the new phylogeny to resolve long-standing taxonomic issues in one of the most studied lizard groups in the New World.

Phylogenomic data resolve higher-level relationships within South American Liolaemus lizards

Phylogenomic approaches now generate hundreds of loci representative of the whole genome that can be used for phylogenetic analyses. The South American lizard genus Liolaemus is the most species-rich vertebrate radiation from temperate zones (more than 265 described species), yet most higher-level phylogenetic relationships within Liolaemus remain poorly resolved. In this study, we used 584 nuclear loci collected using targeted sequenced capture to estimate the phylogenetic relationships among 26 species representing the two subgenera within Liolaemus (Eulaemus + Liolaemus), and all major groups within Eulaemus. Previous molecular and morphological-based phylogenetic analyses of Eulaemus based on a limited number of characters resolved few higher-level relationships, although one point of agreement is that the early divergence within Eulaemus corresponds to the lineomaculatus section, followed by the diversification of eight main clades that are strongly supported and recognized. Liolaemus probably experienced relatively rapid divergences during parts of its evolutionary history, and a phylogenomic approach was used to resolve the relationships among the major groups. The new analyses presented here support the division of Liolaemus into two subgenera, and resolve relationships among many of the major clades of Eulaemus with strong support. A Bayesian divergence dating analysis using 44 protein-coding genes provides an estimation of the split of the two Liolaemus subgenera of approximately 19,7 ma (95% HPD = 16,94-23,04), while diversification within Eulaemus started at 15,05 ma (95% HPD = 12,94 - 17,59) among the L. lineomaculatus and the L. montanus series by Mid Miocene. A novel phylogenetic network analyses for SNP data identified two hybridizing edges among different groups of Eulaemus at different points in time. Having a solid phylogenetic hypothesis of the main Eulaemus clades opens new opportunities to test a variety of macroevolutionary questions for this unique radiation.

Anuran forelimb muscle tendinous structures and their relationship with locomotor modes and habitat use

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 together 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.

You Can't Run, but You Can Hide: The Skeleton of the Sand-Swimmer Lizard Calyptommatus leiolepis (Squamata: Gymnophthalmidae)

Squamates exhibit a vast diversity of body plans, which directly determines habitat use and preference. Here the skeleton of the sand-swimmer burrower gymnophthalmid, Calyptommatus leiolepis, is analyzed to investigate how its peculiar fossorial locomotion affects its osteology. Calyptommatus leiolepis is a limb-reduced, short-intermediate tailed lizard. Although there are other studies on its general anatomy, we performed a detailed description of its skeleton. Using high-resolution computer tomography, each bone element within the skeleton was digitally segmented and a detailed description rendered. Anatomical features related to burrowing include the head having a shovel-like snout with a well-developed horizontal soft tissue ridge, nasal cartilages that exclude sand from the nostrils, reduced eyes covered by a brille, lack of forelimbs, extreme reduction of hind limbs, and imbricated scales among others. The genus Calyptommatus has unique features such as a triradiate jugal (with digit-like posterior projections), a reduced pectoral girdle and forelimbs, parasternal processes that interconnect the ribs, and a single digit in the hind limbs. When comparing this species with other gymnophthalmid lizards including, fossorial species, it is clear that Calyptommatus exhibits the highest number of structural modifications within the family. Despite its specialized morphology, it still retains characters that link this genus to other members of Gymnophthalmidae when included in a phylogeny based solely on phenotypic data. Anat Rec, 303:1305-1326, 2020. © 2019 American Association for Anatomy.

Macroevolution of sexual size dimorphism and reproduction-related phenotypic traits in lizards of the Chaco Domain

BACKGROUND

Comparing sexual size dimorphism (SSD) in the light of the phylogenetic hypothesis may help to understand the phenotypic evolution associated with sexual selection (size of whole body and of reproduction-related body parts). Within a macroevolutionary framework, we evaluated the association between the evolution of SSD and the evolution of reproduction-related phenotypic traits, and whether this association has favored female fecundity, considering also variations according to reproductive modes. We focused on the lizard species that inhabit the Chaco Domain since this is a natural unit with a high diversity of species.

RESULTS

The residual SSD was related positively with the residuals of the reproduction-related phenotypic traits that estimate intrasexual selection and with the residuals of inter-limb length and, according to fecundity selection, those residuals were related positively with the residuals of clutch size in oviparous species. Lizards of the Chaco Domain present a high diversity of SSD patterns, probably related to the evolution of reproductive strategies.

CONCLUSIONS

Our findings highlight that the sexual selection may have acted on the whole-body size as well as on the size of body parts related to reproduction. Male and female phenotypes evolutionarily respond to variations in SSD, and an understanding of these patterns is essential for elucidating the processes shaping sexual phenotype diversity from a macroevolutionary perspective.

Environmental determinism, and not interspecific competition, drives morphological variability in Australasian warblers (Acanthizidae)

Interspecific competition is thought to play a key role in determining the coexistence of closely related species within adaptive radiations. Competition for ecological resources can lead to different outcomes from character displacement to, ultimately, competitive exclusion. Accordingly, divergent natural selection should disfavor those species that are the most similar to their competitor in resource use, thereby increasing morphological disparity. Here, we examined ecomorphological variability within an Australo‐Papuan bird radiation, the Acanthizidae, which include both allopatric and sympatric complexes. In addition, we investigated whether morphological similarities between species are related to environmental factors at fine scale (foraging niche) and/or large scale (climate). Contrary to that predicted by the competition hypothesis, we did not find a significant correlation between the morphological similarities found between species and their degree of range overlap. Comparative modeling based on both a priori and data‐driven identification of selective regimes suggested that foraging niche is a poor predictor of morphological variability in acanthizids. By contrast, our results indicate that climatic conditions were an important factor in the formation of morphological variation. We found a significant negative correlation between species scores for PC1 (positively associated to tarsus length and tail length) and both temperature and precipitation, whereas PC2 (positively associated to bill length and wing length) correlated positively with precipitation. In addition, we found that species inhabiting the same region are closer to each other in morphospace than to species outside that region regardless of genus to which they belong or its foraging strategy. Our results indicate that the conservative body form of acanthizids is one that can work under a wide variety of environments (an all‐purpose morphology), and the observed interspecific similarity is probably driven by the common response to environment.

Phylogenetic patterns and correlation of key structures for jumping: bone crests and cross-sectional areas of muscles in Leptodactylus (Anura, Leptodactylidae)

Anurans are characterized by their saltatory mode of locomotion, which is associated with a specific morphology. The coordinated action of the muscles and bones of the pelvic girdle is key to the transmission of the force of the hindlimbs to the axial skeleton during jumping. Two features are critical for optimal locomotory performance: the cross-sectional area of muscle and the bone crest attachment sites. The first character is a proxy of the force exerted by the muscle, whereas the crests are muscle attachments sites related to muscle force. The provisory relationship between these features has previously been identified and bone crest size can be used to infer the magnitude and, therefore, muscle force in fossils records. In this work, we explore the correlation between the cross-sectional area of essential muscles to the jumping mechanism (longissimus dorsi, extensor iliotibialis B, tenuissimus, puboischiofemoralis internus B, coccygeo-sacralis and coccygeo-iliacus) and the bone crests where these muscles are inserted (dorsal tubercle, dorsal crest and urostylar crest) in species of the genus Leptodactylus. This genus, along with other leptodactylids, exhibits a diversity of locomotor modes, including jumping, hopping, swimming and burrowing. We therefore analyzed the morphometric variation in the two features, cross-sectional area and bone crest area, expecting a correlation with different locomotor types. Our results showed: (i) a correlation between the urostylar crest and the cross-sectional area of the related muscles; (ii) that the bone crest surface area of urostyle and ilium and the cross-sectional area of the corresponding muscles can be utilized to infer locomotor faculties in leptodactylid frogs; and (iii) that the evolution of both characters demonstrates a general tendency from lower values in leptodactylid ancestors to higher values in the Leptodactylus genus. The results attest to the importance of the comparison of current ecological and phylogenetic analogues as they allow us to infer functionality and behavior in fossil and extant groups based on skeletal evidence. Phylogenetic patterns in character evolution and their correlation with locomotory types could imply that functional restrictions are also inherited in leptodactylid.

Tendinous framework of anurans reveals an all-purpose morphology

Tendons are directly associated with movement, amplifying power and reducing muscular work. Taking into account habitat and locomotor challenges faced by anurans, we identify the more conspicuous superficial tendons of a neotropical anuran group and investigate their relation to the former factors. We show that tendons can be visualized as an anatomical framework connected through muscles and/or fascia, and describe the most superficial tendinous layer of the postcranium of Leptodactylus latinasus. To analyze the relation between tendon morphology and ecological characters, we test the relative length ratio of 10 tendon-muscle (t-m) elements in 45 leptodactylid species while taking phylogeny into account. We identify the evolutionary model that best explains our variables. Additionally, we optimize t-m ratio values, and the shape of the longissimus dorsi insertion onto a selected phylogeny of the species. Our data show the existence of an all-purpose morphology that seems to have evolved independently of ecology and functional requirements. This is indicated by no significant relation between morphometric data of the analyzed tendons and habitat use or locomotion, a strong phylogenetic component to most of the analyzed variables, and a generalized pattern of intermediate values for ancestral states. Ornstein-Uhlenbeck is the model that best explains most t-m variables, indicating that stabilizing selection or selective optima might be driving shifts in tendon length within Leptodactylidae. Herein, we show the substantial influence that phylogeny has on tendon morphology, demonstrating that a generalized and stable morphological configuration of tendons is adequate to enable versatile locomotor modes and habitat use. This is an attempt to present the tendinous system as a framework to body support in vertebrates, and can be considered a starting point for further ecomorphological research of this anatomical system in anurans.