• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    絲狀真菌Podospora anserina AA11家族裂解多糖單加氧酶基因的鑒定和功能研究

    2023-12-20 09:16:18杜文珍李元敬吳佳玲陳思羽姜亮劉剛謝寧
    遺傳 2023年12期
    關(guān)鍵詞:子囊孢子突變型突變體

    杜文珍,李元敬,吳佳玲,陳思羽,姜亮,劉剛,謝寧

    研究報(bào)告

    絲狀真菌AA11家族裂解多糖單加氧酶基因的鑒定和功能研究

    杜文珍1,李元敬2,吳佳玲1,陳思羽1,姜亮1,劉剛1,謝寧1

    1. 深圳大學(xué)生命與海洋科學(xué)學(xué)院,深圳市微生物基因工程重點(diǎn)實(shí)驗(yàn)室,深圳 518060 2.華南師范大學(xué)材料與新能源學(xué)院,汕尾 516622

    輔助活性蛋白家族(auxiliary activity family,AA family)中的裂解多糖單加氧酶(lytic polysaccharide monooxygenase, LPMO)能催化纖維素、幾丁質(zhì)和淀粉等多種難降解碳水化合物的氧化解聚。盡管目前對(duì)LPMO的酶學(xué)研究較多,但對(duì)基因失活的研究卻鮮有報(bào)道。本研究利用同源重組方法定點(diǎn)敲除絲狀真菌中AA11家族的5個(gè)基因()、()、()、()和(),分別構(gòu)建了單突變體ΔPaLPMO11A (ΔA)、ΔPaLPMO11B (ΔB)、ΔPaLPMO11C (ΔC)、ΔPaLPMO11D (ΔD)和ΔPaLPMO11E (ΔE),然后通過(guò)遺傳雜交構(gòu)建所有多基因突變體。通過(guò)在不同碳源培養(yǎng)基上的表型分析、DAB和NBT染色以及纖維素酶活測(cè)定分析-野生型菌株與突變型菌株在生長(zhǎng)速率、有性生殖、氧化應(yīng)激和纖維素降解能力等方面的差異,揭示基因在菌株的生長(zhǎng)發(fā)育和木質(zhì)纖維素降解過(guò)程中的作用。實(shí)驗(yàn)結(jié)果表明,在不同纖維素碳源上,ΔBΔCΔE、ΔAΔBΔCΔE、ΔAΔCΔDΔE和ΔAΔBΔCΔDΔE突變型菌株的有性生殖能力降低,其余突變型菌株的孢子萌發(fā)效率、生長(zhǎng)速率和生殖能力幾乎沒(méi)有差異。PaLPMO11家族5個(gè)基因的同時(shí)缺失,會(huì)導(dǎo)致菌株利用各種碳源的能力明顯降低、生長(zhǎng)速率降低、孢子萌發(fā)率降低、子實(shí)體數(shù)減少、部分子實(shí)體發(fā)育異常、壽命縮短和降解纖維素的能力顯著下降,但仍有野生型45%以上的總纖維素酶活力。上述結(jié)果表明,基因可能參與的生長(zhǎng)發(fā)育、有性生殖、衰老和纖維素降解過(guò)程。本研究為系統(tǒng)闡述絲狀真菌中木質(zhì)纖維素降解的調(diào)控機(jī)制提供參考。

    絲狀真菌;;;基因敲除;木質(zhì)纖維素降解

    人類(lèi)使用的大多數(shù)燃料和化學(xué)品資源都來(lái)源于煤炭、石油和天然氣等不可再生的化石能源。隨著全球經(jīng)濟(jì)的快速發(fā)展和工業(yè)化規(guī)模的不斷擴(kuò)大,化石能源的過(guò)度消耗,引發(fā)了不可再生化石能源儲(chǔ)量枯竭、環(huán)境污染和全球變暖等一系列問(wèn)題,給全球環(huán)境保護(hù)和可持續(xù)發(fā)展帶來(lái)了嚴(yán)峻的挑戰(zhàn)。此外,不可再生化石燃料資源預(yù)計(jì)在未來(lái)40~50年內(nèi)耗盡[1],因此,勘探和開(kāi)發(fā)可持續(xù)、生態(tài)友好和儲(chǔ)量大的可再生生物質(zhì)資源已成為近年來(lái)的熱點(diǎn)問(wèn)題,有望克服因消耗化石燃料帶來(lái)的環(huán)境污染和能源短缺等問(wèn)題。--木質(zhì)纖維素生物質(zhì)是地球上最豐富的可再生資源之一,它的年產(chǎn)量約為2×1011噸[2],是獲得生物乙醇和替代化石燃料的理想可再生候選者之一。然而,迄今為止,大部分木質(zhì)纖維素資源尚未得到有效地利用,木質(zhì)纖維素資源的利用率甚至未能達(dá)到其總量的5%[3]。木質(zhì)纖維素主要由木質(zhì)素、纖維素和半纖維素組成,其中纖維素含量占比通常最高[4]。纖維素是生物圈中含量最豐富的有機(jī)多糖,是獲取潔凈能源如基于乙醇的生物燃料和高附加值化學(xué)品的主要原料[5]。因而,提高纖維素的有效利用率是提高木質(zhì)纖維素資源利用率的重點(diǎn)。在自然界中,許多生物都能有效地降解利用木質(zhì)纖維素生物質(zhì),其中能分泌碳水化合物活性酶的絲狀真菌被認(rèn)為是關(guān)鍵的初級(jí)降解者[6]。

    裂解多糖單加氧酶(lytic polysaccharide monoo-xy-genase,LPMO)是銅依賴(lài)的氧化還原酶,該酶通過(guò)氧化裂解頑固多糖鏈表面,使多糖易于處理以便進(jìn)一步酶促作用并最終降解[7]。LPMO廣泛存在于真菌、細(xì)菌、病毒和昆蟲(chóng)中,能催化纖維素、幾丁質(zhì)和淀粉等多種難降解頑固碳水化合物的氧化解聚。其在碳水化合物活性酶數(shù)據(jù)庫(kù)(carbohydrate-active enzymes database,CAZy) (http://www.cazy.org/)中被劃分到 AA9~AA11和AA13~AA16這7個(gè)不同的輔助活性家族(auxiliary activity family,AA family)中[8]。不同家族LPMO的氨基酸序列差異較大,但所有LPMO整體拓?fù)浣Y(jié)構(gòu)相似,催化結(jié)構(gòu)域是免疫球蛋白樣的反向平行β-三明治核心結(jié)構(gòu),活性位點(diǎn)位于該較為平坦的核心區(qū)域上,由兩個(gè)高度保守的組氨酸與金屬銅離子螯合形成的T型“組氨酸支架(histidine brace)”構(gòu)成,其中一個(gè)組氨酸是氨基末端殘基,另一個(gè)為側(cè)鏈組氨酸[9,10]。Hemsworth等[11]發(fā)現(xiàn)來(lái)自米曲霉()AA11家族的LPMO(AoAA11)對(duì)晶體幾丁質(zhì)具有活性,表明AoAA11具有與AA9和AA10的LPMO相似的三維結(jié)構(gòu),包括核心的反向平行β-三明治結(jié)構(gòu)和具有保守“組氨酸支架”的銅活性位點(diǎn)。Wang等[12]從富士鐮刀菌()中獲得截?cái)嘈问降腖PMO11新型酶,稱(chēng)為FfAA11。該酶對(duì)α-幾丁質(zhì)和β-幾丁質(zhì)以及龍蝦殼均表現(xiàn)出氧化活性。近年來(lái),Rieder等[13]發(fā)現(xiàn)煙曲霉()中的AA11B不具有假定的幾丁質(zhì)活性,而對(duì)N-乙酰葡糖胺的可溶性低聚物具有催化活性。St?pamo等[14]發(fā)現(xiàn)AfAA11A雖然缺少碳水化合物結(jié)合模塊,但對(duì)α-和β-幾丁質(zhì)仍具有相當(dāng)大的催化能力。結(jié)合以上研究可知,盡管LPMO11被假定具有幾丁質(zhì)酶活性,但AA11家族的部分成員對(duì)幾丁質(zhì)卻沒(méi)有催化活性。因此,LPMO11可能在真菌生理學(xué)中發(fā)揮著除幾丁質(zhì)酶催化活性外的其他作用[14]。

    迄今為止,雖然已經(jīng)有大量關(guān)于LPMO酶學(xué)的研究,但對(duì)于真菌中基因失活的研究卻鮮有報(bào)道。此外,尚未有利用缺失基因的突變體進(jìn)行復(fù)雜多糖如纖維素、淀粉或幾丁質(zhì)的降解研究,因此其功能尚未得到充分的研究。通過(guò)比對(duì)41種常見(jiàn)真菌的基因組發(fā)現(xiàn),是含有最多降解木質(zhì)纖維素蛋白酶類(lèi)基因的子囊真菌[15],又因其生命周期短,易于培養(yǎng)和進(jìn)行分子遺傳操作,因此它是用來(lái)研究木質(zhì)纖維素降解的良好模式生物[16]。盡管LPMOs在真菌界中廣泛存在,但屬于AA11家族的LPMO尚未得到充分研究。中5個(gè)基因功能仍有待進(jìn)一步被揭示與研究,本研究以期為闡明LPMO11在真菌生理學(xué)中發(fā)揮的不同作用提供一定的參考價(jià)值。

    本研究以為實(shí)驗(yàn)菌株,利用同源重組方法分別構(gòu)建5個(gè)基因和缺失單突變型菌株,通過(guò)遺傳雜交獲得所有多基因突變體。通過(guò)突變型菌株與野生型菌株在生長(zhǎng)速率、有性生殖、氧化應(yīng)激和纖維素降解等方面的差異性分析,揭示AA11家族的基因在中的明確作用與功能,為系統(tǒng)闡述絲狀真菌中木質(zhì)纖維素降解的調(diào)控機(jī)制提供參考。

    1 材料與方法

    1.1 材料

    (1) 菌株和質(zhì)粒:野生型菌株(WT)Δmus51::phleoR菌株(用腐草霉素抗性基因代替基因編碼序列所得到的菌株?;蚓哂行迯?fù)真菌DNA突變的功能,敲除可顯著提高同源重組效率[17])、Δmus51:: HygroR菌株、Δmus51:: NourR菌株和Δmus51:: GeneR菌株。質(zhì)粒pBC-Geneticin、pBC-Hygromycin、pBC- Phleomycin和pBC-Nourseothricin分別用于突變體構(gòu)建。

    (2) 培養(yǎng)基:菌株培養(yǎng)條件以及M2基礎(chǔ)培養(yǎng)基等常用培養(yǎng)基的配置方法可在Genome Project (http://podospora.i2bc.paris-saclay.fr/ methods.php)檢索得到。木質(zhì)纖維素培養(yǎng)基:在M2培養(yǎng)基的基礎(chǔ)上,分別用木質(zhì)素(lignin powder)、木屑(wood shaving)和微晶纖維素(avicel)等量取代糊精作為碳源,制備3種不同成分的木質(zhì)纖維素培養(yǎng)基。

    (3) 主要試劑和儀器:抗生素遺傳霉素(geneticin)和腐草霉素(phleomycin)購(gòu)自南京都萊生物技術(shù)有限公司;諾爾斯菌素(nourseothricin)購(gòu)自上海吉至生化科技有限公司;潮霉素(hygromycin)、羥甲基纖維素鈉(CMC-Na)和其他常用試劑均購(gòu)自上海生工生物工程股份有限公司;Go Taq?DNA Polymerase、Ex Taq?DNA Polymerase、PrimeSTAR?Max DNA Polymerase均購(gòu)自TaKaRa;2×Phanta?Max Master Mix(Dye Plus)和Plasmid Mini Kit均由南京諾唯贊生物科技股份有限公司生產(chǎn);Gel Extraction Kit和Cycle-Pure kit購(gòu)自O(shè)MEGA;溶壁酶購(gòu)自Sigma;M5 Fungal Genomic DNA Kit購(gòu)于北京聚合美生物科技有限公司;二氨基聯(lián)苯胺(diaminobenzidine,DAB)顯色試劑盒和氯化硝基四氮唑藍(lán)(nitro blue tetrazolium,NBT)試劑購(gòu)自Solarbio,4-硝基苯基-β-D-吡喃葡萄糖苷(pNPG)和4-硝基苯基-β-D-纖維二糖(pNPC)均由Aladdin生產(chǎn)。

    1.2 PaLPMO11生物信息學(xué)分析

    從基因組數(shù)據(jù)庫(kù)和CAZy數(shù)據(jù)庫(kù)中獲得5個(gè)假定的基因、、、和,分別命名為和,對(duì)應(yīng)的氨基酸序列為PaLPMO11A (XP_001906381.1)、PaLPMO11B (XP_001912736.1)、PaLPMO11C (XP_001911580------.1)、PaLPMO11D (XP_001911666.1)和PaLPMO11E (XP_001909751.1)。在CAZy數(shù)據(jù)庫(kù)和NCBI基因組數(shù)據(jù)庫(kù)(https://www.ncbi.nlm.nih.gov/)中檢索得到7條來(lái)自真菌子囊菌和輪蟲(chóng)的LPMO11蛋白序列。分別為:(BAE61530.1)、(XP 748042.1)、(CCT67099.1)、(EAA31701.1)、(XP 750980.1)、(EAA35129.1)和(UJR20201.1)。通過(guò)分子進(jìn)化遺傳學(xué)分析軟件MEGA 11.0對(duì)12個(gè)蛋白質(zhì)序列進(jìn)行分析建樹(shù),使用muscle算法對(duì)編譯后的序列進(jìn)行比對(duì),并手動(dòng)調(diào)整,采用MEGA 11.0中的鄰近法(neighbor-joining method)構(gòu)建基于這些序列的系統(tǒng)進(jìn)化樹(shù),利用自展法(bootstrap)設(shè)置1000次重復(fù)檢驗(yàn)。使用I-TASSER在線網(wǎng)站(https://zhanglab. ccmb.med.umich.edu/I-TASSER/)和UCSF Chimera等軟件對(duì)5個(gè)目的基因的蛋白三級(jí)結(jié)構(gòu)和保守氨基酸位點(diǎn)進(jìn)行預(yù)測(cè)。

    1.3 構(gòu)建PaLPMO11家族突變體

    運(yùn)用Split-marker[18]方法在WT菌株中分別敲除和五個(gè)目的基因。構(gòu)建AA11家族5個(gè)單突變型菌株技術(shù)原理如圖1所示。常規(guī)PCR分別擴(kuò)增出目的基因上下游的同源臂片段(約1000 bp)以及相對(duì)應(yīng)的抗性標(biāo)記基因片段。其中選用nourseothricin標(biāo)記,選用hygromycin B標(biāo)記,和選用geneticin標(biāo)記,選用phleomycin標(biāo)記。融合PCR對(duì)常規(guī)PCR純化后得到的目的基因上下游的同源臂片段和抗性標(biāo)記基因片段進(jìn)行融合連接,構(gòu)建敲除表達(dá)盒(引物序列如表1所示)。將融合片段轉(zhuǎn)化至ΔPaMus51菌株的原生質(zhì)體后,分別使用含有nourseothricin、hygromycin、geneticin或phleomycin抗性的抗性平板篩選抗性穩(wěn)定的轉(zhuǎn)化子。轉(zhuǎn)化子與WT菌株雜交,篩選出具有目的基因抗性而不具有原生質(zhì)體抗性的小孢子,即獲得可穩(wěn)定遺傳的突變型菌株。PCR和基因測(cè)序驗(yàn)證5個(gè)單突變型菌株是否構(gòu)建成功。使用相應(yīng)的不同交配型菌株進(jìn)行遺傳雜交,雜交后篩選同時(shí)具有相應(yīng)抗性的小孢子,以此構(gòu)建AA11家族的10個(gè)雙重突變型菌株、10個(gè)三重突變型菌株、5個(gè)四重突變型菌株和五重突變型菌株,PCR驗(yàn)證具有相同抗性雙重突變型菌株ΔAΔD和五重突變體菌株是否構(gòu)建成功。

    1.4 表型分析

    (1) M2培養(yǎng)基培養(yǎng):將不同交配型的WT和31個(gè)AA11家族突變型菌株接種在M2固體培養(yǎng)基上,并于光照條件下27℃恒溫培養(yǎng)1周。在培養(yǎng)過(guò)程中觀察并記錄各個(gè)菌株的生長(zhǎng)速率、菌落大小、菌落形態(tài)、色素沉著、子實(shí)體的形成以及子囊孢子的產(chǎn)生、噴發(fā)和萌發(fā)情況,進(jìn)行顯微觀察并記錄子實(shí)體的數(shù)量。每個(gè)實(shí)驗(yàn)重復(fù)3次。

    圖1 構(gòu)建PaLPMO11家族突變體的流程圖

    (2) 木質(zhì)纖維素培養(yǎng)基培養(yǎng):將不同交配型的WT和突變型菌株混勻破碎,分別接種于添加了木質(zhì)素、木屑和微晶纖維素的M0培養(yǎng)基上,于27℃光照條件下培養(yǎng)兩周。在培養(yǎng)過(guò)程中觀察并記錄各個(gè)菌株的生長(zhǎng)速率、菌落大小、菌落形態(tài)、色素沉著、子實(shí)體的形成以及子囊孢子的產(chǎn)生、噴發(fā)和萌發(fā)情況。每個(gè)實(shí)驗(yàn)重復(fù)3次。

    1.5 菌絲ROS水平分析

    將WT、ΔE和ΔAΔBΔCΔDΔE不同交配型菌株混勻破碎后分別接種在M2培養(yǎng)基上,于27℃黑暗培養(yǎng)3天,分別使用DAB染色試劑盒和NBT染色試劑對(duì)菌絲進(jìn)行染色,檢測(cè)菌絲組織中過(guò)氧化物和超氧化物的積累[19]。在過(guò)氧化物酶的存在下,過(guò)氧化氫與DAB反應(yīng)會(huì)生成棕褐色沉淀;NBT在超氧化物酶催化下,與超氧離子反應(yīng)生成藍(lán)紫色沉淀。菌絲染色強(qiáng)度與過(guò)氧化物和超氧化物的積累成正比關(guān)系[20]。置于水平搖床避光孵育3 h,菌絲體出現(xiàn)顏色變化時(shí)除去染色液,使用PBS洗滌3次,終止顯色反應(yīng)。觀察比對(duì)各菌株菌絲染色程度。每個(gè)實(shí)驗(yàn)重復(fù)3次。

    1.6 纖維素酶活測(cè)定

    將不同交配型的WT、ΔE和ΔAΔBΔCΔDΔE菌株分別接種到纖維素誘導(dǎo)培養(yǎng)基上生長(zhǎng),以促進(jìn)纖維素作用酶的分泌。27℃、150 r/min、光照條件下培養(yǎng)3天、5天、7天、9天,取上清液離心后即得粗酶液。參照國(guó)標(biāo)QB2583-2003及Dashtban[21]和Thankappan[22]的方法對(duì)濾紙酶活(filter paper activity,F(xiàn)PA)、內(nèi)切葡萄糖苷酶活(endoglucosidase activity,EG activity)、外切葡聚糖酶活(cellobiohydrolase activity,CBH activity)和β-葡萄糖苷酶活(β-D-glucosidase activity,BG activity)進(jìn)行測(cè)定,具體操作方法如下。每組樣品均設(shè)置3組重復(fù)。

    表1 本研究使用的引物

    1.6.1 FPA和EG activity測(cè)定

    FPA和EG酶活單位(U)定義為:在實(shí)驗(yàn)反應(yīng)條件下,每分鐘釋放1 μmol葡萄糖所需的酶量。纖維素酶可以降解纖維素底物,生成還原性單糖。還原性單糖與3,5-二硝基水楊酸(DNS)沸水浴后發(fā)生顯色反應(yīng),生成棕紅色氨基化合物。在一定范圍內(nèi),反應(yīng)顏色深淺與酶解反應(yīng)單位時(shí)間內(nèi)釋放的還原性單糖濃度呈正比,在540 nm波長(zhǎng)處可以檢測(cè)到其吸光度。根據(jù)葡萄糖標(biāo)準(zhǔn)曲線計(jì)算酶解反應(yīng)釋放的還原性單糖含量,從而計(jì)算酶活。FPA實(shí)驗(yàn)以濾紙為反應(yīng)底物,EG activity實(shí)驗(yàn)以CMC-Na為反應(yīng)底物。將粗酶液與底物混合,50℃水浴1 h,進(jìn)行反應(yīng)。加入DNS后輕輕混合,沸水浴5 min后立即置于冰上冷卻,取200 μL反應(yīng)液于96孔板中,540 nm處測(cè)定吸光度。每個(gè)樣品均設(shè)置空白對(duì)照組。

    1.6.2 CBH activity和BG activity測(cè)定

    CBH和BG酶活單位(U)定義為:在實(shí)驗(yàn)反應(yīng)條件下,每分鐘釋放1 μmol硝基苯酚(pNP)所需的酶量。無(wú)色底物pNPG 和pNPC在β-葡萄糖苷酶和外切葡聚糖酶的催化下,水解成有色物質(zhì)pNP,在405 nm波長(zhǎng)處有最大吸收峰。根據(jù)pNP標(biāo)準(zhǔn)曲線計(jì)算酶解反應(yīng)釋放的pNP含量,從而計(jì)算酶活。將粗酶液與底物混合,50℃水浴30 min,進(jìn)行反應(yīng)。水浴后立即置于冰上冷卻,向各離心管中加入2% Na2CO3終止反應(yīng)。取200 μL反應(yīng)液于96孔板中,405 nm處測(cè)定吸光度。每個(gè)樣品均設(shè)置空白對(duì)照組。

    2 結(jié)果與分析

    2.1 PaLPMO11家族蛋白的生物信息學(xué)分析

    2.1.1 PaLPMO11的系統(tǒng)發(fā)育分析

    真菌米曲霉(4MAI_A)和煙曲霉(7P3U_A)編碼的LPMO11蛋白相似度為39.63%,由此推測(cè)AA11家族LPMO氨基酸----序列差異雖然較大,但構(gòu)成組氨酸支架的兩個(gè)組氨酸和AA11家族獨(dú)有的軸位點(diǎn)酪氨酸側(cè)鏈均嚴(yán)格保守,該酪氨酸側(cè)鏈不直接參與金屬銅離子的配位,而是與其附近的丙氨酸相互作用[11]。5個(gè)基因編碼的蛋白PaLPMO11A、PaLPMO11B、PaLPMO11C、PaLPMO11D和PaLPMO11E與AoAA11蛋白以及AfAA11蛋白均有較高的一致性(均在27.9%以上)。與AoAA11蛋白相似度分別為37.04%、27.93%、35.94%、45.07%和40.64%;與AfAA11蛋白相似度分別為37.43%、49.02%、48.6%、38.39%和49.05%。因此,推測(cè)5個(gè)編碼的蛋白可能是與AoAA11和AfAA11蛋白類(lèi)似的裂解多糖單加氧酶。使用MEGA 11.0對(duì)12個(gè)LPMO11蛋白構(gòu)建系統(tǒng)發(fā)育樹(shù)(圖2A)。系統(tǒng)發(fā)育分析顯示,PaLPMO11D和AoAA11、AfAA11B和FfAA11親緣關(guān)系近;AfAA11A和PaLPMO11C處于同一分支上,表明親緣關(guān)系最近;其次是與PaLPMO11E關(guān)系較近;而PaLPMO11A和PaLPMO11B位于該進(jìn)化樹(shù)的最外側(cè)。由圖2B可知,5個(gè)PaLPMO11蛋白在催化位點(diǎn)核心保守區(qū)域顯示出較高相似度,預(yù)測(cè)兩個(gè)高度保守的組氨酸會(huì)與金屬銅離子螯合成T形的組氨酸支架,從而構(gòu)成催化活性位點(diǎn)[23]。而PaLPMO11B在多序列比對(duì)結(jié)果中顯示缺少軸位點(diǎn)酪氨酸側(cè)鏈,因此在預(yù)測(cè)蛋白三級(jí)結(jié)構(gòu)中,無(wú)法預(yù)測(cè)出其催化活性位點(diǎn)的保守氨基酸所處位置(圖3)。根據(jù)系統(tǒng)發(fā)育分析結(jié)果和多序列比對(duì)結(jié)果顯示,推測(cè)5個(gè)PaLPMO11蛋白可能分屬于AA11家族的不同亞家族。--

    2.1.2 PaLPMO11蛋白的三級(jí)結(jié)構(gòu)模型預(yù)測(cè)

    和基因編碼蛋白的三級(jí)結(jié)構(gòu)與活性位點(diǎn)保守氨基酸預(yù)測(cè)如圖3所示。5個(gè)PaLPMO11蛋白預(yù)測(cè)具有與AoAA11蛋白相似的整體拓?fù)浣Y(jié)構(gòu),均具有典型的反向平行β-三明治核心結(jié)構(gòu)。將除PaLPMO11B外的4個(gè)PaLPMO11預(yù)測(cè)蛋白的催化活性位點(diǎn)的保守氨基酸與AoAA11蛋白的活性位點(diǎn)的保守氨基酸進(jìn)行比較,所有原子的均方根偏差(root mean square deviation,RMSD)在0.77~0.97 ?范圍內(nèi)顯示出較高水平的結(jié)構(gòu)重疊。根據(jù)系統(tǒng)進(jìn)化分析、多序列比對(duì)和蛋白三級(jí)結(jié)構(gòu)模型預(yù)測(cè)結(jié)果推測(cè),的5個(gè)LPMO11蛋白可能是與AoAA11類(lèi)似的蛋白,可以通過(guò)氧化還原反應(yīng)裂解糖苷鍵,從而降解某種頑固多糖[11]。---

    A: PaLPMO11蛋白的系統(tǒng)發(fā)育樹(shù)。AA11蛋白氨基酸序列來(lái)源于CAZy數(shù)據(jù)庫(kù)和NCBI網(wǎng)站,基因庫(kù)序列號(hào)顯示在系統(tǒng)發(fā)育樹(shù)中。B: LPMO11s和PaLPMO11s的部分氨基酸序列比對(duì)。三角形表示參與銅離子配位的兩個(gè)絕對(duì)保守的組氨酸殘基,以及AA11家族特有的軸向酪氨酸殘基。

    2.2 ΔPaLPMO11突變體的構(gòu)建

    2.2.1五個(gè)基因敲除表達(dá)盒的構(gòu)建

    按1.3構(gòu)建PaLPMO11家族突變體所述方法構(gòu)建突變體。利用4種不同的抗性標(biāo)記geneticin、nourseothricin、hygromycin和phleomycin代替5個(gè)基因序列,敲除這5個(gè)基因,并對(duì)5個(gè)單突變體進(jìn)行抗性驗(yàn)證和PCR驗(yàn)證。以WT基因組DNA為模板,PCR擴(kuò)增得到5個(gè)基因上游片段PaLPMO11-5′和下游片段PaLPMO11-3′約1000 bp片段。以質(zhì)粒pBC- Hygromycin、pBC-Nourseothricin、pBC-Geneticin或pBC-Phleomycin DNA為模板,擴(kuò)增獲得HygroR基因片段(3657 bp)、NourR基因片段(1773 bp)、GeneR基因片段(2342 bp)或PhleR基因片段(3100 bp)。將各個(gè)基因的上下游片段和抗性標(biāo)記基因片段進(jìn)行融合PCR擴(kuò)增,分別構(gòu)建的敲除表達(dá)盒-5′-HygroR和HygroR--3′;的敲除表達(dá)盒-5′-NourR和NourR--3′;的敲除表達(dá)盒-5′-GeneRand GeneR--3′;的敲除表達(dá)盒-5′-PhleR和PhleR--3′以及的敲除表達(dá)盒-5′-GeneR和GeneR--3′。5個(gè)基因敲除表達(dá)盒的構(gòu)建結(jié)果見(jiàn)附圖1和附圖2。通過(guò)原生質(zhì)體轉(zhuǎn)化方法,將敲除表達(dá)盒轉(zhuǎn)化到ΔPaMus51菌株的原生質(zhì)體中,篩選具有相應(yīng)抗性的轉(zhuǎn)化子,利用兩對(duì)驗(yàn)證引物(verify__1F和5T、verify__2F和3T)對(duì)篩選出來(lái)的轉(zhuǎn)化子進(jìn)行PCR驗(yàn)證(驗(yàn)證引物序列見(jiàn)表1)并純化,擴(kuò)增得到的同源重組片段與目的序列大小一致,序列比對(duì)結(jié)果與目標(biāo)序列一致。以上結(jié)果均表明,目的基因片段與抗性基因片段在目的敲除基因位點(diǎn)發(fā)生同源重組置換,證實(shí)抗性標(biāo)記基因成功整合在目的基因座中,成功構(gòu)建PaLPMO11家族5個(gè)基因缺失菌株。

    A: 目標(biāo)蛋白三級(jí)結(jié)構(gòu)預(yù)測(cè)。B:PaLPMO11蛋白與AoAA11蛋白的整體疊加。AoAA11蛋白用黃色標(biāo)注,PaLPMO11蛋白用藍(lán)色標(biāo)注。C:PaLPMO11蛋白(藍(lán)色)與AoAA11蛋白(黃色)的活性位點(diǎn)的重疊。

    2.2.2 ΔPaLPMO11突變體的PCR驗(yàn)證

    單突變體構(gòu)建成功后,通過(guò)遺傳雜交構(gòu)建10個(gè)雙突變型菌株ΔAΔB(ΔPaLPMO11AΔΔPaLPMO11B)、ΔAΔC、ΔAΔD、ΔAΔE、ΔBΔC、ΔBΔD、ΔBΔE、ΔCΔD、ΔCΔE和ΔDΔE,10個(gè)三重突變型菌株ΔAΔBΔC、ΔAΔBΔD、ΔAΔBΔE、ΔAΔCΔD、ΔAΔCΔE、ΔAΔDΔE、ΔBΔCΔD、ΔBΔCΔE、ΔBΔDΔE和ΔCΔDΔE,5個(gè)四重突變型菌株ΔAΔBΔCΔD、ΔAΔBΔCΔE、ΔAΔBΔDΔE、ΔAΔCΔDΔE和ΔBΔCΔDΔE以及缺失所有PaLPMO11家族基因的五重突變體菌株ΔAΔBΔCΔDΔE。具有相同抗性雙重突變型菌株ΔAΔD和五重突變體菌株ΔAΔBΔCΔDΔE均通過(guò)PCR驗(yàn)證。ΔPaLPMO11突變體的PCR驗(yàn)證結(jié)果見(jiàn)附圖3。

    2.3 PaLPMO11家族基因的缺失影響P. anserina利用碳源的能力

    2.3.1 菌株在M2培養(yǎng)基上的生長(zhǎng)發(fā)育

    M2培養(yǎng)基是以糊精作為唯一碳源、能滿(mǎn)足生長(zhǎng)和發(fā)育的基礎(chǔ)培養(yǎng)基。將WT和突變型菌株接種到M2培養(yǎng)基上,27℃培養(yǎng)1周。除ΔAΔBΔCΔDΔE外的突變型菌株均不受影響,沒(méi)有發(fā)現(xiàn)WT和其他突變型菌株在菌落大小、形態(tài)和色素沉著、子實(shí)體形成以及子囊孢子產(chǎn)生、噴發(fā)和萌發(fā)方面存在明顯差異(數(shù)據(jù)未顯示)。而ΔAΔBΔCΔDΔE出現(xiàn)色素沉著的時(shí)間和子實(shí)體生成時(shí)間與其他菌株出現(xiàn)差異,且菌落大小、菌絲密度、子實(shí)體形成數(shù)量及子囊孢子噴發(fā)數(shù)量都遠(yuǎn)低于WT (圖4)。與WT相比,ΔAΔBΔCΔDΔE的生長(zhǎng)有最顯著的延遲,且提前進(jìn)入衰老狀態(tài)。培養(yǎng)第7天起,ΔAΔBΔCΔDΔE開(kāi)始停止生長(zhǎng)并伴隨明顯的色素沉著。結(jié)果表明,PaLPMO11家族單個(gè)基因在的生長(zhǎng)發(fā)育過(guò)程中似乎沒(méi)有發(fā)揮關(guān)鍵作用,但5個(gè)基因同時(shí)敲除,會(huì)引起菌株生長(zhǎng)發(fā)育和有性生殖明顯延緩并導(dǎo)致菌株提前衰老。

    2.3.2 菌株在木質(zhì)素培養(yǎng)基上的生長(zhǎng)發(fā)育

    能夠高效利用木質(zhì)纖維素作為碳源,通過(guò)分泌胞外酶使其能夠有效地滲透和降解植物生物質(zhì)[24]。為了確定PaLPMO11家族基因在木質(zhì)纖維素利用中的作用,將WT與31個(gè)突變型菌株進(jìn)行表型分析比較。用木質(zhì)素、木屑和微晶纖維素等量代替糊精作為唯一的碳源,評(píng)估突變體利用木質(zhì)纖維素的能力。

    在以純木質(zhì)素作為唯一碳源的培養(yǎng)基上培養(yǎng)時(shí),WT菌絲呈紡絲狀生長(zhǎng),在培養(yǎng)基上只形成菌絲體,而不生成子實(shí)體[25]。除ΔAΔBΔCΔDΔE外,其余菌株生長(zhǎng)速度很快,菌絲體在培養(yǎng)6天后幾乎覆蓋了整個(gè)培養(yǎng)基。在木質(zhì)素平板上,ΔAΔBΔCΔDΔE的生長(zhǎng)直徑與菌絲密度明顯減少(圖5),在培養(yǎng)5天時(shí),生長(zhǎng)直徑比WT減少了19%;在培養(yǎng)7天時(shí),直徑減少了近31%。結(jié)果表明,ΔAΔBΔCΔDΔE利用木質(zhì)素的能力降低。除ΔAΔBΔCΔDΔE外,其他突變型菌株的生長(zhǎng)與WT幾乎沒(méi)有差異。這表明,5個(gè)基因之間可能在降解木質(zhì)素以獲取營(yíng)養(yǎng)物質(zhì)的過(guò)程中具有補(bǔ)償效應(yīng),在PaLPMO11家族基因一個(gè)或多個(gè)基因缺失的情況下,菌株的生長(zhǎng)狀態(tài)并沒(méi)有顯著改變,而PaLPMO11家族的基因全部缺失時(shí),菌株的生長(zhǎng)和菌絲量都會(huì)相應(yīng)地下降減少。

    2.3.3 菌株在木屑培養(yǎng)基上的生長(zhǎng)發(fā)育

    在以木屑為唯一碳源的培養(yǎng)基上培養(yǎng)時(shí),除ΔAΔBΔCΔDΔE外,WT和其他突變型菌株的子實(shí)體在培養(yǎng)6天后開(kāi)始發(fā)育,數(shù)量眾多,并在2天后噴發(fā)具有活性的子囊孢子。且除ΔAΔBΔCΔDΔE以外,其他突變型菌株的子實(shí)體萌發(fā)時(shí)間、生長(zhǎng)速率和菌絲密度與WT并無(wú)明顯差異。而5個(gè)基因的同時(shí)缺失會(huì)導(dǎo)致菌株生長(zhǎng)速率延緩,子囊孢子的形成時(shí)間延遲,子實(shí)體的產(chǎn)生顯著減少,這與M2基礎(chǔ)培養(yǎng)基上的觀察結(jié)果一致。

    圖4 WT與ΔAΔBΔCΔDΔE菌株在M2培養(yǎng)基上的生長(zhǎng)發(fā)育情況

    將不同交配型菌株混合破碎,取5 μL菌懸液滴在M2培養(yǎng)基中央,27℃光照培養(yǎng)3天(d)后,拍照并記錄菌株生長(zhǎng)直徑。野生型與其他突變株間生長(zhǎng)直徑并無(wú)顯著差異。檢驗(yàn):*<0.05,**<0.01,***<0.001。

    圖5 WT與ΔAΔBΔCΔDΔE菌株在木質(zhì)素培養(yǎng)基上的生長(zhǎng)發(fā)育情況

    將不同交配型菌株混合破碎,取5 μL菌懸液滴在木質(zhì)素培養(yǎng)基中央,27℃光照培養(yǎng)3天(d)后,拍照并記錄菌株生長(zhǎng)直徑。圖中數(shù)字表示突變型菌株與野生型菌株相比,生長(zhǎng)直徑減少的百分比。野生型與其他突變株間生長(zhǎng)直徑并無(wú)顯著差異。檢驗(yàn):*<0.05,**<0.01,***<0.001。

    圖6 WT與ΔAΔBΔCΔDΔE菌株在木屑培養(yǎng)基上的生長(zhǎng)發(fā)育情況

    將不同交配型菌株混合破碎,取5 μL菌懸液滴在木屑培養(yǎng)基中央,27℃光照培養(yǎng)9天(d)后,拍照并觀察培養(yǎng)基上及培養(yǎng)皿蓋上的子實(shí)體密度。子實(shí)體為可見(jiàn)的黑色小點(diǎn)。野生型與其他突變株間子實(shí)體數(shù)量并無(wú)明顯差異。

    絲狀真菌的子實(shí)體是區(qū)分營(yíng)養(yǎng)生殖和有性生殖的重要生殖器官[26];在真菌生命周期中,子實(shí)體的形成是最為復(fù)雜的發(fā)育過(guò)程之一,需要消耗大量的能量和營(yíng)養(yǎng)物質(zhì)[27,28]。觀察噴發(fā)到培養(yǎng)基蓋子上的子囊孢子密度(圖7A),發(fā)現(xiàn)敲除基因的突變型菌株噴發(fā)的孢子數(shù)明顯少于WT和其他突變體,ΔAΔBΔCΔDΔE的產(chǎn)孢數(shù)是最少的(圖片未完全展示)。為了量化菌株產(chǎn)孢數(shù),對(duì)WT和敲除基因的16個(gè)突變型菌株進(jìn)行子實(shí)體數(shù)量統(tǒng)計(jì)與顯微鏡觀察(圖7B)。每個(gè)菌株均設(shè)置了3組重復(fù)。由圖7D可知,ΔBΔCΔE子實(shí)體數(shù)量為5739,ΔAΔBΔCΔE子實(shí)體數(shù)量為5296,ΔAΔCΔDΔE子實(shí)體數(shù)量為4812,ΔAΔBΔCΔDΔE子實(shí)體數(shù)量為379,與WT的子實(shí)體數(shù)量7216相比,分別減少了20%、27%、33%和95%,均存在顯著差異。統(tǒng)計(jì)上述菌株子囊孢子的萌發(fā)率,WT菌株子囊孢子萌發(fā)率為95.5%,ΔBΔCΔE子囊孢子萌發(fā)率為92%,ΔAΔBΔCΔE子囊孢子萌發(fā)率為81.5%,ΔAΔCΔDΔE子囊孢子萌發(fā)率為70%,ΔAΔBΔCΔDΔE子囊孢子萌發(fā)率為51.5% (圖7E)。子實(shí)體是真菌體內(nèi)的產(chǎn)孢結(jié)構(gòu),孢子在子實(shí)體內(nèi)形成并于成熟后噴發(fā),是真菌進(jìn)行繁殖的主要途徑[29]。ΔAΔBΔCΔDΔE菌株進(jìn)入有性生殖階段后,會(huì)生成部分異常的子實(shí)體,部分子實(shí)體無(wú)內(nèi)容物,不含有活性的子囊孢子(圖7C),因而無(wú)法正常排出子囊孢子,導(dǎo)致噴發(fā)到培養(yǎng)基蓋子上的子囊孢子密度遠(yuǎn)小于其他菌株。實(shí)驗(yàn)結(jié)果顯示,基因缺失的突變體的子囊孢子萌發(fā)率和子實(shí)體數(shù)量均呈現(xiàn)下降趨勢(shì),說(shuō)明基因可能與子囊孢子萌發(fā)和子實(shí)體形成有關(guān);除ΔAΔBΔCΔDΔE外,WT和其余突變體的生長(zhǎng)速率和菌絲密度無(wú)明顯差異。

    圖7 WT和部分PaLPMO11E基因缺失突變型菌株的子囊孢子密度圖、顯微觀察、子實(shí)體數(shù)量統(tǒng)計(jì)與子囊孢子萌發(fā)率統(tǒng)計(jì)

    A:野生型和部分基因缺失突變型菌株的子囊孢子密度圖。菌株于木屑培養(yǎng)基上培養(yǎng)8天后會(huì)將孢子噴發(fā)至培養(yǎng)皿蓋子上,對(duì)接有子囊孢子的培養(yǎng)皿蓋拍照。B:野生型和突變型菌株于倒置顯微鏡下的觀察結(jié)果。不同交配型菌株于M2培養(yǎng)基27℃光照培養(yǎng)5天后,置于倒置顯微鏡下觀察子實(shí)體形態(tài)與數(shù)量。標(biāo)尺:500 μm。C:WT及ΔΔBΔCΔDΔE子實(shí)體放大圖。在M2培養(yǎng)基27℃光照培養(yǎng)5天后,挑取子實(shí)體置于載玻片上,用蓋玻片蓋上碾壓,置于倒置顯微鏡下觀察并拍照。標(biāo)尺:100 μm。D:子實(shí)體數(shù)量統(tǒng)計(jì)。不同交配型菌株在M2培養(yǎng)基上27℃光照7天后,在顯微鏡下對(duì)野生型和突變型菌株的子實(shí)體數(shù)量進(jìn)行計(jì)數(shù)統(tǒng)計(jì)。E:子囊孢子萌發(fā)率統(tǒng)計(jì)。將不同交配型菌株在M2培養(yǎng)基上進(jìn)行雜交,培養(yǎng)7天后開(kāi)始噴發(fā)子囊孢子,用接種孢子培養(yǎng)基接住噴發(fā)出來(lái)的子囊孢子。每個(gè)菌株挑取200個(gè)子囊孢子于萌發(fā)孢子培養(yǎng)基上黑暗培養(yǎng)1天后,統(tǒng)計(jì)子囊孢子的萌發(fā)率。檢驗(yàn): *<0.05,**<0.01,***<0.001。

    這表明基因可能參與了有性生殖過(guò)程,而PaLPMO11家族5個(gè)基因的同時(shí)缺失,會(huì)導(dǎo)致菌株?duì)I養(yǎng)生長(zhǎng)和有性生殖過(guò)程進(jìn)一步受損,不僅表現(xiàn)為菌落的生長(zhǎng)直徑明顯降低和子實(shí)體數(shù)量大幅減少,還表現(xiàn)為生成異常子實(shí)體,導(dǎo)致菌株的低育性。由此推測(cè)上述突變型菌株可能無(wú)法從木屑中正常獲取足夠的碳源,利用營(yíng)養(yǎng)物質(zhì)的能力受損,導(dǎo)致菌株?duì)I養(yǎng)生長(zhǎng)期間能量?jī)?chǔ)備不足,無(wú)法正常供應(yīng)有性生殖過(guò)程的能量需求,從而造成菌株有性生殖過(guò)程部分受損,其中ΔAΔBΔCΔDΔE的表型最為明顯。

    2.3.4 菌株在微晶纖維素平板上的生長(zhǎng)發(fā)育

    在以微晶纖維素為唯一碳源的培養(yǎng)基培養(yǎng)時(shí),野生型菌株生成的子實(shí)體以接種點(diǎn)為中心呈環(huán)狀分布[20],通過(guò)計(jì)算培養(yǎng)皿上的成熟子實(shí)體數(shù)量來(lái)衡量菌株的生育力。由圖8A可知,敲除基因的突變型菌株形成更大更彌散的子實(shí)體環(huán),子實(shí)體分布呈現(xiàn)發(fā)散狀,且子實(shí)體數(shù)量少于WT,這一點(diǎn)與菌株在木屑平板上生長(zhǎng)得到的結(jié)果一致。這一表型在ΔAΔBΔCΔDΔE菌株上尤為明顯,ΔAΔBΔCΔDΔE的子實(shí)體數(shù)量與WT相比下降了19倍。在其他突變型菌株中沒(méi)有觀察到這種表型,它們所生成的子實(shí)體環(huán)與WT無(wú)太大差異,均在培養(yǎng)皿中央生成成熟的子實(shí)體環(huán)。實(shí)驗(yàn)結(jié)果表明,基因可能在的有性生殖過(guò)程中發(fā)揮了作用,基因缺失菌株會(huì)形成一個(gè)更大更彌散的子實(shí)體環(huán),子實(shí)體密度減少(圖8B)。由實(shí)驗(yàn)結(jié)果推測(cè),基因缺失菌株可能無(wú)法充分利用微晶纖維素碳源以維持自身生長(zhǎng)發(fā)育所需,導(dǎo)致有性生殖過(guò)程在一定程度上受到影響,因而生成的子實(shí)體數(shù)量減少。PaLPMO11家族基因的同時(shí)敲除,會(huì)對(duì)的生長(zhǎng)發(fā)育和有性繁殖造成更大程度的影響,導(dǎo)致菌株生長(zhǎng)發(fā)育明顯延緩,生成空心子實(shí)體,子實(shí)體數(shù)量與子囊孢子萌發(fā)率顯著下降。而PaLPMO11家族基因間可能發(fā)揮相互代償?shù)淖饔?,一個(gè)基因或多個(gè)基因缺失時(shí),PaLPMO11家族的其他基因可能發(fā)揮類(lèi)似的功能,以代償其他基因的缺失,從而保證菌株可以維持正常的生理過(guò)程。

    2.4 過(guò)量生成的過(guò)氧化氫可能加速了菌株的衰老

    活性氧(ROS)是細(xì)胞有氧代謝過(guò)程中產(chǎn)生的副產(chǎn)物,主要由超氧陰離子自由基和過(guò)氧化氫組成[30]。生物體內(nèi)存在一系列抗氧化系統(tǒng)使ROS的產(chǎn)生和消除處于平衡狀態(tài),當(dāng)ROS在生物體內(nèi)過(guò)量產(chǎn)生時(shí),會(huì)開(kāi)始攻擊細(xì)胞內(nèi)脂質(zhì)、蛋白質(zhì)和DNA,導(dǎo)致細(xì)胞損傷和凋亡[30,31]。對(duì)WT、ΔE和ΔAΔBΔCΔDΔE菌株菌絲的ROS水平進(jìn)行分析,染色結(jié)果如圖9所示。NBT染色觀察發(fā)現(xiàn),WT和ΔE菌絲的藍(lán)紫色顯色最深,且WT和ΔE之間菌絲染色并無(wú)明顯差異,而ΔAΔBΔCΔDΔE藍(lán)紫色染色最淺,說(shuō)明ΔAΔBΔCΔDΔE超氧化物的分泌量明顯減少。DAB染色觀察發(fā)現(xiàn),除培養(yǎng)第1天外,DAB染色結(jié)果與NBT染色結(jié)果一致,WT和ΔE染色程度無(wú)差異,均呈現(xiàn)明顯的棕褐色,ΔAΔBΔCΔDΔE菌絲染色程度最淺。培養(yǎng)1天時(shí),ΔAΔBΔCΔDΔE菌絲密度低于WT,但其DAB染色程度卻深于WT,結(jié)果表明ΔAΔBΔCΔDΔE菌株可能在培養(yǎng)第1天時(shí)分泌了大量的過(guò)氧化氫。由上述結(jié)果推測(cè),單突變型菌株體內(nèi)可能存在其他相似功能的基因能彌補(bǔ)其缺失,因而其染色結(jié)果與WT無(wú)明顯差異。再結(jié)合在M2和木質(zhì)素平板上的表型分析結(jié)果,ΔAΔBΔCΔDΔE菌株均出現(xiàn)提前衰老的現(xiàn)象。由此推測(cè),ΔAΔBΔCΔDΔE衰老的潛在機(jī)制可能是由于菌株生長(zhǎng)前期體內(nèi)積累了過(guò)量過(guò)氧化氫,而過(guò)量的過(guò)氧化氫可能會(huì)導(dǎo)致菌株進(jìn)入過(guò)度氧化應(yīng)激狀態(tài),導(dǎo)致分子氧化損傷增加,從而促進(jìn)了突變菌株的衰老過(guò)程,導(dǎo)致突變型菌株提前進(jìn)入衰老狀態(tài)。

    2.5 PaLPMO11基因參與纖維素降解過(guò)程

    纖維素可以用于乙醇生產(chǎn)等工業(yè)生物煉制應(yīng)用,而纖維素的降解利用卻受到如物化特性、結(jié)構(gòu)特點(diǎn)以及組成成分等因素的阻礙[32,33]。因此,提高木質(zhì)纖維素的降解效率、減少工業(yè)生產(chǎn)成本對(duì)于木質(zhì)纖維素的工業(yè)應(yīng)用具有重要意義。濾紙酶活是包括內(nèi)切葡聚糖酶、外切葡聚糖酶和β-葡萄糖苷酶三種纖維素降解酶在內(nèi)的總纖維素酶活[34]。單突變體ΔE的四種纖維素酶活力與WT差異不大,而ΔAΔBΔCΔDΔE的酶活降低最為顯著。培養(yǎng)7天時(shí),ΔAΔBΔCΔDΔE的FPA酶活是WT酶活的53%,EG酶活是WT酶活的61%,CBH酶活是WT酶活的61%,BG酶活是WT酶活的32%;菌株培養(yǎng)9天時(shí),ΔAΔBΔCΔDΔE的FPA活性、EG活性、CBH活性和BG活性分別是WT酶活的45%、58%、36%和35% (圖10)。實(shí)驗(yàn)結(jié)果表明,單個(gè)基因的缺失,對(duì)總纖維素酶活影響不明顯,而5個(gè)基因的同時(shí)缺失,總纖維素酶活顯著下降,但卻仍有野生型45%以上的活力,說(shuō)明基因之間可能存在功能冗余,且5個(gè)基因功能可能與纖維素降解有關(guān)。

    圖8 WT與突變型菌株在微晶纖維素培養(yǎng)基上的生長(zhǎng)發(fā)育情況與對(duì)比圖

    子實(shí)體為可見(jiàn)的黑色小點(diǎn)。將不同交配型菌株混合破碎,取5 μL菌懸液滴在微晶纖維素培養(yǎng)基中央,培養(yǎng)8天后拍照。除基因缺失突變型菌株外,未發(fā)現(xiàn)野生型與其他突變株間子實(shí)體數(shù)量與分布存在顯著差異。A:野生型和基因缺失突變型菌株于微晶纖維素培養(yǎng)基培養(yǎng)的結(jié)果圖。B:野生型和突變型菌株放大對(duì)比圖。

    3 討論

    目前僅有來(lái)自米曲霉、富士鐮刀-菌和煙曲霉的LPMO11蛋白被純化和表征[11~13]。所有已知真菌的LPMO均具有保守的中央反向平行β-折疊核心結(jié)構(gòu),而暴露在平坦核心結(jié)構(gòu)表面的金屬銅離子會(huì)與兩個(gè)高度保守的組氨酸提供的3個(gè)氮原子配位,形成T形的組氨酸支架,從而構(gòu)成催化活性位點(diǎn)[10,35]。全基因組序列中存在5個(gè)可能編碼LPMO11蛋白序列的基因,分別為和。進(jìn)行生物信息學(xué)分析發(fā)現(xiàn),5個(gè)PaLPMO11蛋白均含有構(gòu)成催化活性位點(diǎn)的兩個(gè)高度保守的組氨酸殘基和β-三明治折疊核心結(jié)構(gòu),而PaLPMO11B缺乏AA11家族獨(dú)有的軸位點(diǎn)酪氨酸側(cè)鏈。因此推測(cè)5個(gè)PaLPMO11蛋白可能屬于AA11家族的不同亞家族,可以通過(guò)氧化還原反應(yīng)裂解糖苷鍵,從而降解某種難降解多糖。

    圖9 DAB和NBT染色分析野生型菌株與突變型菌株菌絲ROS水平

    將不同交配型菌株混合破碎,取5 μL菌懸液滴在M2培養(yǎng)基中央。黑暗培養(yǎng)1、2、3天(d)后,用DAB和NBT對(duì)菌株菌絲進(jìn)行染色。

    將100 μL菌懸液接種于50 mL 微晶纖維素液體培養(yǎng)基中,150 r/min恒溫?fù)u床培養(yǎng)3、5、7、9天(d),提取粗酶液,測(cè)定不同菌株的纖維素酶活。檢驗(yàn):*<0.05,**<0.01,***<0.001。

    多項(xiàng)研究表明,是一種壽命有限的絲狀真菌,其衰老的表型為:菌株生長(zhǎng)速率和生殖能力的下降、氣生菌絲形成的減少以及生長(zhǎng)停滯會(huì)伴隨著明顯的色素沉著[36,37]。本研究也發(fā)現(xiàn)ΔAΔBΔCΔDΔE在M2和木質(zhì)素培養(yǎng)基上均表現(xiàn)出異常的菌落形態(tài)和顯著的生長(zhǎng)減少,與WT的蓬松菌落表面不同,ΔAΔBΔCΔDΔE菌株的環(huán)狀菌落僅形成稀疏的氣生菌絲;生殖能力顯著下降,僅在內(nèi)圈生成少數(shù)子實(shí)體;培養(yǎng)5天后,ΔAΔBΔCΔDΔE加速了外圈色素的過(guò)度積累(圖4,圖5)。子實(shí)體是真菌體內(nèi)的產(chǎn)孢結(jié)構(gòu),孢子在子實(shí)體內(nèi)形成并于成熟后噴發(fā),是真菌繁殖的主要途徑,該過(guò)程需要消耗大量的能量和營(yíng)養(yǎng)物質(zhì)[26,29]。如圖6、圖7和圖8所示,基因可能與子實(shí)體分布、形成和子囊孢子萌發(fā)有關(guān)。同時(shí)敲除PaLPMO11家族的5個(gè)基因,會(huì)導(dǎo)致菌株子實(shí)體數(shù)量明顯減少和生成部分異常子實(shí)體,進(jìn)而造成有性生殖部分受損。ΔAΔBΔCΔDΔE菌株子實(shí)體數(shù)相較于WT減少95%,子囊孢子萌發(fā)率僅為51.5%,說(shuō)明PaLPMO11家族基因可能是控制菌株有性生殖的關(guān)鍵基因。適量的過(guò)氧化氫可以作為一種信號(hào)分子,發(fā)揮一定的生理作用,但高濃度的過(guò)氧化氫可能會(huì)引發(fā)氧化應(yīng)激,導(dǎo)致程序性細(xì)胞死亡[38]。根據(jù)前人的研究報(bào)道,用400 μmol/L的過(guò)氧化氫會(huì)誘導(dǎo)H9C2細(xì)胞發(fā)生氧化應(yīng)激,導(dǎo)致細(xì)胞活力下降[39];在中,WT菌株體內(nèi)過(guò)氧化氫釋放量會(huì)隨著培養(yǎng)時(shí)間的延長(zhǎng)而不斷增加,菌株培養(yǎng)后期釋放的過(guò)氧化氫量約是培養(yǎng)前期的3倍,證實(shí)了衰老與氧化應(yīng)激的增加有關(guān)[40];在水稻中,基因的突變會(huì)導(dǎo)致過(guò)氧化氫的積累,從而導(dǎo)致葉片過(guò)早衰老[41]。在M2培養(yǎng)基上培養(yǎng)的第1天,ΔAΔBΔCΔDΔE菌株產(chǎn)生了大量過(guò)氧化氫(圖9),再結(jié)合前面表型分析所得結(jié)果推測(cè),過(guò)量的過(guò)氧化氫可能引發(fā)了ΔAΔBΔCΔDΔE菌株過(guò)度氧化應(yīng)激狀態(tài),增強(qiáng)了自由基對(duì)細(xì)胞的氧化損傷,從而促進(jìn)了突變菌株的衰老過(guò)程。綜上,中5個(gè)基因的完全缺失會(huì)導(dǎo)致菌株出現(xiàn)嚴(yán)重的生長(zhǎng)缺陷、形態(tài)變化和壽命縮短等情況。因而,PaLPMO11家族基因在絲狀真菌的生長(zhǎng)發(fā)育和有性生殖過(guò)程中可能是必不可少的。根據(jù)實(shí)驗(yàn)結(jié)果推測(cè)PaLPMO11家族基因功能冗余,5個(gè)基因間可能相互協(xié)調(diào)調(diào)控碳源利用與菌株的生長(zhǎng)發(fā)育、有性生殖和衰老,一個(gè)或多個(gè)基因缺失時(shí),其他基因可能在一定程度上能替代或補(bǔ)充其功能,使能夠應(yīng)對(duì)AA11家族單個(gè)或多個(gè)基因的缺失。在41種常見(jiàn)的真菌基因組中,是含有最多降解木質(zhì)纖維素蛋白酶類(lèi)基因的絲狀真菌[15],具有強(qiáng)大的木質(zhì)纖維降解能力,能夠?qū)⒛举|(zhì)纖維素轉(zhuǎn)化成如生物乙醇等不同的高附加值產(chǎn)品[42]。將作為真菌降解木質(zhì)纖維素的模式生物進(jìn)行研究,發(fā)現(xiàn)一些新的木質(zhì)纖維素降解相關(guān)基因和酶,闡明這些酶和基因在真菌中的作用機(jī)制,有望為深入了解真菌降解復(fù)雜多糖、提高生物能源生產(chǎn)效率等領(lǐng)域的研究提供參考[16,43]。如圖10所示,單突變體的濾紙酶活、內(nèi)切葡聚糖酶活、外切葡聚糖酶活和β-葡萄糖苷酶活與WT沒(méi)有顯著差異,而5個(gè)基因的同時(shí)缺失,會(huì)導(dǎo)致菌株總纖維素酶活顯著下降,但卻仍有野生型45%以上的酶活力。在中,基因缺失突變體降解纖維素的效率比WT提高了近兩倍[44];和漆酶基因亞家族1的突變體降解纖維素的能力明顯下降[25,45];證實(shí)了它們均參與了纖維素降解過(guò)程。在本研究中,出乎意料的是,PaLPMO11家族基因家族貢獻(xiàn)了接近一半的纖維素酶活力。根據(jù)結(jié)果推測(cè),基因之間可能存在功能冗余,且PaLPMO11家族的5個(gè)基因功能可能與纖維素降解有關(guān),PaLPMO11家族的裂解多糖單加氧酶的完全失活可能會(huì)改變的纖維素降解機(jī)制。

    綜上所述,中AA11家族的裂解多糖單加氧酶基因參與菌株的生長(zhǎng)發(fā)育、有性生殖、衰老和纖維素降解過(guò)程。本研究結(jié)果為系統(tǒng)闡述絲狀真菌中木質(zhì)纖維素降解的調(diào)控機(jī)制提供參考。

    附加材料見(jiàn)文章電子版www.chinagene.cn。

    [1] Tran MH, Yu JH, Lee EY. Microwave-assisted two-step liquefaction of acetone-soluble lignin of silvergrass saccharification residue for production of biopolyol and biopolyurethane., 2021, 13(9): 1491.

    [2] Francois JM, Alkim C, Morin N. Engineering microbial pathways for production of bio-based chemicals from lignocellulosic sugars: current status and perspectives., 2020, 13: 118.

    [3] Zhang ZR, Song JL, Han BX. Catalytic transformation of lignocellulose into chemicals and fuel products in ionic liquids., 2017, 117(10): 6834–80.

    [4] Upton BM, Kasko AM. Strategies for the conversion of lignin to high-value polymeric materials: review and perspective., 2016, 116(4): 2275–306.

    [5] Lyu LT, Chu YD, Zhang SF, Zhang Y, Huang QT, Wang S, Zhao ZK. Engineering the oleaginous yeastfor improved resistance against inhibitors in biomass hydrolysates., 2021, 9: 768934.

    [6] Couturier M, Navarro D, Chevret D, Henrissat B, Piumi F, Ruiz-Due?as FJ, Martinez AT, Grigoriev IV, Riley R, Lipzen A, Berrin JG, Master ER, Rosso MN. Enhanced degradation of softwood versus hardwood by the white-rot fungus., 2015, 8: 216.

    [7] Xia W, Xu XX, Qian LC, Shi PJ, Bai YG, Luo HY, Ma R, Yao B. Engineering a highly active thermophilic β-glucosidase to enhance its pH stability and saccharification performance., 2016, 9: 147.

    [8] Calderaro F, Keser M, Akeroyd M, Bevers LE, Eijsink VGH, Várnai A, van den Berg MA. Characterization of an AA9 LPMO from,LPMO9B, under industrially relevant lignocellulose saccharifi-cation conditions., 2020, 13(1): 195.

    [9] Theibich YA, Sauer SPA, Leggio LL, Hedeg?rd ED. Estimating the accuracy of calculated electron paramagnetic resonance hyperfine couplings for a lytic polysaccharide monooxygenase., 2021, 19: 555–67.

    [10] Span EA, Marletta MA. The framework of polysaccharide monooxygenase structure and chemistry., 2015, 35: 93–99.

    [11] Hemsworth GR, Henrissat B, Davies GJ, Walton PH. Discovery and characterization of a new family of lytic polysaccharide monooxygenases., 2014, 10(2): 122–126.

    [12] Wang DM, Li J, Salazar-Alvarez G, McKee LS, Srivastava V, Sellberg JA, Bulone V, Hsieh YSY. Production of functionalised chitins assisted by fungal lytic polysaccharide monooxygenase., 2018, 20(9): 2091–2100.

    [13] Rieder L, Petrovi? D, V?ljam?e P, Eijsink VGH, S?rlie M. Kinetic characterization of a putatively chitin-active LPMO reveals a preference for soluble substrates and absence of monooxygenase activity., 2021, 11(18): 11685–11695.

    [14] St?pamo FG, R?hr ?K, Mekasha S, Petrovi? DM, Várnai A, Eijsink VGH. Characterization of a lytic polysa-ccharide monooxygenase fromshows functional variation among family AA11 fungal LPMOs., 2021, 297(6): 101421.

    [15] Levasseur A, Drula E, Lombard V, Coutinho PM, Henrissat B. Expansion of the enzymatic repertoire of the CAZy database to integrate auxiliary redox enzymes., 2013, 6(1): 41.

    [16] Dicko M, Ferrari R, Tangthirasunun N, Gautier V, Lalanne C, Lamari F, Silar P. Lignin degradation and its use in signaling development by the coprophilous ascomycete, 2020, 6(4): 278.

    [17] El-Khoury R, Sellem CH, Coppin E, Boivin A, Maas MFPM, Debuchy R, Sainsard-Chanet A. Gene deletion and allelic replacement in the filamentous fungus., 2008, 53(4): 249–258.

    [18] Goswami RS. Targeted gene replacement in fungi using a split-marker approach., 2012, 835: 255–269.

    [19] Yan LY, Zhang HJ, Zheng YQ, Cong YQ, Liu CT, Fan F, Zheng C, Yuan GL, Pan G, Yuan DY, Duan MJ. Transcription factor OsMADS25 improves rice tolerance to cold stress., 2021, 43(11): 1078– 1087.閆凌月,張豪健,鄭雨晴,叢韞起,劉次桃,樊帆,鄭鋮,袁貴龍,潘根,袁定陽(yáng),段美娟. 轉(zhuǎn)錄因子OsMADS25提高水稻對(duì)低溫的耐受性. 遺傳,2021, 43(11): 1078–1087.

    [20] Tangthirasunun N, Navarro D, Garajova S, Chevret D, Tong LCH, Gautier V, Hyde KD, Silar P, Berrin JG. Inactivation of cellobiose dehydrogenases modifies the cellulose degradation mechanism of, 2017, 83(2): e02716–16.

    [21] Dashtban M, Maki M, Leung KT, Mao CQ, Qin WS. Cellulase activities in biomass conversion: measurement methods and comparison., 2010, 30(4): 302–309.

    [22] Thankappan S, Kandasamy S, Joshi B, Sorokina KN, Taran OP, Uthandi S. Bioprospecting thermophilic glycosyl hydrolases, from hot springs of Himachal Pradesh, for biomass valorization., 2018, 8(1): 168.

    [23] Bennati-Granier C, Garajova S, Champion C, Grisel S, Haon M, Zhou SM, Fanuel M, Ropartz D, Rogniaux H, Gimbert I, Record E, Berrin JG. Substrate specificity and regioselectivity of fungal AA9 lytic polysaccharide monooxygenases secreted by., 2015, 8: 90.

    [24] Demoor A, Silar P, Brun S. Appressorium: the breakth-rough in., 2019, 5(3): 72.

    [25] Xie N, Ruprich-Robert G, Silar P, Chapeland-Leclerc F. Bilirubin oxidase-like proteins from: promising thermostable enzymes for application in transformation of plant biomass., 2015, 17(3): 866–875.

    [26] Fu C, Thielhelm TP, Heitman J. Unisexual reproduction promotes competition for mating partners in the global human fungal pathogen., 2019, 15(9): e1008394.

    [27] Wang M, Gu BL, Huang J, Jiang S, Chen YJ, Yin YL, Pan YF, Yu GJ, Li YM, Wong BHC, Liang Y, Sun H. Transcriptome and proteome exploration to provide a resource for the study of., 2013, 8(2): e56686.

    [28] Orban A, Weber A, Herzog R, Hennicke F, Rühl M. Transcriptome of different fruiting stages in the cultivated mushroomsuggests a complex regulation of fruiting and reveals enzymes putatively involved in fungal oxylipin biosynthesis., 2021, 22(1): 324.

    [29] Schumacher DI, Lütkenhaus R, Altegoer F, Teichert I, Kück U, Nowrousian M. The transcription factor PRO44 and the histone chaperone ASF1 regulate distinct aspects of multicellular development in the filamentous fungus., 2018, 19(1): 112.

    [30] Ismail HF, Hashim Z, Soon WT, Rahman NSA, Zainudin AN, Majid FAA. Comparative study of herbal plants on the phenolic and flavonoid content, antioxidant activities and toxicity on cells and zebrafish embryo., 2017, 7(4): 452–465.

    [31] Zhang HR, Li BX, Sun ZM, Zhou H, Zhang SS. Integration of intracellular telomerase monitoring by electrochemiluminescence technology and targeted cancer therapy by reactive oxygen species., 2017, 8(12): 8025–8029.

    [32] Kalyani D, Tiwari MK, Li JL, Kim SC, Kalia VC, Kang YC, Lee JK. A highly efficient recombinant laccase from the yeastand its application in the hydrolysis of biomass., 2015, 10(3): e0120156.

    [33] Abraham RE, Wong CS, Puri M. Enrichment of cellulosic waste hemp () hurd into non-toxic microfibres., 2016, 9(7): 562.

    [34] Zhang XJ, Qu YB, Qin YQ. Expression and chromatin structures of cellulolytic enzyme gene regulated by heterochromatin protein 1., 2016, 9: 206.

    [35] Laurent CVFP, Sun PC, Scheiblbrandner S, Csarman F, Cannazza P, Frommhagen M, van Berkel WJH, Oostenbrink C, Kabel MA, Ludwig R. Influence of lytic polysaccharide monooxygenase active site segments on activity and affinity., 2019, 20(24): 6219.

    [36] Silar P.: from laboratory to biotechnology. Genomics of soil- and plant-associated fungi. 2013, 283–309.

    [37] Bernhardt D, Hamann A, Osiewacz HD. The role of mitochondria in fungal aging., 2014, 22: 1–7.

    [38] Pellavio G, Sommi P, Anselmi-Tamburini U, DeMichelis MP, Coniglio S, Laforenza U. Cerium oxide nanoparticles regulate oxidative stress in hela cells by increasing the aquaporin-mediated hydrogen peroxide permeability., 2022, 23(18): 10837.

    [39] Gao H, Liu XP, Tian KM, Meng YC, Yu CC, Peng YF. Insight into the protective effect of salidroside against H2O2-induced injury in H9C2 cells., 2021, 2021: 1060271.

    [40] Wiemer M, Osiewacz HD. Effect of paraquat-induced oxidative stress on gene expression and aging of the filamentous ascomycete., 2014, 1(7): 225–240.

    [41] Wu LW, Ren DY, Hu SK, Li GM, Dong GJ, Jiang L, Hu XM, Ye WJ, Cui YT, Zhu L, Hu J, Zhang GH, Gao ZY, Zeng DL, Qian Q, Guo LB. Down-regulation of a nicotinate phosphoribosyltransferase gene,, leads to withered leaf tips., 2016, 171(2): 1085–1098.

    [42] Bian HY, Wu XX, Luo J, Qiao YZ, Fang GG, Dai HQ. Valorization of alkaline peroxide mechanical pulp by metal chloride-assisted hydrotropic pretreatment for enzymatic saccharification and cellulose nanofibrillation., 2019, 11(2): 331.

    [43] Couturier M, Tangthirasunun N, Ning X, Brun S, Gautier V, Bennati-Granier C, Silar P, Berrin JG. Plant biomass degrading ability of the coprophilic ascomycete fungus., 2016, 34(5): 976–983.

    [44] Brun S, Malagnac F, Bidard F, Lalucque H, Silar P. Functions and regulation of the Nox family in the filamentous fungus: a new role in cellulose degradation., 2009, 74(2): 480–496.

    [45] Xie N, Chapeland-Leclerc F, Silar P, Ruprich-Robert G. Systematic gene deletions evidences that laccases are involved in several stages of wood degradation in the filamentous fungus., 2014, 16(1): 141–161.

    Identification and functional study of AA11 family polysaccharide monooxygenase genes in filamentous fungus

    Wenzhen Du1, Yuanjing Li2, Jialing Wu1, Siyu Chen1, Liang Jiang1, Gang Liu1, Ning Xie1

    The lytic polysaccharide monooxygenase (LPMO) in the auxiliary active protein family (AA family) catalyzes the oxidative depolymerization of various refractory carbohydrates including cellulose, chitin and starch. While accumulating studies investigate the enzymology of LPMO, the research on the inactivation ofgenes has been rarely explored. In this study, fivegenes()()()()()of the AA11 family in the filamentous funguswere knocked out by homologous recombination. Single mutants ΔPaLPMO11A (ΔA), ΔPaLPMO11B (ΔB), ΔPaLPMO11C (ΔC), ΔPaLPMO11D (ΔD) and ΔPaLPMO11E (ΔE) were constructed, and then all polygenic mutants were constructed via genetic crosses. The differences in the growth rate and sexual reproduction between wild type and mutant strains were observed on different carbon source media. The alteration of oxidative stress and cellulose degradation ability were found on DAB and NBT staining and cellulase activity determination. These results implicated thatgenes play a key role in the growth, development, and lignocellulose degradation of. The results showed that the spore germination efficiency, growth rate and reproductive capacity of mutant strains including ΔBΔCΔE, ΔAΔBΔCΔE, ΔAΔCΔDΔE and ΔAΔBΔCΔDΔE was significantly decreased on different cellulose carbon sources and the remaining strains have no difference. The reduced utilization of various carbon sources, the growth rate, the spore germination rate, the number of fruiting bodies, the normal fruiting bodies, the shortened life span and the ability to degrade cellulose were found in strains which all five genes in the PaLPMO11 family were deleted. However, the strain still had 45% cellulase activity compared to wild type. These results suggest thatgenes may be involved in the growth and development, sexual reproduction, senescence and cellulose degradation of. This study provides information for systematically elucidating the regulatory mechanism of lignocellulose degradation in filamentous fungus.

    filamentous fungi;;; gene knockout; lignocellulose degradation

    2023-08-21;

    2023-10-28;

    2023-11-10

    國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(編號(hào):2021YFA0910800),廣東省自然科學(xué)基金項(xiàng)目(編號(hào):2021A1515012166,2021A1515012118)和國(guó)家自然科學(xué)基金項(xiàng)目(編號(hào):31601014,22078199)資助[Supported by the National Key R&D Program of China (No. 2021YFA0910800), the Basic and Applied Basic Research Fund of Guangdong Province (Nos. 2021A1515012166, 2021A1515012118), and the National Natural Science Foundation of China (Nos. 31601014, 22078199)]

    杜文珍,碩士研究生,專(zhuān)業(yè)方向:微生物遺傳學(xué)。E-mail: 2100251038@email.szu.edu.cn

    謝寧,博士,副教授,博士生導(dǎo)師,研究方向:微生物遺傳學(xué)。E-mail: ning.xie@szu.edu.cn

    10.16288/j.yczz.23-223

    (責(zé)任編委: 劉鋼)

    猜你喜歡
    子囊孢子突變型突變體
    破壁方式對(duì)冠突散囊菌有性繁殖體提取物抑菌活性的影響
    冬蟲(chóng)夏草菌室內(nèi)侵染蝙蝠蛾幼蟲(chóng)的代謝分析
    冬蟲(chóng)夏草菌室內(nèi)侵染蝙蝠蛾幼蟲(chóng)分析
    菩提樹(shù)黑痣病菌生物學(xué)特性及病害發(fā)生規(guī)律
    CLIC1及其點(diǎn)突變體與Sedlin蛋白的共定位研究
    擬南芥干旱敏感突變體篩選及其干旱脅迫響應(yīng)機(jī)制探究
    Survivin D53A突變體對(duì)宮頸癌細(xì)胞增殖和凋亡的影響
    表皮生長(zhǎng)因子受體非突變型非小細(xì)胞肺癌分子靶治療有效1病例報(bào)道及相關(guān)文獻(xiàn)復(fù)習(xí)
    CD41-42突變型β地中海貧血重組載體pEGFP-C2-CD41-42的構(gòu)建及其穩(wěn)定轉(zhuǎn)染HeLa細(xì)胞模型的建立
    突變型PUMA(S10A)對(duì)Hela細(xì)胞的凋亡作用及其分子機(jī)制
    亚洲人成伊人成综合网2020| 女人被狂操c到高潮| 熟女少妇亚洲综合色aaa.| 一级毛片精品| 天天躁夜夜躁狠狠躁躁| 丰满的人妻完整版| 亚洲aⅴ乱码一区二区在线播放 | 成熟少妇高潮喷水视频| 亚洲av片天天在线观看| 日本五十路高清| 亚洲久久久国产精品| 18美女黄网站色大片免费观看| 在线看三级毛片| 久久国产亚洲av麻豆专区| 精品少妇一区二区三区视频日本电影| 色综合亚洲欧美另类图片| 日日干狠狠操夜夜爽| 日韩一卡2卡3卡4卡2021年| 亚洲 欧美一区二区三区| 一级黄色大片毛片| 国产蜜桃级精品一区二区三区| 国产高清有码在线观看视频 | 特大巨黑吊av在线直播 | 色精品久久人妻99蜜桃| 老司机福利观看| 久久久久免费精品人妻一区二区 | 免费看a级黄色片| 黄色 视频免费看| 国产国语露脸激情在线看| 日韩一卡2卡3卡4卡2021年| 无人区码免费观看不卡| 在线av久久热| 特大巨黑吊av在线直播 | 美女午夜性视频免费| 久久中文字幕一级| 免费在线观看亚洲国产| 国产伦在线观看视频一区| 久久欧美精品欧美久久欧美| a级毛片a级免费在线| 99精品久久久久人妻精品| 黄色a级毛片大全视频| 一区二区三区激情视频| 老汉色av国产亚洲站长工具| 欧美乱妇无乱码| 亚洲天堂国产精品一区在线| 不卡一级毛片| 变态另类丝袜制服| 精品一区二区三区四区五区乱码| 真人做人爱边吃奶动态| 日韩成人在线观看一区二区三区| 欧美性猛交╳xxx乱大交人| 精品一区二区三区四区五区乱码| 免费在线观看视频国产中文字幕亚洲| 亚洲第一av免费看| 99国产精品一区二区蜜桃av| 国产高清videossex| 亚洲精品一卡2卡三卡4卡5卡| 国产精品香港三级国产av潘金莲| www.www免费av| 正在播放国产对白刺激| 妹子高潮喷水视频| 视频区欧美日本亚洲| 国产成人一区二区三区免费视频网站| 国产精品亚洲美女久久久| 国产三级在线视频| 69av精品久久久久久| 午夜福利成人在线免费观看| 欧美不卡视频在线免费观看 | 最近最新中文字幕大全免费视频| 亚洲第一电影网av| 丝袜在线中文字幕| 天天添夜夜摸| 很黄的视频免费| 人人妻人人看人人澡| 男人舔女人的私密视频| 丝袜美腿诱惑在线| 熟女少妇亚洲综合色aaa.| 午夜激情av网站| 草草在线视频免费看| 母亲3免费完整高清在线观看| 亚洲avbb在线观看| 少妇 在线观看| 一进一出好大好爽视频| 最近在线观看免费完整版| 一级毛片女人18水好多| 一进一出抽搐gif免费好疼| 久久香蕉国产精品| 非洲黑人性xxxx精品又粗又长| 在线国产一区二区在线| av欧美777| 亚洲精华国产精华精| 欧美日韩亚洲国产一区二区在线观看| 亚洲熟女毛片儿| 日韩欧美三级三区| 国产视频内射| 亚洲国产精品合色在线| 欧美激情 高清一区二区三区| 成年人黄色毛片网站| 国产激情欧美一区二区| 伊人久久大香线蕉亚洲五| 91成年电影在线观看| 99久久精品国产亚洲精品| 91九色精品人成在线观看| 91成人精品电影| 51午夜福利影视在线观看| 亚洲精品国产一区二区精华液| 成年版毛片免费区| 久久精品国产清高在天天线| 免费在线观看黄色视频的| 人人妻人人看人人澡| 无遮挡黄片免费观看| 最近最新中文字幕大全免费视频| 伦理电影免费视频| 久久久久久人人人人人| 一个人免费在线观看的高清视频| 亚洲国产中文字幕在线视频| 在线观看免费视频日本深夜| videosex国产| 亚洲欧美精品综合久久99| 国产精品99久久99久久久不卡| 午夜免费观看网址| 精品无人区乱码1区二区| 久久久久免费精品人妻一区二区 | 中文字幕精品亚洲无线码一区 | 亚洲精品美女久久av网站| 免费在线观看黄色视频的| 国产亚洲精品av在线| 亚洲欧美日韩无卡精品| 亚洲黑人精品在线| 精品第一国产精品| 精品少妇一区二区三区视频日本电影| 日韩三级视频一区二区三区| 国产亚洲欧美98| 又黄又爽又免费观看的视频| 亚洲第一青青草原| 99精品久久久久人妻精品| 国产亚洲精品av在线| 99re在线观看精品视频| 在线观看免费日韩欧美大片| 好看av亚洲va欧美ⅴa在| 欧美在线黄色| 在线观看免费日韩欧美大片| 无限看片的www在线观看| 国内揄拍国产精品人妻在线 | 国产黄片美女视频| av有码第一页| 婷婷丁香在线五月| 亚洲成人久久性| 亚洲av电影在线进入| 国产v大片淫在线免费观看| 久久国产精品影院| 天天添夜夜摸| 国产人伦9x9x在线观看| 男女午夜视频在线观看| 中国美女看黄片| 50天的宝宝边吃奶边哭怎么回事| 丁香欧美五月| 亚洲国产精品成人综合色| 国产人伦9x9x在线观看| 午夜福利18| 日韩av在线大香蕉| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产精品999在线| 大型av网站在线播放| 精品午夜福利视频在线观看一区| 免费一级毛片在线播放高清视频| 久久精品人妻少妇| 91av网站免费观看| 我的亚洲天堂| 欧美日本亚洲视频在线播放| av超薄肉色丝袜交足视频| 黄片播放在线免费| 国产野战对白在线观看| 禁无遮挡网站| 一本大道久久a久久精品| 久热这里只有精品99| 午夜免费激情av| 中文字幕高清在线视频| 巨乳人妻的诱惑在线观看| 在线观看午夜福利视频| 一级毛片精品| 一区二区三区高清视频在线| 亚洲熟妇中文字幕五十中出| 日韩av在线大香蕉| 老司机午夜十八禁免费视频| 淫妇啪啪啪对白视频| 在线观看免费视频日本深夜| 黄色视频,在线免费观看| 熟女少妇亚洲综合色aaa.| 婷婷精品国产亚洲av| av在线播放免费不卡| 天天添夜夜摸| 制服诱惑二区| 黄色丝袜av网址大全| 久久久久免费精品人妻一区二区 | 欧美在线一区亚洲| 免费一级毛片在线播放高清视频| 久久99热这里只有精品18| 亚洲成人久久爱视频| 男男h啪啪无遮挡| 国产精品综合久久久久久久免费| 中文字幕人成人乱码亚洲影| 禁无遮挡网站| 亚洲国产精品合色在线| 19禁男女啪啪无遮挡网站| 成年女人毛片免费观看观看9| 变态另类成人亚洲欧美熟女| av有码第一页| 日本三级黄在线观看| 久久精品91蜜桃| 女性生殖器流出的白浆| 91大片在线观看| 亚洲精品一区av在线观看| 午夜福利在线观看吧| 国产1区2区3区精品| 国产成人系列免费观看| 久9热在线精品视频| 日韩国内少妇激情av| 丰满的人妻完整版| 特大巨黑吊av在线直播 | 国产真实乱freesex| 亚洲自偷自拍图片 自拍| 男女床上黄色一级片免费看| 国产精品久久电影中文字幕| 亚洲自拍偷在线| 禁无遮挡网站| 白带黄色成豆腐渣| videosex国产| 50天的宝宝边吃奶边哭怎么回事| 久久99热这里只有精品18| 免费在线观看亚洲国产| 亚洲国产欧美日韩在线播放| 免费看十八禁软件| 久久久久国产精品人妻aⅴ院| 亚洲精品国产精品久久久不卡| 熟女少妇亚洲综合色aaa.| 亚洲全国av大片| 夜夜爽天天搞| 国产精品亚洲美女久久久| 亚洲成av片中文字幕在线观看| 亚洲第一欧美日韩一区二区三区| 欧美又色又爽又黄视频| 精品欧美国产一区二区三| 国产精品乱码一区二三区的特点| videosex国产| 怎么达到女性高潮| 宅男免费午夜| 成熟少妇高潮喷水视频| 亚洲美女黄片视频| 一级作爱视频免费观看| 熟女电影av网| 99国产精品一区二区蜜桃av| 亚洲一区二区三区不卡视频| 免费高清视频大片| 欧美中文日本在线观看视频| 老司机在亚洲福利影院| 国产不卡一卡二| 黄频高清免费视频| 天堂影院成人在线观看| 夜夜看夜夜爽夜夜摸| 88av欧美| 久久这里只有精品19| 99国产精品99久久久久| 国产成人系列免费观看| 母亲3免费完整高清在线观看| 麻豆一二三区av精品| 国产精品久久久人人做人人爽| 人人妻人人澡人人看| 在线观看免费视频日本深夜| 人成视频在线观看免费观看| 不卡av一区二区三区| 在线天堂中文资源库| 免费在线观看成人毛片| 又黄又粗又硬又大视频| 午夜福利欧美成人| 欧美三级亚洲精品| 欧美黑人欧美精品刺激| 在线天堂中文资源库| 老鸭窝网址在线观看| 一本精品99久久精品77| 999精品在线视频| 亚洲熟妇熟女久久| 国产国语露脸激情在线看| 中文字幕高清在线视频| 日韩成人在线观看一区二区三区| 亚洲激情在线av| 久久天躁狠狠躁夜夜2o2o| 在线播放国产精品三级| 黄色成人免费大全| 18禁国产床啪视频网站| tocl精华| 又黄又粗又硬又大视频| 一a级毛片在线观看| 国产精品久久视频播放| 久久热在线av| 日韩欧美 国产精品| 天天添夜夜摸| 特大巨黑吊av在线直播 | 亚洲欧美一区二区三区黑人| 国产精品久久电影中文字幕| 亚洲人成网站在线播放欧美日韩| 一二三四社区在线视频社区8| 禁无遮挡网站| 日日夜夜操网爽| 国产av在哪里看| 黄色视频,在线免费观看| 特大巨黑吊av在线直播 | 亚洲色图 男人天堂 中文字幕| 色老头精品视频在线观看| 极品教师在线免费播放| 精品卡一卡二卡四卡免费| 天堂√8在线中文| 天堂影院成人在线观看| 香蕉av资源在线| 精品久久久久久成人av| 巨乳人妻的诱惑在线观看| 夜夜夜夜夜久久久久| 久久久久久九九精品二区国产 | 精品久久久久久久久久免费视频| 免费av毛片视频| 精品久久久久久成人av| 亚洲精品久久成人aⅴ小说| 99久久国产精品久久久| 女人爽到高潮嗷嗷叫在线视频| 亚洲欧美日韩无卡精品| 最近在线观看免费完整版| 亚洲人成网站在线播放欧美日韩| 丝袜人妻中文字幕| 亚洲午夜精品一区,二区,三区| 免费高清在线观看日韩| 中文资源天堂在线| 极品教师在线免费播放| 国产精品亚洲美女久久久| or卡值多少钱| 好看av亚洲va欧美ⅴa在| 91成年电影在线观看| 国产精品野战在线观看| 欧美性猛交黑人性爽| 高潮久久久久久久久久久不卡| 亚洲精品久久成人aⅴ小说| 国产v大片淫在线免费观看| 久久人人精品亚洲av| 很黄的视频免费| 自线自在国产av| 亚洲一卡2卡3卡4卡5卡精品中文| 妹子高潮喷水视频| 午夜精品在线福利| 国产99白浆流出| 国产亚洲av嫩草精品影院| 熟女电影av网| 青草久久国产| 免费在线观看影片大全网站| 在线观看日韩欧美| 亚洲人成网站在线播放欧美日韩| 香蕉丝袜av| 看片在线看免费视频| 999久久久国产精品视频| 国产精品 国内视频| 夜夜躁狠狠躁天天躁| 在线视频色国产色| 99精品久久久久人妻精品| 色老头精品视频在线观看| 久久久水蜜桃国产精品网| 校园春色视频在线观看| 国产成人欧美| 人成视频在线观看免费观看| 久久久久久久久中文| 丰满的人妻完整版| 夜夜夜夜夜久久久久| 在线观看www视频免费| АⅤ资源中文在线天堂| 亚洲av电影不卡..在线观看| 成人av一区二区三区在线看| 亚洲国产精品成人综合色| 人人妻人人看人人澡| 宅男免费午夜| 精品日产1卡2卡| 亚洲成人久久性| 亚洲精华国产精华精| 国产伦在线观看视频一区| АⅤ资源中文在线天堂| 丰满人妻熟妇乱又伦精品不卡| 日韩欧美一区视频在线观看| 狂野欧美激情性xxxx| 日韩欧美 国产精品| 久久久久九九精品影院| 欧美性长视频在线观看| 亚洲av日韩精品久久久久久密| 亚洲第一电影网av| 此物有八面人人有两片| 给我免费播放毛片高清在线观看| 成年版毛片免费区| 久久中文字幕一级| 最近最新中文字幕大全电影3 | a级毛片a级免费在线| 色尼玛亚洲综合影院| 人人妻人人澡欧美一区二区| 两性夫妻黄色片| 国产精品乱码一区二三区的特点| 99在线人妻在线中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 欧美+亚洲+日韩+国产| 午夜福利视频1000在线观看| www日本在线高清视频| 露出奶头的视频| 国产又黄又爽又无遮挡在线| 无遮挡黄片免费观看| av福利片在线| 国产视频一区二区在线看| 51午夜福利影视在线观看| 中文字幕高清在线视频| 一二三四社区在线视频社区8| 日韩欧美在线二视频| 亚洲精品一区av在线观看| 好男人在线观看高清免费视频 | av超薄肉色丝袜交足视频| 国产亚洲精品av在线| 女生性感内裤真人,穿戴方法视频| 色av中文字幕| 人人妻,人人澡人人爽秒播| 嫁个100分男人电影在线观看| 欧美最黄视频在线播放免费| 国产成人系列免费观看| 日本一区二区免费在线视频| 变态另类成人亚洲欧美熟女| 俄罗斯特黄特色一大片| 精品国产乱子伦一区二区三区| 99re在线观看精品视频| 日本撒尿小便嘘嘘汇集6| 国产99白浆流出| 一级片免费观看大全| 午夜日韩欧美国产| 亚洲激情在线av| 色av中文字幕| 午夜福利18| 中文字幕精品亚洲无线码一区 | 亚洲国产看品久久| 亚洲成人久久爱视频| 妹子高潮喷水视频| 1024视频免费在线观看| 91九色精品人成在线观看| 国产av不卡久久| 精品第一国产精品| 婷婷六月久久综合丁香| 国产97色在线日韩免费| 99久久久亚洲精品蜜臀av| 中文字幕精品亚洲无线码一区 | 中文字幕人成人乱码亚洲影| 欧美在线一区亚洲| 亚洲av中文字字幕乱码综合 | 国产久久久一区二区三区| 亚洲精品在线观看二区| 人人妻人人澡人人看| 日韩欧美 国产精品| 一级作爱视频免费观看| 身体一侧抽搐| av超薄肉色丝袜交足视频| 无限看片的www在线观看| 90打野战视频偷拍视频| 亚洲欧美日韩高清在线视频| 草草在线视频免费看| 91麻豆av在线| 欧美日本亚洲视频在线播放| 麻豆成人午夜福利视频| 久久香蕉精品热| 麻豆av在线久日| 女性被躁到高潮视频| 国产亚洲精品久久久久5区| 黑人操中国人逼视频| 久久久水蜜桃国产精品网| 黄片播放在线免费| 1024香蕉在线观看| 成人18禁在线播放| avwww免费| 中出人妻视频一区二区| 一区二区日韩欧美中文字幕| 欧美最黄视频在线播放免费| 男女之事视频高清在线观看| 女同久久另类99精品国产91| 三级毛片av免费| 久久久久久人人人人人| 真人一进一出gif抽搐免费| 亚洲中文字幕日韩| 色播亚洲综合网| 欧美日韩瑟瑟在线播放| 日本a在线网址| 国产免费男女视频| 色综合站精品国产| 黄色毛片三级朝国网站| 变态另类丝袜制服| 可以在线观看毛片的网站| 久久精品影院6| 怎么达到女性高潮| 天堂√8在线中文| 高潮久久久久久久久久久不卡| 他把我摸到了高潮在线观看| 亚洲avbb在线观看| 国产一区在线观看成人免费| 少妇的丰满在线观看| 午夜免费成人在线视频| bbb黄色大片| 精品电影一区二区在线| av福利片在线| 久久久久久人人人人人| 欧美日韩一级在线毛片| 午夜视频精品福利| av福利片在线| 国内少妇人妻偷人精品xxx网站 | av在线播放免费不卡| 国产亚洲欧美精品永久| 18禁黄网站禁片免费观看直播| 人人妻人人澡人人看| 婷婷精品国产亚洲av在线| 欧美中文综合在线视频| 琪琪午夜伦伦电影理论片6080| 在线十欧美十亚洲十日本专区| 久久九九热精品免费| 欧美日本亚洲视频在线播放| 最新美女视频免费是黄的| 天天添夜夜摸| 国产亚洲精品久久久久久毛片| 亚洲成人久久爱视频| 在线观看一区二区三区| 国产激情偷乱视频一区二区| 男女视频在线观看网站免费 | 亚洲av美国av| 女人爽到高潮嗷嗷叫在线视频| 国产伦人伦偷精品视频| 亚洲一区中文字幕在线| 亚洲黑人精品在线| 精品少妇一区二区三区视频日本电影| 久久精品91无色码中文字幕| 国语自产精品视频在线第100页| 免费搜索国产男女视频| 一级a爱视频在线免费观看| 18美女黄网站色大片免费观看| 一级a爱视频在线免费观看| 99久久无色码亚洲精品果冻| 日本熟妇午夜| 51午夜福利影视在线观看| 欧美国产精品va在线观看不卡| 久久中文字幕一级| 国产成+人综合+亚洲专区| 香蕉丝袜av| 成人永久免费在线观看视频| av视频在线观看入口| 久热爱精品视频在线9| 国产精华一区二区三区| 国产成人欧美| 久久久精品国产亚洲av高清涩受| 日韩欧美国产一区二区入口| 老司机午夜十八禁免费视频| 久久国产精品男人的天堂亚洲| 久久午夜亚洲精品久久| 怎么达到女性高潮| АⅤ资源中文在线天堂| 国产精品电影一区二区三区| 国产精品精品国产色婷婷| 黑人操中国人逼视频| 哪里可以看免费的av片| 人人澡人人妻人| 国产精品九九99| 99久久99久久久精品蜜桃| 日本a在线网址| 久久久久久久久中文| 满18在线观看网站| 动漫黄色视频在线观看| 悠悠久久av| 别揉我奶头~嗯~啊~动态视频| 色在线成人网| 精品久久久久久成人av| 日本精品一区二区三区蜜桃| 国产熟女xx| 日本 欧美在线| 国产午夜福利久久久久久| 免费在线观看亚洲国产| 一a级毛片在线观看| 黑人欧美特级aaaaaa片| or卡值多少钱| av电影中文网址| 亚洲黑人精品在线| 亚洲精品在线美女| 久久久精品国产亚洲av高清涩受| 美女大奶头视频| www.熟女人妻精品国产| 黄片小视频在线播放| 国产精品野战在线观看| 色综合站精品国产| 久久国产亚洲av麻豆专区| 90打野战视频偷拍视频| 亚洲av片天天在线观看| 亚洲无线在线观看| 久久久久免费精品人妻一区二区 | 热99re8久久精品国产| 日韩成人在线观看一区二区三区| 日本在线视频免费播放| 国产精品一区二区三区四区久久 | 黄色 视频免费看| 99久久精品国产亚洲精品| 男女做爰动态图高潮gif福利片| 观看免费一级毛片| 亚洲中文av在线| 香蕉久久夜色| av视频在线观看入口| 亚洲成人免费电影在线观看| 日本一区二区免费在线视频| 美女大奶头视频| 最近最新中文字幕大全电影3 | 国产亚洲欧美98| 黄片大片在线免费观看| 亚洲第一av免费看| 看黄色毛片网站| 欧美乱妇无乱码| 久热这里只有精品99| 一本大道久久a久久精品|