首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Mapping the genetic architecture of forest tree traits is important in order to understand the evolutionary forces that have shaped these traits and to facilitate the development of genomic-based breeding strategies. We examined the number, size, and distribution of allelic effects influencing eight types of traits using 30 published mapping studies (linkage and association mapping) in forest trees. The sizes of allelic effects, measured as the phenotypic variance explained, generally showed a severely right-skewed distribution. We estimated the numbers of underlying causal effects (n qtl ) for different trait categories by improving a method previously developed by Otto and Jones (Genetics 156:2093–2107, 2000). Estimates of n qtl based on association mapping studies were generally higher (median at 643) than those based on linkage mapping (median at 33). Comparisons with simulated linkage and association mapping data suggested that the lower n qtl estimates for the linkage mapping studies could partly be explained by fewer causal loci segregating within the full-sib family populations normally used, but also by the cosegregation of causal loci due to limited recombination. Disease resistance estimates based on linkage mapping studies had the lowest median of four underlying effects, while growth traits based on association mapping had about 580 effects. Theoretically, the capture of 50% of the genetic variation would thus require a population size of about 200 for disease resistance in linkage mapping, while growth traits in association mapping would require about 25,000. The adequacy and reliability of the improved method was successfully verified by applying it to the simulated data.  相似文献   

2.
Genetic linkage maps are indispensable tools in genetic and genomic studies. Recent development of genotyping-by-sequencing (GBS) methods holds great promise for constructing high-resolution linkage maps in organisms lacking extensive genomic resources. In the present study, linkage mapping was conducted for a bivalve mollusc (Chlamys farreri) using a newly developed GBS method—2b-restriction site-associated DNA (2b-RAD). Genome survey sequencing was performed to generate a preliminary reference genome that was utilized to facilitate linkage and quantitative trait locus (QTL) mapping in C. farreri. A high-resolution linkage map was constructed with a marker density (3806) that has, to our knowledge, never been achieved in any other molluscs. The linkage map covered nearly the whole genome (99.5%) with a resolution of 0.41 cM. QTL mapping and association analysis congruously revealed two growth-related QTLs and one potential sex-determination region. An important candidate QTL gene named PROP1, which functions in the regulation of growth hormone production in vertebrates, was identified from the growth-related QTL region detected on the linkage group LG3. We demonstrate that this linkage map can serve as an important platform for improving genome assembly and unifying multiple genomic resources. Our study, therefore, exemplifies how to build up an integrative genomic framework in a non-model organism.  相似文献   

3.

Background

The main goal of our study was to investigate the implementation, prospects, and limits of marker imputation for quantitative genetic studies contrasting map-independent and map-dependent algorithms. We used a diversity panel consisting of 372 European elite wheat (Triticum aestivum L.) varieties, which had been genotyped with SNP arrays, and performed intensive simulation studies.

Results

Our results clearly showed that imputation accuracy was substantially higher for map-dependent compared to map-independent methods. The accuracy of marker imputation depended strongly on the linkage disequilibrium between the markers in the reference panel and the markers to be imputed. For the decay of linkage disequilibrium present in European wheat, we concluded that around 45,000 markers are needed for low cost, low-density marker profiling. This will facilitate high imputation accuracy, also for rare alleles. Genomic selection and diversity studies profited only marginally from imputing missing values. In contrast, the power of association mapping increased substantially when missing values were imputed.

Conclusions

Imputing missing values is especially of interest for an economic implementation of association mapping in breeding populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1366-y) contains supplementary material, which is available to authorized users.  相似文献   

4.
Head smut, caused by the fungus Sphacelotheca reiliana (Kühn) Clint, is a devastating threat to maize production. In this study, QTL mapping of head smut resistance was performed using a recombinant inbred line (RIL) population from a cross between a resistant line “QI319” and a susceptible line “Huangzaosi” (HZS) with a genetic map constructed from genotyping-by-sequencing (GBS) data and composed of 1638 bin markers. Two head smut resistance QTL were identified, located on Chromosome 2 (q2.09HR) and Chromosome 5 (q5.03HR), q2.09HR is co-localized with a previously reported QTL for head smut resistance, and the effect of q5.03HR has been validated in backcross populations. It was also observed that pyramiding the resistant alleles of both QTL enhanced the level of resistance to head smut. A genome-wide association study (GWAS) using 277 diverse inbred lines was processed to validate the mapped QTL and to identify additional head smut resistance associations. A total of 58 associated SNPs were detected, which were distributed in 31 independent regions. SNPs with significant association to head smut resistance were detected within the q2.09HR and q5.03HR regions, confirming the linkage mapping results. It was also observed that both additive and epistastic effects determine the genetic architecture of head smut resistance in maize. As shown in this study, the combined strategy of linkage mapping and association analysis is a powerful approach in QTL dissection for disease resistance in maize.  相似文献   

