首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
为探究高低海拔斑头雁肺脏差异表达基因和差异代谢通路的变化,基于Illumina NovaSeq 6000测序平台对斑头雁肺脏组织进行转录组测序,分析差异表达基因,利用荧光定量PCR对部分差异表达基因进行验证,进一步对差异表达基因进行GO功能富集分析和KEGG信号通路富集分析。结果:与高海拔组斑头雁肺脏相比,筛选出差异表达基因84个,其中上调差异表达基因51个,下调差异表达基因33个。经GO功能注释后富集到21个显著性GO功能,"生物过程"显著富集到11个条目,涉及代谢过程以及内源性刺激的应答等。"细胞组成"显著富集到1个条目,涉及到细胞外区域。"分子功能"显著富集到9个条目,涉及到酶活性、激素、转录因子、钠离子跨膜转运等。注释到KEGG的显著性富集通路有3条,涉及到糖胺聚糖降解通路、丝裂原活化蛋白激酶信号通路以及丙氨酸、天冬氨酸和谷氨酸代谢通路。本研究结果丰富了斑头雁的基因资源,为该物种海拔适应相关基因的功能研究奠定了基础。  相似文献   

2.
旨在对美国短毛黑水貂快速生长过程中转录组差异表达进行分析。本试验将美国短毛黑水貂45、90日龄各3只健康公貂的胸肌组织作为试验材料,利用Illumina HiSeq~(TM)2500高通量测序平台建立转录组文库,对两个生长阶段差异显著性的基因进行GO功能富集和KEGG Pathway分析,寻找调控水貂快速生长期的相关候选基因和代谢通路,并利用qRT-PCR对转录组测序结果进行验证。结果表明,以P0.05,|log_2FC|1为条件筛选到279个显著差异基因,对这些差异基因进行GO和KEGG分析,筛选到与肌肉生长相关显著富集GO条目有66条、相关差异表达基因44个,其中上调20个,下调24个;与肌肉生长相关显著富集KEGG通路12条,如p53信号通路、PPAR信号通路、FOXO信号通路等。这些差异基因可能在水貂的快速生长过程中起到调控作用,可以作为后续研究水貂生长发育的候选基因。本研究结果为探究水貂生长发育调控机制及培育大体型水貂品种奠定基础。  相似文献   

3.
18头肥育猪按饲粮均分为2组,添加红花籽油组和不添加油脂组(对照组)。试验期60d结束,全部屠宰,分别采集肝脏组织样品,进行高通量转录组测序,找出2处理组间的差异表达基因和差异代谢通路。利用Illumina HiSeqTM2500高通量RNA-seq测序技术对2组肥育猪肝脏RNA测序,使用TopHat2软件将测序得到的reads序列与猪(Sscrofa10.2)参考基因组序列比对,找出差异表达基因,在Nr、GO和KEGG数据库中进行功能注释、富集分析和聚类分析。结果显示,红花籽油组和对照组肝脏中差异表达基因共有1 114个,与对照组相比,红花籽油组表达上调的基因有579个,下调的基因有535个。GO功能分类注释到细胞组成、生物学过程和分子功能数据库中的差异表达基因分别有902,903,857个。注释到KEGG通路中差异表达基因414个,显著富集通路4个(P0.05)。  相似文献   

4.
试验选取具有典型内脏型痛风特征的雏鹅5只和相同日龄及生长环境下的健康雏鹅6只,利用RNA-seq技术分析其肾组织转录表达差异,旨在了解雏鹅痛风发生的分子机制。结果表明:痛风雏鹅肾组织转录表达显著区别于健康雏鹅。试验共筛选出差异表达基因1 749个,上调基因501个,下调基因1 248个。GO(Gene Ontology)和KEGG(Kyoto Encyclopedia of Genes and Genomes)分析发现,这些差异基因主要具有趋化因子活性、T细胞受体结合、抗原结合等分子功能,参与免疫反应、趋化因子介导信号通路、炎症反应等生物过程,并主要富集于细胞因子受体相互作用、肠道IgA免疫网络、细胞黏附等信号通路。  相似文献   

