首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 234 毫秒
1.
Cases of geographically restricted co‐occurring sister taxa are rare and may point to potential divergence with gene flow. The two bat species Murina gracilis and Murina recondita are both endemic to Taiwan and are putative sister species. To test for nonallopatric divergence and gene flow in these taxa, we generated sequences using Sanger and next‐generation sequencing, and combined these with microsatellite data for coalescent‐based analyses. MtDNA phylogenies supported the reciprocally monophyletic sister relationship between M. gracilis and M. recondita; however, clustering of microsatellite genotypes revealed several cases of species admixture suggesting possible introgression. Sequencing of microsatellite flanking regions revealed that admixture signatures stemmed from microsatellite allele homoplasy rather than recent introgressive hybridization, and also uncovered an unexpected sister relationship between M. recondita and the continental species Murina eleryi, to the exclusion of M. gracilis. To dissect the basis of these conflicts between ncDNA and mtDNA, we analysed sequences from 10 anonymous ncDNA loci with *beast and isolation‐with‐migration and found two distinct clades of M. eleryi, one of which was sister to M. recondita. We conclude that Taiwan was colonized by the ancestor of M. gracilis first, followed by the ancestor of M. recondita after a period of allopatric divergence. After colonization, the mitochondrial genome of M. recondita was replaced by that of the resident M. gracilis. This study illustrates how apparent signatures of sympatric divergence can arise from complex histories of allopatric divergence, colonization and hybridization, thus highlighting the need for rigorous analyses to distinguish between such scenarios.  相似文献   

2.
Historical events, habitat preferences, and geographic barriers might result in distinct genetic patterns in insular versus mainland populations. Comparison between these two biogeographic systems provides an opportunity to investigate the relative role of isolation in phylogeographic patterns and to elucidate the importance of evolution and demographic history in population structure. Herein, we use a genotype‐by‐sequencing approach (GBS) to explore population structure within three species of mastiff bats (Molossus molossus, M. coibensis, and M. milleri), which represent different ecological histories and geographical distributions in the genus. We tested the hypotheses that oceanic straits serve as barriers to dispersal in Caribbean bats and that isolated island populations are more likely to experience genetic drift and bottlenecks in comparison with highly connected ones, thus leading to different phylogeographic patterns. We show that population structures vary according to general habitat preferences, levels of population isolation, and historical fluctuations in climate. In our dataset, mainland geographic barriers played only a small role in isolation of lineages. However, oceanic straits posed a partial barrier to the dispersal for some populations within some species (M. milleri), but do not seem to disrupt gene flow in others (M. molossus). Lineages on distant islands undergo genetic bottlenecks more frequently than island lineages closer to the mainland, which have a greater exchange of haplotypes.  相似文献   

3.
The level of linkage disequilibrium (LD) is a major factor to determine DNA polymorphism pattern of a population and to construct high‐resolution maps useful in localizing and gene cloning of complicated traits. Here, we investigated LD level of three soybean populations with different genetic backgrounds and taxonomic status of G. gracilis by comparing the DNA polymorphism patterns of four high‐diversity single‐copy nuclear genes. A total of 152, 22, and 77 accessions of G. soja, G. gracilis, and G. max were observed. The results indicated that G. max retained only 75.3 (π) and 39% (θ) of the nucleotide polymorphism found in G. soja. Four gene loci evolved according to neutrality in both G. max and G. gracilis populations, and three gene loci evolved according to neutrality in G. soja population by Tajima's and Fu and Li's test. However, one gene locus deviated from neutrality by Fu and Li's test in the G. soja population. Further, medial level of LD (average r2 = 0.2426) was found in intragene in G. max and G. gracilis populations, but unexpected low level of LD (r2 ≤ 0.0539) was found in G. soja population. Significant genetic differentiation was detected between G. max and G. soja populations and also between G. max and G. gracilis populations; however, nonsignificant genetic differentiation was found between G. gracilis and G. soja populations. The results suggest that LD level depends on genetic background of soybean population, and implicit that G. gracilis should be regarded as the variant of G. soja, not as an independent species.  相似文献   

