中国结核分枝杆菌群体基因组的分析

中国结核分枝杆菌群体基因组的分析

作者:师大云端图书馆 时间:2015-09-19 分类:参考文献 喜欢:2777
师大云端图书馆

【摘要】结核分枝杆菌作为结核病的病原菌,自上世纪五十年代以来,不断发现有效的抗结核菌药物,使古老的结核病流行得到了一定的控制,但是由于不少国家对结核病的忽视,艾滋病的并发感染和抗生素的滥用等,结核病死灰复燃。在中国、印度等地出现了越来越多的耐药结核分枝杆菌,作为中国盛行的lineage2型别也被认为具有更强耐药性、更强毒力的现代型结核菌,这也使得对结核流行病学的、耐药相关的、lineage2型别相关的研究成为关注热点。本实验室通过在中国地区广泛的收集临床耐药结核菌株等,运用高通量二代全基因组测序,对161株结核分枝杆菌进行生物信息分析,从基因组水平挖掘结核分枝杆菌的进化过程、型别特征和耐药相关区域等。本研究的主要结果如下:1.通过与NCBI数据库中22株已知分型的结核分枝杆菌复合群构建系统发育树,确定了收集的161株中国临床菌株分别属于lineage2(122株)、lineage3(2株)和lineage4(37株)。在群体结构方面重点对主要流行的lineage2和lineage4展开分析,从分类比较菌株中突变位点数量发现,结核菌株型别间的差异要比耐药与敏感菌株间的差异要更具有统计学意义,即不同型别间存在明显的突变位点数量的差异。2.基于推测的祖先序列,统计属于lineage2型别相关的位点为295个,lineage4型别相关的突变位点为167个。从G:C>T:A的颠换与G:C>A:T转换的比值中得知,属于lineage2型别相关的SNPs为氧化损伤导致的突变比例较lineage4要大,而氧化损伤又被报道同结核菌的耐药性相关,因此lineage2基因组易发生氧化损伤的突变也为lineage2易突变和易耐药的原因提供条可能线索。另外,我们发现,相对于lineage4型别相关nsSNPs,lineage2型别相关的nsSNPs存在一定数量的影响复制,重组和修复相关的和转录相关蛋白功能的非同义突变(如recD,Rv0922,dinG,uvrC和nth等),这类突变正好也可以用来作为解释lineage2相较于lineage4的易突变和耐多药性的潜在突变机制。因此,这两者型别相关的突变位点除了可作为将来高分辨率分型技术的分子标记,对于研究lineage2的易耐药性、较高突变率和毒力具有一定的价值。3.在依据系统发育树推测结核分枝杆菌的起源和进化史时,我们发现结核分枝杆菌复合群的分化过程同现代人类的迁徙有着密切的联系。追踪结核分枝杆菌复合群的共同祖先,同样是七万年前左右走出非洲,并依此分化非洲分枝lineage5-6(6万年前左右),环印度洋分枝lineage1(5.8万年前左右)、欧美分枝lineage4(3.7万年前左右),并最终3.4万年前左右分化出中亚分枝lineage3和东亚分枝lineage2。4.研究我国主要型别lineage2和lineage4在中国境内的可能历史传播轨迹时,我们推测亚型L2.1、L2.2、L2.3和L4.1、L4.2均是在一万年前左右(新石器时代)发生了扩张和分化出来。其中,L2亚型中相对古老的L2.1很可能同东亚人早期的“南线”从中国的南部地区出现有南至北的传播。而L2.2、L2.3、L4.1和L4.2更可能与东亚人的“北线”迁徙路径有关,从中亚或西伯利亚等区域传入中国并扩散,尤其是最“年轻”的L2.3的共同祖先(5.2千年前-1万年前)很可能由于所处的年代(华夏五千年)或本身的进化优势(更强的毒力、适应性),使L2.3成为如今中国地区的主要流行的型别。5.在对临床结核分枝杆菌的耐药关联分析中,首先根据进化树去除与系统发育树相关的或者是发生同义突变的SNPs,然后通过泊松分布(P<0.05)和正态分布(P<0.01)计算,从剩余的突变基因/间隔区中筛选得到一批多突变位点、耐药菌株中突变比率高的候选耐药相关基因(85个)/间隔区(32个),其中包含已知报道过的14个耐药相关基因和4个耐药相关间隔区。由于临床上有联合用药的习惯,初步统计分析发现该药物与药物之间的关联造成突变频率高的那些已知耐药基因与八种药物都表现出一定关联性,造成关联分析的干扰。6.针对已知报道耐药基因的进一步相关性分析(OR值和Fisher精确检验)发现,生长必需基因突变位点具有一定的偏好性,如rpoB的密码子531和526突变,katG的315突变,这类集中突变的位点很可能是因为对生长必需基因带来的影响小,带给菌株的适应性代价低。相比之下,非必需基因ethA、pncA则表现为耐药菌株中突变的位点较分散、种类多,无明显的特异突变位点;已知报道的inhA和furA基因本身的突变则与异烟肼的耐药性关联性较弱,然而与异烟肼结构上类似的乙硫异烟胺的耐药性则与inhA基因和ethA基因的突变表现出了强关联性;作为同样作用于蛋白质合成的一线药物链霉素和二线药物卷曲霉素、卡那霉素相关的耐药基因则所差异,耐链霉素的菌株以rpsL的突变为主,而编码核糖体16SrRNA的rrs上1400位置的突变则更容易导致卷曲霉素和卡那霉素的耐药性;embA和gyrB则分别在乙胺丁醇和氧氟沙星中耐药性关联远不如embB和gyrA强。7.在挖掘候选耐药相关基因与药物之间的关系中发现,基因水平与耐药机制的关联性要比原本认为的要复杂,基于STITCH数据库分析,除去已知报道的直接药物靶基因突变,还包括间接的耐药补偿相关、细胞壁相关通路(fadD,pks和mmpL家族等)的基因;利用dN/dS分类计算,发现候选的耐药相关基因在耐药菌株中受到强烈的正向选择(平均dN/dS>1.0,其中敏感菌株的全基因组平均dN/dS仅0.430);85个候选耐药基因中含24个结核分枝杆菌生长必需基因,这也很可能与长期的药物刺激和强烈正向选择压力有关;耐药的候选间隔区存在通过sRNA或者新的启动子区域突变来进行调控基因表达的可能。综上所述,本研究通过整合中国地区临床菌株的测序结果,找到lineage2与lineage4型别相关的SNPs和一批耐药候选区域,展示了结核菌同人类长期共进化、迁徙的关系以及在中国的历史传播轨迹,挖掘了一些新结核分枝杆菌的可能耐药机制。本研究为更好预测结核未来的可能传播模式并制定防控、跟踪、诊断和治疗方案提供帮助,同样也为多组学的结核分枝杆菌研究奠定基础。
【作者】林楠;
【导师】汪世华;毕利军;
【作者基本信息】福建农林大学,微生物学,2014,博士
【关键词】全基因组;结核分枝杆菌;单核苷酸多态性;耐药性;进化起源;

