Seascape genomics of the pink abalone (Haliotis corrugata): An insight into a cross-border species in the northeast Pacific coast.

Seascape genomics gives insight into the geographic and environmental factors shaping local adaptations. It improves the understanding of the potential effects of climate change, which is relevant to provide the basis for the international management of fishery resources. The pink abalone (Haliotis corrugata) is distributed from California, United States to Baja California Sur, Mexico, exposed to a latitudinal environmental gradient in the California Current System. Management of the pink abalone contrasts between Mexico and the United States; Mexico has an active fishery organized in four administrative areas, while the United States has kept the fishery in permanent closure since 1996. However, the impact of environmental factors on genetic variation along the species distribution remains unknown, and understanding this relationship is crucial for effective spatial management strategies. This study aims to investigate the neutral and adaptive genomic structure of H. corrugata. A total of 203 samples from 13 locations were processed using ddRADseq, and covering the species' distribution. Overall, 2,231 neutral, nine potentially adaptive and three genomic-environmental association loci were detected. The neutral structure identified two groups: 1) California, United States and 2) Baja California Peninsula, México. In addition, the adaptive structure analysis also detected two groups with genetic divergence observed at Punta Eugenia. Notably, the seawater temperature significantly correlated with the northern group (temperate) and the southern (warmer) group. This study is a valuable foundation for future research and conservation initiatives, emphasizing the importance of considering neutral and adaptive genetic factors when developing management strategies for marine species.

Differences in the genomic diversity, structure, and inbreeding patterns in wild and managed populations of Agave potatorum Zucc. used in the production of Tobalá mezcal in Southern Mexico.

Agave potatorum Zucc. locally known as Tobalá, is an important species for mezcal production. It is a perennial species that takes 10 to 15 years to reach reproductive age. Because of high demand of Tobalá mezcal and the slow maturation of the plants, its wild populations have been under intense anthropogenic pressure. The main objective of this study was to estimate the genome-wide diversity in A. potatorum and determine if the type of management has had any effect on its diversity, inbreeding and structure. We analyzed 174 individuals (105 wild, 42 cultivated and 27 from nurseries) from 34 sites with a reduced representation genomic method (ddRADseq), using 14,875 SNPs. The diversity measured as expected heterozygosity was higher in the nursery and wild plants than in cultivated samples. We did not find private alleles in the cultivated and nursery plants, which indicates that the individuals under management recently derived from wild populations, which was supported by higher gene flow estimated from wild populations to the managed plants. We found low but positive levels of inbreeding (FIS = 0.082), probably related to isolation of the populations. We detected low genetic differentiation among populations (FST = 0.0796), with positive and significant isolation by distance. The population genetic structure in the species seems to be related to elevation and ecology, with higher gene flow among populations in less fragmented areas. We detected an outlier locus related to the recognition of pollen, which is also relevant to self-incompatibility protein (SI). Due to seed harvest and long generation time, the loss of diversity in A. potatorum has been gradual and artificial selection and incipient management have not yet caused drastic differences between cultivated and wild plants. Also, we described an agroecological alternative to the uncontrolled extraction of wild individuals.

Population genetic structure of the maize weevil, Sitophilus zeamais, in southern Mexico.

The maize weevil, Sitophilus zeamais, is a ubiquitous pest of maize and other cereal crops worldwide and remains a threat to food security in subsistence communities. Few population genetic studies have been conducted on the maize weevil, but those that exist have shown that there is very little genetic differentiation between geographically dispersed populations and that it is likely the species has experienced a recent range expansion within the last few hundred years. While the previous studies found little genetic structure, they relied primarily on mitochondrial and nuclear microsatellite markers for their analyses. It is possible that more fine-scaled population genetic structure exists due to local adaptation, the biological limits of natural species dispersal, and the isolated nature of subsistence farming communities. In contrast to previous studies, here, we utilized genome-wide single nucleotide polymorphism data to evaluate the genetic population structure of the maize weevil from the southern and coastal Mexican states of Oaxaca and Chiapas. We employed strict SNP filtering to manage large next generation sequencing lane effects and this study is the first to find fine-scale genetic population structure in the maize weevil. Here, we show that although there continues to be gene flow between populations of maize weevil, that fine-scale genetic structure exists. It is possible that this structure is shaped by local adaptation of the insects, the movement and trade of maize by humans in the region, geographic barriers to gene flow, or a combination of these factors.

