首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
乳头数性状是猪最重要的繁殖性状之一,与养猪业的经济效益密切相关。本试验以Illunima Porcine SNP 60K芯片对22头可乐猪进行基因分型,通过Plink 1.07软件,基于混合线性模型对左乳头数、右乳头数和总乳头数性状进行全基因组关联分析。经过严格的质量控制和多重检验之后,共鉴定出全基因组水平潜在关联的SNPs位点4个(P<2.06E-5);染色体水平显著关联位点3个,潜在关联位点18个;在关联SNP位点可能连锁的上、下游1 cM区间内检索到304个基因,其中Wnt及Fgf信号通路中的候选基因BTRC、FGF5、FGF8和BMP3、RASGEF1B、HMGB3可能影响猪的乳头数性状或繁殖性状。通过全基因组关联分析策略发掘出的关联SNP位点及潜在候选基因,将为可乐猪的选育保种奠定基础。  相似文献   

2.
旨在鉴定荣昌猪初产繁殖性状的重要变异位点和基因,为荣昌猪繁殖性状的遗传改良提供重要的分子标记和基因资源。本研究选取429头荣昌母猪进行猪50K芯片基因分型,经过质量控制和基因型填充后,保留35 046个SNPs用于分析。采用主成分分析法研究群体结构,利用混合线性模型(mixed-linear model, MLM)将出生年、出生月作为固定效应,将主成分值作为协变量对总产仔数、活产仔数、死胎数和初生窝重性状进行全基因组关联分析(GWAS)。结果显示,在全基因组显著水平上鉴定出2个影响荣昌猪初生窝重的SNPs和1个影响荣昌猪死胎数的SNP;在潜在显著水平上鉴定到5个影响荣昌猪总产仔数的SNPs, 3个影响荣昌猪活产仔数的SNPs和10个影响荣昌猪死胎数的SNPs。通过全基因组关联分析筛选到1个显著的SNP(SSC17:57 315 180 bp)同时影响荣昌猪总产仔数、活产仔数和初生窝重,1个显著的SNP(SSC1:279 214 647 bp)同时影响荣昌猪活产仔数和总产仔数,暗示基因在不同性状间具有一因多效性。本研究根据候选基因的相关分子生物学功能,确定BMP7基因为影响荣昌猪总产仔数...  相似文献   

3.
识别野生动物群体内潜在影响动物表型变异的相关基因是进化遗传学研究的主要目的,而动物毛色是研究动物被毛表型形成遗传机制的最佳模型之一。应用Illumina公司提供的猪60 k SNP基因芯片对选取的62只不同被毛表型的野猪个体进行基因分型,利用SNP分析结果,通过对照全基因组关联分析(GWAS)识别影响野猪被毛表型差异的相关变异。结果表明,识别了6个与野猪被毛表型相关的基因组变异区域,分别位于SSC1(ALGA0001794,ASGA0006416)、SSC2(ASGA0011559)、SSC6(H3GA0018683)、SSC7(ASGA0035535)和SSC14(ASGA0060641);最显著相关的SNP(ALGA0001794)位于猪1号染色体上(SSC1)的27 899 596-27 899 696 bp区间(P=2.96×10-(-5))。该研究初步鉴定了6个与野猪毛色性状相关的易感位点,为进一步研究野猪不同毛色性状的形成机制提供了基础。  相似文献   

4.
旨在通过对奶绵羊产奶性状的全基因组关联分析(GWAS),寻找和定位与奶绵羊产奶性状相关的遗传标记和功能基因。本研究以135只戴瑞奶绵羊为试验材料,基于全基因组重测序技术,通过SAMTOOLS进行单核苷酸多态性(SNP)检测,使用PLINK v1.90进行质控,利用GEMMA v 0.98.1的混合线性模型对质控结果进行奶绵羊产奶性状相关的全基因组关联分析。结果表明,有1个SNP与泌乳后90天日均产奶量在全基因组范围内显著相关,8个SNPs达到潜在显著关联,相关的候选基因包括TRNAQ-CUG-2、LOC114117240、ACADLMYL1、CHD6、SLCO3A1;有2个SNPs与150天日均产奶量达潜在显著关联,相关的候选基因包括PRMT6、RNF180;有2个SNPs与泌乳周期达潜在显著关联,相关的候选基因包括PRMT6、TRNAW-CCA-68、TRNAS-GGA-61。进一步基因功能分析推测,ACADLSLCO3A1可能是影响奶绵羊产奶性状的候选基因。本研究为奶绵羊产奶性状的分子机制解析提供了一定的基础,为我国奶绵羊新品种培育提供了一定的理论参考。  相似文献   

