首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 312 毫秒
1.
Multiple disease resistance has important implications for plant fitness, given the selection pressure that many pathogens exert directly on natural plant populations and indirectly via crop improvement programs. Evidence of a locus conditioning resistance to multiple pathogens was found in bin 1.06 of the maize genome with the allele from inbred line “Tx303” conditioning quantitative resistance to northern leaf blight (NLB) and qualitative resistance to Stewart’s wilt. To dissect the genetic basis of resistance in this region and to refine candidate gene hypotheses, we mapped resistance to the two diseases. Both resistance phenotypes were localized to overlapping regions, with the Stewart’s wilt interval refined to a 95.9-kb segment containing three genes and the NLB interval to a 3.60-Mb segment containing 117 genes. Regions of the introgression showed little to no recombination, suggesting structural differences between the inbred lines Tx303 and “B73,” the parents of the fine-mapping population. We examined copy number variation across the region using next-generation sequencing data, and found large variation in read depth in Tx303 across the region relative to the reference genome of B73. In the fine-mapping region, association mapping for NLB implicated candidate genes, including a putative zinc finger and pan1. We tested mutant alleles and found that pan1 is a susceptibility gene for NLB and Stewart’s wilt. Our data strongly suggest that structural variation plays an important role in resistance conditioned by this region, and pan1, a gene conditioning susceptibility for NLB, may underlie the QTL.  相似文献   

2.
Pavlidis P  Metzler D  Stephan W 《Genetics》2012,192(1):225-239
We study the trajectory of an allele that affects a polygenic trait selected toward a phenotypic optimum. Furthermore, conditioning on this trajectory we analyze the effect of the selected mutation on linked neutral variation. We examine the well-characterized two-locus two-allele model but we also provide results for diallelic models with up to eight loci. First, when the optimum phenotype is that of the double heterozygote in a two-locus model, and there is no dominance or epistasis of effects on the trait, the trajectories of selected mutations rarely reach fixation; instead, a polymorphic equilibrium at both loci is approached. Whether a polymorphic equilibrium is reached (rather than fixation at both loci) depends on the intensity of selection and the relative distances to the optimum of the homozygotes at each locus. Furthermore, if both loci have similar effects on the trait, fixation of an allele at a given locus is less likely when it starts at low frequency and the other locus is polymorphic (with alleles at intermediate frequencies). Weaker selection increases the probability of fixation of the studied allele, as the polymorphic equilibrium is less stable in this case. When we do not require the double heterozygote to be at the optimum we find that the polymorphic equilibrium is more difficult to reach, and fixation becomes more likely. Second, increasing the number of loci decreases the probability of fixation, because adaptation to the optimum is possible by various combinations of alleles. Summaries of the genealogy (height, total length, and imbalance) and of sequence polymorphism (number of polymorphisms, frequency spectrum, and haplotype structure) next to a selected locus depend on the frequency that the selected mutation approaches at equilibrium. We conclude that multilocus response to selection may in some cases prevent selective sweeps from being completed, as described in previous studies, but that conditions causing this to happen strongly depend on the genetic architecture of the trait, and that fixation of selected mutations is likely in many instances.  相似文献   

3.
With a small effective population size, random genetic drift is more important than selection in determining the fate of new alleles. Small populations therefore accumulate deleterious mutations. Left unchecked, the effect of these fixed alleles is to reduce the reproductive capacity of a species, eventually to the point of extinction. New beneficial mutations, if fixed by selection, can restore some of this lost fitness. This paper derives the overall change in fitness due to fixation of new deleterious and beneficial alleles, as a function of the distribution of effects of new mutations and the effective population size. There is a critical effective size below which a population will on average decline in fitness, but above which beneficial mutations allow the population to persist. With reasonable estimates of the relevant parameters, this critical effective size is likely to be a few hundred. Furthermore, sexual selection can act to reduce the fixation probability of deleterious new mutations and increase the probability of fixing new beneficial mutations. Sexual selection can therefore reduce the risk of extinction of small populations.  相似文献   

