《International Journal of Molecular Sciences》：Genome-Wide Identification and Bioinformatics Analysis of the FAD Gene Family in Walnut (Juglans regia L.)
编辑推荐：
脂肪酸去饱和酶（fatty acid desaturase，FAD）是植物中合成不饱和脂肪酸的核心催化酶，深刻影响植物生长、发育及对多种环境胁迫的适应性。核桃（Juglans regia L.）是重要的木本油料树种，其种仁富含不饱和脂肪酸。系统鉴定核桃FAD基
脂肪酸去饱和酶（fatty acid desaturase，FAD）是植物中合成不饱和脂肪酸的核心催化酶，深刻影响植物生长、发育及对多种环境胁迫的适应性。核桃（Juglans regia L.）是重要的木本油料树种，其种仁富含不饱和脂肪酸。系统鉴定核桃FAD基因家族并分析其功能，对揭示核桃不饱和脂肪酸代谢的分子机制具有重要意义。基于核桃全基因组数据，研究人员采用同源比对与隐马尔可夫模型（hidden Markov model，HMM）搜索方法鉴定JrFAD基因家族成员，随后利用多种生物信息学工具系统分析其结构特征、进化扩张机制、表达调控与功能。共鉴定出21个JrFAD基因家族成员，划分为5个亚家族，家族基因不均匀分布于9条染色体上。全基因组加倍（whole genome duplication，WGD）/片段重复是主要扩张方式，重复基因对经历强烈纯化选择。家族基因启动子序列富含响应光照、植物激素及多种胁迫的调控元件。表达模式分析显示，JrFAD3.1与JrFAD2.3在核桃种仁油脂快速积累期特异性高表达。本研究阐明了核桃FAD基因家族的组成与进化特征，为深入解析其在脂质代谢调控中的功能机制提供有用信息，并为高不饱和脂肪酸核桃品种遗传改良鉴定了潜在候选基因资源。
研究背景方面，脂肪酸及其衍生脂质是植物细胞重要的储能分子与生物膜结构组分，参与生长发育、信号转导与胁迫响应。不饱和脂肪酸含量与组成直接影响膜脂结构、膜流动性与稳定性，并参与植物对低温、高温、干旱、盐胁迫及病原侵染的适应。FAD是关键催化酶，按亚细胞定位分为可溶性FAD（以硬脂酰-ACP去饱和酶FAB2/SAD为代表，定位于质体基质，催化硬脂酸C18:0向油酸C18:1转化）与膜结合FAD（含FAD4、FAD2/FAD6、FAD3/FAD7/FAD8、ADS/SLD/DES等亚家族，分别定位于内质网或质体/叶绿体膜系统，催化油酸至亚油酸C18:2、亚油酸至α-亚麻酸C18:3等去饱和步骤）。核桃种仁脂肪含量60%–70%，油中约90%为不饱和脂肪酸（含亚油酸、α-亚麻酸等必需脂肪酸），但此前缺乏全基因组水平JrFAD家族系统鉴定。研究人员基于核桃栽培种‘Liaoning 1’的JRELIAO1参考基因组（550.89 Mb，BUSCO完整性98.9%，注释35294个基因），开展JrFAD家族全基因组鉴定与多维生物信息学分析，旨在阐明其组成、进化与表达特征，为核桃油脂品质形成的分子机制解析与遗传改良提供基础。
主要关键技术方法方面，研究人员以‘Liaoning 1’核桃为材料，采集种仁、青皮、壳组织及花后40、60、80、100、120天（DAF）五个主发育阶段种仁（另补50、67、74、87、94、110 DAF中间点），每阶段/组织三生物学重复，液氮速冻存于-80℃。基因组鉴定联合BLASTP（E≤1×10-5）与HMM（Pfam的FA_desaturase PF00487、FA_desaturase_2 PF03405、TMEM189 PF10520，E≤1×10-5）并经NCBI CDD校验；蛋白理化与亚细胞定位用ExPASy ProtParam与CELLO；系统发育用MUSCLE比对与MEGA11最大似然法（1000次bootstrap）；保守基序MEME（最多10个）、结构域Pfam/CDD、基因结构由GFF3提取经TBtools-II展示；染色体内分布与组内共线性用MCScanX，种间共线性比对拟南芥、水稻、大豆、毛果杨；Ka/Ks用Duplicate基因对计算；启动子2000 bp用PlantCARE预测顺式作用元件；表达分析基于实验室既往转录组（DESeq2归一化与差异分析、FPKM热图）；脂质组学用MTBE萃取-UHPLC-Orbitrap MS/MS负离子模式非靶向检测；加权基因共表达网络分析（WGCNA）用21586个表达基因软阈值18构建模块；基因-代谢物关联用Pearson相关与Benjamini-Hochberg校正（FDR<0.05）。
研究结果部分，2.1节鉴定与理化性质：经同源比对与HMM鉴定21个JrFAD成员（JrFAB2.1–5、JrFAD2.1–3、JrFAD3.1–3、JrFAD4、JrFAD6、JrFAD7、JrFAD8、JrADS、JrSLD1–3、JrDES1–2），蛋白长307–455 aa，分子量34.68–51.77 kDa，理论等电点6.10–9.95（18个偏碱），14个不稳定指数<40为稳定蛋白，多数亲水；CELLO预测FAB2在叶绿体、FAD2/3在内质网、SLD/DES/ADS在质膜，但与已知存在偏差需实验验证；二级结构以α-螺旋（34.37%–63.24%）与无规卷曲（31.31%–56.26%）为主。
2.2节系统发育：联用核桃、拟南芥、水稻、毛果杨、山核桃构建ML树，21个JrFAD归5大支（FAB2、FAD4、ADS/SLD/DES、FAD2/FAD6即ω-6、FAD3/FAD7/FAD8即ω-3），亚家族内核桃与毛果杨、山核桃聚得更近，亚家族节点bootstrap≥95%，支持分类保守性。
2.3节基序与结构：同亚家族基序/结构相似，FAB2均含Motif 9/6/1/7且具FA_desaturase_2域与EENRH、DEKRH组氨酸基序；FAD2单外显子含Motif 8/2/4与FA_desaturase+DUF3474；FAD6仅FA_desaturase域且含phase 0/1/2三种剪接相；FAD3/7/8多外显子含Motif 8/2/3（除JrFAD3.2）；SLD含Cyt-b5+FA_desaturase，DES含Sphingolipid Delta4-desaturase域，ADS仅FA_desaturase，FAD4独含Lipid_desat/TMEM189域。
2.4节染色体分布：21基因散于9条染色体，Chr01含4个最多，Chr14上JrFAD3.2与JrFAD3.3相邻呈串联重复，其余分散提示多机制扩张。
2.5节启动子元件：2000 bp上游富集光响应Box 4（66处）、G-box（62处），激素响应CGTCA/TGACG（MeJA）、ABRE（ABA）各约50处，胁迫响应ARE（厌氧，42处）、LTR（低温，17处），发育相关CAT-box（分生组织）、RY（种子特异）等，提示多信号调控。
2.6节重复与共线性：重复类型中WGD/片段重复10个为主，串联重复3个（JrFAB2.4、JrFAD3.2、JrFAD3.3），分散重复6个，单拷贝2个（JrFAD4、JrFAD6）；组内7对共线性（FAD3/7/8亚族最多），Ka/Ks 0.0958–0.1471均<1示纯化选择；种间与大豆（40对）、毛果杨（33对）共线性最强，拟南芥16对、水稻13对。
2.7节表达与关联：种仁发育中JrFAD3.1呈早低-中升-晚降，80–100 DAF（87 DAF达峰FPKM均值4809.35）快速积油期高表达且种仁偏好（果皮/壳<1/70）；JrFAD2.3在100 DAF达峰；JrFAB2.2、JrSLD1基础稳定表达；JrFAB2.1/JrFAD2.1在皮/壳偏高。WGCNA青色模块350基因与S5显著正相关（r=0.76），富集油体发生、脂质储存；JrFAD3.1/JrFAD2.3与JrABI3-1、JrWRI1-1及JrKASI/JrKASII/JrKAR、JrOle、Jr11S共表达。LC-MS特征中棕榈酸/C18:2在S3最高、硬脂酸S4升、油酸S1–S3高后S4降、C18:3在S2最高；105项基因-代谢物检验FDR<0.05者3项：JrSLD3与C18:3正相关（r=0.782），JrFAD2.1、JrFAD3.2与硬脂酸负相关（r=-0.762、-0.750）。
讨论部分总结：研究人员指出FAD家族大小具物种特异性（水稻20、大豆29、小麦68、核桃21），WGD/片段重复与纯化选择主导JrFAD扩张，ω-3亚族重复对扩张贡献大；启动子光/激素/胁迫元件富集与已知FAD参与膜脂代谢和胁迫适应一致但需实验验证；JrFAD3.1内核仁特异高表达与FAD3催化C18:2→C18:3生化功能吻合，JrFAD2.3峰稍晚，二者可能顺序参与多不饱和酸合成；非核仁组织高表达成员提示组织分化；关联分析仅为统计假设。结论部分翻译如下：本研究在核桃基因组中系统鉴定21个FAD基因家族成员并归5亚族，同亚族基序/域/结构保守而亚族间功能分化；基因散居9染色体，WGD/片段重复为主扩张力，重复对强纯化选择；启动子富光/激素/胁迫响应元件；JrFAD3.1与JrFAD2.3在种仁快速积油期特异高表达，为不饱和脂肪酸生物合成转录水平候选基因。研究为解析核桃FAD家族生物学功能及脂质代谢调控提供基础，未来需qRT-PCR、异源表达、基因编辑及群体自然变异与油质关联分析以支撑核桃油品改良性状改良。