Biological distances between Late Holocene hunter-gatherers from the lower Paraná wetland (Argentina).

INTRODUCTION: This study assesses the biological affinity between Late Holocene hunter-gatherers that inhabited the region of the lower Paraná wetland (Argentina) at the end of the Late Holocene (1800-700 years BP). METHODS: Cranial nonmetric traits were registered in 159 adult individuals of both sexes. Biodistances were estimated between environmental and archaeological units through Smith's Mean Measure of Divergence, Mahalanobis distance, multidimensional scaling, Fst index, and R-matrix. Biological distances concerning temporal and geographic distances between the archaeological sites were also evaluated through Mantel and partial Mantel tests. RESULTS: The environmental units showed statistically significant biodistances according to a distance spatial pattern but low genetic differentiation between groups. Additionally, some archaeological sites with median radiocarbon dates around 700 years BP also presented statistically significant biodistances. CONCLUSIONS: Gene flow and genetic drift could be responsible for the observed biodistance results. These processes gave rise to different subpopulations of hunter-gatherers with a shared material culture, whose territorial limits could have been reinforced around 700 years BP or earlier.

Landscape genetics across the Andes mountains: Environmental variation drives genetic divergence in the leaf-cutting ant Atta cephalotes.

Distinguishing among the mechanisms underlying the spatial distribution of genetic variation resulting from the environmental or physical barriers from those arising due to simple geographic distance is challenging in complex landscapes. The Andean uplift represents one of the most heterogeneous habitats where multiple mechanisms may interact, confounding their relative roles. We explore this broad question in the leaf-cutting ant Atta cephalotes, a species that is distributed across the Andes mountains, using nuclear microsatellite markers and mtCOI gene sequences. We investigate spatial genetic divergence across the western range of the northern Andes in Colombia by testing the relative role of alternative scenarios of population divergence, including isolation by geographic distance (IBD), climatic conditions (IBE), and the physical barriers presented by the Andes mountains (IBB). Our results reveal substantial genetic differentiation among A. cephalotes populations for both types of markers, but only nuclear divergence followed a hierarchical pattern with multiple models of genetic divergence imposed by the western range. Model selection showed that the IBD, IBE (temperature and precipitation), and IBB (Andes mountains) models, often proposed as individual drivers of genetic divergence, interact, and explain up to 33% of the genetic divergence in A. cephalotes. The IBE model remained significant after accounting for IBD, suggesting that environmental factors play a more prominent role than IBB. These factors, in combination with the idiosyncratic dispersal patterns of ants, appear to determine the hierarchical patterns of gene flow. This study enriches our understanding of the forces shaping population divergence in complex habitat landscapes.

Categories of variables in analysis of genetic diversity in S1 progenies of Psidium guajava.

Crossing and developing inbred lines have been promising options for guava breeding programs. The purpose of this study was to evaluate the genetic divergence among genotypes of S1 inbred guava families by means of the Gower's technique and the Ward-MLM methodology, to verify the correlation and relative contribution of traits, as well as to identify descriptors with minimum efficiency for this species. The experiment was implemented at the Estação Experimental da Ilha Barra do Pomba, in the municipality of Itaocara, RJ, Brazil. A randomized block design with 18 inbred families, three replicates, and ten plants per plot was used for the experimental design. After 19 months from the planting of the experiment, the 61 earliest and most productive genotypes (individual plants) were evaluated. For this purpose, 29 descriptors were evaluated, of which fifteen were qualitative and fourteen, quantitative. The characteristics required to obtain the distance matrix were analyzed based on the Gower algorithm, and a comparative cluster between the dendrograms of the morphoagronomic variables was achieved from this matrix. Lastly, the Ward-MLM procedure was applied to form the clusters of inbred families. By using all 29 descriptors, greater efficiency was achieved in cluster discrimination. Hence, according to the results identified, it is not possible to indicate minimum descriptors for the culture. Using the Ward-MLM method, the descriptors that most contributed to the divergence among the genotypes were fruit flesh mass, fruit weight, fruit diameter, fruit flesh thickness, fruit placental mass, and fruit length. The most divergent genotypes can be recommended for further crosses or self-pollinations to develop new lines in the guava breeding program of UENF.