4.
Genome-wide association studies (GWAS) have generated sufficient data to assess the role of selection in shaping allelic diversity of disease-associated SNPs. Negative selection against disease risk variants is expected to reduce their frequencies making them overrepresented in the group of minor (<50%) alleles. Indeed, we found that the overall proportion of risk alleles was higher among alleles with frequency <50% (minor alleles) compared to that in the group of major alleles. We hypothesized that negative selection may have different effects on environment (or lifestyle)-dependent versus environment (or lifestyle)-independent diseases. We used an environment/lifestyle index (ELI) to assess influence of environmental/lifestyle factors on disease etiology. ELI was defined as the number of publications mentioning “environment” or “lifestyle” AND disease per 1,000 disease-mentioning publications. We found that the frequency distributions of the risk alleles for the diseases with strong environmental/lifestyle components follow the distribution expected under a selectively neutral model, while frequency distributions of the risk alleles for the diseases with weak environmental/lifestyle influences is shifted to the lower values indicating effects of negative selection. We hypothesized that previously selectively neutral variants become risk alleles when environment changes. The hypothesis of ancestrally neutral, currently disadvantageous risk-associated alleles predicts that the distribution of risk alleles for the environment/lifestyle dependent diseases will follow a neutral model since natural selection has not had enough time to influence allele frequencies. The results of our analysis suggest that prediction of SNP functionality based on the level of evolutionary conservation may not be useful for SNPs associated with environment/lifestyle dependent diseases.  相似文献   

5.
Frequency-dependent selection (FDS) remains a common heuristic explanation for the maintenance of genetic variation in natural populations. The pairwise-interaction model (PIM) is a well-studied general model of frequency-dependent selection, which assumes that a genotype’s fitness is a function of within-population intergenotypic interactions. Previous theoretical work indicated that this type of model is able to sustain large numbers of alleles at a single locus when it incorporates recurrent mutation. These studies, however, have ignored the impact of the distribution of fitness effects of new mutations on the dynamics and end results of polymorphism construction. We suggest that a natural way to model mutation would be to assume mutant fitness is related to the fitness of the parental allele, i.e., the existing allele from which the mutant arose. Here we examine the numbers and distributions of fitnesses and alleles produced by construction under the PIM with mutation from parental alleles and the impacts on such measures due to different methods of generating mutant fitnesses. We find that, in comparison with previous results, generating mutants from existing alleles lowers the average number of alleles likely to be observed in a system subject to FDS, but produces polymorphisms that are highly stable and have realistic allele-frequency distributions.  相似文献   

6.
Kiyoshi Ezawa  Hideki Innan 《Genetics》2013,194(3):709-719
The population genetic behavior of mutations in sperm genes is theoretically investigated. We modeled the processes at two levels. One is the standard population genetic process, in which the population allele frequencies change generation by generation, depending on the difference in selective advantages. The other is the sperm competition during each genetic transmission from one generation to the next generation. For the sperm competition process, we formulate the situation where a huge number of sperm with alleles A and B, produced by a single heterozygous male, compete to fertilize a single egg. This “minimal model” demonstrates that a very slight difference in sperm performance amounts to quite a large difference between the alleles’ winning probabilities. By incorporating this effect of paternity-sharing sperm competition into the standard population genetic process, we show that fierce sperm competition can enhance the fixation probability of a mutation with a very small phenotypic effect at the single-sperm level, suggesting a contribution of sperm competition to rapid amino acid substitutions in haploid-expressed sperm genes. Considering recent genome-wide demonstrations that a substantial fraction of the mammalian sperm genes are haploid expressed, our model could provide a potential explanation of rapid evolution of sperm genes with a wide variety of functions (as long as they are expressed in the haploid phase). Another advantage of our model is that it is applicable to a wide range of species, irrespective of whether the species is externally fertilizing, polygamous, or monogamous. The theoretical result was applied to mammalian data to estimate the selection intensity on nonsynonymous mutations in sperm genes.  相似文献   

