Local adaptation in mainland anole lizards: Integrating population history and genome-environment associations

Notes from the taxonomic disaster zone: Evolutionary drivers of intractable species boundaries in an Australian lizard clade (Scincidae: Ctenotus)

Genomic-scale datasets, sophisticated analytical techniques, and conceptual advances have disproportionately failed to resolve species boundaries in some groups relative to others. To understand the processes that underlie taxonomic intractability, we dissect the speciation history of an Australian lizard clade that arguably represents a "worst-case" scenario for species delimitation within vertebrates: the Ctenotus inornatus species group, a clade beset with decoupled genetic and phenotypic breaks, uncertain geographic ranges, and parallelism in purportedly diagnostic morphological characters. We sampled hundreds of localities to generate a genomic perspective on population divergence, structure, and admixture. Our results revealed rampant paraphyly of nominate taxa in the group, with lineages that are either morphologically cryptic or polytypic. Isolation-by-distance patterns reflect spatially continuous differentiation among certain pairs of putative species, yet genetic and geographic distances are decoupled in other pairs. Comparisons of mitochondrial and nuclear gene trees, tests of nuclear introgression, and historical demographic modelling identified gene flow between divergent candidate species. Levels of admixture are decoupled from phylogenetic relatedness; gene flow is often higher between sympatric species than between parapatric populations of the same species. Such idiosyncratic patterns of introgression contribute to species boundaries that are fuzzy while also varying in fuzziness. Our results suggest that "taxonomic disaster zones" like the C. inornatus species group result from spatial variation in the porosity of species boundaries and the resulting patterns of genetic and phenotypic variation. This study raises questions about the origin and persistence of hybridizing species and highlights the unique insights provided by taxa that have long eluded straightforward taxonomic categorization.

Local adaptation of Dromiciops marsupials (Microbiotheriidae) from southern South America: Implications for species management facing climate change

The two species of the microbiotheriid marsupial genus Dromiciops (Dromiciops bozinovici: "Panchos's monito del monte" and Dromiciops gliroides: "monito del monte") exhibit a marked latitudinal genetic differentiation. Nevertheless, it is unclear whether this differentiation results from neutral processes or can be explained, to some extent, by local adaptation to different environmental conditions. Here, we used an SNP panel gathered by Rad-seq and searched for footprints of local adaptation (putative loci under selection) by exploring genetic associations with environmental variables in the two species of Dromiciops in Chilean and Argentinean populations. We applied three methods for detecting outlier SNPs and two genotype-environment associations approaches to quantify associations between allelic frequencies and environmental variables. Both species display strong genetic structure. D. bozinovici exhibited three distinct genetic groups, marking the first report of such structuring in this species using SNPs. In contrast, D. gliroides displayed four genetic clusters, consistent with previous studies. Both species exhibited an association of their genetic structure with environmental variables. D. bozinovici exhibited significant associations of allelic frequencies with elevation, precipitation during the warmest periods, and seasonality in the thermal regime. For D. gliroides, genetic variation appeared to be associated with more variables than D. bozinovici, including precipitation and temperature-related variables, isothermality, and elevation. All the outlier SNPs were mapped to the D. gliroides reference genome to explore if they fell within functionally known genes. These results represent a necessary first step toward identifying the genome regions that harbor genes associated with climate adaptations in Dromiciops. Notably, we identified genes involved in various functions, including carbohydrate synthesis (ALG8), muscle and neuronal regulation (MEF2D), and stress responses (PTGES3). Ultimately, this study contributes valuable insights that can inform targeted conservation strategies aimed at preserving the genetic diversity of Dromiciops in the face of environmental challenges.

Genomic signatures of climate adaptation in bank voles

Evidence for divergent selection and adaptive variation across the landscape can provide insight into a species' ability to adapt to different environments. However, despite recent advances in genomics, it remains difficult to detect the footprints of climate-mediated selection in natural populations. Here, we analysed ddRAD sequencing data (21,892 SNPs) in conjunction with geographic climate variation to search for signatures of adaptive differentiation in twelve populations of the bank vole (Clethrionomys glareolus) distributed across Europe. To identify the loci subject to selection associated with climate variation, we applied multiple genotype-environment association methods, two univariate and one multivariate, and controlled for the effect of population structure. In total, we identified 213 candidate loci for adaptation, 74 of which were located within genes. In particular, we identified signatures of selection in candidate genes with functions related to lipid metabolism and the immune system. Using the results of redundancy analysis, we demonstrated that population history and climate have joint effects on the genetic variation in the pan-European metapopulation. Furthermore, by examining only candidate loci, we found that annual mean temperature is an important factor shaping adaptive genetic variation in the bank vole. By combining landscape genomic approaches, our study sheds light on genome-wide adaptive differentiation and the spatial distribution of variants underlying adaptive variation influenced by local climate in bank voles.

