Local adaptation to environmental heterogeneity across a landscape can result in population divergence and formation of lineages. On Guadeloupe Island, the active volcano, La Grande Soufri & egrave;re, peaks at 1460 m a.s.l., with rainforest at low elevations transitioning to humid savannahs at high elevations. Two endemic sister species of Eleutherodactylus frogs are co-distributed across this habitat gradient, and previous studies have reported phenotypic differences between lowland and high-elevation populations in each species, in addition to potential ongoing hybridization between the species at high-elevation sites. Here we generate mitochondrial DNA and nuclear DNA genomic data along the elevational transect to quantify population genetic structure, provide historical context for the diversification of these island endemics, and identify potential bottlenecks attributable to the eruptive history of the volcano. We find that both taxa exhibit population clusters that correspond to low- and high-elevation localities; however, genetic divergence is not associated with climate variables or geographical distance. The timing of divergence between the species is estimated at similar to 3.75 Mya; demographic models indicate low levels of migration between the species after divergence, and we find that ongoing hybridization is likely to be limited. Finally, we find moderate heterozygosity across populations, suggesting that they were minimally impacted by recent volcanic activity. A version of this abstract translated to French can be found in the Supplementary Data. Une version de ce resume traduite en fran & ccedil;ais est disponible dans les Donnees Supplementaires.

The introduction of invasive species to islands can be devastating to local biota. While many invasive species’successful colonization of new habitats can be attributed to features like high fecundity or low mate selectivity,parthenogenetic species pose a unique threat in their ability to reproduce asexually at high rates. The Smooth-scaledTegulet (Gymnophthalmus underwoodi), first described from Barbados in the 1950s, has since been documented as anintroduced species across additional Lesser and Greater Antillean islands as well as northwestern Brazil and Venezuela.Herein we report multiple individuals from Puerto Rico, including the first genetic records, which indicate that theindividuals on this island are most closely related to populations recorded from Brazil and Montserrat. Additionalgenetic information is necessary to determine the timing and migration of this species across the Caribbean. Regardless,the continued geographic expansion of this parthenogen across the Caribbean could present a threat to local fauna.

New World mabuyine skinks are a diverse radiation of morphologically cryptic lizards with unique reproductive biologies. Recent studies examining population-level data (morphological, ecological, and genomic) have uncovered novel biodiversity and phenotypes, including the description of dozens of new species and insights into the evolution of their highly complex placental structures. Beyond the potential for this diverse group to serve as a model for the evolution of viviparity in lizards, much of the taxonomic diversity is concentrated in regions experiencing increasing environmental instability from climate and anthropogenic change. Consequently, a better understanding of genome structure and diversity will be an important tool in the adaptive management and conservation of this group. Skinks endemic to Caribbean islands are particularly vulnerable to global change with several species already considered likely extinct and several remaining species either endangered or threatened. Combining PacBio long-read sequencing, Hi-C, and RNAseq data, here we present the first genomic resources for this group by describing new chromosome-level reference genomes for the Puerto Rican Skink Spondylurus nitidus and the Culebra Skink S. culebrae. Results indicate two high quality genomes, both ∼1.4 Gb, assembled nearly telomere to telomere with complete mitochondrion assembly and annotation.

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.

The effects of genetic introgression on species boundaries and how they affect species' integrity and persistence over evolutionary time have received increased attention. The increasing availability of genomic data has revealed contrasting patterns of gene flow across genomic regions, which impose challenges to inferences of evolutionary relationships and of patterns of genetic admixture across lineages. By characterizing patterns of variation across thousands of genomic loci in a widespread complex of true toads (Rhinella), we assess the true extent of genetic introgression across species thought to hybridize to extreme degrees based on natural history observations and multilocus analyses. Comprehensive geographic sampling of five large-ranged Neotropical taxa revealed multiple distinct evolutionary lineages that span large geographic areas and, at times, distinct biomes. The inferred major clades and genetic clusters largely correspond to currently recognized taxa; however, we also found evidence of cryptic diversity within taxa. While previous phylogenetic studies revealed extensive mitonuclear discordance, our genetic clustering analyses uncovered several admixed individuals within major genetic groups. Accordingly, historical demographic analyses supported that the evolutionary history of these toads involved cross-taxon gene flow both at ancient and recent times. Lastly, ABBA-BABA tests revealed widespread allele sharing across species boundaries, a pattern that can be confidently attributed to genetic introgression as opposed to incomplete lineage sorting. These results confirm previous assertions that the evolutionary history of Rhinella was characterized by various levels of hybridization even across environmentally heterogeneous regions, posing exciting questions about what factors prevent complete fusion of diverging yet highly interdependent evolutionary lineages.

