首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Whereas detection and positioning of genes that affect quantitative traits (quantitative trait loci (QTL)) using linkage mapping uses only information from recombinants in the genotyped generations, linkage disequilibrium (LD) mapping uses historical recombinants. Thus, whereas linkage mapping requires large family sizes to detect and accurately position QTL, LD mapping is more dependent on the number of families sampled from the population. In commercial Atlantic salmon breeding programmes, only a small number of individuals per family are routinely phenotyped for traits such as disease resistance and meat colour. In this paper, we assess the power and accuracy of combined linkage disequilibrium linkage analysis (LDLA) to detect QTL in the commercial population using simulation. When 15 half-sib sire families (each sire mated to 30 dams, each dam with 10 progeny) were sampled from the population for genotyping, we were able to detect a QTL explaining 10% of the phenotypic variance in 85% of replicates and position this QTL within 3 cM of the true position in 70% of replicates. When recombination was absent in males, a feature of the salmon genome, power to detect QTL increased; however, the accuracy of positioning the QTL was decreased. By increasing the number of sire families sampled from the population to be genotyped to 30, we were able to increase both the proportion of QTL detected and correctly positioned (even with no recombination in males). QTL with much smaller effect could also be detected. The results suggest that even with the existing recording structure in commercial salmon breeding programmes, there is considerable power to detect and accurately position QTL using LDLA.  相似文献   

2.
Li H  Bradbury P  Ersoz E  Buckler ES  Wang J 《PloS one》2011,6(3):e17573

Background

Nested association mapping (NAM) is a novel genetic mating design that combines the advantages of linkage analysis and association mapping. This design provides opportunities to study the inheritance of complex traits, but also requires more advanced statistical methods. In this paper, we present the detailed algorithm of a QTL linkage mapping method suitable for genetic populations derived from NAM designs. This method is called joint inclusive composite interval mapping (JICIM). Simulations were designed on the detected QTL in a maize NAM population and an Arabidopsis NAM population so as to evaluate the efficiency of the NAM design and the JICIM method.

Principal Findings

Fifty-two QTL were identified in the maize population, explaining 89% of the phenotypic variance of days to silking, and nine QTL were identified in the Arabidopsis population, explaining 83% of the phenotypic variance of flowering time. Simulations indicated that the detection power of these identified QTL was consistently high, especially for large-effect QTL. For rare QTL having significant effects in only one family, the power of correct detection within the 5 cM support interval was around 80% for 1-day effect QTL in the maize population, and for 3-day effect QTL in the Arabidopsis population. For smaller-effect QTL, the power diminished, e.g., it was around 50% for maize QTL with an effect of 0.5 day. When QTL were linked at a distance of 5 cM, the likelihood of mapping them as two distinct QTL was about 70% in the maize population. When the linkage distance was 1 cM, they were more likely mapped as one single QTL at an intermediary position.

Conclusions

Because it takes advantage of the large genetic variation among parental lines and the large population size, NAM is a powerful multiple-cross design for complex trait dissection. JICIM is an efficient and specialty method for the joint QTL linkage mapping of genetic populations derived from the NAM design.  相似文献   

3.
Stella A  Boettcher PJ 《Genetics》2004,166(1):341-350
Simulation was used to evaluate the performance of different selective genotyping strategies when using linkage disequilibrium across large half-sib families to position a QTL within a previously defined genomic region. Strategies examined included standard selective genotyping and different approaches of discordant and concordant sib selection applied to arbitrary or selected families. Strategies were compared as a function of effect and frequency of QTL alleles, heritability, and phenotypic expression of the trait. Large half-sib families were simulated for 100 generations and 2% of the population was genotyped in the final generation. Simple ANOVA was applied and the marker with the greatest F-value was considered the most likely QTL position. For traits with continuous phenotypes, genotyping the most divergent pairs of half-sibs from all families was the best strategy in general, but standard selective genotyping was somewhat more precise when heritability was low. When the phenotype was distributed in ordered categories, discordant sib selection was the optimal approach for positioning QTL for traits with high heritability and concordant sib selection was the best approach when genetic effects were small. Genotyping of a few selected sibs from many families was generally more efficient than genotyping many individuals from a few highly selected sires.  相似文献   

