Exploring seascape genetics and kinship in the reef sponge Stylissa carteri in the Red Sea

A main goal of population geneticists is to study patterns of gene flow to gain a better understanding of the population structure in a given organism. To date most efforts have been focused on studying gene flow at either broad scales to identify barriers to gene flow and isolation by distance or at fine spatial scales in order to gain inferences regarding reproduction and local dispersal. Few studies have measured connectivity at multiple spatial scales and have utilized novel tools to test the influence of both environment and geography on shaping gene flow in an organism. Here a seascape genetics approach was used to gain insight regarding geographic and ecological barriers to gene flow of a common reef sponge, Stylissa carteri in the Red Sea. Furthermore, a small‐scale (<1 km) analysis was also conducted to infer reproductive potential in this organism. At the broad scale, we found that sponge connectivity is not structured by geography alone, but rather, genetic isolation in the southern Red Sea correlates strongly with environmental heterogeneity. At the scale of a 50‐m transect, spatial autocorrelation analyses and estimates of full‐siblings revealed that there is no deviation from random mating. However, at slightly larger scales (100–200 m) encompassing multiple transects at a given site, a greater proportion of full‐siblings was found within sites versus among sites in a given location suggesting that mating and/or dispersal are constrained to some extent at this spatial scale. This study adds to the growing body of literature suggesting that environmental and ecological variables play a major role in the genetic structure of marine invertebrate populations.     

Living on the edge: the role of geography and environment in structuring genetic variation in the southernmost populations of a tropical oak

Understanding the factors determining genetic diversity and structure in peripheral populations is a long‐standing goal of evolutionary biogeography, yet little empirical information is available for tropical species. In this study, we combine information from nuclear microsatellite markers and niche modelling to analyse the factors structuring genetic variation across the southernmost populations of the tropical oak Quercus segoviensis. First, we tested the hypothesis that genetic variability decreases with population isolation and increases with local habitat suitability and stability since the Last Glacial Maximum (LGM). Second, we employed a recently developed multiple matrix regression with randomisation (MMRR) approach to study the factors associated with genetic divergence among the studied populations and test the relative contribution of environmental and geographic isolation to contemporary patterns of genetic differentiation. We found that genetic diversity was negatively correlated with average genetic differentiation with other populations, indicating that isolation and limited gene flow have contributed to erode genetic variability in some populations. Considering the relatively small size of the study area (<120 km), analyses of genetic structure indicate a remarkable inter‐population genetic differentiation. Environmental dissimilarity and differences in current and past climate niche suitability and their additive effects were not associated with genetic differentiation after controlling for geographic distance, indicating that local climate does not contribute to explain spatial patterns of genetic structure. Overall, our data indicate that geographic isolation, but not current or past climate, is the main factor determining contemporary patterns of genetic diversity and structure within the southernmost peripheral populations of this tropical oak.     

Drivers of genetic differentiation in a generalist insect‐pollinated herb across spatial scales

The isolation‐by‐distance model (IBD) predicts that genetic differentiation among populations increases with geographic distance. Yet, empirical studies show that a variety of ecological, topographic and historical factors may override the effect of geographic distance on genetic variation. This may particularly apply to species with narrow but highly heterogeneous distribution ranges, such as those occurring along elevational gradients. Using nine SSR markers, we study the genetic differentiation of the montane pollination‐generalist herb, Erysimum mediohispanicum. Because the effects of any given factor may depend on the geographic scale considered, we investigate the contribution of different environmental and historical factors at three different spatial scales. We evaluate five competing models that put forward the role of geographic distance, local environmental factors [biotic interactions (IBEb) and climatic variables (IBEa)], landscape resistance (IBR) and phylogeographic patterns (IBP), respectively. We find significant IBD regardless of the spatial scale and the genetic distance estimator considered. However, IBEa and IBP also play a prominent role in shaping genetic differentiation patterns at the larger spatial scales, and IBR is significant at the fine spatial scale. Overall, our results highlight the importance of combining different estimators, statistical approaches and spatial scales to disentangle the relative importance of the various ecological factors contributing to the shaping of genetic divergence patterns in natural populations.     

Low but contrasting neutral genetic differentiation shaped by winter temperature in European great tits

