首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 78 毫秒
1.

Background and Aims

It is widely accepted that hydraulic failure due to xylem embolism is a key factor contributing to drought-induced mortality in trees. In the present study, an attempt is made to disentangle phenotypic plasticity from genetic variation in hydraulic traits across the entire distribution area of a tree species to detect adaptation to local environments.

Methods

A series of traits related to hydraulics (vulnerability to cavitation and hydraulic conductivity in branches), growth performance and leaf mass per area were assessed in eight Pinus canariensis populations growing in two common gardens under contrasting environments. In addition, the neutral genetic variability (FST) and the genetic differentiation of phenotypic variation (QST) were compared in order to identify the evolutionary forces acting on these traits.

Key Results

The variability for hydraulic traits was largely due to phenotypic plasticity. Nevertheless, the vulnerability to cavitation displayed a significant genetic variability (approx. 5 % of the explained variation), and a significant genetic × environment interaction (between 5 and 19 % of the explained variation). The strong correlation between vulnerability to cavitation and survival in the xeric common garden (r = –0·81; P < 0·05) suggests a role for the former in the adaptation to xeric environments. Populations from drier sites and higher temperature seasonality were less vulnerable to cavitation than those growing at mesic sites. No trade-off between xylem safety and efficiency was detected. QST of parameters of the vulnerability curve (0·365 for P50 and the slope of the vulnerability curve and 0·452 for P88) differed substantially from FST (0·091), indicating divergent selection. In contrast, genetic drift alone was found to be sufficient to explain patterns of differentiation for xylem efficiency and growth.

Conclusions

The ability of P. canariensis to inhabit a wide range of ecosystems seemed to be associated with high phenotypic plasticity and some degree of local adaptations of xylem and leaf traits. Resistance to cavitation conferred adaptive potential for this species to adapt successfully to xeric conditions.  相似文献   

2.
Genetic adaptation and phenotypic plasticity are two ways in which organisms can adapt to local environmental conditions. We examined genetic and plastic variation in gill and brain size among swamp (low oxygen; hypoxic) and river (normal oxygen; normoxic) populations of an African cichlid fish, Pseudocrenilabrus multicolor victoriae. Larger gills and smaller brains should be advantageous when oxygen is low, and we hypothesized that the relative contribution of local genetic adaptation vs. phenotypic plasticity should be related to potential for dispersal between environments (because of gene flow’s constraint on local genetic adaptation). We conducted a laboratory‐rearing experiment, with broods from multiple populations raised under high‐oxygen and low‐oxygen conditions. We found that most of the variation in gill size was because of plasticity. However, both plastic and genetic effects on brain mass were detected, as were genetic effects on brain mass plasticity. F1 offspring from populations with the highest potential for dispersal between environments had characteristically smaller and more plastic brains. This phenotypic pattern might be adaptive in the face of gene flow, if smaller brains and increased plasticity confer higher average fitness across environment types.  相似文献   

3.
Adaptation to local environmental conditions and the range dynamics of populations can influence evolutionary divergence along environmental gradients. Thus, it is important to investigate patterns of both phenotypic and genetic variations among populations to reveal the respective roles of these two types of factors in driving population differentiation. Here, we test for evidence of phenotypic and genetic structure across populations of a passerine bird (Zosterops borbonicus) distributed along a steep elevational gradient on the island of Réunion. Using 11 microsatellite loci screened in 401 individuals from 18 localities distributed along the gradient, we found that genetic differentiation occurred at two spatial levels: (i) between two main population groups corresponding to highland and lowland areas, respectively, and (ii) within each of these two groups. In contrast, several morphological traits varied gradually along the gradient. Comparison of neutral genetic differentiation (FST) and phenotypic differentiation (PST) showed that PST largely exceeds FST at several morphological traits, which is consistent with a role for local adaptation in driving morphological divergence along the gradient. Overall, our results revealed an area of secondary contact midway up the gradient between two major, cryptic, population groups likely diverged in allopatry. Remarkably, local adaptation has shaped phenotypic differentiation irrespective of population history, resulting in different patterns of variation along the elevational gradient. Our findings underscore the importance of understanding both historical and selective factors when trying to explain variation along environmental gradients.  相似文献   

