Genetic spatial structure in a butterfly metapopulation correlates better with past than present demographic structure

The Glanville fritillary butterfly ( Melitaea cinxia ) has been studied in the Åland Islands in Finland since 1991, where it occurs as a classic metapopulation in a large network of 4000 dry meadows. Much ecological work has been conducted on this species, but population genetic studies have been hampered by paucity of suitable genetic markers. Here, using single nucleotide polymorphisms and microsatellites developed for the Glanville fritillary, we examine the correspondence between the demographic and genetic spatial structures. Given the dynamic nature of the metapopulation, the current genetic spatial structure may bear a signal of past changes in population sizes and past patterns of gene flow rather than reflect the current demographic structure or landscape structure. We analyse this question with demographic data for 10 years, using the Rand index to assess the similarity between the genetic, demographic, and landscape spatial structures. Our results show that the current genetic spatial structure is better explained by the past rather than by the current demographic spatial structure or by the spatial configuration of the habitat in the landscape. Furthermore, current genetic diversity is significantly explained by past metapopulation sizes. The time lag between major demographic events and change in the genetic spatial structure and diversity has implications for the study of spatial dynamics.     

Spatial and temporal patterns of neutral and adaptive genetic variation in the endangered African wild dog (Lycaon pictus)

Deciphering patterns of genetic variation within a species is essential for understanding population structure, local adaptation and differences in diversity between populations. Whilst neutrally evolving genetic markers can be used to elucidate demographic processes and genetic structure, they are not subject to selection and therefore are not informative about patterns of adaptive variation. As such, assessments of pertinent adaptive loci, such as the immunity genes of the major histocompatibility complex (MHC), are increasingly being incorporated into genetic studies. In this study, we combined neutral (microsatellite, mtDNA) and adaptive (MHC class II DLA‐DRB1 locus) markers to elucidate the factors influencing patterns of genetic variation in the African wild dog (Lycaon pictus); an endangered canid that has suffered extensive declines in distribution and abundance. Our genetic analyses found all extant wild dog populations to be relatively small (N_e < 30). Furthermore, through coalescent modelling, we detected a genetic signature of a recent and substantial demographic decline, which correlates with human expansion, but contrasts with findings in some other African mammals. We found strong structuring of wild dog populations, indicating the negative influence of extensive habitat fragmentation and loss of gene flow between habitat patches. Across populations, we found that the spatial and temporal structure of microsatellite diversity and MHC diversity were correlated and strongly influenced by demographic stability and population size, indicating the effects of genetic drift in these small populations. Despite this correlation, we detected signatures of selection at the MHC, implying that selection has not been completely overwhelmed by genetic drift.     

Landscape genetics in a changing world: disentangling historical and contemporary influences and inferring change

Landscape genetics seeks to determine the effect of landscape features on gene flow and genetic structure. Often, such analyses are intended to inform conservation and management. However, depending on the many factors that influence the time to reach equilibrium, genetic structure may more strongly represent past rather than contemporary landscapes. This well‐known lag between current demographic processes and population genetic structure often makes it challenging to interpret how contemporary landscapes and anthropogenic activity shape gene flow. Here, we review the theoretical framework for factors that influence time lags, summarize approaches to address this temporal disconnect in landscape genetic studies, and evaluate ways to make inferences about landscape change and its effects on species using genetic data alone or in combination with other data. Those approaches include comparing correlation of genetic structure with historical versus contemporary landscapes, using molecular markers with different rates of evolution, contrasting metrics of genetic structure and gene flow that reflect population genetic processes operating at different temporal scales, comparing historical and contemporary samples, combining genetic data with contemporary estimates of species distribution or movement, and controlling for phylogeographic history. We recommend using simulated data sets to explore time lags in genetic structure, and argue that time lags should be explicitly considered both when designing and interpreting landscape genetic studies. We conclude that the time lag problem can be exploited to strengthen inferences about recent landscape changes and to establish conservation baselines, particularly when genetic data are combined with other data.     

Autosomal, mitochondrial, and Y chromosome DNA variation in Finland: evidence for a male-specific bottleneck

