对棉花降解的微生物群落基因组学进行生物勘探：新型糖苷水解酶的挖掘

1、文章信息文章历史：2016年3月30日发送2016年7月20日修订2016年7月22日接受2016年7月25日网上发布摘要关键词：糖昔水解酶微生物群落宏基因组学天然棉花表达对棉花降解的微生物群落基因组学进行生物勘探：新型糖昔水解酶的挖掘GuoxiuZhang,PeiLiu,LeiZhang,WeiWei,XuedongWang,DongzhiWei,WeiWang*生物反应器国家重点实验室，新世界生物技术研究所，华东理工大学，上海200237,中国高效表达的糖昔水解酶(GHaseS)急需用于工业生产过程中将植物生物量有效降解成可发酵性糖。目前的研究表明涉及粗棉生物降解的没谱的功能特征。细菌与底

2、物的物理接触对于纤维素的高效水解是必须的。噬纤维菌目通过16SrRNA分析鉴定为普遍存在的菌落，在棉花生物降解中扮演重要的角色。从宏基因组数据中确定了2058中糖昔水解酶同工酶，其中16个在大肠杆菌中成功表达。其中四种酶对对硝基苯基-B-D-吐喃木糖昔有活性，四种对对硝基苯基-B-D-吐喃葡萄糖昔有活性，两种具有对硝基苯基-B-D-葡糖甘酸抗性，一种对昆布素有活性，三种具有对硝基茉基-B-D-氨基葡糖昔抗性，一种对竣甲基纤维素有活性，一种对对硝基茉基-B-D-甘露糖昔有活性。宏基因组为新型生物降解酶的挖掘提供良好的资源。被克隆的16种糖昔水解酶在生物转化和生物产品的生产方面有潜在的应用。棉花生

3、物降解的酶谱的功能特征和糖昔水解酶的分析为酶的合成与分解提供见解及植物生物降解的途径。1 .引言组成高等植物细胞壁主要成分的纤维素和半纤维素是地球上最丰富的可再生碳源(Kobayashietal.,2016)。自然环境中，植物体被微生物菌落产生的纤维素酶降解。糖背水解酶是分布广泛的能够水解糖昔键的酶，主要参与植物体的降解(DaviesandHenrissat,1995)。纤维素是水不溶性的线性聚合物，由重复的3-D-口比喃葡萄糖单元组成。纤维素的微生物降解涉及纤维素水解酶之间的复杂协同作用。普遍知道有三种类型的纤维素酶，包括内切葡聚糖酶，外切葡聚糖酶和3-葡糖甘酶协同涉及纤维素的完全水解(Ly

4、ndetal.,2002)。植物细胞壁中含量仅次于纤维素的半纤维素，是包含木聚糖、阿拉伯木聚糖、甘露聚糖、半乳甘露聚糖、葡糖甘露聚糖等的异质体(Girioetal.,2010)。半纤维素酶如木聚糖酶、甘露聚糖酶、3-木糖聚美、3-甘露聚糖酶、a-阿拉伯吠喃糖甘酶和a-半乳糖甘酶，也在纤维素解聚过程中扮演着水解去除半纤维素成分的角色。在过去十年的纸浆制造，食品和饲料加工以及生物燃料生产中，对纤维素酶和半纤维素酶的应用越来越感兴趣(JuturuandWu,2014)。纤维素酶和木聚糖酶已被用作纤维性能的生物修饰，从而改善纸浆制造的排水和透气性(Dienesetal.,2004)。苹果酸酶(纤维素酶

5、、木聚糖酶、果胶酶)已被用作果汁蔬菜汁的提取和澄清(deCarvalhoetal.,2008)。当前，对使用含高水平纤维素酶和半纤维素酶的酶来改善青贮饲料的利用、谷物脱壳和更好的饲料乳化是极其有意义的(JuturuandWu,2014)。使用纤维素酶来水解木质纤维材料中的纤维素是很吸引人的，因为酶解聚作用是生态友好的。将木质纤维素分解为简单糖类吸引了越来越多的生物能源工业。然而，廉价纤维素酶和半纤维素酶的大规模生产仍是一大瓶颈。为了满足工业过程中纤维素酶和半纤维素酶的新兴需求，立即需要具有更好性能的微生物纤维素酶和半纤维素酶来降低工业生产成本(Himmeletal.,2007;Juturuan

6、dWu,2014,2013,2014)。在自然界中，植物体被复杂的微生物过程分解。植物细胞壁-降解装置在好氧和厌氧微生物体有明显区别(ShallomandShoham,2003)。已经报道了以游离酶或多蛋白复合物的形式分泌这些装置的微生物。有趣的是，某些好氧滑动细菌有很强的纤维素降解能力，如噬纤维菌、生抱嗜纤维菌表现独特的纤维素降解活性，它能在低营养条件下迅速水解结晶纤维素，但是胞外纤维素降解活性十分低(Zhuetal.,2010)。直到现在我们对于基因的组分、协同相互作用以及微生物菌落的植物体降解途径的认识仍是肤浅的。我们需要更深入的研究去更好的理解酶解、微生物解以及植物细胞壁相关的发酵过程