4.
Knowing how microevolutionary processes, such as genetic drift and natural selection, shape variation in adaptive traits is strategic for conservation measures. One way to estimate local adaptation is to compare divergences in quantitative traits (QST) and neutral loci (FST). Therefore, we have assessed the pattern of phenotypic and molecular genetic divergence among natural subpopulations of the fruit tree Eugenia dysenterica DC. A provenance and progeny test was performed to assess the quantitative traits of the subpopulations collected in a wide distribution area of the species in the Brazilian Cerrado. The sampled environments are in a biodiversity hotspot with heterogeneous soil and climate conditions. By associating quantitative trait variation in initial seedling development with neutral microsatellite marker variation, we tested the local adaptation of the traits by the QSTFST contrast. Genetic drift was prevalent in the phenotypic differentiation among the subpopulations, although the traits seedling emergence time and root green mass, which are relevant for adaptation to the Cerrado climate, showed signs of uniform selection. Our results suggest that E. dysenterica has a spatial genetic structure divided into two large groups, separated by a line that divides the Cerrado biome in a southwestern to northeastern direction. This structure must be taken into account for managing E. dysenterica genetic resources both for conservation and breeding purposes.  相似文献   

5.
Gene flow is often considered to be one of the main factors that constrains local adaptation in a heterogeneous environment. However, gene flow may also lead to the evolution of phenotypic plasticity. We investigated the effect of gene flow on local adaptation and phenotypic plasticity in development time in island populations of the common frog Rana temporaria which breed in pools that differ in drying regimes. This was done by investigating associations between traits (measured in a common garden experiment) and selective factors (pool drying regimes and gene flow from other populations inhabiting different environments) by regression analyses and by comparing pairwise FST values (obtained from microsatellite analyses) with pairwise QST values. We found that the degree of phenotypic plasticity was positively correlated with gene flow from other populations inhabiting different environments (among‐island environmental heterogeneity), as well as with local environmental heterogeneity within each population. Furthermore, local adaptation, manifested in the correlation between development time and the degree of pool drying on the islands, appears to have been caused by divergent selection pressures. The local adaptation in development time and phenotypic plasticity is quite remarkable, because the populations are young (less than 300 generations) and substantial gene flow is present among islands.  相似文献   

6.
Both environmental and genetic influences can result in phenotypic variation. Quantifying the relative contributions of local adaptation and phenotypic plasticity to phenotypes is key to understanding the effect of environmental variation on populations. Identifying the selective pressures that drive divergence is an important, but often lacking, next step. High gene flow between high‐ and low‐altitude common frog (Rana temporaria) breeding sites has previously been demonstrated in Scotland. The aim of this study was to assess whether local adaptation occurs in the face of high gene flow and to identify potential environmental selection pressures that drive adaptation. Phenotypic variation in larval traits was quantified in R. temporaria from paired high‐ and low‐altitude sites using three common temperature treatments. Local adaptation was assessed using QSTFST analyses, and quantitative phenotypic divergence was related to environmental parameters using Mantel tests. Although evidence of local adaptation was found for all traits measured, only variation in larval period and growth rate was consistent with adaptation to altitude. Moreover, this was only evident in the three mountains with the highest high‐altitude sites. This variation was correlated with mean summer and winter temperatures, suggesting that temperature parameters are potentially strong selective pressures maintaining local adaptation, despite high gene flow.  相似文献   

7.
E Luquet  J-P Léna  C Miaud  S Plénet 《Heredity》2015,114(1):69-79
Variation in the environment can induce different patterns of genetic and phenotypic differentiation among populations. Both neutral processes and selection can influence phenotypic differentiation. Altitudinal phenotypic variation is of particular interest in disentangling the interplay between neutral processes and selection in the dynamics of local adaptation processes but remains little explored. We conducted a common garden experiment to study the phenotypic divergence in larval life-history traits among nine populations of the common toad (Bufo bufo) along an altitudinal gradient in France. We further used correlation among population pairwise estimates of quantitative trait (QST) and neutral genetic divergence (FST from neutral microsatellite markers), as well as altitudinal difference, to estimate the relative role of divergent selection and neutral genetic processes in phenotypic divergence. We provided evidence for a neutral genetic differentiation resulting from both isolation by distance and difference in altitude. We found evidence for phenotypic divergence along the altitudinal gradient (faster development, lower growth rate and smaller metamorphic size). The correlation between pairwise QSTs–FSTs and altitude differences suggested that this phenotypic differentiation was most likely driven by altitude-mediated selection rather than by neutral genetic processes. Moreover, we found different divergence patterns for larval traits, suggesting that different selective agents may act on these traits and/or selection on one trait may constrain the evolution on another through genetic correlation. Our study highlighted the need to design more integrative studies on the common toad to unravel the underlying processes of phenotypic divergence and its selective agents in the context of environmental clines.  相似文献   