4.
Meuwissen TH  Goddard ME 《Genetics》2000,155(1):421-430
A multimarker linkage disequilibrium mapping method was developed for the fine mapping of quantitative trait loci (QTL) using a dense marker map. The method compares the expected covariances between haplotype effects given a postulated QTL position to the covariances that are found in the data. The expected covariances between the haplotype effects are proportional to the probability that the QTL position is identical by descent (IBD) given the marker haplotype information, which is calculated using the genedropping method. Simulation results showed that a QTL was correctly positioned within a region of 3, 1.5, or 0.75 cM in 70, 62, and 68%, respectively, of the replicates using markers spaced at intervals of 1, 0.5, and 0.25 cM, respectively. These results were rather insensitive to the number of generations since the QTL occurred and to the effective population size, except that 10 generations yielded rather poor estimates of the QTL position. The position estimates of this multimarker disequilibrium mapping method were more accurate than those from a single marker transmission disequilibrium test. A general approach for identifying QTL is suggested, where several stages of disequilibrium mapping are used with increasingly dense marker spacing.  相似文献   

5.
Abstract

The development of dense linkage maps consisting of highly polymorphic loci for livestock species is technically feasible. However, linkage mapping experiments are expensive as they involve many animals and marker typings per animal. To minimize costs of developing linkage maps for livestock species, optimizing designs for mapping studies is necessary. This study provides a general framework for evaluating the efficiency of designs for reference families consisting of two‐ or three‐ generation full‐sib or half‐sib families selected from a segregating population. The influence of number of families, number of offspring per family, family structure (either half‐sib or full‐sib) and marker polymorphism is determined. Evaluation is done for two markers with a recombination rate of .20 and for a marker and a dominant single gene with a recombination rate of .20. Two evaluation criteria are used: expected maximum lod score for detection of linkage and accuracy of an estimated recombination rate defined as probability that the true recombination rate is in an interval around the estimated recombination rate. First, for several designs the contribution of reference families to expected maximum lod score and accuracy is given. Second, the required number of families in a design to obtain a certain value for the evaluation criteria is calculated when number of offspring per family, family structure and marker polymorphism are specified. The required numbers increase when designs are optimized not only for expected maximum lod score but also for accuracy. The required number of animals to map a dominant single gene is very large. Therefore, a set of reference families should be designed for strictly mapping marker loci. Examples illustrate how tabulated results can be generalized to determine the values for a wide range of designs containing two‐ or three‐generation full‐sib or half‐sib families.  相似文献   

6.
A previous genome scan that was conducted in Spanish Churra sheep identified a significant quantitative trait locus (QTL) for milk protein percentage (PP) on chromosome 3 (OAR3), between markers KD103 and OARVH34. The aim of this study was to replicate these results and to refine the mapped position of this QTL. To accomplish this goal, we analysed 14 new half‐sib families of Spanish Churra sheep including 1661 ewes from 29 different flocks. These animals were genotyped for 21 microsatellite markers mapping to OAR3. In addition to a classical linkage analysis (LA), a combined linkage disequilibrium and linkage analysis (LDLA) was performed with the aim of enhancing the resolution of the QTL mapping. The LA that was performed in this sheep population identified the presence of a highly significant QTL for PP near marker KD103 (Pc < 0.001; Pexp < 0.001). The phenotypic variance that was owing to the QTL was 2.74%. Two segregating families for the target QTL were identified in this population with QTL effect estimates of 0.47 and 0.95 SD. The LDLA identified the same QTL as the previous analyses with a high level of statistical significance (P = 9.184 E‐11) and narrowed the confidence interval (CI) to a 13 cM region. These results confirm the segregation of the previously identified OAR3 QTL that influences PP in Spanish Churra sheep. Future research will aim to increase the marker density across the refined CI and to analyse the corresponding candidate genes to identify the allelic variant or variants that underlie this genetic effect.  相似文献   