5.
Soybean white mold (SWM), caused by Sclerotinia sclerotiorum ((Lib.) W. Phillips), is currently considered to be the second most important cause of soybean yield loss due to disease. Research is needed to identify SWM‐resistant germplasm and gain a better understanding of the genetic and molecular basis of SWM resistance in soybean. Stem pigmentation after treatment with oxaloacetic acid is an effective indicator of resistance to SWM. A total of 128 recombinant inbred lines (RILs) derived from a cross of ‘Maple Arrow’ (partial resistant to SWM) and ‘Hefeng 25’ (susceptible) and 330 diverse soybean cultivars were screened for the soluble pigment concentration of their stems, which were treated with oxalic acid. Four quantitative trait loci (QTLs) underlying soluble pigment concentration were detected by linkage mapping of the RILs. Three hundred and thirty soybean cultivars were sequenced using the whole‐genome encompassing approach and 25 179 single‐nucleotide polymorphisms (SNPs) were detected for the fine mapping of SWM resistance genes by genome‐wide association studies. Three out of five SNP markers representing a linkage disequilibrium (LD) block and a single locus on chromosome 13 (Gm13) were significantly associated with the soluble pigment content of stems. Three more SNPs that represented three minor QTLs for the soluble pigment content of stems were identified on another three chromosomes by association mapping. A major locus with the largest effect on Gm13 was found both by linkage and association mapping. Four potential candidate genes involved in disease response or the anthocyanin biosynthesis pathway were identified at the locus near the significant SNPs (<60 kbp). The beneficial allele and candidate genes should be useful in soybean breeding for improving resistance to SWM.  相似文献   

6.
Over the last 20 years, studies of Caenorhabditis elegans natural diversity have demonstrated the power of quantitative genetic approaches to reveal the evolutionary, ecological, and genetic factors that shape traits. These studies complement the use of the laboratory-adapted strain N2 and enable additional discoveries not possible using only one genetic background. In this chapter, we describe how to perform quantitative genetic studies in Caenorhabditis, with an emphasis on C. elegans. These approaches use correlations between genotype and phenotype across populations of genetically diverse individuals to discover the genetic causes of phenotypic variation. We present methods that use linkage, near-isogenic lines, association, and bulk-segregant mapping, and we describe the advantages and disadvantages of each approach. The power of C. elegans quantitative genetic mapping is best shown in the ability to connect phenotypic differences to specific genes and variants. We will present methods to narrow genomic regions to candidate genes and then tests to identify the gene or variant involved in a quantitative trait. The same features that make C. elegans a preeminent experimental model animal contribute to its exceptional value as a tool to understand natural phenotypic variation.  相似文献   

7.
Turnip mosaic virus (TuMV) is the major virus infecting Brassica crops. A dominant gene, TuRB01, that confers extreme resistance to some isolates of TuMV on Brassica napus (oilseed rape), has been mapped genetically. The mapping employed a set of doubled-haploid lines extracted from a population used previously to develop a reference RFLP map of the B. napus genome. The positioning of TuRB01 on linkage group N6 of the B. napus A–genome indicated that the gene probably originated from Brassica rapa. Resistance phenotypes were confirmed by indirect plate-trapped antigen ELISA using a monoclonal antibody raised against TuMV. The specificity of TuRB01 was determined using a wide range of TuMV isolates, including representatives of the European and American/Taiwanese pathotyping systems. Some isolates of TuMV that did not normally infect B. napus plants possessing TuRB01 produced mutant viruses able to overcome the action of the resistance gene. TuRB01 is the first gene for host resistance to TuMV to be mapped in a Brassica crop. A second locus, TuRB02, that appeared to control the degree of susceptibility to the TuMV isolate CHN 1 in a quantitative manner, was identified on the C-genome linkage group N14. The mapping of other complementary genes and the selective combining of such genes, using marker-assisted breeding, will make durable resistance to TuMV a realisable breeding objective. Received: 14 December 1998 / Accepted: 10 April 1999  相似文献   