7.
We study the genetic basis of adaptation in a moving optimum model, in which the optimal value for a quantitative trait increases over time at a constant rate. We first analyze a one-locus two-allele model with recurrent mutation, for which we derive accurate analytical approximations for (i) the time at which a previously deleterious allele becomes beneficial, (ii) the waiting time for a successful new mutation, and (iii) the time the mutant allele needs to reach fixation. On the basis of these results, we show that the shortest total time to fixation is for alleles with intermediate phenotypic effect. We derive an approximation for this “optimal” effect, and we show that it depends in a simple way on a composite parameter, which integrates the ecological parameters and the genetic architecture of the trait. In a second step, we use stochastic computer simulations of a multilocus model to study the order in which mutant alleles with different effects go to fixation. In agreement with the one-locus results, alleles with intermediate effect tend to become fixed earlier than those with either small or large effects. However, the effect size of the fastest mutations differs from the one predicted in the one-locus model. We show how these differences can be explained by two specific effects of multilocus genetics. Finally, we discuss our results in the light of three relevant timescales acting in the system—the environmental, mutation, and fixation timescales—which define three parameter regimes leading to qualitative differences in the adaptive substitution pattern.  相似文献   

8.
The accumulation of beneficial mutations on competing genetic backgrounds in rapidly adapting populations has a striking impact on evolutionary dynamics. This effect, known as clonal interference, causes erratic fluctuations in the frequencies of observed mutations, randomizes the fixation times of successful mutations, and leaves distinct signatures on patterns of genetic variation. Here, we show how this form of “genetic draft” affects the forward-time dynamics of site frequencies in rapidly adapting asexual populations. We calculate the probability that mutations at individual sites shift in frequency over a characteristic timescale, extending Gillespie’s original model of draft to the case where many strongly selected beneficial mutations segregate simultaneously. We then derive the sojourn time of mutant alleles, the expected fixation time of successful mutants, and the site frequency spectrum of beneficial and neutral mutations. Finally, we show how this form of draft affects inferences in the McDonald–Kreitman test and how it relates to recent observations that some aspects of genetic diversity are described by the Bolthausen–Sznitman coalescent in the limit of very rapid adaptation.  相似文献   

9.
There are three basic paradigms of classical conditioning: delay, trace and context conditioning where presentation of a conditioned stimulus (CS) or a context typically predicts an unconditioned stimulus (US). In delay conditioning CS and US normally coterminate, whereas in trace conditioning an interval of time exists between CS termination and US onset. The modeling of trace conditioning is a rather difficult computational problem and is a challenge to the behavior and connectionist approaches mainly due to a time gap between CS and US. To account for trace conditioning, Pavlov (Conditioned reflexes: an investigation of the physiological activity of the cerebral cortex, Oxford University Press, London, 1927) postulated the existence of a stimulus “trace” in the nervous system. Meanwhile, there exist many other options for solving this association problem. There are several excellent reviews of computational models of classical conditioning but none has thus far been devoted to trace conditioning. Eight representative models of trace conditioning aimed at building a prospective model are being reviewed below in a brief form. As a result, one of them, comprising the most important features of its predecessors, can be suggested as a real candidate for a unified model of trace conditioning.  相似文献   

10.
The Effects of Overdominance on Linkage in a Multilocus System   总被引:3,自引:2,他引:1       下载免费PDF全文
Computer simulations were performed with overdominant multiple alleles among tightly linked multiple loci under a multiplicative fitness model. The quantity X2/N(n — 1) was introduced as a new measure of linkage disequilibrium which, unlike previously available measures, can be applied to multiple allele models, where N is the sample size, and n is the number of alleles at the locus possessing fewest alleles. Simulations showed that (1) With multiple (three or four) alleles, the approach to stable disequilibrium is slower and the amount of disequilibrium established is weaker than in a two allele system. (2) The number of complementary chromosomes is a function of number of alleles and of population size. (3) As population size increases, the rate of the approach to stable disequilibrium is slower. (4) There is an optimum selection coefficient which minimizes the transient fixation probability of alleles when linkage is present. (5) The absence of linkage disequilibrium is in most cases not a practical method of testing the hypothesis of balancing selection of genetic polymorphisms because it depends strongly on population size in determining linkage disequilibria.  相似文献   

11.
We study the evolution of a pair of competing behavioural alleles in a structured population when there are non-additive or ‘synergistic’ fitness effects. Under a form of weak selection and with a simple symmetry condition between a pair of competing alleles, Tarnita et al. provide a surprisingly simple condition for one allele to dominate the other. Their condition can be obtained from an analysis of a corresponding simpler model in which fitness effects are additive. Their result uses an average measure of selective advantage where the average is taken over the long-term—that is, over all possible allele frequencies—and this precludes consideration of any frequency dependence the allelic fitness might exhibit. However, in a considerable body of work with non-additive fitness effects—for example, hawk–dove and prisoner''s dilemma games—frequency dependence plays an essential role in the establishment of conditions for a stable allele-frequency equilibrium. Here, we present a frequency-dependent generalization of their result that provides an expression for allelic fitness at any given allele frequency p. We use an inclusive fitness approach and provide two examples for an infinite structured population. We illustrate our results with an analysis of the hawk–dove game.  相似文献   

