MPB：中科院生态环境中心邓晔等-细胞内融合基因技术 (epicPCR) 测定功能类群多样性

为进一步提高《微生物组实验手册》稿件质量，本项目新增大众评审环节。文章在通过同行评审后，采用公众号推送方式分享全文，任何人均可在线提交修改意见。公众号格式显示略有问题，建议电脑端点击文末阅读原文下载PDF审稿。在线文档(https://kdocs.cn/l/cL8RRqHIL)大众评审页面登记姓名、单位和行号索引的修改建议。修改意见的征集截止时间为推文发布后的72小时，文章将会结合有建设性的修改意见进一步修改后获得DOI在线发表，同时根据贡献程度列为审稿人或致谢。感谢广大同行提出宝贵意见。细胞内融合基因技术 (epicPCR) 测定功能类群多样性Diversity of Functional Microbes Detected by epicPCR Technology沈文丽1，王尚2 *，邓晔2, 31山东大学海洋研究院，山东大学，青岛，266237；2中国科学院环境生物技术重点实验室，中国科学院生态环境研究中心，北京，100085；3中国科学院大学资源与环境学院，中国科学院大学，北京，100049*通讯作者邮箱：[e-mail]摘要：细胞内融合基因技术是一项结合了单细胞分离、融合PCR、巢式PCR和测序技术的新型分子技术，不仅可以在未培养的单细胞中将功能基因与系统发育标记基因 (如16S rRNA基因) 连接在一起，而且通量高，成本低。该技术首先通过稀释和涡旋将样品分散成单细胞体系，然后对拟研究的系统发育基因和功能基因进行融合扩增，随后进行巢式PCR，建库并进行高通量测序。细胞内融合基因技术既能够解决扩增子测序物种信息与功能信息脱节的困境，又能够克服大规模宏基因组测序成本高、分析难且研究对象仅是优势物种的局限。关键词：细胞内融合基因技术，功能基因，融合PCR，巢式PCR材料与试剂阶段一：生成聚丙烯酰胺凝珠1.1.5 ml、2 ml、50 ml离心管2.PCR小管3.35 μm细胞筛4.0.2 μm滤膜5.各种型号移液枪头6.双蒸水7.过硫酸铵 (St. Louis, MO, USA, Sigma, catalog number: A3678)8.N,N'-双 (丙烯酰) 胱胺BAC (St. Louis, MO, USA, Sigma)9.丙烯酰胺 (St. Louis, MO, USA, Sigma)10.Tris-HCl (pH 7.5)11.氯化钾12.矿物油 (St. Louis, MO, USA, Sigma, catalog number: M8410)13.司盘80 (St. Louis, MO, USA, Sigma, catalog number: 85548)14.吐温80 (St. Louis, MO, USA, Sigma, catalog number: P8192)15.四甲基乙二胺 (TEMED) (St. Louis, MO, USA, Sigma)16.二乙醚 (St. Louis, MO, USA, Sigma)17.蛋白酶K (St. Louis, MO, USA, Sigma)18.溶菌酶 (St. Louis, MO, USA, Sigma)阶段二：融合PCR1.2 mm玻璃珠2.2 ml小管3.各种型号移液枪头4.PCR小管5.ABIL EM 906.Triton X-1007.矿物油 (St. Louis, MO, USA, Sigma, catalog number: M8410)8.双蒸水9.5× Phusion HF 缓冲液10.超保真酶Phusion Hot Start Flex DNA Polymerase (NEB)11.50 mM MgCl212.10 mM dNTPs13.正向引物F114.反向引物R215.融合引物F2-R116.牛血清白蛋白BSA (molecular biology grade, NEB, Ipswich, MA, USA)17.吐温20 (St. Louis, MO, USA, Sigma)18.乙二胺四乙酸EDTA (St. Louis, MO, USA, Sigma)阶段三：破乳化释放聚丙烯酰胺凝珠1.1.5 ml、2 ml、50 ml离心管2.各种型号移液枪头3.双蒸水4.二乙醚 (St. Louis, MO, USA, Sigma)5.乙酸乙酯 (St. Louis, MO, USA, Sigma)6.AMPure XP磁珠 (Beckman Coulter, Danvers, MA, USA)7.无水乙醇阶段四：巢式PCR1.1.5 ml、2 ml、50 ml离心管2.PCR小管3.各种型号移液枪头4.双蒸水5.超保真酶Phusion Hot Start Flex DNA Polymerase (NEB)6.5× Phusion HF 缓冲液7.50 mM MgCl28.10 mM dNTPs9.正向引物F310.反向引物R311.阻断引物blockF12.阻断引物blockR13.SYBR Green I 核酸染色试剂盒 (10,000 X, Invitrogen, Waltham, MA, USA)表1. 本项目中融合硫酸盐还原菌的功能基因dsrB和16S rRNA所使用的引物信息引物名称引物序列 (5’-3’)F1 (dsrB-F1)GTGTAGCAGTTACCGCAR1-F2 (dsrB-R1_519R)GWATTACCGCGGCKGCTGTGCCTSAAYATGTGYGGYGR2 (1492)GGTTACCTTGTTACGACTTF3 (dsrB-F3)VAGVATSGCGATRTCGGAR3 (E786R)GGACTACHVGGGTWTCTAATBlockR (U519R-block10)TTTTTTTTTTGWATTACCGCGGCKGCTG/3SpC3/BlockF (U519R-block10)TTTTTTTTTTCAGCMGCCGCGGTAATWC/3SpC3/溶液配方1.丙烯酰胺溶液12% 丙烯酰胺0.32% BAC2. 1×TK 缓冲液20 mM Tris HCl60 mM KCl3. STT 乳化油 (14天内可用，每次使用前需要颠倒混匀)4.5% 司盘800.4%吐温800.05%Triton X-100v/v溶于矿物油中4. ABIL 乳化油 (常温储存)4% ABIL EM900.05%Triton X-100v/v溶于矿物油中仪器设备阶段一：生成聚丙烯酰胺凝珠1.1.5 ml台式高速离心机2.各种型号移液枪3.涡旋仪4.超声波清洗仪5.超净操作台6.通风橱阶段二：融合PCR1.超净操作台2.涡旋仪3.PCR仪4.PCR管离心仪5.1.5 ml台式高速离心机阶段三：破乳化释放聚丙烯酰胺凝珠1.1.5 ml台式高速离心机2.各种型号移液枪3.通风橱4.磁珠清洗磁力架5.超净工作台阶段四：破乳化释放聚丙烯酰胺凝珠1.PCR仪2.qPCR仪3.超净操作台4.各种型号移液枪5.PCR管离心仪实验步骤一、生成聚丙烯酰胺凝珠1.细胞悬浮液准备通过ATP等方法对样品细胞悬浮液的浓度进行定量检测，稀释得到30 μl包含1,400万个细胞的样品悬浮液。2.制备聚丙烯酰胺凝珠2.1混合30 μl细胞悬浮液 (1,400万个细胞) 、200 μl丙烯酰胺溶液 (12%丙烯酰胺和0.32% BAC) 和25 μl过硫酸铵溶液 (10%) ，轻轻涡旋混匀。2.2在2 ml圆底离心管中，加入STT乳化油600 μl，然后加入2.1混合所得的水性悬液 (体积255 µl) ，3,000 rpm涡旋30 s，产生约5亿个液滴 (按照液滴平均直径10 μm计算) 。注：STT乳化油由司盘80、吐温80、聚乙二醇单辛基苯基醚和矿物油配置完成后14天内可用，每次使用前需要颠倒混匀。2.3加入四甲基乙二胺 (TEMED) 25 μl，催化聚合反应，3,000 rpm涡旋30 s。2.4乳化液静置90 min聚合，生成如图1所示的聚丙烯酰胺凝珠。3.二乙醚提取聚丙烯酰胺凝珠3.1在2.4乳化液中加入水饱和二乙醚800 μl，立即轻弹颠倒混匀，形成可见沉淀；用移液枪移走乳化油和乙醚的混合物，添加1 ml超纯水，颠倒离心管混匀。3.2将所有样品转移至新离心管中，离心12,000 x g 30 s，形成上层油层、中间水油混合层和下层聚丙烯酰胺bead。3.3采用移液枪移取上层油层后，加入超纯水洗涤，颠倒混匀，继续用移液枪移取油层；重复加入超纯水洗涤至不再有油层形成 (约5次) ；在最后一次水洗不再形成油层时，采用移液枪移取上层水。3.4注入1 ml 1× TK缓冲液，重悬浮珠子于缓冲液中；用35 µm细胞过滤器过滤，转移到1.5 ml离心管中，过滤完成的聚丙烯酰胺凝珠保存于4 °C。

图1. (a) 下层白色胶质物为生成的包裹细胞的聚丙烯酰胺凝珠；(b) 荧光显微镜下聚丙烯酰胺凝珠及其中包裹的单细胞微生物4.细胞溶解 (可选)添加0.8%的溶菌酶 (35,000 U/µl) ，于37 °C过夜培养。离心，移除上清；重悬浮在1× TK缓冲液中。用20%蛋白酶K (1 mg/ml) 和0.8%TritonX-100处理，37 °C培养30 min，然后95 °C培养10 min，重悬浮在1× TK缓冲液中。二、融合PCR1.配制PCR反应混合体系，各组分的体积为1.1倍样品数。试剂体积无菌水1 μL5×Phusion HF 缓冲液20 μL50 mM MgCl22 μL10 mM dNTPs2.5 μL正向引物F1 (10 μM)10 μL反向引物R2 (10 μM)10 μL融合引物R1-F2 (1 μM)1 μLBSA0.5 μL吐温-200.2 μLPhusion Hot Start Flex8μL2.每个样各分装55.2 μl混合液到2 ml 离心管中，加入45 μl合成的聚丙烯酰胺凝珠，混合均匀。加入900 μl ABIL油，3,000 rpm涡旋1 min，充分振荡混匀。3.分装60 μl反应混合液到PCR小管中，每个样约分装16管。然后进行PCR：94 °C 30 s, [94 °C 5 s, 52 °C 30 s, 72 °C 30 s] 33个循环, 72 °C 5 min, 4 °C ∞。4.反应结束后，立马将融合产物收集到1.5 ml离心管中 (将融合产物跑电泳，通常看不到任何产物条带，如图2所示)。在每个收集的样品中，加入终浓度为1 mM的EDTA，该样品可在4 °C保存过夜。融合产物跑水平电泳通常看不到任何条带，如图2所示。

图2. 融合PCR后产物跑胶图三、破乳化释放聚丙烯酰胺凝珠1.破坏ABIL乳化体系，释放凝珠1.1将阶段二生成的融合PCR产物在常温下13,000 x g 离心5 min，弃上层油相。1.2加入1 ml水饱和二乙醚至每个样品中， (同体积水-乙醚混合液，混匀过程中轻起瓶盖放气) ，轻轻涡旋混合，13,000 x g 离心1 min，弃上层液相，重复该步骤一次。1.3同样的方法，用水饱和乙酸乙酯清洗2次。1.4放在通风橱中吹10 min，使残留的二乙醚尽可能挥发干净，约收集100-150 μl产物。该产物可在4 °C保存数小时，在-20 °C保存过夜。2.磁珠清洗PCR产物2.1将磁珠充分摇匀，至于室温中平衡至室温 (10 min左右即可，最多30 min足够) ，每100 μl DNA样品中加入85.5 μl磁珠于1.5 ml离心管中，轻轻混匀，在室温中孵育13 min，使DNA充分结合到磁珠上。2.2将离心管置于磁铁上，分离磁珠2 min，移除悬液 (不要将离心管取下) 。2.3用0.5 ml 70%乙醇清洗磁珠2次 (不要将离心管取下) 。2.4离心管保留在磁铁上，打开离心管盖子，风干15-30 min，直至管壁上没有液滴，磁珠表面干燥。2.5将离心管从磁铁上取下，加入40 μl的缓冲液或无菌水，轻轻混匀，室温中孵育7 min。2.6将离心管置于磁铁上2 min，收集35-40 μl悬液，保存在1.5 ml离心管中。四、巢式PCR1.巢式qPCR (可选)1.1配制反应混合液，各组分的体积为1.1倍样品数试剂体积无菌水7.125 μL5×Phusion HF 缓冲液5 μL10 mM dNTPs0.5 μL正向引物F3 (10 μM)10 μL反向引物R3 (10 μM)10 μLBlockF (32 μM)2.5 μLBlockR (32 μM)2.5 μLPhusion Hot Start Flex0.25 μL100× SYBR Green I0.125 μL1.2每个样分装23 μl混合液至PCR小管中，加入2 μl纯化的融合PCR产物以及水 (作为阴性对照) 。1.3轻弹反应液使混匀，将PCR管置于PCR管离心仪上离心片刻。进行PCR过程：98 °C 30 s, [98 °C 5 s, 52 °C 30 s, 72 °C 30 s] 40个循环, 72 °C 5 min, 4 °C ∞。1.4根据qPCR Ct值评估不同样本的最少循环数，可通过稀释高浓度样本，使不同样本的循环数一样。该PCR产物可在4 °C保存数小时或者-20 °C过夜保存。2.巢式PCR2.1配制反应混合液，各组分的体积为4× 1.1倍样品数试剂体积5× Phusion HF 缓冲液5 μL10 mM dNTPs0.5 μL正向引物F3 (3 μM)2.5 μL反向引物R3 (3 μM)2.5 μLBlockF (32 μM)2.5 μLBlockR (32 μM)2.5 μLPhusion Hot Start Flex0.25 μL2.2每个样分装63 μl混合液至PCR小管中，加入37 μl纯化的融合PCR产物 (根据qPCR结果进行稀释) ，轻轻混匀，然后以25 μl体积进行分装。2.3进行PCR过程：PCR：98 °C 30 s, [98 °C 5 s, 52 °C 30 s, 72 °C 30 s] 40个循环 (循环数根据qPCR结果确定) , 72 °C 5 min, 4 °C ∞。3.将同一样品的平行样收集到1.5 ml离心管中，可以先用5 μl样品跑水平凝胶电泳，观察生成目标产物情况 (如图3所示)。若生成目标产物，将样品跑胶然后割胶纯化，再用磁珠清洗 (方法同阶段三步骤2) ，最后送测序。

图3. 巢式PCR产物跑胶图 (其中黄色箭头所指为产物条带)结果与分析【可选】本课题组采用细胞内融合基因技术 (epicPCR) 对采集自十个西藏盐湖底泥的硫酸盐还原菌群落的系统分类和多样性进行了鉴定，并对影响该群落结构的环境因子进行了鉴别。研究发现了十个新的硫酸盐还原菌门，并由此表明，环境中尚有大量的未知硫酸盐还原菌群类别。可参见如下结果：1.通过对融合的16S rRNA加dsrB基因的16S rRNA扩增子和epicPCR产物的测序，共检测到10个湖泊的微生物总群落的12,519个OTU和硫酸盐还原菌群的883个SRP OTUs (表5) 。在微生物总群落和硫酸盐还原菌类群中，变形杆菌、厚壁菌和拟杆菌是最主要的门，包括了最大比例的OTUs代表 (图4) 。整个微生物群落和SRP亚群落的α-多样性呈显著正相关 (R2=0.443，P<0.05) ，β-多样性部分重叠 (图5) ，揭示了两者间的关系。表5. 对西藏盐湖中微生物总群落和硫酸盐还原菌群落的门，属，以及OTUs水平的结构比较

图4. 微生物总群落和硫酸盐还原菌群落的OTUs的门的水平分类。外饼为微生物总群落，内饼为硫酸盐还原菌群落。

图5. 微生物总群落和硫酸盐还原菌的α-多样性 (a) 与β-多样性 (b) 的对比2.按照epicPCR最初创始文章的方法，筛选出了120个高丰度的SRP-OTU (至少一个样本中为1%) ，并对其构建了进化树以观察进化分化 (图6) 。这些高丰度的SRP-OTU与17个描述的门相关，其中只有7个广为认可的SRP门。类似于微生物总群落的大多数OTUs，大多数SRP-OTUs只存在于某个特定的湖泊中， 证明了在西藏盐湖中微生物群落和SRP亚群落的分布的地方特殊性。

图6. 120个高丰度的硫酸盐还原菌的OTUs的进化关系。加黑字体为参考序列，蓝色字体为广泛分布的OTUs (每一个湖中均有高丰度存在) ，黑色字体为只以高丰度存在于某一个湖中的OTUs，红色字体为高丰度存在于两个以上的湖，但并不在每个湖中都有的OTUs。进化枝上蓝色圆形大小表示自展值的大小 (70-100) 。红色进化枝代表属于新发现的十个硫酸盐还原菌门的OTUs (图下方 * 标记的菌门) 。黑色柱形代表该OTUs相对丰度最高的湖中的丰度。彩色条形标记门的水平的分类。致谢感谢自然科学基金重大研究计划培育项目 (项目号91851106) 以及国家自然科学基金青年基金项目 (项目号32001092) 的经费支持！感谢MIT研究团队Spencer等人，该实验方案参考自Spencer等人于2016年发表在ISME Journl上的文章以及他们提供的详细的实验流程，在此表示特别感谢！参考文献1.Spencer S J, Tamminen M V, Preheim S P, Guo M T, Briggs A W, Brito I L, D A W, Pitkanen L K, Vigneault F, Juhani Virta M P, Alm E J. Massively parallel sequencing of single cells by epicPCR links functional genes with phylogenetic markers. ISME J, 2016, 10: 427-436.2.Qin H, Wang S, Feng K, He Z, Virta M P J, Hou W, Dong H, Deng Y. Unraveling the diversity of sedimentary sulfate-reducing prokaryotes (SRP) across Tibetan saline lakes using epicPCR. Microbiome, 2019, 7: 71.