结核分枝杆菌全基因组测序
文库制备方法、测序和生物信息学分析技术发展迅速，且成本逐渐减低，为二代测序应用于病原菌感染性疾病诊断和分子流行病学分析带来了革命性的进展。
随着测序技术的日益成熟和生物信息学的发展，全基因组测序和分析已被广泛应用于结核分枝杆菌的分子检测（早期辅助诊断）、抗结核药物敏感性的预测（指导临床用药和制定科学的治疗方案），以及对地区传播模式的分子流行病学研究（有助于发现和控制疫情）中，为解决结核病的早期诊断和耐药问题，以及有效控制疫情传播提供了很大的帮助。
晶诺专业致力于结核病相关的基础和应用研究。针对结核分枝杆菌的基因组特征，建立了针对性、标准化的基因组提取和质控评估平台，已有几万株结核分枝杆菌基因组数据生物信息学分析经验。
下面详细介绍一下晶诺的结核测序、分析平台：
一、质控评估平台
为了从源头上保证测序数据的准确性、可靠性，对样品提取、DNA检测、建库、测序每一个生产步骤都严格把控，从源头上确保了高质量数据的产出。
（1）DNA提取：
基因组提取我们采用改良版的CTAB法，由经验丰富的专业人员处理，以保证DNA的高质量。
（2）建库前质控：
浓度检测：Qubit荧光光度计
完整性检测：琼脂糖凝胶电泳
根据浓度、体积、总量、电泳结果来综合评判样本质量。
（3）文库及测序数据质控：
文库构建完成后，先使用Qubit 2.0进行初步定量，稀释文库，随后使用Agilent 2100对文库的插入片段进行检测，插入片段大小符合预期后，使用Q-PCR方法对文库的有效浓度进行准确定量分析，多个环节并重以保证文库质量。文库检测合格后，按照有效浓度及目标下机数据量的需求将不同文库pooling至flowcell，cBOT成簇后使用illumina高通量测序平台NovaSeq 6000进行测序，原始数据量产出1G以上，测序平均深度可达200X以上。
二、结核菌全基因组测序数据分析平台
（1）耐药分析
基因型预测是基于异烟肼（ahpC、inhA、fabG1和katG）、利福平（rpoB）、乙胺丁醇（embA、embB和embC）和吡嗪酰胺（pncA）等耐药性相关基因内或上游的突变。根据已知的结核耐药突变数据库，预测结核分枝杆菌对常用抗结核药物（异烟肼、利福平、乙胺丁醇、吡嗪酰胺、链霉素、氟喹诺酮、氨基糖苷类、阿米卡星、卡那霉素、卷曲霉素、乙硫异烟胺、对氨基水杨酸、环丝氨酸、利奈唑胺、贝达喹啉、氯法齐明、迪拉马尼）的敏感性，全基因组测序涵盖的基因组位点数量几乎不受限制，但可能无法检测到具有复杂基础基因相互作用的耐药机制。
耐药分析示例：
图1
图2
图1给出了样本菌株的谱系、敏感、耐多药或广泛耐药，以及对应药物耐药的突变点信息，其中-为未发现已知对应的药物耐药基因型。图2为耐药基因的突变频率，突变点在样本中可能有些菌株突变、有些未突变。
注：lineage1为环印度洋分枝、lineage2为东亚分枝（北京系）、lineage3为中亚东非印度分枝、lineage4为欧美分枝、lineage5和lineage6为主要流行于西非的非洲分枝，M. africanum West Africa 1 和 M. africanum West Africa 2、Lineage7 埃塞俄比亚分枝。
（2）遗传距离分析——近期传播关系
传统分子流行病学方法一般使用IS6110 RFLP、spoligotype、VNTR等基因分型方法鉴定结核病传播链，但使用全基因组测序的方法相较于基因分型具有更多的优势，如全基因组测序具有更高的分辨率，全基因组测序是基于基因组间的SNP差异数来确定两两菌株间是否具有近期传播关系，而传统基因分型只是基于部分基因进行区分，结果可能造成基因型相同的菌株存在不同的传播源和传播路径。另一方面，传统基因分型很难准确鉴别传播发生的路径、顺序及是否存在多个传播源，而通过全基因组数据分析可解决此类问题。
遗传距离分析示例：
该图为每两个菌株间的SNP差异数矩阵，一般认为SNP<12个，两菌株间具有近期传播关系，可结合流行病学调查鉴定地区结核病的传播链。
此项分析在流行病学调查研究中更具有明显优势，随着高通量测序的普遍应用及成本优势，有望取代传统的基因分型方法，成为分子流行病学的金标准。
除了以上常规分析内容，晶诺新增个性化分析内容——构建系统发育树，并将以上的数据内容进行汇总图形化显示，更清晰区分各菌株间关系、药物的敏感性、谱系及耐药分类。
系统进化树示例：
图中线长短用于表示两个菌株间的遗传距离，谱系、耐药分类使用颜色进行了区分，药物敏感性使用了空实圈进行标注，黑色圆圈表示耐受对应药物，白色圆圈表示对应药物敏感。
除此之外，也可定制其它高级分析内容，欢迎咨询！
分枝杆菌基因组提取样品制备方法
一．样品培养
1.固体培养
吸取一定体积的菌液（根据菌体浓度调整吸取体积）至固体培养基划线（单线法、四区法均可），锡箔纸包裹，37°C倒置培养，直至克隆直径为0.5-1 cm大小。
2.液体培养
挑取单克隆或吸取一定体积的菌液（根据菌体浓度调整吸取体积）至液体培养基中，培养基体积不超过容器总体积的1/5或特殊容器规定的最大培养体积。37°C静置培养，直至菌液浓度达到OD600=0.8-1。
二. 样品灭活
1.菌落处理方法
挑取黄豆粒大小的菌落至TE溶液（25 mM Tris-HCl (pH 8.0)，10 mM EDTA）中，尽量避免挑取任何固体培养基（混有培养基会影响后续基因组提取效果），TE溶液体积需淹没菌落。金属浴80℃加热30分钟。灭活后的样品立即取出并冻存于-80℃冰箱。寄送前请勿反复冻融。
2.菌液处理方法
吸取1mL浓度达到OD600=0.8-1新鲜菌液至螺纹管中，13000 rpm离心10分钟。弃掉上清，加入500 μL的TE溶液（25 mM Tris-HCl (pH 8.0)，10 mM EDTA）。金属浴80℃加热30分钟。灭活后的样品立即取出并冻存于-80℃冰箱。寄送前请勿反复冻融。
| 菌株分类 | 序号 | 菌株名称 | 灭活条件 |
| 细菌 | 1 | 沙门氏菌 | 80℃水浴15min |
| 2 | 志贺氏菌 | 70℃水浴15min | |
| 3 | 单核细胞增生李斯特氏菌 | 80℃水浴15min | |
| 4 | 副溶血性弧菌 | 60℃水浴30min | |
| 5 | 结核分枝杆菌 | 80℃水浴30min | |
| 6 | 唐菖蒲伯克霍尔德 | 80℃水浴30min | |
| 7 | 肺炎链球菌 | 60℃水浴30min | |
| 8 | 流感嗜血杆菌 | 60℃水浴30min | |
| 9 | 霍乱弧菌 | 60℃水浴30min | |
| 10 | 创伤弧菌 | 60℃水浴30min | |
| 11 | 非结核分枝杆菌 | 80℃水浴30min | |
| 12 | 大肠埃希氏菌 | 70℃水浴15min | |
| 13 | 鲍曼不动杆菌 | 80℃水浴30min | |
| 14 | 幽门螺杆菌 | 80℃水浴30min | |
| 15 | 肺炎克雷伯菌 | 60℃水浴30min | |
| 16 | 支原体 | 60℃水浴30min | |
| 17 | 铜绿假单胞菌 | 60℃水浴60min | |
| 18 | 沙门氏菌 | 80℃水浴15min | |
| 19 | 金黄色葡萄球菌 | 80℃水浴30min | |
| 20 | 化脓链球菌 | 60℃水浴30min | |
| 真菌 | 1 | 白假丝酵母菌(白色念珠菌) | 60℃水浴60min |
| 2 | 马尔尼菲篮状菌 | 60℃水浴60min | |
| 3 | 新生隐球菌 | 60℃水浴10min |
以上均为最佳灭活条件，如温度过高、时间过长，或没有加TE buffer，则可能影响提取效果。
请勿将干菌体直接进行灭活，根据晶诺实验室的操作标准流程，干灭的样本会被认为仍然存在风险，如仍需继续项目，我们将启动生物安全实验室流程并产生费用，请注意！
如需TE Buffer，我们有100 mL和500 mL规格的产品，请与我们联系。
产品名称：TE缓冲液
产品货号：60340K
三．样品运输
样品在运输途中避免温度的变化而出现DNA降解，推荐10-15kg干冰运输。
【温馨小提示】
*冻存盒可用皮筋或胶带缠好，如盒内有空隙可用干净的纸填充，避免运输过程中样本晃动使盖子打开。
*如用自封袋，请保证袋子不易破损和袋口密封，避免样本掉落。可再套一两个袋子。
*请勿直接将样本管置于干冰中。
测序样本命名标准
1、样本编号只能由阿拉伯数字、字母、- 组成。例如A001-1是可以识别的编号。
2、不要用数字0作为开头。
3、样本编号不能包含中文、空格、/\ . : * ? ” < > | ① ②……等特殊字符，如有，分析时我们会用拼音替换中文，数字替换特殊字符。
分享一些有关结核分枝杆菌全基因组测序分方面的综述文章，希望能对结核基因组测序感兴趣的您有所帮助。
文献分享：
1、张洁,任怡宣,潘丽萍,张宗德.全基因组测序在结核分枝杆菌研究中的应用[J].中国防痨杂志，2020年7月第42卷第7期.
链接：http://www.zgflzz.cn/CN/10.3969/j.issn.1000-6621.2020.07.017
2、徐鹏,甘明宇,高谦.二代测序技术在结核分枝杆菌研究中的应用进展[J].微生物与感染，2015年2月25日，10（1）：54-60.
链接：http://jmi.fudan.edu.cn/CN/abstract/abstract491.shtml
