2019年3-基因工程原理ds-测序技术.ppt

上传人:上海哈登 文档编号:2820567 上传时间:2019-05-22 格式:PPT 页数:91 大小:7.08MB
返回 下载 相关 举报
2019年3-基因工程原理ds-测序技术.ppt_第1页
第1页 / 共91页
2019年3-基因工程原理ds-测序技术.ppt_第2页
第2页 / 共91页
2019年3-基因工程原理ds-测序技术.ppt_第3页
第3页 / 共91页
2019年3-基因工程原理ds-测序技术.ppt_第4页
第4页 / 共91页
2019年3-基因工程原理ds-测序技术.ppt_第5页
第5页 / 共91页
点击查看更多>>
资源描述

《2019年3-基因工程原理ds-测序技术.ppt》由会员分享,可在线阅读,更多相关《2019年3-基因工程原理ds-测序技术.ppt(91页珍藏版)》请在三一文库上搜索。

1、基因工程原理,何光源 2016-3-10,第三章 测序技术,第一节 DNA测序策略 第二节 传统的DNA序列分析法 第三节 第二代测序技术简介 第四节 第三代测序技术,2019/5/22,造福人类的HGP,人类基因组计划(Human Genome ProjectHGP)于1990年正式启动,其主要目标有:识别人类DNA中所有基因(超过10万个);测定组成人类DNA的30亿碱基对的序列;将这些信息储存到数据库中;开发出有关数据分析工具;致力于解决该计划可能引发的伦理、法律和社会问题。 人类基因组计划是当代生命科学一项伟大的科学工程,它奠定了21世纪生命科学发展和现代医药生物技术产业化的基础,具有

2、科学上的巨大意义和商业上的巨大价值。,2019/5/22,第一节 DNA测序策略,DNA的测序策略因待测DNA分子性质、测序方法和测序目的不同而有所不同。在测序之前,必须根据待测序列区的长度,所要求的测序精确度以及现有设施来制定测序总策略。 一、已知序列的确证性测序策略 二、未知序列的从头测序策略,2019/5/22,一、已知序列的确证性测序策略,确证性测序(confirmatory sequencing):对已知的核酸序列进行鉴定和证实; 例如:临床分子诊断中,对有遗传倾向的病例检测相关基因有无突变; 不仅有助于临床诊断,还有助于探索有关疾病的发病机理及基因突变引起的表型变化。,2019/5

3、/22,一、已知序列的确证性测序策略,多数情况下,确证性测序的DNA片段较小,一套反应即可获得双链DNA其中一条链上局部区域的核苷酸序列,通常只需对亚克隆于M13噬菌体或噬菌粒载体上的一段合适的限制酶切片段进行测序; 待测DNA片段位于通用引物的测序范围之内时,可按克隆位点两侧的载体序列设计和合成寡核苷酸片段作为“通用引物”; 否则,最好的方法是合成一段长度为18-20核苷酸的寡核苷酸引物,与距离待测区约50-100核苷酸的序列互补。,2019/5/22,(一)随机测序法,鸟枪法shotgun strategy 利用超声处理、核酸酶(DNase)切割或限制性酶切割,将长链DNA随机断裂成适于测

4、序的片段,把这些DNA片段装入适当载体,建立亚克隆文库,含有子片段的亚克隆扩增后分别进行DNA序列测定,然后将序列重叠排列,确定待测DNA的完整序列; 该法所获得的序列由许多序列累积产生,结果准确; 但由于测序的亚克隆是随机挑选出来,某些序列可能多次重复测定,而一些序列一直不出现,通常要测定比实际靶DNA所含核苷酸多47倍的序列;如果靶DNA含较多的重复序列,计算机将难以进行分析。,2019/5/22,(二)定向测序法,定向测序法(directed sequencing):指从待测DNA片段的某一端开始按顺序进行测序,直至将待测DNA的序列全部测完; 该类方法具有方向性,不会出现大量重复测定的