Amazonian colonization from the Atlantic Forest: New perspectives on the connections of South American tropical forests

An open and dry vegetation belt separates Amazonia (AM) and the Atlantic Forest (AF). Evidence from palaeoclimatic and phylogenetic studies suggests past connections between these forests during cycles of increased humidity through the formation of forest corridors. The distinctive northern AF avifauna is known to have affinities both with AM and the southern AF. Still, the extent of how these two regions contributed to the assemblage of this avifauna remains poorly understood. Using historical demographic analyses and comparative phylogeography based on sub-genomic genetic sampling, we assessed how past connections between AM and AF led to shared vicariance and colonization events in four avian AF endemic taxa. Our results supported the occurrence of humid forest corridors promoting the contact between AF and AM populations and suggested two vicariant events and two colonization events from AF to AM. Population divergences were mostly non-synchronous and occurred multiple times during the Pleistocene. Historical gene flow was prevalent across study groups, supporting migration flows after the initial separation between AM and AF - a pattern previously unknown in birds between these regions. Idiosyncratic histories and divergent demographic syndromes suggest that organisms' responses to climate-driven habitat shifts broadly depend on their ecological attributes. This study strengthened our knowledge of past connections between AM and AF and provided demographic scenarios amenable for testing in other groups of co-distributed organisms.

Evidence of introgression, ecological divergence and adaptation in Asterias sea stars

Hybrid zones are important windows into the evolutionary dynamics of populations, revealing how processes like introgression and adaptation structure population genomic variation. Importantly, they are useful for understanding speciation and how species respond to their environments. Here, we investigate two closely related sea star species, Asterias rubens and A. forbesi, distributed along rocky European and North American coastlines of the North Atlantic, and use genome-wide molecular markers to infer the distribution of genomic variation within and between species in this group. Using genomic data and environmental niche modelling, we document hybridization occurring between northern New England and the southern Canadian Maritimes. We investigate the factors that maintain this hybrid zone, as well as the environmental variables that putatively drive selection within and between species. We find that the two species differ in their environmental niche breadth; Asterias forbesi displays a relatively narrow environmental niche while conversely, A. rubens has a wider niche breadth. Species distribution models accurately predict hybrids to occur within environmental niche overlap, thereby suggesting environmental selection plays an important role in the maintenance of the hybrid zone. Our results imply that the distribution of genomic variation in North Atlantic sea stars is influenced by the environment, which will be crucial to consider as the climate changes.

Testing assertions of widespread introgressive hybridization in a clade of neotropical toads with low mate selectivity (Rhinella granulosa species group)

Discordance between different genomic regions, often identified through multilocus sequencing of selected markers, presents particular difficulties in identifying historical processes which drive species diversity and boundaries. Mechanisms causing discordance, such as incomplete lineage sorting or introgression due to interspecific hybridization, are better identified based on population-level genomic datasets. In the toads of the Rhinella granulosa species group, patterns of mito-nuclear discordance and potential hybridization have been reported by several studies. However, these patterns were proposed based on few loci, such that alternative mechanisms behind gene-tree heterogeneity cannot be ruled out. Using genome-wide ddRADseq loci from a subset of species within this clade, we found only partial concordance between currently recognized species-level taxon boundaries and patterns of genetic structure. While most taxa within the R. granulosa group correspond to clades, genetic clustering analyses sometimes grouped distinct taxonomic units into a single cluster. Moreover, levels of admixture between inferred clusters were limited and restricted to a single taxon pair which is best explained by incomplete lineage sorting as opposed to introgressive hybridization, according to D-statistics results. These findings contradict previous assertions of widespread cryptic diversity and gene flow within the R. granulosa clade. Lastly, our analyses suggest that diversification events within the Rhinella granulosa group mostly dated back to the early Pliocene, being generally younger than species divergences in other closely related clades that present high levels of cross-species gene flow. This finding uniquely contradicts common assertions that this young clade of toads exhibits interspecific hybridization.

Population genomic analyses support sympatric origins of parapatric morphs in a salamander

