岳倩倩，趙永霞，王先桂，周禮紅

（1.茅臺(tái)學(xué)院 實(shí)驗(yàn)實(shí)訓(xùn)教學(xué)中心，貴州 遵義 564507；2.遵義醫(yī)科大學(xué) 藥學(xué)院，貴州 遵義 563003；3.貴州大學(xué) 生命科學(xué)學(xué)院，貴州 貴陽 550025）

紅曲霉（Monascus）作為一種藥食用絲狀真菌，其次級代謝產(chǎn)物包括天然紅曲色素、γ-氨基丁酸、莫納可林（Monacolin）K、桔霉素等[1-2]。紅曲色素作為一種聚酮類代謝產(chǎn)物，現(xiàn)已有110多種被發(fā)現(xiàn)[3]，根據(jù)其溶解性可分為醇溶性和水溶性2種[4]。紅曲色素具有食品著色劑的生物安全性和適用性[5]，已被用于研制具有功能特性的調(diào)味品[6]、酒類[7]、肉制品[8]和面制品[9]等。除此之外，紅曲色素具有羥基、醚鍵、羰基等官能團(tuán)和共軛結(jié)構(gòu)，具有抗突變[10]、抗癌[11]、抗肥胖[12-13]、抗氧化[14]、降血壓[15]和抗阿爾茲海默癥[4]等多種生物活性，使得紅曲色素在功能性食品、制藥等方面有更廣的應(yīng)用前景。

氮源作為重要的能源物質(zhì)，對真菌生長發(fā)育過程中大分子物質(zhì)有調(diào)控作用，已有關(guān)于構(gòu)巢曲霉（Aspergillus nidulans）、粗糙脈胞菌（Neurosporacrassa）中氮代謝調(diào)控機(jī)制的研究[16-17]。岳倩倩等[18]考察不同的無機(jī)氮源對紅曲霉菌絲生長和產(chǎn)色素的調(diào)控情況發(fā)現(xiàn)，銨態(tài)氮對紅曲霉生長和產(chǎn)色素有促進(jìn)作用，亞硝態(tài)氮?jiǎng)t有抑制作用，但其氮代謝調(diào)控機(jī)制尚不清楚。

隨著第二代測序技術(shù)的出現(xiàn)，組學(xué)技術(shù)應(yīng)運(yùn)而生。到目前為止，主要的組學(xué)方法包括基因組學(xué)、轉(zhuǎn)錄組學(xué)、蛋白質(zhì)組學(xué)、代謝組學(xué)等[19-20]。其中，轉(zhuǎn)錄組學(xué)因其具有靈敏度高、費(fèi)用低、針對性強(qiáng)、對其篩出的遺傳標(biāo)記可直接進(jìn)行應(yīng)用等優(yōu)點(diǎn)而率先取得進(jìn)步[21]。因其可以用于深度挖掘新基因、確定代謝途徑及進(jìn)化分析[22-23]，已被廣泛應(yīng)用于生物學(xué)、醫(yī)學(xué)、農(nóng)學(xué)等領(lǐng)域。如稻曲病菌有性生殖中光誘導(dǎo)基因的轉(zhuǎn)錄組分析[24]、大青轉(zhuǎn)錄組測序及生物信息學(xué)分析[25]等。為深入研究氮源對紅曲霉生長和產(chǎn)色素的調(diào)控機(jī)制，本研究分別以銨態(tài)氮、亞硝態(tài)氮為氮源培養(yǎng)安卡紅曲霉（Monascus anka）菌株GZUM-25，采用Illumina二代測序平臺(tái)對紅曲霉GZUM-25進(jìn)行轉(zhuǎn)錄組測序分析，分析不同氮源條件下調(diào)控紅曲霉菌絲生長和產(chǎn)色素的差異表達(dá)基因，并對基因進(jìn)行GO功能富集及KEGG通路富集分析，以期為后續(xù)紅曲霉基因的深度挖掘提供信息資源和理論參考。

1 材料與方法

1.1 材料與試劑

1.1.1 菌種

