首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Patterns of local adaptation at fine spatial scales are central to understanding how evolution proceeds, and are essential to the effective management of economically and ecologically important forest tree species. Here, we employ single and multilocus analyses of genetic data (= 116 231 SNPs) to describe signatures of fine‐scale adaptation within eight whitebark pine (Pinus albicaulis Engelm.) populations across the local extent of the environmentally heterogeneous Lake Tahoe Basin, USA. We show that despite highly shared genetic variation (FST = 0.0069), there is strong evidence for adaptation to the rain shadow experienced across the eastern Sierra Nevada. Specifically, we build upon evidence from a common garden study and find that allele frequencies of loci associated with four phenotypes (mean = 236 SNPs), 18 environmental variables (mean = 99 SNPs), and those detected through genetic differentiation (n = 110 SNPs) exhibit significantly higher signals of selection (covariance of allele frequencies) than could be expected to arise, given the data. We also provide evidence that this covariance tracks environmental measures related to soil water availability through subtle allele frequency shifts across populations. Our results replicate empirical support for theoretical expectations of local adaptation for populations exhibiting strong gene flow and high selective pressures and suggest that ongoing adaptation of many P. albicaulis populations within the Lake Tahoe Basin will not be constrained by the lack of genetic variation. Even so, some populations exhibit low levels of heritability for the traits presumed to be related to fitness. These instances could be used to prioritize management to maintain adaptive potential. Overall, we suggest that established practices regarding whitebark pine conservation be maintained, with the additional context of fine‐scale adaptation.  相似文献   

2.
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.  相似文献   

3.
Striking genetic structure among marine populations at small spatial scales is becoming evident with extensive molecular studies. Such observations suggest isolation at small scales may play an important role in forming patterns of genetic diversity within species. Isolation‐by‐distance, isolation‐by‐environment and historical priority effects are umbrella terms for a suite of processes that underlie genetic structure, but their relative importance at different spatial and temporal scales remains elusive. Here, we use marine lakes in Indonesia to assess genetic structure and assess the relative roles of the processes in shaping genetic differentiation in populations of a bivalve mussel (Brachidontes sp.). Marine lakes are landlocked waterbodies of similar age (6,000–10,000 years), but with heterogeneous environments and varying degrees of connection to the sea. Using a population genomic approach (double‐digest restriction‐site‐associated DNA sequencing), we show strong genetic structuring across populations (range FST: 0.07–0.24) and find limited gene flow through admixture plots. At large spatial scales (>1,400 km), a clear isolation‐by‐distance pattern was detected. At smaller spatial scales (<200 km), this pattern is maintained, but accompanied by an association of genetic divergence with degree of connection. We hypothesize that (incomplete) dispersal barriers can cause initial isolation, allowing priority effects to give the numerical advantage necessary to initiate strong genetic structure. Priority effects may be strengthened by local adaptation, which the data may corroborate by showing a high correlation between mussel genotypes and temperature. Our study indicates an often‐neglected role of (evolution‐mediated) priority effects in shaping population divergence.  相似文献   

4.
Detecting and estimating the degree of genetic differentiation among populations of highly mobile marine fish having pelagic larval stages is challenging because their effective population sizes can be large, and thus, little genetic drift and differentiation is expected in neutral genomic sites. However, genomic sites subject to directional selection stemming from variation in local environmental conditions can still show substantial genetic differentiation, yet these signatures can be hard to detect with low‐throughput approaches. Using a pooled RAD‐seq approach, we investigated genomewide patterns of genetic variability and differentiation within and among 20 populations of Atlantic herring in the Baltic Sea (and adjacent Atlantic sites), where previous low‐throughput studies and/or studies based on few populations have found limited evidence for genetic differentiation. Stringent quality control was applied in the filtering of 1 791 254 SNPs, resulting in a final data set of 68 182 polymorphic loci. Clear differentiation was identified between Atlantic and Baltic populations in many genomic sites, while differentiation within the Baltic Sea area was weaker and geographically less structured. However, outlier analyses – whether including all populations or only those within the Baltic Sea – uncovered hundreds of directionally selected loci in which variability was associated with either salinity, temperature or both. Hence, our results support the view that although the degree of genetic differentiation among Baltic Sea herring populations is low, there are many genomic regions showing elevated divergence, apparently as a response to temperature‐ and salinity‐related natural selection. As such, the results add to the increasing evidence of local adaptation in highly mobile marine organisms, and those in the young Baltic Sea in particular.  相似文献   