8.
There is currently tremendous interest in the possibility of using genome-wide association mapping to identify genes responsible for natural variation, particularly for human disease susceptibility. The model plant Arabidopsis thaliana is in many ways an ideal candidate for such studies, because it is a highly selfing hermaphrodite. As a result, the species largely exists as a collection of naturally occurring inbred lines, or accessions, which can be genotyped once and phenotyped repeatedly. Furthermore, linkage disequilibrium in such a species will be much more extensive than in a comparable outcrossing species. We tested the feasibility of genome-wide association mapping in A. thaliana by searching for associations with flowering time and pathogen resistance in a sample of 95 accessions for which genome-wide polymorphism data were available. In spite of an extremely high rate of false positives due to population structure, we were able to identify known major genes for all phenotypes tested, thus demonstrating the potential of genome-wide association mapping in A. thaliana and other species with similar patterns of variation. The rate of false positives differed strongly between traits, with more clinal traits showing the highest rate. However, the false positive rates were always substantial regardless of the trait, highlighting the necessity of an appropriate genomic control in association studies.  相似文献   

9.

Key message

Using association and linkage mapping, two SNP markers closely linked to the SBWMV resistance gene on chromosome 5D were identified and can be used to select the gene in breeding.

Abstract

Soil-borne wheat mosaic virus (SBWMV) disease is a serious viral disease of winter wheat growing areas worldwide. SBWMV infection can significantly reduce grain yield up to 80 %. Developing resistant wheat cultivars is the only feasible strategy to reduce the losses. In this study, wheat Infinium iSelect Beadchips with 9 K wheat SNPs were used to genotype an association mapping population of 205 wheat accessions. Six new SNPs from two genes were identified to be significantly associated with the gene for SBWMV resistance on chromosome 5D. The SNPs and Xgwm469, an SSR marker that has been reported to be associated with the gene, were mapped close to the gene using F6-derived recombinant inbred lines from the cross between a resistant parent ‘Heyne’ and a susceptible parent ‘Trego’. Two representative SNPs, wsnp_CAP11_c209_198467 and wsnp_JD_c4438_5568170, from the two linked genes in wheat were converted into KBioscience Competitive Allele-Specific Polymerase assays and can be easily used in marker-assisted selection to improve wheat resistance to SBWMV in breeding.  相似文献   

10.
Grain size is one of the key traits that determines the quality of Basmati rice from the consumers’ as well as the traders’ point of view. Though many genes governing grain size have been identified in indica and japonica, little work has been done in Basmati rice. The present study aims at dissection of a QTL region governing grain size traits in Basmati employing association and linkage mapping approaches. Association mapping revealed that three markers, i.e., RM 6024 (grain breadth), RM1237 and RM18582 (grain length-breadth ratio), which cover 889 kb in the targeted QTL region have been significantly associated with grain size traits. Using linkage mapping, the targeted QTL region has been further delimited to a physical distance of 268 kb that comprises 24 annotated genes. The gene expression analysis of parents, revealed 19 genes differentially expressing within the QTL. Of them, 15 genes showed high expression in Basmati370, while four were expressed in Jaya, and whereas five genes did not show any differential expression between parents. Among differentially expressed genes, a highly expressed gene in Basmati370, Os05g0374200 (Cytokinin dehydrogenase 1 precursor) seems to be involved in accumulation of cytokinins, thus affecting the grain size. Therefore, our findings demonstrated that by complimenting association and linkage mapping, it is likely to dissect a QTL governing grain size traits in Basmati rice and also the QTL could be a potential target for marker-assisted breeding and map-based cloning studies.  相似文献   

