The effect of a barrier to gene flow on patterns of geographic variation.

Explicit formulae are given for the effects of a barrier to gene flow on random fluctuations in allele frequency; these formulae can also be seen as generating functions for the distribution of coalescence times. The formulae are derived using a continuous diffusion approximation, which is accurate over all but very small spatial scales. The continuous approximation is confirmed by comparison with the exact solution to the stepping stone model. In both one and two spatial dimensions, the variance of fluctuations in allele frequencies increases near the barrier; when the barrier is very strong, the variance doubles. However, the effect on fluctuations close to the barrier is much greater when the population is spread over two spatial dimensions than when it occupies a linear, one-dimensional habitat: barriers of strength comparable with the dispersal range (B approximately equal to sigma) can have an appreciable effect in two dimensions, whereas only barriers with strength comparable with the characteristic scale (B approximately equal to L=sigma/sqrt{2mu}) are significant in one dimension (mu is the rate of mutation or long-range dispersal). Thus, in a two-dimensional population, barriers to gene flow can be detected through their effect on the spatial pattern of genetic marker alleles.     

A conceptual framework for the spatial analysis of landscape genetic data

Understanding how landscape heterogeneity constrains gene flow and the spread of adaptive genetic variation is important for biological conservation given current global change. However, the integration of population genetics, landscape ecology and spatial statistics remains an interdisciplinary challenge at the levels of concepts and methods. We present a conceptual framework to relate the spatial distribution of genetic variation to the processes of gene flow and adaptation as regulated by spatial heterogeneity of the environment, while explicitly considering the spatial and temporal dynamics of landscapes, organisms and their genes. When selecting the appropriate analytical methods, it is necessary to consider the effects of multiple processes and the nature of population genetic data. Our framework relates key landscape genetics questions to four levels of analysis: (i) node-based methods, which model the spatial distribution of alleles at sampling locations (nodes) from local site characteristics; these methods are suitable for modeling adaptive genetic variation while accounting for the presence of spatial autocorrelation. (ii) Link-based methods, which model the probability of gene flow between two patches (link) and relate neutral molecular marker data to landscape heterogeneity; these methods are suitable for modeling neutral genetic variation but are subject to inferential problems, which may be alleviated by reducing links based on a network model of the population. (iii) Neighborhood-based methods, which model the connectivity of a focal patch with all other patches in its local neighborhood; these methods provide a link to metapopulation theory and landscape connectivity modeling and may allow the integration of node- and link-based information, but applications in landscape genetics are still limited. (iv) Boundary-based methods, which delineate genetically homogeneous populations and infer the location of genetic boundaries; these methods are suitable for testing for barrier effects of landscape features in a hypothesis-testing framework. We conclude that the power to detect the effect of landscape heterogeneity on the spatial distribution of genetic variation can be increased by explicit consideration of underlying assumptions and choice of an appropriate analytical approach depending on the research question.     

Gene surfing in expanding populations

Large scale genomic surveys are partly motivated by the idea that the neutral genetic variation of a population may be used to reconstruct its migration history. However, our ability to trace back the colonization pathways of a species from their genetic footprints is limited by our understanding of the genetic consequences of a range expansion. Here, we study, by means of simulations and analytical methods, the neutral dynamics of gene frequencies in an asexual population undergoing a continual range expansion in one dimension. During such a colonization period, lineages can fix at the wave front by means of a "surfing" mechanism [Edmonds, C.A., Lillie, A.S., Cavalli-Sforza, L.L., 2004. Mutations arising in the wave front of an expanding population. Proc. Natl. Acad. Sci. 101, 975-979]. We quantify this phenomenon in terms of (i) the spatial distribution of lineages that reach fixation and, closely related, (ii) the continual loss of genetic diversity (heterozygosity) at the wave front, characterizing the approach to fixation. Our stochastic simulations show that an effective population size can be assigned to the wave that controls the (observable) gradient in heterozygosity left behind the colonization process. This effective population size is markedly higher in the presence of cooperation between individuals ("pushed waves") than when individuals proliferate independently ("pulled waves"), and increases only sub-linearly with deme size. To explain these and other findings, we develop a versatile analytical approach, based on the physics of reaction-diffusion systems, that yields simple predictions for any deterministic population dynamics. Our analytical theory compares well with the simulation results for pushed waves, but is less accurate in the case of pulled waves when stochastic fluctuations in the tip of the wave are important.     

