首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
目的:筛选比较大鼠C6胶质瘤细胞与正常星型胶质细胞中的差异表达基因。方法:利用Illumina Hi Seq4000技术对大鼠C6胶质瘤细胞和正常星形胶质细胞进行转录组测序并对测序结果进行功能注释。以|log2 Fold Change|≥1和q值0.05为标准筛选差异表达基因,并通过GO和KEGG数据库分析其功能和参与的信号通路。结果:较之正常星形胶质细胞,大鼠C6胶质瘤细胞中共筛选得到4363个差异表达基因,其中1603个基因表达上调,2760个基因表达下调;GO富集结果表明:4109个差异表达基因共映射到906个GO term,主要与GTP酶活性正性调节、基因表达正性调节、蛋白结合和钙离子结合等功能相关;KEGG富集结果表明:1535个差异表达基因共参与106条通路,富集最多的是PI3K-Akt信号通路和癌症发生通路。结论:成功获取大鼠C6胶质瘤细胞与正常胶质细胞的基因表达谱,并对其差异基因进行GO和KEGG注释,结果有助于胶质瘤发生发展的研究。  相似文献   

2.
目的应用转录组测序方法, 分析北京地区人呼吸道合胞病毒(respiratory syncytial virus, RSV)A亚型优势流行基因型ON1地方株感染A549细胞后差异表达基因, 为RSV防治提供潜在靶点。方法选用已经过全基因组测序确定为RSV A亚型ON1基因型的地方株(61397-ON1)感染A549细胞, 提取总mRNA, 通过转录组测序筛选出与未感染的A549细胞为对照的差异表达基因, 对其进行GO分析、KEGG通路分析, 同时随机选择6个差异表达倍数大于2倍的基因进行qRT-PCR验证。结果以未感染的A549细胞为对照, 筛选出1 632个差异表达基因, 其中807个基因表达上调, 825个基因表达下调。差异基因主要参与细胞因子反应以及MAPK级联反应正向调控等免疫应答相关生物过程, 并在MAPK信号通路、NOD样受体信号通路、p53信号通路、TNF信号通路、IL-17信号通路及NF-κB信号通路发生了富集。选择的6个差异表达基因qRT-PCR验证结果与转录组数据趋势一致。结论 RSV A亚型ON1基因型毒株感染A549细胞后的差异表达基因主要参与细胞因子应答及免疫相...  相似文献   

3.
目的:分析富含AT序列的特异性核基质区结合蛋白1(SATB1)表达改变对人鼻咽癌(NPC)细胞中蛋白表达谱特征的影响,并采用生物信息学富集分析差异信号通路。方法:利用SATB1过表达慢病毒和阴性对照慢病毒感染CNE1细胞,获得稳定表达细胞系。提取2组细胞的总蛋白,利用TMT标记蛋白定量技术和串联质谱分析技术筛选出差异表达蛋白,采用RT-qPCR对差异蛋白编码基因在mRNA水平进行定量验证;应用GO数据库对差异表达蛋白进行注释和富集分析,应用KEGG数据库对差异蛋白涉及信号通路进行富集分析。结果:SATB1过表达慢病毒感染CNE1细胞后获得稳定表达细胞系。与对照组相比,过表达SATB1的CNE1细胞鉴别出278个差异表达蛋白,其中115个表达上调,163个表达下调;采用RT-qPCR对10个代表性差异蛋白的编码基因进行mRNA水平的验证,其结果与蛋白组学结果具有一致性。GO数据库分析差异表达蛋白主要参与了细胞过程、单有机体过程、生物调控和代谢过程,发挥生物分子结合和催化等分子功能;细胞组件主要存在于细胞部分、细胞以及细胞器上。KEGG数据库分析显示差异表达蛋白参与和肿瘤关系密切的信号通路...  相似文献   