7.
Natural mating and mass spawning in the European sea bass (Dicentrarchus labrax L., Moronidae, Teleostei) complicate genetic studies and the implementation of selective breeding schemes. We utilized a two‐step experimental design for detecting QTL in mass‐spawning species: 2122 offspring from natural mating between 57 parents (22 males, 34 females and one missing) phenotyped for body weight, eight morphometric traits and cortisol levels, had been previously assigned to parents based on genotypes of 31 DNA microsatellite markers. Five large full‐sib families (five sires and two dams) were selected from the offspring (570 animals), which were genotyped with 67 additional markers. A new genetic map was compiled, specific to our population, but based on the previously published map. QTL mapping was performed with two methods: half‐sib regression analysis (paternal and maternal) and variance component analysis accounting for all family relationships. Two significant QTL were found for body weight on linkage group 4 and 6, six significant QTL for morphometric traits on linkage groups 1B, 4, 6, 7, 15 and 23 and three suggestive QTL for stress response on linkage groups 3, 14 and 23. The QTL explained between 8% and 38% of phenotypic variance. The results are the first step towards identifying genes involved in economically important traits like body weight and stress response in European sea bass.  相似文献   

8.
A novel and robust method for the fine-scale mapping of genes affecting complex traits, which combines linkage and linkage-disequilibrium information, is proposed. Linkage information refers to recombinations within the marker-genotyped generations and linkage disequilibrium to historical recombinations before genotyping started. The identity-by-descent (IBD) probabilities at the quantitative trait locus (QTL) between first generation haplotypes were obtained from the similarity of the marker alleles surrounding the QTL, whereas IBD probabilities at the QTL between later generation haplotypes were obtained by using the markers to trace the inheritance of the QTL. The variance explained by the QTL is estimated by residual maximum likelihood using the correlation structure defined by the IBD probabilities. Unlinked background genes were accounted for by fitting a polygenic variance component. The method was used to fine map a QTL for twinning rate in cattle, previously mapped on chromosome 5 by linkage analysis. The data consisted of large half-sib families, but the method could also handle more complex pedigrees. The likelihood of the putative QTL was very small along most of the chromosome, except for a sharp likelihood peak in the ninth marker bracket, which positioned the QTL within a region <1 cM in the middle part of bovine chromosome 5. The method was expected to be robust against multiple genes affecting the trait, multiple mutations at the QTL, and relatively low marker density.  相似文献   

9.
The objective of this study was to locate quantitative trait loci (QTL) causing variation in birth weight and age of puberty of doe kids in a population of Rayini cashmere goats. Four hundred and thirty kids from five half‐sib families were genotyped for 116 microsatellite markers located on the caprine autosomes. The traits recorded were birth weight of the male and female kids, body weight at puberty, average daily gain from birth to age of puberty and age at puberty of the doe kids. QTL analysis was conducted using the least squares interval mapping approach. Linkage analysis indicated significant QTL for birth weight on Capra hircus chromosomes (CHI) 4, 5, 6, 18 and 21. Five QTL located on CHI 5, 14 and 29 were associated with age at puberty. Across‐family analysis revealed evidence for overlapping QTL affecting birth weight (78 cM), body weight at puberty (72 cM), average daily gain from birth to age of puberty (72 cM) and age at puberty (76 cM) on CHI 5 and overlapping QTL controlling body weight at puberty and age at puberty on CHI 14 at 18–19 cM. The proportion of the phenotypic variance explained by the detected QTL ranged between 7.9% and 14.4%. Confirming some of the previously reported results for birth weight and growth QTL in goats, this study identified more QTL for these traits and is the first report of QTL for onset of puberty in doe kids.  相似文献   

10.
RAPD markers were employed for construction of a linkage map and localization of QTLs for oleic acid level using a set of 94 recombinant inbred lines (RILs) of mustard (Brassica juncea L.) as a mapping population. Only 30% of the 235 random primers used were useful in terms of polymorphism detected and the reproducibility of those patterns. Normal Mendelian segregation was observed for the majority of the 130 markers obtained with 71 informative primers; only 13.1% deviated (P < 0.01) from the expected 1:1 ratio. One-hundred and fourteen markers were assigned to 21 linkage groups (LGs) covering a total length of 790.4 cM with an average distance of 6.93 cM between markers. Two quantitative trait loci (QTL) for oleic acid level were mapped to 14- and 10.6-cM marker intervals on two different LGs. Both loci together explained 32.2% of phenotypic variance. One major QTL explained 28.5% of the trait variance observed in this species.  相似文献   

