Including Fossils in Phylogenetic Climate Reconstructions: A Deep Time Perspective on the Climatic Niche Evolution and Diversification of Spiny Lizards (Sceloporus)

Exploring the causes underlying the latitudinal variation in range sizes: Evidence for Rapoport's rule in spiny lizards (genus Sceloporus)

Species' range size is a fundamental unit of analysis in biodiversity research, given its association with extinction risk and species richness. One of its most notable patterns is its positive relationship with latitude, which has been considered an ecogeographical rule called Rapoport's rule. Despite this rule being confirmed for various taxonomic groups, its validity has been widely discussed and several taxa still lack a formal assessment. Different hypotheses have been proposed to explain their potential mechanisms, with those related to temperature and elevational being the most supported thus far. In this study, we employed two level of analyses (cross-species and assemblage) to investigate the validity of Rapoport's rule in spiny lizards (genus Sceloporus). Additionally, we evaluated four environmental-related hypotheses (minimum temperature, temperature variability, temperature stability since the last glacial maximum, and elevation) posed to explain such pattern, contrasting our results to those patterns expected under a null model of range position. Our results provided support for Rapoport's rule at both levels of analyses, contrasting with null expectations. Consistently, minimum temperature and elevation were the most relevant variables explaining the spatial variation in range size. At the cross-species level, our null simulations revealed that both variables deviated significantly from random expectations. Conversely, at the assemblage level, none of the variables were statistically different from the expected relationships. We discussed the implication of our findings in relation to the ecology and evolution of spiny lizards.

Identifying Climatic Drivers of Hybridization with a New Ancestral Niche Reconstruction Method

Applications of molecular phylogenetic approaches have uncovered evidence of hybridization across numerous clades of life, yet the environmental factors responsible for driving opportunities for hybridization remain obscure. Verbal models implicating geographic range shifts that brought species together during the Pleistocene have often been invoked, but quantitative tests using paleoclimatic data are needed to validate these models. Here, we produce a phylogeny for Heuchereae, a clade of 15 genera and 83 species in Saxifragaceae, with complete sampling of recognized species, using 277 nuclear loci and nearly complete chloroplast genomes. We then employ an improved framework with a coalescent simulation approach to test and confirm previous hybridization hypotheses and identify one new intergeneric hybridization event. Focusing on the North American distribution of Heuchereae, we introduce and implement a newly developed approach to reconstruct potential past distributions for ancestral lineages across all species in the clade and across a paleoclimatic record extending from the late Pliocene. Time calibration based on both nuclear and chloroplast trees recovers a mid- to late-Pleistocene date for most inferred hybridization events, a timeframe concomitant with repeated geographic range restriction into overlapping refugia. Our results indicate an important role for past episodes of climate change, and the contrasting responses of species with differing ecological strategies, in generating novel patterns of range contact among plant communities and therefore new opportunities for hybridization. The new ancestral niche method flexibly models the shape of niche while incorporating diverse sources of uncertainty and will be an important addition to the current comparative methods toolkit. [Ancestral niche reconstruction; hybridization; paleoclimate; pleistocene.].

Integrative Phylogenetics: Tools for Palaeontologists to Explore the Tree of Life

The modern era of analytical and quantitative palaeobiology has only just begun, integrating methods such as morphological and molecular phylogenetics and divergence time estimation, as well as phenotypic and molecular rates of evolution. Calibrating the tree of life to geological time is at the nexus of many disparate disciplines, from palaeontology to molecular systematics and from geochronology to comparative genomics. Creating an evolutionary time scale of the major events that shaped biodiversity is key to all of these fields and draws from each of them. Different methodological approaches and data employed in various disciplines have traditionally made collaborative research efforts difficult among these disciplines. However, the development of new methods is bridging the historical gap between fields, providing a holistic perspective on organismal evolutionary history, integrating all of the available evidence from living and fossil species. Because phylogenies with only extant taxa do not contain enough information to either calibrate the tree of life or fully infer macroevolutionary dynamics, phylogenies should preferably include both extant and extinct taxa, which can only be achieved through the inclusion of phenotypic data. This integrative phylogenetic approach provides ample and novel opportunities for evolutionary biologists to benefit from palaeontological data to help establish an evolutionary time scale and to test core macroevolutionary hypotheses about the drivers of biological diversification across various dimensions of organisms.