4.
目的:分析胚鼠腭裂形成过程中腭胚组织全基因组DNA差异甲基化位点相关基因及基因组差异甲基化水平。方法:使用高效低成本全基因组DNA甲基化检测技术——甲基化修饰依赖性内切酶测序法(MethylRAD-Seq法)对妊娠第13.5天(GD13.5)、GD14.5和GD16.5(n=18)的胚鼠腭裂模型组和对照组腭胚组织进行全基因组DNA甲基化测序,比较基因组不同功能元件(3’端非翻译区、5’端非翻译区、TSS2000、内含子、外显子和基因间区)中甲基化位点的差异分布及甲基化水平差异基因,并对相关基因进行GO和KEGG通路分析。结果:(1)DNA不同元件的甲基化位点的峰值个数和峰值,对照组无明显变化,实验组甲基化整体水平随着时间推移逐渐降低;(2)GO功能富集和KEGG通路富集分析差异甲基化位点的基因及甲基化水平有差异的相关基因,其中与腭裂形成相关的基因为Fgf16、Jarid2、Kdm6a、Nr3c1、Tbx22和Trp53;(3)与腭裂形成相关的信号通路主要有MAPK、Wnt和TGF-β信号通路。结论:DNA甲基化在腭裂形成过程中可能起重要的调控作用。  相似文献   

5.
目的:筛选在结直肠癌远处转移中T细胞淋巴瘤侵袭转移1(T-cell lymphoma invasion and metastasis 1, Tiam1)相关基因并鉴定其表达。方法腹腔注射二甲基肼( dimethylhydrazine, DMH)建立小鼠结直肠癌模型,收集Tiam1转基因小鼠和野生型小鼠的结直肠癌原发灶新鲜组织5对,应用Affymetrix小鼠全基因组表达谱芯片检测两组的差异表达基因,qRT-PCR技术验证部分差异基因以确定结果的可靠性。结果基因芯片检测结果显示差异在3倍以上的基因共794个,有远处转移的Ti-am1转基因小鼠结直肠癌组织中表达上调的基因共400个,下调基因共394个;通过生物信息学分析及GENEMANIA共表达网络构建,筛选3个(PIK3R5、FIGF、RPS6KA6)特异相关的差异基因进行qRT-PCR技术验证,验证结果与基因芯片结果检测相符。结论通过基因芯片技术筛选和鉴定结直肠癌转移中Tiam1相关基因,建立了特异性结直肠癌转移表达谱,为寻找结直肠癌远处转移的分子标志物提供依据。  相似文献   

6.
目的:筛选1型糖尿病(T1DM)的潜在分子标志物并分析HCST参与T1DM的可能机制。方法:从GEO数据库下载T1DM和健康人群的单个核细胞微阵列表达数据,借助GEO2R分析差异基因,同时使用R语言进行火山图和热图的可视化。利用String在线工具对差异基因进行PPI分析,MCODE插件筛选关键子网络,DAVID在线分析工具对子网络基因进行GO和KEGG富集分析;CytoHubba插件对PPI网络进行关键基因筛选。结果:T1DM组与对照组相比共得到71个差异基因,其中包含上调基因22个,下调基因49个;利用MCODE插件从PPI网络中筛选1个关键子网络,得分7.5,共有9个节点和30个互作。子网络的9个基因的GO和KEGG富集分析发现,这9个基因主要参与调节免疫反应、细胞的防御反应、免疫应答及细胞溶解的生物过程,且涉及的主要信号通路是自然杀伤细胞介导的细胞毒作用的信号通路;四种算法得到的关键基因通过Venn分析最终筛选出1个关键基因——HCST。HCST在T1DM组明显低于健康对照组,差异有统计学意义。HCST的ROC曲线分析结果显示曲线下面积为0.983 3,P=0.000 1,差异有统计学意义。结论:HCST可成为诊断T1DM的潜在分子标志物和治疗T1DM的潜在靶点。  相似文献   

7.
目的:筛选和分析大鼠胶质瘤模型与正常星型胶质细胞中的差异表达基因。方法:构建大鼠C6颅内胶质瘤模型,利用Illumina Hi Seq 4000技术对大鼠胶质瘤模型的胶质细胞和正常星形胶质细胞进行转录组测序并对测序结果进行功能注释。以q值0. 05为标准筛选差异表达基因,并通过GO和KEGG数据库分析其功能和参与的信号通路。结果:较之正常星形胶质细胞,大鼠胶质瘤模型组细胞共筛选得到9221个差异表达基因,其中4604个基因表达上调,4617个基因表达下调; GO富集结果表明:这些差异表达基因共映射到575个GO term,与各类结合相关的分子如蛋白结合、离子结合、小分子结合、核苷酸结合等功能富集显著; KEGG富集结果表明:差异表达基因共参与22条通路,富集最多的是癌症发生通路、MAPK通路、胞吞通路和黏着斑通路。结论:成功获取大鼠胶质瘤模型细胞与正常胶质细胞的基因表达谱,差异表达基因富集在癌症发生通路、MAPK通路、胞吞作用通路和黏着斑等生物学通路,这些信号通路可能在胶质瘤发生发展过程中发挥重要作用,对其进一步分析将为临床治疗胶质瘤提供新的靶点。  相似文献   