Comparative analysis of the within-population genetic structure in wild cherry (Prunus avium L.) at the self-incompatibility locus and nuclear microsatellites

Gametophytic self-incompatibility (SI) systems in plants exhibit high polymorphism at the SI controlling S-locus because individuals with rare alleles have a higher probability to successfully pollinate other plants than individuals with more frequent alleles. This process, referred to as frequency-dependent selection, is expected to shape number, frequency distribution, and spatial distribution of self-incompatibility alleles in natural populations. We investigated the genetic diversity and the spatial genetic structure within a Prunus avium population at two contrasting gene loci: nuclear microsatellites and the S-locus. The S-locus revealed a higher diversity (15 alleles) than the eight microsatellites (4-12 alleles). Although the frequency distribution of S-alleles differed significantly from the expected equal distribution, the S-locus showed a higher evenness than the microsatellites (Shannon's evenness index for the S-locus: E = 0.91; for the microsatellites: E = 0.48-0.83). Also, highly significant deviations from neutrality were found for the S-locus whereas only minor deviations were found for two of eight microsatellites. A comparison of the frequency distribution of S-alleles in three age-cohorts revealed no significant differences, suggesting that different levels of selection acting on the S-locus or on S-linked sites might also affect the distribution and dynamics of S-alleles. Autocorrelation analysis revealed a weak but significant spatial genetic structure for the multilocus average of the microsatellites and for the S-locus, but could not ascertain differences in the extent of spatial genetic structure between these locus types. An indirect estimate of gene dispersal, which was obtained to explain this spatial genetic pattern, indicated high levels of gene dispersal within our population (sigma(g) = 106 m). This high gene dispersal, which may be partly due to the self-incompatibility system itself, aids the effective gene flow of the microsatellites, thereby decreasing the contrast between the neutral microsatellites and the S-locus.     

Human neutral alpha-glucosidase C: genetic polymorphism including a "null" allele.

We describe a genetic polymorphism of human neutral alpha-glucosidase C, detected in lymphoid cells by a combination of starch gel electrophoresis and isoelectric focusing. The seven phenotypes observed appear to result from the expression of four different alleles. The distribution of the observed phenotypes fits the expected distribution predicted from calculated gene frequencies in Hardy-Weinberg equilibrium. Family studies are consistent with autosomal inheritance of the gene. The product of one of the alleles is unusual in that it is "silent," with an estimated gene frequency of .174 in an outbred white population. Approximately one-third of the population is heterozygous "null." Homozygosity for the allele has not been associated with any obvious disease state. This is the third example of a "null" allele which has a substantial gene frequency in an outbred population but does not appear to result in disease in the homozygous state.     

Gene and Allelic Genealogies at a Gametophytic Self-Incompatibility Locus

The properties of gene and allelic genealogies at a gametophytic self-incompatibility locus in plants have been investigated analytically and checked against extensive numerical simulations. It is found that, as with overdominant loci, there are two genealogical processes with markedly different time scales. First, functionally distinct allelic lines diverge on an extremely long time scale which is inversely related to the mutation rate to new alleles. These alleles show a genealogical structure which is similar, after an appropriate rescaling of time, to that described by the coalescent process for genes at a neutral locus. Second, gene copies sampled within the same functional allelic line show genealogical relationships similar to neutral gene genealogies but on a much shorter time scale, which is on the same order of magnitude as the harmonic mean of the number of gene copies within an allelic line. These results are discussed in relation to data showing trans-specific polymorphisms for alleles at the gametophytic self-incompatibility locus in the Solanaceae. It is shown that population sizes on the order of 4 X 10(5) and a mutation rate per locus per generation as high as 10(-6) could account for estimated allelic divergence times in this family.     

Mating system as a barrier to gene flow

Understanding mating system as one of reproductive isolating barriers remains important although this barrier is classified in a different sense from behavioral, ecological, and mechanical isolating barriers. Selfing enhances incipient speciation while outcrossing facilitates species integrity. Here, I study how mating system affects gene exchanges between genetically diverging species in a hybrid zone. Results show that a predominant selfing species has a greater barrier to selective gene flow than does a predominant outcrossing species. Barrier to neutral gene flow convexly changes with the selfing rate due to linkage disequilibrium, with a maximum at around intermediate selfing rate. Asymmetric transient or steady‐state barriers to neutral gene flow occur between two sides of a hybrid zone when the neutral gene is affected by its linked selective gene whose alternative alleles are adaptive to heterogeneous habitats. Selfing interacts with both a physical barrier and a density‐dependent ecological regulation (a logarithmic model) to strengthen the barriers to neutral and selective gene flow. This theory helps to interpret incipient speciation driven by selfing or to explain the asymmetric gene flow or unequal genomic mixtures between closely related species caused by their asymmetric mating systems in natural hybrid zones.     

Duplications That Suppress and Deletions That Restore Expression from a Chalcone Synthase Multigene Family

Seed coat color in soybean is determined by four alleles of the classically defined / (inhibitor) locus that controls the presence or absence as well as the spatial distribution of anthocyanin pigments in the seed coat. By analyzing spontaneous mutations of the / locus, we demonstrated that the / locus is a region of chalcone synthase (CHS) gene duplications. Paradoxically, deletions of CHS gene sequences allow higher levels of CHS mRNAs and restore pigmentation to the seed coat. The unusual nature of the / locus suggests that its dominant alleles may represent naturally occurring examples of homology-dependent gene silencing and that the spontaneous deletions erase the gene-silencing phenomena. Specifically, mutations from the dominant ii allele (yellow seed coats with pigmented hila) to the recessive i allele (fully pigmented) can be associated with the absence of a 2.3-kb Hindlll fragment that carries CHS4, a member of the multigene CHS family. Seven independent mutations exhibit deletions in the CHS4 promoter region. The dominant / allele (yellow seed coats) exhibits an extra 12.1-kb Hindlll fragment that hybridizes with both the CHS coding region and CHS1 promoter-specific probes. Mutations of the dominant / allele to the recessive i allele (pigmented seed coats) give rise to 10.4- or 9.6-kb Hindlll CHS fragments that have lost the duplicated CHS1 promoter. Finally, gene expression analysis demonstrated that heterozygous plants (I/i) with yellow seed coats have reduced mRNA levels, indicating that the 12.1-kb Hindlll CHS fragment associated with the dominant / allele inhibits pigmentation in a trans-dominant manner. Moreover, CHS gene-specific expression in seed coats shows that multiple CHS genes are expressed in seed coats.     

Genetic analysis of differentiation among breeding ponds reveals a candidate gene for local adaptation in Rana arvalis

One of the main questions in evolutionary and conservation biology is how geographical and environmental features of the landscape shape neutral and adaptive genetic variation in natural populations. The identification of genomic polymorphisms that account for adaptive variation can aid in finding candidate loci for local adaptation. Consequently, a comparison of spatial patterns in neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection at the landscape scale. Many amphibians breed in wetlands, which differ in environmental conditions and in the degree of isolation, enhancing the potential for local adaptation. We used microsatellite markers to measure genetic differentiation among 17 local populations of Rana arvalis breeding in a network of wetlands. We found that locus RC08604 deviated from neutral expectations, suggesting that it is a good candidate for directional selection. We used a genetic network analysis to show that the allele distribution in this locus is correlated with habitat characteristics, whereas this was not the case at neutral markers that displayed a different allele distribution and population network in the study area. The graph approach illustrated the genomic heterogeneity (neutral loci vs. the candidate locus for directional selection) of gene exchange and genetic divergence among populations under directional selection. Limited gene flow between wetlands was only observed at the candidate genomic region under directional selection. RC08604 is partially located inside an up‐regulated thyroid‐hormone receptor (TRβ) gene coordinating the expression of other genes during metamorphosis and appears to be linked with variation in larval life‐history traits found among R. arvalis populations. We suggest that directional selection on genes coding larval life‐history traits is strong enough to maintain the divergence in these genomic regions, reducing the effective recombination of locally adapted alleles but not in other regions of the genome. Integrating this knowledge into conservation plans at the landscape scale will improve the design of management strategies to preserve adaptive genetic diversity in wetland networks.     

Pulsatile flow in a constricted channel.

A nonuniform channel is used as a simple model of a constricted arterial vessel. Flow patterns have been calculated for pulsatile flow with both sinusoidal and nonsinusoidal flow rates for a range of Reynolds number, Re, and Strouhal number, St. The results show that even for relatively low frequency flows a strong vortex wave will be generated with a complex wall shear stress distribution and peak values much greater than those found in steady or unsteady parallel flow. The vortex wave increases in strength with increasing Re and St, with its total length and wavelength independent of Re but inversely proportional to St. The form of the imposed flow rate is found to have an important effect on the flow and the shear stress distribution.     

Geographic and climatic factors associated with the spatial structure of gene frequencies inCastanea sativa Mill. from Turkey

Palynological evidence points to the north-eastern region of present day Turkey as the most probable area of initial domestication of chestnut (Castanea sativa Mill.), after the Würm glaciation. Genetic data on chestnut are consistent with the above hypothesis. We have attempted a finer analysis of the spatial structure of isozyme polymorphisms in 13 populations from eastern and western Turkey. The main objectives were: (i) to determine if the observed spatial patterns are consistent with a patch-like distribution of genetic isolates, or with a mechanism involving long-range gene flow with admixture; (ii) to ascertain if there is some spatial pattern related to possible selective pressures; (iii) if so, to investigate if such pressures are likely to be due to some climatic factors. Our results suggest that long-range gene flow seems to be the major mechanism involved in genetic differentiation in Turkish chestnut. On the other hand, some alleles show a peculiar spatial pattern which could be associated with macro-climatic variables, in particular with the amount of rainfall.     

Effects of parasite and historic driven selection on the diversity and structure of a MHC-II gene in a small mammal species (<Emphasis Type="Italic">Peromyscus leucopus</Emphasis>) undergoing range expansion

Genetic diversity may decrease from the centre to the margin of a species distribution range due to neutral stochastic processes. Selection may also alter genetic diversity in non-neutral markers, such as genes associated with the immune system. Both neutral processes and selection on the immune system are thus expected to affect the spatial distribution of such markers, but the relative strength of each has been scarcely studied. Here, we compared the diversity of a neutral marker (mitochondrial cytochrome b) and a selected marker (DRB gene from the MHC-II), in eastern-North American populations of white-footed mice (Peromyscus leucopus), a species known for its role of main reservoir of the Lyme disease. We observed distinct phylogeographic patterns with these two markers, which may be the result of selection pressure acting upon the DRB gene. As predicted by the central marginal hypothesis, we observed a loss of neutral genetic diversity toward the margin of the species distribution. A decrease in diversity was also observed for the DRB gene, likely due to genetic drift and positive selection operated by helminth parasites. Such a loss in genetic diversity at the range margin may slow down the ongoing expansion of P. leucopus, by counterbalancing the effect of global warming on the mouse survival at higher latitude.     

Modes of speciation and the neutral theory of biodiversity

Hubbell's neutral theory of biodiversity has generated much debate over the need for niches to explain biodiversity patterns. Discussion of the theory has focused on its neutrality assumption, i.e. the functional equivalence of species in competition and dispersal. Almost no attention has been paid to another critical aspect of the theory, the assumptions on the nature of the speciation process. In the standard version of the neutral theory each individual has a fixed probability to speciate. Hence, the speciation rate of a species is directly proportional to its abundance in the metacommunity. We argue that this assumption is not realistic for most speciation modes because speciation is an emergent property of complex processes at larger spatial and temporal scales and, consequently, speciation rate can either increase or decrease with abundance. Accordingly, the assumption that speciation rate is independent of abundance (each species has a fixed probability to speciate) is a more natural starting point in a neutral theory of biodiversity. Here we present a neutral model based on this assumption and we confront this new model to 20 large data sets of tree communities, expecting the new model to fit the data better than Hubbell's original model. We find, however, that the data sets are much better fitted by Hubbell's original model. This implies that species abundance data can discriminate between different modes of speciation, or, stated otherwise, that the mode of speciation has a large impact on the species abundance distribution. Our model analysis points out new ways to study how biodiversity patterns are shaped by the interplay between evolutionary processes (speciation, extinction) and ecological processes (competition, dispersal).     

Neutral mutation as the source of genetic variation in life history traits

The mechanism underlying the maintenance of adaptive genetic variation is a long-standing question in evolutionary genetics. There are two concepts (mutation-selection balance and balancing selection) which are based on the phenotypic differences between alleles. Mutation - selection balance and balancing selection cannot properly explain the process of gene substitution, i.e. the molecular evolution of quantitative trait loci affecting fitness. I assume that such loci have non-essential functions (small effects on fitness), and that they have the potential to evolve into new functions and acquire new adaptations. Here I show that a high amount of neutral polymorphism at these loci can exist in real populations. Consistent with this, I propose a hypothesis for the maintenance of genetic variation in life history traits which can be efficient for the fixation of alleles with very small selective advantage. The hypothesis is based on neutral polymorphism at quantitative trait loci and both neutral and adaptive gene substitutions. The model of neutral - adaptive conversion (NAC) assumes that neutral alleles are not neutral indefinitely, and that in specific and very rare situations phenotypic (relative fitness) differences between them can appear. In this paper I focus on NAC due to phenotypic plasticity of neutral alleles. The important evolutionary consequence of NAC could be the increased adaptive potential of a population. Loci responsible for adaptation should be fast evolving genes with minimally discernible phenotypic effects, and the recent discovery of genes with such characteristics implicates them as suitable candidates for loci involved in adaptation.     

The Age of a Neutral Mutant Persisting in a Finite Population

Formulae for the mean and the mean square age of a neutral allele which is segregating with frequency x in a population of effective size N(e) have been obtained using the diffusion equation method, for the case of 4N(e)v<1 where v is the mutation rate. It has been shown that the average ages of neutral alleles, even if their frequencies are relatively low, are quite old. For example, a neutral mutant whose current frequency is 10% has the expected age roughly equal to the effective population size N(e) and the standard deviation 1.4N(e) (in generations), assuming that this mutant has increased by random drift from a very low frequency. Also, formulae for the mean "first arrival time" of a neutral mutant to a certain frequency x have been presented. In addition, a new, approximate method has been developed which enables us to obtain the condition under which frequencies of "rare" polymorphic alleles among local populations are expected to be uniform if the alleles are selectively neutral.-It was concluded that exchange of only a few individuals on the average between adjacent colonies per generation is enough to bring about such a uniformity of frequencies.     

High genetic diversity in Sarracenia leucophylla(Sarraceniaceae), a carnivorous wetland herb

Eighteen allozyme loci were used to examine genetic diversity in 10 natural populations of Sarracenia leucophylla Raf., a pitcher plant restricted to the southeastern United States. One ex situ population propagated for restoration in Georgia was also analyzed. S. leucophylla is an insect-pollinated, outcrossing perennial wetland herb that is threatened over much of its geographic range. Fifteen loci (83.3%) were polymorphic, with a mean number of alleles of 3.33. Compared to species having similar life-history traits and to previously analyzed Sarracenia species, S. leucophylla displayed unexpectedly high genetic diversity. For example, genetic diversity within the species (Hes) was 0.224 and mean population genetic diversity (Hep) was 0.183. Although small S. leucophylla populations maintained less genetic diversity than larger ones, these differences were not statistically significant. Nonetheless, this suggests that small populations may have lost rare alleles. Statistically significant genetic differentiation among populations was found (theta = 0.192, P < .01), although it was not atypical considering the species' life-history characteristics. A significant correlation (P < .01) between genetic and geographic distance was found, indicating an isolation-by-distance effect. However, the correlation coefficient for this relationship was low (r = 0.46), suggesting that factors other than gene flow play a prominent role in the geographic distribution of genetic diversity within the species. The ex situ population captured most of the allozyme variation found in its source population.     

Testing the generalized neutrality hypothesis

Various tests of the hypothesis of selective neutrality based on gene frequency are now available. These tests take as null hypothesis the concept of “strict neutrality”: all new mutants are required to be selectively identical to each other. For evolutionary questions, however, (as opposed to those of genetic polymorphism), a wider null hypothesis might be of interest. Since deleterious alleles have essentially no evolutionary importance, one might wish to test the null hypothesis that only neutral or deleterious mutations occur. The principal alternative to this hypothesis is that there exists heterotic selection of some form for some alleles tending to maintain a level of genetic polymorphism higher than that under neutrality. In this paper an assessment is made of the usefulness of a test of strict neutrality first proposed by this author (Ewens, 1972) as a test of null hypothesis of “generalized neutrality,” i.e. that only neutral or deleterious alleles occur. At the same time some remarks will be made about estimation of the fundamental parameter θ defining these processes.     

Human disorders in N-glycosylation and animal models

Genes that cause human disorders in N-linked oligosaccharide biosynthesis have appeared much faster than animal model systems to study them. In most models, a single gene is altered or deleted while other genes and the environment are held constant. Since humans have variable genetic backgrounds and environments, model systems may only partially mimic the actual disorders. Mutations in seven of the 30-40 genes needed for the synthesis and transfer of oligosaccharides from the lipid donor to the nascent protein acceptors in the endoplasmic reticulum cause Type I Congenital Disorders of Glycosylation (CDG). Since all of these gene products ultimately contribute to the same final step, one might suspect that all the diseases would be very similar. However, even patients with mutations in the same gene show considerable phenotypic variability. Modifier, or susceptibility genes in the background likely explain some variations of the "primary" gene chosen for study. Add to this the stress of infections, dietary insufficiencies, and the demands of growth itself. These issues are particularly important during development when the temporal and spatial specific interplay of cell adhesions and signals has only a single opportunity. Multiple hypomorphic alleles of genes in the same pathway may have synergistic effects. Investigators designing model systems to study human glycosylation disorders may want to construct strains with several heterozygous hypomorphic alleles in rate-limiting steps in the glycosylation pathway.     

Selection and Biased Gene Conversion in a Multigene Family: Consequences of Interallelic Bias and Threshold Selection

In a previous paper, I investigated the interactions in a gene family of additive selection and biased gene conversion in a finite population when conversion events are rare. Here I extend my "weak-conversion limit" model by allowing biased interallelic conversion (conversion between alleles at the same locus) of arbitrary frequency and various threshold selection schemes for rare interlocus conversion events. I suggest that it is not unreasonable for gene families to experience threshold fitness functions, and show that certain types of thresholds can greatly constrain the rate at which advantageous alleles are fixed as compared to other fitness schemes, such as additive selection. It is also shown that the double sampling process operating on a gene family in a finite population (sampling over the number of genes in the gene family and over the number of individuals in the population) can have interesting consequences. For selectively neutral alleles that experience interallelic bias, the probability of fixation at each single locus may be essentially neutral, but the cumulative effects on the entire gene family of small departures from neutrality can be significant, especially if the gene family is large. Thus, in some situations, gene families can respond to directional forces that are weak in comparison to drift at single loci.     

The influence of altitude and topography on genetic structure in the long-toed salamander (Ambystoma macrodactulym)

A primary goal of molecular ecology is to understand the influence of abiotic factors on the spatial distribution of genetic variation. Features including altitudinal clines, topography and landscape characteristics affect the proportion of suitable habitat, influence dispersal patterns, and ultimately structure genetic differentiation among populations. We studied the effects of altitude and topography on genetic variation of long-toed salamanders (Ambystoma macrodactylum), a geographically widespread amphibian species throughout northwestern North America. We focused on 10 low altitude sites (< 1200 m) and 11 high-altitude sites in northwestern Montana and determined multilocus genotypes for 549 individuals using seven microsatellite loci. We tested four hypotheses: (1) gene flow is limited between high- and low-altitude sites; and, (2) gene flow is limited among high-altitude sites due to harsh habitat and extreme topographical relief between sites; (3) low-altitude sites exhibit higher among-site gene flow due to frequent flooding events and low altitudinal relief; and (4) there is a negative correlation between altitude and genetic variation. Overall F(ST) values were moderate (0.08611; P < 0.001). Pairwise F(ST) estimates between high and low populations and a population graphing method supported the hypothesis that low-altitude and high-altitude sites, taken together, are genetically differentiated from each other. Also as predicted, gene flow is more prominent among low-altitude sites than high-altitude sites; low-altitude sites had a significantly lower F(ST) (0.03995; P < 0.001) than high altitude sites (F(ST) = 0.10271; P < 0.001). Use of Bayesian analysis of population structure (BAPS) resulted in delineation of 10 genetic groups, two among low-altitude populations and eight among high-altitude populations. In addition, within high altitude populations, basin-level genetic structuring was apparent. A nonequilibrium algorithm for detecting current migration rates supported these population distinctions. Finally, we also found a significant negative correlation between genetic diversity and altitude. These results are consistent with the hypothesis that topography and altitudinal gradients shape the spatial distribution of genetic variation in a species with a broad geographical range and diverse life history. Our study sheds light on which key factors limit dispersal and ultimately species' distributions.