8.
9.

Background and Aims

The successful spread of invasive plants in new environments is often linked to multiple introductions and a diverse gene pool that facilitates local adaptation to variable environmental conditions. For clonal plants, however, phenotypic plasticity may be equally important. Here the primary adaptive strategy in three non-native, clonally reproducing macrophytes (Egeria densa, Elodea canadensis and Lagarosiphon major) in New Zealand freshwaters were examined and an attempt was made to link observed differences in plant morphology to local variation in habitat conditions.

Methods

Field populations with a large phenotypic variety were sampled in a range of lakes and streams with different chemical and physical properties. The phenotypic plasticity of the species before and after cultivation was studied in a common garden growth experiment, and the genetic diversity of these same populations was also quantified.

Key Results

For all three species, greater variation in plant characteristics was found before they were grown in standardized conditions. Moreover, field populations displayed remarkably little genetic variation and there was little interaction between habitat conditions and plant morphological characteristics.

Conclusions

The results indicate that at the current stage of spread into New Zealand, the primary adaptive strategy of these three invasive macrophytes is phenotypic plasticity. However, while limited, the possibility that genetic diversity between populations may facilitate ecotypic differentiation in the future cannot be excluded. These results thus indicate that invasive clonal aquatic plants adapt to new introduced areas by phenotypic plasticity. Inorganic carbon, nitrogen and phosphorous were important in controlling plant size of E. canadensis and L. major, but no other relationships between plant characteristics and habitat conditions were apparent. This implies that within-species differences in plant size can be explained by local nutrient conditions. All together this strongly suggests that invasive clonal aquatic plants adapt to a wide range of habitats in introduced areas by phenotypic plasticity rather than local adaptation.  相似文献   

10.
In ectotherms, variation in life history traits among populations is common and suggests local adaptation. However, geographic variation itself is not a proof for local adaptation, as genetic drift and gene flow may also shape patterns of quantitative variation. We studied local and regional variation in means and phenotypic plasticity of larval life history traits in the common frog Rana temporaria using six populations from central Sweden, breeding in either open‐canopy or partially closed‐canopy ponds. To separate local adaptation from genetic drift, we compared differentiation in quantitative genetic traits (QST) obtained from a common garden experiment with differentiation in presumably neutral microsatellite markers (FST). We found that R. temporaria populations differ in means and plasticities of life history traits in different temperatures at local, and in FST at regional scale. Comparisons of differentiation in quantitative traits and in molecular markers suggested that natural selection was responsible for the divergence in growth and development rates as well as in temperature‐induced plasticity, indicating local adaptation. However, at low temperature, the role of genetic drift could not be separated from selection. Phenotypes were correlated with forest canopy closure, but not with geographical or genetic distance. These results indicate that local adaptation can evolve in the presence of ongoing gene flow among the populations, and that natural selection is strong in this system.  相似文献   

11.
Parasite‐mediated selection varying across time and space in metapopulations is expected to result in host local adaptation and the maintenance of genetic diversity in disease‐related traits. However, nonadaptive processes like migration and extinction‐(re)colonization dynamics might interfere with adaptive evolution. Understanding how adaptive and nonadaptive processes interact to shape genetic variability in life‐history and disease‐related traits can provide important insights into their evolution in subdivided populations. Here we investigate signatures of spatially fluctuating, parasite‐mediated selection in a natural metapopulation of Daphnia magna. Host genotypes from infected and uninfected populations were genotyped at microsatellite markers, and phenotyped for life‐history and disease traits in common garden experiments. Combining phenotypic and genotypic data a QSTFST‐like analysis was conducted to test for signatures of parasite mediated selection. We observed high variation within and among populations for phenotypic traits, but neither an indication of host local adaptation nor a cost of resistance. Infected populations have a higher gene diversity (Hs) than uninfected populations and Hs is strongly positively correlated with fitness. These results suggest a strong parasite effect on reducing population level inbreeding. We discuss how stochastic processes related to frequent extinction‐(re)colonization dynamics as well as host and parasite migration impede the evolution of resistance in the infected populations. We suggest that the genetic and phenotypic patterns of variation are a product of dynamic changes in the host gene pool caused by the interaction of colonization bottlenecks, inbreeding, immigration, hybrid vigor, rare host genotype advantage and parasitism. Our study highlights the effect of the parasite in ameliorating the negative fitness consequences caused by the high drift load in this metapopulation.  相似文献   