Gene flow is usually thought to reduce genetic divergence and impede local adaptation by homogenising gene pools between populations. However, evidence for local adaptation and phenotypic differentiation in highly mobile species, experiencing high levels of gene flow, is emerging. Assessing population genetic structure at different spatial scales is thus a crucial step towards understanding mechanisms underlying intraspecific differentiation and diversification. Here, we studied the population genetic structure of a highly mobile species – the great tit Parus major – at different spatial scales. We analysed 884 individuals from 30 sites across Europe including 10 close‐by sites (< 50 km), using 22 microsatellite markers. Overall we found a low but significant genetic differentiation among sites (F_ST = 0.008). Genetic differentiation was higher, and genetic diversity lower, in south‐western Europe. These regional differences were statistically best explained by winter temperature. Overall, our results suggest that great tits form a single patchy metapopulation across Europe, in which genetic differentiation is independent of geographical distance and gene flow may be regulated by environmental factors via movements related to winter severity. This might have important implications for the evolutionary trajectories of sub‐populations, especially in the context of climate change, and calls for future investigations of local differences in costs and benefits of philopatry at large scales.     

Environmental gradients predict the genetic population structure of a coral reef fish in the Red Sea

The relatively recent fields of terrestrial landscape and marine seascape genetics seek to identify the influence of biophysical habitat features on the spatial genetic structure of populations or individuals. Over the last few years, there has been accumulating evidence for the effect of environmental heterogeneity on patterns of gene flow and connectivity in marine systems. Here, we investigate the population genetic patterns of an anemonefish, Amphiprion bicinctus, along the Saudi Arabian coast of the Red Sea. We collected nearly one thousand samples from 19 locations, spanning approximately 1500 km, and genotyped them at 38 microsatellite loci. Patterns of gene flow appeared to follow a stepping‐stone model along the northern and central Red Sea, which was disrupted by a distinct genetic break at a latitude of approximately 19°N. The Red Sea is characterized by pronounced environmental gradients along its axis, roughly separating the northern and central from the southern basin. Using mean chlorophyll‐a concentrations as a proxy for this gradient, we ran tests of isolation by distance (IBD, R² = 0.52) and isolation by environment (IBE, R² = 0.64), as well as combined models using partial Mantel tests and multiple matrix regression with randomization (MMRR). We found that genetic structure across our sampling sites may be best explained by a combined model of IBD and IBE (Mantel: R² = 0.71, MMRR: R² = 0.86). Our results highlight the potential key role of environmental patchiness in shaping patterns of gene flow in species with pelagic larval dispersal. We support growing calls for the integration of biophysical habitat characteristics into future studies of population genetic structure.     

Are habitat fragmentation,local adaptation and isolation‐by‐distance driving population divergence in wild rice Oryza rufipogon?

Habitat fragmentation weakens the connection between populations and is accompanied with isolation by distance (IBD) and local adaptation (isolation by adaptation, IBA), both leading to genetic divergence between populations. To understand the evolutionary potential of a population and to formulate proper conservation strategies, information on the roles of IBD and IBA in driving population divergence is critical. The putative ancestor of Asian cultivated rice (Oryza sativa) is endangered in China due to habitat loss and fragmentation. We investigated the genetic variation in 11 Chinese Oryza rufipogon populations using 79 microsatellite loci to infer the effects of habitat fragmentation, IBD and IBA on genetic structure. Historical and current gene flows were found to be rare (m_h = 0.0002–0.0013, m_c = 0.007–0.029), indicating IBD and resulting in a high level of population divergence (F_ST = 0.343). High within‐population genetic variation (H_E = 0.377–0.515), relatively large effective population sizes (N_e = 96–158), absence of bottlenecks and limited gene flow were found, demonstrating little impact of recent habitat fragmentation on these populations. Eleven gene‐linked microsatellite loci were identified as outliers, indicating local adaptation. Hierarchical AMOVA and partial Mantel tests indicated that population divergence of Chinese O. rufipogon was significantly correlated with environmental factors, especially habitat temperature. Common garden trials detected a significant adaptive population divergence associated with latitude. Collectively, these findings imply that IBD due to historical rather than recent fragmentation, followed by local adaptation, has driven population divergence in O. rufipogon.     

Local and regional founder effects in lake zooplankton persist after thousands of years despite high dispersal potential