Genetic divergence at species boundaries of the dolphinfish (Coryphaena hippurus) in the Tropical Eastern Pacific.

Background: Marine species constitute commercially important resources, and knowledge about mechanisms that shape phylogeographic patterns and genetic structure provides valuable information for conservation. The dolphinfish, Coryphaena hippurus, is one of the most important species caught in the Tropical Eastern Pacific (TEP). However, the lack of consensus about the existence of genetically differentiated populations in the area has hindered the adoption of management strategies to ensure its viability. Methods: We assessed genetic variation and phylogeographic structure using two mitochondrial genes and 14 nuclear DNA microsatellite loci. Population genetic tools were used to characterize the spatial distribution of genetic variation of C. hippurus in the TEP, evaluate the extent of connectivity between dolphinfish populations, infer potential barriers to gene flow, and test for signals of contemporary and historical demographic expansions. Results: Mitochondrial DNA sequences showed genetic homogeneity across locations in the TEP, as well as a strong signal of population expansion dated to the late Pleistocene. In contrast, nuclear microsatellite markers resolved four genetically distinct groups with a remarked genetic differentiation between the most distant locations, at the northern and southern boundaries of the species' range. High mean genetic diversity was found at all localities (Hs = 0.66-0.81). Notwithstanding, positive F IS and low effective population size (Ne = 77.9-496.4) were also recorded. Conclusions: The distribution of genetic variation could be related to expansion-contraction cycles following seasonal temperature changes at transitional areas, promoting population subdivisions. However, we cannot rule out the effect of oceanographic dynamics to the observed patterns. Although this marine species remains highly abundant despite commercial exploitation, the low Ne values are of conservation concern and must be considered in fishery management plans.

Population genomics of the neotropical palm Copernicia prunifera (Miller) H. E. Moore: Implications for conservation.

Copernicia prunifera (Miller) H. E. Moore is a palm tree native to Brazil. The products obtained from its leaf extracts are a source of income for local families and the agroindustry. Owing to the reduction of natural habitats and the absence of a sustainable management plan, the maintenance of the natural populations of this palm tree has been compromised. Therefore, this study aimed to evaluate the diversity and genetic structure of 14 C. prunifera populations using single nucleotide polymorphisms (SNPs) identified through genotyping-by-sequencing (GBS) to provide information that contributes to the conservation of this species. A total of 1,013 SNP markers were identified, of which 84 loci showed outlier behavior and may reflect responses to natural selection. Overall, the level of genomic diversity was compatible with the biological aspects of this species. The inbreeding coefficient (f) was negative for all populations, indicating excess heterozygotes. Most genetic variations occurred within populations (77.26%), and a positive correlation existed between genetic and geographic distances. The population structure evaluated through discriminant analysis of principal components (DAPC) revealed low genetic differentiation between populations. The results highlight the need for efforts to conserve C. prunifera as well as its distribution range to preserve its global genetic diversity and evolutionary potential.

A new species of Tricorythodes Ulmer from Colombia and additional notes on T. serratus (Baumgardner amp; vila) (Ephemeroptera: Leptohyphidae).

A new species of Tricorythodes is described and illustrated based on nymphs from Colombia. Legs, mouthparts and other parts of the body were photographed through SEM and molecular characterization was performed by amplification of the COImt gene. Tricorythodes carare sp. nov. shows close affinities with T. serratus (Baumgardner vila), both sharing a similar coloration pattern and the presence of chalazae on femora. The new species can be distinguished from T. serratus and all the species of the genus by the following combination of characters: maxillary palp 3-segmented, with apical setae; dorsal surface of median and hind femora with a longitudinal row of 5 to 9 chalazae serrated in the apex, tarsal claws with 4 or 5 marginal denticles and with 13 submarginal denticles; posterior margin of the abdominal terga with serrated chalazae in the apex. SEM images and distance analysis of the partial sequences of the COI gene were useful to discriminate the new species from T. serratus (genetic divergence 1719%). An additional undescribed species of this group is reported from Esmeraldas in Ecuador and for the departments of Sucre and Choc in Colombia, but it is not described here because of the low number of individuals and the lack of DNA integrity. What is known under the name Tricorythodes serratus is suspected to be a species complex.

