Integrating multiple analytical approaches to spatially delineate and characterize genetic population structure: an application to boreal caribou (<Emphasis Type="Italic">Rangifer tarandus caribou</Emphasis>) in central Canada

Individual-based clustering (IBC) methods have become increasingly popular for the characterization and delineation of genetic population units for numerous species. These methods delineate populations based on the genetic assumptions of a breeding unit which may provide a better representation of the behaviour of the species. The increasing use of IBC has resulted in the development of several analytical models all of which vary in their theoretical assumptions to infer genetic population structure. In this paper, we report a comparative strategy utilizing three IBC methods to characterize the spatial genetic structure of the boreal population of woodland caribou (Rangifer tarandus caribou) in central Canada. In addition, we implement both tests for isolation-by-distance (IBD) and frequency-based assignment tests to validate the consensus genetic clusters as defined by IBC. We also compare indirect metrics of genetic diversity and gene flow using both a priori defined herds and the IBC defined populations. Although our results show some concordance between both pre-defined herds and IBC derived genetic clusters, the IBC analyses identified a cluster that was cryptic to observation-based caribou herds and found no difference between several adjacent herds. By comparing multiple IBC methods and integrating both IBD and indirect genetic diversity metrics a posteriori, our strategy provides an effective means to delineate wildlife population structure and accurately assess genetic diversity and connectivity.     

Partitioning drivers of spatial genetic variation for a continuously distributed population of boreal caribou: Implications for management unit delineation

Isolation by distance (IBD) is a natural pattern not readily incorporated into theoretical models nor traditional metrics for differentiating populations, although clinal genetic differentiation can be characteristic of many wildlife species. Landscape features can also drive population structure additive to baseline IBD resulting in differentiation through isolation‐by‐resistance (IBR). We assessed the population genetic structure of boreal caribou across western Canada using nonspatial (STRUCTURE) and spatial (MEMGENE) clustering methods and investigated the relative contribution of IBD and IBR on genetic variation of 1,221 boreal caribou multilocus genotypes across western Canada. We further introduced a novel approach to compare the partitioning of individuals into management units (MU) and assessed levels of genetic connectivity under different MU scenarios. STRUCTURE delineated five genetic clusters while MEMGENE identified finer‐scale differentiation across the study area. IBD was significant and did not differ for males and females both across and among detected genetic clusters. MEMGENE landscape analysis further quantified the proportion of genetic variation contributed by IBD and IBR patterns, allowing for the relative importance of spatial drivers, including roads, water bodies, and wildfires, to be assessed and incorporated into the characterization of population structure for the delineation of MUs. Local population units, as currently delineated in the boreal caribou recovery strategy, do not capture the genetic variation and connectivity of the ecotype across the study area. Here, we provide the tools to assess fine‐scale spatial patterns of genetic variation, partition drivers of genetic variation, and evaluate the best management options for maintaining genetic connectivity. Our approach is highly relevant to vagile wildlife species that are of management and conservation concern and demonstrate varying degrees of IBD and IBR with clinal spatial genetic structure that challenges the delineation of discrete population boundaries.     

Reconstruction of caribou evolutionary history in Western North America and its implications for conservation

The role of Beringia as a refugium and route for trans-continental exchange of fauna during glacial cycles of the past 2million years are well documented; less apparent is its contribution as a significant reservoir of genetic diversity. Using mitochondrial DNA sequences and 14 microsatellite loci, we investigate the phylogeographic history of caribou (Rangifer tarandus) in western North America. Patterns of genetic diversity reveal two distinct groups of caribou. Caribou classified as a Northern group, of Beringian origin, exhibited greater number and variability in mtDNA haplotypes compared to a Southern group originating from refugia south of glacial ice. Results indicate that subspecies R. t. granti of Alaska and R. t. groenlandicus of northern Canada do not constitute distinguishable units at mtDNA or microsatellites, belying their current status as separate subspecies. Additionally, the Northern Mountain ecotype of woodland caribou (presently R. t. caribou) has closer kinship to caribou classified as granti or groenlandicus. Comparisons of mtDNA and microsatellite data suggest that behavioural and ecological specialization is a more recently derived life history characteristic. Notably, microsatellite differentiation among Southern herds is significantly greater, most likely as a result of human-induced landscape fragmentation and genetic drift due to smaller population sizes. These results not only provide important insight into the evolutionary history of northern species such as caribou, but also are important indicators for managers evaluating conservation measures for this threatened species.     

