首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Overlapping runs of homozygosity (ROH islands) shared by the majority of a population are hypothesized to be the result of selection around a target locus. In this study we investigated the impact of selection for coat color within the Noriker horse on autozygosity and ROH patterns. We analyzed overlapping homozygous regions (ROH islands) for gene content in fragments shared by more than 50% of horses. Long‐term assortative mating of chestnut horses and the small effective population size of leopard spotted and tobiano horses resulted in higher mean genome‐wide ROH coverage (SROH) within the range of 237.4–284.2 Mb, whereas for bay, black and roan horses, where rotation mating is commonly applied, lower autozygosity (SROH from 176.4–180.0 Mb) was determined. We identified seven common ROH islands considering all Noriker horses from our dataset. Specific islands were documented for chestnut, leopard spotted, roan and bay horses. The ROH islands contained, among others, genes associated with body size (ZFAT, LASP1 and LCORL/NCAPG), coat color (MC1R in chestnut and the factor PATN1 in leopard spotted horses) and morphogenesis (HOXB cluster in all color strains except leopard spotted horses). This study demonstrates that within a closed population sharing the same founders and ancestors, selection on a single phenotypic trait, in this case coat color, can result in genetic fragmentation affecting levels of autozygosity and distribution of ROH islands and enclosed gene content.  相似文献   

2.
Tomato (Solanum lycopersicum L.) has become a popular model for genetic studies of fruit flavor in the last two decades. In this article we present a study of tomato fruit flavor, including an analysis of the genetic, metabolic and sensorial variation of a collection of contemporary commercial glasshouse tomato cultivars, followed by a validation of the associations found by quantitative trait locus (QTL) analysis of representative biparental segregating populations. This led to the identification of the major sensorial and chemical components determining fruit flavor variation and detection of the underlying QTLs. The high representation of QTL haplotypes in the breeders’ germplasm suggests that there is great potential for applying these QTLs in current breeding programs aimed at improving tomato flavor. A QTL on chromosome 4 was found to affect the levels of the phenylalanine‐derived volatiles (PHEVs) 2‐phenylethanol, phenylacetaldehyde and 1‐nitro‐2‐phenylethane. Fruits of near‐isogenic lines contrasting for this locus and in the composition of PHEVs significantly differed in the perception of fruity and rose‐hip‐like aroma. The PHEV locus was fine mapped, which allowed for the identification of FLORAL4 as a candidate gene for PHEV regulation. Using a gene‐editing‐based (CRISPR‐CAS9) reverse‐genetics approach, FLORAL4 was demonstrated to be the key factor in this QTL affecting PHEV accumulation in tomato fruit.  相似文献   

3.
Multiparental cross designs for mapping quantitative trait loci (QTL) provide an efficient alternative to biparental populations because of their broader genetic basis and potentially higher mapping resolution. We describe the development and deployment of a recombinant inbred line (RIL) population in durum wheat (Triticum turgidum ssp. durum) obtained by crossing four elite cultivars. A linkage map spanning 2664 cM and including 7594 single nucleotide polymorphisms (SNPs) was produced by genotyping 338 RILs. QTL analysis was carried out by both interval mapping on founder haplotype probabilities and SNP bi‐allelic tests for heading date and maturity date, plant height and grain yield from four field experiments. Sixteen QTL were identified across environments and detection methods, including two yield QTL on chromosomes 2BL and 7AS, with the former mapped independently from the photoperiod response gene Ppd‐B1, while the latter overlapped with the vernalization locus VRN‐A3. Additionally, 21 QTL with environment‐specific effects were found. Our results indicated a prevalence of environment‐specific QTL with relatively small effect on the control of grain yield. For all traits, functionally different QTL alleles in terms of direction and size of genetic effect were distributed among parents. We showed that QTL results based on founder haplotypes closely matched functional alleles at known heading date loci. Despite the four founders, only 2.1 different functional haplotypes were estimated per QTL, on average. This durum wheat population provides a mapping resource for detailed genetic dissection of agronomic traits in an elite background typical of breeding programmes.  相似文献   