We reconstructed the genetic structure of a planktonic crustacean Daphnia longispina living in high mountain lakes and ponds in the Pyrenees to investigate whether it was shaped by persistent founder effects originating shortly after the last glacial maximum or by ongoing dispersal and effective migration (gene flow). We found that the genetic structure can largely be explained by a single colonization event following gradual deglaciation of the Pyrenees ~10 000–15 000 years ago. Nuclear genetic diversity declined steeply from southeast to northwest, suggestive of serial colonization of available habitats with advancing deglaciation. The spatial genetic structure suggests that founder effects were major determinants of the present‐day diversity, both at the catchment level and at the level of individual water bodies, further supporting extremely low effective migration rates. This study reveals a prime example of a founder effect that is both long lasting and maintained at small spatial scales. Our data suggest a process of isolation by colonization as a result of strong priority effects and monopolization. We found evidence for the spread of haplotypes with Pyrenean ancestry across the Palaearctic over distances up to 5500 km, although the local genetic structure after colonization was hardly influenced by contemporary dispersal. Finally, our data also suggest that mitochondrial mutation rates in the studied populations were seven times higher than typically assumed. Overall, we show that founder effects can persist for centuries even at small spatial scales at which the potential for dispersal is high.     

Genome‐wide differentiation in closely related populations: the roles of selection and geographic isolation

Population divergence in geographic isolation is due to a combination of factors. Natural and sexual selection may be important in shaping patterns of population differentiation, a pattern referred to as ‘isolation by adaptation’ (IBA). IBA can be complementary to the well‐known pattern of ‘isolation by distance’ (IBD), in which the divergence of closely related populations (via any evolutionary process) is associated with geographic isolation. The barn swallow Hirundo rustica complex comprises six closely related subspecies, where divergent sexual selection is associated with phenotypic differentiation among allopatric populations. To investigate the relative contributions of selection and geographic distance to genome‐wide differentiation, we compared genotypic and phenotypic variation from 350 barn swallows sampled across eight populations (28 pairwise comparisons) from four different subspecies. We report a draft whole‐genome sequence for H. rustica, to which we aligned a set of 9493 single nucleotide polymorphisms (SNPs). Using statistical approaches to control for spatial autocorrelation of phenotypic variables and geographic distance, we find that divergence in traits related to migratory behaviour and sexual signalling, as well as geographic distance, together explain over 70% of genome‐wide divergence among populations. Controlling for IBD, we find 42% of genomewide divergence is attributable to IBA through pairwise differences in traits related to migratory behaviour and sexual signalling alone. By (i) combining these results with prior studies of how selection shapes morphological differentiation and (ii) accounting for spatial autocorrelation, we infer that morphological adaptation plays a large role in shaping population‐level differentiation in this group of closely related populations.     

Spatial scale affects landscape genetic analysis of a wetland grasshopper

Most landscape genetic studies assess the impact of landscape elements on species' dispersal and gene flow. Many of these studies perform their analysis on all possible population pairs in a study area and do not explicitly consider the effects of spatial scale and population network topology on their results. Here, we examined the effects of spatial scale and population network topology on the outcome of a landscape genetic analysis. Additionally, we tested whether the relevant spatial scale of landscape genetic analysis could be defined by population network topology or by isolation‐by‐distance (IBD) patterns. A data set of the wetland grasshopper Stethophyma grossum, collected in a fragmented agricultural landscape, was used to analyse population network topology, IBD patterns and dispersal habitats, using least‐cost transect analysis. Landscape genetic analyses neglecting spatial scale and population network topology resulted in models with low fits, with which a most likely dispersal habitat could not be identified. In contrast, analyses considering spatial scale and population network topology resulted in high model fits by restricting landscape genetic analysis to smaller scales (0–3 km) and neighbouring populations, as represented by a Gabriel graph. These models also successfully identified a likely dispersal habitat of S. grossum. The above results suggest that spatial scale and potentially population network topology should be more explicitly considered in future landscape genetic analyses.     

Closely related octopus species show different spatial genetic structures in response to the Antarctic seascape

