The effects of contemporary processes in maintaining the genetic structure of western song sparrows (Melospiza melodia)

Historic events and contemporary processes work in concert to create and maintain geographically partitioned variation and are instrumental in the generation of biodiversity. We sought to gain a better understanding of how contemporary processes such as movement and isolation influence the genetic structure of widely distributed vagile species such as birds. Song sparrows (Melospiza melodia) in western North America provide a natural system for examining the genetics of populations that have different patterns of geographic isolation and migratory behavior. We examined the population genetics of 576 song sparrows from 23 populations using seven microsatellite loci to assess genetic differentiation among populations and to estimate the effects of drift and immigration (gene flow) on each population. Sedentary, isolated populations were characterized by low levels of immigration and high levels of genetic drift, whereas those populations less isolated displayed signals of high gene flow and little differentiation from other populations. Contemporary dispersal rates from migratory populations, estimated by assignment test, were higher and occurred over larger distances than dispersal from sedentary populations but were also probably too low to counter the effects of drift in most populations. We suggest that geographic isolation and limited gene flow facilitated by migratory behavior are responsible for maintaining observed levels of differentiation among Pacific coastal song sparrow populations.     

Estimates of gene flow, genetic substructure and population heterogeneity in bracken (Pteridium aquilinum)

Bracken [ Pteridium aquilinum (L.) Kuhn] is a cosmopolitan species and is a noxious weed in many areas. Because of its abundance, particularly in Britain, bracken affords an ideal system for investigating various aspects of population genetics and evolution. High mobility of dispersal units (spores) suggests that rates of gene flow among distant populations should be high. Gene flow is a major evolutionary force that influences the genetic structure of populations. To examine the effects of gene flow on population heterogeneity and population substructuring in bracken, starch gel electrophoresis of enzymes was used to provide the necessary genetic database. Allele frequency data at 21 loci were obtained for seven British populations, one Majorcan and one from the eastern United States. A model was employed to estimate the amount of gene flow ( Nm ) at several levels. Gene flow among British populations was extremely high ( Nm = 36.51), one of the highest estimates reported for plants. Among eight European populations gene flow was lower (but still considered high) at Nm = 2.47. Trans-Atlantic gene flow was low ( Nm = 0.0926).
F -statistics were used to assess population heterogeneity and substructuring. The data indicate that, compared with other species, there is very little genetic differentiation among British populations of bracken. Indeed, it appears that the whole island is behaving as a single randommating population. This result is consistent with high levels of gene flow. Only one population (on the Isle of Arran) showed statistically significant genetic substructuring. Habitat heterogeneity on the island and age structure are hypothesized as possible causes of this result.
The data reported here support previous studies demonstrating that bracken is genetically polymorphic and is an outcrossing species.     

Gene flow among small populations of a self-incompatible plant: an interaction between demography and genetics

We assessed the effects of population size and genetic relatedness on rates of pollen gene flow into experimental populations of the insect-pollinated, self-incompatible plant Raphanus sativus. We created synthetic populations of sizes 2, 5, 10, and 20 with three genetic structures (full siblings, half siblings, and unrelated plants). Following pollination in a natural setting, we conducted a simple paternity exclusion analysis using the allozyme genotypes of progeny to measure apparent gene flow and Monte Carlo simulations to estimate total gene flow. Estimates of apparent pollen gene flow rates ranged from 0 to 100% and were similar in rank to estimates of total gene flow. There were significant effects of population size and relatedness on the rate of apparent gene flow, and there were significant population size by relatedness interactions. Populations of size 2 had higher gene flow rates than larger populations, gene flow being negatively associated with the level of cross-compatibility (as measured by hand pollinations). Gene flow into populations of size 2 was also negatively associated with the distance to the nearest population of size 10 or 20. These results suggest that interactions among demography (population size), genetics (cross-compatibility), and ecology (pollinator behavior) are important influences on pollen gene flow rates into small plant populations.     