4.
The molecular analysis of genes influencing human height has been notoriously difficult. Genome-wide association studies (GWAS) for height in humans based on tens of thousands to hundreds of thousands of samples so far revealed ~200 loci for human height explaining only 20% of the heritability. In domestic animals isolated populations with a greatly reduced genetic heterogeneity facilitate a more efficient analysis of complex traits. We performed a genome-wide association study on 1,077 Franches-Montagnes (FM) horses using ~40,000 SNPs. Our study revealed two QTL for height at withers on chromosomes 3 and 9. The association signal on chromosome 3 is close to the LCORL/NCAPG genes. The association signal on chromosome 9 is close to the ZFAT gene. Both loci have already been shown to influence height in humans. Interestingly, there are very large intergenic regions at the association signals. The two detected QTL together explain ~18.2% of the heritable variation of height in horses. However, another large fraction of the variance for height in horses results from ECA 1 (11.0%), although the association analysis did not reveal significantly associated SNPs on this chromosome. The QTL region on ECA 3 associated with height at withers was also significantly associated with wither height, conformation of legs, ventral border of mandible, correctness of gaits, and expression of the head. The region on ECA 9 associated with height at withers was also associated with wither height, length of croup and length of back. In addition to these two QTL regions on ECA 3 and ECA 9 we detected another QTL on ECA 6 for correctness of gaits. Our study highlights the value of domestic animal populations for the genetic analysis of complex traits.  相似文献   

5.
A previous analysis of an F2/Backcross Charolais × Holstein cross population identified the presence of a highly significant QTL on chromosome 6 (BTA6) affecting the proportion of bone in the carcass. Two closely linked QTL affected birth weight (BW) and body length at birth (BBL). In this report, the marker density around the QTL on BTA6 was increased, adding four additional microsatellite markers across the chromosome and 46 SNPs within the target QTL confidence interval. Of the SNPs, 26 were in positional candidate genes and the remaining 20 provided an even distribution of markers in the target QTL region. As a bone‐related trait, the sum of the bone weight for all the left fore‐ and hindquarter joints of the carcass was analysed. We also studied the BW and BBL. Analyses of the data substantially reduced the QTL confidence interval. No strong evidence was found that the QTL for the three traits studied are different, and we conclude that the results are consistent with a single pleiotropic QTL influencing the three traits, with the largest effects on the proportion of bone in the carcass. The analyses also suggest that none of the SNPs tested is the sole causative variant of the QTL effects. Specifically, the SNP in the NCAPG gene previously reported as a causal mutation for foetal growth and carcass traits in other cattle populations was excluded as the causal mutation for the QTL reported here. Polymorphisms located in other previously identified candidate genes including SPP1, ABCG2, IBSP, MEPE and PPARGC1A were also excluded. The results suggest that SNP51_BTA‐119876 is the polymorphism in strongest linkage disequilibrium with the causal mutation(s). Further research is required to identify the causal variant(s) associated with this bone‐related QTL.  相似文献   

6.
Recently, the SERPINA6 gene encoding corticosteroid‐binding globulin (CBG) has been proposed as a candidate gene for a quantitative trait locus (QTL) affecting cortisol level on pig chromosome 7. The QTL was repeatedly detected in different lines, including a Piétrain × (German Landrace × German Large White) cross (PiF1) and purebred German Landrace (LR). In this study, we investigated whether the known non‐synonymous polymorphisms c.44G>T, c.622C>T, c.770C>T, c.793G>A, c.832G>A and c.919G>A of SERPINA6 are sufficient to explain the QTL in these two populations. Our investigations revealed that SNPs c.44G>T, c.622C>T, c.793G>A and c.919G>A are associated with cortisol level in PiF1 (< 0.01). Haplotype analysis showed that these associations are largely attributable to differences between a major haplotype carrying SNPs c.793G>A and c.919G>A and a haplotype carrying SNPs c.44G>T and c.622C>T. Furthermore, some SNPs, particularly c.44G>T and c.622C>T and the carrier haplotype, showed association with meat quality traits including pH and conductivity (< 0.05). In LR, the non‐synonymous SNPs segregate at very low frequency (<5%) and/or show only weak association with cortisol level (SNPs c.832G>A and c.919G>A; < 0.05). These findings suggest that the non‐synonymous SNPs are not sufficient to explain the QTL across different breeds. Therefore, we examined whether the expression of SERPINA6 is affected by cis‐regulatory polymorphisms in liver, the major organ for CBG production. We found allelic expression imbalance of SERPINA6, which suggests that its expression is indeed affected by genetic variation in cis‐acting elements. This represents candidate causal variation for future studies of the molecular background of the QTL.  相似文献   

