聚蘋果酸（PMLA）是以蘋果酸為唯一單體的均聚高分子聚合物，屬于聚酯類聚合物，具有高生物相容性、高水溶性、生物可吸收性、化學(xué)可衍生性、可降解性和無免疫原性等多種優(yōu)良性能，在生物醫(yī)藥、食品和生物材料領(lǐng)域具有潛在的應(yīng)用前景。研究表明產(chǎn)黑色素短梗霉（Aureobasidium melanogenum）是一種具有較強(qiáng)產(chǎn)PMLA能力的類酵母真菌，對該菌種進(jìn)行基因組測序及組裝可為改造菌種，提高產(chǎn)量提供理論依據(jù)。

1 通過二代測序數(shù)據(jù)進(jìn)行基因組預(yù)組裝

基于Illumina NovaSeq 6000測序平臺(tái)對產(chǎn)黑色素短梗霉基因組進(jìn)行測序，共得到35.85 Gb×2的雙端測序結(jié)果。且該測序平臺(tái)得到的raw reads具有較高的質(zhì)量評分（Q＞35）。之后，選取質(zhì)控后clean reads進(jìn)行后續(xù)的組裝實(shí)驗(yàn)。

不同k-mer長度對產(chǎn)黑色素短梗霉clean reads進(jìn)行組裝結(jié)果見表1，組裝結(jié)果顯示，產(chǎn)黑色素短梗霉基因組大小約為44 Mb。其中，在k-mer長度為127時(shí)得到N50值為908694，L50值為14，基因組覆蓋度為121.2×。基于該結(jié)果可以判斷該菌基因組長度理論值在44 Mb左右；且隨著k-mer長度的增加，N50值增大，L50值減小；若繼續(xù)增大k-mer值會(huì)進(jìn)一步優(yōu)化組裝參數(shù)。但由于SPAdes軟件進(jìn)行組裝的k-mer值最大為127，這可能是考慮到繼續(xù)加大k-mer值對測序深度要求較高，從而提高測序成本。

2 三代測序組裝結(jié)果及基因組結(jié)構(gòu)注釋

基于PacBio平臺(tái)的單分子測序共產(chǎn)生46 Gb大小的bam基因組測序文件，轉(zhuǎn)換為fasta文件后，根據(jù)二代組裝結(jié)果設(shè)置基因組大小為44 Mb進(jìn)行組裝。各組裝軟件組裝結(jié)果如表2所示，其中，選用Canu組裝的最優(yōu)結(jié)果通過quickmerge軟件與其他組裝結(jié)果進(jìn)行融合，并結(jié)合二代測序文件進(jìn)行基因組糾正（polish），在刪除重復(fù)contigs后的最終組裝結(jié)果見表3。

結(jié)合轉(zhuǎn)錄組測序文件進(jìn)行基因組結(jié)構(gòu)注釋共注釋出15684 個(gè)基因，并找到基因編碼區(qū)與氨基酸預(yù)測區(qū)，因結(jié)合轉(zhuǎn)錄組測序文件進(jìn)行結(jié)構(gòu)注釋，這些基因可能包含有可變剪切和重復(fù)注釋的結(jié)構(gòu)，會(huì)增加注釋出的基因數(shù)目。因此，進(jìn)行功能注釋和基因名稱注釋后需刪除重復(fù)的基因名；最終，共獲得6202 個(gè)基因注釋結(jié)果。

3 GO、KEGG、COG以及antiSMASH次級(jí)代謝注釋結(jié)果

對注釋出的6202 個(gè)基因進(jìn)行GO、KEGG與COG數(shù)據(jù)庫注釋，結(jié)果如圖1所示。其中，COG注釋結(jié)果顯示大部分基因與碳水化合物轉(zhuǎn)運(yùn)及代謝、氨基酸轉(zhuǎn)運(yùn)代謝、轉(zhuǎn)錄后修飾、RNA加工及修飾有關(guān)。KEGG注釋結(jié)果顯示大部分基因所處代謝通路與核糖體、過氧化物體、RNA轉(zhuǎn)運(yùn)有關(guān)。GO注釋結(jié)果顯示大部分基因與RNA、過氧化物體以及線粒體有關(guān)。最終，antiSMASH次級(jí)代謝物預(yù)測結(jié)果共發(fā)現(xiàn)4 個(gè)非核糖體肽合成酶（NRPS）基因簇、5 個(gè)一類聚酮合酶（pks）基因簇、3 個(gè)β-內(nèi)酯類合成基因簇以及7 個(gè)萜類合成基因簇，其中1 個(gè)一類pks基因簇和黑色素合成有關(guān)（相關(guān)性100%），1 個(gè)一類pks基因簇與黑麥酮酸類化合物合成有關(guān)（相關(guān)性18%）。