5、现象; 主要包括:嵌套缺失法和引物步入法。,2019/5/22,(二)定向测序法,1嵌套缺失法 嵌套缺失法(nested deletion):利用工具酶酶解待测DNA,只要控制消化时间,即可产生一组从同一端缺失的长度不一的互套缺失突变体,测定其序列,通过相互重叠部分将相邻片段的序列彼此拼接,如此重叠互套,获得待测DNA的全序列; 产生互套缺失突变体的主要工具酶:核酸外切酶III、核酸酶Bal 31、DNA酶I和T4 DNA聚合酶。,2019/5/22,(二)定向测序法,用外切酶构建互套缺失突变体测序方法示意图,2019/5/22,(二)定向测序法,2引物步入法 引物步入法(primer wal

6、king ):从待测DNA片段的3端开始,利用载体上的序列设计第一次反应的引物,测定一段DNA序列,然后依次根据前一次测序结果,设计下一次反应引物,循序渐进测序,直至完成,得到靶DNA的全部序列。,2019/5/22,(二)定向测序法,引物步入法测序原理示意图,2019/5/22,Key Genomics Technologies,1975 - Southern DNA hybridization technique 1977 - Sangers chain-termination and Maxam、Gilberts chemical DNA sequencing methods 1980

7、- Automated in situ oligonucleotide synthesis instrument 1985 - Mulliss discovery of PCR at Cetus 1992 - Affymetrix (Fodors group) first gene-chip 1995 - ABIs first automated DNA sequencer 2006 - 2nd generation DNA sequencer on market 2007 and beyond Single molecule sequencing techniques,2019/5/22,测

8、序技术的发展,双脱氧末端终止法 (Sanger 测序法) 1970s 同位素标记,手工 1980s 荧光标记,自动 1990s 毛细管电泳 合成测序法(第二代测序) 焦磷酸测序(Pyrosequencing, Roche/454) 合成测序(Sequencing-By-Synthesis, Illumina/Solexa) 连接测序(Sequencing-By-Ligation, ABI/SOLiD) 单分子测序技术(第三代测序) Helicos Pacific Biosciences Oxford Nanopore,2019/5/22,第二节 传统的DNA序列分析法,Sanger双脱氧链终止

9、法 Maxam-Gilbert DNA化学修饰法 DNA杂交测序 毛细管电泳 DNA序列分析的自动化,2019/5/22,1970年发明了第一种DNA测序方法,2019/5/22,Sanger双脱氧链终止法,Sanger于1977年在加减法测序的基础上创建了双脱氧链末端合成终止法(chain termination method),简称Sanger法、双脱氧法或酶法。,2019/5/22,一、基本原理,单链DNA模板、DNA合成引物、DNA聚合酶以单链的DNA为模板,合成出准确的DNA互补链 能够利用2,3-双脱氧核苷三磷酸作底物,参入到寡核苷酸链的3-末端,从而终止DNA链的延长,2019/

10、5/22,双脱氧链终止法的测序基础是以双脱氧核苷三磷酸为测序反应的链终止剂。,2019/5/22,2019/5/22,2019/5/22,在4组相互独立的测序反应体系中,分别加入四种不同ddNTP,其余成分相同。调整每个测序反应中dNTP与ddNTP的比例,使引物的延伸在对应于待测模板DNA每个可能掺入的位置都有可能发生终止。 产生4组分别终止于互补链3末端的每一个A、G、C和T位置上的一系列长度的核苷酸链。通过高分辨率变性聚丙烯酰胺凝胶电泳和放射自显影检测,直接读出新合成DNA链的序列。,2019/5/22,2019/5/22,双脱氧链末端终止法测序原理示意图,2019/5/22,二、测序反

11、应体系的组成,(一)待测模板 1单链模板DNA Sanger法经典测序反应-将待测DNA片段克隆于M13mp载体中,得到单链的DNA模板,再按Sanger法进行测序。 2双链模板DNA 将待测的DNA片段克隆到质粒载体上,得到含待测DNA克隆片段的双链质粒,经碱变性处理,与寡核苷酸引物序列退火,按Sanger法进行测序;,2019/5/22,Sanger双脱氧-M13体系,M13-噬菌体载体,可得到单链DNA序列 Sanger双脱氧-pUC体系 pUC-质粒载体,利用双链质粒DNA模板的双脱氧DNA序列分析法,2019/5/22,2019/5/22,2019/5/22,2019/5/22,20