8.
目的 利用生物信息学方法分析溃疡性结肠炎(UC)的枢纽基因和关键通路。方法 通过基因表达数据库(GEO)下载UC表达谱芯片GSE134025。利用R语言的Limma包筛选UC组与正常组肠黏膜细胞差异表达基因,对这些基因进行GO和KEGG分析;利用STRING数据库构建蛋白互作网络(PPI)并将结果导入Cytoscape筛选出核心基因;利用KEGG mapper绘制核心基因所在的信号通路图。结果 筛选出190个差异表达基因,其中147个上调,43个下调。GO分析结果表明,差异表达基因主要涉及炎症反应、细胞增殖的正调控等生物学过程,主要富集于细胞膜、质膜等细胞成分,具有肝素结合、生长因子等分子功能。KEGG分析显示差异基因主要富集于趋化因子信号通路,细胞因子受体相互作用等相关信号通路。从PPI网络中筛选出10个枢纽基因:白细胞介素6(IL-6)、CXC趋化因子配体8(CXCL8)、CXCL10、CXCL1、CXCL9、膜联素A1(ANXA1)、IL-1β、CC趋化因子配体20(CCL20)、CXCL2、CXCL11,其中多个基因位于IL-17调控的信号通路下游。结论 发现了10个与UC发病...  相似文献   

9.
目的 探究多囊卵巢综合征(PCOS)骨骼肌氧化应激差异基因,识别PCOS患者骨骼肌关键的氧化应激基因,并预测基因及转录因子调控网络。方法 以PCOS患者的骨骼肌为研究对象,筛选GEO数据库数据集GSE8157和GSE6798,通过R软件获得差异表达的氧化应激相关基因。然后对差异表达的基因进行基因(GO)富集分析、KEGG信号通路富集分析及转录因子的预测。结果 在26例PCOS患者和26例对照女性,共鉴定出7个氧化应激差异基因(PDX1、SLC1A1、SMPD3、PAWR、CHEK2、WRN、S100A8)及4个转录因子(NFKB1、STAT3、YY1、FOXC1)。GO和KEGG富集分析表明,有多个与氧化应激及其相关的富集术语,在IL-17信号通路,p53信号通路等KEGG通路富集。结论 我们鉴定出的7个PCOS骨骼肌潜在的氧化应激相关基因和4个转录因子,可能通过调节骨骼肌氧化应激在PCOS的发生发展中发挥重要作用。  相似文献   

10.
目的:对再生障碍性贫血(再障)患者进行外周血全外显子测序,观察基因突变的类型及突变率以探讨再障的发病机制。方法:采集19例再障患者的外周血并提取单个核细胞进行全外显子测序,对所得基因进行GO和KEGG富集分析。数据采用SPSS 22.0进行统计分析,以P<0.05为存在显著差异。结果:测得19例再障患者共有2 006个基因突变,参照GeneCards及OMIM数据库及患者测序结果筛选出49个基因,绘制患者基因突变图谱,发现MUC5B基因可能与再障发病相关,其突变率达100%,其余基因如ERCC6、SYNE1和WFS1等突变率均超过40%,亦可能参与再障的形成。在GO分析中其主要与细胞运动、细胞骨架及组成、蛋白结合等功能相关。KEGG通路分析中,主要是在Hippo通路、ECM-受体相互作用通路等地方富集。结论:MUC5B基因可能是造成T细胞比例失调,进而导致再障的关键基因。此外,Hippo通路可能通过调节T细胞比例、端粒酶活性等参与再障的发病机制。  相似文献   