4.
Aim The purpose of this paper was to study the patterns of genetic variation, demographic history, haplotype relationships and potential location of diversity centres of two closely related species, Machilus thunbergii and Machilus kusanoi. Location The phylogeography of M. thunbergii and M. kusanoi was examined by sampling 110 and 106 individuals from 25 and 16 sampling sites, respectively, across their distributional range in Taiwan. Machilus thunbergii is distributed on the Asian mainland, South Korea, southern Japan, the Ryukyus, Taiwan and the Philippines, whereas M. kusanoi is endemic to Taiwan. These two species are closely related, and both are widely distributed in Taiwan but occupy different altitudinal zones and habitats. Methods The range‐wide variation of M. thunbergii and M. kusanoi in Taiwan was studied using chloroplast DNA (cpDNA) variations. A haplotype network was constructed with the computer program tcs . Nested clade analysis was conducted with the computer program ceodis , and various parameters of genetic diversity were calculated and neutrality tested by the computer program Dna SP. Population differentiation was estimated using the programs arlequin and hapstep . The contribution of the populations to gene diversity and to allelic richness was calculated using the software contrib . The level of divergence for each population from the remaining populations was calculated as the mean value of pairwise FST for each population against the rest of the populations. Results Extremely low levels of genetic differentiation were found for both species. This result suggested that these two species probably survived in multiple relict refugia with different population sizes throughout the island during low‐temperature periods of the Pleistocene. In addition, nested clade analysis (NCA) of cpDNA haplotypes indicated that restricted gene flow with isolation‐by‐distance characterized the recolonization after the Pleistocene by Tashueshan and Shiouhluan populations of M. thunbergii in the north‐central area west of the Central Mountain Range (CMR). In contrast, NCA analysis indicated that a major diversity centre on the southern tip of the island (Kending population) and contiguous range expansion characterized the recolonization by M. kusanoi of northern areas along the east side of the CMR. The major diversity centres found for the two species examined were further supported by the results of the mean FST for individual populations in comparison with other populations, and of the contribution of the divergence component to the total diversity. Main conclusions This research supports the multiple relict refugia hypothesis for both species investigated. Populations of M. thunbergii at Shiouhluan and Tashueshan in the north‐central area west of the CMR represent a diversity centre currently expanding its size. A diversity centre at the southern‐edge population of M. kusanoi, and a contiguous range expansion from Kending, were found. These results indicate that the M. thunbergii populations at Tashueshan and Shiouhluan and the M. kusanoi population at Kending, and even Soukar, are evolutionarily significant units for conservation programmes.  相似文献   

5.
Little is known about the natural history, biology, and population genetic structure of the Hardhead Silverside, Atherinomorus stipes, a small schooling fish found around islands throughout the Caribbean. Our field observations of A. stipes in the cays of Belize and the Florida Keys found that populations tend to be in close association with the shoreline in mangrove habitats. Due to this potential island‐based population structuring, A. stipes represents an ideal system to examine questions about gene flow and isolation by distance at different geographic scales. For this study, the mitochondrial gene nd2 was amplified from 394 individuals collected from seven different Belizean Cays (N = 175) and eight different Floridian Keys (N = 219). Results show surprisingly high haplotype diversity both within and between island‐groups, as well as a high prevalence of unique haplotypes within each island population. The results are consistent with models that require gene flow among populations as well as in situ evolution of rare haplotypes. There was no evidence for an isolation by distance model. The nd2 gene tree consists of two well‐supported monophyletic groups: a Belizean‐type clade and a Floridian‐type clade, indicating potential species‐level differentiation.  相似文献   

6.
1. Previous studies have identified lowland areas as barriers to gene flow (dispersal) between distinct mountain ranges in montane species of aquatic insects. In this study, we investigated the population genetic structure of two closely related Atalophlebia (mayfly) species inhabiting lowland areas of south‐east Queensland, Australia, with the expectation of widespread gene flow throughout the low‐altitude environment and associated homogeneous genetic structure. 2. In particular, we asked whether species with lower‐altitude distributions demonstrate greater spatial distribution of mtDNA (COI) alleles than the upland species studied previously. This pattern would be expected if good dispersal ability is associated with population persistence in these drought‐prone habitats. 3. The two species demonstrated contrasting genetic population structure. Atalophlebia sp. AV13 D revealed strong population structure, with populations on each side of the low‐altitude area isolated from each other for a long time (c.350 kya), and the presence of an isolation‐by‐distance pattern over relatively small geographical distances (<40 km). In contrast, Atalophlebia sp. AV13 A was panmictic at the scale investigated (≤160 km), with no history of past population fragmentation. 4. Examination of sample distribution along the altitudinal gradient reveals that Atalophlebia sp. AV13 D may have a more upland distribution (associated with greater habitat stability) than previously supposed, while Atalophlebia sp. AV13 A inhabits more xeric lowland areas, where freshwater habitats are less stable. We consequently hypothesise that these contrasting genetic population structures result from differences in habitat stability along the altitudinal gradient, only species with good dispersal ability being able to persist in unstable habitats. These findings may be applicable to other regions of the globe where habitat instability is associated with altitudinal gradients.  相似文献   