In numerous clades, divergent sister species have largely non-overlapping geographic ranges. This pattern presumably arises because species diverged in allopatry or parapatry, prior to a subsequent contact. Here, we provide population-genomic evidence for the opposite scenario: previously sympatric ecotypes that have spatially separated into divergent monomorphic populations over large geographic scales (reverse sympatric scenario). We analyzed a North American salamander (Plethodon cinereus) with two color morphs that are broadly sympatric: striped (redback) and unstriped (leadback). Sympatric morphs can show considerable divergence in other traits, and many Plethodon species are fixed for a single morph. Long Island (New York) is unusual in having many pure redback and leadback populations that are spatially separated, with pure redback populations in the west and pure leadbacks in the east. Previous work showed that these pure-morph populations were genetically, morphologically, and ecologically divergent. Here, we performed a coalescent-based analysis of new data from 88,696 single-nucleotide polymorphisms to address the origins of these populations. This analysis strongly supports the monophyly of Long Island populations and their subsequent divergence into pure redback and pure leadback populations. Taken together, these results suggest that the formerly sympatric mainland morphs separated into parapatric populations on Long Island, reversing the conventional speciation scenario.

Population genetic differentiation and genomic signatures of adaptation to climate in an abundant lizard

Species distributed across climatic gradients will typically experience spatial variation in selection, but gene flow can prevent such selection from causing population genetic differentiation and local adaptation. Here, we studied genomic variation of 415 individuals across 34 populations of the common wall lizard (Podarcis muralis) in central Italy. This species is highly abundant throughout this region and populations belong to a single genetic lineage, yet there is extensive phenotypic variation across climatic regimes. We used redundancy analysis to, first, quantify the effect of climate and geography on population genomic variation in this region and, second, to test if climate consistently sorts specific alleles across the landscape. Climate explained 5% of the population genomic variation across the landscape, about half of which was collinear with geography. Linear models and redundancy analyses identified loci that were significantly differentiated across climatic regimes. These loci were distributed across the genome and physically associated with genes putatively involved in thermal tolerance, regulation of temperature-dependent metabolism and reproductive activity, and body colouration. Together, these findings suggest that climate can exercise sufficient selection in lizards to promote genetic differentiation across the landscape in spite of high gene flow.

Genome-scale data reveal deep lineage divergence and a complex demographic history in the Texas horned lizard (Phrynosoma cornutum) throughout the southwestern and central US

The southwestern and central United States serve as an ideal region to test alternative hypotheses regarding biotic diversification. Genomic data can now be combined with sophisticated computational models to quantify the impacts of paleoclimate change, geographic features, and habitat heterogeneity on spatial patterns of genetic diversity. In this study, we combine thousands of genotyping-by-sequencing (GBS) loci with mtDNA sequences (ND1) from the Texas horned lizard (Phrynosoma cornutum) to quantify relative support for different catalysts of diversification. Phylogenetic and clustering analyses of the GBS data indicate support for at least three primary populations. The spatial distribution of populations appears concordant with habitat type, with desert populations in AZ and NM showing the largest genetic divergence from the remaining populations. The mtDNA data also support a divergent desert population, but other relationships differ and suggest mtDNA introgression. Genotype–environment association with bioclimatic variables supports divergence along precipitation gradients more than along temperature gradients. Demographic analyses support a complex history, with introgression and gene flow playing an important role during diversification. Bayesian multispecies coalescent analyses with introgression (MSci) analyses also suggest that gene flow occurred between populations. Paleo-species distribution models support two southern refugia that geographically correspond to contemporary lineages. We find that divergence times are underestimated and population sizes are overestimated when introgression occurred and is ignored in coalescent analyses, and furthermore, inference of ancient introgression events and demographic history is sensitive to inclusion of a single recently admixed sample. Our analyses cannot refute the riverine barrier or glacial refugia hypotheses. Results also suggest that populations are continuing to diverge along habitat gradients. Finally, the strong evidence of admixture, gene flow, and mtDNA introgression among populations suggests that P. cornutum should be considered a single widespread species under the General Lineage Species Concept.

Distinct interspecific and intraspecific vulnerability of coastal species to global change