Population structure over a broad spatial scale driven by nonanthropogenic factors in a wide‐ranging migratory mammal,Alaskan caribou

Wide‐ranging mammals face significant conservation threats, and knowledge of the spatial scale of population structure and its drivers is needed to understand processes that maintain diversity in these species. We analysed DNA from 655 Alaskan caribou (Rangifer tarandus granti) from 20 herds that vary in population size, used 19 microsatellite loci to document genetic diversity and differentiation in Alaskan caribou, and examined the extent to which genetic differentiation was associated with hypothesized drivers of population subdivision including landscape features, population size and ecotype. We found that Alaskan caribou are subdivided into two hierarchically structured clusters: one group on the Alaska Peninsula containing discrete herds and one large group on the Mainland lacking differentiation between many herds. Population size, geographic distance, migratory ecotype and the Kvichak River at the nexus of the Alaska Peninsula were associated with genetic differentiation. Contrary to previous hypotheses, small Mainland herds were often differentiated genetically from large interconnected herds nearby, and genetic drift coupled with reduced gene flow may explain this pattern. Our results raise the possibility that behaviour helps to maintain genetic differentiation between some herds of different ecotypes. Alaskan caribou show remarkably high diversity and low differentiation over a broad geographic scale. These results increase information for the conservation of caribou and other migratory mammals threatened by population reductions and landscape barriers and may be broadly applicable to understanding the spatial scale and ecological drivers of population structure in widespread species.     

Population structure of caribou in an ice‐bound archipelago

Aim

Archipelagos provide ideal natural systems for testing the effects of isolation and fragmentation of habitats on the genetic makeup of populations—an important consideration, given that many insular species are of conservation concern. Two theories predominate: Island Biogeography Theory (IBT) posits that proximity to the mainland drives the potential for migrants and gene flow. The Central Marginal Hypothesis (CMH) predicts that island populations at the periphery of a species range may experience low gene flow, small population size and high rates of genetic drift. We investigated population genetic structure, genetic diversity and key drivers of diversity for Arctic island‐dwelling caribou (Rangifer tarandus). Our aim was to inform intraspecific units for conservation and decipher how IBT and CMH could act in an archipelago where isolation is highly variable due to sea ice and open water.

Location

Canadian Arctic Archipelago, Canada (Latitude, 55–82°N; Longitude, 61–123°W).

Methods

We genotyped 447 caribou at 16 microsatellite loci; these caribou represented two subspecies (R. t. groenlandicus, R. t. pearyi) and three designatable units. We used hierarchical Bayesian clustering and ordination to determine genetic groups. We evaluated the influence of ecological and geographic variables on genetic diversity using linear mixed‐effects models and compared diversity among mainland and island herds.

Results

Bayesian clustering revealed nine genetic clusters with differentiation among and within caribou subspecies. Genetic differentiation was explained predominantly by isolation‐by‐distance across all caribou, even at the scale of subspecies. Island caribou were less genetically diverse than mainland herds; individual heterozygosity was negatively correlated with distance‐to‐mainland and the extent of autumn ice‐free coastline and positively correlated with unglaciated island size.

Main conclusions

Our findings underscore the importance of hierarchical analysis when investigating genetic population structure. Genetic diversity and its key drivers lend support to both IBT and CMH and highlight the pending threat of climate change for Arctic island caribou.

     

Conservation implications of significant population differentiation in an endangered estuarine seahorse

The spatial distribution of a species’ genetic diversity can provide insights into underlying evolutionary, ecological and environmental processes, and can contribute information towards the delineation of conservation units. The Knysna seahorse, Hippocampus capensis, is endangered and occurs in only three estuaries on the warm-temperate south coast of South Africa: Knsyna, Keurbooms and Swartvlei. Population sizes in the latter two estuaries have been very small for a prolonged period of time, and the populations residing in them may thus benefit from translocations as a means of increasing population sizes and possibly also genetic diversity. However, information on whether these three estuaries constitute distinct conservation units that warrant separate management is presently lacking. Here, we used genetic information from mitochondrial (control region) and nuclear microsatellite loci to assess the genetic diversity and spatial structure across the three estuaries, and also whether translocations should be included in the management plan for the Knysna seahorse. Although each population had a unique combination of alleles, and clustering methods identified the Swartvlei Estuary as being distinct from the others, levels of genetic admixture were high, and there was no evidence for reciprocal monophyly that would indicate that each estuary has a unique demographic history. On these grounds, we suggest recognising the three populations as a single evolutionarily significant unit (ESU), and encourage translocations between them to ensure the species’ long-term survival.     