5.
In order to identify the molecular markers related to alive litter size of Bama Xiang pigs,the genome-wide association study (GWAS) was used to map and screen the candidate genes affecting the alive litter size trait.Ear tissue samples of 297 Bama Xiang pigs with multiple parity records were collected,and DNA was extracted and genotyped by porcine 50K SNP beadchip.After quality control and genotype imputation,the alive litter size of Bama Xiang pigs were GWAS by Tassel.The results showed that the average number born alive per litter of Bama Xiang pigs increased gradually with the increasing of parity in the range of 1-9 parities.A total of 32 816 SNPs were obtained after quality control and filtration.8 SNPs related to alive litter size of Bama Xiang pigs were screened by genome-wide association analysis,which were significant at genome or chromosome level.Based on the enrichment analysis of the coding genes in the region between 500 kb upstream and downstream of the associated significant SNP loci,and the QTL regions and gene functions related to porcine reproductive traits,4 genes (CAPZB,MSH3,CITED2 and HSD17B7) were finally identified to be candidate genes related to alive litter size of Bama Xiang pigs.  相似文献   

6.
为寻找与巴马香猪产活仔数相关的分子标记,试验利用全基因组关联分析(GWAS)定位并筛选了影响产活仔数性状的候选基因,采集297头具有多胎产仔记录的巴马香猪耳组织样品,提取DNA并利用猪50K SNP芯片进行基因分型,分型结果经质控与基因型填充后,使用Tassel软件对巴马香猪产活仔数性状进行全基因组关联分析。结果显示,巴马香猪平均窝产活仔数在1~9胎内随着胎次增加逐渐升高;经质控过滤后共获得32 816个SNPs位点,利用全基因组关联分析共筛选到8个与巴马香猪产活仔数相关的SNPs位点,分别在基因组或染色体水平达到显著;对关联显著SNP位点上下游500 kb内的编码基因进行富集分析,并依据猪繁殖性状相关QTL区域及基因功能,最终筛选到4个基因(CAPZB、MSH3、CITED2和HSD17B7)作为影响巴马香猪产活仔数的候选基因。  相似文献   

7.
本实验旨在探究猪丝裂原活化蛋白激酶激酶6(Mitogen-Activated Protein Kinase Kinase 6,MAP2K6)多态性与猪重要经济性状的关系,为猪育种提供新的标记资源。实验选用法系大白猪为研究对象,采用质谱技术对猪MAP2K6基因的SNPs位点进行基因分型,同时利用线性混合模型分析单个标记位点与猪重要经济性状的相关性,共检测出4个SNPs,通过与生长育肥和乳头数性状的关联分析表明,rs345304630(SNP1)、rs325278117(SNP2)和rs332017877(SNP3)与达100 kg体重日龄和右乳头数显著相关,rs325752048(SNP4)与眼肌面积显著相关;连锁不平衡结果表明,MAP2K6基因SNP1、SNP2和SNP3处于强连锁不平衡状态。结果提示,在育种进程中,通过对优质基因型的选择来缩短母猪达100 kg体重日龄,增加母猪的泌乳和带仔能力,对提高哺乳仔猪的存活率和断奶重具有重要意义,也可为猪生长育肥和乳头数性状的分子标记的挖掘和利用提供一定的理论依据。  相似文献   