12.
Pseudogenization is a widespread phenomenon in genome evolution, and it has been proposed to serve as an engine of evolutionary change, especially during human origins (the “less-is-more” hypothesis). However, there has been no comprehensive analysis of human-specific pseudogenes. Furthermore, it is unclear whether pseudogenization itself can be selectively favored and thus play an active role in human evolution. Here we conduct a comparative genomic analysis and a literature survey to identify 80 nonprocessed pseudogenes that were inactivated in the human lineage after its separation from the chimpanzee lineage. Many functions are involved among these genes, with chemoreception and immune response being outstandingly overrepresented, suggesting potential species-specific features in these aspects of human physiology. To explore the possibility of adaptive pseudogenization, we focus on CASPASE12, a cysteinyl aspartate proteinase participating in inflammatory and innate immune response to endotoxins. We provide population genetic evidence that the nearly complete fixation of a null allele at CASPASE12 has been driven by positive selection, probably because the null allele confers protection from severe sepsis. We estimate that the selective advantage of the null allele is about 0.9% and the pseudogenization started shortly before the out-of-Africa migration of modern humans. Interestingly, two other genes related to sepsis were also pseudogenized in humans, possibly by selection. These adaptive gene losses might have occurred because of changes in our environment or genetic background that altered the threat from or response to sepsis. The identification and analysis of human-specific pseudogenes open the door for understanding the roles of gene losses in human origins, and the demonstration that gene loss itself can be adaptive supports and extends the “less-is-more” hypothesis.  相似文献   

13.
Central to Wright's shifting-balance theory is the idea that genetic drift and selection in systems with gene interaction can lead to the formation of “adaptive gene complexes.” The theory of genetic drift has been well developed over the last 60 years; however, nearly all of this theory is based on the assumption that only additive gene effects are acting. Wright's theory was developed recognizing that there was a “universality of interaction effects,” which implies that additive theory may not be adequate to describe the process of differentiation that Wright was considering. The concept of an adaptive gene complex implies that an allele that is favored by individual selection in one deme may be removed by selection in another deme. In quantitative genetic terms, the average effects of an allele relative to other alleles changes from deme to deme. The model presented here examines the variance in local breeding values (LBVs) of a single individual and the covariance in the LBVs of a pair of individuals mated in the same deme relative to when they are mated in different demes. Local breeding value is a measure of the average effects of the alleles that make up that individual in a particular deme. I show that when there are only additive effects the covariance between the LBVs of individuals equals the variance in the LBV of an individual. As the amount of epistasis in the ancestral population increases, the variance in the LBV of an individual increases and the covariance between the LBVs of a pair of individuals decreases. The divergence in these two values is a measure of the extent to which the LBV of an individual varies independently of the LBVs of other individuals. When this value is large, it means that the relative ordering of the average effects of alleles will change from deme to deme. These results confirm an important component of Wright's shifting-balance theory: When there is gene interaction, genetic drift can lead to the reordering of the average effects of alleles and when coupled with selection this will lead to the formation of the adaptive gene complexes.  相似文献   

14.
A recent commentary raised concerns about aspects of the model and assumptions used in a previous study which demonstrated that selection can favor chromosomal alleles that confer higher plasmid donation rates. Here, the authors of that previous study respond to the concerns raised.