The high prevalence of rare genetic diseases in Finland has been attributed to a founder effect some 2,000 years ago. However, this hypothesis has not been supported from mtDNA sequence and autosomal microsatellite data which indicate high levels of gene diversity. Here we have identified genetic evidence for a population bottleneck by examining variable microsatellite loci on the nonrecombining portion of Y chromosomes from Finland and four populations from Europe and the Americas. Sequence data from segment I of the control region (HVS-1) of mtDNA (360 bases) and 20 autosomal dinucleotide repeat markers were also analyzed. Partitions of genetic variance within and between populations revealed significant levels of Y-chromosome differentiation between populations. Phylogenetic and diversity analyses revealed divergent Finnish Y-haplotype clades and significantly lower Y-haplotype diversity among Finns as compared to other populations. Surprisingly, Finnish Y-haplotype diversity was even lower than the Native American populations. These results provide support for the Finnish bottleneck hypothesis. Evidence for two separate founding Finnish Y-chromosome lineages was also observed from the Y-chromosome phylogeny. A limited number of closely related founding males may have contributed to the low number of paternal lineages in the Finnish population. In contrast, high levels of genetic diversity for mtDNA and autosomal STRs may be the result of sex-biased gene flow and recent immigration to urban areas from established internal isolates within Finland.     

Spatial population structure in a patchily distributed beetle

The dynamics and evolution of populations will critically depend on their spatial structure. Hence, a recent emphasis on one particular type of structure--the metapopulation concept of Levins--can only be justified by empirical assessment of spatial population structures in a wide range of organisms. This paper focuses on Aphodius fossor, a dung beetle specialized on cattle pastures. An agricultural database was used to locate nearly 50 000 local populations of A. fossor in Finland. Several independent methods were then used to quantify key processes in this vast population system. Allozyme markers and mitochondrial DNA (mtDNA) sequences were applied to examine genetic differentiation of local populations and to derive indirect estimates of gene flow. These estimates were compared to values expected on the basis of direct observations of dispersing individuals and assessments of local effective population size. Molecular markers revealed striking genetic homogeneity in A. fossor. Differentiation was only evident in mtDNA haplotype frequencies between the isolated Aland islands and the Finnish mainland. Thus, indirect estimates of gene flow agreed with direct observations that local effective population size in A. fossor is large (hundreds of individuals), and that in each generation, a substantial fraction (approximately one-fifth) of the individuals move between populations. Large local population size, extreme haplotype diversity and a high regional incidence of A. fossor all testify against recurrent population turnover. Taken together, these results provide strong evidence that the whole mainland population of A. fossor is better described as one large 'patchy population', with substantial movement between relatively persistent local populations, than as a classical metapopulation.     

Population structure of wild bananas, Musa balbisiana, in China determined by SSR fingerprinting and cpDNA PCR-RFLP

Both demographic history and dispersal mechanisms influence the apportionment of genetic diversity among plant populations across geographical regions. In this study, phylogeography and population structure of wild banana, Musa balbisiana, one of the progenitors of cultivated bananas and plantains in China were investigated by an analysis of genetic diversity of simple sequence repeat (SSR) fingerprint markers and cpDNA PCR-RFLP. A chloroplast DNA (cpDNA) genealogy of 21 haplotypes identified two major clades, which correspond to two geographical regions separated by the Beijiang and Xijiang rivers, suggesting a history of vicariance. Significant genetic differentiation was detected among populations with cpDNA markers, a result consistent with limited seed dispersal in wild banana mediated by foraging of rodents. Nuclear SSR data also revealed significant geographical structuring in banana populations. In western China, however, there was no detected phylogeograpahical pattern, possibly due to frequent pollen flow via fruit bats. In contrast, populations east of the Beijiang River and the population of Hainan Island, where long-range soaring pollinators are absent, are genetically distinct. Colonization-extinction processes may have influenced the evolution of Musa populations, which have a metapopulation structure and are connected by migrating individuals. Effective gene flow via pollen, estimated from the nuclear SSR data, is 3.65 times greater than gene flow via seed, estimated from cpDNA data. Chloroplast and nuclear DNAs provide different insights into phylogeographical patterns of wild banana populations and, taken together, can inform conservation practices.     

Genetic diversity, structure, and size of an endangered brown bear population threatened by highway construction in the Pindos Mountains, Greece