5.
Host–parasite co‐evolution can lead to genetic differentiation among isolated host–parasite populations and local adaptation between parasites and their hosts. However, tests of local adaptation rarely consider multiple fitness‐related traits although focus on a single component of fitness can be misleading. Here, we concomitantly examined genetic structure and co‐divergence patterns of the trematode Coitocaecum parvum and its crustacean host Paracalliope fluviatilis among isolated populations using the mitochondrial cytochrome oxidase I gene (COI). We then performed experimental cross‐infections between two genetically divergent host–parasite populations. Both hosts and parasites displayed genetic differentiation among populations, although genetic structure was less pronounced in the parasite. Data also supported a co‐divergence scenario between C. parvum and P. fluviatilis potentially related to local co‐adaptation. Results from cross‐infections indicated that some parasite lineages seemed to be locally adapted to their sympatric (home) hosts in which they achieved higher infection and survival rates than in allopatric (away) amphipods. However, local, intrinsic host and parasite characteristics (host behavioural or immunological resistance to infections, parasite infectivity or growth rate) also influenced patterns of host–parasite interactions. For example, overall host vulnerability to C. parvum varied between populations, regardless of parasite origin (local vs. foreign), potentially swamping apparent local co‐adaptation effects. Furthermore, local adaptation effects seemed trait specific; different components of parasite fitness (infection and survival rates, growth) responded differently to cross‐infections. Overall, data show that genetic differentiation is not inevitably coupled with local adaptation, and that the latter must be interpreted with caution in a multi‐trait context.  相似文献   

6.
Local adaptation at range edges influences species’ distributions and how they respond to environmental change. However, the factors that affect adaptation, including gene flow and local selection pressures, are likely to vary across different types of range edge. We performed a reciprocal transplant experiment to investigate local adaptation in populations of Plantago lanceolata and P. major from central locations in their European range and from their latitudinal and elevation range edges (in northern Scandinavia and Swiss Alps, respectively). We also characterized patterns of genetic diversity and differentiation in populations using molecular markers. Range‐centre plants of P. major were adapted to conditions at the range centre, but performed similarly to range‐edge plants when grown at the range edges. There was no evidence for local adaptation when comparing central and edge populations of P. lanceolata. However, plants of both species from high elevation were locally adapted when compared with plants from high latitude, although the reverse was not true. This asymmetry was associated with greater genetic diversity and less genetic differentiation over the elevation gradient than over the latitudinal gradient. Our results suggest that adaptation in some range‐edge populations could increase their performance following climate change. However, responses are likely to differ along elevation and latitudinal gradients, with adaptation more likely at high‐elevation. Furthermore, based upon these results, we suggest that gene flow is unlikely to constrain adaptation in range‐edge populations of these species.  相似文献   

7.
The distribution, spatial pattern and population dynamics of a species can be influenced by differences in the environment across its range. Spatial variation in climatic conditions can cause local populations to undergo disruptive selection and ultimately result in local adaptation. However, local adaptation can be constrained by gene flow and may favour resident individuals over migrants—both are factors critical to the assessment of invasion potential. The Natal fruit fly (Ceratitis rosa) is a major agricultural pest in Africa with a history of island invasions, although its range is largely restricted to south east Africa. Across Africa, C. rosa is genetically structured into two clusters (R1 and R2), with these clusters occurring sympatrically in the north of South Africa. The spatial distribution of these genotypic clusters remains unexamined despite their importance for understanding the pest's invasion potential. Here, C. rosa, sampled from 22 South African locations, were genotyped at 11 polymorphic microsatellite loci and assessed morphologically using geometric morphometric wing shape analyses to investigate patterns of population structure and determine connectedness of pest‐occupied sites. Our results show little to no intraspecific (population) differentiation, high population connectivity, high effective population sizes and only one morphological type (R2) within South Africa. The absence of the R1 morphotype at sites where it was previously found may be a consequence of differences in thermal niches of the two morphotypes. Overall, our results suggest high invasion potential of this species, that area‐wide pest management should be undertaken on a country‐wide scale, and that border control is critical to preventing further invasions.  相似文献   