12.
Evaluating the relative importance of neutral and adaptive processes as determinants of population differentiation across environments is a central theme of evolutionary biology. We applied the QSTFST comparison flanked by a direct test for local adaptation to infer the role of climate‐driven selection and gene flow in population differentiation of an annual grass Avena sterilis in two distinct parts of the species range, edge and interior, which represent two globally different climates, desert and Mediterranean. In a multiyear reciprocal transplant experiment, the plants of desert and Mediterranean origin demonstrated home advantage, and population differentiation in several phenotypic traits related to reproduction exceeded neutral predictions, as determined by comparisons of QST values with theoretical FST distributions. Thus, variation in these traits likely resulted from local adaptation to desert and Mediterranean environments. The two separate common garden experiments conducted with different experimental design revealed that two population comparisons, in contrast to multi‐population comparisons, are likely to detect population differences in virtually every trait, but many of these differences reflect effects of local rather than regional environment. We detected a general reduction in neutral (SSR) genetic variation but not in adaptive quantitative trait variation in peripheral desert as compared with Mediterranean core populations. On the other hand, the molecular data indicated intensive gene flow from the Mediterranean core towards desert periphery. Although species range position in our study (edge vs. interior) was confounded with climate (desert vs. Mediterranean), the results suggest that the gene flow from the species core does not have negative consequences for either performance of the peripheral plants or their adaptive potential.  相似文献   

13.
During the process of ecological speciation, reproductive isolation results from divergent natural selection and leads to a positive correlation between genetic divergence and adaptive phenotypic divergence, that is, isolation by adaptation (IBA). In natural populations, phenotypic differentiation is often autocorrelated with geographic distance, making IBA difficult to distinguish from the neutral expectation of isolation by distance (IBD). We examined these two alternatives in a dramatic case of clinal phenotypic variation in an Andean songbird, the Line‐cheeked Spinetail (Cranioleuca antisiensis). At its geographic extremes, this species shows a near threefold difference in body mass (11.5 to 31.0 g) with marked plumage differences. We analysed phenotypic, environmental and genetic data (5,154 SNPs) from 172 individuals and 19 populations sampled along its linear distribution in the Andes. We found that body mass was tightly correlated with environmental temperature, consistent with local adaptation as per Bergmann's rule. Using a PSTFST analysis, we found additional support for natural selection driving body mass differentiation, but these results could also be explained by environment‐mediated phenotypic plasticity. When we assessed the relative support for patterns of IBA and IBD using variance partitioning, we found that IBD was the best explanation for genetic differentiation along the cline. Adaptive phenotypic or environmental divergence can reduce gene flow, a pattern interpreted as evidence of ecological speciation's role in diversification. Our results provide a counterexample to this interpretation. Despite conditions conducive to ecological speciation, our results suggest that dramatic size and environmental differentiation within C. antisiensis are not limiting gene flow.  相似文献   

14.
Phenotypic plasticity can increase tolerance to heterogeneous environments but the elevations and slopes of reaction norms are often population specific. Disruption of locally adapted reaction norms through outcrossing can lower individual viability. Here, we sampled five genetically distinct populations of brown trout (Salmo trutta) from within a river network, crossed them in a full-factorial design, and challenged the embryos with the opportunistic pathogen Pseudomonas fluorescens. By virtue of our design, we were able to disentangle effects of genetic crossing distance from sire and dam effects on early life-history traits. While pathogen infection did not increase mortality, it was associated with delayed hatching of smaller larvae with reduced yolk sac reserves. We found no evidence of a relationship between genetic distance (W, FST) and the expression of early-life history traits. Moreover, hybrids did not differ in phenotypic means or reaction norms in comparison to offspring from within-population crosses. Heritable variation in early life-history traits was found to remain stable across the control and pathogen environments. Our findings show that outcrossing within a rather narrow geographical scale can have neutral effects on F1 hybrid viability at the embryonic stage, i.e. at a stage when environmental and genetic effects on phenotypes are usually large.  相似文献   

