首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到16条相似文献,搜索用时 412 毫秒
1.
WGCNA鉴定奶山羊妊娠至泌乳期乳腺发育关键基因   总被引:2,自引:0,他引:2  
高慧杰  郑惠玲 《畜牧兽医学报》2020,51(11):2679-2688
旨在通过加权基因共表达网络分析(weighted gene co-expression network analysis,WGCNA)筛选奶山羊不同生理阶段乳腺发育的关键基因。本研究从GEO数据库中下载奶山羊不同生理阶段(妊娠46天、70天、90天、110天和产后40天)的乳腺组织微阵列数据集GSE14008,使用R语言的WGCNA包对数据进行共表达分析。将得到的模块与生理阶段进行关联分析,选择目标模块,并根据连接度选出枢纽基因。使用DAVID网站对模块进行富集分析后,使用String网站构建模块的蛋白互作网络,并使用Cytoscape软件得到核心基因,最终与枢纽基因取交集得到目标基因。对18个样本的8 443个基因进行加权基因共表达分析,得到30个模块,并选出4个与不同生理阶段相关的目标模块及每个模块的30个枢纽基因,同时也得到4个模块的蛋白互作网络及每个网络的20个核心基因。最终,4个模块共得到13个与乳腺发育相关的目标基因(UQCR、RGL2、NOTCH1、PTBP1、PPP5C、FZR1、UBE2L3、TNF、MAT2A、ITGB2、GPR18、JAK1、CSN2)。本研究通过WGCNA、GO富集分析和PPI网络等生物信息学技术,阐明了不同生理时期乳腺发育的关键过程,及在这些过程中起关键作用的基因,这为进一步研究乳腺发育机制提供了新的思路和线索。  相似文献   

2.
旨在基于RNA-Seq技术对塔里木马鹿毛色相关基因进行筛选及分析。采用Illumina Hi Seq TM2000测序平台对塔里木马鹿和天山马鹿的皮肤组织进行转录组测序,所得序列经质控、组装后比对到NR、Swiss-Prot、COG、KOG、KEGG、GO和Pfam数据库中注释,并对差异表达基因进行筛选、功能注释和富集分析。结果表明,测序获得25 038个有注释信息的Unigenes,比对分析显示,塔里木马鹿与天山马鹿有922个差异表达基因,其中上调表达基因495个,下调表达基因427个;GO功能富集分析结果显示,568个差异表达基因富集到61个GO条目上,分别参与了生物学过程、细胞组分及分子功能;KEGG代谢通路富集分析发现,在差异表达基因中富集最显著的代谢通路是ECM-受体相互作用。利用实时荧光定量PCR(qRT-PCR)方法分析与塔里木马鹿毛色相关的7个候选基因的转录水平变化来验证转录组测序结果的准确性和可靠性,这些基因的表达趋势与转录组测序结果相一致。ECM-受体相互作用、蛋白质消化与吸收、PI3K-Akt信号通路及与黑色素合成相关的酪氨酸等通路可能与塔里木马鹿的毛色有关;候选基因MITFGgt1、VDRPTPRFCⅡTAARPC5L、POMC等可能在塔里木马鹿毛色形成过程中发挥重要作用。本研究结果为今后塔里木马鹿毛色相关基因的分子调控机制方面及挖掘潜在的新基因提供了丰富的试验数据。  相似文献   