8.
Defining the scale of connectivity among marine populations and identifying the barriers to gene flow are tasks of fundamental importance for understanding the genetic structure of populations and for the design of marine reserves. Here, we investigated the population genetic structure at three spatial scales of the red gorgonian Paramuricea clavata (Cnidaria, Octocorallia), a key species dwelling in the coralligenous assemblages of the Mediterranean Sea. Colonies of P. clavata were collected from 39 locations across the Mediterranean Sea from Morocco to Turkey and analysed using microsatellite loci. Within three regions (Medes, Marseille and North Corsica), sampling was obtained from multiple locations and at different depths. Three different approaches (measures of genetic differentiation, Bayesian clustering and spatially explicit maximum‐difference algorithm) were used to determine the pattern of genetic structure. We identified genetic breaks in the spatial distribution of genetic diversity, which were concordant with oceanographic conditions in the Mediterranean Sea. We revealed a high level of genetic differentiation among populations and a pattern of isolation by distance across the studied area and within the three regions, underlining short effective larval dispersal in this species. We observed genetic differentiation among populations in the same locality dwelling at different depths, which may be explained by local oceanographic conditions and which may allow a process of local adaptation of the populations to their environment. We discuss the implications of our results for the conservation of the species, which is exposed to various threats.  相似文献   

9.
Ecological parameters vary in space, and the resulting heterogeneity of selective forces can drive adaptive population divergence. Clinal variation represents a classical model to study the interplay of gene flow and selection in the dynamics of this local adaptation process. Although geographic variation in phenotypic traits in discrete populations could be remainders of past adaptation, maintenance of adaptive clinal variation requires recurrent selection. Clinal variation in genetically determined traits is generally attributed to adaptation of different genotypes to local conditions along an environmental gradient, although it can as well arise from neutral processes. Here, we investigated whether selection accounts for the strong clinal variation observed in a highly heritable pheomelanin‐based color trait in the European barn owl by comparing spatial differentiation of color and of neutral genes among populations. Barn owl's coloration varies continuously from white in southwestern Europe to reddish‐brown in northeastern Europe. A very low differentiation at neutral genetic markers suggests that substantial gene flow occurs among populations. The persistence of pronounced color differentiation despite this strong gene flow is consistent with the hypothesis that selection is the primary force maintaining color variation among European populations. Therefore, the color cline is most likely the result of local adaptation.  相似文献   

10.
11.
1. Determined by landscape structure as well as dispersal‐related traits of species, connectivity influences various key aspects of population biology, ranging from population persistence to genetic structure and diversity. Here, we investigated differences in small‐scale connectivity in terms of gene flow between populations of two ecologically important invertebrates with contrasting dispersal‐related traits: an amphipod (Gammarus fossarum) with a purely aquatic life cycle and a mayfly (Baetis rhodani) with a terrestrial adult stage. 2. We used highly polymorphic markers to estimate genetic differentiation between populations of both species within a Swiss pre‐alpine catchment and compared these results to the broader‐scale genetic structure within the Rhine drainage. Landscape genetic approaches were used to test for correlations of genetic and geographical structures and in‐stream barrier effects. 3. We found overall very weak genetic structure in populations of B. rhodani. In contrast, G. fossarum showed strong genetic differentiation, even at spatial scales of a few kilometres, and a clear pattern of isolation by distance. Genetic diversity decreased from downstream towards upstream populations of G. fossarum, suggesting asymmetric gene flow. Correlation of genetic structure with landscape topography was more pronounced in the amphipod. Our study also indicates that G. fossarum might be capable of dispersing overland in headwater regions and of crossing small in‐stream barriers. 4. We speculate that differences in dispersal capacity but also habitat specialisation and potentially the extent of local adaptation could be responsible for the differences in genetic differentiation found between the two species. These results highlight the importance of taking into account dispersal‐related traits when planning management and conservation strategies.  相似文献   