Lack of genetic differentiation in yellowfin tuna has conservation implications in the Eastern Pacific Ocean.

Yellowfin tuna, Thunnus albacares, is an important global fishery and of particular importance in the Eastern Pacific Ocean (EPO). According to the 2019 Inter-American Tropical Tuna Commission (IATTC) assessment, yellowfin tuna within the EPO is a single stock, and is being managed as one stock. However, previous studies indicate site fidelity, or limited home ranges, of yellowfin tuna which suggests the potential for multiple yellowfin tuna stocks within the EPO, which was supported by a population genetic study using microsatellites. If numerous stocks are present, management at the wrong spatial scales could cause the loss of minor yellowfin tuna populations in the EPO. In this study we used double digestion RADseq to assess the genetic structure of yellowfin tuna in the EPO. A total of 164 yellowfin tuna from Cabo San Lucas, México, and the Galápagos Islands and Santa Elena, Ecuador, were analysed using 18,011 single nucleotide polymorphisms. Limited genetic differentiation (FST = 0.00058-0.00328) observed among the sampling locations (México, Ecuador, Peru, and within Ecuador) is consistent with presence of a single yellowfin tuna population within the EPO. Our findings are consistent with the IATTC assessment and provide further evidence of the need for transboundary cooperation for the successful management of this important fishery throughout the EPO.

New Insight into the Genome-Wide Diversity and Admixture of Six Colombian Sheep Populations.