【参考文献】
[1]金富贵.中国人口过度和影响消费需求2010-2060[D].复旦大学,区域经济学,2012,硕士.
[2]俞霖.政府监管视野下的食品安全问题研究[D].江西财经大学,行政管理,2013,硕士.
[3]边馥苓,王金鑫.论数字城市工程及其技术体系[J].武汉大学学报(信息科学版),2004,12:1045-1049.
[4]柳维林.miRNAs参与ApoE(-/-)小鼠Toll样受体4免疫作用机制的干预研究[D].湖南师范大学,运动人体科学,2011,硕士.
[5]尹嘉誉.辽河文化美学特征与艺术设计创新取向研究[D].沈阳师范大学,美术学,2013,硕士.
[6]周新原.马克思“人的本质”理论的方法论研究[D].华侨大学,哲学,2014,硕士.
[7]施高萍.多柔体系统动力学模型的降阶研究[D].浙江工业大学,固体力学,2004,硕士.
[8]黄娟.Fe_(16)N_2的物性研究[D].华东师范大学,纳米物理学,2013,硕士.
[9]王璇.基于羧基咪唑鎓两性离子配体的羧基—叠氮混桥配位聚合物的合成、结构与磁性[D].华东师范大学,无机化学,2013,硕士.
[10]陈星运,贺江平,舒远杰.聚合物/石墨导电纳米复合材料研究进展[J].化工新型材料,2010,10:12-14+64.
[11]熊建英.荧光定量PCR快速检测登革1型病毒[D].南方医科大学,病原生物学,2012,硕士.
[12]党朵.秦艽花抗乙酰胆碱酯酶活性成分及伊贝母提取物生物活性研究[D].兰州理工大学,食品科学,2013,硕士.
[13]杨长河.重庆市乡镇基层党员干部队伍建设研究[D].西南大学,马克思主义中国化研究,2014,硕士.
[14]汪恒毅.低功耗Gm-C带通滤波器的研究与设计[D].杭州电子科技大学,微电子与固体电子学,2014,硕士.
[15]衣德龙,郑慕军.VERSATEC静电绘图仪绘图中出现白色条带的原因及排除方法[J].石油地球物理勘探,1994,S2:234-236+238.
[16]常亮.顺铂对小鼠耳蜗钙蛋白酶表达的影响[D].辽宁医学院,生理学,2012,硕士.
[17]李娅.鄂尔多斯盆地铁边城地区延长组长6储层综合评价[D].西安石油大学,油气田开发工程,2012,硕士.
[18]宋印海.直流发电机短路后极性改变的分析[J].冶金自动化,1992,01:52-53.
[19]程潜.考虑列车纵向作用的高架车站动力分析[D].北京交通大学,2014.
[20]杨克文.城市隧道穿越人防空洞关键技术研究[D].北京交通大学,桥梁与隧道工程,2013,硕士.
[21]肖培亮.酶催化不对称合成手性咖啡酸衍生物的研究[D].吉林大学,生物化学与分子生物学,2013,硕士.
[22]谭海伦.论《黑暗的心》中的日神冲动与酒神冲动[D].湖南师范大学,欧洲语言文学,2013,硕士.
[23]田峰.国有土地使用权的权利冲突与解决[D].吉林大学,民商法学,2004,硕士.
[24]张杰一.基于P6的大科学工程项目管理信息系统的设计与开发[D].哈尔滨工业大学,机械工程,2013,硕士.
[25]郭佳琦.股指期货与现货市场联动性及引导关系研究[D].山西财经大学,金融学,2013,硕士.
[26]李义明,王舒驰.从“控网”、“用网”到“生活在网中”——媒介化社会的网络思想政治教育创新[J].思想教育研究,2014,07:53-55.
[27]扶小红.洞庭湖区土地利用/覆盖变化及对洪涝灾害易损性影响分析[D].湖南师范大学,土地资源管理,2014,硕士.
[28]王金磊.垃圾渗滤液回灌的最佳物理吸附条件研究[D].大连海事大学,环境工程,2014,硕士.
[29]田海姣,王铁龙,王颖.垂直轴风力发电机发展概述[J].应用能源技术,2006,11:22-27.
[30]李志敏.起爆药静电响应规律与安全设计[D].北京理工大学,2014.
[31]洪威.全量程谱域光学相干层析成像方法研究[D].浙江大学,光学工程,2013,硕士.
[32]郭翌.心脏占位性疾病的超声图像分析与处理方法研究[D].复旦大学,医学电子学,2013,博士.
[33]谢朝霞.基于组件的智能HRMS研究与开发[D].华北工学院,2003.
[34]易虹天地电脑技术服务(北京)有限公司潘维民博士.数据仓库会成为CRM的灵魂吗?[N].中国计算机报,2001-03-01A15.
[35]相皛.中国建设银行信用卡业务中欺诈风险的防控研究[D].云南大学,工商管理,2012,硕士.
[36]彭贺明.基于工序要求的夹具定位方案规划软件的研制[D].杭州电子科技大学,机械设计及理论,2013,硕士.
[37]郑婷婷.强化层叠式垂直流人工湿地脱氮效果的研究[D].扬州大学,市政工程,2012,硕士.
[38]卓蕊潋.新型便携式多道伽马能谱仪的研制[D].中国地质大学(北京),电子与通信工程,2013,硕士.
[39]李坤.新资源食品人参主要成分及其质量标准研究[D].吉林大学,生物工程,2014,硕士.
[40]常丽君.首秦制氧站设备故障排除及可靠性管理[D].燕山大学,机械工程,2012,硕士.
[41]张军华,吕宁,田连玉,雷凌.地震资料去噪方法综合评述[J].石油地球物理勘探,2005,S1:121-127+138.
[42]程红波.黑龙江移动通信公司供应链发展规划综合管理研究[D].吉林大学,项目管理,2014,硕士.
[43]刘卫.预加固对软弱围岩隧道掌子面稳定性的影响研究[D].北京交通大学,2013.
[44]史文东,张立昌,尹克敏,孙海宁.地质统计地震反演方法及其在民丰断裂带封堵中的应用研究[J].石油地球物理勘探,2002,04:407-411+432.
[45]谢纪刚,赵立彬.融资约束、现金持有量与并购支付方式[J].北京交通大学学报(社会科学版),2014,03:.
[46]李会波.内蒙古师范大学体育学院足球课程资源的利用现状和对策研究[D].内蒙古师范大学,体育教育训练学,2012,硕士.
[47]徐东超.(S)-N-(2-苄氧基—丙烯基)-肼基甲酸叔丁酯的合成与表征[D].大连工业大学,材料学,2012,硕士.
[48]唐国富.山水画语境对水彩风景画创作的影响[D].湖南师范大学,美术学,2014,硕士.
[49]王福桢.北方浓香型白酒发酵微生物多样性分布模式解析[D].内蒙古大学,生物工程,2014,硕士.
[50]林天发.基于非均匀性的沥青路面质量控制研究[D].长安大学,交通运输工程(专业学位),2014,硕士.

相关推荐
更多