3.
旨在通过WGCNA技术解析转录组和代谢组数据,探究驴肉嫩度调控机制。试验动物采用24~36月龄的健康雌性广灵驴(平均体重236.10 kg),将驴肉剪切力和肌内脂肪含量作为表型数据,以每个样品3个重复进行表型数据测定。本试验基于前期研究的具有显著剪切力和肌内脂肪含量差异的14个广灵驴背最长肌样本的转录组和代谢组测序数据,运用WGCNA技术筛选与驴肉嫩度相关的基因及代谢物并进行转录组与代谢组联合分析,解析嫩度相关基因与代谢物。结果表明,利用WGCNA技术通过|r|≥0.5以及P≤0.05筛选标准得到3个与嫩度相关的关键基因模块Greenyellow、Darkgrey、Darkgreen以及2个关键代谢物模块Brown、Yellow。对关键基因模块进行GO富集分析发现,模块内基因主要在甘油磷脂的生物合成、脂质氧化、脂肪酸β-氧化、细胞大分子分解代谢过程、肌肉器官发育、钙离子结合等GO功能上富集。存在于关键模块内的基因及代谢物经KEGG功能富集分析发现,其大多集中在精氨酸和脯氨酸代谢、Wnt信号通路、蛋白质消化吸收、脂肪酸代谢、TCA循环、胰高血糖素信号通路、甘油磷脂代谢、嘌呤代谢、β-丙氨酸代谢等通路上。联合分析表明,丙氨酸、天冬氨酸和谷氨酸代谢、精氨酸和脯氨酸代谢、β-丙氨酸代谢以及PPAR信号通路可能调控驴肉嫩度。WGCNA及联合KEGG共富集分析筛选到的丙氨酸、天冬氨酸和谷氨酸代谢、精氨酸和脯氨酸代谢、β-丙氨酸代谢以及PPAR信号通路可能对调控驴肉嫩度有重要作用;而GAD1、PPATNIT2AGMATCARNS1、ACOXL以及腺苷酸基琥珀酸、L-脯氨酸、L-谷氨酸、肌酸、高肌肽、肌肽、泛酸、(9S)-羟基十八碳二烯酸则可能是影响驴肉嫩度的候选基因及代谢物。本试验可为今后广灵驴肉质嫩度的分子调控与改良育种提供一定的理论基础。  相似文献   

4.
中国美利奴羊胚胎骨骼肌发育的加权基因共表达网络分析   总被引:2,自引:1,他引:1  
旨在对绵羊胚胎骨骼肌lncRNAs(long non-coding RNA)进行鉴定分析,以阐明其在肌纤维类型转换与肌纤维增粗过程中的调控机制。本研究选取体重相近的成年中国美利奴母羊进行同期发情和人工授精,通过全转录组测序技术对其妊娠第85(D85N)、105(D105N)和135天(D135N)的胎儿背最长肌组织进行测序,设置D85N vs D105N、D105N vs D135N和D85N vs D135N 3个比较组,通过比较筛选出显著差异表达的lncRNA与mRNA。利用加权基因共表达网络分析(weighted gene co-expression network analysis,WGCNA)方法构建共表达模块,使用DAVID在线工具和R-package进行GO和KEGG富集分析以找到与肌肉发育相关的模块。最后从目标模块中筛选出高连通度的lncRNAs和mRNAs,通过它们与miRNAs间的靶向预测关系建立lncRNA-miRNA-mRNA共表达网络。根据WGCNA分析结果,共得到25个模块。功能富集显示,模块中的核心基因主要富集于细胞粘附、Wnt、紧密连接、mTOR、AMPK及ECM-受体相互作用等肌肉发育相关的信号通路,选出模块中连通度高的lncRNAs和mRNAs构建子网络,得到TNNI2、PIP5K1A、PDK4等关键相关基因,预测出MSTRG.3903、MSTRG.10154、MSTRG.1629、MSTRG.10496、MSTRG.9559、MSTRG.10178、MSTRG.10521、MSTRG.3911、MSTRG.4586、MSTRG.7232等10个与肌肉发育、肌肉疾病、细胞增殖相关的lncRNAs。本研究成功构建了肌纤维发育相关的lncRNA-miRNA-mRNA共表达网络,找到多个与妊娠后期胚胎骨骼肌发育相关的潜在候选基因,为深入研究lncRNA在中国美利奴羊胚胎发育过程中骨骼肌的发育调控机制奠定了基础,也为其他家畜骨骼肌发育机制的研究提供了参考和方向。  相似文献   