7.
Environmental gradients have been shown to disrupt gene flow in marine species, yet their influence in structuring populations at depth remains poorly understood. The Cape hakes (Merluccius paradoxus and M. capensis) are demersal species co‐occurring in the Benguela Current system, where decades of intense fishing resulted in severely depleted stocks in the past. Previous studies identified conflicting mtDNA genetic substructuring patterns and thus contrasting evolutionary trajectories for both species. Using 10 microsatellite loci, the control region of mtDNA and employing a seascape genetics approach, we investigated genetic connectivity and the impact of prolonged exploitation in the two species, which are characterized by different patterns of fishing pressure. Three consecutive years were sampled covering the entire distribution (N = 2100 fishes). Despite large estimated population sizes, both species exhibited low levels of contemporary genetic diversity (0.581 < HE < 0.692), implying that fishing has had a significant impact on their genetic composition and evolutionary trajectories. Further, for M. paradoxus, significant temporal, but not spatial, divergence points to the presence of genetic chaotic patchiness. In contrast, M. capensis exhibited a clear latitudinal cline in genetic differentiation between Namibia and South Africa (FST = 0.063, P < 0.05), with low (0.2% per generation) estimates of contemporary gene flow. Seascape analyses reveal an association with bathymetry and upwelling events, suggesting that adaptation to local environmental conditions may drive genetic differentiation in M. capensis. Importantly, our results highlight the need for temporal sampling in disentangling the complex factors that impact population divergence in marine fishes.  相似文献   

8.
1. The highly variable hydrology of dryland rivers has important implications for population dynamics in these systems. In western Queensland, fluctuations in sub‐population size are likely to lead to local bottlenecks and extinctions, increasing the need for connectivity and gene flow to maintain population viability. 2. Using microsatellite markers, we explored evidence for this metapopulation structure in two species of freshwater fish (Maquaria ambigua and Tandanus tandanus) and one crustacean (Macrobrachium australiense) in a sub‐catchment of the upper Murray–Darling Basin, Australia. 3. Overall, we found very weak genetic structure for all three species. Two species (M. ambigua and M. australiense) showed some significant genetic structure that did not correlate with geographic distance. However, decomposed pairwise regression analysis revealed evidence for intense genetic drift at the waterhole scale, suggesting that local bottlenecks are driving what little genetic structure does exist for these species. 4. The results identify the local impact of bottlenecks on genetic diversity, but highlight the importance of gene flow in maintaining population viability in these highly variable systems. As the impacts of bottlenecks are likely to be tempered by gene flow, it is suggested that the maintenance of connectivity is of paramount importance in this dryland system.  相似文献   

9.
Genomewide markers enable us to study genetic differentiation within a species and the factors underlying it at a much higher resolution than before, which advances our understanding of adaptation in organisms. We investigated genomic divergence in Metrosideros polymorpha, a woody species that occupies a wide range of ecological habitats across the Hawaiian Islands and shows remarkable phenotypic variation. Using 1659 single nucleotide polymorphism (SNP) markers annotated with the genome assembly, we examined the population genetic structure and demographic history of nine populations across five elevations and two ages of substrates on Mauna Loa, the island of Hawaii. The nine populations were differentiated into two genetic clusters distributed on the lower and higher elevations and were largely admixed on the middle elevation. Demographic modelling revealed that the two genetic clusters have been maintained in the face of gene flow, and the effective population size of the high‐altitude cluster was much smaller. A FST‐based outlier search among the 1659 SNPs revealed that 34 SNPs (2.05%) were likely to be under divergent selection and the allele frequencies of 21 of them were associated with environmental changes along elevations, such as temperature and precipitation. This study shows a genomic mosaic of M. polymorpha, in which contrasting divergence patterns were found. While most genomic polymorphisms were shared among populations, a small fraction of the genome was significantly differentiated between populations in diverse environments and could be responsible for the dramatic adaptation to a wide range of environments.  相似文献   