5.
旨在分析初生和成年牦牛大脑皮质转录组表达谱,筛选影响大脑皮质发育的候选基因。本研究选取初生(1~7日龄(NYB),n=3)和成年(3~4岁(AYB),n=3)牦牛,放血处死,采集大脑皮质,利用Illumina HiSeq平台对转录组进行测序与分析,筛选得到差异表达基因,在此基础上对差异基因进行GO功能注释、KEGG富集分析和实时荧光定量PCR(qRT-PCR)方法验证。测序结果共分析得到了4 790个显著差异表达的基因(P<0.05),其中包括2 670个上调表达的差异基因,2 120个下调表达的差异基因。随机选取了9个差异基因进行qRT-PCR验证,表达趋势与转录组测序结果一致,表明测序结果可靠。GO功能注释发现其中占比较多的条目为与蛋白合成及ATP利用相关的转录调控、RNA聚合酶II正调控转录、蛋白结合及ATP结合;KEGG富集分析发现占比较多与大脑发育相关的通路为调控轴突发育及重塑的轴突导向通路,与神经胶质细胞增殖相关的MAPK信号通路,与突触传递相关的神经活性配体-受体相互作用、 PI3K-Akt、钙信号及磷脂酶D信号通路,与免疫防御相关的Toll样受体通路,与低氧适应调...  相似文献   

6.
旨在基于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等可能在塔里木马鹿毛色形成过程中发挥重要作用。本研究结果为今后塔里木马鹿毛色相关基因的分子调控机制方面及挖掘潜在的新基因提供了丰富的试验数据。  相似文献   

7.
痛风雏鹅肾组织损伤的RNA-seq分析   总被引:2,自引:2,他引:0  
试验选取具有典型内脏型痛风特征的雏鹅5只和相同日龄及生长环境下的健康雏鹅6只,利用RNA-seq技术分析其肾组织转录表达差异,旨在了解雏鹅痛风发生的分子机制。结果表明:痛风雏鹅肾组织转录表达显著区别于健康雏鹅。试验共筛选出差异表达基因1 749个,上调基因501个,下调基因1 248个。GO(Gene Ontology)和KEGG(Kyoto Encyclopedia of Genes and Genomes)分析发现,这些差异基因主要具有趋化因子活性、T细胞受体结合、抗原结合等分子功能,参与免疫反应、趋化因子介导信号通路、炎症反应等生物过程,并主要富集于细胞因子受体相互作用、肠道IgA免疫网络、细胞黏附等信号通路。  相似文献   

8.
《中国兽医学报》2019,(12):2458-2466
旨在通过对奶牛干奶期和泌乳初期肝脏组织进行转录组测序和生物信息学分析,探明其不同生理阶段的差异表达基因和差异代谢通路。选取305 d产奶量接近、干奶天数一致、怀孕天数接近的中国荷斯坦奶牛3头,分别在干奶期(产前50 d)和泌乳初期(产后15 d)活体采集肝脏组织,利用Illumina Hiseq~(TM)2500高通量RNA-seq测序技术对干奶期和泌乳初期肝脏RNA测序,使用TopHat2软件将测序得到的reads序列与牛(UMD3.1.66)参考基因组序列比对,在GO和KEGG数据库中进行聚类分析和富集分析。结果显示,泌乳初期和干奶期肝脏中差异表达基因共有409个,与干奶期相比,泌乳初期表达上调的基因有235个,下调的基因有174个。GO功能分类注释到生物学过程、细胞组成和分子功能数据库中的GO条目分别有63,41,15条。注释到KEGG的显著富集通路有21条(P0.05)。该研究结果为挖掘奶牛产奶性状的候选基因提供了技术依据。  相似文献   

9.
旨在通过转录组测序技术获得简州大耳羊肌内前体脂肪细胞成脂分化前后的差异表达基因,并经生物信息学分析获得相关信号通路及可能发挥作用的关键功能候选基因。本研究以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个基因可作为功能候选基因。研究结果为阐明肉用山羊肌内脂肪细胞成脂分化的分子调控网络提供基础数据和系统资料。  相似文献   

10.
旨在通过转录组测序技术获得简州大耳羊肌内前体脂肪细胞成脂分化前后的差异表达基因,并经生物信息学分析获得相关信号通路及可能发挥作用的关键功能候选基因。本研究以7日龄的健康简州大耳羊公羊为试验动物(n=3),采用胶原酶消化法分离获得其肌内前体脂肪细胞;利用Illumina平台对肌内前体脂肪细胞和诱导分化5 d的肌内脂肪细胞cDNA样品(n=3)进行高通量测序;以|log_2fold change|0和P0.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、ACACB、ACOT11、ACOX3、APOA1和WISP2在分化前后的山羊肌内脂肪细胞中的表达趋势与RNA-seq结果一致,提示这些基因适合作为下一步研究的功能候选基因。本研究筛选得到简州大耳羊肌内脂肪细胞成脂分化的差异基因,并确认了6个基因可作为功能候选基因。研究结果为阐明肉用山羊肌内脂肪细胞成脂分化的分子调控网络提供基础数据和系统资料。  相似文献   