Determining whether comparable processes drive genetic divergence among marine species is relevant to molecular ecologists and managers alike. Sympatric species with similar life histories might be expected to show comparable patterns of genetic differentiation and a consistent influence of environmental factors in shaping divergence. We used microsatellite loci to quantify genetic differentiation across the Scotia Arc in three species of closely related benthic octopods, Pareledone turqueti, P. charcoti, and Adelieledone polymorpha. The relative importance of environmental factors (latitude, longitude, depth, and temperature) in shaping genetic structure was investigated when significant spatial genetic structure was uncovered. Isolated populations of P. turqueti and A. polymorpha at these species’ range margins were genetically different to samples close to mainland Antarctica; however, these species showed different genetic structures at a regional scale. Samples of P. turqueti from the Antarctic Peninsula, Elephant Island, and Signy Island were genetically different, and this divergence was associated primarily with sample collection depth. By contrast, weak or nonsignificant spatial genetic structure was evident across the Antarctic Peninsula, Elephant Island, and Signy Island region for A. polymorpha, and slight associations between population divergence and temperature or depth (and/or longitude) were detected. Pareledone charcoti has a limited geographic range, but exhibited no genetic differentiation between samples from a small region of the Scotia Arc (Elephant Island and the Antarctic Peninsula). Thus, closely related species with similar life history strategies can display contrasting patterns of genetic differentiation depending on spatial scale; moreover, depth may drive genetic divergence in Southern Ocean benthos.     

Living in two worlds: Evolutionary mechanisms act differently in the native and introduced ranges of an invasive plant

Identifying the factors that influence spatial genetic structure among populations can provide insights into the evolution of invasive plants. In this study, we used the common reed (Phragmites australis), a grass native in Europe and invading North America, to examine the relative importance of geographic, environmental (represented by climate here), and human effects on population genetic structure and its changes during invasion. We collected samples of P. australis from both the invaded North American and native European ranges and used molecular markers to investigate the population genetic structure within and between ranges. We used path analysis to identify the contributions of each of the three factors—geographic, environmental, and human‐related—to the formation of spatial genetic patterns. Genetic differentiation was observed between the introduced and native populations, and their genetic structure in the native and introduced ranges was different. There were strong effects of geography and environment on the genetic structure of populations in the native range, but the human‐related factors manifested through colonization of anthropogenic habitats in the introduced range counteracted the effects of environment. The between‐range genetic differences among populations were mainly explained by the heterogeneous environment between the ranges, with the coefficient 2.6 times higher for the environment than that explained by the geographic distance. Human activities were the primary contributor to the genetic structure of the introduced populations. The significant environmental divergence between ranges and the strong contribution of human activities to the genetic structure in the introduced range suggest that invasive populations of P. australis have evolved to adapt to a different climate and to human‐made habitats in North America.     

Testing the role of ancient and contemporary landscapes on structuring genetic variation in a specialist grasshopper

Understanding the processes underlying spatial patterns of genetic diversity and structure of natural populations is a central topic in evolutionary biogeography. In this study, we combine data on ancient and contemporary landscape composition to get a comprehensive view of the factors shaping genetic variation across the populations of the scrub‐legume grasshopper (Chorthippus binotatus binotatus) from the biogeographically complex region of southeast Iberia. First, we examined geographical patterns of genetic structure and employed an approximate Bayesian computation (ABC) approach to compare different plausible scenarios of population divergence. Second, we used a landscape genetic framework to test for the effects of (1) Late Miocene paleogeography, (2) Pleistocene climate fluctuations, and (3) contemporary topographic complexity on the spatial patterns of population genetic differentiation. Genetic structure and ABC analyses supported the presence of three genetic clusters and a sequential west‐to‐east splitting model that predated the last glacial maximum (LGM, c. 21 Kya). Landscape genetic analyses revealed that population genetic differentiation was primarily shaped by contemporary topographic complexity, but was not explained by any paleogeographic scenario or resistance distances based on climate suitability in the present or during the LGM. Overall, this study emphasizes the need of integrating information on ancient and contemporary landscape composition to get a comprehensive view of their relative importance to explain spatial patterns of genetic variation in organisms inhabiting regions with complex biogeographical histories.     

Spatial and ecological population genetic structures within two island‐endemic Aeonium species of different niche width

The Crassulacean genus Aeonium is a well‐known example for plant species radiation on oceanic archipelagos. However, while allopatric speciation among islands is documented for this genus, the role of intra‐island speciation due to population divergence by topographical isolation or ecological heterogeneity has not yet been addressed. The aim of this study was to investigate intraspecific genetic structures and to identify spatial and ecological drivers of genetic population differentiation on the island scale. We analyzed inter simple sequence repeat variation within two island‐endemic Aeonium species of La Palma: one widespread generalist that covers a large variety of different habitat types (Ae. davidbramwellii) and one narrow ecological specialist (Ae. nobile), in order to assess evolutionary potentials on this island. Gene pool differentiation and genetic diversity patterns were associated with major landscape structures in both species, with phylogeographic implications. However, overall levels of genetic differentiation were low. For the generalist species, outlier loci detection and loci–environment correlation approaches indicated moderate signatures of divergent selection pressures linked to temperature and precipitation variables, while the specialist species missed such patterns. Our data point to incipient differentiation among populations, emphasizing that ecological heterogeneity and topographical structuring within the small scales of an island can foster evolutionary processes. Very likely, such processes have contributed to the radiation of Aeonium on the Canary Islands. There is also support for different evolutionary mechanisms between generalist and specialist species.     