11.
目的通过生物信息学方法探究类风湿关节炎患者的滑膜成纤维细胞差异表达基因及相关信号通路,寻找潜在的类风湿关节炎特异性分子标志物。方法利用R语言limma包等程序方法分析基因芯片GSE21959并筛选差异基因(differentially expressed genes,DEGs),利用DAVID数据库分析DEGs获得其GO富集分析和KEGG信号通路分析的结果。利用STRING数据库构建蛋白互作网络,再将结果导入Cytoscape软件中模块化核心基因并绘制蛋白互作网络图。结果筛选获得了123个差异基因,其中表达上调的基因38个,表达下调的基因85个。GO富集分析表明DEGs主要参与了趋化因子调节、CXCR趋化因子受体结合和血管生成正向调控等生物学过程,KEGG信号通路富集分析主要包括了趋化因子信号通路、Rap1信号通路和血管平滑肌收缩等信号通路。模块化分析获得了7个核心基因分别为:CXCL1、CXCL8、CXCL6、ADRA2A、ADCY8、S1PR1和SAA1。结论通过生物信息学分析获得类风湿关节炎的DEGs、核心基因、生物学过程和信号通路等信息,为探究类风湿关节炎的发病机制、发现诊断标志物和探索新治疗靶点提供理论依据与新的方向。  相似文献   

12.
目的:基于癌症基因组图谱(The Cancer Genome Atlas,TCGA)数据集,探索tspear基因在结直肠癌中的预后作用及机制.方法:自TCGA获取数据,应用R语言分析结直肠癌组织和正常组织间tspear mRNA的表达差异及TSPEAR表达与临床病理特征之间的相关性.采用Kaplan-Meier分析、比例风险回归分析(Cox''s proportional hazards regression model,Cox)和列线图确定tspear对结直肠癌患者总生存率和疾病特异性生存率的预测价值.对tspear 相关的差异表达基因进行基因集富集分析(Gene Set Enrichment Analysis,GSEA)、基因本体富集(Gene Ontology,GO)和信号通路富集(Kyoto Encyclopedia of Genes and Genomes,KEGG).结果:tspear在结直肠癌中的表达显著增加,其表达升高与肿瘤T分期、M分期、残余肿瘤、总生存期和疾病特异性生存期显著相关.单变量Cox回归分析表明tspear是结直肠癌生存率的预测因子.多变量Cox回归分析显示tspear是结直肠癌的独立预后因素.KEGG分析显示tspear相关差异表达基因主要富集于神经活性配体-受体相互作用通路.结论:tspear可能是一种新的结直肠癌独立预后生物标志物.  相似文献   

13.
目的:利用全基因组表达谱芯片筛查与卵巢浆液性囊腺癌发生相关的基因,对在卵巢浆液性囊腺癌发生过程中可能参与的基因间的信号转导通路进行分析。方法:选取癌症基因组图谱(TCGA)数据库中卵巢浆液性囊腺癌的Affymetrix Gene Chip Human Exon 1.0 ST Array数据共16张,分别为卵巢浆液性囊腺癌组8张和正常组8张,筛选出差异表达基因,并进行基因本体(gene ontology,GO)分析和信号通路分析,构建卵巢浆液性囊腺癌相关基因间的信号转导通路,分析网络中具有重要作用的基因。结果:共筛选出1 144个在卵巢癌中差异表达的基因,其中表达上调的基因有747个,表达下调的基因有397个。GO分析得到上调差异基因的显著性功能分析结果362项,下调差异基因的显著性功能分析结果 160项(P0.05)。其中包括与肿瘤发生相关的基因功能有细胞周期、DNA复制、细胞增殖、细胞凋亡、细胞黏附等。信号通路分析得到45个显著上调信号通路和14个显著下调信号通路(P0.05)。其中参与肿瘤发生相关的信号通路主要有细胞周期、P53信号通路、DNA复制、肿瘤中的信号通路、PI3K-Akt信号通路、ECM-receptor信号通路、细胞黏附因子、细胞凋亡等。挑选显著性基因功能和信号通路分析的交集基因229个,构建显著性GO与信号通路基因间信号转导网络。分析发现CDK1、PLK1、MCM3和PGK1这4个基因在卵巢癌的基因调控网络中具有重要作用。结论:卵巢浆液性囊腺癌中有大量差异表达基因,差异表达的基因在多个与肿瘤发生密切相关的信号通路中发挥重要的调控作用。  相似文献   