One of the major negative effects of roads is the creation of barriers to the movement of wildlife, ultimately disconnecting populations and increasing extinction risk. We collected genetic data from a threatened brown bear population in the central part of the Pindos mountain range in northwestern Greece to provide information about this, as yet genetically undescribed, population and to evaluate its status prior to the construction of a major highway. We used noninvasive genetic sampling methods and microsatellite analysis to investigate nuclear genetic diversity, population genetic structure, demographic history, relatedness within the population and estimated effective and total population size. Brown bears in the study area were found to possess a relatively high level of nuclear genetic diversity and low levels of inbreeding; the population did not show any signs of substructuring but seems to have gone through a genetic bottleneck in the recent past. The estimated effective population size was 29, and the total population size estimate obtained by two different methods was 33 and 51 individuals, respectively. Our results indicate a good conservation status of this bear population and provide baseline genetic data for the future evaluation of the effects on bears from the construction of a major highway, for monitoring the genetic status of this and other bear populations in Greece and for assessing gene flow in bear populations in southern Europe.     

Population genetic structure,differentiation, and diversity in Tetrix subulata pygmy grasshoppers: roles of population size and immigration

Genetic diversity within and among populations and species is influenced by complex demographic and evolutionary processes. Despite extensive research, there is no consensus regarding how landscape structure, spatial distribution, gene flow, and population dynamics impact genetic composition of natural populations. Here, we used amplified fragment length polymorphisms (AFLPs) to investigate effects of population size, geographic isolation, immigration, and gene flow on genetic structure, divergence, and diversity in populations of Tetrix subulata pygmy grasshoppers (Orthoptera: Tetrigidae) from 20 sampling locations in southern Sweden. Analyses of 1564 AFLP markers revealed low to moderate levels of genetic diversity (PPL = 59.5–90.1; Hj = 0.23–0.32) within and significant divergence among sampling localities. This suggests that evolution of functional traits in response to divergent selection is possible and that gene flow is restricted. Genetic diversity increased with population size and with increasing proportion of long‐winged phenotypes (a proxy of recent immigration) across populations on the island of Öland, but not on the mainland. Our data further suggested that the open water separating Öland from the mainland acts as a dispersal barrier that restricts migration and leads to genetic divergence among regions. Isolation by distance was evident for short interpopulation distances on the mainland, but gradually disappeared as populations separated by longer distances were included. Results illustrate that integrating ecological and molecular data is key to identifying drivers of population genetic structure in natural populations. Our findings also underscore the importance of landscape structure and spatial sampling scheme for conclusions regarding the role of gene flow and isolation by distance.     

Population genetic structure in Myrtus communis L. in a chronically fragmented landscape in the Mediterranean: can gene flow counteract habitat perturbation?

Ancient managed landscapes provide ideal opportunities to assess the consequences of habitat fragmentation on the patterns of genetic diversity and gene flow in long-lived plant species. Using amplified fragment length polymorphism (AFLP) and allozyme markers, we quantified seed-mediated gene flow and population genetic diversity and structure in 14 populations of Myrtus communis (myrtle), a common endozoochorous shrub species of forest patches in lowland agricultural Mediterranean areas. Overall, allozyme diversity for myrtle was low (P₉₅   =   25%; A   =   1.411; H_e = 0.085) compared to other known populations, and a significant portion of populations (57%) had lower levels of allelic diversity and/or heterozygosity than expected at random, as shown by simulated resampling of the whole diversity of the landscape. We found significant correlations between allozyme variability and population size and patch isolation, but no significant inbreeding in any population. Genetic differentiation among populations for both allozyme and AFLP markers was significant (Φ_ST = 0.144 and Φ_ST = 0.142, respectively) but an isolation-by-distance pattern was not detected. Assignment tests on AFLP data indicated a high immigration rate in the populations ( ca. 20–22%), likely through effective seed dispersal across the landscape by birds and mammals. Our results suggest that genetic isolation is not the automatic outcome of habitat destruction since substantial levels of seed-mediated gene flow are currently detectable. However, even moderate rates of gene flow seem insufficient in this long-lived species to counteract the genetic erosion and differentiation imposed by chronic habitat destruction.     

Range-wide genetic differentiation of Eugenia dysenterica (Myrtaceae) populations in Brazilian Cerrado