10.
Parrotia subaequalis (Hamamelidaceae) is a Tertiary relic species endemic in eastern China. We used inter‐simple sequence repeat (ISSR) markers to access genetic diversity and population genetic structure in natural five populations of P. subaequalis. The levels of genetic diversity were higher at species level (= 0.2031) but lower at population level (= 0.1096). The higher genetic diversity at species levels might be attributed to the accumulation of distinctive genotypes which adapted to the different habitats after Quaternary glaciations. Meanwhile, founder effects on the early stage, and subsequent bottleneck of population regeneration due to its biological characteristics, environmental features, and human activities, seemed to explain the low population levels of genetic diversity. The hierarchical AMOVA revealed high levels (42.60%) of among‐population genetic differentiation, which was in congruence with the high levels of Nei's genetic differentiation index (GST = 0.4629) and limited gene flow (Nm = 0.5801) among the studied populations. Mantel test showed a significant isolation‐by‐distance, indicating that geographic isolation has a significant effect on genetic structure in this species. Unweighted pair‐group method with arithmetic average clustering, PCoA, and Bayesian analyses uniformly recovered groups that matched the geographical distribution of this species. In particular, our results suggest that Yangtze River has served as a natural barrier to gene flow between populations occurred on both riversides. Concerning the management of P. subaequalis, the high genetic differentiation among populations indicates that preserving all five natural populations in situ and collecting enough individuals from these populations for ex situ conservation are necessary.  相似文献   

11.
12.
Resolving the genetic population structure of species inhabiting pristine, high latitude ecosystems can provide novel insights into the post‐glacial, evolutionary processes shaping the distribution of contemporary genetic variation. In this study, we assayed genetic variation in lake trout (Salvelinus namaycush) from Great Bear Lake (GBL), NT and one population outside of this lake (Sandy Lake, NT) at 11 microsatellite loci and the mtDNA control region (d‐loop). Overall, population subdivision was low, but significant (global FST θ = 0.025), and pairwise comparisons indicated that significance was heavily influenced by comparisons between GBL localities and Sandy Lake. Our data indicate that there is no obvious genetic structure among the various basins within GBL (global FST = 0.002) despite the large geographic distances between sampling areas. We found evidence of low levels of contemporary gene flow among arms within GBL, but not between Sandy Lake and GBL. Coalescent analyses suggested that some historical gene flow occurred among arms within GBL and between GBL and Sandy Lake. It appears, therefore, that contemporary (ongoing dispersal and gene flow) and historical (historical gene flow and large founding and present‐day effective population sizes) factors contribute to the lack of neutral genetic structure in GBL. Overall, our results illustrate the importance of history (e.g., post‐glacial colonization) and contemporary dispersal ecology in shaping genetic population structure of Arctic faunas and provide a better understanding of the evolutionary ecology of long‐lived salmonids in pristine, interconnected habitats.  相似文献   

13.
Identifying processes that promote or limit gene flow can help define the ecological and evolutionary history of a species. Furthermore, defining those factors that make up “species boundaries” can provide a definition of the independent evolutionary trajectories of related taxa. For many species, the historic processes that account for their distribution of genetic variation remain unresolved. In this study, we examine the geographic distribution of genetic diversity for two species of Louisiana Irises, Iris brevicaulis and Iris fulva. Specifically, we asked how populations are structured and if population structure coincides with potential barriers to gene flow. We also asked whether there is evidence of hybridization between these two species outside Louisiana hybrid zones. We used a genotyping‐by‐sequencing approach and sampled a large number of single nucleotide polymorphisms across these species' genomes. Two different population assignment methods were used to resolve population structure in I. brevicaulis; however, there was considerably less population structure in I. fulva. We used a species tree approach to infer phylogenies both within and between populations and species. For I. brevicaulis, the geography of the collection locality was reflected in the phylogeny. The I. fulva phylogeny reflected much less structure than detected for I. brevicaulis. Lastly, combining both species into a phylogenetic analysis resolved two of six populations of I. brevicaulis that shared alleles with I. fulva. Taken together, our results suggest major differences in the level and pattern of connectivity among populations of these two Louisiana Iris species.  相似文献   