11.
试验设计了椰子油组和对照组2个处理。采集育成猪肝脏组织,进行转录组高通量测序,找出2种处理间的差异表达基因及差异代谢通路。利用Illumina Hi SeqTM2 500高通量RNA-seq测序技术测序,使用Top Hat2软件将测序得到的reads序列与猪参考基因组(Sscrofa10.2)序列比对,找出差异表达基因,并在Nr、GO和KEGG数据库中进行功能注释、富集分析和聚类分析。结果表明:椰子油组和对照组肝脏中差异表达基因共有487个,与对照组相比,椰子油组表达上调的基因有222个,下调的基因有265个。在差异表达上/下调前15位基因中,有多个基因是参考基因组中没有的新基因和功能未知基因。GO功能分类注释到细胞组成、分子功能和生物学过程数据库中的差异表达基因数分别有398、368个和398个;注释到KEGG通路中差异表达基因数188个。  相似文献   

12.
捻转血矛线虫阿苯达唑耐药株给药前后比较转录组学分析   总被引:1,自引:0,他引:1  
为探明捻转血矛线虫阿苯达唑耐药株给药前后在转录组水平上的差异,本研究采用Illumina Hiseq4000对给药前后的捻转血矛线虫进行转录组测序,筛选得到差异表达基因并通过GO和KEGG数据库对其进行功能注释及富集性分析,并利用荧光定量PCR验证部分差异表达基因。结果显示,共筛选获得851个显著差异表达基因,其中包括584个上调基因和267个下调基因。通过GO功能富集分析显示,有458、418和367个基因分别注释到生物学过程、分子功能和细胞组分三大类;对差异表达基因进行KEGG富集分析显示,有173个差异表达基因参与到75条KEGG通路中,显著富集在核糖体合成、细胞凋亡、过氧化物酶增殖体激活受体(PPAR)等信号通路。本试验初步筛选了阿苯达唑耐药株给药前后的差异表达基因,为深入探索捻转血矛线虫耐药性分子机制、筛选对耐药性检测的分子标记及耐药性早期鉴别诊断方法的建立提供重要数据。  相似文献   

13.
牦牛发情期卵巢比较转录组学研究   总被引:1,自引:0,他引:1  
旨在进一步了解牦牛发情期卵巢的分子机制,解析牦牛繁殖的特殊性。本研究应用RNA-seq技术对牦牛和平原黄牛发情期卵巢进行转录组高通量测序及全基因组差异表达模式比对分析。通过比较分析牦牛和黄牛卵巢转录组数据,共筛选出1 307个差异表达基因,其中661个基因表达量上调和646个基因表达量下调。进一步功能分析表明,这些差异基因涉及多种GO分类及KEGG通路。其中GO分类注释显示,差异基因与细胞粘附、激素调控等生物学过程存在密切关联,同时钙离子结合、阳离子跨膜转运等分子事件表现活跃。KEGG通路分析显示,补体和凝血级联通路的富集水平最高,其次为细胞色素P450相关通路。昼夜节律等一些新型通路,尽管与生殖功能没有明显关联,但也表现出显著富集。本研究首次对比牦牛和黄牛发情期卵巢转录组数据,筛选并分析相关差异基因。该研究结果为进一步阐述牦牛卵巢的基本分子机理提供基础,同时也为全面理解牦牛繁殖特异性提供新的思路。  相似文献   

14.
sRNA伴侣蛋白Hfq敲除条件下沙门菌的转录组分析   总被引:1,自引:0,他引:1  
为探明沙门菌(Salmonella)在sRNA伴侣蛋白Hfq敲除条件下转录组变化情况。本试验采用HiSeq测序平台对沙门菌LT2菌株及其Hfq敲除株进行高通量测序,通过DESeq2差异分析方法筛选敲除sRNA伴侣蛋白Hfq条件下的差异表达基因,对其进行生物信息学GO功能显著性富集分析和Pathway显著性富集分析。结果显示,对照组和试验组分别得到12 753 534条、8 254 308条Clean Reads。通过DESeq2差异分析获得差异基因1 055个,其中表达上调基因516个,表达下调基因539个。GO功能显著性富集分析表明显著差异表达基因主要富集在钴胺素代谢过程、钴胺素生物合成过程、碳水化合物代谢过程等生物过程中。Pathway显著性富集分析表明显著差异表达基因富集到细菌趋化性、丙酸代谢和碳代谢等11个信号通路中。结果表明,本试验应用RNA-seq技术丰富了sRNA伴侣蛋白Hfq调控基因数量,筛选出20个受伴侣蛋白Hfq调控的显著差异表达基因,其中4个基因功能及编码蛋白未知,注释了差异表达基因的生物学功能及信号通路,推断Hfq在细菌对数期主要通过影响营养提供和趋化性等途径,继而影响细菌的正常生理活动,为Hfq协同sRNA的调控机理研究及sRNA靶基因的筛选奠定了基础。  相似文献   

