基于全基因組的紫芝細胞色素P450分析△

曾欣宜，朱英杰，李瀅，孫超*

(1.中國醫學科學院 北京協和醫學院 藥用植物研究所，北京 100193；2.中國中醫科學院 中藥研究所，北京 100700)

·基礎研究·

基于全基因組的紫芝細胞色素P450分析△

曾欣宜1，朱英杰2，李瀅1，孫超1*

(1.中國醫學科學院 北京協和醫學院 藥用植物研究所，北京 100193；2.中國中醫科學院 中藥研究所，北京 100700)

目的：紫芝是多孔菌科靈芝屬的重要傳統中藥。細胞色素P450是一類含血紅素的單加氧酶。本研究的目的是對紫芝基因組編碼的CYP450序列進行比對及系統進化分析。方法：發掘紫芝基因組編碼的所有P450基因，構建紫芝CYP450系統進化樹，比較靈芝屬3個已測序物種中P450基因家族的差異，并對紫芝中可能參與三萜生物合成的P450基因家族進行染色體定位和內含子相位分析。結果：紫芝基因組共編碼228個CYP450基因，其中包括9個假基因，在進化上分為3個分支，分布在41個家族中，其中CYP5359家族最大，包含有53條序列。紫芝所有的P450基因家族在其他兩個靈芝屬真菌中都存在。CYP512和CYP5144家族可能參與靈芝酸生物合成，其中一些家族成員可能是通過基因復制進化而來的。結論：紫芝P450基因家族的系統進化分析與進化機制研究為進一步揭示紫芝P450基因的功能與生理生化作用奠定了基礎。

紫芝；CYP450；靈芝三萜；進化

紫芝Ganodermasinense是名貴傳統中藥材，為擔子菌綱(Basidiomycetes)多孔菌目(Polyporales)多孔菌科(Polyporaceae)靈芝屬Ganodermae真菌[1]。現代藥理學研究表明紫芝具有提高免疫力，抗腫瘤[2-3]，調節血糖血脂[4-7]，延緩衰老[8]等多種功效。紫芝與赤芝G.lucidum作為中藥靈芝的基原物種，已被收錄于《中華人民共和國藥典》2010版[9]。

本課題組相繼完成了赤芝和紫芝基因組精細圖，為研究靈芝有效成分的生物合成及其他生物學問題奠定了基礎[10-11]。靈芝三萜是靈芝的主要活性成分之一，目前已從靈芝中分離得到150多種[12]。靈芝三萜是高度氧化的羊毛甾醇衍生物，推測有多個細胞色素P450單加氧酶(Cytochrome P450 monooxygenases，P450s)參與了靈芝三萜的生物合成過程。P450s是一類能夠在溫和的條件下把底物中反應惰性的碳氫鍵氧化的單加氧酶，是參與次生代謝的重要酶類之一，根據標準的細胞色素P450命名法，目前已鑒定大約19 000種P450s[13]。P450s主要參與的反應有：羥基化、環氧化、脫烷基化、碘化硫化、脫氨基化、脫硫化、脫鹵化和硝基還原作用[14]。本研究對紫芝編碼的所有P450基因進行了分類和進化分析，然后與靈芝屬其它兩個已測序真菌赤芝和Ganodermasp.中的P450基因進行了比較研究，并對可能參與靈芝酸合成的P450家族進行了預測。

1 材料與方法

1.1 序列來源

赤芝和紫芝細胞色素P450數據來源于本課題組紫芝基因組。

1.2 分析軟件及處理方法

1.2.1 系統進化分析 本文中用到的系統進化分析軟件有：分子進化遺傳分析軟件MEGA6.06，用鄰接法(Neighbor Joining，NJ)構建系統發育樹，自舉檢驗法Bootstrap(1000次)檢驗各分支的支持率[15]。

1.2.2 內含子相位分析 CYP512基因和CYP5144基因的內含子相位圖是利用in-house腳本對基因結構注釋文件進行分析得來。

2 結果與分析