History,geography and host use shape genomewide patterns of genetic variation in the redheaded pine sawfly (Neodiprion lecontei)

Divergent host use has long been suspected to drive population differentiation and speciation in plant‐feeding insects. Evaluating the contribution of divergent host use to genetic differentiation can be difficult, however, as dispersal limitation and population structure may also influence patterns of genetic variation. In this study, we use double‐digest restriction‐associated DNA (ddRAD) sequencing to test the hypothesis that divergent host use contributes to genetic differentiation among populations of the redheaded pine sawfly (Neodiprion lecontei), a widespread pest that uses multiple Pinus hosts throughout its range in eastern North America. Because this species has a broad range and specializes on host plants known to have migrated extensively during the Pleistocene, we first assess overall genetic structure using model‐based and model‐free clustering methods and identify three geographically distinct genetic clusters. Next, using a composite‐likelihood approach based on the site frequency spectrum and a novel strategy for maximizing the utility of linked RAD markers, we infer the population topology and date divergence to the Pleistocene. Based on existing knowledge of Pinus refugia, estimated demographic parameters and patterns of diversity among sawfly populations, we propose a Pleistocene divergence scenario for N. lecontei. Finally, using Mantel and partial Mantel tests, we identify a significant relationship between genetic distance and geography in all clusters, and between genetic distance and host use in two of three clusters. Overall, our results indicate that Pleistocene isolation, dispersal limitation and ecological divergence all contribute to genomewide differentiation in this species and support the hypothesis that host use is a common driver of population divergence in host‐specialized insects.     

Widespread evidence for incipient ecological speciation: a meta‐analysis of isolation‐by‐ecology

Ecologically mediated selection has increasingly become recognised as an important driver of speciation. The correlation between neutral genetic differentiation and environmental or phenotypic divergence among populations, to which we collectively refer to as isolation‐by‐ecology (IBE), is an indicator of ecological speciation. In a meta‐analysis framework, we determined the strength and commonality of IBE in nature. On the basis of 106 studies, we calculated a mean effect size of IBE with and without controlling for spatial autocorrelation among populations. Effect sizes were 0.34 (95% CI 0.24–0.42) and 0.26 (95% CI 0.13–0.37), respectively, indicating that an average of 5% of the neutral genetic differentiation among populations was explained purely by ecological contrast. Importantly, spatial autocorrelation reduced IBE correlations for environmental variables, but not for phenotypes. Through simulation, we showed how the influence of isolation‐by‐distance and spatial autocorrelation of ecological variables can result in false positives or underestimated correlations if not accounted for in the IBE model. Collectively, this meta‐analysis showed that ecologically induced genetic divergence is pervasive across time‐scales and taxa, and largely independent of the choice of molecular marker. We discuss the importance of these results in the context of adaptation and ecological speciation and suggest future research avenues.     

Fine‐scale spatial genetic structure of common and declining bumble bees across an agricultural landscape

Land‐use changes have threatened populations of many insect pollinators, including bumble bees. Patterns of dispersal and gene flow are key determinants of species' ability to respond to land‐use change, but have been little investigated at a fine scale (<10 km) in bumble bees. Using microsatellite markers, we determined the fine‐scale spatial genetic structure of populations of four common Bombus species (B. terrestris, B. lapidarius, B. pascuorum and B. hortorum) and one declining species (B. ruderatus) in an agricultural landscape in Southern England, UK. The study landscape contained sown flower patches representing agri‐environment options for pollinators. We found that, as expected, the B. ruderatus population was characterized by relatively low heterozygosity, number of alleles and colony density. Across all species, inbreeding was absent or present but weak (F_IS = 0.01–0.02). Using queen genotypes reconstructed from worker sibships and colony locations estimated from the positions of workers within these sibships, we found that significant isolation by distance was absent in B. lapidarius, B. hortorum and B. ruderatus. In B. terrestris and B. pascuorum, it was present but weak; for example, in these two species, expected relatedness of queens founding colonies 1 m apart was 0.02. These results show that bumble bee populations exhibit low levels of spatial genetic structure at fine spatial scales, most likely because of ongoing gene flow via widespread queen dispersal. In addition, the results demonstrate the potential for agri‐environment scheme conservation measures to facilitate fine‐scale gene flow by creating a more even distribution of suitable habitats across landscapes.     

