首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 78 毫秒
1.
The hermaphroditic nematode Pristionchus pacificus is a model organism with a range of fully developed genetic tools. The species is globally widespread and highly diverse genetically, consisting of four major independent lineages (lineages A, B, C, and D). Despite its young age (~2.1 Ma), volcanic La Réunion Island harbors all four lineages. Ecological and population genetic research studies suggest that this diversity is due to repeated independent island colonizations by P. pacificus. Here, we use model‐based statistical methods to rigorously test hypotheses regarding the evolutionary history of P. pacificus. First, we employ divergence analyses to date diversification events among the four “world” lineages. Next, we examine demographic properties of a subset of four populations (“a”, “b”, “c”, and “d”), present on La Réunion Island. Finally, we use the results of the divergence and demographic analyses to inform a modeling‐based approximate Bayesian computation (ABC) approach, where we test hypotheses about the order and timing of establishment of the Réunion populations. Our dating estimates place the recent common ancestor of P. pacificus lineages at nearly 500,000 generations past. Our demographic analysis supports recent (<150,000 generations) spatial expansion for the island populations, and our ABC approach supports c>a>b>d as the most likely colonization order of the island populations. Collectively, our study comprehensively improves previous inferences about the evolutionary history of P. pacificus.  相似文献   

2.

Background  

The nematode Pristionchus pacificus has originally been developed as a satellite organism for comparison to Caenorhabditis elegans. A 10X coverage of the whole genome of P. pacificus is available, making P. pacificus the first non- Caenorhabditis nematode with a fully sequenced genome. The macroevolutionary comparison between P. pacificus and C. elegans has been complemented by microevolutionary studies of closely related strains and species within the genus Pristionchus. In addition, new understanding of the biology of Pristionchus from field studies, demonstrating a close association with various scarab beetles and the Colorado potato beetle, supports consideration of this nematode in studies of ecosystems. In the course of field studies on four continents more than 1,200 isolates were established from 15,000 beetle specimens representing 18 Pristionchus species. Two remarkable features of the Pristionchus – beetle association are the high species specificity of the interaction and the interception of the beetle's sex communication system for host recognition by the nematodes, as suggested by chemotaxis studies. Evolutionary interpretations of differences in developmental, behavioral and ecological patterns require a phylogenetic framework of the genus Pristionchus.  相似文献   

3.
Pristionchus pacificus has been established as a nematode model system in evolutionary developmental biology and evolutionary ecology. Field studies in North and South America, Asia, Africa and Europe indicated that nematodes of the genus Pristionchus live in association with scarab beetles. Here, we describe the first account of soil‐ and beetle‐associated nematodes on an island setting by investigating the island of Réunion in the Indian Ocean. Réunion has high numbers of endemic insects and is one among several attractive islands for biodiversity studies. Being of volcanic origin, Réunion is 2–3 million years old, making it the youngest of the Mascareigne islands. We show that beetle‐ and soil‐derived nematodes on Réunion are nearly exclusively hermaphroditic, suggesting that selfing is favoured over gonochorism (outcrossing) during island colonization. Among members of four nematode genera observed on Réunion, Pristionchus pacificus was the most prevalent species. A total of 76 isolates, in association with five different scarab beetles, has been obtained for this cosmopolitan nematode. A detailed mitochondrial haplotype analysis indicates that the Réunion isolates of P. pacificus cover all four worldwide clades of the species. This extraordinary haplotype diversity suggests multiple independent invasions, most likely in association with different scarab beetles. Together, we establish Réunion as a case study for nematode island biogeography, in which the analysis of nematode population genetics and population dynamics can provide insight into evolutionary and ecological processes. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 170–179.  相似文献   

