Phylogeography and population structure of the Yunnan snub-nosed monkey (Rhinopithecus bieti) inferred from mitochondrial control region DNA sequence analysis

Rhinopithecus bieti, the Yunnan snub-nosed monkey, is the nonhuman primate with the highest altitudinal distribution and is also one of the 25 most globally endangered primate species. Currently, R. bieti is found in forests between 3000 and 4500 m above sea level, within a narrow area on the Tibetan Plateau between the Yangtze and Mekong rivers, where it is suffering from loss of habitat and shrinking population size (approximately 1500). To assess the genetic diversity within this species, its population structure and to infer its evolutionary history, we sequenced 401 bp of the hypervariable I (HVI) segment from the mitochondrial DNA control region (CR) for 157 individuals from 11 remnant patches throughout the fragmented distribution area. Fifty-two variable sites were observed and 30 haplotypes were defined. Compared with other primate species, R. bieti cannot be regarded as a taxon with low genetic diversity. Phylogenetic analysis partitioned haplotypes into two divergent haplogroups (A and B). Haplotypes from the two mitochondrial clades were found to be mixed in some patches although the distribution of haplotypes displayed local homogeneity, implying a strong population structure within R. bieti. Analysis of molecular variance detected significant differences among the different geographical regions, suggesting that R. bieti should be separated into three management units (MUs) for conservation. Based on our results, it can be hypothesized that the genetic history of R. bieti includes an initial, presumably allopatric divergence between clades A and B 1.0-0.7 million years ago (Ma), which might have been caused by the Late Cenozoic uplift of the Tibetan Plateau, secondary contact after this divergence as a result of a population expansion 0.16-0.05 Ma, and population reduction and habitat fragmentation in the very recent past.     

Relationship between microspatial population genetic structure and habitat heterogeneity in Pomatias elegans (O.F. Müller 1774) (Caenogastropoda, Pomatiasidae)

In the present study the population genetic structure of the terrestrial snail Pomatias elegans was related to habitat structure on a microspatial scale. The genetic variability of 1607 individuals from 51 sampling sites in five different populations in Provence, France, was studied with an allozyme marker using population genetic methods, Mantel tests and spatial autocorrelation techniques were applied to different connectivity networks accounting for the structural features of the landscape. It is suggested that the population structure is, to a large extent, a function of the habitat quality, quantified as population density, and of the spatial arrangement of the habitat in the landscape and not of the geographical distance per se . In fragmented habitats, random genetic drift was the prevailing force for sampling sites separated by a few hundred meters.  © 2002 The Linnean Society of London, Biological Journal of the Linnean Society , 2002, 76 , 565–575.     

Effect of anthropogenic landscape features on population genetic differentiation of Przewalski's gazelle: main role of human settlement