Partitioning the effects of isolation by distance,environment, and physical barriers on genomic divergence between parapatric threespine stickleback

Genetic divergence between populations is shaped by a combination of drift, migration, and selection, yielding patterns of isolation‐by‐distance (IBD) and isolation‐by‐environment (IBE). Unfortunately, IBD and IBE may be confounded when comparing divergence across habitat boundaries. For instance, parapatric lake and stream threespine stickleback (Gasterosteus aculeatus) may have diverged due to selection against migrants (IBE), or mere spatial separation (IBD). To quantitatively partition the strength of IBE and IBD, we used recently developed population genetic software (BEDASSLE) to analyze partial genomic data from three lake‐stream clines on Vancouver Island. We find support for IBD within each of three outlet streams (unlike prior studies of lake‐stream stickleback). In addition, we find evidence for IBE (controlling for geographic distance): the genetic effect of habitat is equivalent to geographic separation of ～1.9 km of IBD. Remarkably, of our three lake‐stream pairs, IBE is strongest where migration between habitats is easiest. Such microgeographic genetic divergence would require exceptionally strong divergent selection, which multiple experiments have failed to detect. Instead, we propose that nonrandom dispersal (e.g., habitat choice) contributes to IBE. Supporting this conclusion, we show that the few migrants between habitats are a nonrandom subset of the phenotype distribution of the source population.     

Niche specificity influences gene flow across fine‐scale habitat mosaics in Succulent Karoo plants

While the tempo of diversification in biodiversity hotspots has received much attention, the spatial scale of diversification has often been overlooked. Addressing this deficiency requires understanding the drivers of population divergence and the spatial scales at which they operate in species‐rich clades and ecosystems. South Africa's Succulent Karoo (SK) hotspot provides an excellent system for such research, being both compact (ca. 110,000 km²) and home to spectacular in‐situ radiations, such as the ruschioid Aizoaceae. Here we use GBS to document genetic structure in two co‐occurring ruschioid species, at both coarse (>10 km) and fine (<500 m) spatial scales. Where Ruschia burtoniae shows strong between‐population genetic differentiation and no gene flow, Conophytum calculus shows weak differentiation, with high levels of admixture suggesting recent or ongoing gene flow. Community analysis and transplant experiments reveal that R. burtoniae occupies a narrow, low‐pH edaphic niche, and at scales of a few hundred metres, areas of elevated genetic turnover correspond to patches of edaphically unsuitable habitat. In contrast, C. calculus occupies a broader niche and exhibits isolation‐by‐distance without a habitat effect. We suggest that edaphic specialisation, coupled with highly restricted seed and pollen dispersal in heterogeneous landscapes, has played a major role in driving rapid diversification at small spatial scales in this system. However, the contrasting patterns in our study species show that these factors do not influence all organisms uniformly, being strongly modulated by lineage‐specific traits that influence both the spatial scale of gene flow and habitat specificity.     

Drivers for genetic structure at different geographic scales for Pacific red snapper (Lutjanus peru) and yellow snapper (Lutjanus argentiventris) in the tropical eastern Pacific

The tropical eastern Pacific (TEP) is a highly dynamic region and a model system to study how habitat discontinuities affect the distribution of shorefishes, particularly for species that display ontogenetic habitat shifts, including snappers (Lutjanidae). To evaluate the genetic structure of the Pacific red snapper (Lutjanus peru) and the yellow snapper (Lutjanus argentiventris) throughout their distribution range along the TEP, 13 and 11 microsatellite loci were analysed, respectively. The genetic diversity of L. peru (N = 446) and L. argentiventris (N = 170) was evaluated in 10 and 5 localities, respectively, showing slightly higher but non-significant values in the Gulf of California for both species. The genetic structure analysis identified the presence of significant genetic structure in both species, but the locations of the identified barriers for the gene flow differed between species. The principal driver for the genetic structure at large scales >2500 km was isolation by distance. At smaller scales (<250 km), the habitat discontinuity for juveniles and adults and the environmental differences throughout the distribution range represented potential barriers to gene flow between populations for both species.