The recombination landscapes of spiny lizards (genus Sceloporus)

Despite playing a critical role in evolutionary processes and outcomes, relatively little is known about rates of recombination in the vast majority of species, including squamate reptiles-the second largest order of extant vertebrates, many species of which serve as important model organisms in evolutionary and ecological studies. This paucity of data has resulted in limited resolution on questions related to the causes and consequences of rate variation between species and populations, the determinants of within-genome rate variation, as well as the general tempo of recombination rate evolution on this branch of the tree of life. In order to address these questions, it is thus necessary to begin broadening our phylogenetic sampling. We here provide the first fine-scale recombination maps for two species of spiny lizards, Sceloporus jarrovii and Sceloporus megalepidurus, which diverged at least 12 Mya. As might be expected from similarities in karyotype, population-scaled recombination landscapes are largely conserved on the broad-scale. At the same time, considerable variation exists at the fine-scale, highlighting the importance of incorporating species-specific recombination maps in future population genomic studies.

A chromosome-level genome assembly for the eastern fence lizard (Sceloporus undulatus), a reptile model for physiological and evolutionary ecology

BACKGROUND

High-quality genomic resources facilitate investigations into behavioral ecology, morphological and physiological adaptations, and the evolution of genomic architecture. Lizards in the genus Sceloporus have a long history as important ecological, evolutionary, and physiological models, making them a valuable target for the development of genomic resources.

FINDINGS

We present a high-quality chromosome-level reference genome assembly, SceUnd1.0 (using 10X Genomics Chromium, HiC, and Pacific Biosciences data), and tissue/developmental stage transcriptomes for the eastern fence lizard, Sceloporus undulatus. We performed synteny analysis with other snake and lizard assemblies to identify broad patterns of chromosome evolution including the fusion of micro- and macrochromosomes. We also used this new assembly to provide improved reference-based genome assemblies for 34 additional Sceloporus species. Finally, we used RNAseq and whole-genome resequencing data to compare 3 assemblies, each representing an increased level of cost and effort: Supernova Assembly with data from 10X Genomics Chromium, HiRise Assembly that added data from HiC, and PBJelly Assembly that added data from Pacific Biosciences sequencing. We found that the Supernova Assembly contained the full genome and was a suitable reference for RNAseq and single-nucleotide polymorphism calling, but the chromosome-level scaffolds provided by the addition of HiC data allowed synteny and whole-genome association mapping analyses. The subsequent addition of PacBio data doubled the contig N50 but provided negligible gains in scaffold length.

CONCLUSIONS

These new genomic resources provide valuable tools for advanced molecular analysis of an organism that has become a model in physiology and evolutionary ecology.

Occurrence data uncover patterns of allopatric divergence and interspecies interactions in the evolutionary history of Sceloporus lizards

As shown from several long-term and time-intensive studies, closely related, sympatric species can impose strong selection on one another, leading to dramatic examples of phenotypic evolution. Here, we use occurrence data to identify clusters of sympatric Sceloporus lizard species and to test whether Sceloporus species tend to coexist with other species that differ in body size, as we would expect when there is competition between sympatric congeners. We found that Sceloporus species can be grouped into 16 unique bioregions. Bioregions that are located at higher latitudes tend to be larger and have fewer species, following Rapoport's rule and the latitudinal diversity gradient. Species richness was positively correlated with the number of biomes and elevation heterogeneity of each bioregion. Additionally, most bioregions show signs of phylogenetic underdispersion, meaning closely related species tend to occur in close geographic proximity. Finally, we found that although Sceloporus species that are similar in body size tend to cluster geographically, small-bodied Sceloporus species are more often in sympatry with larger-bodied Sceloporus species than expected by chance alone, whereas large-bodied species cluster with each other geographically and phylogenetically. These results suggest that community composition in extant Sceloporus species is the result of allopatric evolution, as closely related species move into different biomes, and interspecies interactions, with sympatry between species of different body sizes. Our phyloinformatic approach offers unique and detailed insights into how a clade composed of ecologically and morphologically disparate species are distributed over large geographic space and evolutionary time.