In our original work [1], we demonstrated experimentally that selection can favor chromosomal alleles that confer higher plasmid donation rates, given the plasmid is beneficial and the recipient has an elevated chance of carrying the donor allele (i.e., preferential donation to kin). Our experiments demonstrated this effect via 2 mechanisms of preferential donation: biased conjugation rates and structured populations. We interpreted these results through the lens of kin selection theory (benefits via horizontal gene transfer to kin), supported by simulations and an analytical fitness function model. These results hold importance by outlining that the evolution of plasmid transfer rates (a key aspect of the antibiotic resistance crisis) is not necessarily the sole product of selection on the plasmid itself and forms part of a broader series of papers from our labs investigating the sociomicrobiology of plasmids [24].A new commentary raises concerns over our fitness function model, flagging issues with both the structure of the model and assumptions made in our analysis [5]. We stand by the general conclusions of our work but accept that our fitness function and stated analysis assumptions could be better formulated. Our initial fitness function is heuristic in the sense it was designed to capture general processes acting on the fitness of individuals, dependent on the plasmid and donor allele status—without explicitly modeling the myriad demographic events of dispersal, reproduction, conjugation, and death that result in selective shifts across a metapopulation of cells. Specifically, we captured the “force of infection” faced by an uninfected cell as the product of average plasmid prevalence and average donor allele prevalence in the local patch (pjqj; see commentary for notation details). We agree with the authors that this force of infection is better phrased as the average of the product ((1/N)∑pij qij), in part because this avoids the potential pathology under limit conditions described by the authors, but also because this approach better highlights that the particular social trait in question is an “other only” cooperative trait [6], illustrated by commentary equation [2], where transmission to self and transmission to others are separated. This separation has the important consequence of highlighting that unlike many microbial social traits where benefits accrue to a group (including self), a cooperative plasmid donor trait can only benefit other cells that lack the plasmid. Given established costs of donation (e.g., see figure S2 in our original article), this defines our “donor” behavior as an altruistic trait, which can, therefore, only be favored by selection given nonrandom interactions among individuals (e.g., [7]).Our experimental results outline 2 mechanisms of nonrandom interactions: preferential donation to kin and population structure. Each of these mechanisms will generate positive covariances between focal individual qij and non-self-recipient qj donor allele states (cov(qj, qij) > 0). The pathway via preferential donation to kin (order-of-magnitude differences according to our analyses and more recent measurements among lineages coexisting within natural populations [8]) will also likely generate positive covariances between donor and recipient abilities (cov(sij, qij) > 0). In contrast, to arrive at the result that selection always works against plasmid donor alleles (equation [4]), the commentary makes the assumption that both of the above covariances are zero. We suggest that the additional analyses begun by the authors are an exciting starting point to better map selection on donor alleles, under a broader array of defined assumptions on cell–cell and gene–gene structure, ideally informed by data on structures found in natural bacterial populations.  相似文献   

15.
Self-fertilization is generally seen to be disadvantageous in the long term. It increases genetic drift, which subsequently reduces polymorphism and the efficiency of selection, which also challenges adaptation. However, high selfing rates can increase the fixation probability of recessive beneficial mutations, but existing theory has generally not accounted for the effect of linked sites. Here, we analyze a model for the fixation probability of deleterious mutants that hitchhike with selective sweeps in diploid, partially selfing populations. Approximate analytical solutions show that, conditional on the sweep not being lost by drift, higher inbreeding rates increase the fixation probability of the deleterious allele, due to the resulting reduction in polymorphism and effective recombination. When extending the analysis to consider a distribution of deleterious alleles, as well as the average fitness increase after a sweep, we find that beneficial alleles generally need to be more recessive than the previously assumed dominance threshold (h < 1/2) for selfing to be beneficial from one-locus theory. Our results highlight that recombination aiding the efficiency of selection on multiple loci amplifies the fitness benefits of outcrossing over selfing, compared to results obtained from one-locus theory. This effect additionally increases the parameter range under which obligate outcrossing is beneficial over partial selfing.  相似文献   

16.
Mating type in the Gibberella fujikuroi species complex is controlled by a single locus with two alleles and is usually identified following sexual crosses with standard, female-fertile tester isolates. The mating type alleles have been arbitrarily designated “+” and “−” within each biological species, and the nomenclature is tied to the standard tester strains. We developed a pair of PCR primers that can be used to amplify a unique fragment of one of the mating type alleles (MAT-2) from at least seven of the biological species in this species complex. Based on the amplification pattern, we propose a replacement for the existing, arbitrary +/− terminology that is presently in use. The new terminology is based on DNA sequence similarities between the mating type allele fragments from the biological species of the G. fujikuroi species complex and the corresponding fragments from other filamentous ascomycetes.  相似文献   

