摘要：該研究以‘鐵觀音’茶樹品種的種子為試驗材料,采用轉(zhuǎn)錄組測序技術(shù)分析種子發(fā)育的3個時期(幼果期、膨大期、成熟期)的表達(dá)差異,探究茶樹種子油脂代謝的分子機制。結(jié)果表明:(1)經(jīng)轉(zhuǎn)錄組測序、組裝后共獲得30940581個clean reads,經(jīng)數(shù)據(jù)合并拼接最終得到36951條非冗余Unigene序列,其中28476個Unigene可得到功能注釋;在轉(zhuǎn)錄本中能夠被注釋到GO分類的Unigene有11201條(30.3%),KEGG分析發(fā)現(xiàn)共有17172個基因參與了127個代謝通路。(2)經(jīng)KEGG通路篩選出14條與脂肪酸代謝相關(guān)的通路,且隨著茶籽的發(fā)育,大部分脂肪酸調(diào)控途徑相關(guān)基因呈下調(diào)趨勢,其中上調(diào)基因數(shù)最多的有α-亞麻酸代謝途徑和脂肪酸降解途徑(有17個基因表達(dá)量上調(diào)),下調(diào)基因數(shù)最多的是甘油磷脂代謝途徑(有58個基因表達(dá)量下調(diào));在茶籽發(fā)育幼果期α-亞麻酸代謝途徑中表達(dá)量上調(diào)的基因數(shù)超過表達(dá)量下調(diào)的基因數(shù)。(3)研究發(fā)現(xiàn)茶籽脂肪酸合成相關(guān)的基因涉及14個脂類調(diào)控途徑,共409條差異基因;隨著茶樹種子發(fā)育到成熟期,上調(diào)的差異表達(dá)基因數(shù)量在減少,下調(diào)的差異表達(dá)基因數(shù)量增加,其中α-亞麻酸途徑中的基因PLA2G16、DAD1、pldA、FabF、FabI表達(dá)量上調(diào)顯著,隨后表達(dá)量下調(diào)。(4)qRT-PCR檢測結(jié)果表明,7個茶樹FAD和1個ACP差異表達(dá)基因的水平與轉(zhuǎn)錄組測序結(jié)果基本一致;隨著茶籽的發(fā)育,基因CsFAD7和Δ6-CsFAD從幼果期、果實膨大期至果實成熟期都為差異下調(diào)表達(dá),CsFAD2、CsFAD 6和Δ7-CsFAD為差異上調(diào)表達(dá),CsFAD 8、Δ8-CsFAD和CsACP在幼果期至果實膨大期差異上調(diào)表達(dá),在果實膨大期至果實成熟期差異下調(diào)表達(dá)。