Conservation genetics of the yellow-bellied toad (Bombina variegata): population structure,genetic diversity and landscape effects in an endangered amphibian

Revealing patterns of genetic diversity and barriers for gene flow are key points for successful conservation in endangered species. Methods based on molecular markers are also often used to delineate conservation units such as evolutionary significant units and management units. Here we combine phylo-geographic analyses (based on mtDNA) with population and landscape genetic analyses (based on microsatellites) for the endangered yellow-bellied toad Bombina variegata over a wide distribution range in Germany. Our analyses show that two genetic clusters are present in the study area, a northern and a southern/central one, but that these clusters are not deeply divergent. The genetic data suggest high fragmentation among toad occurrences and consequently low genetic diversity. Genetic diversity and genetic connectivity showed a negative relationship with road densities and urban areas surrounding toad occurrences, indicating that these landscape features act as barriers to gene flow. To preserve a maximum of genetic diversity, we recommend considering both genetic clusters as management units, and to increase gene flow among toad occurrences with the aim of restoring and protecting functional meta-populations within each of the clusters. Several isolated populations with especially low genetic diversity and signs of inbreeding need particular short-term conservation attention to avoid extinction. We also recommend to allow natural gene flow between both clusters but not to use individuals from one cluster for translocation or reintroduction into the other. Our results underscore the utility of molecular tools for species conservation, highlight outcomes of habitat fragmentation onto the genetic structure of an endangered amphibian and reveal particularly threatened populations in need for urgent conservation efforts.

     

Prioritization of landscape connectivity for the conservation of Peary caribou

Adequate connectivity between discontinuous habitat patches is crucial for the persistence of metapopulations across space and time. Loss of landscape connectivity is often a direct result of fragmentation caused by human activities but also can be caused indirectly through anthropogenic climate change. Peary caribou (Rangifer tarandus pearyi) are widely dispersed across the islands of the Canadian Arctic Archipelago and rely on sea ice to move seasonally between island habitats throughout their range. Seasonal connectivity provided by sea ice is necessary to maintain genetic diversity and to facilitate dispersal and recolonization of areas from which caribou have been extirpated. We used least‐cost path analysis and circuit theory to model connectivity across Peary caribou range, and future climate projections to investigate how this connectivity might be affected by a warming climate. Further, we used measures of current flow centrality to estimate the role of High Arctic islands in maintaining connectivity between Peary caribou populations and to identify and prioritize those islands and linkages most important for conservation. Our results suggest that the Bathurst Island complex plays a critical role in facilitating connectivity between Peary caribou populations. Large islands, including Banks, Victoria, and Ellesmere have limited roles in connecting Peary caribou. Without rigorous greenhouse gas emission reductions our projections indicate that by 2100 all connectivity between the more southern Peary caribou populations will be lost for important spring and early‐winter movement periods. Continued connectivity across the Canadian Arctic Archipelago, and possibly Peary caribou persistence, ultimately hinges on global commitments to limit climate change. Our research highlights priority areas where, in addition to emission reductions, conservation efforts to maintain connectivity would be most effective.     

Spatial genetic structure of <Emphasis Type="Italic">Lissotriton helveticus L</Emphasis>. following the restoration of a forest ponds network

Preserving amphibian genetic diversity through ecological restoration and conservation actions is a major challenge since their populations are declining worldwide. We studied the genetic diversity and spatial genetic structure of the palmate newt (Lissotriton helveticus) 2 years after the restoration of a pond network in northwestern France with the aim of reconstructing fine-scale genetic structure and patterns of colonization. We sampled newts from 29 forest ponds including both restored and non-degraded reference ponds, and genotyped 391 individuals at 12 microsatellite loci. We used two Bayesian clustering methods to spatially delineate genetic clusters, and we also detected potential recent migrants within the network. All ponds showed low levels of observed heterozygosity (Ho?=?0.534) and a mean F _IS of 0.251, possibly indicating a Wahlund or bottleneck effect. Pairwise F _ST suggested limited evidence of genetic differentiation among ponds. Within the pond network, we identified 3 to 4 genetic clusters. Combined with the detection of migrants, the results suggest an increase in gene flow within the restored pond network and that a high number of migrants came from the reference ponds. Our findings indicate an unexpected high dispersal ability for this small-bodied species. Overall, the absence of population structure represents a positive beginning for the restoration project. It also emphasizes the importance of spatial design in restoring a pond network and that such genetic data and methods should be used to monitor amphibians in restored habitats.     