Characterising and predicting species responses to anthropogenic global change is one of the key challenges in contemporary ecology and conservation. The sensitivity of marine species to climate change is increasingly being described with forecasted species distributions, yet these rarely account for population level processes such as genomic variation and local adaptation. This study compares inter- and intraspecific patterns of biological composition to determine how vulnerability to climate change, and its environmental drivers, vary across species and populations. We compare species trajectories for three ecologically important southern African marine invertebrates at two time points in the future, both at the species level, with correlative species distribution models, and at the population level, with gradient forest models. Reported range shifts are species-specific and include both predicted range gains and losses. Forecasted species responses to climate change are strongly influenced by changes in a suite of environmental variables, from sea surface salinity and sea surface temperature, to minimum air temperature. Our results further suggest a mismatch between future habitat suitability (where species can remain in their ecological niche) and genomic vulnerability (where populations retain their genomic composition), highlighting the inter- and intraspecific variability in species' sensitivity to global change. Overall, this study demonstrates the importance of considering species and population level climatic vulnerability when proactively managing coastal marine ecosystems in the Anthropocene.

The effects of climate and demographic history in shaping genomic variation across populations of the Desert Horned Lizard (Phrynosoma platyrhinos)

Species often experience spatial environmental heterogeneity across their range, and populations may exhibit signatures of adaptation to local environmental characteristics. Other population genetic processes, such as migration and genetic drift, can impede the effects of local adaptation. Genetic drift in particular can have a pronounced effect on population genetic structure during large-scale geographic expansions, where a series of founder effects leads to decreases in genetic variation in the direction of the expansion. Here, we explore the genetic diversity of a desert lizard that occupies a wide range of environmental conditions and that has experienced post-glacial expansion northwards along two colonization routes. Based on our analyses of a large SNP data set, we find evidence that both climate and demographic history have shaped the genetic structure of populations. Pronounced genetic differentiation was evident between populations occupying cold versus hot deserts, and we detected numerous loci with significant associations with climate. The genetic signal of founder effects, however, is still present in the genomes of the recently expanded populations, which comprise subsets of genetic variation found in the southern populations.

Speciation and secondary contact in a fossorial island endemic, the São Tomé caecilian

A period of isolation in allopatry typically precedes local adaptation and subsequent divergence among lineages. Alternatively, locally adapted phenotypes may arise and persist in the face of gene flow, resulting in strong correlations between ecologically-relevant phenotypic variation and corresponding environmental gradients. Quantifying genetic, ecological, and phenotypic divergence in such lineages can provide insights into the abiotic and biotic mechanisms that structure populations and drive the accumulation of phenotypic and taxonomic diversity. Low-vagility organisms whose distributions span ephemeral geographic barriers present the ideal evolutionary context within which to address these questions. Here, we combine genetic (mtDNA and genome-wide SNPs) and phenotypic data to investigate the divergence history of caecilians (Amphibia: Gymnophiona) endemic to the oceanic island of São Tomé in the Gulf of Guinea archipelago. Consistent with a previous mtDNA study, we find two phenotypically and genetically distinct lineages that occur along a north-to-south axis with extensive admixture in the centre of the island. Demographic modelling supports divergence in allopatry (~300 kya) followed by secondary contact (~95 kya). Consequently, in contrast to a morphological study that interpreted latitudinal phenotypic variation in these caecilians as a cline within a single widespread species, our analyses suggest a history of allopatric lineage divergence and subsequent hybridization that may have blurred species boundaries. We propose that late Pleistocene volcanic activity favoured allopatric divergence between these lineages with local adaptation to climate maintaining a stable hybrid zone in the centre of São Tomé Island. Our study joins a growing number of systems demonstrating lineage divergence on volcanic islands with stark environmental transitions across small geographic distances.

Evolutionary drivers of sexual signal variation in Amazon Slender Anoles

Phenotypic variation among populations, as seen in the signaling traits of many species, provides an opportunity to test whether similar factors generate repeated phenotypic patterns in different parts of a species' range. We investigated whether genetic divergence, abiotic gradients, and sympatry with closely related species explain variation in the dewlap colors of Amazon Slender Anoles, Anolis fuscoauratus. To this aim, we characterized dewlap diversity in the field with respect to population genetic structure and evolutionary relationships, assessed whether dewlap phenotypes are associated with climate or landscape variables, and tested for nonrandom associations in the distributions of A. fuscoauratus phenotypes and sympatric Anolis species. We found that dewlap colors vary among but not within sites in A. fuscoauratus. Regional genetic clusters included multiple phenotypes, while populations with similar dewlaps were often distantly related. Phenotypes did not segregate in environmental space, providing no support for optimized signal transmission at a local scale. Instead, we found a negative association between certain phenotypes and sympatric Anolis species with similar dewlap color attributes, suggesting that interactions with closely related species promoted dewlap divergence among A. fuscoauratus populations. Amazon Slender Anoles emerge as a promising system to address questions about parallel trait evolution and the contribution of signaling traits to speciation.