11.
A set of 148 modern spring barley cultivars was explored for the extent of linkage disequilibrium (LD) between genes governing traits and nearby marker alleles. Associations of agronomically relevant traits (days to heading, plant height), resistance traits (leaf rust, barley yellow dwarf virus (BYD)), and morphological traits (rachilla hair length, lodicule size) with AFLP markers and SSR markers were found. Known major genes and QTLs were confirmed, but also new putative QTLs were found. The LD mapping clearly indicated the common occurrence of Rph3, a gene for hypersensitivity resistance against Puccinia hordei, and also confirmed the QTL Rphq2 for prolonging latency period of P. hordei in seedlings. We also found strong indication for a hitherto not reported gene for resistance or tolerance to BYD on chromosome 2, linked to SSR marker HVM054. Our conclusion is that LD mapping is a valuable additional tool in the search for applicable marker associations with major genes and QTLs. Electronic supplementary material Electronic supplementary material is available for this article at and accessible for authorised users.  相似文献   

12.

Key Message

Rapid LD decay in wild emmer population from Israel allows high-resolution association mapping. Known and putative new stripe rust resistance genes were found.

Abstract

Genome-wide association mapping (GWAM) is becoming an important tool for the discovery and mapping of loci underlying trait variation in crops, but in the wild relatives of crops the use of GWAM has been limited. Critical factors for the use of GWAM are the levels of linkage disequilibrium (LD) and genetic diversity in mapped populations, particularly in those of self-pollinating species. Here, we report LD estimation in a population of 128 accessions of self-pollinating wild emmer, Triticum turgidum ssp. dicoccoides, the progenitor of cultivated wheat, collected in Israel. LD decayed fast along wild emmer chromosomes and reached the background level within 1 cM. We employed GWAM for the discovery and mapping of genes for resistance to three isolates of Puccinia striiformis, the causative agent of wheat stripe rust. The wild emmer population was genotyped with the wheat iSelect assay including 8643 gene-associated SNP markers (wheat 9K Infinium) of which 2,278 were polymorphic. The significance of association between stripe rust resistance and each of the polymorphic SNP was tested using mixed linear model implemented in EMMA software. The model produced satisfactory results and uncovered four significant associations on chromosome arms 1BS, 1BL and 3AL. The locus on 1BS was located in a region known to contain stripe rust resistance genes. These results show that GWAM is an effective strategy for gene discovery and mapping in wild emmer that will accelerate the utilization of this genetic resource in wheat breeding.  相似文献   

13.
Hybrid speciation represents a relatively rapid form of diversification. Early models of homoploid hybrid speciation suggested that reproductive isolation between the hybrid species and progenitors primarily resulted from karyotypic differences between the species. However, genic incompatibilities and ecological divergence may also be responsible for isolation. Iris nelsonii is an example of a homoploid hybrid species that is likely isolated from its progenitors primarily by strong prezygotic isolation, including habitat divergence, floral isolation and post-pollination prezygotic barriers. Here, we used linkage mapping and quantitative trait locus (QTL) mapping approaches to investigate genomic collinearity and the genetic architecture of floral differences between I. nelsonii and one of its progenitor species I. hexagona. The linkage map produced from this cross is highly collinear with another linkage map produced between I. fulva and I. brevicaulis (the two other species shown to have contributed to the genomic makeup of I. nelsonii), suggesting that karyotypic differences do not contribute substantially to isolation in this homoploid hybrid species. Similar to other studies of the genetic architecture of floral characteristics, at least one QTL was found that explained >20% variance in each color trait, while minor QTLs were detected for each morphological trait. These QTLs will serve as hypotheses for regions under selection by pollinators.  相似文献   

14.
Culex pipiens pallens is the most abundant Culex mosquito species in northern China and is an important vector of bancroftian filariasis and West Nile virus. Deltamethrin is an insecticide that is widely used for mosquito control, however resistance to this and other insecticides has become a major challenge in the control of vector-borne diseases that appear to be inherited quantitatively. Furthermore, the genetic basis of insecticide resistance remains poorly understood. In this study, quantitative trait loci (QTL) mapping of resistance to deltamethrin was conducted in F2 intercross segregation populations using bulked segregation analysis (BSA) and amplified fragment length polymorphism markers (AFLP) in Culex pipiens pallens. A genetic linkage map covering 381 cM was constructed and a total of seven QTL responsible for resistance to deltamethrin were detected by composite interval mapping (CIM), which explained 95% of the phenotypic variance. The major QTL in linkage group 2 accounted for 62% of the variance and is worthy of further study. 12 AFLP markers in the map were cloned and the genomic locations of these marker sequences were determined by applying the Basic Local Alignment Search Tool (BLAST) tool to the genome sequence of the closely related Culex quinquefasciatus. Our results suggest that resistance to deltamethrin is a quantitative trait under the control of a major QTL in Culex pipiens pallens. Cloning of related AFLP markers confirm the potential utility for anchoring the genetic map to the physical map. The results provide insight into the genetic architecture of the trait.  相似文献   