4 產(chǎn)黑色素短梗霉菌體透射電鏡結(jié)果

如圖2所示，在低產(chǎn)組中菌株具有較大細(xì)胞核（N）以及周圍存在的具有雙層膜結(jié)構(gòu)線粒體（M），在高產(chǎn)組中出現(xiàn)了類似乙醛酸循環(huán)體的圓形結(jié)構(gòu)，提示黑色素短梗霉中可能存在有乙醛酸循環(huán)途徑。

5 PMLA合成相關(guān)基因結(jié)構(gòu)預(yù)測

在對基因結(jié)構(gòu)與基因名注釋后，找到PCKA、MASY的蛋白質(zhì)序列，其晶體結(jié)構(gòu)預(yù)測結(jié)果如圖3、4所示，其中，PCKA、MASY與PDB上傳晶體結(jié)構(gòu)（PCKA：1YLH，MASY：3CUZ）進(jìn)行比對，結(jié)果顯示預(yù)測的蛋白結(jié)構(gòu)與PDB數(shù)據(jù)庫中序列比對結(jié)果基本一致，氨基酸比對結(jié)果中序列與其晶體結(jié)構(gòu)中小分子配體的結(jié)合位點(diǎn)也具有較高的一致性。同時(shí)蛋白的保守作用位點(diǎn)如PCKA中86R、140V、146G、287G、288D、289D，MASY中440C、275C、276G、277R、278W在比對結(jié)果中一致。

結(jié)合基因組及轉(zhuǎn)錄組測序結(jié)果可以組裝并注釋出質(zhì)量較高的基因組。本研究通過三代測序（Pacbio sequel II平臺(tái)）、二代測序（Illumina NovaSeq 6000）平臺(tái)對產(chǎn)黑色素短梗霉基因組進(jìn)行測序，因三代組裝需要預(yù)估基因組大小，因此，首先通過二代數(shù)據(jù)進(jìn)行基因組預(yù)組裝共得到44 Mbp基因組大小。之后考察了不同組裝軟件對產(chǎn)黑色素短梗霉基因組的組裝效果，在一般的默認(rèn)選項(xiàng)下，Canu軟件得到了較好的組裝效果。對該組裝結(jié)果通過二代數(shù)據(jù)修正并去重后得到包含26 個(gè)contigs、N50為2204220、GC值為50.09%、大小為42 Mb的較高質(zhì)量基因組組裝結(jié)果。通過轉(zhuǎn)錄組數(shù)據(jù)對該組裝結(jié)果進(jìn)行結(jié)構(gòu)注釋，共找到6202 個(gè)基因。產(chǎn)黑色短梗霉屬于出芽短梗霉的亞種，其在不同環(huán)境中也會(huì)呈現(xiàn)出酵母狀與菌絲體狀的不同形態(tài)。將基因組組裝與注釋結(jié)果與酵母和絲狀真菌的模式生物基因組進(jìn)行比較；其中釀酒酵母（Saccharomyces cerevisiae）基因組大小為12.15 Mb，編碼6016 個(gè)蛋白，GC含量為38.15%；構(gòu)巢曲霉基因組大小為30.30 Mb，編碼10008 個(gè)蛋白，GC含量為50.10%。因此，產(chǎn)黑色素短梗霉基因組組成更偏向構(gòu)巢曲霉，且本實(shí)驗(yàn)室先前通過無參轉(zhuǎn)錄組注釋發(fā)現(xiàn)很多與構(gòu)巢曲霉同源的基因。因構(gòu)巢曲霉是絲狀真菌的模式生物，該結(jié)果提示產(chǎn)黑色素短梗霉可能同樣適用構(gòu)巢曲霉的分子轉(zhuǎn)化方法。