12、19/5/22,二、测序反应体系的组成,(二)引物 “通用”引物:在DNA聚合酶催化的测序反应中需要测序引物。不论是单链DNA模板,还是双链DNA模板,都可通过使用克隆位点两侧的载体序列互补的“通用”引物; 通用引物长度:一般为1530个核苷酸。,2019/5/22,二、测序反应体系的组成,(三)DNA聚合酶 1大肠杆菌DNA聚合酶I大片段(Klenow片段): Sanger法最早使用的酶,具有53聚合酶和35外切酶活性,催化链延伸反应的能力较差,用它进行测序常产生较高的本底和假带; 2测序酶(sequenase):经过修饰的T7噬菌体DNA聚合酶,消除了35外切酶活性。酶活非常稳定,很高的链

13、延伸能力和极快的聚合反应速度,测定较长DNA的首选酶; 3Taq DNA聚合酶:PCR反应关键酶,在95的高温下保持稳定,可在高温下进行反应,该酶用于测序具有很好的链延伸性能,能克服富含GC序列的模板形成自身二级结构对测序的影响。,2019/5/22,二、测序反应体系的组成,(四)放射性核素标记的dNTP -32P-dNTP或-35S-dNTP。目前通常采用-35S-dNTP。,2019/5/22,二、测序反应体系的组成,(五)测序产物的凝胶电泳及识读 能否将测序反应中产生的各种不同长度的DNA片段进行有效分离是序列分析成败的关键。,2019/5/22,三、Sanger双脱氧链终止法的特点,1

14、.快速简便,一次反应可测500个以上的碱基序列; 2.荧光染料代替放射性核素标记,使该法便于实现自动化分析,能够满足绝大多数DNA样品序列分析的需要。,2019/5/22,Maxam-Gilbert DNA化学修饰法,1977年哈佛大学A.M. Maxam和W. Gilbertf发明 原理: 化学试剂处理末端放射性标记的DNA片段,造成碱基的特异性切割产生一组具有不同长度的DNA链混合物经凝胶电泳分离和放射自显影,读出待测DNA片段的核苷酸序列,2019/5/22,由于这种化学切割反应的特异性是由碱基的修饰作用决定的,所以切割反应必定是定量的 化学切割反应的试剂: 肼 hydrazine 联氨

15、 NH2-NH2 硫酸二甲酯 dimethylsulphate 六氢吡啶,2019/5/22,碱性条件下,肼与胸腺嘧啶(T)和胞嘧啶( C )作用 通过六氢吡啶作用,使两个磷酸分子从糖片段上释放出来,导致在核苷酸位置上发生DNA的断裂 盐存在的条件下,肼同T的反应被抑制,只发生胞嘧啶碱基特异的切割反应 由此区别C和C+T两种反应,2019/5/22,硫酸二甲酯((CH3O)2SO4) 一种碱性的化学试剂 作用于DNA碱基环中的氮原子,使之甲基化 在中性的pH值环境中,可导致配糖键发生水解,使去碱基的糖-磷酸键十分微弱 碱性条件下磷酸分子从DNA链上脱落,造成链的断裂 鸟嘌呤(G)的N7腺嘌呤(

16、A)的N3 六氢吡啶作用,从脱氧核糖上移去2个磷酸分子,使DNA链在甲基化的鸟嘌呤G位点断裂,2019/5/22,2019/5/22,CS载体系统:末端标记载体 核酸内切酶Th111 特点: 产生5单碱基的突出末端,便于标记; Th111位点十分稀少; 5突出末端可以是G或A,也可以是T或C,选择性强; 在载体中具有两个Th111 位点,可对克隆在载体上的片段任何一段作选择性标记,2019/5/22,化学修饰法的优点: 不需要体外酶切; 只要有末端标记,无论单双链都可用来测序。 限制因素:测序胶的分辨率,2019/5/22,DNA杂交测序,原理:如果一段短的DNA探针能够与较长的靶DNA片段杂

17、交,并形成完全的双链分子,可据此推断靶DNA序列相对应的互补序列 步骤: 将待测的靶DNA分子与一组已知核苷酸序列的寡核苷酸探针进行杂交 比较分析与待测DNA杂交的探针之间的碱基重叠关系,据此推算靶DNA的核苷酸序列,2019/5/22,2019/5/22,两种操作方式: 1. 将不同的寡核苷酸与固定在滤膜上的靶DNA序列样品进行杂交 2. 应用寡核苷酸矩阵芯片的DNA杂交测序法,2019/5/22,杂交测序的应用,1.检测靶DNA的单碱基突变 2. 用于不同片段之间的序列比较分析 (同源性序列检测) 3. 用表达序列标签(EST)的矩阵芯片,检测不同类型细胞中或不同生长发育状态下的细胞中特定