8.
单倍型标记与数量性状基因座(quantitative trait loci,QTL)之间具有较强的连锁不平衡(linkage disequilibrium,LD)关系,在基因定位和因果突变鉴定方面具有较高的应用价值。为了评估单倍型标记在基因组研究中的作用,本研究在华西牛资源群体中,选取该群体于2008—2021年间屠宰的共计1 478头平均月龄为24个月的个体进行研究,其中公牛1 333头,母牛145头。利用770K高密度芯片数据,基于LD阈值(r2>0.3)及固定单核苷酸多态(single nucleotide polymorphism,SNP)个数(5个连续SNP)两种方法进行单倍型构建,分别采用单位点SNP标记和两种单倍型标记共3种标记,基于GCTA的混合线性模型(mixed linear model,MLM),开展宰前活重(LW)和屠宰率(DP)等屠宰性状的全基因组关联分析(genome-wide association study,GWAS),定位影响屠宰性状的显著(P<0.05) SNPs、单倍型块和候选基因,同时比较3种标记的GWAS结果,评估3种标记的优劣。结果显示,3种标记在全基因组范围内共找到16个的显著SNPs及单倍型区域,主要分布于1、5、6、14、16、17和28号染色体上,同时鉴定到FAM184B、PPM1K、LCORL、RIMS2等10个与屠宰性状相关的候选基因,其中,基于SNP标记方法鉴定到的3个候选基因,在利用基于单倍型标记的方法中也鉴定到,且单倍型鉴定到的显著性位点或区域大多位于基因内部。在两种单倍型构建方法中,与基于固定SNP个数构建单倍型进行GWAS相比,基于LD阈值的构建方法鉴定到了更多候选基因。本研究结果表明,以单倍型开展GWAS可以综合考虑SNP位点间连锁关系,能较好地揭示复杂性状的遗传结构。  相似文献   

9.
旨在揭示白羽肉鸡孵化性状的遗传基础。本试验以白羽肉鸡A、B两个公鸡群体为素材,测定A群体5个世代(556只)、B群体3个世代(398只)的40周龄种蛋受精率和孵化率,并利用55K SNP芯片对两个群体共954个健康个体进行基因分型。基于系谱和基因型信息计算的A群体的受精率遗传力分别为0.21和0.14。利用混合线性模型对受精率和孵化率进行全基因组关联分析(genome-wide association study, GWAS)。全基因组关联分析共关联到12个显著SNPs,分布在1、9、11、12、20号染色体上,其中有9个SNPs与受精率显著关联,有3个SNPs与孵化率显著关联,对A群体受精率高、低组各4只公鸡睾丸组织进行转录组测序,实时荧光定量PCR验证两个差异表达基因。全基因组关联分析筛选出9个与受精率相关的候选基因COPG1、ADAMTS18、FABP2、CDH8、SLCO4A1、ATP13A3、NELL2、VEZT和IGHMBP2;3个与孵化率相关的候选基因TMEM、SCO1和MYH1A。转录组测序与荧光定量PCR验证了两个与受精率相关的候选基因HMGCLL1和COA6。本研究...  相似文献   

10.
Genome-wide association studies (GWAS) have introduced an influential tool in the search for quantitative trait loci (QTL) influencing economically important traits in sheep. To identify QTL associated with greasy fleece weight, a GWAS with 50 K single nucleotide polymorphisms (SNPs) was performed in a Baluchi sheep population. Association with greasy fleece weights was tested using the software Plink. The results of our GWAS provided three novel SNP markers and candidate genes associated with greasy fleece weight. A total of three chromosome-wide significant associations were detected for SNP on chromosomes 17 and 20 affecting greasy fleece weight across the four shearing. One of the significant SNP markers was located within ovine known genes namely FAM101A. Further investigation of these identified regions in validation studies will facilitate the identification of strong candidate genes for wool production in sheep.  相似文献   

11.
The aim of this study was to detect the copy number variation (CNV) in the genome of Bama Xiang pigs and investigate the effect of marker density on the efficiency and accuracy of CNV detection. PennCNV and R-Gada were employed to detect CNVs using the 1.4M high-density SNP chip data of 319 (160 hogs and 159 gilts) Bama Xiang pigs, and the CNV region (CNVR) was constructed by merging overlapping CNVs. Only the CNVR with higher frequency than 5% was verified by the genome-wide association study (GWAS). Finally, according to the marker densities, a certain number of SNPs were evenly extracted, and the effect of marker density on CNV detection efficiency and accuracy was explored. There were 6 327 CNVs detected by PennCNV and 3 489 CNVs detected by R-Gada, which made up of 795 and 340 CNVRs, respectively, including 226 CNVRs identified by both programs. Among the 226 CNVRs, the shortest was 3.98 kb, the longest was 1 297.78 kb, and their total length was 33.27 Mb, of which 102 (45%) overlapped the CNVRs reported previously. Among the 795 CNVRs detected by PennCNV, 135 had a higher frequency than 5%, 20 of which had been verified by GWAS, and the verification rate was 15%. With the SNP density increasing, the efficiency and accuracy of CNV detection were increased, especially for the small size CNVs. A CNVR sketch of Bama Xiang pigs had been drawn using 1.4M SNP chips, which was helpful to identify CNVRs associated with important economic traits in the future. At the same time, we revealed the positive effect of marker density on the efficiency and accuracy of CNV detection, and the results provided a reference of choosing marker density for the follow-up research of CNV detection.  相似文献   