安卡紅曲霉（Monascus anka）GZUM-25：本實(shí)驗(yàn)室分離保藏。

1.1.2 培養(yǎng)基

沙氏瓊脂改良培養(yǎng)基[25]：麥芽糖5 g/L，蛋白胨10 g/L，酵母膏5 g/L，瓊脂20 g/L，葡萄糖20 g/L，pH自然，121 ℃高壓滅菌25 min。液體培養(yǎng)基中不添加瓊脂。

察氏培養(yǎng)基[25]：NaNO33g/L，K2HPO41g/L，MgSO·47H2O 0.5 g/L，KCl 0.5 g/L，F(xiàn)eSO40.01 g/L，蔗糖30 g/L，瓊脂20 g/L，121 ℃高壓滅菌25 min。

調(diào)控培養(yǎng)基：向察氏培養(yǎng)基中分別添加1%的NaNO2、（NH4）2SO4兩種不同的無機(jī)氮源，用于紅曲霉菌株培養(yǎng)。

1.1.3 化學(xué)試劑

麥芽糖、葡萄糖、酵母膏、蛋白胨（均為生化試劑）：奧博星生物技術(shù)有限公司；NaNO2、K2HPO4、MgSO4·7H2O、KCl、FeSO4（均為分析純）：重慶北碚精細(xì)化工廠；（NH4）2SO4（分析純）：重慶川江化工有限公司。其他試劑均為國產(chǎn)分析純。

1.2 儀器與設(shè)備

LDZX-50KBS立式高壓蒸汽滅鍋：上海申安醫(yī)療器械廠；SW-CJ-1F超凈工作臺(tái)：蘇州凈化設(shè)備有限公司；DH-360A電熱恒溫培養(yǎng)箱：北京科偉永興儀器有限公司；THZ-82A數(shù)顯氣浴恒溫振蕩器：金壇市晶玻實(shí)驗(yàn)儀器廠；Trizol總核糖核酸（ribonucleic acid，RNA）抽提試劑盒、TruseqTMRNA sample prep試劑盒、Dynal M280磁珠：美國Invitrogen公司；NanoDrop2000：美國Thermo Scientific公司。

1.3 試驗(yàn)方法

1.3.1 紅曲霉樣品制備

將安卡紅曲霉菌株GZMU-25接種到沙氏瓊脂改良培養(yǎng)基上，30 ℃恒溫培養(yǎng)10 d。將已活化的安卡紅曲霉菌株GZMU-25用10 mL無菌水洗下后，將孢子懸液分別轉(zhuǎn)移到含1%NaNO2和1%（NH4）2SO4的調(diào)控培養(yǎng)基中，分別編號為MP1、MP2，30 ℃、230 r/min條件下振蕩培養(yǎng)7 d。培養(yǎng)結(jié)束后，將菌絲用無菌鑷子挑出，用液氮速凍保存于-80 ℃冰箱。

1.3.2 總RNA的提取

采用總RNA提取試劑盒提取經(jīng)液氮冷凍保藏的菌絲球總RNA，通過NanoDrop2000對所提RNA的濃度和純度進(jìn)行檢測，采用瓊脂糖凝膠電泳檢測RNA完整性，Agilent 2100測定RNA完整值（RNA integrity number，RIN），挑選滿足測序要求的RNA樣本進(jìn)行轉(zhuǎn)錄組測序。

1.3.3 文庫的構(gòu)建和Illummina Hiseq測序

利用帶有Oligo（dT）的磁珠與poly A進(jìn)行腺嘌呤（adenine，A）-胸腺嘧啶（thymine，T）配對，從樣本MP1、MP2總RNA中分離出信使核糖核酸（messenger ribonucleic acid，mRNA），加入fragmentation buffer，將mRNA隨機(jī)斷裂成堿基長度為200 bp左右的小片段。在逆轉(zhuǎn)錄酶的作用下，加入六堿基隨機(jī)引物，以mRNA為模板反轉(zhuǎn)錄合成第一鏈互補(bǔ)脫氧核糖核酸（complementary deoxyribonucleic acid，cDNA），隨后進(jìn)行第二鏈合成，形成穩(wěn)定的雙鏈結(jié)構(gòu)。對雙鏈cDNA進(jìn)行末端修飾、片段大小選擇，最后選擇富集合格的文庫進(jìn)行Illummina Hiseq測序。