18、基因的表达状况 4. 检测传统的测序技术所得DNA片段的核苷酸序列数据,2019/5/22,毛细管电泳,以高压直流电场为驱动力,在毛细管内使带电粒子按淌度或分配系数进行分离的一种电泳技术。 具有分辨率高、重现性好、灵敏度高、快速和易于实现自动化等优点。 常用于DNA测序的毛细管电泳形式有以下几种: (一)毛细管凝胶电泳 (二)非凝胶基质毛细管电泳 (三)阵列毛细管电泳,2019/5/22,(一)毛细管凝胶电泳CGE,Capillary gel electrophoresis 将平板电泳的凝胶移到毛细管中做支持物进行电泳,由于电场强度比板凝胶电泳大大提高,使分离时间缩短约25倍; 主要以聚丙烯酰

19、胺凝胶作筛分介质,它黏度大、抗对流,能减少溶质的扩散,有极好的分离效果,除用于DNA序列分析外,还可用于DNA片段的分离和PCR产物的分析。 问题:制备凝胶柱费力;凝胶形成和电泳期间产生的气泡影响分离;使用过程中焦耳热对分离介质的降解使毛细管的寿命缩短等。,2019/5/22,(二)非凝胶基质毛细管电泳 为了避免CGE存在的问题,非凝胶基质毛细管电泳用于DNA分析的研究越来越多; 非凝胶基质:线性聚丙烯酰胺和纤维素衍生物 非凝胶基质在DNA测序中可以更换,使毛细管的寿命延迟,在优化条件下可获得高效及高速的分离效果。,2019/5/22,(三)阵列毛细管电泳CAE,Capillary array

20、 electrophoresis 将毛细管电泳与板凝胶电泳的优势相结合,采用毛细管凝胶电泳的装置,将多支毛细管并列进行分离与检测,以板凝胶电泳的优势弥补了毛细管凝胶电泳的不足,可一次检测多个样品; 阵列毛细管电泳散热效率高,适于高电场强度电泳,是一种高速、高通量的测序法; 使用非凝胶基质,毛细管壁可以抑制横向扩散,具有易于实现自动化的优点。,2019/5/22,DNA序列分析的自动化,激光测序法终止标记系统 用4种不同的荧光染料标记不同的ddNTP; 测序反应在同一反应管内进行,不必分成4管; 反应产物按终止位置的碱基不同其3末端带有不同的荧光基团,被激发后产生不同的荧光; 将反应产物加样于凝

21、胶的同一加样孔,电泳分离后,经过DNA测序仪分析系统识别,将检测将信号不断传送到计算机,通过软件分析,自动读出待测DNA的全部核苷酸序列。,2019/5/22,激光测序法终止标记系统测序原理示意图,2019/5/22,测序仪原理,1.荧光标记 BigDye Terminators-美国应用生物系统公司专利,四色荧光分别标记的起链终止剂作用的四种ddNTPs; 一个反应管中完成循环测序反应,通过一个电泳道电泳即可完成测序任务; 传统双脱氧链终止法:四个反应管循环测序反应,四个泳道电泳检测,2019/5/22,2.循环测序,2019/5/22,377型遗传分析仪,310型遗传分析仪,3100遗传分

22、析仪,3730遗传分析仪,2019/5/22,第三节 第二代测序技术简介,第二代测序技术,2019/5/22,Next Generation Sequencing,2019/5/22,表 目前使用最广泛的三大第二代测序平台测序能力统计信息(2010年年初数据),Solexa和SOLiD配对末端测序所需时间和产出是单末端的两倍,454的配对末端和单末端差异在于建库方法,所需时间和测序量不变。ABI SOLiD包含两张芯片,这里的数据是一张芯片的量。,2019/5/22,Roche 454 焦磷酸测序,原理:循环芯片测序法cyclic-array sequencing 对布满DNA样品的芯片重复进