The analysis of geographic patterns in genetic variation has been one of the most important current tools to understand ecological and evolutionary processes underlying population structure. However, inferring such processes from population data may be misleading if biased geographic samples are analyzed. Here we expand previous analyses of Eugenia dysenterica population structure in Brazilian Cerrado, analyzing a larger number of populations distributed throughout a broader geographic region covering most of species' range. We provide new estimates of genetic diversity and population structure based on SSR markers from both neutral and genic regions, using several cluster and ordination techniques. These analyzes reveal a continuous northwestern-southeastern gradient in population differentiation, and not two distinct clusters of populations as suggested in some previous studies. This more comprehensive analysis also reinforces that a simple process of stochastic differentiation do not explain the observed patterns. Moreover, we conclude that explanations for population differentiation may focus on why genetic diversity decreases toward southeastern populations and not necessarily on barriers and interruption of gene flow creating regional patterns of population differentiation.     

Fine-scale population structure and dispersal in Biomphalaria glabrata, the intermediate snail host of Schistosoma mansoni, in Venezuela

Biomphalaria glabrata is the main intermediate host of Schistosoma mansoni in America and one of the most intensely studied species of freshwater snails, yet very little is known about its population biology. Here, we used seven highly polymorphic microsatellite loci to analyse genetic diversity in the Valencia lake basin, which represents the core of the endemic area for schistosomiasis in Venezuela. Populations were sampled at short spatial scale (a few kilometres), both inside the lake and in ponds or rivers near the lake. Our results indicate that B. glabrata essentially cross-fertilizes, with little variation in selfing rates among populations. Our markers detected considerable genetic variation, with an average heterozygosity of 0.60. More diversity per population was found within than outside the lake, suggesting an influence of connectivity among populations on the levels of genetic diversity. A marked population structure was detected and lake populations were less structured than other populations. Most individuals were assigned to their population of origin using an assignment test. No strong demographic signal (e.g. bottleneck) was detected, though lake populations are likely to experience bottlenecks more frequently than the other populations analysed. Differences in gene flow therefore seem to play an important role in population differentiation and in the restoring of genetic diversity in demographically unstable populations.     

Y chromosome haplotype distribution of brown bears (Ursus arctos) in Northern Europe provides insight into population history and recovery

High‐resolution, male‐inherited Y‐chromosomal markers are a useful tool for population genetic analyses of wildlife species, but to date have only been applied in this context to relatively few species besides humans. Using nine Y‐chromosomal STRs and three Y‐chromosomal single nucleotide polymorphism markers (Y‐SNPs), we studied whether male gene flow was important for the recent recovery of the brown bear (Ursus arctos) in Northern Europe, where the species declined dramatically in numbers and geographical distribution during the last centuries but is expanding now. We found 36 haplotypes in 443 male extant brown bears from Sweden, Norway, Finland and northwestern Russia. In 14 individuals from southern Norway from 1780 to 1920, we found two Y chromosome haplotypes present in the extant population as well as four Y chromosome haplotypes not present among the modern samples. Our results suggested major differences in genetic connectivity, diversity and structure between the eastern and the western populations in Northern Europe. In the west, our results indicated that the recovered population originated from only four male lineages, displaying pronounced spatial structuring suggestive of large‐scale population size increase under limited male gene flow within the western subpopulation. In the east, we found a contrasting pattern, with high haplotype diversity and admixture. This first population genetic analysis of male brown bears shows conclusively that male gene flow was not the main force of population recovery.     

Low genetic differentiation in a sedentary bird: house sparrow population genetics in a contiguous landscape

The house sparrow Passer domesticus has been declining in abundance in many localities, including Finland. We studied the genetic diversity and differentiation of the house sparrow populations across Finland in the 1980s, at the onset of the species'' decline in abundance. We genotyped 472 adult males (the less dispersive sex) from 13 locations in Finland (covering a range of 400 × 800 km) and one in Sweden (Stockholm) for 13 polymorphic microsatellite markers. Our analysis of Finnish ringing records showed that natal dispersal distances are limited (90% <16 km), which confirmed earlier finding from other countries. The Finnish populations were panmictic, and genetically very homogeneous and the limited dispersal was sufficiently large to maintain their connectivity. However, all Finnish populations differed significantly from the Stockholm population, even though direct geographical distance to it was often smaller than among Finnish populations. Hence, the open sea between Finland and Sweden appears to form a dispersal barrier for this species, whereas dispersal is much less constrained across the Finnish mainland (which lacks geographical barriers). Our findings provide a benchmark for conservation biologists and emphasize the influence of landscape structure on gene flow.     