Reconstructing recent divergence: evaluating nonequilibrium population structure in New Zealand chinook salmon

Newly established or perturbed populations are often the focus of conservation concerns but they pose special challenges for population genetics because drift?migration equilibrium is unlikely. To advance our understanding of the evolution of such populations, we investigated structure and gene flow among populations of chinook salmon that formed via natural straying following introduction to New Zealand in the early 1900s. We examined 11 microsatellite loci from samples collected in several sites and years to address two questions: (i) what population differentiation has arisen in the ≈ 30 generations since salmon were introduced to New Zealand, relative to temporal variation within populations; and (ii) what are the approximate effective population sizes and amounts of gene flow in these populations? These questions are routinely addressed in studies of indigenous populations, but less often in the case of new populations and rarely with consideration of equilibrium assumptions. We show that despite the recent introduction, continued gene flow and high temporal variability among samples, detectable population structure has arisen among the New Zealand populations, consistent with their colonization pattern and isolation by geographical distance. Furthermore, we use simple individual‐based simulations and estimates of effective population sizes to estimate the effective gene flow among drainages under likely nonequilibrium conditions. Similar methodology may be broadly applicable to other studies of population structure and phenotypic evolution under similar nonequilibrium, high gene flow conditions.     

Evaluating Genetic Diversity for the Conservation of the Threatened Galapagos Endemic Calandrinia galapagosa (Portulacaceae)

Calandrinia galapagosa (Portulacaceae) is a threatened plant endemic to San Cristobal Island, Galapagos. The species is threatened by goat herbivory: less than 1300 individuals remain in nine populations, three of which have been fenced for protection. It is not known whether these three populations are genetically representative of the species at large or whether the three distinct color morphs exhibited by the species represent different taxonomic units. We measured the frequency of 113 polymorphic amplified fragment length polymorphism markers in 189 specimens representing three different phenotypes from all nine populations to relate the genetic structure of the species to both its geological history and morphological variation. The analysis revealed high diversity within populations and low, but significant, differentiation between populations. In addition, no genetic diversification was observed across color morphs. We conclude that (1) regular gene flow has historically occurred between populations, and (2) individual color morphs are not separate taxa. We identify another population that should be fenced to ensure maximal protection of genetic variation. Our research serves as a case study of using population molecular genetics for optimizing protection strategies for rare plants in tropical archipelagoes affected by introduced herbivores.     

Simulating natural selection in landscape genetics