11.
We analysed a QTL affecting milk yield (MY), milk protein yield (PY) and milk fat yield (FY) in the dual purpose cattle breed Fleckvieh on BTA5. Twenty-six microsatellite markers covering 135 cM were selected to analyse nine half-sib families containing 605 sons in a granddaughter design. We thereby assigned two new markers to the public linkage map using the CRI-MAP program. Phenotypic records were daughter yield deviations (DYD) originating from the routinely performed genetic evaluations of breeding animals. To determine the position of the QTL, three different approaches were applied: interval mapping (IM), linkage analysis by variance component analysis (LAVC), and combined linkage disequilibrium (LD) and linkage (LDL) analysis. All three methods mapped the QTL in the same marker interval ( BM2830-ETH152 ) with the greatest test-statistic value at 118, 119.33 and 119.33 cM respectively. The positive QTL allele simultaneously increases DYD in the first lactation by 272 kg milk, 7.1 kg milk protein and 7.0 kg milk fat. Although the mapping accuracy and the significance of a QTL effect increased from IM over LAVC to LDL, the confidence interval was large (13, 20 and 24 cM for FY, MY and PY respectively) for the positional cloning of the causal gene. The estimated averages of pair wise marker LD with a distance <5 cM were low (0.107) and reflect the large effective population size of the Fleckvieh subpopulation analysed. This low level of LD suggests a need for increase in marker density in following fine mapping steps.  相似文献   

12.
四倍体栽培棉种产量和纤维品质性状的QTL定位   总被引:29,自引:1,他引:28  
陆地棉和海岛棉是两个不同的四倍体栽培种 ,但在生产上各有其特点 ,陆地棉丰产性强 ,海岛棉纤维品质优良 ,利用其种间杂交群体定位产量和品质性状的QTL ,对于分子标记辅助的海岛棉优质纤维向陆地棉转移很有意义。以SSR和RAPD为分子标记 ,陆地棉与海岛棉杂种 (邯郸 2 0 8×Pima90 )F2 群体为作图群体 ,构建了一张含 12 6个标记的遗传图谱 ,包括 6 8个SSR标记和 5 8个RAPD标记 ,可分为 2 9个连锁群 ,标记间平均距离为 13 7cM ,总长1717 0cM ,覆盖棉花总基因组约 34 34% ;以遗传图 12 6个标记为基础 ,对F2 :3 家系符合正态分布的 10个农艺性状及纤维品质性状进行全基因组QTL扫描 ,结果发现 2 9个QTL分别与产量和品质性状有关。其中与衣指、籽指、皮棉产量、子棉产量、衣分等产量性状相关的QTL分别有 1、3、5、6和 1个 ,与纤维长度、整齐度、强度、伸长率和马克隆值等品质性状相关的QTL分别有 2、4、2、4和 1个。各QTL解释的变异量在 12 4 2 %~ 47 0 1%之间。其中比强度有关的 2个QTL能够解释的表型变异率分别为 34 15 %和 13 86 %。  相似文献   

13.
In a simulation study different designs for a pure line pig population were compared for efficiency of mapping QTL using the variance component method. Phenotypes affected by a Mendelian QTL, a paternally expressed QTL, a maternally expressed QTL or by a QTL without an effect were simulated. In all alternative designs 960 progeny were phenotyped. Given the limited number of animals there is an optimum between the number of families and the family size. Estimation of Mendelian and parentally expressed QTL is more efficient in a design with large family sizes. Too small a number of sires should be avoided to minimize chances of sires to be non-segregating. When a large number of families is used, the number of haplotypes increases which reduces the accuracy of estimating the QTL effect and thereby reduces the power to show a significant QTL and to correctly position the QTL. Dense maps allow for smaller family size due to exploitation of LD-information. Given the different possible modes of inheritance of the QTL using 8 to16 boars, two litters per dam was optimal with respect to determining significance and correct location of the QTL for a data set consisting of 960 progeny. The variance component method combining linkage disequilibrium and linkage analysis seems to be an appropriate choice to analyze data sets which vary in marker density and which contain complex family structures.  相似文献   