7.
A dataset consisting of 787 animals with high‐density SNP chip genotypes (346 774 SNPs) and 939 animals with medium‐density SNP chip genotypes (33 828 SNPs) from eight indigenous Swiss sheep breeds was analyzed to characterize population structure, quantify genomic inbreeding based on runs of homozygosity and identify selection signatures. In concordance with the recent known history of these breeds, the highest genetic diversity was observed in Engadine Red sheep and the lowest in Valais Blacknose sheep. Correlation between FPED and FROH was around 0.50 and thereby lower than that found in similar studies in cattle. Mean FROH estimates from medium‐density data and HD data were highly correlated (0.95). Signatures of selection and candidate gene analysis revealed that the most prominent signatures of selection were found in the proximity of genes associated with body size (NCAPG, LCORL, LAP3, SPP1, PLAG1, ALOX12, TP53), litter size (SPP1), milk production (ABCG2, SPP1), coat color (KIT, ASIP, TBX3) and horn status (RXFP2). For the Valais Blacknose sheep, the private signatures in proximity of genes/QTL influencing body size, coat color and fatty acid composition were confirmed based on runs of homozygosity analysis. These private signatures underline the genetic uniqueness of the Valais Blacknose sheep breed. In conclusion, we identified differences in the genetic make‐up of Swiss sheep breeds, and we present relevant candidate genes responsible for breed differentiation in locally adapted breeds.  相似文献   

8.
Based on RNA-seq analysis, we recently found that the bovine NCAPG (non-SMC condensin I complex, subunit G) gene is differentially expressed during development of the longissimus muscle. In the present study, we validated this result and, using quantitative real-time PCR analysis, identified two adjacent genes, LCORL and DCAF16, that are more abundant in fetal muscle tissue; further analysis of tissue-specific expression patterns indicated high abundance of NCAPG in muscle. Since no polymorphisms were detected in a previous study of Qinchuan cattle, we extended our investigation to examine the occurrence of single-nucleotide polymorphisms (SNPs) in the NCAPG gene. Three SNPs, i.e., one located in the intron region (g47747: T > G), a synonymous mutation (g52535: A > G) and a missense mutation (g53208: T > G) that leads to a change in the amino acid of interest (pIle442Met), were detected in a population of Qinchuan beef cattle (n = 300). Association analysis showed that these SNPs were significantly associated with the growth traits of Qinchuan beef cattle. Our results indicate that the bovine NCAPG gene may be involved in the development of the longissimus muscle. These polymorphisms in the NCAPG gene may be useful for marker-assisted selection of optimal body size in Qinchuan beef cattle.  相似文献   

9.
Common wild rice (Oryza rufipogon Griff.) is an important genetic reservoir for rice improvement. We investigated a quantitative trait locus (QTL), qGP5‐1, which is related to plant height, leaf size and panicle architecture, using a set of introgression lines of O. rufipogon in the background of the Indica cultivar Guichao2 (Oryza sativa L.). We cloned and characterized qGP5‐1 and confirmed that the newly identified gene OsEBS (enhancing biomass and spikelet number) increased plant height, leaf size and spikelet number per panicle, leading to an increase in total grain yield per plant. Our results showed that the increased size of vegetative organs in OsEBS‐expressed plants was enormously caused by increasing cell number. Sequence alignment showed that OsEBS protein contains a region with high similarity to the N‐terminal conserved ATPase domain of Hsp70, but it lacks the C‐terminal regions of the peptide‐binding domain and the C‐terminal lid. More results indicated that OsEBS gene did not have typical characteristics of Hsp70 in this study. Furthermore, Arabidopsis (Arabidopsis thaliana) transformed with OsEBS showed a similar phenotype to OsEBS‐transgenic rice, indicating a conserved function of OsEBS among plant species. Together, we report the cloning and characterization of OsEBS, a new QTL that controls rice biomass and spikelet number, through map‐based cloning, and it may have utility in improving grain yield in rice.  相似文献   