5.
对生长速度长期的高压选择导致了快大型肉鸡与优质型地方鸡显著的表型差异,本研究旨在研究这种表型差异的分子遗传机理。采用Agilent鸡全基因组表达谱芯片对达到性成熟的优质型清远麻鸡(112d)和快大型科宝肉鸡(42d)比目鱼肌中与肌纤维发育和类型组成相关的候选基因及信号通路进行系统筛查。结果,芯片分析共筛选到差异倍数在2倍及以上的基因1 318个,以科宝肉鸡作为参照,清远麻鸡中上调基因501个,下调基因817个,主要涉及到肌肉发育、能量代谢、脂质代谢等生物学过程。基于KEGG Pathway分析发现,差异基因除了参与肌纤维发育和分化相关信号通路(如Hedgehog信号通路和Ca2+信号通路)外,Wnt信号通路,mTOR信号通路,MAPK信号通路,ErbB信号通路、JAK-STAT信号通路等一些跟能量代谢相关的信号通路也被富集为显著的信号通路。与能量代谢相关的信号通路和与肌肉发育相关的信号通路相互作用形成一个调控网络从而影响肌纤维的发育,筛选出了20个在肌纤维生长发育及分化过程中可能具有重要影响的候选基因,但这些差异表达基因在肌纤维的发育和分化过程中的作用尚待深入研究。  相似文献   

6.
【目的】通过对1和60日龄广西麻鸡腿肌进行转录组测序(RNA-Seq),筛选出广西麻鸡生长发育关键基因。【方法】选择1和60日龄健康母鸡各3只,分别采集腿肌组织,利用Illumina NovaseqTM 6000平台进行mRNA转录组测序,利用DESeq2软件进行差异表达基因的分析,并对差异表达基因进行GO功能和KEGG通路富集分析,进一步筛选出生长发育相关基因,并用实时荧光定量PCR验证转录组测序结果。【结果】转录组测序结果显示,6个样本分别获得34 296 198~41 647 642条clean reads。与1日龄腿肌相比,60日龄腿肌共获得2 304个差异表达基因,其中998个基因上调,1 306个基因下调。GO功能富集分析发现,共获得富集条目572条,其中生物过程381条,细胞组分70条,分子功能121条。KEGG通路富集分析发现,共获得34个显著富集的信号通路,其中心肌收缩、MAPK信号通路、紧密连接、肌动蛋白细胞骨架信号通路、心肌细胞中的肾上腺素能信号传导等与生长发育相关。通过GO功能和KEGG通路富集分析获得5个生长发育相关基因,分别为肌球蛋白重链10(MYH10)、肌球蛋白链15(MYH15)、成纤维细胞生长因子10(FGF10)、成纤维细胞生长因子16(FGF16)、肌肉生长抑制素(GDF8)基因,其中MYH10、MYH15、FGF10为上调差异表达基因,FGF16和GDF8为下调差异表达基因。所选差异表达基因实时荧光定量PCR结果与转录组测序的基因表达水平一致,说明转录组测序结果可靠。【结论】本试验通过对广西麻鸡1和60日龄2个不同生长阶段进行转录组测序分析,获得MYH10、MYH15、FGF10、FGF16、GDF8共5个与生长发育相关的关键基因,为后续进一步探讨广西麻鸡生长发育的分子调控机制提供参考。  相似文献   

7.
本研究对课题组前期在京星黄鸡转录组研究中筛选到的与肌内脂肪(IMF)差异沉积相关的14个候选基因进行验证,检测其在中外两个鸡种群体中与胸肌IMF沉积的关联性。以98日龄慢速型地方鸡京星黄鸡和42日龄快速型白羽科宝肉鸡胸肌组织为素材,通过胸肌甘油三酯(TG)含量区分高低表型组,并检测候选基因在组间的基因表达差异。结果表明,在京星黄鸡胸肌TG高、低组间ATP结合盒亚家族B成员8(ABCB8)、脂联素(ADIPOQ)、第6号染色体开放阅读框65(BEND6)、CD74分子(CD74)、核糖基5-磷酸转移酶(FKTN)、组蛋白乙酰转移酶1(HAT1)、硫酸乙酰肝素-氨基葡萄糖3-磺基转移酶5(HS3ST5)、介体复合物亚基4(MED4)、肿瘤坏死因子超家族成员8(TNFSF8)和TNFAIP3相互作用蛋白1(TNIP1)共10个基因表达差异显著(P<0.05);在科宝肉鸡TG高、低组间ADIPOQ、BEND6、FKTN、HAT1、HS3ST5、MED4和TNIP1共7个基因表达差异显著(P<0.05);ADIPOQ、FKTN、HAT1、HS3ST5、MED4和TNIP1共6个基因在两个品种中差异表达趋势一致(P<0.05)。本研究提供了鸡IMF沉积相关新候选基因,为IMF分子调控机理研究和相关分子标记筛选研究奠定了良好的基础。  相似文献   