Genomic signatures of drift and selection driven by predation and human pressure in an insular lizard

Genomic divergence was studied in 10 small insular populations of the endangered Balearic Islands lizard (Podarcis lilfordi) using double digest restriction-site associated DNA sequencing. The objectives were to establish levels of divergence among populations, investigate the impact of population size on genetic variability and to evaluate the role of different environmental factors on local adaptation. Analyses of 72,846 SNPs supported a highly differentiated genetic structure, being the populations with the lowest population size (Porros, Foradada and Esclatasang islets) the most divergent, indicative of greater genetic drift. Outlier tests identified ~ 2% of loci as candidates for selection. Genomic divergence-Enviroment Association analyses were performed using redundancy analyses based on SNPs putatively under selection, detecting predation and human pressure as the environmental variables with the greatest explanatory power. Geographical distributions of populations and environmental factors appear to be fundamental drivers of divergence. These results support the combined role of genetic drift and divergent selection in shaping the genetic structure of these endemic island lizard populations.

Multi-model seascape genomics identifies distinct environmental drivers of selection among sympatric marine species

BACKGROUND

As global change and anthropogenic pressures continue to increase, conservation and management increasingly needs to consider species' potential to adapt to novel environmental conditions. Therefore, it is imperative to characterise the main selective forces acting on ecosystems, and how these may influence the evolutionary potential of populations and species. Using a multi-model seascape genomics approach, we compare putative environmental drivers of selection in three sympatric southern African marine invertebrates with contrasting ecology and life histories: Cape urchin (Parechinus angulosus), Common shore crab (Cyclograpsus punctatus), and Granular limpet (Scutellastra granularis).

RESULTS

Using pooled (Pool-seq), restriction-site associated DNA sequencing (RAD-seq), and seven outlier detection methods, we characterise genomic variation between populations along a strong biogeographical gradient. Of the three species, only S. granularis showed significant isolation-by-distance, and isolation-by-environment driven by sea surface temperatures (SST). In contrast, sea surface salinity (SSS) and range in air temperature correlated more strongly with genomic variation in C. punctatus and P. angulosus. Differences were also found in genomic structuring between the three species, with outlier loci contributing to two clusters in the East and West Coasts for S. granularis and P. angulosus, but not for C. punctatus.

CONCLUSION

The findings illustrate distinct evolutionary potential across species, suggesting that species-specific habitat requirements and responses to environmental stresses may be better predictors of evolutionary patterns than the strong environmental gradients within the region. We also found large discrepancies between outlier detection methodologies, and thus offer a novel multi-model approach to identifying the principal environmental selection forces acting on species. Overall, this work highlights how adding a comparative approach to seascape genomics (both with multiple models and species) can elucidate the intricate evolutionary responses of ecosystems to global change.

Molecular phylogenetic inference of the howler monkey radiation (Primates: Alouatta)

Howler monkeys (Alouatta), comprising between nine and 14 species and ranging from southern Mexico to northern Argentina, are the most widely distributed platyrrhines. Previous phylogenetic studies of howlers have used chromosomal and morphological characters and a limited number of molecular markers; however, branching patterns conflict between studies or remain unresolved. We performed a new phylogenetic analysis of Alouatta using both concatenated and coalescent-based species tree approaches based on 14 unlinked non-coding intergenic nuclear regions. Our taxon sampling included five of the seven South American species (Alouatta caraya, Alouatta belzebul, Alouatta guariba, Alouatta seniculus, Alouatta sara) and the two recognized species from Mesoamerica (Alouatta pigra, Alouatta palliata). Similarly to previous studies, our phylogenies supported a Mesoamerican clade and a South American clade. For the South American howlers, both methods recovered the Atlantic Forest endemic A. guariba as sister to all remaining South American species, albeit with moderate support. Moreover, we found no support for the previously proposed sister relationship between A. guariba and A. belzebul. For the first time, a clade composed of A. sara and A. caraya was identified. The relationships among the other South American howlers, however, were not fully supported. Our estimates for divergence times within Alouatta are generally older compared to estimates in earlier studies. However, they conform to recent studies proposing a Miocene age for the Isthmus of Panama and for the uplift of the northern Andes. Our results also point to an early genetic isolation of A. guariba in the Atlantic Forest, in agreement with the hypothesis of biotic exchange across South American rain forests in the Miocene. Collectively, these findings contribute to a better understanding of the diversification processes among howler monkey species; however, they also suggest that further comprehension of the evolutionary history of the Alouatta radiation will rely on broadened taxonomic, geographic, and genomic sampling.