15.
Predictions about the fate of species or populations under climate change scenarios typically neglect adaptive evolution and phenotypic plasticity, the two major mechanisms by which organisms can adapt to changing local conditions. As a consequence, we have little understanding of the scope for organisms to track changing environments by in situ adaptation. Here, we use a detailed individual-specific long-term population study of great tits (Parus major) breeding in Wytham Woods, Oxford, UK to parameterise a mechanistic model and thus directly estimate the rate of environmental change to which in situ adaptation is possible. Using the effect of changes in early spring temperature on temporal synchrony between birds and a critical food resource, we focus in particular on the contribution of phenotypic plasticity to population persistence. Despite using conservative estimates for evolutionary and reproductive potential, our results suggest little risk of population extinction under projected local temperature change; however, this conclusion relies heavily on the extent to which phenotypic plasticity tracks the changing environment. Extrapolating the model to a broad range of life histories in birds suggests that the importance of phenotypic plasticity for adjustment to projected rates of temperature change increases with slower life histories, owing to lower evolutionary potential. Understanding the determinants and constraints on phenotypic plasticity in natural populations is thus crucial for characterising the risks that rapidly changing environments pose for the persistence of such populations.  相似文献   

16.
Gene flow between phenotypically divergent populations can disrupt local adaptation or, alternatively, may stimulate adaptive evolution by increasing genetic variation. We capitalised on historical Trinidadian guppy transplant experiments to test the phenotypic effects of increased gene flow caused by replicated introductions of adaptively divergent guppies, which were translocated from high‐ to low‐predation environments. We sampled two native populations prior to the onset of gene flow, six historic introduction sites, introduction sources and multiple downstream points in each basin. Extensive gene flow from introductions occurred in all streams, yet adaptive phenotypic divergence across a gradient in predation level was maintained. Descendants of guppies from a high‐predation source site showed high phenotypic similarity with native low‐predation guppies in as few as ~12 generations after gene flow, likely through a combination of adaptive evolution and phenotypic plasticity. Our results demonstrate that locally adapted phenotypes can be maintained despite extensive gene flow from divergent populations.  相似文献   

17.
Adaptive plasticity is expected to be important when the grain of environmental variation is encompassed in offspring dispersal distance. We investigated patterns of local adaptation, selection and plasticity in an association of plant morphology with fine-scale habitat shifts from oak canopy understory to adjacent grassland habitat in Claytonia perfoliata. Populations from beneath the canopy of oak trees were >90 % broad leaved and large seeded, while plants from adjacent grassland habitat were >90 % linear-leaved and small seeded. In a 2-year study, we used reciprocal transplants and phenotypic selection analysis to investigate local adaptation, selection, plasticity and maternal effects in this trait-environment association. Transgenerational effects were studied by planting offspring of inbred maternal families grown in both environments across the same environments in the second year. Reciprocal transplants revealed local adaptation to habitat type: broad-leaved forms had higher fitness in oak understory and linear-leaved plants had higher fitness in open grassland habitat. Phenotypic selection analyses indicated selection for narrower leaves and lower SLA in open habitat, and selection for broad leaves and intermediate values of SLA in understory. Both plant morphs exhibited plastic responses in traits in the same direction as selection on traits (narrower leaves and lower SLA in open habitat) suggesting that plasticity is adaptive. We detected an adaptive transgenerational effect in which maternal environment influenced offspring fitness; offspring of grassland-reared plants had higher fitness than understory-reared plants when grown in grassland. We did not detect costs of plasticity, but did find a positive association between leaf shape plasticity and fitness in linear-leaved plants in grassland habitat. Together, these findings indicate that fixed differences in trait values corresponding to selection across habitat contribute to local adaptation, but that plasticity and maternal environmental effects may be favored through promotion of survival across heterogeneous environments.  相似文献   