Significance of caribou (Rangifer tarandus) ecotypes from a molecular genetics viewpoint

Three caribou ecotypes are present in easternNorth America: the mountain caribou which isfound south of the St. Lawrence River, thebarren-ground caribou which calves in thetundra, and in between, the forest-dwellingecotype which lives all year long in the borealforest. Blood and muscle samples were collectedfrom seven populations and characterized ateight microsatellite loci to test thehypotheses that forest-dwelling andbarren-ground ecotypes constitute a singlemetapopulation and that geographical isolationresults in reduced genetic diversity. The meannumber of alleles per locus, allelic richness,and observed and expected heterozygositydeclined from north to south and were thesmallest in isolated forest-dwellingpopulations. Correspondence analysis showedthree groups of samples corresponding to thethree ecotypes. Gene flow estimates weremoderate or high among all forest-dwellingpopulations and particularly between those <200 km apart. Our results suggest that thethree caribou ecotypes represent three distinctgenetic entities and that the forest-dwellingpopulations in the continuous range form ametapopulation. Genetic diversity was lower inisolated populations but does not seem to be ofimmediate concern for conservation. We proposethat management strategies should favorincrease in caribou numbers in order to avoidextinction due to stochastic events and tomaintain local biodiversity. In the continuousrange, conservation strategies of cariboupopulations must be planned on a large scale tomaintain occasional exchanges betweenpopulations, thus preserving genetic diversity.     

Cryptic structure and niche divergence within threatened Galápagos giant tortoises from southern Isabela Island

Although Galápagos giant tortoises are an icon for both human-mediated biodiversity losses and conservation management successes, populations of two species on southern Isabela Island (Chelonoidis guntheri, and C. vicina) remain threatened by hunting and persistence of feral animals. Conservation management of these tortoises has been hampered by lack of clarity regarding their taxonomy, ecological and morphological diversity, and the spatial distribution of evolutionarily significant units that may exist. Analyses of 16 microsatellite loci did not group samples according to current taxonomy. Instead, three (rather than two) genetic clusters were revealed. We show that the three regions of southern Isabela associated with these genetic clusters are significantly different in their ecological niches, which could suggest that ecological divergence may have shaped patterns of genetic differentiation in these tortoises. Furthermore, results suggest limited recent gene flow among sampled localities and between each of the three regions associated with genetic clusters. We discuss the need for further research on the ecological factors shaping the genetic and morphological diversity of southern Isabela tortoises. We suggest that current strategies whereby populations are managed separately are warranted pending further study, but due to mixed ancestry we recommend that Cerro Paloma tortoises be excluded from management programs.     

Spatial differences in genetic diversity and northward migration suggest genetic erosion along the boreal caribou southern range limit and continued range retraction

With increasing human activities and associated landscape changes, distributions of terrestrial mammals become fragmented. These changes in distribution are often associated with reduced population sizes and loss of genetic connectivity and diversity (i.e., genetic erosion) which may further diminish a species' ability to respond to changing environmental conditions and lead to local population extinctions. We studied threatened boreal caribou (Rangifer tarandus caribou) populations across their distribution in Ontario/Manitoba (Canada) to assess changes in genetic diversity and connectivity in areas of high and low anthropogenic activity. Using data from >1,000 caribou and nine microsatellite loci, we assessed population genetic structure, genetic diversity, and recent migration rates using a combination of network and population genetic analyses. We used Bayesian clustering analyses to identify population genetic structure and explored spatial and temporal variation in those patterns by assembling networks based on R_ST and F_ST as historical and contemporary genetic edge distances, respectively. The Bayesian clustering analyses identified broad‐scale patterns of genetic structure and closely aligned with the R_ST network. The F_ST network revealed substantial contemporary genetic differentiation, particularly in areas presenting contemporary anthropogenic disturbances and habitat fragmentation. In general, relatively lower genetic diversity and greater genetic differentiation were detected along the southern range limit, differing from areas in the northern parts of the distribution. Moreover, estimation of migration rates suggested a northward movement of animals away from the southern range limit. The patterns of genetic erosion revealed in our study suggest ongoing range retraction of boreal caribou in central Canada.     