Comparison of environmental inference approaches for ecometric analyses: Using hypsodonty to estimate precipitation

Ecometrics is the study of community-level functional trait-environment relationships. We use ecometric analyses to estimate paleoenvironment and to investigate community-level functional changes through time.We evaluate four methods that have been used or have the potential to be used in ecometric analyses for estimating paleoenvironment to determine whether there have been systematic differences in paleoenvironmental estimation due to choice of the estimation method. Specifically, we evaluated linear regression, polynomial regression, nearest neighbor, and maximum-likelihood methods to explore the predictive ability of the relationship for a well-known ecometric dataset of mammalian herbivore hypsodonty metrics (molar tooth crown to root height ratio) and annual precipitation. Each method was applied to 43 Pleistocene fossil sites and compared to annual precipitation from global climate models. Sites were categorized as glacial or interglacial, and paleoprecipitation estimates were compared to the appropriate model.Estimation methods produce results that are highly correlated with log precipitation and estimates from the other methods (p < 0.001). Differences between estimated precipitation and observed precipitation are not significantly different across the four methods, but maximum likelihood produces the most accurate estimates of precipitation. When applied to paleontological sites, paleoprecipitation estimates align more closely with glacial global climate models than with interglacial models regardless of the age of the site.Each method has constraints that are important to consider when designing ecometric analyses to avoid misinterpretations when ecometric relationships are applied to the paleontological record. We show interglacial fauna estimates of paleoprecipitation more closely match glacial global climate models. This is likely because of the anthropogenic effects on community reassembly in the Holocene.

Volatile fatty acid and aldehyde abundances evolve with behavior and habitat temperature in Sceloporus lizards

Animal signals evolve by striking a balance between the need to convey information through particular habitats and the limitations of what types of signals can most easily be produced and perceived. Here, we present new results from field measures of undisturbed behavior and biochemical analyses of scent marks from 12 species of Sceloporus lizards to explore whether evolutionary changes in chemical composition are better predicted by measures of species behavior, particularly those associated with visual displays, chemoreception, and locomotion, or by measures of habitat climate (precipitation and temperature). We found that more active lizard species used fewer compounds in their volatile scent marks, perhaps conveying less specific information about individual and species identity. Scent marks from more active lizard species also had higher proportions of saturated fatty acids, and the evolution of these compounds has been tracking the phylogeny closely as we would expect for a metabolic byproduct. In contrast, the proportions of unsaturated fatty acids were better explained by evolutionary shifts in habitat temperature (and not precipitation), with species in warmer climates using almost no volatile unsaturated fatty acids. The proportion of aldehydes was explained by both behavior and environment, decreasing with behavioral activity and increasing with habitat temperature. Our results highlight the evolutionary flexibility of complex chemical signals, with different chemical compounds responding to different elements of the selective landscape over evolutionary time.

Integrating Agriculture and Ecosystems to Find Suitable Adaptations to Climate Change

Climate change is altering agricultural production and ecosystems around the world. Future projections indicate that additional change is expected in the coming decades, forcing individuals and communities to respond and adapt. Current research efforts typically examine climate change effects and possible adaptations but fail to integrate agriculture and ecosystems. This failure to jointly consider these systems and associated externalities may underestimate climate change impacts or cause adaptation implementation surprises, such as causing adaptation status of some groups or ecosystems to be worsened. This work describes and motivates reasons why ecosystems and agriculture adaptation require an integrated analytical approach. Synthesis of current literature and examples from Texas are used to explain concepts and current challenges. Texas is chosen because of its high agricultural output that is produced in close interrelationship with the surrounding semi-arid ecosystem. We conclude that future effect and adaptation analyses would be wise to jointly consider ecosystems and agriculture. Existing paradigms and useful methodology can be transplanted from the sustainable agriculture and ecosystem service literature to explore alternatives for climate adaptation and incentivization of private agriculturalists and consumers. Researchers are encouraged to adopt integrated modeling as a means to avoid implementation challenges and surprises when formulating and implementing adaptation.