Development of genetic diversity,differentiation and structure over 500 years in four ponderosa pine populations

Population history plays an important role in shaping contemporary levels of genetic variation and geographic structure. This is especially true in small, isolated range‐margin populations, where effects of inbreeding, genetic drift and gene flow may be more pronounced than in large continuous populations. Effects of landscape fragmentation and isolation distance may have implications for persistence of range‐margin populations if they are demographic sinks. We studied four small, disjunct populations of ponderosa pine over a 500‐year period. We coupled demographic data obtained through dendroecological methods with microsatellite data to discern how and when contemporary levels of allelic diversity, among and within‐population levels of differentiation, and geographic structure, arose. Alleles accumulated rapidly following initial colonization, demonstrating proportionally high levels of gene flow into the populations. At population sizes of approximately 100 individuals, allele accumulation saturated. Levels of genetic differentiation among populations (F_ST and Jost's D_est) and diversity within populations (F_IS) remained stable through time. There was no evidence of geographic genetic structure at any time in the populations' history. Proportionally, high gene flow in the early stages of population growth resulted in rapid accumulation of alleles and quickly created relatively homogenous genetic patterns among populations. Our study demonstrates that contemporary levels of genetic diversity were formed quickly and early in population development. How contemporary genetic diversity accumulates over time is a key facet of understanding population growth and development. This is especially relevant given the extent and speed at which species ranges are predicted to shift in the coming century.     

Birth of a hotspot of intraspecific genetic diversity: notes from the underground

Hotspots of intraspecific diversity have been observed in most species, often within areas of putative Pleistocene refugia. They have thus mostly been viewed as the outcome of prolonged stability of large populations within the refugia. However, recent evidence has suggested that several other microevolutionary processes could also be involved in their formation. Here, we investigate the contribution of these processes to current range-wide patterns of genetic diversity in the Italian endemic mole Talpa romana, using both nuclear (30 allozyme loci) and mitochondrial markers (cytochrome b sequences). Southern populations of this species showed an allozyme variation that is amongst the highest observed in small mammals (most populations had an expected heterozygosity of 0.10 or above), which was particularly unexpected for a subterranean species. Population genetic, phylogeographic and historical demographic analyses indicated that T. romana populations repeatedly underwent allopatric differentiations followed by secondary admixture within the refugial range in southern Italy. A prolonged demographic stability was reliably inferred from the mitochondrial DNA data only for a population group located north and east of the Calabrian peninsula, showing comparatively lower levels of allozyme variability, and lacking evidence of secondary admixture with other groups. Thus, our results point to the admixture between differentiated lineages as the main cause of the higher levels of diversity of refugial populations. When compared with the Pleistocene evolutionary history recently inferred for species from both the same and other geographic regions, these results suggest the need for a reappraisal of the role of gene exchange in the formation of intraspecific hotspots of genetic diversity.     

Genetic drift and rapid evolution of viviparity in insular fire salamanders (Salamandra salamandra)

Continental islands offer an excellent opportunity to investigate adaptive processes and to time microevolutionary changes that precede macroevolutionary events. We performed a population genetic study of the fire salamander (Salamandra salamandra), a species that displays unique intraspecific diversity of reproductive strategies, to address the microevolutionary processes leading to phenotypic and genetic differentiation of island, coastal and interior populations. We used eight microsatellite markers to estimate genetic diversity, population structure and demographic parameters in viviparous insular populations and ovoviviparous coastal and interior populations. Our results show considerable genetic differentiation (F(ST) range: 0.06-0.27), and no clear signs of gene flow among populations, except between the large and admixed interior populations. We find no support for island colonization by rafting or intentional/accidental anthropogenic introductions, indicating that rising sea levels were responsible for isolation of the island populations approximately 9000 years ago. Our study provides evidence of rapid genetic differentiation between island and coastal populations, and rapid evolution of viviparity driven by climatic selective pressures on island populations, geographic isolation with genetic drift, or a combination of these factors. Studies of these viviparous island populations in early stages of divergence help us better understand the microevolutionary processes involved in rapid phenotypic shifts.     

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.     

It's about time: divergence, demography, and the evolution of developmental modes in marine invertebrates

Differences in larval developmental mode are predicted to affect ecological and evolutionary processes ranging from gene flow and population bottlenecks to rates of population recovery from anthropogenic disturbance and capacity for local adaptation. The most powerful tests of these predictions use comparisons among species to ask how phylogeographic patterns are correlated with the evolution and loss of prolonged planktonic larval development. An important and largely untested assumption of these studies is that interspecific differences in population genetic structure are mainly caused by differences in dispersal and gene flow (rather than by differences in divergence times among populations or changes in effective population sizes), and that species with similar patterns of spatial genetic variation have similar underlying temporal demographic histories. Teasing apart these temporal and spatial patterns is important for understanding the causes and consequences of evolutionary changes in larval developmental mode. New analytical methods that use the coalescent history of allelic diversity can reveal these temporal patterns, test the strength of traditional population-genetic explanations for variation in spatial structure based on differences in dispersal, and identify strongly supported alternative explanations for spatial structure based on demographic history rather than on gene flow alone. We briefly review some of these recent analytical developments, and show their potential for refining ideas about the correspondence between the evolution of larval developmental mode, population demographic history, and spatial genetic variation.     

Population genetic structure of Bellamya aeruginosa (Mollusca: Gastropoda: Viviparidae) in China: weak divergence across large geographic distances

Bellamya aeruginosa is a widely distributed Chinese freshwater snail that is heavily harvested, and its natural habitats are under severe threat due to fragmentation and loss. We were interested whether the large geographic distances between populations and habitat fragmentation have led to population differentiation and reduced genetic diversity in the species. To estimate the genetic diversity and population structure of B. aeruginosa, 277 individuals from 12 populations throughout its distribution range across China were sampled: two populations were sampled from the Yellow River system, eight populations from the Yangtze River system, and two populations from isolated plateau lakes. We used seven microsatellite loci and mitochondrial cytochrome oxidase I sequences to estimate population genetic parameters and test for demographic fluctuations. Our results showed that (1) the genetic diversity of B. aeruginosa was high for both markers in most of the studied populations and effective population sizes appear to be large, (2) only very low and mostly nonsignificant levels of genetic differentiation existed among the 12 populations, gene flow was generally high, and (3) relatively weak geographic structure was detected despite large geographic distances between populations. Further, no isolation by linear or stream distance was found among populations within the Yangtze River system and no signs of population bottlenecks were detected. Gene flow occurred even between far distant populations, possibly as a result of passive dispersal during flooding events, zoochoric dispersal, and/or anthropogenic translocations explaining the lack of stronger differentiation across large geographic distances. The high genetic diversity of B. aeruginosa and the weak population differentiation are likely the results of strong gene flow facilitated by passive dispersal and large population sizes suggesting that the species currently is not of conservation concern.     

Molecular phylogeography and population evolution analysis of <Emphasis Type="Italic">Citrus ichangensis</Emphasis> (Rutaceae)

Ichang papeda (Citrus ichangensis), a wild and endemic perennial plant in Rutaceae, is characterized by the existence of wild and natural populations in southwestern and middle-west China. We analyzed a total of 231 individuals across 16 natural populations using chloroplast SSR markers, nuclear SSR markers, and single-copy nuclear genes. Standard population genetic analyses as well as Bayesian and maximum likelihood models were used to clarify the genetic diversity, population differentiation, barriers to gene flow, bottleneck events, isolation by distance, history migration, demographic history among populations, and phylogeny evolution. The chloroplast and nuclear genome analyses revealed a low level of genetic diversity in C. ichangensis. Clear signals of recent bottlenecks and strong patterns of isolation by distance were detected among different subpopulations, indicating a low extent of historical gene flow for this species and that genetic drift would occur after population differentiation. Bayesian clustering analyses revealed a clear pattern of genetic structure, with one cluster spanning the potential refugia in Wuling Mountains and Ta-pa Mountains, and other two clusters covering a more limited distribution range. The demographic history also supported the scenario that two isolated clusters originated in parallel from the genetic diversity center. Taxonomically, Ichang papeda may be a member of subgenus Citrus. Owing to the complicated topography, the mountainous regions and the Yangtze River have provided long-term stable habitats for C. ichangensis and acted as main barriers for its expansion, which might facilitate the process of speciation. Statistical population models and genetic data indicated strong genetic structure in C. ichangensis, which might result from the restricted gene flow, genetic drift, and population bottlenecks.