15.
为初步鉴定并挖掘出猪戊型肝炎病毒(Hepatitis E virus,HEV)ORF3蛋白影响HepG2细胞核黄素代谢信号通路的lncRNA-mRNA调控网络,本试验通过构建腺病毒过表达载体,制备高滴度过表达腺病毒,介导猪 HEV-ORF3在HepG2细胞中实现过表达。Western blotting检测ORF3蛋白过表达成功后,运用lncRNA高通量组学测序,筛选出差异表达的lncRNA并进行靶向差异基因预测,对lncRNA靶向差异基因进行GO功能和KEGG通路富集分析,初步鉴定出与核黄素代谢信号通路相关的lncRNA-mRNA调控网络。Western blotting结果显示,在约为12 ku处出现目的条带,说明成功实现腺病毒介导猪HEV-ORF3在HepG2细胞中过表达。lncRNA高通量组学测序结果显示,共发现102个显著差异表达lncRNAs表达量上调,80个显著差异表达lncRNAs表达量下调。GO功能和KEGG通路富集分析显示,初步挖掘出lncRNA(MSTRG.13995.2)和lncRNA(MSTRG.1960.1)可能是与核黄素代谢信号通路相关的显著差异表达lncRNA,分别通过顺式调控其靶向基因APC-5和FLAD1来影响核黄素代谢信号通路。本试验初步鉴定出lncRNA(MSTRG.13995.2)-APC5和lncRNA(MSTRG.1960.1)-FLAD1可能是影响HepG2细胞核黄素代谢信号通路的lncRNA-mRNA调控网络。  相似文献   

16.
为鉴定禽源多杀性巴氏杆菌(avian Pasteurella multocida,Pm)四环素(Tetracycline,Tet)敏感株与耐药株(MIC=32μg/mL)在转录组水平上的差异,本研究采用RNA-seq技术对Tet敏感株与耐药株进行双向测序,通过GO、KEGG和ARDB数据库对表达差异基因进行注释和富集性分析。结果显示,共筛选出442个差异表达基因,其中包括103个上调基因和339个下调基因。通过GO功能富集分析显示差异表达基因多数富集于核糖体结构成分、rRNA结合、跨膜运输和离子转运等生物过程;对差异表达基因进行信号通路富集性分析,大多数差异表达基因分布在ABC转运蛋白与核糖体代谢途径中;对差异表达基因进行耐药基因注释(ARDB),发现其中的耐药基因主要为核糖体保护蛋白和外排系统两类。本研究表明禽源Pm对四环素耐药的耐药机制可能与细胞膜上ABC转运蛋白的过量表达以及核糖体保护蛋白有关。  相似文献   

17.
旨在对蒙古马高负荷运动训练前后转录组差异表达进行分析。本试验对6匹进行了4个月高负荷运动训练的蒙古马进行研究,分别在训练前后两个时期采集肌肉样品,利用二代测序技术对肌肉样品建立转录组文库,对发生显著变化的差异基因进行GO功能富集和KEGG Pathway分析,寻找训练前后发生显著变化的代谢通路及相关候选基因。结果表明,在转录组筛选到的1 102个显著差异表达基因中,有299个在训练后发生下调,803个发生上调。这些差异基因被注释到3 398条GO term上。其中与运动学相关的生物学过程有:肌肉结构发育、心血管系统发育、循环系统发育、肌肉组织发育、肌肉器官发育、肌肉收缩和肌肉细胞分化等。对差异表达基因进行KEGG pathway分析,结果表明,差异基因被富集到如扩张型心肌病、肥厚型心肌病、心肌收缩、钙信号通路、肌动蛋白骨架调节等与运动相关的193条通路中。本研究结果为探究蒙古马强耐力的基本分子机理提供理论基础,同时也为理解蒙古马运动学特性提供新的思路。  相似文献   