Knowledge of how contemporary and historical factors drive patterns of genetic structure across geographic space can shed light on the processes underlying diversification. This approach is especially fruitful in studies of widespread species or species clades that occur across multiple environmental conditions and biomes. In the Neotropics, specifically, molecular data from widespread vertebrate species have revealed high levels of lineage diversity and spatial genetic structure - yet studies that explore the possible correlates of local structure patterns are lacking. We investigate the distribution of lineage diversity within two widespread South American skink species complexes of the genus Mabuya. We characterize genetic structure and diversity in these widely ranged lizards, and identify potential geographic and environmental correlates, to shed light on the processes that promote lineage diversification across the heterogeneous landscapes which they occupy. In both groups, we found mitochondrial lineages to be spatially structured along the coastal forests and the savannas of Brazil. These mtDNA patterns are, however, not shared with those inferred from nuclear DNA markers. The geographic location of major mitochondrial genetic breaks is consistent with those of other taxa, suggesting common responses to former landscape change in eastern South America, particularly along geological faults. Genetic differentiation is correlated with environmental turnover and geographic separation in one, but not in the other, group of skinks. Compared to other studies of similarly widely distributed organisms, the link between spatial environmental gradients and genetic differentiation is not as strong or consistent, suggesting a more complex history underlying current phylogeographic patterns. Our genetic data indicate the existence of yet undescribed diversity in wide-ranging lizards, and the value of phylogenetic and phylogeographic studies of similarly understudied species.

Aim: Species with wide distributions spanning the African Guinean and Congolian rain forests are often composed of genetically distinct populations or cryptic species with geographic distributions that mirror the locations of the remaining forest habitats. We used phylogeographic inference and demographic model testing to evaluate diversification models in a widespread rain forest species, the African foam-nest treefrog Chiromantis rufescens. Location: Guinean and Congolian rain forests, West and Central Africa. Taxon: Chiromantis rufescens. Methods: We collected mitochondrial DNA (mtDNA) and single-nucleotide polymorphism (SNP) data for 130 samples of C. rufescens. After estimating population structure and inferring species trees using coalescent methods, we tested demographic models to evaluate alternative population divergence histories that varied with respect to gene flow, population size change and periods of isolation and secondary contact. Species distribution models were used to identify the regions of climatic stability that could have served as forest refugia since the last interglacial. Results: Population structure within C. rufescens resembles the major biogeographic regions of the Guinean and Congolian forests. Coalescent-based phylogenetic analyses provide strong support for an early divergence between the western Upper Guinean forest and the remaining populations. Demographic inferences support diversification models with gene flow and population size changes even in cases where contemporary populations are currently allopatric, which provides support for forest refugia and barrier models. Species distribution models suggest that forest refugia were available for each of the populations throughout the Pleistocene. Main conclusions: Considering historical demography is essential for understanding population diversification, especially in complex landscapes such as those found in the Guineo-Congolian forest. Population demographic inferences help connect the patterns of genetic variation to diversification model predictions. The diversification history of C. rufescens was shaped by a variety of processes, including vicariance from river barriers, forest fragmentation and adaptive evolution along environmental gradients.