Paleoclimatic evolution as the main driver of current genomic diversity in the widespread and polymorphic Neotropical songbird Arremon taciturnus

Several factors have been proposed as drivers of species diversification in the Neotropics, including environmental heterogeneity, the development of drainage systems and historical changes in forest distribution due to climatic oscillations. Here, we investigate which drivers contributed to the evolutionary history and current patterns of diversity of a polymorphic songbird (Arremon taciturnus) that is widely distributed in Amazonian and Atlantic forests as well as in Cerrado gallery and seasonally-dry forests. We use genomic, phenotypic and habitat heterogeneity data coupled with climatic niche modelling. Results suggest the evolutionary history of the species is mainly related to paleoclimatic changes, although changes in the strength of the Amazon river as a barrier to dispersal, current habitat heterogeneity and geographic distance were also relevant. We propose an ancestral distribution in the Guyana Shield, and recent colonization of areas south of the Amazon river at ~380 to 166 kya, and expansion of the distribution to southern Amazonia, Cerrado and the Atlantic Forest. Since then, populations south of the Amazon River have been subjected to cycles of isolation and possibly secondary contact due to climatic changes that affected habitat heterogeneity and population connectivity. Most Amazonian rivers are not associated with long lasting isolation of populations, but some might act as secondary barriers, susceptible to crossing under specific climatic conditions. Morphological variation, while stable in some parts of the distribution, is not a reliable indicator of genetic structure or phylogenetic relationships.

Diversification history of clown tree frogs in Neotropical rainforests (Anura, Hylidae, Dendropsophus leucophyllatus group)

General consensus emphasizes that no single biological process can explain the patterns of species' distributions and diversification in the Neotropics. Instead, the interplay of several processes across space and time must be taken into account. Here we investigated the phylogenetic relationships and biogeographic history of tree frogs in the Dendropsophus leucophyllatus species group (Amphibia: Hylidae), which is distributed across Amazonia and the Atlantic rainforests. Using Next Generation Sequencing (NGS) and double digest restriction-site associated DNA (ddRADseq), we inferred phylogenetic relationships, species limits, and temporal and geographic patterns of diversification relative to the history of these biomes. Our results indicate that the D. leucophyllatus species group includes at least 14 independent lineages, which are currently arranged into ten described species. Therefore, a significant portion of species in the group are still unnamed. Different processes were associated to the group diversification history. For instance, the Andes uplift likely caused allopatric speciation for Cis-Andean species, whereas it may also be responsible for changes in the Amazonian landscape triggering parapatric speciation by local adaptation to ecological factors. Meanwhile, Atlantic Forest ancestors unable to cross the dry diagonal biomes after rainforest's retraction, evolved in isolation into different species. Diversification in the group began in the early Miocene, when connections between Atlantic Forest and the Andes (Pacific Dominion) by way of a south corridor were possible. The historical scenario in Amazonia, characterized by several speciation events and habitat heterogeneity, helped promoting diversification, resulting in the highest species diversity for the group. This marked species diversification did not happen in Atlantic Forest, where speciation is very recent (late Pliocene and Pleistocene), despite its remarkable climatic heterogeneity.

Selection at behavioural, developmental and metabolic genes is associated with the northward expansion of a successful tropical colonizer

What makes a species able to colonize novel environments? This question is key to understand the dynamics of adaptive radiations and ecological niche shifts, but the mechanisms that underlie expansion into novel habitats remain poorly understood at a genomic scale. Lizards from the genus Anolis are typically tropical, and the green anole (Anolis carolinensis) constitutes an exception since it expanded into temperate North America from subtropical Florida. Thus, we used the green anole as a model to investigate signatures of selection associated with colonization of a new environment, namely temperate North America. To this end, we analysed 29 whole-genome sequences, covering the entire native range of the species. We used a combination of recent methods to quantify both positive and balancing selection in northern populations, including F_ST outlier methods, machine learning and ancestral recombination graphs. We naively scanned for genes of interest and assessed the overlap between multiple tests. Strikingly, we identified many genes involved in behaviour, suggesting that the recent successful colonization of northern environments may have been linked to behavioural shifts as well as physiological adaptation. Using a candidate genes strategy, we determined that genes involved in response to cold or behaviour displayed more frequently signals of selection, while controlling for local recombination rate, gene clustering and gene length. In addition, we found signatures of balancing selection at immune genes in all investigated genetic groups, but also at genes involved in neuronal and anatomical development.