7、。一.30.据估计，大约有4-6X10原核生物栖息在地球上，每克土壤约有4000种(Curtisetal.,2002;CurtisandSloan,2004)。虽然土壤微生物群落非常复杂，但它们被证明是探索新型植物生物降解生物体和酶类的有吸引力的来源。然而，大约95%-99.9%的微生物还未用标准实验室技术培养出来(Amannetal.,1990)。作为绕过基于培养方法的限制，宏基因组学正在成为直接研究微生物群落物种多样性、探索新型生物催化剂、从环境样品中分析其生物化学活性的有力途径(LorenzandEck,2005;Lietal.,2009)。为了探索新型生物降解酶，IlluminaHis

8、eq2000(IlluminaInc.,SanDiego,USA)常常调查土壤微生物群体降解天然棉花的能力，尤其着力于研究纤维素酶和半纤维素酶。为获得优势社区成员分布的代表性估计，使用IlluminaMiSeq平台进行16SrRNAU序。CAZy(/)数据库进行BLASTX佥索，并鉴定了含有2058个候选基因的优势糖背水解酶(GH家族，并进一步研究了基因克隆、蛋白表达和粗酶表征。对土壤微生物群体的测序来挖掘新型生物降解酶会提高生物量转化和生物产物生产的技术。2 .方法2.1 .土壤采样及加工从中国山东临沂棉花地收集到的土壤样品储存在冰箱中运送至实验室，大约1

9、00克亚样品储存在-80C。为了增加生物量降解物的丰度并加快发现新的糖甘水解酶基因的过程，将样品悬浮液的稀释液在含有粗棉作为主要碳源的Mandels培养基(pH5.5)中的振荡烧瓶中在30c下培养(MandelsandAndreotti,1987)。悬浮3夜每7天传代培养一次，不含微生物悬浮液的培养基作为对照。从棉田传播土壤样品42次后，使用数码相机进行棉花分解的数字成像。微生物悬浮液和上清液的糖甘水解酶活性可通过p-硝基苯基底物或者多糖底物测定。2.2 基因组DNA分离通过富集培养的发酵液体培养基在10000g、4c下离心10min。SedimentswerecollectedforDNAe

10、xtractionusingtheFastDNASPINKitforsoil(MPBiomed-icals,Solon,OH,USA%用0.8%的凝胶电泳观察单体DNA提取物，DNA的浓度和提取物的纯度是由显微分光光度法所确定的(NanoDrop-1000,ThermoScientific,Willmington,DE,USA)。最后，已经质检的DN群品用于霰弹枪配对端库构建和随后的高通量测序。2.3 MiSeq焦磷酸测序与16S数据分析从基因组DNA中扩增得到的细菌16SrRNA基因的V4不变区用通用引物Y1和丫2对真菌来说，引物ITS休DITS2ITS1进行扩增。PCRi程依据于先前的报道

11、(Zhaoetal.,2014)。成功扩增后，PCR物在Personalbio(上海，中国)使用IlluminaMiSeq(USA)进行末端焦磷酸序列测序。使用分割库中的默认参数对序列进行质量控制(Caporasoetal.,2010)。在除去条形码和引物序列后，剩下的序列以97%序列同一性的级别聚类成操作分类单位(OTU)。使用RD明类器对OTU代表序列进行分类(Wangetal.,2007)(补充材料S1所得到的OTUI用于确定分类学相对丰度和随后的“-多样性的统计分析2.4HiSeq焦磷酸测序宏基因组测序是在中国上海的Personalbio用IlluminaHiseq2000完成的(Il

12、luminaInc.,SanDiego,CA,USA已经质检的DNA任意断裂成300-400bp的小片段。然后将含有标记的接头的小片段构建到DNA文库中，然后直接进行配对末端测序。测序策略是PE101+8+101周期的指标(配对末端测序，8bp指标序列对应可读101bp)。大约可从样品中产生3.3Gb的宏基因组数据。原始数据由FastQC®行初步评估(用于高通量序列数据的质量控制工具)。2.5 宏基因组分析根据原始数据的质量，通过DynamicTrim.pl过滤序列，过滤低质量碱基，然后通过质量过滤筛选序列(Coxetal.,2010)。然后，除去小于75bp的序列，并将天鹅绒1.2

13、.10用于拼接序列以获得的最终重叠群(ZerbinoandBirney,2008)。最终组装的重叠群体已提交给MG-RAST(/,accessionnumber:4557667.3)来完成基因功能注释(与GenBank比对)和注释分类(由NOG和SEEDS系统注释)(最大值是10-5,最小值是70%)。将糖昔水解酶序列的重叠群进一步提取并上载到MG-RAST(accessionnumber:4560530.3)。最后，糖背水解酶序列的全长由ORFfinder工具(/projects/gorf

14、/)、BLASTX(http:/Blast.cgi)和NNPP(http:/seqtools/promoter.html)获得(补充材料S2)。2.6克隆并表达全长糖昔水解酶基因总共24个全长候选基因使用引物从宏基因组DNA进行PCFT增，并根据每个候选序列进行设计。扩增片段随后用pEASY-Uni无缝克隆和装配套件(中国北京感受态细胞)根据制造商的协议克隆到表达载体pET22b中。通过DNA测序确认每个质粒，并将其导入大肠杆菌宿主BL21(DE3)或Transetta(DE3)(Trans-Gen,Beijing,C

15、hina)用于蛋白质表达，大肠杆菌重组子在有50g/mL的氨茉青霉素的T断养基中37c下培养直到OD600nm=0.5-0.6,然后培养物用异丙基-3-D-硫代半乳糖昔(IPTG0.05M最终浓度)18c下诱导8h。细菌收获后，重悬于1/6培养基体积的裂解平衡洗涤(LEW)冲液(50mMsodiumphosphate,pH7.0,300mMNaCl,10mM2-mercaptoethanol,10%TritonX-100),中，然后用超声波破胞。细胞裂解物离心并且用SDS-PAG检测上清液与沉淀物。2.7 酶活性检测诱导细胞在100mMNaCl的50mM磷酸缓冲?夜(pH4topH10)中收获

16、和重悬。超声波破碎后，在37c下，在40科l含有100mMNaCl的50mM磷酸盐缓冲液(pH4至pH10)和50科l0.5%特异性底物(一式三份)中测试106细胞裂解物。培养合适的一段时间后，加入100dL的二硝基水杨酸并在99c下加热10min是反应终止(Miller,1959)。水解产物测540nm处的吸光值。用pNP进行的酶测试如下进行一式三份：10dL的细胞裂解物在10科的5mMpNP底物和30dL的50mM磷酸缓冲液(pH413H10)中培养。加入100dL的1M碳酸钠停止试验。水解产物对硝基苯基在405nm处测吸光值。候选蛋白的最适pH在最适pH范围中确定。纯葡萄糖、木糖、半乳糖

17、醛酸和pNP用于制作标准曲线。底物从SigmaAldrich公司(St,Louis,MO,USA购买得到。3.结果与讨论3.1 天然棉花降解和酶活性检测为了寻找相对不足的基因，富集可以增加其克隆的可能性，并加快发现新基因的过程。在原始土壤样品悬浮液的作用下，棉花完全降解需要2个月。从棉田传播土壤样品42次后，粗棉的降解速度越来越快。7天内，絮状棉完全降解，液体培养基颜色从光线变为棕褐色(Figure1)。对于控制变量，絮状棉花没有显着变化。此外，为了分析腐解过程，用对硝基苯Figr1.Irnj-singofEDn1td；left:rindderddedcollovibinmus:right.T

18、h干microbijilsuspensian.was,&utxuhuiedfor42面仃7djysjc30CinshakingflanksinMjndelsmedium.(pM5.5)withcrudrcollaii邦mjjr>rcjrhansimjic?,the-mediiumwithoutmicrobiillsuspensionasrpn(Tflia基底物和多糖底物来检测微生物裂解物和沉淀物的糖甘水解酶的活性。即使活性很低，微生物裂解物仍能水解pNP-口比喃葡糖、pNP-口比喃木聚糖、pNP氨基葡糖酸、pNP-口比喃甘露糖甘、pNP名千维二糖糖昔和pNP-葡糖甘酸(Fig.2

19、A)。酶在广泛的pH(frompH5topH8).具有高活性。水解pNP-葡糖甘酸的最适pH在8;而对于pNP派生物最适pH是6.微生物裂解物对竣甲基纤维素和昆布素、a,3-1-3:1-6葡萄糖连接的多糖也展现活性(Fig.2B)。然而，活性很低。类似的是，CMCNa仅被微生物裂解物水解并且对昆布素无反应。3.2 土壤微生物群落组成分析去除条形码和引物序列后，产生578830.12o.wo.sO.G0.4a,2iua音理名卜密oA质科强i次*一峭闺呻g川而出”“in勺pyrin曲期也bioshie*pNP宙iMCiLroiiia0.0.rlS6pHF*工Evahiatlond("雨j

20、ciavltyof匕mkr-nbljl.siisprnElonagansl-p-nicropbrrylwgirjiejmd|»5J"lurideimdHrHesMdbnernipllA：在时产时机正0安匕y&using,斯力irmMimylde*MhA电工a：Eiuyrnitic占石2u51nM划uc<h±-nl-JiiningpalyiarctundH.DatadrjwnfrombioAnmlrvpli口tr，：ri-3;Errortun可阅读序列。典型的序列在分类学上用RDP分析，鉴定了97%序列相似的214操作分类单元。！16S序列分析表明土壤

21、微生物群体的多样性。阿尔法多样性由赵，ACE辛普森，香农和Coverage五个指标估计（补充材料S3）,99.89%被完全覆盖。正如我们在宏基因组分析中发现的，Proteobacteria和Bacteroidetes代表了最丰富的微生物门（分别为72%和17%）,其次是Firmicutes（10%）（Fig.3）。Proteobacteria占全社区的72%,由Y-变形杆菌和3-变形杆菌占主导地位（分另1J为41%和26.1%）。Gamma-Proteobacteria和3-Proteobacteria主要由Burkholderiales（26%）,Pseudomonadales（19.7%）

22、和Xanthomonadales（14.2%）的成员组成。伯氏假单胞菌，假单胞菌属和黄单胞菌属主要以大麦科植物（22.3%）,不动杆菌（11）和嗜麦芽窄食单胞菌（11.4%）为主。细菌类主要由细胞吞噬体（14.3%）主导，主要由Dy属（Dyadobacter）（14.2%）表示。占总社区10%的Firmicutes以芽抱杆菌属为主（9.99%）。最后，还观察到一群仅占社区0.4%的真核生物，表明真菌在粗棉分解中不起主要作用。3.3 序列分析通常来说，使用其间隔为10bp的酶和该反应相容的缓冲液，否则双酶切应该顺序表达。当前研究表明，OM-6和OK-5碱性蛋白酶基因的扩增产物和pET21a用Ba

23、mH1和SalI双酶切（Fig.2）3.2.3限制性分析的克隆确认阳性菌株的确认基于质粒载体的氨茉青霉素抗性的参照特性。而且，含有蛋白酶基因的阳性克隆子的选择基于插入物的释放，并且我们观察到在琼脂糖凝胶上OM-6和OK-5的位置分别显示为大约1.2kb和0.5kb处。在克隆的每一步过程中基于生物信息学的工具CLCWorkbench被用来预测结果以及规范化方案。在克隆到大肠杆菌BL21之前，选择大肠杆菌DH5a菌株用于转化，因为对于目标DNA的起始克隆它是一个方便的宿主，因其具有高转化效率以及高质量DNA的产生。DH5a不包含T7RN咪合酶的基因；这就使得在非表达条件下构建重组质粒显得不合理。当

24、我们使用包含T7乳糖操纵子的质粒来表达时，宿主菌株中不包含用来目的基因表达的pLysS宿主菌株蛋白酶不足对于我们近期研究时满意的，以保证菌株不仅促进功能性产物的形成，而且确保没有假阳性结果的几率。3.3重组蛋白酶的表达将组装的重叠群体提交给MG-RAST（/）。由于人造重复读取（ADR）未被检测到，只有一个序列未通过QC水线。共406,305个序列（包括190,880,944个基数对平均长度为469bps的对象）通过质量控制。约98.3%（399,443）序列共产生461,990个预测蛋白质编码区，其中321,789个（69.7%的特征）被使

25、用至少一个蛋白质数据库（M5NR）,其余的140,201（30.3%的特征）与蛋白质数据库（孤儿）没有显着的相似性。其余251,195个功能（321,789个注释功能的78.1%）被分配到功能类别。分析结果显示在MG-RAST0站上（登录号：4557667.3）。如图1所示。S1（补充材料S4）,功能类别的百分比通过采集组装的重叠群构建。结果表明，最丰富的类别是基于聚类的子系统，其中包括没有特定/已知任务的功能耦合证据的基因。高度丰富的类别遵循碳水化合物，氨基酸，蛋白质和DNA/RNA代谢的预期顺序，也在MG-RAS近艮务器中的其他可比较的基因组中也类似地观察到。另一方面，几种蛋白质与光合作用

26、，次生代谢，休眠和抱子形成有关。致谢S.D.Gohe酢常感谢印度新德里UGC的一等奖学金，这项工作是由印度拉杰果德Saurashtra大学赞助。参考文献1 M.S.Dodia,H.G,Bhimani,C.M,Rawal,R.H,Joshi,S.P,Singh,BioresourceTechnol.99(2008)6223-6227.2 M.S.Dodia,C.M.Rawal,H.G.Bhimani,R.H.Joshi,S.K.Khare,S.P.Singh,J.Ind.Microbiol.Biotechnol.35(2008)121-131.S.D.Gohel,S.P.Singh/Intern

27、ationalJournalofBiologicalMacromolecules50(2012)664-6716713 A.Gupta,I.Roy,R.K.Patel,S.P.Singh,S.K.Khare,M.N.Gupta,J.Chromatogr.A1075(2005)103-108.4 U.Boominadhan,R.Rajakumar,P.K.V.Sivakumaar,M.Melvin,Bot.Res.Int.2(2)(2009)8387.5 S.Ramesh,M.Rajesh,N.Mathivanan,BioprocessBiosyst.Eng.32(2009)91800.6 L.

28、V.A.Reddy,J.Chem.Technol.Biotechnol.8(11)(2008)15261533.7 P.C?alik,N.Kalender,H.Tunc?er,EnzymeMicrob.Tech.32(2003)706720.8 J.T.Thumar,S.P.Singh,Braz.J.Microbiol.38(2007)19.9 M.J.Page,S.L.Wong,J.Heiqwtt,N.C.Strynadka,R.T.MacGillivray,Biochemistry42(30)(2003)90609066.10 P.L.Jorgensen,M.Tangney,P.E.Ped

29、ersen,S.Hastrup,B.Diderichsen,S.T.Jorgensen,Appl.Environ.Microbiol.66(2)(2000)825827.11 H.Tsujibo,T.Kubota,M.Yamamoto,K.Miyamoto,Y.Inamori,Appl.Environ.Microbiol.69(2003)894900.12 J.M.Fernandez-Abalos,V.Reviejo,M.Diaz,S.Rodriguez,F.Leal,R.L.Santamaia,Microbiology149(7)(2003)162232.13 S.Bakhtiar,J.Ve

30、vodova,R.Hatti-Kaul,X.D.Su,ActaCrystallgr.DBiol.Crystallogr.59(3)(2002)529531.14 P.M.Carolina,G.Augusto,D.Castro-Ochoa,A.Farres,Appl.Microbiol.Biotechnol.78(2008)603-612.15 B.Yan,X.Chen,X.Hou,H.He,B.Zhou,Y.Zhang,Extremophiles13(2009)725733.16 X.Ni,L.Yue,Z.Chi,Z.Li,X.Wang,C.Madzak,Mar.Biotechnol.11(2

31、009)8189.17 M.K.Purohit,S.P.Singh,Int.J.Biol.Macromol.49(2011)103112.18 M.K.Purohit,S.P.Singh,Lett.Appl.Microbiol.49(3)(2009)338344.19 B.Hagihara,TheEnzymes,vol.4,AcademicPress,NewYork,1958,p.224.20 J.Sambrook,E.F.Fritsch,T.Maniatis,MolecularCloning.ALaboratoryManual,ColdSpringHarborLaboratoryPress,

32、ColdSpring,Harbor,NewYork,USA,1989.21 J.Daintith,Workbench,ADictionaryofComputing,2004.22 A.Balan,C.S.deSouza,A.Moutran,R.C.Ferreira,C.S.Franco,C.H.Ramos,L.C.deSouzaFerreira,ProteinExpr.Purif.43(2005)103110.23 L.Liu,R.Wang,W.He1,F.He,G.Huang,ActaBiochim.Biophys.42(2010)585592.24 F.Rajamohan,C.R.Engs

33、trom,T.J.Denton,L.A.Engen,I.Kourinov,F.M.Uckun,ProteinExpr.Purif.16(1999)359368.25 S.Y.Kwon,C.S.An,J.R.Liu,S.S.Kwak,H.S.Lee,J.K.Kim,K.H.Paek,Mol.Cells10(2000)812.26 E.Pittaluga,A.Poma,A.Tucci,L.Spano,J.Biotechnol.117(2005)263266.27 M.J.Luo,W.X.Liu,X.Y.Yang,Y.Xu,F.Yan,P.Huang,F.Chen,Russ.J.PlantPhysiol.54(2007)202206.28 S.G