12.
旨在检测巴马香猪基因组上的拷贝数变异(CNV),并探究标记密度对于CNV检测效率和准确率的影响。本研究利用319头巴马香猪(其中阉公猪160头和母猪159头)1.4M高密度SNP芯片的数据,采用PennCNV和R-Gada两种软件进行CNVs检测;然后通过重叠CNV融合法,构建拷贝数变异区域(CNVR),并用全基因组关联分析(GWAS)对频率大于5%的CNVR进行验证;最后根据不同的标记密度,均匀抽取一定数目的SNPs来探究标记密度对CNV检测效率和准确性的影响。结果,PennCNV和R-Gada软件分别检测到6 327和3 489个CNVs,分别构成795和340个CNVRs,其中226个为共同CNVRs。在这226个共同CNVRs中,最短的为3.98 kb,最长的为1 297.78 kb,总长度为33.27 Mb,其中102个(45%)与前人报道的CNVRs重叠。在PennCNV检出的795个CNVRs中,有135个频率大于5%,其中20个得到GWAS验证,验证率为15%。随着SNP密度的逐渐增加,CNV的检测效率和检测准确性不断提高,尤其是小片段CNVs的检测效率。本研究利用1.4M SNP芯片的数据,通过PennCNV和R-Gada软件绘制巴马香猪CNVR的草图,为将来鉴别与重要经济性状相关的CNVRs奠定了基础。同时,揭示了标记密度对CNV检测效率和准确性有正面影响,为后续CNV研究选择合适的标记密度提供了一定的参考。  相似文献   

13.
绵羊体重性状全基因组关联分析   总被引:1,自引:1,他引:0  
Motivated by mining major candidate genes across Ovine genome, the present study is to perform genome-wide association studies(GWAS) to detect genes associated with body weight traits. Using Illumina OvineSNP50 BeadChip, we performed a GWA study in 329 purebred sheep phenotyped for 6 body weight traits(birth weight, weaning weight, 6-month weight, pre-weaning gain, post-weaning gain, daily weight gain). Statistics and data analysis were based on TASSEL program,mixed linear model and the latest Ovis_aries_v3.1 genome sequence (released October 2012). The results indicated that 10 SNPs consistently reached genome-wise significant level for post-weaning gain and 22 SNPs reached chromosome-wise significant level for other body weight traits. The SNPs were within (MEF2B,RFXANK,et al) or close to some ovine genes, which were thought to be the most important candidate genes associated with body weight traits. The results will contribute to identify candidate genes for ovine body weight traits, and facilitate the potential utilization of genes involved production traits in sheep in future.  相似文献   

14.
The lipid content and fatty acid (FA) profile have an important impact in human health as well as in the technological transformation and nutritional and organoleptic quality of meat. A genome-wide association study (GWAS) on 144 backcross pigs (25% Iberian × 75% Landrace) was performed for 32 traits associated with intramuscular FA composition and indices of FA metabolism. The GWAS was carried out using Qxpak 5.0 and the genotyping information obtained from the Porcine SNP60K BeadChip (Illumina Inc., San Diego, CA). Signals of significant association considering a false- discovery rate (q-value < 0.05) were observed in 15 of the 32 analyzed traits, and a total of 813 trait-associated SNP (TAS), distributed in 43 chromosomal intervals on almost all autosomes, were annotated. According to the clustering analysis based on functional classification, several of the annotated genes are related to FA composition and lipid metabolism. Some interesting positional concordances among TAS and previously reported QTL for FA compositions and/or other lipid traits were also found. These common genomic regions for different traits suggest pleiotropic effects for FA composition and were found primarily on SSC4, SSC8, and SSC16. These results contribute to our understanding of the complex genetic basis of FA composition and FA metabolism.  相似文献   