2.1 紫芝CYP450分析

紫芝P450有228條序列，9條假P450序列，分布在41個家族中，其中CYP5359家族最為龐大，有53條序列均歸屬到該家族，其次為39條序列的CYP5150家族和31條序列的CYP512家族(圖1)。其中最長的3條CYP450蛋白序列分別是：CYP6005B1(1102aa)、CYP6005A1(1024aa)和CYP63A37(614aa)，而最短的條序列僅200多氨基酸：CYP5136D12P(221aa)和CYP5349A3(218aa)。對紫芝CYP450的系統進化分析，發現可以分成3枝(圖2)，其中最大的一枝共有132個CYP家族，包括：CYP5150(39)、CYP512(31)、CYP5035(25)、CYP5136(9)、CYP63(7)、CYP537(2)、CYP5140(2)、CYP5364(2)以及僅1條序列的CYP5151、CYP5138、CYP5137、CYP51、CYP5341、CYP5340和CYP5141。擁有88個CYP家族是第二大枝，這枝有著紫芝CYP家族中數目最多的CYP5359家族，共53條序列，其他家族數目均較少，為6條序列的CYP5037，以及4條序列的CYP5144家族和CYP5148家族，兩外還有3條序列的CYP5348家族和2條序列的CYP5058家族，以及僅1條序列的CYP502、CYP5351、CYP5158、CYP5065、CYP5360、CYP5152、CYP5349 、CYP5361、CYP5362、CYP5365、CYP53、CYP61、CYP642、CYP5347、CYP5366和CYP5139。第三枝所含CYP家族最少，共4個家族，包括7條序列：CYP505(3)、CYP6005(2)、CYP5141(1)和CYP5156(1)。在這枝中，CYP6005B1和CYP6005A1有著紫芝CYP450中最長的序列。

2.2 靈芝屬真菌P450分類及家族比較

對基因組水平的靈芝屬真菌P4 50分析發現(圖1)，靈芝屬的三個真菌的P450數目龐大，家族類別較多，雖比植物少，但是在真菌中，是較為豐富的。在這三種靈芝屬真菌中，相較于其他家族，CYP63、CYP512、CYP5035、CYP5037、CYP5136、CYP5139、CYP5144、CYP5150和CYP5359的數目是較多的，而CYP53、CYP61、CYP502、CYP642、CYP5065、CYP5137、CYP5138、CYP5151、CYP5152、CYP5156、CYP5158、CYP5347、CYP5349、CYP5351、CYP5360、CYP5362、CYP5365和CYP5366這些CYP家族在三種靈芝屬真菌中的數目均為一。其中紫芝CYP450家族最為龐大，并且紫芝CYP512、CYP5035、CYP5150和CYP5359均明顯多于其他兩個靈芝屬真菌(圖2)。而在G.sp中，CYP63、CYP5136、CYP5144和CYP5148是較多的。在赤芝中，僅CYP5356是其獨有的P450家族，而其他家族與另外兩種靈芝屬在數量上均較為一致，且并不多于其他兩種。

2.3 可能參與紫芝酸生物合成的CYP450分析

研究表明，CYP512和CYP5144在動物中對甾體激素、睪酮等激素的結構進行修飾，而在赤芝中CYP512基因與CYP5144基因與羊毛甾醇合酶(LSS)共表達[11]，因此預測這兩個P450家族可能參與紫芝三萜的生物合成。

圖1 紫芝CYP450家族系統進化樹

圖2 靈芝屬真菌P450家族分類比較(圖中G.sp數據來源于Khajamohiddin[16])

