结核分枝杆菌北京型家族菌株的微进化与遗传多样性研究

【摘要】结核病(tuberculosis,TB)是严重危害人类健康的传染性疾病。WHO的统计数字显示,2011年全球共有新发结核病例数870万,包括31万耐多药结核病例(Multi-drugresistanttuberculosis,MDR-TB),因结核病死亡的病例数高达140万。中国是22个结核病高负担国家之一,病例总数居全球第二,耐药病例数世界第一。结核分枝杆菌北京型家族菌株(下称北京菌株)是我国结核病的主要致病菌,该菌株大多分布于东南亚地区,但随着亚洲人口的迁移,北京菌株已遍布全球,是一百多年来扩张最快的结核菌家族。课题组前期通过复制、重组、修复相关的3R基因(Replication,Recombination,andRepair)单核苷酸多态性位点(SingleNucleotidePolymorphism,SNPs)分型技术,将北京菌株进一步细分为7种亚型,并发现最近分化(5,000-6,000年前)的Bmyc10亚型在我国人群中的分布占主要地位(约80%),而种系发生关系上与之相邻的Bmyc26亚型在人群中的比例不足1%。这两个种系关系相近的结核菌亚型为什么会出现如此悬殊的流行差异,是本论文研究的第一个问题。为了揭示Bmyc10菌株在人群中广泛分布的分子机制,我们根据全国6县级地区1448株结核菌的3R基因SNPs以及16位点数目可变串联重复序列(VariableNumberofTandemRepeats,VNTR)基因分型结果,挑选了35株代表5种北京菌株亚型的临床菌株用作高通量全基因组测序。联合4株在全球范围发现并已公布基因组序列的北京菌株和结核菌标准菌株H37Rv,我们构建了北京菌株各亚型的种系关系发生树,以此鉴定Bmyc10亚型特异的SNPs位点和插入缺失序列(InsertionsandDeletions,Indels)。种系发生树的结构显示,目前引起广泛关注的北京菌株均属于Bmyc10亚型。对该亚型突变基因的功能相似性分析结果提示,这些基因在细胞定位上偏好于细胞膜和细胞壁等与宿主频繁接触的区域。通过位点保守性、正选择分析等方法,我们鉴定出11处可能对基因功能有重要影响的遗传突变。参考结核菌在134种应激条件下的转录组学数据,我们发现诱导这些突变基因上调的应激条件为自由脂肪酸以及反应氧/氮化合物损伤,这与结核菌在宿主巨噬细胞内的生存环境十分相似。因此,我们的分析结果提示,Bmyc10亚型特异的遗传突变很可能以多个突变基因协同作用的方式,增强了该亚型菌株对宿主巨噬细胞内应激条件的适应性。这些适应性突变的积累,为此后Bmyc10亚型的扩张奠定了重要的遗传基础。流行病学研究提示,耐多药结核病的产生与传播是结核病在某些地区广泛流行的主要原因,然而,结核菌获得性耐药的形成机制目前尚不清楚。在本论文的第二章,我们将研究方向转向了Bmyc10菌株在临床治疗过程中的微进化过程,其中结核菌获得性耐药的形成过程是我们重点研究的问题。为了研究在治疗压力下,结核菌从全敏感发展为耐多药结核的微进化过程,我们从上海疾病预防控制中心(下称上海疾控中心)选取了3例结核病患者的7份痰菌样本作全基因组测序。经过以似然比检验(LikelihoodRatioTest,LRT)为核心的分析流程,我们发现每例测序样本具有8-41处频率大于5%的不固定突变,提示体内的结核菌群体具有丰富的遗传多样性。接着,在对耐药相关基因更为细致的分析中,我们发现多达4-5种耐药突变在耐药形成过程中被筛选,这说明结核菌有能力选择多种耐药突变应对多药联用的抗结核治疗方案。最终,随着治疗的持续进行,只有一种耐药突变被保留,提示不同的耐药菌亚群之间存在相互竞争(克隆干扰)的关系。结核菌克隆干扰的发现,强调了早期诊断、及时治疗对预防获得性耐药的重要意义。此外,我们还发现了19处在治疗过程中频率发生显著变化的基因组突变,其中的14处突变可能增强了耐药结核菌对体内环境的适应能力。结核菌体内遗传多样性的发现,对其适应性进化的理论研究具有重要意义,然而,这一发现的临床应用价值仍有待探索。第二章的研究中,我们发现有效的治疗方案,使痰菌样本中的特征性突变减少,而当治疗方案无效时,高通量测序检测到的突变数量增加,这提示痰菌样本的遗传多样性可能与结核病患者的疗效相关。随着测序成本的不断下降,全基因组测序有望替代当前的痰液检测手段,成为指导临床用药,监测治疗效果的常规检验方法。定量监测结核菌遗传多样性变化并探索与治疗转归的联系,因此成为第三章研究的重点问题。为了初步探讨治疗效果与痰菌样本遗传多样性的联系,我们从河南省胸科医院的菌种库中选取了1例耐多药复治病例在治疗前,以及治疗后第2、4、6、8、16周的痰菌样本,并对这些样本进行了1,000倍以上的高深度测序。治疗记录显示,在药敏测试结果不详的前4周里,患者接受了无效的治疗方案,而第5周后二线抗结核药阿米卡星的应用,使痰液菌载量略为下降。与此相应,我们在治疗前4周的3例痰菌样本中发现特征性突变从33持续上升至50SNPs/103倍测序深度,而第6周的样本中,突变数量又骤降至25SNPs/103倍测序深度,下降趋势一直延续至第16周。这一结果提示,“治疗无效,多样性上升；治疗有效,多样性下降”的科学假设很可能是成立的。接着,通过鉴定治疗过程中频率发生显著性变化的基因组突变,我们发现仅有一处同义突变的频率持续上升。虽然没有改变蛋白质一级结构,但该突变使编码天冬氨酸密码子的偏好性提升了2倍,说明同义突变可能通过改变密码子使用偏好性影响了突变蛋白的表达水平。最后,我们发现阿米卡星主要的杀菌效应只维持了不到1周的时间。由于治疗期间并没有阿米卡星耐药相关突变出现,我们推测结核菌的表型耐药是治疗失败的主要原因。这一结果提示,高通量技术也可以用于结核菌表型耐药的检测和评估,为将来的个性化治疗提供用药指导。综上所述,本论文通过高通量测序技术,研究了北京菌株Bmyc10亚型在不同时空尺度下的进化过程。通过对痰菌样本遗传多样性的定量监测,我们还初步探讨了高通量测序对指导用药,评估疗效的可行性。本研究的结论不仅对结核菌微进化的理论研究有参考价值,也有望为高通量测序技术的临床应用指出一条新的方向。
【作者】孙刚；
【导师】高谦；
【作者基本信息】复旦大学，病原生物学，2013，博士
【关键词】结核分枝杆菌；北京型家族；高通量测序；遗传多样性；