4.
Phylogeographical studies have shown that some shallow‐water marine organisms, such as certain coral reef fishes, lack spatial population structure at oceanic scales, despite vast distances of pelagic habitat between reefs and other dispersal barriers. However, whether these dispersive widespread taxa constitute long‐term panmictic populations across their species ranges remains unknown. Conventional phylogeographical inferences frequently fail to distinguish between long‐term panmixia and metapopulations connected by gene flow. Moreover, marine organisms have notoriously large effective population sizes that confound population structure detection. Therefore, at what spatial scale marine populations experience independent evolutionary trajectories and ultimately species divergence is still unclear. Here, we present a phylogeographical study of a cosmopolitan Indo‐Pacific coral reef fish Naso hexacanthus and its sister species Naso caesius, using two mtDNA and two nDNA markers. The purpose of this study was two‐fold: first, to test for broad‐scale panmixia in N. hexacanthus by fitting the data to various phylogeographical models within a Bayesian statistical framework, and second, to explore patterns of genetic divergence between the two broadly sympatric species. We report that N. hexacanthus shows little population structure across the Indo‐Pacific and a range‐wide, long‐term panmictic population model best fit the data. Hence, this species presently comprises a single evolutionary unit across much of the tropical Indian and Pacific Oceans. Naso hexacanthus and N. caesius were not reciprocally monophyletic in the mtDNA markers but showed varying degrees of population level divergence in the two nuclear introns. Overall, patterns are consistent with secondary introgression following a period of isolation, which may be attributed to oceanographic conditions of the mid to late Pleistocene, when these two species appear to have diverged.  相似文献   

5.
Characterizing the genetic basis of among‐species variation in lifespan is a major goal of evolutionary gerontology research, but the very feature that defines separate species – the inability to interbreed – makes achieving this goal impractical, if not impossible, for most taxa. Pristionchus nematodes provide an intriguing system for tackling this problem, as female lifespan varies among species that can be crossed to form viable (although infertile) hybrids. By conducting reciprocal crosses among three species – two dioecious (long‐lived Pristionchus exspectatus and short‐lived Pristionchus arcanus) and one androdioecious (short‐lived Pristionchus pacificus) – we found that female lifespan was long for all hybrids, consistent with the hypothesis that the relatively short lifespans seen for P. pacificus hermaphrodites and P. arcanus females are caused by independent, recessive alleles that are masked in hybrid genomes. Cross‐direction had a small effect on survivorship for crosses involving P. exspectatus, indicating that nuclear–mitochondrial interactions may also influence Pristionchus longevity. Our findings suggest that long lifespan in P. exspectatus reflects the realization of an ancestral potential for extended longevity in the P. pacificus species complex. This work demonstrates the utility of interspecific hybrids for ageing research and provides a foundation for future work on the genetic architecture of interspecific lifespan variation.  相似文献   

6.
7.
Ribosomal RNA genes have long been a favoured locus in phylogenetic and metabarcoding studies. Within a genome, rRNA loci are organized as tandem repeated arrays and the copies are homogenized through the process of concerted evolution. However, some level of rRNA variation (intragenomic polymorphism) is known to persist and be maintained in the genomes of many species. In nematode worms, the extent of rRNA polymorphism (RP) across species and the evolutionary and life history factors that contribute to the maintenance of intragenomic RP is largely unknown. Here, we present an extensive analysis across 30 terrestrial nematode species representing a range of free‐living and parasitic taxa isolated worldwide. Our results indicate that RP is common and widespread, ribosome function appears to be maintained despite mutational changes, and intragenomic variants are stable in the genome and neutrally evolving. However, levels of variation were varied widely across rRNA locus and species, with some taxa observed to lack RP entirely. Higher levels of RP were significantly correlated with shorter generation time and high reproductive rates, and population‐level factors may play a role in the geographic and phylogenetic structuring of rRNA variants observed in genera such as Rotylenchulus and Pratylenchus. Although RP did not dramatically impact the clustering and recovery of taxa in mock metabarcoding analyses, the present study has significant implications for global biodiversity estimates of nematode species derived from environmental rRNA amplicon studies, as well as our understanding of the evolutionary and ecological factors shaping genetic diversity across the nematode Tree of Life.  相似文献   

8.
Social insects are the target of numerous pathogens. This is because the high density of closely‐related individuals frequently interacting with each other enhances the transmission and establishment of pathogens. This high selective pressure results in the rapid evolution of immune genes, which might be counteracted by a reduced effective population size (Ne) lowering the effectiveness of selection. We tested the effect of Ne on the evolutionary rate of an important immune gene for the antimicrobial peptide Hymenoptaecin in two common central European bumblebee species: Bombus terrestris and Bombus lapidarius. Both species are similar in their biology and are expected to be under similar selective pressures because pathogen prevalence does not differ between species. However, previous studies indicated a higher Ne in B. terrestris compared to B. lapidarius. We found high intraspecific variability in the coding sequence but low variability for silent polymorphisms in B. lapidarius. Estimates of long‐ and short‐term Ne were three‐ to four‐fold higher Ne in B. terrestris, although the species did not differ in census population sizes. The difference in Ne might result in less efficient selection and suboptimal adaptation of immune genes (e.g. hymenoptaecin) in B. lapidarius, and thus this species might become less resistant and more tolerant, turning into a superspreader of diseases.  相似文献   