15.
The racing line of Quarter Horses is characterized by great sprinting speed over short distances on straight tracks. To perform selection on racing horses, the speed index (SI) and conformation traits are often used. A genome-wide association study (GWAS) on 112 racing Quarter Horses was performed for the SI trait, and markers and genes associated were reported. The GWAS was carried out using the Qxpak.5 software and the genotyping data obtained from the Equine SNP50 BeadChip. A total of eight significant single-nucleotide polymorphisms (SNPs) (P < .0001; Q = 0.25) distributed on Equus caballus autosomal chromosomes 2, 4, 10, 18, and 27 were found and served as markers for genomic regions mined for candidate genes associated with SI. For candidate gene, annotation was considered 100 kb windows upstream and downstream to each important SNP. The highlighted genes were GRM8, GRIK2, NEB, ANK1, and KAT6A because their function could be related to racing performance. Future studies should consider a validation study with an independent population, and sequencing of these candidate genes should be done to identify causal alleles.  相似文献   

16.
【目的】 挖掘影响地方鸡体尺性状的有效SNP位点及功能基因, 给儋州鸡育种工作提供有效的数据基础和理论支撑。【方法】 共采集200只儋州鸡血样并提取基因组DNA, 利用10×全基因组重测序技术获得全基因组SNP标记并对试验个体基因型进行分型。使用EMMAX软件基于混合线性模型对70日龄的儋州鸡体尺性状(胫长、胫围、体斜长、胸宽、髋骨宽、胸深、龙骨长)进行全基因组关联分析。【结果】 共发现与胫长性状和胫围性状基因组水平显著相关的SNPs位点有12和8个, 与胫长性状相关SNPs分别定位于1、2、4和8号染色体上; 与胫围性状相关的SNPs定位于2、4、8和13号染色体上。预测与胫长相关的候选基因为KCNA1、TPK1、EZH2、FSTL5和AMY2A基因, 与胫围相关的候选基因为TPK1、FSTL5、AMY2ATGFBILECT2和IL-9。通过KEGG通路分析和GO注释发现, 8个基因参与钾离子跨膜转运、硫胺素新陈代谢、细胞增殖、钙离子结合、骨骼肌卫星细胞维持与骨骼肌再生、细胞受体相互作用、生长因子活性等生物学进程。【结论】 本研究发现了20个与儋州鸡体尺性状关联的SNPs位点, 并筛选到8个目标性状候选基因, 为儋州鸡育种提供候选的分子标记, 为地方鸡标记辅助选择提供新的思路。  相似文献   

17.
旨在通过对中国肉用西门塔尔牛纵向体重性状的全基因组关联分析(genome-wide association study, GWAS),定位与肉牛生长发育性状显著关联的候选基因。本研究利用808头中国肉用西门塔尔牛公牛0、6、12、18月龄的纵向体重数据,采用Gompertz、Logistic和Brody 3种非线性模型拟合个体的体重预测模型,估计参数A(成熟体重)、b(达到最大生长率的时间)和K(成熟率),然后以参数值为表型,BovineHD (770K) 芯片数据质控后剩余671 991 个SNPs,利用GAPIT进行关联分析,结合基因注释筛选影响肉牛发育的候选基因。选取拟合度最高的Gompertz模型(R2=0.954)确定相应参数估计值,GWAS共筛选到了9、49和7个显著的SNPs分别与A、b和K关联,且主要分布在2、3、7、9、11、14、22和25号染色体上。基因注释结果发现,PLIN3、KCNS3、ANGPTL2和ALPL与生长发育过程相关,其中KCNS3被认为是肌内脂肪含量的候选基因;IGF-1、TMCO1、PRKAG3和SHISA9影响肌肉发育过程,其中IGF-1被报道为生长发育过程的核心调控元件;ASPH基因参与调控肉牛胴体发育和肉质性状。本研究利用体重预测模型的参数估计值作为表型进行GWAS分析,定位到了一些与生长发育性状相关的候选基因,为其他纵向数据的研究提供了参考,也为调控肉牛生长发育进而提高产肉量的育种工作提供了新的候选分子标记。  相似文献   