8.
旨在探究中卫山羊羊毛弯曲随生长发生变化的分子机制,挖掘不同发育时期影响羊毛弯曲度的关键基因。试验采集了3只中卫山羊出生后45日龄(弯曲毛)和365日龄(直毛)的皮肤组织,进行转录组测序(RNA-seq),使用WGCNA与GSEA分析找出Hub基因。以P-value<0.05和|log2FoldChange|≥1作为差异表达基因筛选的标准,45日龄为对照组,365日龄为试验组共筛选得到1 252个显著差异表达基因,包括812个表达上调基因和440个下调基因。基于WGCNA方法分析共得到14个模块,其中黄色和绿松石模块与被毛弯曲表型相关。利用String和Cytoscape构建网络,筛选到20个影响羊毛弯曲度的核心基因。从GSEA结果中筛选的被显著富集通路Hedgehog signaling pathway、JAK/STAT signaling pathway、TGF/BETA signaling pathway等参与了毛囊发育调控。综合KEGG、WGCNA、分子网络构建、GSEA分析结果,共同筛选得到基因CCL27、IL7、WNT2,推断这3个基因在调控毛囊发育中起到了关键的调控作用。本研究为进一步阐明动物毛发弯曲的分子机制提供了理论基础。  相似文献   

9.
旨在通过对北京油鸡(BY)和“广明2号”白羽肉鸡B系(广明白鸡,GM)在热应激和正常饲养条件下的脾脏转录组进行分析,鉴定两鸡种在应对热应激时的差异表达基因和信号通路,解析不同鸡种耐热性能差异的分子调控机制。试验动物采用相同饲养条件下的25日龄北京油鸡和广明白鸡,测定异嗜性粒细胞与淋巴细胞的比值(H/L)、超氧化物歧化酶(SOD)和总抗氧化能力(T-AOC)等表型数据,采集脾脏组织进行转录组测序(RNA-Seq)。基于RNA-Seq数据,结合表型进行加权基因共表达网络分析(WGCNA),筛选出与性状相关性较高的模块与基因。脾脏转录组测序结果表明,热应激组相比于对照组,北京油鸡中共鉴定到313个差异表达基因,其中169个上调表达,144个下调表达;广明白鸡中鉴定到235个差异表达基因,其中152个上调表达,83个下调表达。利用WGCNA分析,在北京油鸡中筛选出2个与H/L相关性较高的模块,在广明白鸡中筛选出4个与H/L和T-AOC强相关的模块。通过筛选模块中的Hub基因,发现两品种中均存在TRIM29基因,说明该基因可能在抗热应激方面发挥重要作用。本研究通过转录组分析,揭示了不同鸡种在热应...  相似文献   