15.
Due to its high nutritional value, chickpea is one of the most important and cost-effective legumes for human diet. Nutrient elements, such as Cu, P, K have numerous essential functions for the human metabolism. In this study, association mapping of loci linked to the seed Cu, P and K concentrations were performed on a population consisting of 107 Cicer reticulatum and 73 Cicer arietinum individuals in four environments (two locations x two years). A total of 121,840 SNPs were genotyped across 180 individuals by GBS analysis. The association mapping between the SNP markers and the seed Cu, P, K concentrations were identified and eight SNPs were found to be significantly associated with variations in three nutrient elements in more than two environments This research suggests that association mapping is a useful methodology for the identification of loci controlling the Cu, P and K uptake in chickpea seeds for further association mapping, molecular breeding, and marker-assisted selection and plant breeding studies and provides a broader understanding of the relationship between the investigated Cicer species and the effects of environmental conditions.  相似文献   

16.
Mungbean yellow mosaic India virus (MYMIV) is a major constraint on mungbean production in South and Southeast Asia. The virus belongs to the genus Begomovirus, causing yellow mosaic disease and subsequently yield loss up to 75–100 %. The present study employed F2 and BC1F1 populations derived from a cross between susceptible (BARImung 1; BM1) and resistant (BARImung 6; BM6) mungbeans to identify quantitative trait loci (QTLs) associated with resistance to MYMIV. Resistance to the virus was evaluated using F2:3 and BC1F1:2 populations under field conditions in two locations in Bangladesh in 2012. A total of 1,165 simple sequence repeat markers from different legumes were used to detect the polymorphism between BM1 and BM6. Sixty-one polymorphic markers were used to construct a linkage map comprising 11 linkage groups. Composite interval mapping consistently identified two major QTLs, qMYMIV2 on linkage group 2 and qMYMIV7 on linkage group 7, conferring the resistance in both F2 and BC1F1 populations. qMYMIV2 and qMYMIV7 accounted for 31.42–37.60 and 29.07–47.36 %, respectively, of the disease score variation, depending on populations and locations. At both loci, the resistant alleles were contributed by the parent BM6. qMYMIV2 appeared to be common to a major QTL for MYMIV resistance in mungbean reported previously, while qMYMIV7 is a new QTL for the resistance. The markers linked to the QTLs in this study are useful in marker-assisted breeding for development of mungbean varieties resistant to MYMIV.  相似文献   

17.

Background

Sudden death syndrome (SDS) is a serious threat to soybean production that can be managed with host plant resistance. To dissect the genetic architecture of quantitative resistance to the disease in soybean, two independent association panels of elite soybean cultivars, consisting of 392 and 300 unique accessions, respectively, were evaluated for SDS resistance in multiple environments and years. The two association panels were genotyped with 52,041 and 5,361 single nucleotide polymorphisms (SNPs), respectively. Genome-wide association mapping was carried out using a mixed linear model that accounted for population structure and cryptic relatedness.

Result

A total of 20 loci underlying SDS resistance were identified in the two independent studies, including 7 loci localized in previously mapped QTL intervals and 13 novel loci. One strong peak of association on chromosome 18, associated with all disease assessment criteria across the two panels, spanned a physical region of 1.2 Mb around a previously cloned SDS resistance gene (GmRLK18-1) in locus Rfs2. An additional variant independently associated with SDS resistance was also found in this genomic region. Other peaks were within, or close to, sequences annotated as homologous to genes previously shown to be involved in plant disease resistance. The identified loci explained an average of 54.5% of the phenotypic variance measured by different disease assessment criteria.

Conclusions