23、行基于DNA的聚合酶反应(模板变性、引物退火杂交及延伸)以及荧光序列读取反应。 优势:操作更简易、费用更低廉、广泛应用。,2019/5/22,原 理,荧光素酶,硫酸化酶,焦磷酸,荧光素,氧化荧光素,可见光,2019/5/22,2019/5/22,测序原理,体外构建好的两端带接头的单链DNA文库; 单链DNA文库在一个油包水的乳液环境中进行PCR扩增 由于起始阶段模板浓度非常低,因此大部分包含磁珠的乳液滴都只含有一种模板; 经PCR扩增,每个磁珠只连有一种模板的扩增产物 打破为乳液滴,收集磁珠,加入芯片中;,2019/5/22,测序原理,经乳液PCR法扩增后携带有大量模板的磁珠被置于芯片上的微孔

24、中; 焦磷酸法测序 每一轮反应都会掺入一个核苷酸,同时释放一个焦磷酸; 焦磷酸和腺苷酰硫酸在硫酸化酶的催化下生成ATP; 荧光素酶在ATP参与下将荧光素转化为氧化荧光素,同时发出荧光被检测器检测到。,2019/5/22,工作流程,文库制备,乳液PCR,焦磷酸测序,加入含酶小微珠,芯片制备,2019/5/22,工作流程,一、样品输入并片段化 基因组DNA、BAC等被打断成300800 bp的片段; 对于小分子非编码RNA或者PCR扩增产物,直接使用 二、文库制备 将A和B接头(3和5端具有特异性)连接到DNA片段上; 接头将用于后续的纯化、扩增和测序步骤; 具有A、B接头的单链DNA片段组成了样

25、品文库。,2019/5/22,工作流程,三、一个DNA片段连接一个磁珠 单链DNA文库被固定在特别设计的DNA捕获磁珠上,每一个磁珠携带了一个独特的单链DNA片段; 磁珠结合的文库被扩增试剂乳化,形成油包水的混合物,成为只包含一个磁珠和一个独特片段的微反应器。 四、乳液PCR扩增 每个独特的片段在自己的微反应器里进行独立的扩增,没有其他的竞争性或者污染性序列的影响; 整个片段文库的扩增平行进行; 每一个片段扩增后产生几百万个相同的拷贝; 乳液混合物被打破,扩增的片段结合在磁珠上。,2019/5/22,工作流程,五、测序 携带DNA的捕获磁珠(20m)随后放入PTP板(Pico Titer Pl

26、ate)的微孔(29 m)中进行后继的测序反应; 反应释放出的光信号实时被仪器配置的高灵敏度CCD捕获到,有一个碱基和测序模板进行配对,就会捕获到一分子的光信号; 由此一一对应,可以准确快速地确定待测模板的碱基序列 六、数据分析 454测序系统可以在10小时的运行当中获得100多万个读长,读取超过4-6亿个碱基信息; 提供两种不同的生物信息学工具对测序数据进行分析,适用于不同的应用:达400 MB的从头拼接和任何大小基因组的重测序。,2019/5/22,454 sequencing: Emulsion PCR (emPCR),Generation of millions of clonally

27、 amplified templates on each bead No cloning and colony picking,2019/5/22,454 sequencing: Deposition of DNA beads into the PicoTiterPlate,2019/5/22,454测序系统图示,A为液体试剂供应装置 B为反应池 C为光线检测成像系统和计算机控制系统,2019/5/22,测序质量,454测序技术的准确率在99%以上; 其主要限制来自同聚物,也就是相同碱基的连续掺入,如AAA或GGG; 由于没有终止元件来阻止单个循环的连续掺入,同聚物的长度需要从信号强度中推断出

28、来。这个过程就可能产生误差; 454测序平台的主要错误类型是插入-缺失。,2019/5/22,图a为高通量鸟枪Sanger测序法策略 图b为鸟枪循环芯片测序法策略,2019/5/22,454测序仪的先行者地位 Leamon、Rothberg等人撰写的一篇介绍该技术的论文被引用了570多次 100多篇经过同行审议的关于人类遗传学、代谢组学、生态学、进化学以及古生物学的论文(peer-reviewed publications)都是使用454测序仪开展的研究,2019/5/22,一个4Mb基因组和3个6Mb基因组测序 传统的Sanger测序法:需要好几个月的时间, 454测序仪,一位实验人员,包括