Oscillayers: A dataset for the study of climatic oscillations over Plio-Pleistocene time-scales at high spatial-temporal resolution

MOTIVATION

In order to understand how species evolutionarily responded to Plio-Pleistocene climate oscillations (e.g. in terms of speciation, extinction, migration and adaptation), it is first important to have a good understanding of those past climate changes per se. This, however, is currently limited due to the lack of global-scale climatic datasets with high temporal resolution spanning the Plio-Pleistocene. To fill this gap, I here present Oscillayers, a global-scale and region-specific bioclim dataset, facilitating the study of climatic oscillations during the last 5.4 million years at high spatial (2.5 arc-minutes) and temporal (10 kyr time periods) resolution. This dataset builds upon interpolated anomalies (Δ layers) between bioclim layers of the present and the Last Glacial Maximum (LGM) that are scaled relative to the Plio-Pleistocene global mean temperature curve, derived from benthic stable oxygen isotope ratios, to generate bioclim variables for 539 time periods. Evaluation of the scaled, interpolated estimates of palaeo-climates generated for the Holocene, Last Interglacial and Pliocene showed good agreement with independent general circulation models (GCMs) for respective time periods in terms of pattern correlation and absolute differences. Oscillayers thus provides a new tool for studying spatial-temporal patterns of evolutionary and ecological processes at high temporal and spatial resolution.

MAIN TYPES OF VARIABLE CONTAINED

Nineteen bioclim variables for time periods throughout the Plio-Pleistocene. Input data and R script to recreate all 19 bioclim variables.

SPATIAL LOCATION AND GRAIN

Global at 2.5 arc-minutes (4.65 x 4.65 = 21.62 km2 at the equator).

TIME PERIOD AND GRAIN

The last 5.4 million years. The grain is 10 kyr (= 539 time periods).

LEVEL OF MEASUREMENT

Data are for terrestrial climates (excluding Antarctica) taking sea level changes into account.

SOFTWARE FORMAT

All data are available as ASCII grid files.

Thermal niche evolution across replicated Anolis lizard adaptive radiations

Elucidating how ecological and evolutionary mechanisms interact to produce and maintain biodiversity is a fundamental problem in evolutionary ecology. Here, we focus on how physiological evolution affects performance and species coexistence along the thermal niche axis in replicated radiations of Anolis lizards best known for resource partitioning based on morphological divergence. We find repeated divergence in thermal physiology within these radiations, and that this divergence significantly affects performance within natural thermal environments. Morphologically similar species that co-occur invariably differ in their thermal physiology, providing evidence that physiological divergence facilitates species coexistence within anole communities. Despite repeated divergence, phylogenetic comparative analyses indicate that physiological traits have evolved more slowly than key morphological traits related to the structural niche. Phylogenetic analyses also reveal that physiological divergence is correlated with divergence in broad-scale habitat climatic features commonly used to estimate thermal niche evolution, but that the latter incompletely predicts variation in the former. We provide comprehensive evidence for repeated adaptive evolution of physiological divergence within Anolis adaptive radiations, including the complementary roles of physiological and morphological divergence in promoting community-level diversity. We recommend greater integration of performance-based traits into analyses of climatic niche evolution, as they facilitate a more complete understanding of the phenotypic and ecological consequences of climatic divergence.

Repeated evolution of viviparity in phrynosomatid lizards constrained interspecific diversification in some life-history traits

In vertebrates, viviparity has evolved independently multiple times, apparently increasing morphological diversification and speciation rates as a consequence. We tested whether the evolution of viviparity has also increased diversification of life-history traits by estimating evolutionary rates of lizards from the North American family Phrynosomatidae. Using modern phylogenetic comparative methods, we compared these rates between oviparous and viviparous species, and found no support for this hypothesis. Instead, we found higher evolutionary rates for oviparous species in some life-history traits. Our results suggest that the evolution of viviparity may have constrained rather than facilitated evolution of life histories.