14.
Cadmium (Cd) is a widespread soil pollutant and poses a significant threat to human health via the food chain. Large phenotypic variations in Cd concentration of radish roots and shoots have been observed. However, the genetic and molecular mechanisms of Cd accumulation in radish remain to be elucidated. In this study, a genetic linkage map was constructed using an F(2) mapping population derived from a cross between a high Cd-accumulating cultivar NAU-Dysx and a low Cd-accumulating cultivar NAU-Yh. The linkage map consisted of 523 SRAP, RAPD, SSR, ISSR, RAMP, and RGA markers and had a total length of 1,678.2 cM with a mean distance of 3.4 cM between two markers. All mapped markers distributed on nine linkage groups (LGs) having sizes between 134.7 and 236.8 cM. Four quantitative trait loci (QTLs) for root Cd accumulation were mapped on LGs 1, 4, 6, and 9, which accounted for 9.86 to 48.64 % of all phenotypic variance. Two QTLs associated with shoot Cd accumulation were detected on LG1 and 3, which accounted for 17.08 and 29.53 % of phenotypic variance, respectively. A major-effect QTL, qRCd9 (QTL for root Cd accumulation on LG9), was identified on LG 9 flanked by NAUrp011_754 and EM5me6_286 markers with a high LOD value of 23.6, which accounted for 48.64 % of the total phenotypic variance in Cd accumulation of F(2) lines. The results indicated that qRCd9 is a novel QTL responsible for controlling root Cd accumulation in radish, and the identification of specific molecular markers tightly linked to the major QTL could be further applied for marker-assisted selection (MAS) in low-Cd content radish breeding program.  相似文献   

15.
路明  周芳  谢传晓  李明顺  徐云碧  张世煌 《遗传》2007,29(9):1131-1138
为了增加单位面积产量, 玉米育种者已经开始了更密植更紧凑株型的选育。叶夹角和叶向值是评价玉米株型的重要指标。本研究以掖478×丹340的500个F2单株为作图群体, 构建了具有138个位点的SSR标记连锁图谱, 图谱总长度为1 394.9 cM, 平均间距10.1 cM。利用397个F2:3家系对叶夹角和叶向值进行QTL定位分析, 结果表明: 叶夹角和叶向值分别检测到6和8个QTL, 累计解释表型变异41.0%和60.8%, 单个QTL的贡献率在2.9%~13.6%之间。与叶夹角和叶向值有关的基因主要作用方式为加性和部分显性。此外两个性状共检测到9对上位性互作位点, 表明上位性互作在叶夹角和叶向值的遗传中也起较重要的作用。  相似文献   

16.
Results from power studies for linkage detection have led to many ongoing and planned collections of phenotypically extreme nuclear families. Given the great expense of collecting these families and the imminent availability of a dense diallelic marker map, the families are likely to be used in allelic-association as well as linkage studies. However, optimal selection strategies for linkage may not be equally powerful for association. We examine the power to detect linkage disequilibrium for quantitative traits after phenotypic selection. The results encompass six selection strategies that are in widespread use, including single selection (two designs), affected sib pairs, concordant and discordant pairs, and the extreme-concordant and -discordant design. Selection of sibships on the basis of one extreme proband with high or low trait scores provides as much power as discordant sib pairs but requires the screening and phenotyping of substantially fewer initial families from which to select. Analysis of the role of allele frequencies within each selection design indicates that common trait alleles generally offer the most power, but similarities between the marker- and trait-allele frequencies are much more important than the trait-locus frequency alone. Some of the most widespread selection designs, such as single selection, yield power gains only when both the marker and quantitative trait loci (QTL) are relatively rare in the population. In contrast, discordant pairs and the extreme-proband design provide power for the broadest range of QTL-marker-allele frequency differences. Overall, proband selection from either tail provides the best balance of power, robustness, and simplicity of ascertainment for family-based association analysis.  相似文献   

17.

Background

Information for mapping of quantitative trait loci (QTL) comes from two sources: linkage disequilibrium (non-random association of allele states) and cosegregation (non-random association of allele origin). Information from LD can be captured by modeling conditional means and variances at the QTL given marker information. Similarly, information from cosegregation can be captured by modeling conditional covariances. Here, we consider a Bayesian model based on gene frequency (BGF) where both conditional means and variances are modeled as a function of the conditional gene frequencies at the QTL. The parameters in this model include these gene frequencies, additive effect of the QTL, its location, and the residual variance. Bayesian methodology was used to estimate these parameters. The priors used were: logit-normal for gene frequencies, normal for the additive effect, uniform for location, and inverse chi-square for the residual variance. Computer simulation was used to compare the power to detect and accuracy to map QTL by this method with those from least squares analysis using a regression model (LSR).