10.
旨在通过转录组测序技术获得简州大耳羊肌内前体脂肪细胞成脂分化前后的差异表达基因,并经生物信息学分析获得相关信号通路及可能发挥作用的关键功能候选基因。本研究以7日龄的健康简州大耳羊公羊为试验动物(n=3),采用胶原酶消化法分离获得其肌内前体脂肪细胞;利用Illumina平台对肌内前体脂肪细胞和诱导分化5 d的肌内脂肪细胞cDNA样品(n=3)进行高通量测序;以|log2fold change|>0和P<0.05为阈值筛选获得差异表达基因,并利用clusterProfiler R包对其进行GO功能和KEGG通路富集;最后利用实时荧光定量PCR(quantitative real-time PCR,qRT-PCR)技术检测功能候选基因在细胞分化前后的表达量变化。结果显示,共获得差异表达基因7 916个,其中4 143个为表达上调基因,主要富集到氧化磷酸化、核糖体合成和三羧酸循环通路等305条通路;3 773个为表达下调基因,且主要富集到甲状腺激素信号通路等303条通路;差异基因GO功能注释中52.8%为生物过程、13.4%是细胞组成和33.8%是分子功能。qRT-PCR结果表明,UCP3、ACACBACOT11、ACOX3、APOA1和WISP2在分化前后的山羊肌内脂肪细胞中的表达趋势与RNA-seq结果一致,提示这些基因适合作为下一步研究的功能候选基因。本研究筛选得到简州大耳羊肌内脂肪细胞成脂分化的差异基因,并确认了6个基因可作为功能候选基因。研究结果为阐明肉用山羊肌内脂肪细胞成脂分化的分子调控网络提供基础数据和系统资料。  相似文献   

11.
The purpose of this study was to select the key genes of mammary gland development at different physiological stages in dairy goats by weighted gene co-expression network analysis(WGCNA). GSE14008 mammary gland tissue microarray data set of dairy goats at different physiological stages (pregnancy 46 days, 70 days, 90 days, 110 days and 40 days postpartum) was downloaded from GEO database, and co-expression analysis was carried out using WGCNA package of R language. The target modules were selected by correlation analysis between the modules and physiological stages, and the hub genes were selected according to the connectivity degree. After enrichment analysis of the modules using DAVID website, the protein-protein interaction network of the modules was constructed using String website and the core genes were obtained using Cytoscape software, and finally the target genes were obtained from intersection of core genes and hub genes. A total of 8 443 genes from 18 samples were analyzed for co-expression of weighted genes. The 30 modules were obtained, and 4 target modules related to different physiological stages and 30 hub genes of each module were selected. Meanwhile, protein-protein interaction network of 4 modules and 20 core genes of each network were also obtained. Finally, 13 target genes related to mammary gland development were obtained from the 4 modules, which were UQCR, RGL2, NOTCH1, PTBP1, PPP5C, FZR1, UBE2L3, TNF, MAT2A, ITGB2, GPR18, JAK1 and CSN2 genes. The key processes of mammary gland development at different physiological stages and the genes that played the key role in these processes were elucidated by WGCNA, GO enrichment analysis, PPI network and other bioinformatics techniques, which provided a new idea and clue for the further research on the mechanism of mammary gland development.  相似文献   

12.
旨在通过转录组测序技术筛选出共轭亚油酸(conjugated linoleic acids,CLA)影响边鸡胸肌脂类代谢的相关基因,为日粮中CLA调控边鸡脂类代谢的作用机制提供理论依据。本试验选取90日龄健康边鸡180只,随机分为5组,每组3个重复,每组日粮中分别添加不同比例的CLA 0%(对照组)、0.5%、1.0%、1.5%和2.0%,预饲期1周,正饲期6周。饲养试验结束后采集胸肌组织进行转录组测序,使用DESeq2软件对测序数据进行差异表达分析,对差异表达基因进行GO功能和KEGG通路富集分析,筛选出与胸肌脂类代谢相关的差异表达基因,利用qPCR对差异表达基因进行验证。通过转录组测序数据分析,共筛选出1 229个差异表达基因,其中594个基因上调,635个基因下调。随机选取9个差异表达基因进行qPCR验证,其相对表达量变化趋势与测序结果一致。GO功能分析发现,差异表达基因主要集中在生物学过程中的细胞过程、单一生物过程、生物调节和代谢过程。KEGG通路富集分析发现,各试验组中差异表达基因显著富集在不同的信号通路中,主要富集在ECM-受体相互作用、黏着斑、吞噬体、细胞凋亡、肌动蛋白细胞骨架调节和细胞因子-细胞因子受体相互作用等通路中。通过GO和KEGG功能注释分析,共筛选出MCAT、APOA1、PTGDS、ALDH3A2、PLTP、FABP3、FOXO1、UCP3和FABP4等18个主要的参与脂类代谢相关的差异表达基因,可能在调控边鸡胸肌脂类代谢中发挥重要的作用。本研究通过转录组测序筛选出CLA影响边鸡胸肌脂类代谢过程中的主要调控基因,发现了日粮中CLA影响差异表达基因作用的主要生物学过程和信号通路,为今后研究CLA影响边鸡胸肌脂类代谢的分子作用机制奠定了基础。  相似文献   