Environmental and spatial characterisation of an unknown fauna using DNA sequencing – an example with Himalayan Hydropsychidae (Insecta: Trichoptera)

相似文献   

Population size and major valleys explain microsatellite variation better than taxonomic units for caribou in western Canada

Identifying conservation units below the species level is becoming increasingly important, particularly when limited resources necessitate prioritization for conservation among such units. This problem is exemplified with caribou, a mammal with a circum-Arctic distribution that is exposed to a broad spectrum of ecological conditions, but is also declining in many parts of its range. We used microsatellite markers to evaluate the suitability of existing intra-specific taxonomic designations to act as population units for conservation and contrasted this with landscape features that were independent of taxonomy. We also quantified the relationship between genetic differentiation and subpopulation size, a factor that has been under-represented in landscape genetic research. Our data set included three subspecies and three ecotypes of caribou that varied in population size by five orders of magnitude. Our results indicated that genetic structure did not correspond to existing taxonomic designation, particularly at the level of ecotype. Instead, we found that major valleys and population size were the strongest factors associated with substructure. There was a negative exponential relationship between population size and F(ST) between pairs of adjacent subpopulations, suggesting that genetic drift was the mechanism causing the structure among the smallest subpopulations. A genetic assignment test revealed that movement among subpopulations was a fraction of the level needed to stabilize smaller subpopulations, indicating little chance for demographic rescue. Such results may be broadly applicable to landscape genetic studies, because population size and corresponding rates of drift have the potential to confound interpretations of landscape effects on population structure.     

Linking genetic and ecological differentiation in an ungulate with a circumpolar distribution

Genetic differentiation among populations may arise from the disruption of gene flow due to local adaptation to distinct environments and/or neutral accumulation of mutations and genetic drift resulted from geographical isolation. Quantifying the role of these processes in determining the genetic structure of natural populations remains challenging. Here, we analyze the relative contribution of isolation‐by‐resistance (IBR), isolation‐by‐environment (IBE), genetic drift and historical isolation in allopatry during Pleistocene glacial cycles on shaping patterns of genetic differentiation in caribou/reindeer populations Rangifer tarandus across the entire distribution range of the species. Our study integrates analyses at range‐wide and regional scales to partial out the effects of historical and contemporary isolation mechanisms. At the circumpolar scale, our results indicate that genetic differentiation is predominantly explained by IBR and historical isolation. At a regional scale, we found that IBR, IBE and population size significantly explained the spatial distribution of genetic variation among populations belonging to the Euro‐Beringian lineage within North America. In contrast, genetic differentiation among populations within the North American lineage was predominantly explained by IBR and population size, but not IBE. We also found discrepancies between genetic and ecotype designation across the Holarctic species distribution range. Overall, these results indicate that multiple isolating mechanisms have played roles in shaping the spatial distribution of genetic variation across the distribution range of a large mammal with high potential for gene flow. Considering multiple spatial scales and simultaneously testing a comprehensive suite of potential isolating mechanisms, our study contributes to understand the ecological and evolutionary processes underlying organism–landscape interactions.     

Extreme levels of hidden diversity in microscopic animals (Rotifera) revealed by DNA taxonomy

Knowledge and estimates of species richness at all scales are biased both by our understanding of the evolutionary processes shaping diversity and by the methods used to delineate the basic units of diversity. DNA taxonomy shows that diversity may be underestimated by traditional taxonomy, especially for microscopic animals. The effects of such hidden diversity are usually overlooked in ecological studies. Here, we estimate hidden diversity in bdelloid rotifers, a group of microscopic animals. We analyse cryptic diversity using a coalescent approach to infer taxonomical units from phylogenetic trees. Cryptic diversity was measured for eight traditional species of bdelloid rotifers and the results compared to that of the monogonont rotifer Brachionus plicatilis species complex, which is well studied and for which cross-mating experiments have been performed to explicitly define some of the species boundaries. A taxonomic inflation of up to 34 potential cryptic taxa was found in bdelloids. Cryptic taxa within each traditional species may be spatially isolated, but do not have narrower ecological niches. The species deemed as generalists exhibit the highest cryptic diversity. Cryptic diversity based on molecular characterization is commonly found in animals; nevertheless, the amount of cryptic diversity in bdelloids is much higher than in other groups analysed so far, maybe because of their peculiar parthenogenetic reproduction, other than microscopic size. We discuss this hypothesis in the light of the available empirical evidence from other groups of microscopic invertebrates, such as tardigrades and mites, which share size, habitat heterogeneity, potential for dispersal, and/or parthenogenetic reproduction.     