14.
Geological history of oceanic islands can have a profound effect on the evolutionary history of insular flora, especially in complex islands such as Tenerife in the Canary Islands. Tenerife results from the secondary connection of three paleo‐islands by a central volcano, and other geological events that further shaped it. This geological history has been shown to influence the phylogenetic history of several taxa, including genus Micromeria (Lamiaceae). Screening 15 microsatellite markers in 289 individuals representing the eight species of Micromeria present in Tenerife, this study aims to assess the genetic diversity and structure of these species and its relation with the geological events on the island. In addition, we evaluate the extent of hybridization among species and discuss its influence on the speciation process. We found that the species restricted to the paleo‐islands present lower levels of genetic diversity but the highest levels of genetic differentiation suggesting that their ranges might have contracted over time. The two most widespread species in the island, M. hyssopifolia and M. varia, present the highest genetic diversity levels and a genetic structure that seems correlated with the geological composition of the island. Samples from M. hyssopifolia from the oldest paleo‐island, Adeje, appear as distinct while samples from M. varia segregate into two main clusters corresponding to the paleo‐islands of Anaga and Teno. Evidence of hybridization and intraspecific migration between species was found. We argue that species boundaries would be retained despite hybridization in response to the habitat's specific conditions causing postzygotic isolation and preserving morphological differentiation.  相似文献   

15.
Patterns of genetic differentiation among taxa at early stages of divergence provide an opportunity to make inferences about the history of speciation. Here, we conduct a survey of DNA‐sequence polymorphism and divergence at loci on the autosomes, X chromosome, Y chromosome and mitochondrial DNA in samples of Mus domesticus, M. musculus and M. castaneus. We analyzed our data under a divergence with gene flow model and estimate that the effective population size of M. castaneus is 200 000–400 000, of M. domesticus is 100 000–200 000 and of M. musculus is 60 000–120 000. These data also suggest that these species started to diverge approximately 500 000 years ago. Consistent with this recent divergence, we observed considerable variation in the genealogical patterns among loci. For some loci, all alleles within each species formed a monophyletic group, while at other loci, species were intermingled on the phylogeny of alleles. This intermingling probably reflects both incomplete lineage sorting and gene flow after divergence. Likelihood ratio tests rejected a strict allopatric model with no gene flow in comparisons between each pair of species. Gene flow was asymmetric: no gene flow was detected into M. domesticus, while significant gene flow was detected into both M. castaneus and M. musculus. Finally, most of the gene flow occurred at autosomal loci, resulting in a significantly higher ratio of fixed differences to polymorphisms at the X and Y chromosomes relative to autosomes in some comparisons, or just the X chromosome in others, emphasizing the important role of the sex chromosomes in general and the X chromosome in particular in speciation.  相似文献   

16.
Identifying drivers of dispersal limitation and genetic differentiation is a key goal in biogeography. We examine patterns of population connectivity and genetic diversity using restriction site‐associated DNA sequencing (RADseq) in two bumble bee species, Bombus vosnesenskii and Bombus bifarius, across latitude and altitude in mountain ranges from California, Oregon and Washington, U.S.A. Bombus vosnesenskii, which occurs across a broader elevational range at most latitudes, exhibits little population structure while B. bifarius, which occupies a relatively narrow higher elevation niche across most latitudes, exhibits much stronger population differentiation, although gene flow in both species is best explained by isolation with environmental niche resistance. A relationship between elevational habitat breadth and genetic diversity is also apparent, with B. vosnesenskii exhibiting relatively consistent levels of genetic diversity across its range, while B. bifarius has reduced genetic diversity at low latitudes, where it is restricted to high‐elevation habitat. The results of this study highlight the importance of the intersect between elevational range and habitat suitability in influencing population connectivity and suggest that future climate warming will have a fragmenting effect even on populations that are presently well connected, as they track their thermal niches upward in montane systems.  相似文献   