搜索CYP512和CYP5144基因后，對其染色體位置圖(圖3A)進行分析發現，這兩個家族的基因分布在2號、3號、4號、5號、7號、10號、11號和12號染色體，而1號染色體、6號染色體、8號染色體和9號染色體上均沒有分布。在5號染色體的0-500kb的位置共分布18個基因(圖3B)，并且分布位置是臨近的，因此推測可能存在重復基因的現象，對此構建了N-J樹和內含子相位圖(圖3C)。分析后可以發現，GS09913和GS09914可能互為重復基因，從內含子相位圖來看，極有可能為某段外顯子的“移動”所導致的基因重復，在N-J樹中這兩個基因進化關系也十分接近。GS06038和GS06037在2號染色體上位置接近，內含子相位也一致，由此可推測這兩個基因可能為重復基因。同樣地，GS04141和GS04140、GS07489和GS07486、GS03132和GS03141以及GS11489和GS11434也可能是互為重復基因。

圖3 CYP512基因和CYP5144基因分析

3 討論

紫芝CYP450共有228條P450序列，9條假P450s，分布在41個家族中，CYP5359家族(53 P450s)數目最多，其次為CYP5150家族(39 P450s)和CYP512家族(31 P450s)。對紫芝CYP450的系統進化可分為三枝，其中132個P450聚在數目最多的一枝，88個P450聚成的是第二枝，其中包括紫芝中數目最多的CYP5359家族，第三枝的數目較少，僅7個P450，聚在該枝的是紫芝中P450序列最長的幾條。而對三種靈芝屬真菌P450比較分析發現，在真菌中，其P450數目是比較龐大的，且CYP63、CYP512、CYP5035、CYP5037、CYP5136、CYP5139、CYP5144、CYP5150和CYP5359這些P450家族在三種靈芝中的數目都是較多的。值得關注的是，在三種靈芝屬CYP家族分析中，每種靈芝的基因組對P450家族數目的偏好是不同的。在紫芝中，CYP512、CYP5035、CYP5150和CYP5359這四個家族是其所偏好的，家族數目明顯多于其他兩種靈芝。而在G.sp中，CYP63、CYP5136、CYP5144和CYP5148是其偏好的。與另外兩種靈芝相比，赤芝并沒有其特有偏好的P450家族。

據報道，P450基因常與參與次生代謝的基因組成基因簇，共同調控靈芝次生代謝。在紫芝中，同樣地，預測CYP512家族和CYP5144家族極有可能與LSS(羊毛甾醇合酶)共表達，參與到紫芝三萜的生物合成中。進而對這兩個家族的基因進行分析，基因分布在染色體的位置圖表明其中有18個基因分布在5號染色體上，且位于100～600 Kb的位置。進一部分析基因的內含子相位和系統進化分析后發現，可能存在多組基因重復現象，分別為GS09913和GS09914、GS06038和GS06037、GS04141和GS04140、GS07489和GS07486、GS03132和GS03141以及GS11489和GS11434。

靈芝三萜是靈芝中主要的活性成分，它們是由共同的底物環狀骨架羊毛甾醇生成的，而CYP450是這一類底物的重要修飾酶家族。因此，通過對紫芝中CYP450的組成、分類及靈芝屬內的比較分析，對CYP450的功能驗證及解析紫芝中三萜類成分的生物合成奠定了基礎。為中藥靈芝資源的保護及合理開發利用，以及道地藥材的挖掘提供理論指導。

[1] 林志彬.靈芝的現代研究[M].北京：北京大學醫學出版社，2007.25-198.

[2] Liu J Q，Wang C F，Li Y，et al.Isolation and bioactivity evaluation of terpenoids from the medicinal fungusGanodermasinense[J].PlantaMed，2012，78(4)：368-376.

[3] Yue G G，Fung K P，Tse G M，et al.Comparative studies of various ganoderma species and their different parts with regard to their antitumor and immunomodulating activities in vitro[J].JAlternComplement，2006，12(8)：777-789.

[4] Sato N，Ma C M，Komatsu K，et al.Triterpene-Farnesyl hydroquinone conjugates fromGanodermasinense[J].Nat Prod，2009，72(5)：958-961.

[5] Yue G G，Chan B C，Han X Q，et al.Immunomodulatory Activities ofGanodermasinensePolysaccharides in Human Immune Cells[J].NutrCancer，2013，65(5)：765-774.