1.3.4 序列分析和功能注釋

通過Illummina Hiseq測序得到包含測序接頭序列、低質(zhì)量讀段、N（無法確定堿基信息）率較高序列及長度過短序列的原始測序數(shù)據(jù)，采用軟件SeqPrep、Sickle對原始數(shù)據(jù)進(jìn)行閱讀框（reads）中接頭序列的去除、3'端低質(zhì)量堿基的修剪及含N比率＞10%reads的去除。利用bowtie將修剪后的測序數(shù)據(jù)比對到基因組上，檢測外顯子之間的可變剪切信息，并將最終得到的序列與Nr、GO、KEGG數(shù)據(jù)庫進(jìn)行對比，并對其功能進(jìn)行注釋。

2 結(jié)果與分析

2.1 紅曲霉轉(zhuǎn)錄組測序與轉(zhuǎn)錄分析

對原始數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，結(jié)果見表1。由表1可知，從樣品MP1中得到的閱讀框（reads）總數(shù)和堿基總數(shù)分別為62 127 568條和9 381 262 768 bp；從樣品MP2中得到的reads總數(shù)和堿基總數(shù)分別為52 725 002條和7 961 475 302 bp。對原始數(shù)據(jù)進(jìn)行質(zhì)量修剪后，樣品MP1得到的Clean reads和Clean堿基總數(shù)分別為61 233 790條和9 145 916 108 bp；樣品MP2得到的Clean reads 和Clean 堿基總數(shù)分別為51 890 476條和7 746 137 365 bp。樣品MP1和MP2的鳥嘌呤（guanine，G）胞嘧啶（cytosine，C）含量分別為53.29%、53.10%，Q20堿基比例分別為98.59%、98.53%，Q30堿基比例分別為95.78%、95.64%。通過與參考基因組對比后，樣品MP1和MP2得到的Uniquely mapped總數(shù)分別為45 194 295、35 139 204，占比分別為73.81%、67.72%。結(jié)果表明，該轉(zhuǎn)錄組樣本測序質(zhì)量較高，可進(jìn)行下一步樣品中差異基因的分析。

表1 轉(zhuǎn)錄組測序結(jié)果統(tǒng)計(jì)Table 1 Statistics of transcriptome sequencing analysis results

2.2 紅曲霉樣本間差異基因的分析

使用edgeR軟件，根據(jù)gene read count數(shù)據(jù)進(jìn)行差異表達(dá)計(jì)算，按照錯(cuò)誤發(fā)現(xiàn)率（false discoveryrate，F(xiàn)DR）＜0.05且差異倍數(shù)（fold change，F(xiàn)C）對數(shù)值的絕對值（｜log2FC｜）≥1為標(biāo)準(zhǔn)，對樣品差異基因進(jìn)行篩選，結(jié)果見圖1。

圖1 紅曲霉樣本間差異表達(dá)基因散點(diǎn)圖Fig.1 Scatter plot of differentially expressed genes between Monascus samples

由圖1可知，樣本MP1與MP2中共存在2 117個(gè)基因，其中存在顯著差異基因有95個(gè)，37個(gè)基因?qū)儆陲@著上調(diào)，58個(gè)基因?qū)儆陲@著下調(diào)。

2.3 紅曲霉樣本間差異表達(dá)基因的GO功能富集分析