Creole sheep represent a strategic genetic resource for populations living in marginal areas under financial restrictions on the American continent. Six Colombian sheep breeds (two wool (BCL-Boyacá and NCL-Nariño, 12 and 14 samples) and four hair (OPCE-Ethiopian, 54 samples; OPCS-Sudan, 74 samples; OPCP-Pelibeuy, 59 samples; OPCW-Wayúu, 24 samples) were genotyped using the Illumina Ovine SNP50 BeadChip. Data was also included from international 44 breeds from International Sheep Genomics Consortium (ISGC) and from data published in previous a previous work on the Caribbean and African breeds. Although geographically separated, wool (NCL, BCL) and hair types (OPCE, OPCS, OPCW) presented little genetic differentiation (FST 0.05) at a global level but several groups of animals separated suggesting local clustering due to geographical isolation. The OPCP underwent a recent crossing with Mexican Pelibuey, explaining its differentiation. Findings in this work such as the proximity to West African Djallonké (WAD) and Barbados Black Belly (BBB), suggest different introductions of African type animals from the Caribbean region on a pre-existing genetic basis formed by animals deriving from the first importations coming from Europe in colonial times. As expected, Colombian wool breeds showed, in particular in Admixture software results, a greater genomic component in common with European breeds and in particular with Iberian ones (Churra). This study provides a basis for future research into the genetic diversity within and between the Colombian sheep breeds analysed, and scientific data for policy decisions on Farm Animal Genetic Resources (FAnGR).

Two new species of geckos of the genus Phyllopezus Peters, 1878 (Squamata: Gekkota: Phyllodactylidae) from northeastern Brazil.

We describe two new species of Brazilian geckos of the genus Phyllopezus based on morphological and molecular data. The first species is currently known from a relictual Cerrado enclavecampos rupestres, in the mountains of the Serra do Espinhao in the state of Bahia. The second species is known from northeastern Atlantic Forest and transitional areas with Caatinga biome in the state of Alagoas. The two new species are sister taxa and together are the sister clade to the remaining species in the Phyllopezus pollicaris species complex. These new species can be morphologically distinguished from their congeners by meristic and morphometric characters, in addition to color pattern and genetic differentiation.

The taxonomic status of Cryptobatrachus frogs (Anura: Hemiphractidae) from the Serrana del Perij.

The Neotropical frog genus Cryptobatrachus includes six currently recognized species distributed throughout the northern Andes in Colombia and Venezuela. Cryptobatrachus conditus, C. pedroruizi, and C. remotus were described from the Colombian and Venezuelan slopes of the Serrana del Perij. Due to the great morphological similarity among these species, we re-assess their taxonomic status based on morphological and molecular data from types, topotypes, and specimens from localities referred to in the original descriptions of these species. Morphometric comparisons showed that all these species are virtually indistinguishable, although some subtle qualitative differences in morphological traits distinguish C. conditus from the other species. Phylogenetic analyses of DNA sequences found that C. remotus and C. pedroruizi are not reciprocally monophyletic and exhibit low genetic divergence ( 1 %). Therefore, C. remotus should be considered a junior synonym of C. pedroruizi. This work stresses the importance of comprehensive taxon sampling in poorly explored areas, especially between neighbouring countries.

Signatures of natural selection in morphological quantitative traits in Argentinean populations of Senegalia gilliesii (Fabaceae).

In order to elucidate the role of evolutionary forces in shaping the variation of quantitative traits in Senegalia gilliesii we evaluate seven phenotypic traits in three Argentinean populations, two of them sharing environmental and vegetation type conditions, and a third one ecologically differentiated from the former. The phenotypic traits were compared with molecular markers. Here, we search for signatures of selection by means of the comparison PST-FST . We assessed if the averages of the seven phenotypic traits were different among populations by means of ANOVA and we performed discriminant analysis of principal components (DAPC) for both morphological and molecular data. The ANOVA showed significant results only for two traits. For all foliar traits and two spine traits, the PST-FST comparison suggested the occurrence of stabilizing selection. The DAPC obtained from AFLP data showed three well defined groups of populations; when the same analysis was conducted with morphological data the scatterplot showed high overlapping among individuals and could not separate the populations. Overall, our findings suggest a prominent role of stabilizing selection in all foliar traits and stipular spine length. These results could be extrapolated to other tropical and subtropical acacias. Further studies are needed to analyse the mechanisms underlying genetic differentiation in natural populations of S. gilliesii, find its relationship with eco-geographical variables.

Too simple models may predict the island rule for the wrong reasons.

Biddick & Burns (2021) proposed a null/neutral model that reproduces the island rule as a product of random drift. We agree that it is unnecessary to assume adaptive processes driving island dwarfing or gigantism, but several flaws make their approach unrealistic and thus unsuitable as a stochastic model for evolutionary size changes.

Neutral evolution test of the spike protein of SARS-CoV-2 and its implications in the binding to ACE2.

As the SARS-CoV-2 has spread and the pandemic has dragged on, the virus continued to evolve rapidly resulting in the emergence of new highly transmissible variants that can be of public health concern. The evolutionary mechanisms that drove this rapid diversity are not well understood but neutral evolution should open the first insight. The neutral theory of evolution states that most mutations in the nucleic acid sequences are random and they can be fixed or disappear by purifying selection. Herein, we performed a neutrality test to better understand the selective pressures exerted over SARS-CoV-2 spike protein from homologue proteins of Betacoronavirus, as well as to the spikes from human clinical isolates of the virus. Specifically, Tyr and Asn have higher occurrence rates on the Receptor Binding Domain (RBD) and in the overall sequence of spike proteins of Betacoronavirus, whereas His and Arg have lower occurrence rates. The in vivo evolutionary phenomenon of SARS-CoV-2 shows that Glu, Lys, Phe, and Val have the highest probability of occurrence in the emergent viral particles. Amino acids that have higher occurrence than the expected by the neutral control, are favorable and are fixed in the sequence while the ones that have lower occurrence than expected, influence the stability and/or functionality of the protein. Our results show that most unique mutations either for SARS-CoV-2 or its variants of health concern are under selective pressures, which could be related either to the evasion of the immune system, increasing the virus' fitness or altering protein - protein interactions with host proteins. We explored the consequences of those selected mutations in the structure and protein - protein interaction with the receptor. Altogether all these forces have shaped the spike protein and the continually evolving variants.

Green, yellow or black? Genetic differentiation and adaptation signatures in a highly migratory marine turtle.

Marine species may exhibit genetic structure accompanied by phenotypic differentiation related to adaptation despite their high mobility. Two shape-based morphotypes have been identified for the green turtle (Chelonia mydas) in the Pacific Ocean: the south-central/western or yellow turtle and north-central/eastern or black turtle. The genetic differentiation between these morphotypes and the adaptation of the black turtle to environmentally contrasting conditions of the eastern Pacific region has remained a mystery for decades. Here we addressed both questions using a reduced-representation genome approach (Dartseq; 9473 neutral SNPs) and identifying candidate outlier loci (67 outlier SNPs) of biological relevance between shape-based morphotypes from eight Pacific foraging grounds (n = 158). Our results support genetic divergence between morphotypes, probably arising from strong natal homing behaviour. Genes and enriched biological functions linked to thermoregulation, hypoxia, melanism, morphogenesis, osmoregulation, diet and reproduction were found to be outliers for differentiation, providing evidence for adaptation of C. mydas to the eastern Pacific region and suggesting independent evolutionary trajectories of the shape-based morphotypes. Our findings support the evolutionary distinctness of the enigmatic black turtle and contribute to the adaptive research and conservation genomics of a long-lived and highly mobile vertebrate.

Evidence of local adaptation despite strong drift in a Neotropical patchily distributed bromeliad.

Both genetic drift and divergent selection are predicted to be drivers of population differentiation across patchy habitats, but the extent to which these forces act on natural populations to shape traits is strongly affected by species' ecological features. In this study, we infer the genomic structure of Pitcairnia lanuginosa, a widespread herbaceous perennial plant with a patchy distribution. We sampled populations in the Brazilian Cerrado and the Central Andean Yungas and discovered and genotyped SNP markers using double-digest restriction-site associated DNA sequencing. In addition, we analyzed ecophysiological traits obtained from a common garden experiment and compared patterns of phenotypic and genetic divergence (PST-FST comparisons) in a subset of populations from the Cerrado. Our results from molecular analyses pointed to extremely low genetic diversity and a remarkable population differentiation, supporting a major role of genetic drift. Approximately 0.3% of genotyped SNPs were flagged as differentiation outliers by at least two distinct methods, and Bayesian generalized linear mixed models revealed a signature of isolation by environment in addition to isolation by distance for high-differentiation outlier SNPs among the Cerrado populations. PST-FST comparisons suggested divergent selection on two ecophysiological traits linked to drought tolerance. We showed that these traits vary among populations, although without any particular macro-spatial pattern, suggesting local adaptation to differences in micro-habitats. Our study shows that selection might be a relevant force, particularly for traits involved in drought stress, even for populations experiencing strong drift, which improves our knowledge on eco-evolutionary processes acting on non-continuously distributed species.

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.

Idiosyncratic responses to drivers of genetic differentiation in the complex landscapes of Isthmian Central America.

Isthmian Central America (ICA) is one of the most biodiverse regions in the world, hosting an exceptionally high number of species per unit area. ICA was formed <25 million years ago and, consequently, its biotic assemblage is relatively young and derived from both colonization and in situ diversification. Despite intensive taxonomic work on the local fauna, the potential forces driving genetic divergences and ultimately speciation in ICA remain poorly studied. Here, we used a landscape genetics approach to test whether isolation by distance, topography, habitat suitability, or environment drive the genetic diversity of the regional frog assemblage. To this end, we combined data on landscape features and mitochondrial DNA sequence variation for nine codistributed amphibian species with disparate life histories. In five species, we found that at least one of the factors tested explained patterns of genetic divergence. However, rather than finding a general pattern, our results revealed idiosyncratic responses to historical and ecological processes, indicating that intrinsic life-history characteristics may determine the effect of different drivers of isolation on genetic divergence in ICA. Our work also suggests that the convergence of several factors promoting isolation among populations over a heterogeneous landscape might maximize genetic differentiation, despite short geographical distances. In conclusion, abiotic factors and geographical features have differentially affected the genetic diversity across the regional frog assemblage. Much more complex models (i.e., considering multiple drivers), beyond simple vicariance of Caribbean and Pacific lineages, are needed to better understand the evolutionary history of ICA's diverse biotas.