The accumulation of biodiversity in tropical forests can occur through multiple allopatric and parapatric models of diversification, including forest refugia, riverine barriers and ecological gradients. Considerable debate surrounds the major diversification process, particularly in the West African Lower Guinea forests, which contain a complex geographic arrangement of topographic features and historical refugia. We used genomic data to investigate alternative mechanisms of diversification in the Gaboon forest frog, Scotobleps gabonicus, by first identifying population structure and then performing demographic model selection and spatially explicit analyses. We found that a majority of population divergences are best explained by allopatric models consistent with the forest refugia hypothesis and involve divergence in isolation with subsequent expansion and gene flow. These population divergences occurred simultaneously and conform to predictions based on climatically stable regions inferred through ecological niche modelling. Although forest refugia played a prominent role in the intraspecific diversification of S.gabonicus, we also find evidence for potential interactions between landscape features and historical refugia, including major rivers and elevational barriers such as the Cameroonian Volcanic Line. We outline the advantages of using genomewide variation in a model-testing framework to distinguish between alternative allopatric hypotheses, and the pitfalls of limited geographic and molecular sampling. Although phylogeographic patterns are often species-specific and related to life-history traits, additional comparative studies incorporating genomic data are necessary for separating shared historical processes from idiosyncratic responses to environmental, climatic and geological influences on diversification.

Shifts in the geographic distribution of habitats over time can promote dispersal and vicariance, thereby influencing large‐scale biogeographic patterns and ecological processes. An example is that of transient corridors of suitable habitat across disjunct but ecologically similar regions, which have been associated with climate change over time. Such connections likely played a role in the assembly of tropical communities, especially within the highly diverse Amazonian and Atlantic rainforests of South America. Although these forests are presently separated by open and dry ecosystems, paleoclimatic and phylogenetic evidence suggest that they have been transiently connected in the past. However, little is known about the timing, magnitude and the distribution of former forest connections. We employ sequence data at multiple loci from three codistributed arboreal lizards (Anolis punctatus, Anolis ortonii and Polychrus marmoratus) to infer the phylogenetic relationships among Amazonian and Atlantic Forest populations and to test alternative historical demographic scenarios of colonization and vicariance using coalescent simulations and approximate Bayesian computation (ABC). Data from the better‐sampled Anolis species support colonization of the Atlantic Forest from eastern Amazonia. Hierarchical ABC indicates that the three species colonized the Atlantic Forest synchronously during the mid‐Pleistocene. We find support of population bottlenecks associated with founder events in the two Anolis, but not in P. marmoratus, consistently with their distinct ecological tolerances. Our findings support that climatic fluctuations provided key opportunities for dispersal and forest colonization in eastern South America through the cessation of environmental barriers. Evidence of species‐specific histories strengthens assertions that biological attributes play a role in responses to shared environmental change.

Phylogeographic endemism, the degree to which the history of recently evolved lineages is spatially restricted, reflects fundamental evolutionary processes such as cryptic divergence, adaptation and biological responses to environmental heterogeneity. Attempts to explain the extraordinary diversity of the tropics, which often includes deep phylogeographic structure, frequently invoke interactions of climate variability across space, time and topography. To evaluate historical versus contemporary drivers of phylogeographic endemism in a tropical system, we analyse the effects of current and past climatic variation on the genetic diversity of 25 vertebrates in the Brazilian Atlantic rainforest. We identify two divergent bioclimatic domains within the forest and high turnover around the Rio Doce. Independent modelling of these domains demonstrates that endemism patterns are subject to different climatic drivers. Past climate dynamics, specifically areas of relative stability, predict phylogeographic endemism in the north. Conversely, contemporary climatic heterogeneity better explains endemism in the south. These results accord with recent speleothem and fossil pollen studies, suggesting that climatic variability through the last 250 kyr impacted the northern and the southern forests differently. Incorporating sub-regional differences in climate dynamics will enhance our ability to understand those processes shaping high phylogeographic and species endemism, in the Neotropics and beyond.

3 Introduction 4 Materials & Methods 8 Sampling and Molecular Data 8 Phylogenetic and STRUCTURE Analyses 8 Environmental Analyses 10 Multiple Matrix Regression with Randomization (MMRR) 10 Results 12 Phylogenetic and STRUCTURE Analyses 12 Environmental Analyses and MMRR 14 Discussion 16 Acknowledgements 18 References 19 Figures 35 Tables 41 Supplementary Material 42 Methods 42 Figures 43 Tables 47