9.
Pristionchus pacificus, recently established as a model organism in evolutionary biology, is a cosmopolitan nematode that has a necromenic association with scarab beetles. The diverse array of host beetle species and habitat types occupied by P. pacificus make it a good model for investigating local adaptation to novel environments. Presence of P. pacificus on La Réunion Island, a young volcanic island with a dynamic geological history and a wide variety of ecozones, facilitates such investigation in an island biogeographic setting. Microsatellite data from 20 markers and 223 strains and mitochondrial sequence data from 272 strains reveal rich genetic diversity among La Réunion P. pacificus isolates, shaped by differentially timed introductions from diverse sources and in association with different beetle species. Distinctions between volcanic zones and between arid western and wet eastern climatic zones have likely limited westward dispersal of recently colonized lineages and maintained a genetic distinction between eastern and western clades. The highly selfing lifestyle of P. pacificus contributes to the strong fine‐scale population structure detected, with each beetle host harbouring strongly differentiated assemblages of strains. Periodic out‐crossing generates admixture between genetically diverse lineages, creating a diverse array of allelic combinations likely to increase the evolutionary potential of the species and facilitate adaptation to local environments and beetle hosts.  相似文献   

10.
Parasitoid wasps represent a large proportion of hymenopteran species. They have complex evolutionary histories and are important biocontrol agents. To advance parasitoid research, a combination of Illumina short‐read, PacBio long‐read and Hi‐C scaffolding technologies was used to develop a high‐quality chromosome‐level genome assembly for Pteromalus puparum, which is an important pupal endoparasitoid of caterpillar pests. The chromosome‐level assembly has aided in studies of venom and detoxification genes. The assembled genome size is 338 Mb with a contig N50 of 38.7 kb and a scaffold N50 of 1.16 Mb. Hi‐C analysis assembled scaffolds onto five chromosomes and raised the scaffold N50 to 65.8 Mb, with more than 96% of assembled bases located on chromosomes. Gene annotation was assisted by RNA sequencing for the two sexes and four different life stages. Analysis detected 98% of the BUSCO (Benchmarking Universal Single‐Copy Orthologs) gene set, supporting a high‐quality assembly and annotation. In total, 40.1% (135.6 Mb) of the assembly is composed of repetitive sequences, and 14,946 protein‐coding genes were identified. Although venom genes play important roles in parasitoid biology, their spatial distribution on chromosomes was poorly understood. Mapping has revealed venom gene tandem arrays for serine proteases, pancreatic lipase‐related proteins and kynurenine–oxoglutarate transaminases, which have amplified in the P. puparum lineage after divergence from its common ancestor with Nasonia vitripennis. In addition, there is a large expansion of P450 genes in P. puparum. These examples illustrate how chromosome‐level genome assembly can provide a valuable resource for molecular, evolutionary and biocontrol studies of parasitoid wasps.  相似文献   

11.
The nematode Caenorhabditis elegans is a central laboratory model system in almost all biological disciplines, yet its natural life history and population biology are largely unexplored. Such information is essential for in‐depth understanding of the nematode's biology because its natural ecology provides the context, in which its traits and the underlying molecular mechanisms evolved. We characterized natural phenotypic and genetic variation among North German C. elegans isolates. We used the unique opportunity to compare samples collected 10 years apart from the same compost heap and additionally included recent samples for this and a second site, collected across a 1.5‐year period. Our analysis revealed significant population genetic differentiation between locations, across the 10‐year time period, but for only one location a trend across the shorter time frame. Significant variation was similarly found for phenotypic traits of likely importance in nature, such as choice behavior and population growth in the presence of pathogens or naturally associated bacteria. Phenotypic variation was significantly influenced by C. elegans genotype, time of isolation, and sampling site. The here studied C. elegans isolates may provide a valuable, genetically variable resource for future dissection of naturally relevant gene functions.  相似文献   