18.
为鉴定与宁海黄鸡和广西黄鸡繁殖性状相关的分子标记及候选基因,本研究利用600K SNP芯片技术对鸡繁殖性状进行了全基因组关联分析(GWAS)。结果显示:4个单核苷酸多态性位点(SNPs)与繁殖性状的相关性达到了Bonferroni校正的5%全基因组显著性水平,其中,rs313221983与死胚蛋数相关,rs314844182与死胚蛋数及受精蛋孵化率相关,rs14758703与受精率相关,rs312721292与入孵蛋孵化率和受精蛋孵化率相关;在每个显著的SNP位点上下游5 kb范围内进行检测,共发现3个可能与死胚蛋数、受精率、入孵蛋孵化率和受精蛋孵化率相关近端基因,即唾液酸转移酶1(ST8SIA1)基因、神经外胚层皮质1(ENC1)基因和LOC101750905。本研究为地方品种鸡标记辅助选择和基因组选择提供了理论依据。  相似文献   

19.
To investigate single nucleotide polymorphism (SNP) loci associated with yearling wool traits of fine-wool sheep for optimizing marker-assisted selection and dissection of the genetic architecture of wool traits, we conducted a genome-wide association study (GWAS) based on the fixed and random model circulating probability unification (FarmCPU) for yearling staple length (YSL), yearling mean fiber diameter (YFD), yearling greasy fleece weight (YGFW), and yearling clean fleece rate (YCFR) by using the whole-genome re-sequenced data (totaling 577 sheep) from the following four fine-wool sheep breeds in China: Alpine Merino sheep (AMS), Chinese Merino sheep (CMS), Qinghai fine-wool sheep (QHS), and Aohan fine-wool sheep (AHS). A total of 16 SNPs were detected above the genome-wise significant threshold (P = 5.45E-09), and 79 SNPs were located above the suggestive significance threshold (P = 5.00E-07) from the GWAS results. For YFD and YGFW traits, 7 and 9 SNPs reached the genome-wise significance thresholds, whereas 10 and 12 SNPs reached the suggestive significance threshold, respectively. For YSL and YCFR traits, none of the SNPs reached the genome-wise significance thresholds, whereas 57 SNPs exceeded the suggestive significance threshold. We recorded 14 genes located at the region of ±50-kb near the genome-wise significant SNPs and 59 genes located at the region of ±50-kb near the suggestive significant SNPs. Meanwhile, we used the Average Information Restricted Maximum likelihood algorithm (AI-REML) in the “HIBLUP” package to estimate the heritability and variance components of the four desired yearling wool traits. The estimated heritability values (h2) of YSL, YFD, YGFW, and YCFR were 0.6208, 0.7460, 0.6758, and 0.5559, respectively. We noted that the genetic parameters in this study can be used for fine-wool sheep breeding. The newly detected significant SNPs and the newly identified candidate genes in this study would enhance our understanding of yearling wool formation, and significant SNPs can be applied to genome selection in fine-wool sheep breeding.  相似文献   

20.
试验旨在研究MAP3K5基因多态性与杜洛克猪生长、饲料利用性状的关联性。利用Illumina SNP60芯片测序结果比较发现了MAP3K5基因的4个SNPs位点(Ssc1:30769583 A>C;Ssc1:30781169 A>G;Ssc1:30940839 A>G;Ssc1:30962276 G>A)。统计获得MAP3K5基因的4个SNPs位点的基因型频率及等位基因频率,发现4个SNPs为低度-中度的多态突变位点。对MAP3K5基因多态性与生长、饲料利用性状的关联性进行分析,结果表明,杜洛克猪Ssc1:30769583 A>C位点CC基因型个体与AA基因型个体相比剩余采食量(RFI)性状降低了133.08 g/d(P<0.05);Ssc1:30781169 A>G位点的GG基因型个体与AG基因型个体相比RFI性状降低了116.18 g/d(P<0.05),ADFI性状降低了0.23 kg/d(P<0.05);Ssc1:30940839 A>G位点的GG基因型个体与AG基因型个体相比饲料转化率(FCR)性状降低了0.10%(P<0.05);Ssc1:30962276 G>A位点的AA基因型个体与GG基因型个体相比ADG性状降低了0.04 kg/d(P<0.05)。试验初步认为MAP3K5基因4个SNPs位点的多态性对杜洛克猪的RFI、ADFI、FCR和ADG具有一定影响。  相似文献   

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

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