29、样品制备等步骤在内所用的时间仅需要一周 使用454测序仪还避免了传统测序方法中细菌克隆阶段可能出现的错误,获得了高质量的测序结果 发现了导致结核分枝杆菌对R207910产生抗药性的两个点突变位点,此研究成果使得人类最近的40年内第一次找到了特异性治疗结核病的药物,2019/5/22,2007年6月,James Watson的基因组序列登录到了GenBank数据库当中,这是第一次使用非Sanger测序法获得了人类个体基因组序列,且第一次将个人基因组序列公之于众 整个测序在两个月之内完成,花费不到100万美元,仅占耗时10年之久的人类基因组计划使用经费的千分之一,Venter基因组计划费用的百分之

30、一,2019/5/22,454测序仪最初的技术参数(每次可以获得两千万碱基序列,测序长度100bp,准确率96%) James Watson测序时的技术参数(一亿碱基序列,250bp,99%) 个体基因组测序的费用由100,000美元降低到10,000美元,继而降低到1,000美元甚至更低,2019/5/22,Illumina Solexa 合成测序,Illumina Solexa 合成测序 基本原理,2019/5/22,Clonal Single Molecule Arrays 单分子克隆,Prepare DNA fragments,Ligate adapters,Attach single

31、 molecules to surface Amplify to form clusters,1000 molecules per 1 m cluster 1000 clusters per 100 m square 40 million clusters per experiment,2019/5/22,Reversible Terminator Chemistry 可逆终止反应,All 4 labelled nucleotides in 1 reaction,2019/5/22,Sequencing-by-Synthesis (SBS),Cycle 1: Add sequencing re

32、agents First base incorporated Remove unincorporated bases Detect signal Cycle 2-n: Add sequencing reagents and repeat,1、每轮测序反应加入四种带有荧光标记的dNTP,末端带有可以被去除的阻断基团 2、每轮反应只能整合一个核苷酸,仪器读取相应的荧光信号 3、信号读取结束,用化学方法去除阻断基团,进行下一轮测序反应,2019/5/22,Base calling from the raw data,The identity of each base of a cluster is

33、read off from sequential images 根据每个点每轮反应读取的荧光信号序列,转换成相应的DNA序列,2019/5/22,Solexa 测序 Workflow,2019/5/22,第四节 第三代测序技术,第三代测序技术又称单分子测序技术,DNA在测序时,不需要经过PCR扩增,可以实现对每一条DNA分子的单独测序。 主要分为两类: 单分子荧光测序技术 新型纳米孔测序技术,2019/5/22,单分子荧光测序技术,脱氧核苷酸用荧光标记,显微镜可以实时记录荧光的强度变化。当荧光标记的脱氧核苷酸被掺入DNA链的时候,它的荧光就同时能在DNA链上探测到。当它与DNA链形成化学键的时

34、候,它的荧光基团就被DNA聚合酶切除,荧光消失。这种荧光标记的脱氧核苷酸不会影响DNA聚合酶的活性,并且在荧光被切除之后,合成的DNA链和天然的DNA链完全一样。,2019/5/22,新型纳米孔测序技术,新型纳米孔测序法(nanopore sequencing)是采用电泳技术,借助电泳驱动单个分子逐一通过纳米孔 来实现测序的。由于纳米孔的直径非常细小,仅允许单个核酸聚合物通过,而ATCG单个碱基的带电性质不一样,每一个核苷酸通过孔道时都能出现一种特定形式的电流改变,通过电信号的差异就能检测出通过的碱基类别,从而实现测序。,2019/5/22,第三代测序技术的特点,实现了DNA聚合酶内在自身的反

35、应速度,一秒可以测10个碱基,测序速度是化学法测序的2万倍; 实现了DNA聚合酶内在自身的延续性,一个反应就可以测非常长的序列。二代测序现在可以测到上百个碱基,但是三代测序现在就可以测几千个碱基; 精度非常高,达到99.99%;,2019/5/22,第三代测序技术的特点,直接测RNA的序列,RNA的直接测序,将大大降低体外逆转录产生的系统误差; 可以直接测甲基化的DNA序列,根据DNA聚合酶在甲基化位点与非甲基化时停顿的时间不同,可以判断模板是否甲基化。,2019/5/22,思考题,阐述DNA的测序方法有哪些? 二代测序技术包括哪些类型?各有什么优势?,2019/5/22,2019/5/22,Thank You!,

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 其他


经营许可证编号:宁ICP备18001539号-1