18.
本文旨在通过转录组测序和生物学信息分析探究在鸡的成骨发育过程中可能发挥重要作用的调控基因和生物学通路。从18胚龄明星肉鸡的胫骨获取骨髓间充质干细胞(BMSC),并进行体外培养与染色鉴定。试验分为试验组和对照组,每组3个生物学重复,试验组将BMSC进行体外成骨诱导分化21 d,对照组不进行处理。分别提取试验组和对照组细胞样品的总RNA,进行转录组分析。结果表明:每个样品平均获得了23 M clean reads数,大约20 M clean reads被比对到鸡的参考基因组上。相比于对照组细胞,试验组筛选到2 715个差异表达基因,包含1 437个上调的差异基因,1 278个下调的差异基因。GO富集分析表明,差异基因主要富集在骨骼系统发育、骨矿化、骨化以及整合素结合等相关功能;KEGG通路主要富集在细胞的焦点黏着、物质合成及能量代谢以及细胞外基质受体相互作用等相关通路。在20条骨发育相关的GO条目或KEGG通路中总共富集到182个基因,其中上调的基因81个,下调的基因101个。通过查找NCBI、Web of Science等相关数据库以及文献,对BMSC成骨分化的过程中发挥重要作用的基因进...  相似文献   

19.
本研究旨在利用RNA-seq技术对不同季节槟榔江水牛全血mRNA进行分析,找出不同季节精液品质候选基因。选取5头生长状况相近的3~5岁槟榔江水牛,进行春、夏、秋、冬4个季节的精液品质(鲜精活力、精子畸形率、射精量、精子密度)检测。根据精液品质,选取差异极显著的春(BL)、夏(SBL)两个季节的水牛样品进行RNA-seq转录组测序分析。颈静脉采集全血,抽提mRNA并组内混合制样,利用RNA-seq技术分析不同季节槟榔江水牛全血mRNA差异表达。根据差异倍数1.5,P0.05的标准进行筛选,得到差异基因,再利用DAVID6.7软件对差异基因进行GO和KEGG分析,通过RT-PCR验证基因芯片结果的准确性。结果表明:差异表达基因84个,其中上调基因56个,下调基因28个;GO分析发现,这些差异基因与应激反应、免疫反应有关;KEGG分析表明,差异基因显著富集到肿瘤坏死因子(TNF)和T细胞受体信号通路。其中,神经营养因子受体类似物死亡结构域(NRADD,又叫热休克蛋白因子1(HSF1))、原癌基因FOS和JUN,与精液品质相关性强,可作为进一步研究的靶基因。  相似文献   

20.
试验旨在探讨过氧化物酶体增殖物激活受体γ(peroxisome proliferators-activated receptors,PPARγ)激动因子(环格列酮)和胰岛素对延边黄牛前脂肪细胞脂质代谢的影响及可能的调控机制。试验分为对照组(CON,5% FBS)、胰岛素处理组(I组,5% FBS+10 mg/L胰岛素)、环格列酮处理组(C组,5% FBS+10 mg/L环格列酮)、胰岛素+环格列酮处理组(IC组,5% FBS+10 mg/L胰岛素+10 mg/L环格列酮)。各组处理144 h,通过RNA-Seq技术对其进行转录本测序,差异表达基因进行GO功能注释和KEGG富集分析。差异表达基因分析结果显示,与对照组相比,IC组差异表达基因数量为1 764个,比I和C组分别多276和569个,3个试验组间共有371个差异基因;IC组上调基因数量为974个,分别比I和C组多140和249个,IC组下调基因为790个,I和C组下调基因分别为654和470个。差异基因的GO分析表明,各处理组的差异基因参与了细胞过程、代谢过程以及生物过程的调节;在分子功能上,主要是分子功能调节剂、转运活性、转录调节活性等;在细胞成分上,大多数差异基因被富集到细胞器、细胞膜及胞外区等。KEGG通路分析发现,差异表达基因主要富集在FoxO、PPAR、TGF-β、p53等信号通路。综上所述,胰岛素和环格列酮单独处理与二者共同处理对延边黄牛前脂肪细胞分化的调控机制有所差别,其中二者共同处理对分化的促进作用更加明显。本研究结果为研究环格列酮与胰岛素在调节脂肪生成中的功能和分子机制提供了依据。  相似文献   

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

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