18.
Sitobion avenae (F.) can survive on various plants in the Poaceae, which may select for highly plastic genotypes. But phenotypic plasticity was often thought to be non-genetic, and of little evolutionary significance historically, and many problems related to adaptive plasticity, its genetic basis and natural selection for plasticity have not been well documented. To address these questions, clones of S. avenae were collected from three plants, and their phenotypic plasticity under alternative environments was evaluated. Our results demonstrated that nearly all tested life-history traits showed significant plastic changes for certain S. avenae clones with the total developmental time of nymphs and fecundity tending to have relatively higher plasticity for most clones. Overall, the level of plasticity for S. avenae clones’ life-history traits was unexpectedly low. The factor ‘clone’ alone explained 27.7–62.3% of the total variance for trait plasticities. The heritability of plasticity was shown to be significant in nearly all the cases. Many significant genetic correlations were found between trait plasticities with a majority of them being positive. Therefore, it is evident that life-history trait plasticity involved was genetically based. There was a high degree of variation in selection coefficients for life-history trait plasticity of different S. avenae clones. Phenotypic plasticity for barley clones, but not for oat or wheat clones, was frequently found to be under significant selection. The directional selection of alternative environments appeared to act to decrease the plasticity of S. avenae clones in most cases. G-matrix comparisons showed significant differences between S. avenae clones, as well as quite a few negative covariances (i.e., trade-offs) between trait plasticities. Genetic basis and evolutionary significance of life-history trait plasticity were discussed.  相似文献   

19.
Genetic differentiation along environmental clines is often observed as a result of interplay between gene flow and natural selection. In order to understand the relative roles of these processes in shaping this differentiation, we designed a study in which we used two approaches that have not previously been combined, the Q STF ST comparison and crossbreeding. We examined (1) interpopulation phenotypic and genetic (AFLP) variation, and (2) performance of interpopulation hybrids in a common annual Senecio glaucus. Fitness of interpopulation hybrids (F1 and F2) was tested under simulated population natural conditions in terms of aridity and analyzed for a relationship with (1) spatial distance and (2) environmental differences (amount of annual rainfall). While phenotypic variation corresponded to the clinal changes in aridity along population locations, viz. narrower and longer leaves, longer leaf outgrowths and advanced flowering in more arid environments, the F ST < 0.1 calculated from AFLP data suggested intensive interpopulation gene flow, with little if any contribution of genetic drift. Performance of hybrids in simulated natural environments revealed heterosis in F1, but a hybrid breakdown in F2 generation. These effects were related to both the spatial distance between hybrid parents and their population rainfall differences. The detected clinal phenotypic variation and outbreeding depression in F2 strongly support presence of aridity-induced clinal natural selection, which is matched by the observed Q ST ≫ F ST. From this we conclude that Q STF ST comparison can detect effect of diversifying selection when patterns of phenotypic variation across sampled locations can be reliably predicted from environmental variation.  相似文献   

20.
Understanding how natural selection generates and maintains adaptive genetic diversity in heterogeneous environments is key to predicting the evolutionary response of populations to rapid environmental change. Detecting selection in complex spatial environments remains challenging, especially for threatened species where the effects of strong genetic drift may overwhelm signatures of selection. We carried out a basinwide riverscape genomic analysis in the threatened southern pygmy perch (Nannoperca australis), an ecological specialist with low dispersal potential. High‐resolution environmental data and 5162 high‐quality filtered SNPs were used to clarify spatial population structure and to assess footprints of selection associated with a steep hydroclimatic gradient and with human disturbance across the naturally and anthropogenically fragmented Murray–Darling Basin (Australia). Our approach included FST outlier tests to define neutral loci, and a combination of spatially explicit genotype–environment association analyses to identify candidate adaptive loci while controlling for the effects of landscape structure and shared population history. We found low levels of genetic diversity and strong neutral population structure consistent with expectations based on spatial stream hierarchy and life history. In contrast, variables related to precipitation and temperature appeared as the most important environmental surrogates for putatively adaptive genetic variation at both regional and local scales. Human disturbance also influenced the variation in candidate loci for adaptation, but only at a local scale. Our study contributes to understanding of adaptive evolution along naturally and anthropogenically fragmented ecosystems. It also offers a tangible example of the potential contributions of landscape genomics for informing in situ and ex situ conservation management of biodiversity.  相似文献   

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

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