12.
Potential declines in native pollinator communities and increased reliance on pollinator‐dependent crops have raised concerns about native pollinator conservation and dispersal across human‐altered landscapes. Bumble bees are one of the most effective native pollinators and are often the first to be extirpated in human‐altered habitats, yet little is known about how bumble bees move across fine spatial scales and what landscapes promote or limit their gene flow. In this study, we examine regional genetic differentiation and fine‐scale relatedness patterns of the yellow‐faced bumble bee, Bombus vosnesenskii, to investigate how current and historic habitat composition impact gene flow. We conducted our study across a landscape mosaic of natural, agricultural and urban/suburban habitats, and we show that B. vosnesenskii exhibits low but significant levels of differentiation across the study system (FST = 0.019, Dest = 0.049). Most importantly, we reveal significant relationships between pairwise FST and resistance models created from contemporary land use maps. Specifically, B. vosnesenskii gene flow is most limited by commercial, industrial and transportation‐related impervious cover. Finally, our fine‐scale analysis reveals significant but declining relatedness between individuals at the 1–9 km spatial scale, most likely due to local queen dispersal. Overall, our results indicate that B. vosnesenskii exhibits considerable local dispersal and that regional gene flow is significantly limited by impervious cover associated with urbanization.  相似文献   

13.
The spatial genetic structure of populations is strongly influenced by current and historical patterns of gene flow and drift, which in the simplest case, is limited by geographic distance. We examined the microspatial genetic structure within 33 populations of song sparrows (Melospiza melodia) which included eight subspecies located across coastal areas in southern British Columbia (BC) and California. We also examined the effect of water barriers and local density estimates on genetic structuring. Across both regions, positive genetic structure was detectable at distances of less than 10 km. Genetic divergence was highest in Californian subspecies, perhaps due to reduced gene flow across sub-specific contact zones. In BC, populations distributed across islands displayed greater genetic structuring over similar spatial scales than those across mainland sites, supporting the prediction that water barriers reduce gene flow in this species. Our results confirm both the expectation for fine-scale genetic structure in these generally sedentary subspecies, and the role of landscape features in generating geographic variation in genetic structure.  相似文献   

14.
Aim Climatic changes and fluctuations in the past have strongly influenced the distribution of animal and plant species. Such fluctuations are also reflected in the patterns of genetic diversity on both local and global scales. The genetic pattern of the pearly heath butterfly, Coenonympha arcania, was used to evaluate the genetic differentiation of isolated (in north‐western Europe), peripheral (in north‐eastern Europe) and central (in southern Europe) populations in the context of post‐glacial distributional changes of the species. Location Europe (Sweden, Germany, the Baltic states, Italy, Slovenia, Hungary, Romania, Bulgaria). Thus, samples were collected from large parts of the species’ distribution representing the three categories mentioned above. Methods We analysed 18 loci of 569 individuals from 28 populations by allozyme electrophoresis. We used both individual‐based and population‐based analyses, including F‐statistics, various clustering methods and Markov chain Monte Carlo simulations. Results All loci, except Fum, were polymorphic. The mean FST for all samples was 0.18. The mean genetic distance among populations was 0.046. Two major genetic lineages were distinguished. Populations from the centre of the distributional range in southern Europe and the northern periphery of the distributional range differed significantly in their level of genetic variability. The central populations of south‐eastern Europe showed high levels of genetic diversity and no differentiation among populations. Main conclusions Most probably the two major genetic lineages evolved during glacial isolation in two disjunct Mediterranean refugia. The lack of genetic differentiation across south‐eastern Europe implies a continuous Würm ice age distribution in this area, thus supporting the functional existence of steppe forests throughout this region. The peripheral‐isolated populations in Sweden seem to have suffered from one or more severe bottlenecks, resulting in substantial genetic impoverishment. The peripheral‐connected eastern Baltic populations, on the other hand, are affected by post‐glacial and possibly recurrent gene flow from more central parts of the distribution.  相似文献   