使用Goatools軟件對紅曲霉樣品間差異表達(dá)的基因進(jìn)行GO功能顯著性富集，利用GO數(shù)據(jù)庫，以P≤0.05為篩選標(biāo)準(zhǔn)，按照基因參與的生物學(xué)過程（biological process，BP）、構(gòu)成細(xì)胞組分（cellular component，CC）、實(shí)現(xiàn)的分子功能（molecular function，MF）等將MP1、MP2樣品存在的差異表達(dá)基因進(jìn)行比對分析，結(jié)果見圖2。由圖2可知，在樣本MP1和MP2中，差異基因參與了生物學(xué)過程、細(xì)胞組分和分子功能。將95個(gè)差異基因分別映射到46個(gè)代謝通路中，在分子功能中，主要富集通路為鐵離子連接（iron ion binding）和單氧化酶活性（monooxygenase activity），均有5個(gè)差異基因參與，均占差異表達(dá)基因總數(shù)的5.26%。其次是氧化還原酶活性（oxidoreductase activity）、血紅素結(jié)合（heme binding）等。在細(xì)胞組分中，主要富集通路為葡萄糖苷酶II復(fù)合物（glucosidase II complex）、細(xì)胞表面（cell surface），參與的差異表達(dá)基因數(shù)均為1個(gè)，均占總差異表達(dá)基因總數(shù)的1%。在生物學(xué)過程中，主要富集通路為葡萄糖代謝過程（sucrose metabolic process）、淀粉代謝過程（starch metabolic process）、細(xì)胞碳水化合物代謝過程（cellular carbohydrate metabolic process），參與差異表達(dá)的基因數(shù)分別有3個(gè)、3個(gè)、4個(gè)，占總差異表達(dá)基因總數(shù)的3.16%、3.16%、4.21%。

圖2 紅曲霉樣本間差異表達(dá)基因的GO功能富集結(jié)果Fig.2 GO functional enrichment results of differentially expressed genes between Monascus samples

2.4 紅曲霉樣本間差異表達(dá)基因的KEGG通路富集分析

以P≤0.05作為篩選差異基因在KEGG通路顯著富集的標(biāo)準(zhǔn)，使用軟件KOBAS對紅曲霉樣本間差異基因進(jìn)行KEGG通路富集分析，結(jié)果見圖3。由圖3可知，在樣品MP1和MP2中，95個(gè)差異基因中有11個(gè)基因注釋到KEGG 17條通路中。包括磷酸肌醇代謝（inositol phosphate metabolism）、氨基苯基酸酯降解（aminobenzoate degradation）、淀粉和蔗糖代謝（starch and sucrose metabolism）、氯代環(huán)己烷和氯苯降解（chlorocyclohexane and chlorobenzene degradation）、緊密連接（tight junction）、纈氨酸、亮氨酸和異亮氨酸生物合成（valine，leucine and isoleucine biosynthesis）、酪氨酸代謝（tyrosine metabolism）等。在以上差異表達(dá)基因富集通路中，磷酸肌醇代謝（inositol phosphate metabolism）、淀粉和蔗糖代謝（starch and sucrose metabolism）、苯乙烯降解（styrenedegradation）以及氨基苯基酸酯降解（aminobenzoate degradation）差異基因參與數(shù)目最多，均為2個(gè)，可推測以上4條通路在紅曲霉菌絲生長和色素產(chǎn)生方面存在重要作用。

圖3 紅曲霉樣本間差異表達(dá)基因的KEGG通路富集分析結(jié)果Fig.3 KEGG pathway enrichment analysis results of differentially expressed genes between Monascus samples

3 結(jié)論

以銨態(tài)氮、亞硝態(tài)氮為氮源培養(yǎng)安卡紅曲霉（Monascus anka）菌株GZUM-25，采用Illumina二代測序平臺(tái)對紅曲霉GZUM-25進(jìn)行轉(zhuǎn)錄組測序分析發(fā)現(xiàn)，銨態(tài)氮和亞硝態(tài)氮培養(yǎng)下的安卡紅曲霉GZUM-25基因組中有95個(gè)差異表達(dá)基因，其中，上調(diào)基因有37個(gè)，下調(diào)基因有58個(gè)。通過GO功能富集及KEGG通路富集分析發(fā)現(xiàn)，差異基因主要富集的功能為生物學(xué)過程、細(xì)胞組分和分子功能，富集的通路為磷酸肌醇代謝、淀粉和蔗糖代謝、苯乙烯降解、氨基苯基酸酯降解。本研究為紅曲霉相關(guān)功能基因的挖掘提供了理論支持。