14.
目的 运用生物信息学方法筛选子宫内膜异位症的差异表达基因,系统探讨子宫内膜异位症的病因病理机制。方法 通过GEO数据库搜索子宫内膜异位症基因芯片,联合Perl语言及R语言分析异位子宫内膜组织与正常组织的差异表达基因;利用STRING数据库及Cytoscape软件绘制PPI网络图;通过g:profiler数据库对关键差异基因进行GO和KEGG富集分析。结果 共纳入2套子宫内膜异位症基因芯片,筛选出共同差异表达基因308个(178个上调,130个下调),剔除degree值为1和2的基因后得到关键差异表达基因170个(上调91个,下调79个),其中差异最显著的基因为:TAGLN(上调)、EpCAM(下调);PPI网络图中得到degree值最大的差异基因为:CDK1。差异表达基因GO生物过程富集分析得出细胞周期改变、细胞黏附、增殖异常等与EMs发生发展密切相关,KEGG通路富集分析发现包括ECM受体相互作用通路、局灶性黏着斑、PI3K-Akt信号通路等在内的8条通路参与其中。结论 子宫内膜异位症的病因病理机制可能与TAGLN、EpCAM、CDK1基因表达异常及其涉及的表观遗传修饰改变、上皮-间...  相似文献   

15.
目的:探讨结直肠癌细胞中脾酪氨酸激酶基因甲基化和表达的关系。方法:应用亚硫酸盐修饰测序、甲基化特异性聚合酶链反应和蛋白印迹技术检测结直肠癌细胞脾酪氨酸激酶的甲基化状态以及表达情况;荧光素酶报告分析法研究启动子区域CpG岛的甲基化与启动子活性的关系;甲基化转移酶抑制剂处理脾酪氨酸激酶甲基化失表达的结直肠癌细胞株,观察处理前后细胞内脾酪氨酸激酶基因甲基化状态和表达情况。结果:(1)23个结直肠癌细胞中,9个细胞启动子发生甲基化而失去蛋白质表达;其余则正常表达,甲基化发生率为39.2%;(2)9个甲基化的细胞中,7个存在微卫星不稳定;而14个未发生甲基化的细胞中,仅有4个存在微卫星不稳定。二者之间的差异显著(P<0.05);(3)脾酪氨酸激酶启动子全长和未甲基化启动子荧光素酶的活性分别是甲基化组的4.5和4.7倍;5-Aza-CdR可恢复甲基化启动子的活性;(4)5-Aza-CdR可去甲基化而使脾酪氨酸激酶基因重新表达,而且具有时间依从性。结论:结直肠癌细胞中,启动子区域的甲基化导致Syk基因丧失表达,5-Aza-CdR可以去甲基化而恢复脾酪氨酸激酶基因的表达。  相似文献   

16.
目的 探究HIV感染者接受高效抗逆转录病毒治疗失败后的外周血单核细胞基因表达差异和信号通路的表达情况.方法 从GEO数据库中下载了数据集GSE52900,进行差异化基因的筛选,并通过GO、KEGG和PPI等网络分析确定关键基因和重要生物学通路.结果 从高效抗病毒治疗失败的外周血单核细胞基因表达谱中按照P<0.05和|log2(FC)|>1筛选出79个表达差异化基因,包括55个表达上调基因和24个表达下调基因.PPI网络鉴定出5个关键基因:CD4、CCL5、CXCR4、ITGAL、C1 QB.结论 GO和KEGG分析发现差异化基因主要富集在免疫应答通路,提示免疫因素在抗病毒治疗中发挥了重要作用.  相似文献   