12.
By constraining gene flow, group living and natal philopatry can result in fine‐scale genetic structure. Although the genetic structure of some group‐living lizards has been characterised, studies are few compared with those for group‐living bird and mammal species. The Egerniinae group of lizards exhibits a high diversity of social structures, making it a useful group for comparative studies of genetic structure across a broader range of social taxa. A well‐studied member of Egerniinae is Egernia stokesii, a lizard that forms long‐term pair bonds and stable social groups and exhibits natal philopatry and limited dispersal. Evidence exists for consistent E. stokesii social structure across seven close but disconnected rocky outcrops within a 40 × 10 km area. We used summary statistics, analysis of molecular variance, Bayesian clustering, and discriminant analysis of principal components to assess if E. stokesii exhibit a consistent pattern of fine‐scale genetic structure across the same seven outcrops. Due to E. stokesii social structure and constrained dispersal, we predicted significant genetic structuring – based on microsatellite markers – among outcrops. We found significant fine‐scale genetic structuring and evidence for two genetic clusters. We discuss features of E. stokesii biology and ecology that may explain our findings. Some rocky outcrops, and some social groups, contained lizards from both genetic clusters. An examination of the composition of mixed cluster social groups did not detect any notable patterns. Therefore, further work is necessary to identify how the observed patterns may have arisen. Future investigations in E. stokesii and other group‐living lizard species are likely to contribute greatly to our understanding of the genetic consequences of group living.  相似文献   

13.
Parasites and hosts interact across both micro‐ and macroevolutionary scales where congruence among their phylogeographic and phylogenetic structures may be observed. Within southern Africa, the four‐striped mouse genus, Rhabdomys, is parasitized by the ectoparasitic sucking louse, Polyplax arvicanthis. Molecular data recently suggested the presence of two cryptic species within P. arvicanthis that are sympatrically distributed across the distributions of four putative Rhabdomys species. We tested the hypotheses of phylogeographic congruence and cophylogeny among the two parasite lineages and the four host taxa, utilizing mitochondrial and nuclear sequence data. Despite the documented host‐specificity of P. arvicanthis, limited phylogeographic correspondence and nonsignificant cophylogeny was observed. Instead, the parasite–host evolutionary history is characterized by limited codivergence and several duplication, sorting and host‐switching events. Despite the elevated mutational rates found for P. arvicanthis, the spatial genetic structure was not more pronounced in the parasite lineages compared with the hosts. These findings may be partly attributed to larger effective population sizes of the parasite lineages, the vagility and social behaviour of Rhabdomys, and the lack of host‐specificity observed in areas of host sympatry. Further, the patterns of genetic divergence within parasite and host lineages may also be largely attributed to historical biogeographic changes (expansion‐contraction cycles). It is thus evident that the association between P. arvicanthis and Rhabdomys has been shaped by the synergistic effects of parasite traits, host‐related factors and biogeography over evolutionary time.  相似文献   

14.
The leopard coral grouper, Plectropomus leopardus, belonging to the family Epinephelinae, is a carnivorous coral reef fish widely distributed in tropical and subtropical waters of the Indo‐Pacific. Due to its appealing body appearance and delicious taste, P. leopardus has become a popular commercial fish for aquaculture in many countries. However, the lack of genomic and molecular resources for P. leopardus has hindered study of its biology and genomic breeding programmes. Here we report the de novo sequencing and assembly of the P. leopardus genome using a combination of 10 × Genomics, high‐throughput chromosome conformation capture (Hi‐C) and PacBio long‐read sequencing technologies. The genome assembly has a total length of 881.55 Mb with a scaffold N50 of 34.15 Mb, consisting of 24 pseudochromosome scaffolds. busco analysis showed that 97.2% of the conserved single‐copy genes were retrieved, indicating the assembly was almost entire. We predicted 25,248 protein‐coding genes, among which 96.5% were functionally annotated. Comparative genomic analyses revealed that gene family expansions in P. leopardus were associated with immune‐related pathways. In addition, we identified 5,178,453 single nucleotide polymorphisms based on genome resequencing of 54 individuals. The P. leopardus genome and genomic variation data provide valuable genomic resources for studies of its genetics, evolution and biology. In particular, it is expected to benefit the development of genomic breeding programmes in the farming industry.  相似文献   