[6] Han X Q，Chan B C，Dong C，et al.Isolation,structure characterization，and immunomodulating activity of a hyperbranched polysaccharide from the fruiting bodies ofGanodermasinense[J].JAgricFoodChem，2012，60(17)：4276-4281.

[7] Wu H T，Lu F H，Su Y C，et al.Invivoandinvitroanti-tumor effects of fungal extracts[J].Molecules，2014，19(2)：2546-2556.

[8] Paterson R R M.Ganoderma-A therapeutic fungal biofactory[J].Phytochemistry，2006，67(18)：1985-2001

[9] 國家藥典委員會.中華人民共和國藥典：一部[S].北京：中國醫藥科技出版社，2010：205.

[10] Chen S L，Xu J，Liu C，et al.Genome sequence of the model medicinal mushroomGanodermalucidum[J].Nat.Commun.，2012(3)：913.

[11] Zhu Y J，Xu J，Sun C，et al.Chromosome-level genome maP provides insights into diverse defense mechanisms in the medicinal fungusGanodermasinense[J].Sci.Rep.，2015(5)：11087.

[12] 李曄，朱忠敏，姚渭溪，等.靈芝三萜類化合物的研究進展[J].中國中藥雜志，2012，37(2)：165-171.

[13] Nelson D R.Progress in tracing the evolutionary paths of cytochrome P450[J].ActaBiochim.Biophys.Sin.，2011，1814(1)：14-18.

[14] Sono M，Roach M P，Coulter E D，et al.Hemecontaining oxygenases[J].Chem.Rev.，1996，96(7)：2841-2888.

[15] Tamura K，Stecher G，Peterson D，et al.MEGA6：molecular evolutionary genetics analysis version 6.0.[J]Mol.Biol.Evol.，2013，30(12)：2725-2729.

[16] Khajamohiddin S，David R.，Robert R，et al.Genomewide annotation and comparative genomics of cytochrome P450 monooxygebases (P450s) in the polypore speciesBjerkanderaadusta，Gandermasp.andPhlebiabrevispora[J].Mycologia，2013，105(6)：1445-1455.

ResearchofCytochromeP450inGanodermasinenseBasedonAnalysisofGenomicData

ZENG Xinyi1，ZHUYingjie2，LIYing1，SUNChao1*

(1.InstituteofMedicinalPlantDevelopment，ChinaAcademyofMedicalSciences&PekingUnionMedicalCollege，Beijing100193，China；2.InstituteofChineseMateriaMedica，ChinaAcademyofChineseMedicalSciences，Beijing100700，China)

Objective：Ganodermasinense，belonging to the Polyporales，Ganodermae，is one of the most well-known medicinal mushrooms in the world.Cytochrome P450s are monooxygenases containing heme.Our study aims to understand the evolution of P450 families inG.sinense.Methods：We found all P450 genes from theG.sinensegenome and constructed their phylogenetic tree.The P450s fromG.sinensewere compared with those from the other two sequenced fungi from Ganoderma genus.We also analyzed the chromosome location and intron phases of the P450 genes possibly involved in the biosynthesis of triterpenoids.Results：TheG.sinensegenome encoded a total of 228 P450 genes，including 9 pseudogenes.All of them were divided into three clans and 41 families，among which，CYP5359 is the biggest family containing 53 members.All the P450 families inG.sinensealso exist in the other two sequenced Ganoderma species.CYP512 and CYP5144 families were supposed to be involved in the biosynthesis of triterpenoids and some of their members were generated by gene duplications.Conclusion：We analyzed the phylogenetics of P450s inG.sinenseand explored their possible evolutional mechanism，which will lay the foundation for further study on the function and physiological activities of these P450s.

Ganodermasinense；CYP450；triterpenoid；evolution

2016-04-06)

國家自然科學基金委面上項目(81273485)；國家自然科學基金委面上項目(81573704)

*

孫超，研究員，研究方向：功能基因組學；Tel：(010) 57833197，E-mail：csun@implad.ac.cn

10.13313/j.issn.1673-4890.2016.12.005