Implications of recurrent disturbance for genetic diversity

Exploring interactions between ecological disturbance, species’ abundances and community composition provides critical insights for ecological dynamics. While disturbance is also potentially an important driver of landscape genetic patterns, the mechanisms by which these patterns may arise by selective and neutral processes are not well‐understood. We used simulation to evaluate the relative importance of disturbance regime components, and their interaction with demographic and dispersal processes, on the distribution of genetic diversity across landscapes. We investigated genetic impacts of variation in key components of disturbance regimes and spatial patterns that are likely to respond to climate change and land management, including disturbance size, frequency, and severity. The influence of disturbance was mediated by dispersal distance and, to a limited extent, by birth rate. Nevertheless, all three disturbance regime components strongly influenced spatial and temporal patterns of genetic diversity within subpopulations, and were associated with changes in genetic structure. Furthermore, disturbance‐induced changes in temporal population dynamics and the spatial distribution of populations across the landscape resulted in disrupted isolation by distance patterns among populations. Our results show that forecast changes in disturbance regimes have the potential to cause major changes to the distribution of genetic diversity within and among populations. We highlight likely scenarios under which future changes to disturbance size, severity, or frequency will have the strongest impacts on population genetic patterns. In addition, our results have implications for the inference of biological processes from genetic data, because the effects of dispersal on genetic patterns were strongly mediated by disturbance regimes.     

Reconceptualising the beta diversity‐environmental heterogeneity relationship in running water systems

相似文献   

Genetic diagnosis of a rare myrmecochorous species,Plagiorhegma dubium (Berberidaceae): Historical genetic bottlenecks and strong spatial structures among populations

Distribution of genetic variation over time and space is relevant to demographic histories and tightly linked to ecological disturbances as well as evolutionary potential of an organism. Therefore, understanding the pattern of genetic diversity is a primary step in conservation and management projects for rare and threatened plant species. We used eight microsatellite markers to examine the level of genetic diversity, spatial structure, and demographic history of Plagiorhegma dubium, a rare myrmecochorous herb, populations sampled across northeast Asia and Siberia. We found low within‐population genetic variation associated with historical bottlenecks. Although pairwise F_ST values were not much higher than the ones found in similar life form species, STRUCTURE and PCoA revealed a clear broadscale spatial pattern of genetic structure. Bayesian clustering (best K = 6) and PCoA identified three populations that are distinctive from neighboring populations in the Korean peninsula, which suggests potential units for conservation and management plans in Korea. MIGRATE‐N and BAYESASS showed that both contemporary (0.003–0.045) and historical migration rates (2 × e^?5?4.6 × e^?4) were low. Our findings provide a good example, where genetic considerations should be integrated for conservation and management plans of rare and threatened species.     

Impact of Geography and Climate on the Genetic Differentiation of the Subtropical Pine Pinus yunnanensis

Southwest China is a biodiversity hotspot characterized by complex topography, heterogeneous regional climates and rich flora. The processes and driving factors underlying this hotspot remain to be explicitly tested across taxa to gain a general understanding of the evolution of biodiversity and speciation in the region. In this study, we examined the role played by historically neutral processes, geography and environment in producing the current genetic diversity of the subtropical pine Pinus yunnanensis. We used genetic and ecological methods to investigate the patterns of genetic differentiation and ecological niche divergence across the distribution range of this species. We found both continuous genetic differentiation over the majority of its range, and discrete isolated local clusters. The discrete differentiation between two genetic groups in the west and east peripheries is consistent with niche divergence and geographical isolation of these groups. In the central area of the species’ range, population structure was shaped mainly by neutral processes and geography rather than by ecological selection. These results show that geographical and environmental factors together created stronger and more discrete genetic differentiation than isolation by distance alone, and illustrate the importance of ecological factors in forming or maintaining genetic divergence across a complex landscape. Our findings differ from other phylogenetic studies that identified the historical drainage system in the region as the primary factor shaping population structure, and highlight the heterogeneous contributions that geography and environment have made to genetic diversity among taxa in southwest China.