17.
Scanning genomes for loci with high levels of population differentiation has become a standard of population genetics. FST outlier loci are most often interpreted as signatures of local selection, but outliers might arise for many other reasons too often left unexplored. Here, we tried to identify further the history and genetic basis underlying strong differentiation at FST outlier loci in a marine mussel. A genome scan of genetic differentiation has been conducted between Atlantic and Mediterranean populations of Mytilus galloprovincialis. The differentiation was low overall (FST = 0.03), but seven loci (2%) were strong FST outliers. We then analysed DNA sequence polymorphism at two outlier loci. The genetic structure proved to be the consequence of differential introgression of alleles from the sister‐hybridizing species Mytilus edulis. Surprisingly, the Mediterranean population was the most introgressed at these two loci, although the contact zone between the two species is nowadays localized along the Atlantic coasts of France and the British Isles. A historical contact between M. edulis and Mediterranean M. galloprovincialis should have happened during glacial periods. It proved difficult to disentangle two hypotheses: (i) introgression was adaptive, implying edulis alleles have been favoured in Mediterranean populations, or (ii) the genetic architecture of the barrier to edulis gene flow is different between the two M. galloprovincialis backgrounds. Five of the seven outliers between M. galloprovincialis populations were also outliers between M. edulis and Atlantic M. galloprovincialis, which would support the latter hypothesis. Differential introgression across semi‐permeable barriers to gene flow is a neglected scenario to interpret outlying loci that may prove more widespread than anticipated.  相似文献   

18.
19.
The Red‐headed Wood Pigeon Columba janthina nitens is endemic to the Ogasawara Islands, an oceanic island chain located 1000 km south of the main islands of Japan. The subspecies is at high risk of extinction because of its small population size and restricted habitat range. We undertook genetic analyses of this pigeon using sequences of a portion of the mitochondrial control region and five microsatellite markers to estimate the genetic characteristics of two wild populations from the Bonin and Volcano Islands, as well as one captive breeding population. The genetic diversity of the wild individuals was exceptionally low in both the mitochondria (nucleotide diversity = 0.00105) and at the microsatellite (3.2 alleles per locus and HE = 0.12) loci. Higher numbers of microsatellite genotypes were observed in the Volcano Islands population than in the Bonin Islands population, which may be because of the relatively low impact of human disturbance. The most common mitochondrial haplotypes and microsatellite alleles observed in the two wild populations were completely fixed in the captive population. Our results suggest that the genetic diversity of the captive population needs to be increased. However, introduction of a wild individual into a captive population can lead to a decreased genetic diversity in the wild population and therefore should be done with caution. The genetic differentiation between the Bonin and the Volcano island groups was low, and the populations of the two island groups should be regarded as a single evolutionarily significant unit. However, special consideration is required for habitat conservation in the Volcano Islands, which may be functioning as a sanctuary for the Red‐headed Wood Pigeon. For the long‐term conservation of threatened bird species that live on remote oceanic islands, determination of management units considering gene flow caused by their flying capacity and maintenance of genetically suitable wild and captive populations are essential.  相似文献   

20.
Estimating population connectivity and species' abilities to disperse across the landscape is crucial for understanding the long‐term persistence of species in changing environments. Surprisingly, few landscape genetic studies focused on tropical regions despite the alarming extinction rates within these ecosystems. Here, we compared the influence of landscape features on the distribution of genetic variation of an Afromontane frog, Amietia wittei, with that of its more broadly distributed lowland congener, Amietia angolensis, on Mt. Kilimanjaro, Tanzania. We predicted high gene flow in the montane species with movements enhanced through terrestrial habitats of the continuous rainforest. In contrast, dispersal might be restricted to aquatic corridors and reduced by anthropogenic disturbance in the lowland species. We found high gene flow in A. wittei relative to other montane amphibians. Nonetheless, gene flow was lower than in the lowland species which showed little population structure. Least‐cost path analysis suggested that dispersal is facilitated by stream networks in both species, but different landscape features were identified to influence connectivity among populations. Contrary to a previous study, gene flow in the lowland species was negatively correlated with the presence of human settlements. Also, genetic subdivision in A. wittei did not coincide with specific physical barriers as in other landscape genetic studies, suggesting that factors other than topography may contribute to population divergence. Overall, these results highlight the importance of a comparative landscape genetic approach for assessing the influence of the landscape matrix on population connectivity, particularly because nonintuitive results can alter the course of conservation and management.  相似文献   

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

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