17.
背景:铁死亡与骨关节炎发生、发展密切相关,但具体特征基因及调控机制尚不清楚。目的:运用WGCNA及多种机器学习方法识别骨关节炎铁死亡特征基因及免疫浸润分析。方法:从GEO数据库下载骨关节炎相关数据集,同时在Ferr Db网站中获取铁死亡相关基因,采用R语言对骨关节炎数据集进行批次校正、提取骨关节炎铁死亡基因并进行差异分析,对差异基因进行GO功能及KEGG信号通路分析;同时运用WGCNA分析及机器学习(随机森林、LASSO回归及SVM-RFE分析)筛选骨关节炎铁死亡特征基因,并进行体外细胞实验,将软骨细胞分为正常组和骨关节炎组,运用数据集及q PCR验证表达并行相关免疫浸润分析。结果与结论:(1)经批次校正及PCA分析获得骨关节炎基因12 548个,同时获得铁死亡基因484个,进而得到24个骨关节炎铁死亡差异基因;(2)GO分析主要涉及对氧化应激反应、对有机磷反应等生物过程;涉及细胞顶端、顶端质膜等细胞组分;涉及血红素结合、四吡咯结合等分子功能;(3)KEGG分析显示,骨关节炎铁死亡差异基因与白细胞介素17信号通路、肿瘤坏因子信号通路等信号通路有关;(4)运用WGCNA分析及机器学习筛选...  相似文献   

18.
目的探索模拟失重对小鼠骨髓来源巨噬细胞基因表达的影响。方法采用尾悬吊法构建模拟失重小鼠模型,对照组不做尾吊处理。模拟失重28 d后提取小鼠骨髓来源巨噬细胞并将其分别诱导为M0型和M1型巨噬细胞。随后,提取细胞总RNA进行转录组测序及分析,并通过实时定量RT-PCR验证重要基因的差异表达。结果转录组测序结果显示,模拟失重组与对照组之间的M0型巨噬细胞共有807个基因差异表达;模拟失重组与对照组之间的M1型巨噬细胞共有898个基因差异表达。GO分析表明,M0型巨噬细胞差异表达基因主要集中于免疫应答、趋化等生物学过程中;M1型巨噬细胞差异基因主要富集于炎症反应、细胞迁移的正调控和单核细胞趋化等生物学过程中。KEGG通路分析发现M0型巨噬细胞的差异基因主要涉及趋化因子、细胞因子-细胞因子受体相互作用等信号通路;M1型巨噬细胞的差异基因主要涉及造血细胞谱系、流体剪切应力等信号通路。进一步分析发现,尾吊组与对照组在M0或M1型巨噬细胞中的差异基因均在细胞因子-细胞因子受体相互作用信号通路中富集,其中尾吊组的M0与M1型巨噬细胞中调节单核细胞/巨噬细胞迁移和浸润的关键趋化因子Ccl2均显著高表达,并进一步通过RT-PCR实验得到验证。结论模拟失重可显著影响小鼠巨噬细胞中以Ccl2为代表的与炎症发生和加重密切相关的基因表达水平,这些改变的基因富集于多个生物学过程,可能对巨噬细胞的黏附、迁移等功能产生影响。  相似文献   

19.
目的:通过对癌症基因组图谱(TCGA)数据库中大样本结直肠癌组织基因表达谱RNA测序(RNA-Seq)数据进行生物信息学分析,挖掘与结直肠癌预后密切相关的新基因,并对筛选出的新基因进行验证和功能研究。方法:从TCGA数据库下载结直肠癌RNA-Seq数据筛选差异表达基因,采用R-survival包对差异表达基因进行生存分析,筛选出与结直肠癌预后相关的基因,并从中挑选尚未在肿瘤研究中报道的表达上调最显著的新基因进行验证。采用RT-qPCR检测新基因在人结肠癌细胞、人正常结肠上皮细胞、人结直肠癌组织及配对的癌旁组织中表达水平的变化。采用慢病毒载体构建稳定沉默新基因表达的结肠癌细胞,以转染阴性对照慢病毒的结肠癌细胞为对照,采用CCK-8和Transwell实验研究稳定沉默新基因表达对结肠癌细胞活力、迁移和侵袭的影响。结果:(1)筛选出结直肠癌差异表达基因4 017个,Kaplan-Meier生存分析显示69个基因与结直肠癌患者预后差密切相关(P<0. 01),其中表达上调基因36个,这36个基因中有11个尚未在肿瘤研究中被报道。(2)表达上调倍数前3位的基因为CCDC78、PGGHG及T...  相似文献   

20.
为探讨外周血T细胞及双阴性T细胞(double-negative T cell,DNT)在结直肠癌(colorectal cancer,CRC)患者中的变化及临床意义,收集93例初次就诊的CRC患者(其中55例行CRC根治术)外周血,同时收集54例体检中心健康者的外周血.采用FACS检测T细胞亚群及DNT比例,分析CR...  相似文献   

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

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