17.
Female mammals are functional mosaics of their parental X-linked gene expression due to X chromosome inactivation (XCI). This process inactivates one copy of the X chromosome in each cell during embryogenesis and that state is maintained clonally through mitosis. In mice, the choice of which parental X chromosome remains active is determined by the X chromosome controlling element (Xce), which has been mapped to a 176-kb candidate interval. A series of functional Xce alleles has been characterized or inferred for classical inbred strains based on biased, or skewed, inactivation of the parental X chromosomes in crosses between strains. To further explore the function structure basis and location of the Xce, we measured allele-specific expression of X-linked genes in a large population of F1 females generated from Collaborative Cross (CC) strains. Using published sequence data and applying a Bayesian “Pólya urn” model of XCI skew, we report two major findings. First, inter-individual variability in XCI suggests mouse epiblasts contain on average 20–30 cells contributing to brain. Second, CC founder strain NOD/ShiLtJ has a novel and unique functional allele, Xceg, that is the weakest in the Xce allelic series. Despite phylogenetic analysis confirming that NOD/ShiLtJ carries a haplotype almost identical to the well-characterized C57BL/6J (Xceb), we observed unexpected patterns of XCI skewing in females carrying the NOD/ShiLtJ haplotype within the Xce. Copy number variation is common at the Xce locus and we conclude that the observed allelic series is a product of independent and recurring duplications shared between weak Xce alleles.  相似文献   

18.
In many species, genomic data have revealed pervasive adaptive evolution indicated by the fixation of beneficial alleles. However, when selection pressures are highly variable along a species'' range or through time adaptive alleles may persist at intermediate frequencies for long periods. So called “balanced polymorphisms” have long been understood to be an important component of standing genetic variation, yet direct evidence of the strength of balancing selection and the stability and prevalence of balanced polymorphisms has remained elusive. We hypothesized that environmental fluctuations among seasons in a North American orchard would impose temporally variable selection on Drosophila melanogaster that would drive repeatable adaptive oscillations at balanced polymorphisms. We identified hundreds of polymorphisms whose frequency oscillates among seasons and argue that these loci are subject to strong, temporally variable selection. We show that these polymorphisms respond to acute and persistent changes in climate and are associated in predictable ways with seasonally variable phenotypes. In addition, our results suggest that adaptively oscillating polymorphisms are likely millions of years old, with some possibly predating the divergence between D. melanogaster and D. simulans. Taken together, our results are consistent with a model of balancing selection wherein rapid temporal fluctuations in climate over generational time promotes adaptive genetic diversity at loci underlying polygenic variation in fitness related phenotypes.  相似文献   

19.
Three days of fear conditioning that combines tactile stimulation of a row of facial vibrissae (conditioned stimulus, CS) with a tail shock (unconditioned stimulus, UCS) expands the representation of “trained” vibrissae, which can be demonstrated by labeling with 2-deoxyglucose in layer IV of the barrel cortex. We have also shown that functional reorganization of the primary somatosensory cortex (S1) increases GABAergic markers in the hollows of “trained” barrels of the adult mouse. This study investigated how whisker-shock conditioning (CS+UCS) affected the expression of puncta of a high-affinity GABA plasma membrane transporter GAT-1 in the barrel cortex of mice 24 h after associative learning paradigm. We found that whisker-shock conditioning (CS+UCS) led to increase expression of neuronal and astroglial GAT-1 puncta in the “trained” row compared to controls: Pseudoconditioned, CS-only, UCS-only and Naïve animals. These findings suggest that fear conditioning specifically induces activation of systems regulating cellular levels of the inhibitory neurotransmitter GABA.  相似文献   

20.
Defects in mismatch repair (MMR) genes result in a mutator phenotype by inducing microsatellite instability (MI), a characteristic of hereditary nonpolyposis colorectal cancers (HNPCC) and a subset of sporadic colon tumors. Present models describing the mechanism by which germ line mutations in MMR genes predispose kindreds to HNPCC suggest a “two-hit” inactivation of both alleles of a particular MMR gene. Here we present experimental evidence that a nonsense mutation at codon 134 of the hPMS2 gene is sufficient to reduce MMR and induce MI in cells containing a wild-type hPMS2 allele. These results have significant implications for understanding the relationship between mutagenesis and carcinogenesis and the ability to generate mammalian cells with mutator phenotypes.  相似文献   

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

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