This study identified multiple novel loci and refined the map locations of known loci related to SDS resistance. These insights into the genetic basis of SDS resistance can now be used to further enhance durable resistance to SDS in soybean. Additionally, the associations identified here provide a basis for further efforts to pinpoint causal variants and to clarify how the implicated genes affect SDS resistance in soybean.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-809) contains supplementary material, which is available to authorized users.  相似文献   

18.
Different hypotheses on the genetic control of the resistance to the plum pox virus (PPV) have been reported in apricot, but there was a lack of agreement about the number of loci involved. In recent years, apricot genetic maps have been constructed from progenies derived from ‘Stark Early Orange’ or ‘Goldrich’, two main sources of resistance, three of these including the mapping of the PPV resistance loci. As the location of the locus was not precisely established, we mapped the PPV resistance loci using interval mapping (IM), composite interval mapping (CIM), and the Kruskal–Wallis non-parametric test in the F1 progeny derived from a cross between the susceptible cv. ‘Polonais’ and ’Stark Early Orange’. Four genomic regions were identified as being involved in PPV resistance. One of these mapped to the upper region of linkage group 1 of ‘Stark Early Orange’, and accounted for 56% of the phenotypic variation. Its location was similar to the one previously identified in ‘Goldrich’ and Prunus davidiana. In addition, a gene strongly associated to these major quantitative trait loci (QTL) was found to be related to PPV infection. Two putative QTLs were detected on linkage groups 3 of ‘Polonais’ and 5 of both ‘Polonais’ and ‘Stark Early Orange’ with both parametric and non-parametric methods at logarithm of odds (LOD) scores slightly above the detection threshold. The last QTL was only detected in the early stage of the infection. PPV resistance is, thus, controlled by a major dominant factor located on linkage group 1. The hypothesis of recessive factors with lower effect is discussed.  相似文献   

19.
Field evolved resistance of insect populations to Bacillus thuringiensis (Bt) crystalline (Cry) toxins expressed by crop plants has resulted in reduced control of insect feeding damage to field crops, and threatens the sustainability of Bt transgenic technologies. A single quantitative trait locus (QTL) that determines resistance in Ostrinia nubilalis larvae capable of surviving on reproductive stage transgenic corn that express the Bt Cry1Fa toxin was previously mapped to linkage group 12 (LG12) in a backcross pedigree. Fine mapping with high-throughput single nucleotide polymorphism (SNP) anchor markers, a candidate ABC transporter (abcc2) marker, and de novo mutations predicted from a genotyping-by-sequencing (GBS) data redefined a 268.8 cM LG12. The single QTL on LG12 spanned an approximate 46.1 cM region, in which marker 02302.286 and abcc2 were ≤2.81 cM, and the GBS marker 697 was an estimated 1.89 cM distant from the causal genetic factor. This positional mapping data showed that an O. nubilalis genome region encoding an abcc2 transporter is in proximity to a single QTL involved in the inheritance of Cry1F resistance, and will assist in the future identification the mutation(s) involved with this phenotype.  相似文献   

20.
Viruses vectored by the aphid Amphorophora agathonica cause decline in black raspberry plant health resulting in a shortened life and poor fruit quality of the infected plantings. New aphid resistant cultivars could increase the longevity of plantings providing growers and processors with consistent fruit production. Recent exploration of the native range of black raspberry identified three sources of aphid resistance: Ag4 from Ontario (ON), Canada, Ag5 from Maine (ME), and a third source from Michigan (MI) with no formal designation. The objectives of this study were to assess segregation of these three sources of aphid resistance in populations with single and combined sources and develop markers that can identify each source of resistance. A genetic linkage map constructed for ORUS 4305 placed the ON aphid resistance locus on Rubus linkage group (RLG) 6. Segregation ratios in populations with single and combined sources, and linkage mapping in two populations (ORUS 4304 and ORUS 4812) segregating for the Ag5 and MI sources, respectively, indicated that these three sources of resistance are each conferred by single dominant genes/alleles that are linked on RLG6. Confirmation of marker association in 16 validation populations identified four markers that could be used to predict resistance; however, none could distinguish between the ON and MI sources. These four markers may be useful for screening populations to enrich the field-planted progeny for aphid resistance. Fine mapping of the resistance loci is needed to develop functional markers at each of the resistance loci to enable pyramiding and durable aphid resistance.  相似文献   

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

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