10.
Cho IC  Park HB  Yoo CK  Lee GJ  Lim HT  Lee JB  Jung EJ  Ko MS  Lee JH  Jeon JT 《Animal genetics》2011,42(6):621-626
Haematological traits play important roles in disease resistance and defence functions. The objective of this study was to locate quantitative trait loci (QTL) and the associated positional candidate genes influencing haematological traits in an F2 intercross between Landrace and Korean native pigs. Eight blood‐related traits (six erythrocyte traits, one leucocyte trait and one platelet trait) were measured in 816 F2 progeny. All experimental animals were genotyped with 173 informative microsatellite markers located throughout the pig genome. We report that nine chromosomes harboured QTL for the baseline blood parameters: genomic regions on SSC 1, 4, 5, 6, 8, 9, 11, 13 and 17. Eight of twenty identified QTL reached genome‐wide significance. In addition, we evaluated the KIT locus, an obvious candidate gene locus affecting variation in blood‐related traits. Using dense single nucleotide polymorphism marker data on SSC 8 and the marker‐assisted association test, the strong association of the KIT locus with blood phenotypes was confirmed. In conclusion, our study identified both previously reported and novel QTL affecting baseline haematological parameters in pigs. Additionally, the positional candidate genes identified here could play an important role in elucidating the genetic architecture of haematological phenotype variation in swine and in humans.  相似文献   

11.
Plants employ a diverse set of defense mechanisms to mediate interactions with insects and fungi. These relationships can leave lasting impacts on host plant genome structure such as rapid expansion of gene families through tandem duplication. These genomic signatures provide important clues about the complexities of plant/biotic stress interactions and evolution. We used a pseudo‐backcross hybrid family to identify quantitative trait loci (QTL) controlling associations between Populus trees and several common Populus diseases and insects. Using whole‐genome sequences from each parent, we identified candidate genes that may mediate these interactions. Candidates were partially validated using mass spectrometry to identify corresponding QTL for defensive compounds. We detected significant QTL for two interacting fungal pathogens and three insects. The QTL intervals contained candidate genes potentially involved in physical and chemical mechanisms of host–plant resistance and susceptibility. In particular, we identified adjoining QTLs for a phenolic glycoside and Phyllocolpa sawfly abundance. There was also significant enrichment of recent tandem duplications in the genomic intervals of the native parent, but not the exotic parent. Tandem gene duplication may be an important mechanism for rapid response to biotic stressors, enabling trees with long juvenile periods to reach maturity despite many coevolving biotic stressors.  相似文献   

12.
Milk production traits, such as 305‐day milk yield (305MY), have been under direct selection to improve production in dairy cows. Over the past 50 years, the average milk yield has nearly doubled, and over 56% of the increase is attributable to genetic improvement. As such, additional improvements in milk yield are still possible as new loci are identified. The objectives of this study were to detect SNPs and gene sets associated with 305MY in order to identify new candidate genes contributing to variation in milk production. A population of 781 primiparous Holstein cows from six central Washington dairies with records of 305MY and energy corrected milk were used to perform a genome‐wide association analysis (GWAA) using the Illumina BovineHD BeadChip (777 962 SNPs) to identify QTL associated with 305MY (< 1.0 × 10?5). A gene set enrichment analysis with SNP data (GSEA‐SNP) was performed to identify gene sets (normalized enrichment score > 3.0) and leading edge genes (LEGs) influencing 305MY. The GWAA identified three QTL comprising 34 SNPs and 30 positional candidate genes. In the GSEA‐SNP, five gene sets with 58 unique and 24 shared LEGs contributed to 305MY. Identification of QTL and LEGs associated with 305MY can provide additional targets for genomic selection to continue to improve 305MY in dairy cattle.  相似文献   

13.
14.
PCR–RFLP was applied to analyse polymorphisms within the MRF4 and heart fatty acid-binding protein (H-FABP) gene for correlation studies with growth traits in three-month-old Qinchuan (QQ), Qinchuan × Limousin (LQ) and Qinchuan × Red Angus (AQ) cattle. The results showed that 874 bp PCR products of MRF4 digested with XbaI and 2,075 bp PCR products of H-FABP digested with HaeIII were polymorphic in the three populations. Moreover, the frequencies of allele A at MRF4 locus and allele B at H-FABP locus in the QQ, AQ, and LQ populations were 0.8358/0.8888/0.8273 and 0.8358/0.7500/0.8195 respectively. Allele A at MRF4 locus and allele B at H-FABP locus were dominant in the three populations. No statistically significant differences in growth traits were observed among the genotypes of the all three populations at H-FABP locus. However, the association of MRF4 polymorphism with growth traits was then determined in all three populations. The body weight, withers height, heart girth and height at hip cross of individuals with genotype AA were higher than those with genotype AB or BB (P < 0.05). Therefore, we suggest that the MRF4 gene may function in the control or expression of growth traits, particularly body weight, withers height, heart girth and height at hip cross.  相似文献   