13.
【目的】 挖掘影响地方鸡体尺性状的有效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个目标性状候选基因, 为儋州鸡育种提供候选的分子标记, 为地方鸡标记辅助选择提供新的思路。  相似文献   

14.
为研究肌细胞生成素(myogenin,MyoG)和胰岛素样生长因子1(insulin-like growth factor 1,IGF-1)基因表达量对家禽生长发育的影响,试验以体型差异较大的盐津乌骨鸡和大围山微型鸡为材料,分别测定300日龄鸡的体重、体尺和屠宰性状来进行比较分析,检测胸肌、腿肌和肝脏中MyoGIGF-1基因mRNA的相对表达量,将生长、屠宰性状与基因的相对表达量进行两两相关性分析。结果显示,盐津乌骨鸡体重、体尺和屠宰性状等指标均高于大围山微型鸡,同时盐津乌骨鸡胸肌、腿肌和肝脏中MyoGIGF-1基因的表达量也高于大围山微型鸡,表达量趋势为:胸肌 > 腿肌 > 肝脏;相关性分析结果显示,两个品种鸡胸肌、腿肌和肝脏中MyoGIGF-1基因的表达量与体重、体尺和屠宰性状等指标呈显著正相关(P<0.05),其中MyoG基因的相对表达量与除胸深外的所有性状均呈极显著相关(P<0.01),而IGF-1基因的相对表达量与体重、体斜长、胸深、胸骨长均呈显著相关(P<0.05)。因此,MyoG、IGF-1基因可作为盐津乌骨鸡和大围山微型鸡的生长性状候选基因来进行选育。  相似文献   

15.
为深入了解家禽肉品质性状形成的分子机制,本实验测定金茅黑鸡胸肌的肉品质,利用RNA-seq技术进行转录组测序,分析转录组基因表达量与肉品质性状间的关联性。结果显示:与金茅黑鸡14周龄肉色显著相关的基因5个,pH57个,剪切力33个,系水力17个;PLEKHS1、MUC等基因为影响肉色的候选基因,LHX9、SPDEF、GRXCR1等基因为系水力的候选基因,NPY、POMC、FGF20、MGAT4C和GUK1基因为影响肌肉pH的候选基因,EDAR基因和HS6ST3基因为肌肉剪切力的主要候选基因;功能分析发现这些基因在调节肌肉纤维、肌间脂肪沉积、脂类代谢、糖代谢过程和磷酸二酯酶水解等方面发挥重要作用,进而影响胸肌肌肉品质的形成。  相似文献   

16.
试验采用PCR-RFLP方法检测脂滴包被蛋白(perilipin,PLIN)基因内含子6在济宁百日鸡、莱芜黑鸡等4个地方鸡种和1个培育品系中的遗传多态性,分析了多态位点不同基因型与鸡胴体及脂肪性状的相关性。结果发现,在5个供试群体中检测到1个多态位点,测序证实为新发现的鸡PLIN基因2 467 bp处G→A突变,该突变位点在供试群体中检测到3种基因型:A1A1、A1A2和A2A2,2个等位基因:A1和A2。等位基因A1在所有供试群体中均表现为优势等位基因。关联分析结果表明,PLIN基因2 467 bp位点对鸡部分胴体性状和脂肪性状影响显著(P<0.05)。多重比较结果表明,A1A1基因型个体活体重、屠体重、全净膛重、腹脂重和腹脂率均显著高于A2A2基因型个体(P<0.05)。A1A2基因型个体的胸肌肌内脂肪含量显著高于A1A1和A2A2基因型个体(P<0.05)。研究结果表明,PLIN基因对鸡胴体及脂肪性状有一定影响。  相似文献   

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

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