【参考文献】
[1]钱坤荣.A医院作业成本法下的医疗项目成本管理研究[D].苏州大学,工商管理（专业学位）,2013,硕士.
[2]张赟.山西省对外贸易对产业结构升级转型影响的实证研究[D].兰州商学院,国际贸易学,2014,硕士.
[3]胡钢,刘哲,徐小平,高瑞.像素级图像融合技术的研究与进展[J].计算机应用研究,2008,03:650-655.
[4]胡泱,陈刚.一种有效的基于网格和密度的聚类分析算法[J].计算机应用,2003,12:64-67.
[5]王燕.我国证券业系统性风险研究[D].山西财经大学,金融学,2014,硕士.
[6]付彩欣.WCGDV-Hop无线传感器网络定位算法研究[D].吉林大学,电子与通信工程,2014,硕士.
[7]韦京,杜林林,孙明志.地铁车站PBA工法小导洞施工方案优化研究[J].建筑技术,2014,09:796-799.
[8]谢鲲鹏.化学混凝—生物接触氧化处理高浓度油田采油废水的工艺条件研究[D].辽宁师范大学,细胞生物学,2003,硕士.
[9]李琳.元认知策略在高中生英语阅读理解中应用情况的调查研究[D].沈阳师范大学,外国语言学及应用语言学,2013,硕士.
[10]刘明.创新村级组织党的建设问题研究[D].安徽大学,公共管理,2012,硕士.
[11]R·M·奥尔福德,K·R·凯利,斯文·特雷特尔,R.W.沃德.地震模拟用的有限差分法[J].石油地球物理勘探,1976,S1:28-29.
[12]杜俊卿.氯钾离子共体诱导黄瓜对霜霉病抗性及其酶机理的研究[D].内蒙古农业大学,蔬菜学,2004,硕士.
[13]王萍.运用数据挖掘技术预测客户购买倾向——方法与实证研究[J].情报科学,2005,05:738-741.
[14]徐振伟.基于SOA的烟草物流调度与服务平台的研究[D].浙江理工大学,计算机技术,2012,硕士.
[15]唐国维,申静波,赵建民,李井辉.知识发现和数据挖掘及其在油田生产辅助决策中的应用[J].自动化技术与应用,2003,04:58-60+69.
[16]郑申白.采用差动变压器实现辊缝测量的方法[J].冶金自动化,1996,03:41-43.
[17]耿静漪.负载型二氧化钛光催化降解染料的研究[D].苏州大学,2012.
[18]邹红.同伴互助学习法在ICU护士核心能力培养中的应用[D].南方医科大学,护理学,2012,硕士.
[19]尹承荣.半煤岩巷柔模混凝土沿空留设技术研究[D].西安科技大学,岩土工程,2013,硕士.
[20]翁启献,段裘龙.〈八〉主机选型与轴系设计[J].江苏船舶.1990(01)
[21]李东胜,史作民,冯秋红,刘峰.中国东部南北样带暖温带区栎属树种叶片形态性状对气候条件的响应[J].植物生态学报,2013,09:793-802.
[22]胡昌平.面向用户的资源整合与服务平台建设战略——国家可持续发展中的图书情报事业战略分析(2)[J].中国图书馆学报,2005,02:5-9+24.
[23]赵艳.我国证券公司内部控制问题研究及系统设计[D].对外经济贸易大学,工商管理,2004,硕士.
[24]王扶林.我院第四届科研成果、优秀著作、优秀论文评审工作结束[J].石油管理干部学院学报.1999(01)
[25]方月娥1肖雄亮2.工程机械专业教学模式改革探索与实践[J].中国校外教育.
[26]李芳.基于多智能体的地铁站台乘客紧急疏散模型[D].北京交通大学,2015.
[27]寇蕾.基于职业生涯规划的连锁酒店集团内部营销管理研究[D].沈阳师范大学,旅游管理,2013,硕士.
[28]郑莉.碳纳米管及其复合材料修饰电极的电催化作用研究[D].西北大学,2009.
[29]岳智忠,陈建军,刘永昌.三维综合解释技术在长沙岭构造成藏条件研究中的应用[J].石油地球物理勘探,2002,S1:97-100+221.
[30]张晋春,李广峰.基于数据仓库的站段计划经营决策支持系统的研究[J].铁道运输与经济,2002,08:39-41+47.
[31]曹娅玲.基于纳米材料构建电致化学发光免疫传感器的研究[D].西南大学,分析化学,2013,硕士.
[32]谭锦荣.供应链环境下医药分销企业转型研究[D].华中师范大学,工商管理（专业学位）,2013,硕士.
[33]曹春梅.跳板跳水辅助训练系统的设计与开发[D].清华大学,2004.
[34]李兆楼.基于环糊精的功能性有机超分子凝胶新体系研究[D].山东大学,有机化学,2014,博士.
[35]卢文锦.新媒体冲击下我国纸媒转型策略探析[D].四川省社会科学院,新闻学,2014,硕士.
[36]陆建兵.情态班学生的自我效能感、学习动机与学习成绩的相关与干预[D].山东师范大学,教育（专业学位）,2013,硕士.
[37]褚阳.我国石油行业垄断法律规制研究[D].郑州大学,经济法学,2013,硕士.
[38]沈强.书籍设计的形态品质塑造方式探析[J].大众文艺,2013,06:92-93.
[39]赵旸.乳腺癌MRI表现与ER、PR、HER-2、ki-67相关性的研究[D].吉林大学,临床医学,2013,硕士.
[40]王小会.表面功能化的磁性荧光纳米颗粒的制备、表征及其在胰腺癌靶向双模态成像中的应用研究[D].苏州大学,影像医学与核医学（专业学位）,2014,硕士.
[41]寇琦明.济南市楼盘广告语的语用策略研究[D].曲阜师范大学,语言学及应用语言学,2013,硕士.
[42]娄家星.基于单片机的超低功耗智能遥控车位锁的设计与实现[D].山东大学,计算机软件与理论,2012,硕士.
[43]刘铭.纳米氧化物的盐助溶液燃烧合成与表征[D].南昌大学,2014.
[44]葛珊珊.老年突发性耳聋患者急性应激反应的护理干预研究[D].山西医科大学,耳鼻咽喉科学,2013,硕士.
[45]胡鹏飞.中小学生学习“楷书”对硬笔书写的促进作用[D].曲阜师范大学,美术学,2013,硕士.
[46]董柏举.丝素蛋白纳米纤维—纳米金复合材料的制备及性能研究[D].苏州大学,纺织材料与纺织品设计,2014,硕士.
[47]孙露.计及侧风的某跑车高速气动稳定性研究[D].湖南大学,车辆工程,2011,硕士.
[48]孙青.大豆肽对大鼠血脂代谢的影响及其抗动脉粥样硬化机理研究[D].山东大学,营养与食品卫生学,2013,硕士.
[49]孙安超.刘琦声乐作品创作特征与演唱研究[D].天津音乐学院,音乐与舞蹈学,2014,硕士.
[50]张忠福.农村中小学布局调整问题研究[D].华东师范大学,教育经济与管理,2003,硕士.