15.
Mature body size is genetically correlated with growth rate, an important economic trait in the sheep industry. Mature body size has been studied extensively in humans as well as cattle and other domestic animal populations but not in sheep. Six‐hundred and sixteen ewes, across 22 breeds, were measured for 28 linear measurements representing various skeletal parts. PCA from these measures generated principal components 1 and 2 which represented 66 and 7% of the phenotypic variation respectively. Two‐hundred and twenty sheep were genotyped on the Illumina Ovine HD beadchip for a GWAS investigating mature body size and linear body measurements. Forty‐six (Bonferroni P < 0.05) SNP associations across 14 chromosomes were identified utilizing principal component 1, representing overall body size, revealing mature body size to have fewer loci of large effect than other domestic species such as dogs and horses. Genome‐wide associations for individual linear measures identified major quantitative trait loci for withers height and ear length. Withers height was associated (Bonferroni P < 0.05) with 12 SNPs across six chromosomes whereas ear length was associated with a single locus on chromosome 3, containing MSRB3. This analysis identified several loci known to be associated with mature body size in other species such as NCAPG, LCORL, and HMGA2. Mature body size is more polygenic in sheep than other domesticated species, making the development of genomic selection for the trait the most efficient option for maintaining or reducing mature body size in sheep.  相似文献   

16.
Genomic studies have been used to identify genes underlying many important plant secondary metabolic pathways. However, genes for salicinoid phenolic glycosides (SPGs)—ecologically important compounds with significant commercial, cultural, and medicinal applications—remain largely undescribed. We used a linkage map derived from a full‐sib population of hybrid cottonwoods (Populus spp.) to search for quantitative trait loci (QTL) for the SPGs salicortin and HCH‐salicortin. SSR markers and primer sequences were used to anchor the map to the V3.0 P. trichocarpa genome. We discovered 21 QTL for the two traits, including a major QTL for HCH‐salicortin (R2 = .52) that colocated with a QTL for salicortin on chr12. Using the V3.0 Populus genome sequence, we identified 2,983 annotated genes and 1,480 genes of unknown function within our QTL intervals. We note ten candidate genes of interest, including a BAHD‐type acyltransferase that has been potentially linked to PopulusSPGs. Our results complement other recent studies in Populus with implications for gene discovery and the evolution of defensive chemistry in a model genus. To our knowledge, this is the first study to use a full‐sib mapping population to identify QTL intervals and gene lists associated with SPGs.  相似文献   

17.
The impact of environmental factors, such as diet, and the genetic basis of autoimmune pancreatitis (AIP) are largely unknown. Here, we used an experimental murine AIP model to identify the contribution of diet to AIP development, as well as to fine‐map AIP‐associated genes in outbred mice prone to develop the disease. For this purpose, we fed mice of an autoimmune‐prone intercross line (AIL) three different diets (control, calorie‐reduced and western diet) for 6 months, at which point the mice were genotyped and phenotyped for AIP. Overall, 269 out of 734 mice (36.6%) developed AIP with signs of parenchymal destruction, equally affecting mice of both sexes. AIP prevalence and severity were reduced by approximately 50% in mice held under caloric restriction compared to those fed control or western diet. We identified a quantitative trait locus (QTL) on chromosome 4 to be associated with AIP, which is located within a previously reported QTL. This association does not change when considering diet or sex as an additional variable for the mapping. Using whole‐genome sequences of the AIL founder strains, we resolved this QTL to a single candidate gene, namely Map3k7. Expression of Map3k7 was largely restricted to islet cells as well as lymphocytes found in the exocrine pancreas of mice with AIP. Our studies suggest a major impact of diet on AIP. Furthermore, we identify Map3k7 as a novel susceptibility gene for experimental AIP. Both findings warrant clinical translation.  相似文献   