15.
The green anole (Anolis carolinensis) is a lizard widespread throughout the southeastern United States and is a model organism for the study of reproductive behavior, physiology, neural biology, and genomics. Previous phylogeographic studies of A. carolinensis using mitochondrial DNA and small numbers of nuclear loci identified conflicting and poorly supported relationships among geographically structured clades; these inconsistencies preclude confident use of A. carolinensis evolutionary history in association with morphological, physiological, or reproductive biology studies among sampling localities and necessitate increased effort to resolve evolutionary relationships among natural populations. Here, we used anchored hybrid enrichment of hundreds of genetic markers across the genome of A. carolinensis and identified five strongly supported phylogeographic groups. Using multiple analyses, we produced a fully resolved species tree, investigated relative support for each lineage across all gene trees, and identified mito‐nuclear discordance when comparing our results to previous studies. We found fixed differences in only one clade—southern Florida restricted to the Everglades region—while most polymorphisms were shared between lineages. The southern Florida group likely diverged from other populations during the Pliocene, with all other diversification during the Pleistocene. Multiple lines of support, including phylogenetic relationships, a latitudinal gradient in genetic diversity, and relatively more stable long‐term population sizes in southern phylogeographic groups, indicate that diversification in A. carolinensis occurred northward from southern Florida.  相似文献   

16.
All organisms must interact with their environment, responding in behavioral, chemical, and other ways to various stimuli throughout their life cycles. Characterizing traits that directly represent an organism's ability to sense and react to their environment provides useful insight into the evolution of life‐history strategies. One such trait for the nematode Pristionchus pacificus, chemosensation, is involved in navigation to beetle hosts. Essential for the survival of the nematode, chemosensory behavior may be subject to variation as nematodes discriminate among chemical cues to complete their life cycle. We examine this hypothesis using natural isolates of P. pacificus from La Réunion Island. We select strains from a variety of La Réunion beetle hosts and geographic locations and examine their chemoattraction response toward organic compounds, beetle washes, and live beetles. We find that nematodes show significant differences in their response to various chemicals and are able to chemotax to live beetles in a novel assay. Further, strains can discriminate among different cues, showing more similar responses toward beetle washes than to organic compounds in cluster analyses. However, we find that variance in chemoattraction response is not significantly associated with temperature, location, or beetle host. Rather, strains show a more concerted response toward compounds they most likely directly encounter in the wild. We suggest that divergence in odor‐guided behavior in P. pacificus may therefore have an important ecological component.  相似文献   

17.
Soybean cyst nematode (SCN, Heterodera glycines) is a major pest of soybean that is spreading across major soybean production regions worldwide. Increased SCN virulence has recently been observed in both the United States and China. However, no study has reported a genome assembly for H. glycines at the chromosome scale. Herein, the first chromosome‐level reference genome of X12, an unusual SCN race with high infection ability, is presented. Using whole‐genome shotgun (WGS) sequencing, Pacific Biosciences (PacBio) sequencing, Illumina paired‐end sequencing, 10X Genomics linked reads and high‐throughput chromatin conformation capture (Hi‐C) genome scaffolding techniques, a 141.01‐megabase (Mb) assembled genome was obtained with scaffold and contig N50 sizes of 16.27 Mb and 330.54 kilobases (kb), respectively. The assembly showed high integrity and quality, with over 90% of Illumina reads mapped to the genome. The assembly quality was evaluated using Core Eukaryotic Genes Mapping Approach and Benchmarking Universal Single‐Copy Orthologs. A total of 11,882 genes were predicted using de novo, homolog and RNAseq data generated from eggs, second‐stage juveniles (J2), third‐stage juveniles (J3) and fourth‐stage juveniles (J4) of X12, and 79.0% of homologous sequences were annotated in the genome. These high‐quality X12 genome data will provide valuable resources for research in a broad range of areas, including fundamental nematode biology, SCN–plant interactions and co‐evolution, and also contribute to the development of technology for overall SCN management.  相似文献   