Results

To simplify the analysis, data from unrelated individuals in a purebred population were simulated, where only LD information contributes to map the QTL. LD was simulated in a chromosomal segment of 1 cM with one QTL by random mating in a population of size 500 for 1000 generations and in a population of size 100 for 50 generations. The comparison was studied under a range of conditions, which included SNP density of 0.1, 0.05 or 0.02 cM, sample size of 500 or 1000, and phenotypic variance explained by QTL of 2 or 5%. Both 1 and 2-SNP models were considered. Power to detect the QTL for the BGF, ranged from 0.4 to 0.99, and close or equal to the power of the regression using least squares (LSR). Precision to map QTL position of BGF, quantified by the mean absolute error, ranged from 0.11 to 0.21 cM for BGF, and was better than the precision of LSR, which ranged from 0.12 to 0.25 cM.

Conclusions

In conclusion given a high SNP density, the gene frequency model can be used to map QTL with considerable accuracy even within a 1 cM region.  相似文献   

18.
Recent studies report a surprisingly high degree of marker-to-marker linkage disequilibrium (LD) in ruminant livestock populations. This has important implications for QTL mapping and marker-assisted selection. This study evaluated LD between microsatellite markers in a number of breeding populations of layer chickens using the standardized chi-square (chi(2')) measure. The results show appreciable LD among markers separated by up to 5 cM, decreasing rapidly with increased separation between markers. The LD within 5 cM was strongly conserved across generations and differed among chromosomal regions. Using marker-to-marker LD as an indication for marker-QTL LD, a genome scan of markers spaced 2 cM apart at moderate power would have good chances of uncovering most QTL segregating in these populations. However, of markers showing significant trait associations, only 57% are expected to be within 5 cM of the responsible QTL, and the remainder will be up to 20 cM away. Thus, high-resolution LD mapping of QTL will require dense marker genotyping across the region of interest to allow for interval mapping of the QTL.  相似文献   

19.
The power to separate the variance of a quantitative trait locus (QTL) from the polygenic variance is determined by the variability of genes identical by descent (IBD) at the QTL. This variability may increase with inbreeding. Selfing, the most extreme form of inbreeding, increases the variability of the IBD value shared by siblings, and thus has a higher efficiency for QTL mapping than random mating. In self-incompatible organisms, sib mating is the closest form of inbreeding. Similar to selfing, sib mating may also increase the power of QTL detection relative to random mating. In this study, we develop an IBD-based method under sib mating designs for QTL mapping. The efficiency of sib mating is then compared with random mating. Monte Carlo simulations show that sib mating designs notably increase the power for QTL detection. When power is intermediate, the power to detect a QTL using full-sib mating is, on average, 7% higher than under random mating. In addition, the IBD-based method proposed in this paper can be used to combine data from multiple families. As a result, the estimated QTL parameters can be applied to a wide statistical inference space relating to the entire reference population. This revised version was published online in July 2006 with corrections to the Cover Date.  相似文献   

20.
玉米株高和穗位高的QTL定位   总被引:20,自引:0,他引:20  
杨晓军  路明  张世煌  周芳  曲延英  谢传晓 《遗传》2008,30(11):1477-1486
摘要: 用玉米自交系掖478和丹340构建了397个F2:3家系群体, 利用双亲间多态的150个共显性SSR标记绘制分子连锁图谱, 图谱总长度1 478.7 cM, 标记间平均距离10.0 cM。在5种环境下对株高和穗位高性状进行鉴定, 复合区间作图法检测到21个株高QTL和25个穗位高QTL。于第1和5染色体的umc2025-umc1035及umc1822-bnlg1118区域检测到平均贡献率分别为12.2%和14.9%的株高QTL。于第3和5染色体的phi029-umc1102及phi109188-bnlg1118区域检测到平均贡献率达到10.2%和22.8%的穗位高QTL。第5染色体的Bin5.05-5.07区域可能存在控制株高和穗位高的主效QTL。株高和穗位高的基因作用方式主要是加性和部分显性效应。文章还分析了群体大小及试验环境对株高和穗位高QTL定位结果的影响  相似文献   

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

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