18.
Fruit neck length (FNL) is an important quality trait in cucumber because it directly affects its market value. However, its genetic basis remains largely unknown. We identified a candidate gene for FNL in cucumber using a next‐generation sequencing‐based bulked segregant analysis in F2 populations, derived from a cross between Jin5‐508 (long necked) and YN (short necked). A quantitative trait locus (QTL) on chromosome 7, Fnl7.1, was identified through a genome‐wide comparison of single nucleotide polymorphisms between long and short FNL F2 pools, and it was confirmed by traditional QTL mapping in multiple environments. Fine genetic mapping, sequences alignment and gene expression analysis revealed that CsFnl7.1 was the most likely candidate Fnl7.1 locus, which encodes a late embryogenesis abundant protein. The increased expression of CsFnl7.1 in long‐necked Jin5‐508 may be attributed to mutations in the promoter region upstream of the gene body. The function of CsFnl7.1 in FNL control was confirmed by its overexpression in transgenic cucumbers. CsFnl7.1 regulates fruit neck development by modulating cell expansion. Probably, this is achieved through the direct protein–protein interactions between CsFnl7.1 and a dynamin‐related protein CsDRP6 and a germin‐like protein CsGLP1. Geographical distribution differences of the FNL phenotype were found among the different cucumber types. The East Asian and Eurasian cucumber accessions were highly enriched with the long‐necked and short‐necked phenotypes, respectively. A further phylogenetic analysis revealed that the Fnl7.1 locus might have originated from India. Thus, these data support that the CsFnl7.1 has an important role in increasing cucumber FNL.  相似文献   

19.
A region on rat Chromosome (Chr) 2 of the Dahl salt-sensitive rat (S) was shown previously to contain a quantitative trait locus (QTL) for blood pressure (BP). This was achieved first by linkage, followed by the use of congenic strains. A congenic strain, designated S.MNS-D2Mit6/Adh, contained a segment of Chr 2 from the Milan Normotensive (MNS) rat in the S genetic background. Since the region containing the QTL was roughly 80 cM in size, a further reduction was needed towards the positional or candidate gene cloning. Currently, two congenic substrains were made from the original strain S.MNS-D2Mit6/Adh. One of these two substrains showed a BP-lowering effect, whereas the other substrain did not. Deducing the segment not shared in the two substrains, the BP QTL has to be present in a chromosome region of roughly 5.7 cM between the marker D2Rat303 and the locus for the neutroendopeptidase gene (Nep). Nep is not included within the segment. This region does not seem to contain any candidate genes well known for the BP control. Thus, the final identification of the QTL will most likely lead to the discovery of a brand new gene for the BP regulation. Received: 14 December 2000 / Accepted: 18 January 2001  相似文献   

20.
Better understanding of the mechanisms underlying interindividual variation in stress responses and their links with production traits is a key issue for sustainable animal breeding. In this study, we searched for quantitative trait loci (QTL) controlling the magnitude of the plasma cortisol stress response and compared them to body size traits in five F2 full‐sib families issued from two rainbow trout lines divergently selected for high or low post‐confinement plasma cortisol level. Approximately 1000 F2 individuals were individually tagged and exposed to two successive acute confinement challenges (1 month interval). Post‐stress plasma cortisol concentrations were determined for each fish. A medium density genome scan was carried out (268 markers, overall marker spacing less than 10 cM). QTL detection was performed using qtlmap software, based on an interval mapping method ( http://www.inra.fr/qtlmap ). Overall, QTL of medium individual effects on cortisol responsiveness (<10% of phenotypic variance) were detected on 18 chromosomes, strongly supporting the hypothesis that control of the trait is polygenic. Although a core array of QTL controlled cortisol concentrations at both challenges, several QTL seemed challenge specific, suggesting that responses to the first and to a subsequent exposure to the confinement stressor are distinct traits sharing only part of their genetic control. Chromosomal location of the steroidogenic acute regulatory protein (STAR) makes it a good potential candidate gene for one of the QTL. Finally, comparison of body size traits QTL (weight, length and body conformation) with cortisol‐associated QTL did not support evidence for negative genetic relationships between the two types of traits.  相似文献   

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

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