15.
Global climate change poses a serious threat to the future health of coral reef ecosystems. This calls for management strategies that are focused on maximizing the evolutionary potential of coral reefs. Fundamental to this is an accurate understanding of the spatial genetic structure in dominant reef‐building coral species. In this study, we apply a genotyping‐by‐sequencing approach to investigate genome‐wide patterns of genetic diversity, gene flow, and local adaptation in a reef‐building coral, Pocillopora damicornis, across 10 degrees of latitude and a transition from temperate to tropical waters. We identified strong patterns of differentiation and reduced genetic diversity in high‐latitude populations. In addition, genome‐wide scans for selection identified a number of outlier loci putatively under directional selection with homology to proteins previously known to be involved in heat tolerance in corals and associated with processes such as photoprotection, protein degradation, and immunity. This study provides genomic evidence for both restricted gene flow and local adaptation in a widely distributed coral species, and highlights the potential vulnerability of leading‐edge populations to rapid environmental change as they are locally adapted, reproductively isolated, and have reduced levels of genetic diversity.  相似文献   

16.
Identifying processes underlying the genetic and morphological differences among populations is a central question of evolutionary biology. Forest trees typically contain high levels of neutral genetic variation, and genetic differences are often correlated with geographic distance between populations [isolation by distance (IBD)] or are due to historic vicariance events [isolation by colonization (IBC)]. In contrast, morphological differences are largely due to local adaptation. Here, we examined genetic (microsatellite) and morphological (from a common garden experiment) variation in Populus nigra L., European black poplar, collected from 13 sites across western Europe and grown in a common garden in Belgium. Significant genetic differentiation was observed, with populations from France displaying greater admixture than the distinct Spanish and central European gene pools, consistent with previously described glacial refugia (IBC). Many quantitative traits displayed a bimodal distribution, approximately corresponding to small‐leaf and large‐leaf ecotypes. Examination of nine climatic variables revealed the sampling locations to have diverse climates, and although the correlation between morphological and climatic differences was significant, the pattern was not consistent with strict local adaptation. Partial Mantel tests based on multivariate summary statistics identified significant residual correlation in comparisons of small‐leaf to large‐leaf ecotypes, and within the small‐leaf samples, but not within large‐leaf ecotypes, indicating that variation within the small‐leaf morphotype in particular may be adaptive. Some small‐leaf populations experience climates very similar to those in large‐leaf sites. We conclude that adaptive differentiation and persistent IBC acted in combination to produce the genetic and morphological patterns observed in P. nigra.  相似文献   

17.
A comprehensive study of the phylogeography and population genetics of the largest wild artiodactyl in the arid and cold‐temperate South American environments, the guanaco (Lama guanicoe) was conducted. Patterns of molecular genetic structure were described using 514 bp of mtDNA sequence and 14 biparentally inherited microsatellite markers from 314 samples. These individuals originated from 17 localities throughout the current distribution across Peru, Bolivia, Argentina and Chile. This confirmed well‐defined genetic differentiation and subspecies designation of populations geographically separated to the northwest (L. g. cacsilensis) and southeast (L. g. guanicoe) of the central Andes plateau. However, these populations are not completely isolated, as shown by admixture prevalent throughout a limited contact zone, and a strong signal of expansion from north to south in the beginning of the Holocene. Microsatellite analyses differentiated three northwestern and 4–5 southeastern populations, suggesting patterns of genetic contact among these populations. Possible genetic refuges were identified, as were source‐sink patterns of gene flow at historical and recent time scales. Conservation and management of guanaco should be implemented with an understanding of these local population dynamics while also considering the preservation of broader adaptive variation and evolutionary processes.  相似文献   