18.
Pheromones are among the most important sexual signals used by organisms throughout the animal kingdom. However, few are identified in vertebrates, leaving the evolutionary mechanisms underlying vertebrate pheromones poorly understood. Pre‐existing biases in receivers' perceptual systems shape visual and auditory signalling systems, but studies on how receiver biases influence the evolution of pheromone communication remain sparse. The lamprey Petromyzon marinus uses a relatively well‐understood suite of pheromones and offers a unique opportunity to study the evolution of vertebrate pheromone communication. Previous studies indicate that male signalling with the mating pheromone 3‐keto petromyzonol sulphate (3kPZS) may exploit a nonsexual attraction to juvenile‐released 3kPZS that guides migration into productive rearing habitat. Here, we infer the distribution of male signalling with 3kPZS using a phylogenetic comparison comprising six of 10 genera and two of three families. Our results indicate that only P. marinus and Ichthyomyzon castaneus release 3kPZS at high rates. Olfactory and behavioural assays with P. marinus, I. castaneus and a subset of three other species that do not use 3kPZS as a sexual signal indicate that male signalling might have driven the evolution of female adaptations to detect 3kPZS with specific olfactory mechanisms and respond to 3kPZS with targeted attraction relevant during mate search. We postulate that 3kPZS communication evolved independently in I. castaneus and P. marinus, but cannot eliminate the alternative that other species lost 3kPZS communication. Regardless, our results represent a rare macroevolutionary investigation of a vertebrate pheromone and provide insight into the evolutionary mechanisms underlying pheromone communication.  相似文献   

19.
Understanding the earliest events in speciation remains a major challenge in evolutionary biology. Thus identifying species whose populations are beginning to diverge can provide useful systems to study the process of speciation. Drosophila aldrichi, a cactophilic fruit fly species with a broad distribution in North America, has long been assumed to be a single species owing to its morphological uniformity. While previous reports either of genetic divergence or reproductive isolation among different D. aldrichi strains have hinted at the existence of cryptic species, the evolutionary relationships of this species across its range have not been thoroughly investigated. Here we show that D. aldrichi actually is paraphyletic with respect to its closest relative, Drosophila wheeleri, and that divergent D. aldrichi lineages show complete hybrid male sterility when crossed. Our data support the interpretation that there are at least two species of D. aldrichi, making these flies particularly attractive for studies of speciation in an ecological and geographical context.  相似文献   

20.
Accurate representation of temperature sensitivity (Q10) of soil microbial activity across time is critical for projecting soil CO2 efflux. As microorganisms mediate soil carbon (C) loss via exo‐enzyme activity and respiration, we explore temperature sensitivities of microbial exo‐enzyme activity and respiratory CO2 loss across time and assess mechanisms associated with these potential changes in microbial temperature responses. We collected soils along a latitudinal boreal forest transect with different temperature regimes (long‐term timescale) and exposed these soils to laboratory temperature manipulations at 5, 15, and 25°C for 84 days (short‐term timescale). We quantified temperature sensitivity of microbial activity per g soil and per g microbial biomass at days 9, 34, 55, and 84, and determined bacterial and fungal community structure before the incubation and at days 9 and 84. All biomass‐specific rates exhibited temperature sensitivities resistant to change across short‐ and long‐term timescales (mean Q10 = 2.77 ± 0.25, 2.63 ± 0.26, 1.78 ± 0.26, 2.27 ± 0.25, 3.28 ± 0.44, 2.89 ± 0.55 for β‐glucosidase, N‐acetyl‐β‐d ‐glucosaminidase, leucine amino peptidase, acid phosphatase, cellobiohydrolase, and CO2 efflux, respectively). In contrast, temperature sensitivity of soil mass‐specific rates exhibited either resilience (the Q10 value changed and returned to the original value over time) or resistance to change. Regardless of the microbial flux responses, bacterial and fungal community structure was susceptible to change with temperature, significantly differing with short‐ and long‐term exposure to different temperature regimes. Our results highlight that temperature responses of microbial resource allocation to exo‐enzyme production and associated respiratory CO2 loss per unit biomass can remain invariant across time, and thus, that vulnerability of soil organic C stocks to rising temperatures may persist in the long term. Furthermore, resistant temperature sensitivities of biomass‐specific rates in spite of different community structures imply decoupling of community constituents and the temperature responses of soil microbial activities.  相似文献   

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

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