Anthropogenic landscapes influence evolutionary processes such as population genetic differentiation, however, not every type of landscape features exert the same effect on a species, hence it is necessary to estimate their relative effect for species management and conservation. Przewalski's gazelle (Procapra przewalskii), which inhabits a human-altered area on Qinghai-Tibet Plateau, is one of the most endangered antelope species in the world. Here, we report a landscape genetic study on Przewalski's gazelle. We used skin and fecal samples of 169 wild gazelles collected from nine populations and thirteen microsatellite markers to assess the genetic effect of anthropogenic landscape features on this species. For comparison, the genetic effect of geographical distance and topography were also evaluated. We found significant genetic differentiation, six genetic groups and restricted dispersal pattern in Przewalski's gazelle. Topography, human settlement and road appear to be responsible for observed genetic differentiation as they were significantly correlated with both genetic distance measures [F(ST)/(1-F(ST)) and F'(ST)/(1-F'(ST))] in Mantel tests. IBD (isolation by distance) was also inferred as a significant factor in Mantel tests when genetic distance was measured as F(ST)/(1-F(ST)). However, using partial Mantel tests, AIC(c) calculations, causal modeling and AMOVA analysis, we found that human settlement was the main factor shaping current genetic differentiation among those tested. Altogether, our results reveal the relative influence of geographical distance, topography and three anthropogenic landscape-type on population genetic differentiation of Przewalski's gazelle and provide useful information for conservation measures on this endangered species.     

Environmental factors affecting population level genetic divergence of the striped field mouse (<Emphasis Type="Italic">Apodemus agrarius</Emphasis>) in South Korea

The cognizing of connectivity among small mammal populations across heterogeneous landscapes is complicated due to complex influences of landscape and anthropogenic factors on gene flow. A landscape genetics approach offers inferences on how landscape features drive population structure. Through a landscape genetics approach, we investigated influences of geographical, environmental, and anthropogenic features on populations of Apodemus agrarius, the striped field mouse, the prime vector of hemorrhagic fever by a landscape genetic approach. We identified landscape features that might affect the population structure of striped field mice by analyzing microsatellite markers of 197 striped field mice from 21 populations throughout South Korea. We developed Maximum-likelihood population effects models based on landscape distances and resistance matrices and pairwise F_ST values for meta-populations of striped field mouse. We also conducted Mantel and partial Mantel tests to investigate geographic patterns of genetic similarities. In Mantel and partial Mantel tests, the F_ST was significantly correlated with all three models of movement; movement cost, Euclidian distance and least-cost distance, although the magnitudes of correlations varied. The 4 top-ranked models included three variables; temperature, precipitation and one human disturbance factor (population). We did not attain a significant effect for anthropogenic factors on genetic similarities among populations in the Korean striped field mouse, but we confirmed a significant association for genetic similarity with climatic features (temperature and precipitation).     

Using multiple landscape genetic approaches to test the validity of genetic clusters in a species characterized by an isolation‐by‐distance pattern

Bayesian clustering methods are typically used to identify barriers to gene flow, but they are prone to deduce artificial subdivisions in a study population characterized by an isolation‐by‐distance pattern (IbD). Here we analysed the landscape genetic structure of a population of wild boars (Sus scrofa) from south‐western Germany. Two clustering methods inferred the presence of the same genetic discontinuity. However, the population in question was characterized by a strong IbD pattern. While landscape‐resistance modelling failed to identify landscape features that influenced wild boar movement, partial Mantel tests and multiple regression of distance matrices (MRDMs) suggested that the empirically inferred clusters were separated by a genuine barrier. When simulating random lines bisecting the study area, 60% of the unique barriers represented, according to partial Mantel tests and MRDMs, significant obstacles to gene flow. By contrast, the random‐lines simulation showed that the boundaries of the inferred empirical clusters corresponded to the most important genetic discontinuity in the study area. Given the degree of habitat fragmentation separating the two empirical partitions, it is likely that the clustering programs correctly identified a barrier to gene flow. The differing results between the work published here and other studies suggest that it will be very difficult to draw general conclusions about habitat permeability in wild boar from individual studies.     

Low levels of genetic differentiation with isolation by geography and environment in populations of Drosophila melanogaster from across China

The fruit fly, Drosophila melanogaster, is a model species in evolutionary studies. However, population processes of this species in East Asia are poorly studied. Here we examined the population genetic structure of D. melanogaster across China. There were 14 mitochondrial haplotypes with 10 unique ones out of 23 known from around the globe. Pairwise F_ST values estimated from 15 novel microsatellites ranged from 0 to 0.11, with geographically isolated populations showing the highest level of genetic uniqueness. STRUCTURE analysis identified high levels of admixture at both the individual and population levels. Mantel tests indicated a strong association between genetic distance and geographical distance as well as environmental distance. Full redundancy analysis (RDA) showed that independent effects of environmental conditions and geography accounted for 62.10% and 31.58% of the total explained genetic variance, respectively. When geographic variables were constrained in a partial RDA analysis, the environmental variables bio2 (mean diurnal air temperature range), bio13 (precipitation of the wettest month), and bio15 (precipitation seasonality) were correlated with genetic distance. Our study suggests that demographic history, geographical isolation, and environmental factors have together shaped the population genetic structure of D. melanogaster after its introduction into China.Subject terms: Genetic variation, Population genetics     

Population subdivision in marine environments: the contributions of biogeography, geographical distance and discontinuous habitat to genetic differentiation in a blennioid fish, Axoclinus nigricaudus

The relative importance of factors that may promote genetic differentiation in marine organisms is largely unknown. Here, contributions to population structure from a biogeographic boundary, geographical distance and the distribution of suitable habitat were investigated in Axoclinus nigricaudus, a small subtidal rock-reef fish, throughout its range in the Gulf of California. A 408-bp fragment of the mitochondrial control region was sequenced from 105 individuals. Variation was significantly partitioned between 28 of 36 possible combinations of population pairs. Phylogenetic analyses, hierarchical analyses of variance and a modified Mantel test substantiated a major break between two putative biogeographic regions. This genetic discontinuity coincides with an abrupt change in ecological characteristics, including temperature and salinity, but does not coincide with known oceanographic circulation patterns or any known historic barriers. There was an overall relationship of increasing genetic distance with increasing geographical distance between population pairs, in a manner consistent with isolation-by-distance. A significant habitat-by-geographical-distance interaction term indicated that, for a given geographical distance, populations separated by discontinuous habitat (sand) are more distinct genetically than are populations separated by continuous habitat (rock). In addition, populations separated by deep open waters were more genetically distinct than populations separated by continuous habitat (rock). These results indicate that levels of genetic differentiation among populations of A. nigricaudus cannot be explained by a single factor, but are due to the combined influences of biogeography, geographical distance and availability of suitable habitat.     

Genome-wide association to fine-scale ecological heterogeneity within a continuous population of Biscutella laevigata (Brassicaceae)

Gene flow, drift and selection can be detected through different signatures across the genome and the landscape. Genetic discontinuities along with their correlation to environmental features can be used to tease out isolation-by-distance and isolation-by-time from processes related to selection. Using spatial statistics (spatial autocorrelation methods, canonical correspondence analysis and partial Mantel tests) dealing with genome-wide amplified fragment length polymorphism (AFLP) under unlikely Hardy-Weinberg assumptions, this study investigates 124 individuals within a continuous population of the autopolyploid Biscutella laevigata (Brassicaceae). Fine-scale spatial genetic structure was strong and the mosaic-like distribution of AFLP genotypes was consistently associated with habitat factors, even when controlled for geographical distances. The use of multivariate analyses enabled separation of the factors responsible for the repartition of the genetic variance and revealed a composite effect of isolation by distance, phenological divergence and local adaptation to habitats characterised by different solar radiation regimes. These results suggest that the immigrant inviability barrier facilitated the maintenance of adapted subpopulations to distinct environmental conditions at the local scale.     

Machine learning identifies specific habitats associated with genetic connectivity in Hyla squirella

The goal of this study was to identify and differentiate the influence of multiple habitat types that span a spectrum of suitability for Hyla squirella, a widespread frog species that occurs in a broad range of habitat types. We collected microsatellite data from 675 samples representing 20 localities from the southeastern USA and used machine‐learning methodologies to identify significant habitat features associated with genetic structure. In simulation, we confirm that our machine‐learning algorithm can successfully identify landscape features responsible for generating between‐population genetic differentiation, suggesting that it can be a useful hypothesis‐generating tool for landscape genetics. In our study system, we found that H. squirella were spatially structured and models including specific habitat types (i.e. upland oak forest and urbanization) consistently explained more variation in genetic distance (median pR² = 47.78) than spatial distance alone (median pR² = 23.81). Moreover, we estimate the relative importance that spatial distance, upland oak and urbanized habitat have in explaining genetic structure of H. squirella. We discuss how these habitat types may mechanistically facilitate dispersal in H. squirella. This study provides empirical support for the hypothesis that habitat‐use can be an informative correlate of genetic differentiation, even for species that occur in a wide range of habitats.     

Landscape genetics of the blotched tiger salamander (Ambystoma tigrinum melanostictum)

The field of landscape genetics has great potential to identify habitat features that influence population genetic structure. To identify landscape correlates of genetic differentiation in a quantitative fashion, we developed a novel approach using geographical information systems analysis. We present data on blotched tiger salamanders (Ambystoma tigrinum melanostictum) from 10 sites across the northern range of Yellowstone National Park in Montana and Wyoming, USA. We used eight microsatellite loci to analyse population genetic structure. We tested whether landscape variables, including topographical distance, elevation, wetland likelihood, cover type and number of river and stream crossings, were correlated with genetic subdivision (F(ST)). We then compared five hypothetical dispersal routes with a straight-line distance model using two approaches: (i) partial Mantel tests using Akaike's information criterion scores to evaluate model robustness and (ii) the BIOENV procedure, which uses a Spearman rank correlation to determine the combination of environmental variables that best fits the genetic data. Overall, gene flow appears highly restricted among sites, with a global F(ST) of 0.24. While there is a significant isolation-by-distance pattern, incorporating landscape variables substantially improved the fit of the model (from an r2 of 0.3 to 0.8) explaining genetic differentiation. It appears that gene flow follows a straight-line topographic route, with river crossings and open shrub habitat correlated with lower F(ST) and thus, decreased differentiation, while distance and elevation difference appear to increase differentiation. This study demonstrates a general approach that can be used to determine the influence of landscape variables on population genetic structure.     

人为干扰对神农架川金丝猴连通性及遗传多样性的影响

川金丝猴湖北种群是川金丝猴分布的最东缘,现主要分布于神农架国家公园和巴东自然保护区,其作为单独的管理单元,相对较低的遗传多样性、孤立的遗传状态和较小的种群数量,使得该种群面对环境变化脆弱性高。以神农架为研究区域,首先基于物种分布模型获得川金丝猴移动的阻力图层,利用电路理论分析适宜生境连接度;其次,设置道路和旅游两个人为干扰情景,运用最小费用距离模型与电路理论模型,分析道路和旅游活动对川金丝猴移动路径以及对不同亚群内和亚群间连通性的影响;最后,利用CDPOP软件模拟道路对川金丝猴遗传多样性的影响。结果表明:(1)神农架川金丝猴适宜生境分布较集中,生境斑块间存在多个潜在扩散区域;(2)道路改变了川金丝猴往来于南北生境间的最优路径,旅游干扰对川金丝猴移动路径无明显影响,其干扰主要在于川金丝猴对其具有回避性,回避距离为100 m;(3)道路显著降低金猴岭亚群内的连通性,并显著增加了所有亚群间的移动阻力,旅游活动对亚群内和亚群间的连通性无显著影响;(4)川金丝猴种群在150年内的遗传多样性呈下降趋势,道路的影响将加剧遗传多样性下降的程度。基于以上结果本文提出了针对性的保护建议,对于该物种的保护和其他濒危物种的类似研究具有较强的参考价值和借鉴意义。     

Fine scale population genetic structure of pumas in the Intermountain West

In this study, I examined the population genetic structure of subpopulations of pumas (Puma concolor) in Idaho and surrounding states. Patterns of genetic diversity, population structure, levels of inbreeding, and the relationship between genetic differentiation and dispersal distance within and between 15 subpopulations of pumas were compared. Spatial analyses revealed that the Snake River plain was an important barrier to movement between northern and southern regions of Idaho. In addition, subpopulations south of the Snake River plain exhibited lower levels of genetic diversity, higher levels of inbreeding, and a stronger pattern of isolation by distance relative to subpopulations north of the Snake River plain. Lower levels of diversity and restricted gene flow are likely the result of historically lower population sizes in conjunction with more recent changes in habitat use and available dispersal corridors for movement. The subdivision of puma populations north and south of the Snake River plain, along with the patterns of genetic diversity within regions, indicate that landscape features are affecting the population genetic structure of pumas in Idaho. These results indicate that information about the effects of landscape features on the distribution of genetic diversity should be considered when designing plans for the management and conservation of pumas.     

影响披碱草属植物遗传分化和亲缘关系的地理因素分析

用微卫星标记法对产自中国青藏高原和蒙古高原地区9种披碱草属（Elymus L．）植物的40个居群进行遗传变异和亲缘关系分析,居群问遗传相似系数为0．020～0．957,居群间遗传距离与地理因素关系密切。以垂穗披碱草（Elymus nutctns Griseb．）为例,采用二维和三维Mantel检验,对遗传距离和地理因素间的关系（纬度、经度和海拔）作进一步分析,结果表明,地理位置（纬度和经度）是影响垂穗披碱草居群微卫星遗传差异的最重要因素（R=0．478;P＜0．05）,其次是海拔（R=0．258;P＞0．05）。单一地理位置因素（纬度或经度）与遗传距离的相关性并不明显,但从单一地理位置因素与遗传距离的相关系数、可信度区间来看,纬度对遗传分化的影响程度大于经度。     

Patterning of geographic variation in Middle Pleistocene Homo frontal bone morphology

A quantitative assessment of the frontal bone morphology of a sample of Middle Pleistocene hominins was undertaken in order to address questions regarding their population structure and evolutionary history. Outline tracings of the frontal bones of forty-seven fossil crania were obtained, and size-standardized measurements were then computed using an Elliptical Fourier analysis of these tracings. Principal component scores of the Fourier harmonic amplitudes were derived and served as a quantitative representation of the morphology of the frontal bone. Morphological, geographical, and temporal distance matrices were then constructed between each pair of fossils. A partial Mantel matrix correlation test was performed between morphological and geographical distance matrices, controlling for temporal distance, in order to determine if the pattern of geographical differentiation in features of the frontal bone of mid-Pleistocene Homo followed that of an isolation-by-distance model of population structure. The results of the partial Mantel tests indicate that the overall patterning of differentiation in the features of the frontal bone cannot best be explained by a population structure shaped by isolation-by-distance. Additionally, various aspects of the frontal bone quantified here follow different patterns of geographical differentiation, suggesting that a mosaic pattern of evolution holds true for characters within one cranial region and not just for those between regions.     

Landscape resistance and habitat combine to provide an optimal model of genetic structure and connectivity at the range margin of a small mammal

We evaluated the effect of habitat and landscape characteristics on the population genetic structure of the white‐footed mouse. We develop a new approach that uses numerical optimization to define a model that combines site differences and landscape resistance to explain the genetic differentiation between mouse populations inhabiting forest patches in southern Québec. We used ecological distance computed from resistance surfaces with Circuitscape to infer the effect of the landscape matrix on gene flow. We calculated site differences using a site index of habitat characteristics. A model that combined site differences and resistance distances explained a high proportion of the variance in genetic differentiation and outperformed models that used geographical distance alone. Urban and agriculture‐related land uses were, respectively, the most and the least resistant landscape features influencing gene flow. Our method detected the effect of rivers and highways as highly resistant linear barriers. The density of grass and shrubs on the ground best explained the variation in the site index of habitat characteristics. Our model indicates that movement of white‐footed mouse in this region is constrained along routes of low resistance. Our approach can generate models that may improve predictions of future northward range expansion of this small mammal.     

Effects of landscape and history on diversification of a montane, stream-breeding amphibian

Aim The aim of this study was to understand the roles of landscape features in shaping patterns of contemporary and historical genetic diversification among populations of the Andean tree frog (Hypsiboas andinus) across spatial scales. Location Andes mountains, north‐western Argentina, South America. Methods Mitochondrial DNA control region sequences were utilized to assess genetic differentiation among populations and calculate population pair‐wise genetic distances. Three models of movement, namely traditional straight‐line distance and two effective distances based on habitat classification, were examined to determine which of these explained the most variation in pair‐wise population genetic differentiation. The two habitat classifications were based on digital vegetation and hydrology layers that were generated from a 90‐m resolution digital elevation model (DEM) and known relationships between elevation and habitat. Mantel tests were conducted to test for correlations between geographic and genetic distance matrices and to estimate the percentage variation explained by each type of geographic distance. To investigate the location of possible barriers to gene flow, we used Monmonier’s maximum difference algorithm as implemented in barrier 2.2. Results At both geographic scales, effective distances explained more variation in genetic differentiation than did straight‐line distance. The least‐cost distances based on the simple classification performed better than the more detailed habitat classification. We controlled for the effects of historical range fragmentation determined from previous nested clade analyses, and therefore evaluated the effect of different distances on the genetic variation attributable to more recent factors. Effective distances identified populations that were highly divergent as a result of isolation in unsuitable habitats. The proposed locations of barriers to gene flow identified using Monmonier’s maximum difference algorithm corresponded well with earlier analyses and supported findings from our partial Mantel tests. Main conclusions Our results indicate that landscape features have been important in both historical and contemporary genetic structuring of populations of H. andinus at both large and small spatial scales. A landscape genetic perspective offers novel insights not provided by traditional phylogeographic studies: (1) effective distances can better explain patterns of differentiation in populations, especially in heterogeneous landscapes where barriers to dispersal may be common; and (2) least‐cost path analysis can help to identify corridors of movement between populations that are biologically more realistic.     

Landscape composition modulates population genetic structure of Eriosoma lanigerum (Hausmann) on Malus domestica Borkh in central Chile

Landscape genetics have been particularly relevant when assessing the influence of landscape characteristics on the genetic variability and the identification of barriers to gene flow. Linking current practices of area-wide pest management information on pest population genetics and geographical barriers would increase the efficiency of these programs. The woolly apple aphid, Eriosoma lanigerum (Hausmann), an important pest of apple orchards worldwide, was collected on apple trees (Malus domestica Borkh) from different locations in a 400 km north-south transect trough central Chile. In order to determine if there was population structure, diversity and flow were assessed. A total of 215 individuals from these locations were analysed using Inter Simple Sequence Repeat (ISSR) markers. Four ISSR primers generated a total of 114 polymorphic loci. The percentage of molecular variation among locations was 18%. As the algorithm used by structure may be poorly suited for inferring the number of genetic clusters in a data set that has an IBD relationship, the number of genetic clusters in the samples was also analyzed using a Bayesian clustering method implemented in software Baps version 4.14. We inferred the presence of four genetic clusters in the study region. Clustering of individuals followed a pattern explained by some geographical barriers. Using partial Mantel tests, we detected barriers to gene flow other than distance, created by a combination of main rivers and mountains. Although landscape genetics are rarely used in pest management, our results suggest that these tools may be suitable for the design of area-wide pest management programs.     

Streams over mountains: influence of riparian connectivity on gene flow in the Pacific jumping mouse (Zapus trinotatus)

In species affiliated with heterogeneous habitat, we expect gene flow to be restricted due to constraints placed on individual movement by habitat boundaries. This is likely to impact both individual dispersal and connectivity between populations. In this study, a GIS-based landscape genetics approach was used, in combination with fine-scale spatial autocorrelation analysis and the estimation of recent intersubpopulation migration rates, to infer patterns of dispersal and migration in the riparian-affiliated Pacific jumping mouse (Zapus trinotatus). A total of 228 individuals were sampled from nine subpopulations across a system of three rivers and genotyped at eight microsatellite loci. Significant spatial autocorrelation among individuals revealed a pattern of fine-scale spatial genetic structure indicative of limited dispersal. Geographical distances between pairwise subpopulations were defined following four criteria: (i) Euclidean distance, and three landscape-specific distances, (ii) river distance (distance travelled along the river only), (iii) overland distance (similar to Euclidean, but includes elevation), and (iv) habitat-path distance (a least-cost path distance that models movement along habitat pathways). Pairwise Mantel tests were used to test for a correlation between genetic distance and each of the geographical distances. Significant correlations were found between genetic distance and both the overland and habitat-path distances; however, the correlation with habitat-path distance was stronger. Lastly, estimates of recent migration rates revealed that migration occurs not only within drainages but also across large topographic barriers. These results suggest that patterns of dispersal and migration in Pacific jumping mice are largely determined by habitat connectivity.     

A range-wide genetic bottleneck overwhelms contemporary landscape factors and local abundance in shaping genetic patterns of an alpine butterfly (Lepidoptera: Pieridae: Colias behrii)

Spatial and environmental heterogeneity are major factors in structuring species distributions in alpine landscapes. These landscapes have also been affected by glacial advances and retreats, causing alpine taxa to undergo range shifts and demographic changes. These nonequilibrium population dynamics have the potential to obscure the effects of environmental factors on the distribution of genetic variation. Here, we investigate how demographic change and environmental factors influence genetic variation in the alpine butterfly Colias behrii. Data from 14 microsatellite loci provide evidence of bottlenecks in all population samples. We test several alternative models of demography using approximate Bayesian computation (ABC), with the results favouring a model in which a recent bottleneck precedes rapid population growth. Applying independent calibrations to microsatellite loci and a nuclear gene, we estimate that this bottleneck affected both northern and southern populations 531–281 years ago, coinciding with a period of global cooling. Using regression approaches, we attempt to separate the effects of population structure, geographical distance and landscape on patterns of population genetic differentiation. Only 40% of the variation in F_ST is explained by these models, with geographical distance and least‐cost distance among meadow patches selected as the best predictors. Various measures of genetic diversity within populations are also decoupled from estimates of local abundance and habitat patch characteristics. Our results demonstrate that demographic change can have a disproportionate influence on genetic diversity in alpine species, contrasting with other studies that suggest landscape features control contemporary demographic processes in high‐elevation environments.     

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.