18.
Chromosomal inversions have been implicated in facilitating adaptation in the face of high levels of gene flow, but whether chromosomal fusions also have similar potential remains poorly understood. Atlantic salmon are usually characterized by population structure at multiple spatial scales; however, this is not the case for tributaries of the Miramichi River in North America. To resolve genetic relationships between populations in this system and the potential for known chromosomal fusions to contribute to adaptation, we genotyped 728 juvenile salmon using a 50 K SNP array. Consistent with previous work, we report extremely weak overall population structuring (Global FST = 0.004) and failed to support hierarchical structure between the river's two main branches. We provide the first genomic characterization of a previously described polymorphic fusion between chromosomes 8 and 29. Fusion genomic characteristics included high LD, reduced heterozygosity in the fused homokaryotes, and strong divergence between the fused and the unfused rearrangement. Population structure based on fusion karyotype was five times stronger than neutral variation (FST = 0.019), and the frequency of the fusion was associated with summer precipitation supporting a hypothesis that this rearrangement may contribute local adaptation despite weak neutral differentiation. Additionally, both outlier variation among populations and a polygenic framework for characterizing adaptive variation in relation to climate identified a 250‐Kb region of chromosome 9, including the gene six6 that has previously been linked to age‐at‐maturity and run‐timing for this species. Overall, our results indicate that adaptive processes, independent of major river branching, are more important than neutral processes for structuring these populations.  相似文献   

19.
While the tempo of diversification in biodiversity hotspots has received much attention, the spatial scale of diversification has often been overlooked. Addressing this deficiency requires understanding the drivers of population divergence and the spatial scales at which they operate in species‐rich clades and ecosystems. South Africa's Succulent Karoo (SK) hotspot provides an excellent system for such research, being both compact (ca. 110,000 km2) and home to spectacular in‐situ radiations, such as the ruschioid Aizoaceae. Here we use GBS to document genetic structure in two co‐occurring ruschioid species, at both coarse (>10 km) and fine (<500 m) spatial scales. Where Ruschia burtoniae shows strong between‐population genetic differentiation and no gene flow, Conophytum calculus shows weak differentiation, with high levels of admixture suggesting recent or ongoing gene flow. Community analysis and transplant experiments reveal that R. burtoniae occupies a narrow, low‐pH edaphic niche, and at scales of a few hundred metres, areas of elevated genetic turnover correspond to patches of edaphically unsuitable habitat. In contrast, C. calculus occupies a broader niche and exhibits isolation‐by‐distance without a habitat effect. We suggest that edaphic specialisation, coupled with highly restricted seed and pollen dispersal in heterogeneous landscapes, has played a major role in driving rapid diversification at small spatial scales in this system. However, the contrasting patterns in our study species show that these factors do not influence all organisms uniformly, being strongly modulated by lineage‐specific traits that influence both the spatial scale of gene flow and habitat specificity.  相似文献   

20.
Quantifying spatial genetic structure can reveal the relative influences of contemporary and historic factors underlying localized and regional patterns of genetic diversity and gene flow – important considerations for the development of effective conservation efforts. Using 10 polymorphic microsatellite loci, we characterize genetic variation among populations across the range of the Eastern Sand Darter (Ammocrypta pellucida), a small riverine percid that is highly dependent on sandy substrate microhabitats. We tested for fine scale, regional, and historic patterns of genetic structure. As expected, significant differentiation was detected among rivers within drainages and among drainages. At finer scales, an unexpected lack of within‐river genetic structure among fragmented sandy microhabitats suggests that stratified dispersal resulting from unstable sand bar habitat degradation (natural and anthropogenic) may preclude substantial genetic differentiation within rivers. Among‐drainage genetic structure indicates that postglacial (14 kya) drainage connectivity continues to influence contemporary genetic structure among Eastern Sand Darter populations in southern Ontario. These results provide an unexpected contrast to other benthic riverine fish in the Great Lakes drainage and suggest that habitat‐specific fishes, such as the Eastern Sand Darter, can evolve dispersal strategies that overcome fragmented and temporally unstable habitats.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号