Linking landscape effects to key evolutionary processes through individual organism movement and natural selection is essential to provide a foundation for evolutionary landscape genetics. Of particular importance is determining how spatially-explicit, individual-based models differ from classic population genetics and evolutionary ecology models based on ideal panmictic populations in an allopatric setting in their predictions of population structure and frequency of fixation of adaptive alleles. We explore initial applications of a spatially-explicit, individual-based evolutionary landscape genetics program that incorporates all factors--mutation, gene flow, genetic drift and selection--that affect the frequency of an allele in a population. We incorporate natural selection by imposing differential survival rates defined by local relative fitness values on a landscape. Selection coefficients thus can vary not only for genotypes, but also in space as functions of local environmental variability. This simulator enables coupling of gene flow (governed by resistance surfaces), with natural selection (governed by selection surfaces). We validate the individual-based simulations under Wright-Fisher assumptions. We show that under isolation-by-distance processes, there are deviations in the rate of change and equilibrium values of allele frequency. The program provides a valuable tool (cdpop v1.0; http://cel.dbs.umt.edu/software/CDPOP/) for the study of evolutionary landscape genetics that allows explicit evaluation of the interactions between gene flow and selection in complex landscapes.     

Extinction-colonization dynamics structure genetic variation of spotted sunfish (Lepomis punctatus) in the Florida Everglades

The population genetics of aquatic animals in the Florida Everglades may be strongly influenced by extinction and colonization dynamics. We combined analyses of allozyme and microsatellite loci to test the hypothesis that two levels of population structure are present for spotted sunfish (Pisces: Centrarchidae: Lepomis punctatus) inhabiting the Everglades. We hypothesized that annual cycles of marsh dry-down increase local-scale genetic variation through a process of local extinction and colonization; we hypothesized that barriers to gene flow by levee/canal systems create a second, regional level of genetic variation. In 1996 and 1997, we sampled spotted sunfish from 11 Everglades sites that were distributed in three regions separated by levees. We documented patterns of genetic variation at 7 polymorphic allozyme loci and 5 polymorphic microsatellite loci. Most genetic variation was present among local populations, according to both types of genetic markers. Furthermore, samples from marsh sites were heterogeneous, while those from canals were not. These data supported our hypothesis that dry-down events and local population dynamics in the marsh have a significant effect on population genetic structure of spotted sunfish. We found no support for our hypothesis that water-management structures superimpose a second level of genetic structure on this species, possibly because canals obscure historical structure by facilitating gene flow or because the complete canal system has been in place for fewer than 20 generations of this species. Our data suggests a continent-island (canal-marsh) structure of populations with high gene flow among regions and recurrent mixing in marshes from canal and creek habitats.     

Landscape genetic connectivity in European wildcat (Felis silvestris silvestris): a matter of food,shelters and demographic status of populations

Understanding landscape impacts on gene flow is necessary to plan comprehensive management and conservation strategies of both the species of interest and its habitat. Nevertheless, only a few studies have focused on the landscape genetic connectivity of the European wildcat, an umbrella species whose conservation allows the preservation of numerous other species and habitat types. We applied population and landscape genetics approaches, using genotypes at 30 microsatellites from 232 genetically-identified wildcats to determine if, and how, landscape impacted gene flow throughout France. Analyses were performed independently within two population patches: the historical north-eastern patch and the central patch considered as the colonization front. Our results showed that gene flow occurred at large spatial scales but also revealed significant spatial genetic structures within population patches. In both population patches, arable areas, pastures and permanent grasslands and lowly fragmented forested areas were permeable to gene flow, suggesting that shelters and dietary resources are among the most important parameters for French wildcat landscape connectivity, while distance to forest had no detectable effect. Anthropized areas appeared highly resistant in the north-eastern patch but highly permeable in the central patch, suggesting that different behaviours can be observed according to the demographic context in which populations are found. In line with this hypothesis, spatial distribution of genetic variability seemed uneven in the north-eastern patch and more clinal in the central patch. Overall, our results highlighted that European wildcat might be a habitat generalist species and also the importance of performing spatial replication in landscape genetics studies.

     

Population structure and genetic diversity of black redhorse (Moxostoma duquesnei) in a highly fragmented watershed

Dams have the potential to affect population size and connectivity, reduce genetic diversity, and increase genetic differences among isolated riverine fish populations. Previous research has reported adverse effects on the distribution and demographics of black redhorse (Moxostoma duquesnei), a threatened fish species in Canada. However, effects on genetic diversity and population structure are unknown. We used microsatellite DNA markers to assess the number of genetic populations in the Grand River (Ontario) and to test whether dams have resulted in a loss of genetic diversity and increased genetic differentiation among populations. Three hundred and seventy-seven individuals from eight Grand River sites were genotyped at eight microsatellite loci. Measures of genetic diversity were moderately high and not significantly different among populations; strong evidence of recent population bottlenecks was not detected. Pairwise F_ST and exact tests identified weak (global F_ST = 0.011) but statistically significant population structure, although little population structuring was detected using either genetic distances or an individual-based clustering method. Neither geographic distance nor the number of intervening dams were correlated with pairwise differences among populations. Tests for regional equilibrium indicate that Grand River populations were either in equilibrium between gene flow and genetic drift or that gene flow is more influential than drift. While studies on other species have identified strong dam-related effects on genetic diversity and population structure, this study suggests that barrier permeability, river fragment length and the ecological characteristics of affected species can counterbalance dam-related effects.     

Gene flow and effective population size in Lepanthes (Orchidaceae): a case for genetic drift

Genetic drift can play an important role in population differentiation, particularly when effective population sizes are small and gene flow is limited. Such conditions are suspected to be common in the species-rich Orchidaceae. We investigated the likelihood of genetic drift in natural populations of three endemic species of Lepanthes (Orchidaceae) from Puerto Rico. We estimated effective population size, Ne, using three ecologically based methods. Two of the three estimates were based on variance in reproductive potential and the third was based on coalescence time. All estimates of Ne were usually <40% of the standing population size, resulting in values of <20 individuals per population. Based on starch gel electrophoresis of isozymes, Nm estimates suggest restricted gene flow among populations in the range of one or less successful migrant per generation. Genetic differentiation among populations is expected under such conditions from random genetic drift. Indeed we observed high genetic differentiation among populations (L. rubripetala, F_ST> G_ST, 6; θ.248, 0.266, 0.293; L. rupestris, 0.148, 0.169, 0.138; L. eltoroensis, 0.251, 0.219, 0.218, respectively). Genetic drift is likely to be important for population differentiation in Lepanthes as a result of small effective population sizes and restricted gene flow.     

Historical sketch of Slovak Haban (Hutterite) population based on autosomal STR analysis

According to the Hutterite chronicles, the Habans arrived from Austrian Tyrol, Switzerland, and northernmost Italy and stayed in four regions of Slovakia (Soboti?te, Vel'ké Leváre, Moravsky Sv?ty Ján, Tren?ín). There are some communities in western Slovakia that retained their Haban cultural identity and still identify themselves as descendents of the Hutterite population with their own specific customs. Slovak Habans are typical founder population with significant social isolation for which high degree of inbreeding is typical. Present study investigated STR polymorphisms as a powerful genetic tool for population genetic studies. The aim was to perform a comparative, population genetic study based on 15 STR loci widely used in forensic genetics, of the Haban population, the Slovak majority population and the population of Tyrol. We analyzed allele frequencies and other statistical parameters in three selected populations in order to identify groups of specific ethnic origin and establish their genetic relationship. The data set included 110 unrelated Habans and 201 unrelated individuals from the Slovak majority population, as well as allelic frequencies for the population of Austrian Tyrol available in the literature. Population pairwise FST values used as a short term genetic distance between populations showed significant differentiation between the Habans and both reference populations (FST=0.0025 and 0.0042 for comparison with the Slovaks and Austrians, respectively; p<10(-3)). The Slovak Hutterites were demonstrated to be genetically distinct and more closely related to their geographic neighbors than to their historical ancestral population, which may be at least partially explained by gene flow between neighboring Haban and Slovak populations.     

Gene flow among Cydia pomonella (Lepidoptera: Tortricidae) geographic and host populations in South Africa

Information on gene flow among geographic and host populations of C. pomonella (L.) (Lepidoptera: Tortricidae) in South Africa is lacking, despite the importance of these measures for the success of control practices such as chemical control and sterile insect release, which are affected by the amount of gene flow among populations. Therefore, populations collected from nine geographically distant regions in South Africa from apples, pears, and stone fruit were compared using amplified fragment length polymorphism with five selective primer pairs. Results showed that although populations from different hosts were not genetically differentiated, significant evidence for population substructure was apparent between geographic populations. Over local scales, it was possible to distinguish between populations collected from orchards situated <1 km apart. These results suggest that although extensive gene flow occurs among populations from different hosts, gene flow among local geographic C. pomonella populations may be limited and is explained in terms of limited moth flight, the relative isolation of pome fruit production areas, and the absence of wild hosts.     

High levels of gene flow and heterozygote excess characterize Rhizoctonia solani AG-1 IA (Thanatephorus cucumeris) from Texas

To date, much of the genetics of the basidiomycete Thanatephorus cucumeris (anamorph = Rhizoctonia solani) remains unknown. Here, we present a population genetics study using codominant markers to augment laboratory analyses. Seven single-copy nuclear RFLP markers were used to examine 182 isolates of Rhizoctonia solani AG-1 IA collected from six commercial rice fields in Texas. Thirty-six multilocus RFLP genotypes were identified. Population subdivision analyses indicated a high degree of gene flow/migration between the six geographic populations. Tests for Hardy-Weinberg equilibrium (HWE) among the 36 multilocus RFLP genotypes revealed that four of the seven loci did not significantly differ from HWE. Subsequent analysis demonstrated that departures from HWE at the three remaining loci were due to an excess of heterozygotes. Data presented here suggest that R. solani AG-1 IA is actively outbreeding (heterothallic). Possible explanations for heterozygote excess, which was observed at all seven RFLP loci, are discussed.     

Population genetics of the hispid cotton rat (Sigmodon hispidus): patterns of genetic diversity at the major histocompatibility complex

The hispid cotton rat, Sigmodon hispidus, is a common rodent widely distributed across the southern USA and south into South America. To characterize major histocompatibility complex (MHC) diversity in this species and to elucidate large-scale patterns of genetic partitioning, we examined MHC genetic variability within and among 13 localities, including a disjunct population in Arizona and a population from Costa Rica that may represent an undescribed species. We also tested the hypothesis that populations within the USA are at equilibrium with regard to gene flow and genetic drift, resulting in isolation-by-distance. Using single-strand conformation polymorphism (SSCP) analysis we identified 25 alleles from 246 individuals. Gene diversity within populations ranged from 0.000 to 0.908. Analysis of molecular variance (AMOVA) revealed that 83.7% of observed variation was accounted for by within-population diversity and 16.3% was accounted for by among-population divergence. The disjunct population in Arizona was fixed for a single allele. The Costa Rican population was quite divergent based on allelic composition and was the only population with unique alleles. Within the main portion of the geographical distribution of S. hispidus in the USA there was considerable divergence among some populations; however, there was no significant pattern of isolation-by-distance overall (P = 0.090). Based on the significant divergence of the only sampled population to its east, the Mississippi River appears to represent a substantial barrier to gene flow.     

Unconstrained gene flow between populations of a widespread epiphytic lichen Usnea subfloridana (Parmeliaceae,Ascomycota) in Estonia

Few studies have investigated the genetic diversity of populations of common and widespread lichenized fungi using microsatellite markers, especially the relationships between different measures of genetic diversity and environmental heterogeneity. The main aim of our study was to investigate the population genetics of a widespread and mainly clonally reproducing Usnea subfloridana at the landscape scale, focusing on the comparison of lichen populations within hemiboreal forest stands. Particular attention has been paid to the genetic differentiation of lichen populations in two geographically distinct regions in Estonia and the relationships between forest characteristics and measures of genetic diversity. We genotyped 578 Usnea thalli from eleven lichen populations using seven specific fungal microsatellite markers. Measures of genetic diversity (allelic richness, Shannon's information index, Nei's unbiased genetic diversity, clonal diversity, the number of multilocus genotypes, the number of private alleles, and the minimum number of colonization events) were calculated and compared between Usnea populations. Shared haplotypes, gene flow and AMOVA analyses suggest that unconstrained gene flow and exchange of multilocus genotypes exist between the two geographically remote regions in Estonia. Stand age, mean circumference of the host tree, size of forest site and tree species composition did not show any significant influence on allelic richness, Shannon's information index, Nei's unbiased genetic diversity, clonal diversity, the number of private alleles, and the minimum number of colonization events of U. subfloridana populations. Therefore it was concluded that other factors of habitat heterogeneity could probably have a more significant effect on population genetics of U. subfloridana populations.     

Population genetic structure of Bemisia tabaci (Hemiptera: Aleyrodidae) utilizing microsatellite markers

We aimed to characterize the population genetic structure within and among five Bemisia tabaci (Gennadius) populations collected from different host plants and geographic regions by using microssatelites as a molecular marker. Each population was represented by 19 specimens. The host plants and geographic origins of these populations were described as follows: Pop 1: Squash Barreiras (BA); Pop 2: Cotton Barreiras (BA); Pop 3: Soybean Campinas (SP); Pop 4: Tomato Cruz das Almas (BA); and Pop 5: Soybean Rondonópolis (MT). Six polymorphic loci were observed, which discriminated 31 different alleles in the studied populations, with a mean number of alleles per population of 3.30 (2.67 - 4.00). Using Fisher's Exact test, it was observed that at least three populations were in Hardy-Weinberg equilibrium for most of the studied loci (six). The dendrogram (UPGMA) separated populations into groups mainly related to the geographic origin of the samples. Only population 5 differed from the others at a 0.15 distance (74.5% group consistency). The most similar populations were 1 and 2, with a 0.01 distance (65.3%). This is in agreement with their geographic origins and it was not consistent with host specificity. The results suggest considerable gene flow (7.3%) among all whitefly populations and indicate that a better understanding of the gene flow in populations of B. tabaci associated with different hosts is required for the management of this insect.     

Population differentiation and nuclear gene flow in the Dominican anole (Anolis oculatus)

Allele frequency data from nuclear microsatellite loci were used to investigate patterns of nuclear gene flow and population structure in the morphologically variable Dominican anole (Anolis oculatus). All six loci used proved to be highly polymorphic, with an average of 18.8 alleles per locus. Test for Hardy-Weinberg equilibrium revealed small numbers of heterozygote deficiencies at single loci in single populations and consistent patterns of increasingly significant heterozygote deficiency in global tests across populations and loci. No significant relationship between FST and patristic distances estimated from mitochondrial DNA sequences was detected and estimates of FIS were significantly higher in females than in males, indicating that gene flow may be sex-biased and mediated mainly by male migration. A highly significant correlation between linearized FST and loge (geographical distance) indicates that geographical proximity is a significant factor in the genetic structure of A. oculatus populations. However, levels of gene flow between morphologically differentiated parapatric populations are frequently seen to be relatively high. This supports the hypothesis of natural selection being the driving force behind the development and maintenance of morphological variation and shows that adaptive differentiation may be maintained despite the homogenizing influence of gene flow. Generally, the morphologically variable populations of A. oculatus are seen to be poor candidates for in situ speciation, but an exceptional case on the west coast of Dominica indicates that isolation resulting from vicariant events may lead to rapid differentiation at both mitochondrial and nuclear loci. This provides a possible mechanism for anole speciation on other Caribbean islands.     

Population genetics of Wolbachia‐infected,parthenogenetic and uninfected,sexual populations of Tetrastichus coeruleus (Hymenoptera: Eulophidae)

Wolbachia are endosymbiotic bacteria known to manipulate the reproduction of their hosts. These manipulations are expected to have consequences on the population genetics of the host, such as heterozygosity levels, genetic diversity and gene flow. The parasitoid wasp Tetrastichus coeruleus has populations that are infected with parthenogenesis‐inducing Wolbachia and populations that are not infected. We studied the population genetics of T. coeruleus between and within Wolbachia‐infected and uninfected populations, using nuclear microsatellites and mitochondrial DNA. We expected reduced genetic diversity in both DNA types in infected populations. However, migration and gene flow could introduce new DNA variants into populations. We therefore paid special attention to individuals with unexpected (genetic) characteristics. Based on nuclear and mitochondrial DNA, two genetic clusters were evident: a thelytokous cluster containing all Wolbachia‐infected, parthenogenetic populations and an arrhenotokous cluster containing all uninfected, sexual populations. Nuclear and mitochondrial DNA did not exhibit concordant patterns of variation, although there was reduced genetic diversity in infected populations for both DNA types. Within the thelytokous cluster, there was nuclear DNA variation, but no mitochondrial DNA variation. This nuclear DNA variation may be explained by occasional sex between infected females and males, by horizontal transmission of Wolbachia, and/or by novel mutations. Several females from thelytokous populations were uninfected and/or heterozygous for microsatellite loci. These unexpected characteristics may be explained by migration, by inefficient transmission of Wolbachia, by horizontal transmission of Wolbachia, and/or by novel mutations. However, migration has not prevented the build‐up of considerable genetic differentiation between thelytokous and arrhenotokous populations.     

Using molecular markers with high mutation rates to obtain estimates of relative population size and to distinguish the effects of gene flow and mutation: a demonstration using data from endemic Mauritian skinks

We propose a method of analysing genetic data to obtain separate estimates of the size (N(p)) and migration rate (m(p)) for the sampled populations, without precise prior knowledge of mutation rates at each locus ( micro(L)). The effects of migration and mutation can be distinguished because high migration has the effect of reducing genetic differentiation across all loci, whereas a high mutation rate will only affect the locus in question. The method also takes account of any differences between the spectra of immigrant alleles and of new mutant alleles. If the genetic data come from a range of population sizes, and the loci have a range of mutation rates, it is possible to estimate the relative sizes of the different N(p) values, and likewise the m(p) and the micro(L). Microsatellite loci may also be particularly appropriate because loci with a high mutation rate can reach mutation-drift-migration equilibrium more quickly, and because the spectra of mutants arriving in a population can be particularly distinct from the immigrants. We demonstrate this principle using a microsatellite data set from Mauritian skinks. The method identifies low gene flow between a putative new species and populations of its sister species, whereas the differentiation of two other populations is attributed to small population size. These distinct interpretations were not readily apparent from conventional measures of genetic differentiation and gene diversity. When the method is evaluated using simulated data sets, it correctly distinguishes low gene flow from small population size. Loci that are not at mutation-migration-drift equilibrium can distort the parameter estimates slightly. We discuss strategies for detecting and overcoming this effect.     

Effect of landscape features on genetic structure of the goitered gazelle (<Emphasis Type="Italic">Gazella subgutturosa</Emphasis>) in Central Iran

The populations of goitered gazelle suffered significant decline due to natural and anthropogenic factors over the last century. Investigating the effects of barriers on gene flow among the remaining populations is vital for conservation planning. Here we adopted a landscape genetics approach to evaluate the genetic structure of the goitered gazelle in Central Iran and the effects of landscape features on gene flow using 15 polymorphic microsatellite loci. Spatial autocorrelation, isolation by distance (IBD) and isolation by resistance (IBR) models were used to elucidate the effects of landscape features on the genetic structure. Ecological modeling was used to construct landscape permeability and resistance map using 12 ecogeographical variables. Bayesian algorithms revealed three genetically homogeneous groups and restricted dispersal pattern in the six populations. The IBD and spatial autocorrelation revealed a pattern of decreasing relatedness with increasing distance. The distribution of potential habitats was strongly correlated with bioclimatic factors, vegetation type, and elevation. Resistance distances and graph theory were significantly related with variation in genetic structure, suggesting that gazelles are affected by landscape composition. The IBD showed greater impact on genetic structure than IBR. The Mantel and partial Mantel tests indicated low but non-significant effects of anthropogenic barriers on observed genetic structure. We concluded that a combination of geographic distance, landscape resistance, and anthropogenic factors are affecting the genetic structure and gene flow of populations. Future road construction might impede connectivity and gene exchange of populations. Conservation measures on this vulnerable species should consider some isolated population as separate management units.