《Horticulturae》：Genome-Wide Identification and Multi-Project Transcriptome Compendium Profiling of the UDP-Glycosyltransferase Gene Family in Litchi (Litchi chinensis Sonn.)
编辑推荐：
荔枝（Litchi chinensis Sonn.）是一种高值亚热带果树，因其独特风味和营养特性而受到重视。UDP-糖基转移酶（UDP-glycosyltransferases, UGTs）催化多种小分子的糖基化，在果实发育、次生代谢和胁迫响应中发挥关键作用，
荔枝（Litchi chinensis Sonn.）是一种高值亚热带果树，因其独特风味和营养特性而受到重视。UDP-糖基转移酶（UDP-glycosyltransferases, UGTs）催化多种小分子的糖基化，在果实发育、次生代谢和胁迫响应中发挥关键作用，但该基因家族在荔枝中尚未被系统表征。本研究在15条荔枝染色体上鉴定出102个LcUGT基因，并将其划分为10个系统发育亚家族。重复分析表明，全基因组/片段重复和串联重复是主要驱动因素。通过整合跨越不同组织、发育阶段和胁迫条件的212个公开RNA测序（RNA-seq）样本，研究人员构建了LcUGT家族的综合表达图谱，并鉴定了在果实发育过程中具有差异化时空模式的候选基因。共表达网络分析进一步揭示了与假种皮和种子发育相关的组织特异性模块、枢纽UGT及候选上游转录因子，并且在另一个品种中进行的定量实时PCR（quantitative real-time PCR, qRT-PCR）检测显示，其表达趋势与转录组数据一致。该研究首次对荔枝UGT基因家族进行全基因组表征，为果实品质和耐逆性的分子育种提供了候选靶标。
研究背景方面，糖基化是小分子普遍存在的修饰，UDP-糖基转移酶（UDP-glycosyltransferases, UGTs；催化糖基转移的酶）参与果实发育、次生代谢和胁迫响应。荔枝（Litchi chinensis Sonn.）为高值亚热带果树，已有LcUFGT1和LcUGT87A等个别基因研究，但UGT家族缺乏全基因组系统鉴定，完整成员组成、进化关系及候选基因仍不清楚，因此研究人员开展全基因组鉴定和多项目转录组分析。该研究首次系统表征荔枝UGT家族，为功能解析和分子育种提供资源与候选靶标。
关键技术方法（不超过250字）：研究人员基于SCAULchv2.0蛋白组，以PF00201隐马尔可夫模型（hidden Markov model, HMM）和拟南芥UGT BLASTP鉴定成员，并用NCBI CDD/SMART验证结构域；用Muscle/FastTree构建最大似然树，MEME、TBtools分析基序与基因结构；MCScanX和KaKs_Calculator分析复制与Ka/Ks（非同义/同义替换率比值）；PlantCARE预测启动子；Blast2GO进行GO（Gene Ontology，基因本体）/KEGG（Kyoto Encyclopedia of Genes and Genomes，京都基因与基因组百科全书）富集；整合NCBI SRA的8个BioProject、212个RNA-seq样本（果实发育、组织、冷胁迫、甲基紫精（methyl viologen, MV）氧化胁迫、裂果、暗果皮病、转基因和野生比较）统一流程定量TPM（transcripts per million，每百万转录本）；WGCNA（weighted gene co-expression network analysis，加权基因共表达网络分析）/Cytoscape/PlantTFDB构建网络；对‘Feizixiao’假种皮47、55、62、69 dap（days after pollination，授粉后天数）样本进行qRT-PCR验证。
3.1 鉴定与理化性质：通过BLASTP和HMM鉴定102个LcUGT，蛋白长度99–1469 aa，分子量10.78–164.33 kDa，pI 4.70–9.23，不稳定性指数21.76–57.95，脂肪族指数68.94–117.74，GRAVY为?0.448–0.441，多数为负值而呈亲水性。
3.2 系统发育分类、染色体分布与复制分析：与79个AtUGT构建最大似然树后分为10组，E组30个最大，F和80/81各1个；15条染色体分布不均，13号最多15个，4号仅1个，并存在串联阵列；MCScanX将99个基因分为四类：近端重复41.4%、串联重复34.3%、全基因组复制（whole-genome duplication, WGD）/片段重复13.1%、分散重复11.1%；Ka/Ks为0.117–0.452，均小于1，表明纯化选择。
3.3 保守基序、基因结构与蛋白结构域：外显子1–12个，平均约2.5个；鉴定10个保守基序，Motif1和Motif2对应GTB型超家族且几乎存在于全部成员；注释到PLN02448、PLN02534、GT1_Gtf-like和YjiC等结构域，核心cl10013存在，基序和结构保守与分化支持亚家族分类。
3.4 启动子顺式调控元件：在2000 bp上游鉴定16类顺式元件，分为激素响应、胁迫响应、发育相关及光/昼夜响应；光响应和防御/胁迫响应元件最丰富，提示多种信号通路参与调控。
3.5 GO和KEGG功能富集：GO富集于木质素代谢、吲哚丁酸代谢、油菜素甾醇代谢及胁迫响应；分子功能以糖基转移、槲皮素3-O-葡萄糖基转移酶、UDP-糖基转移酶和槲皮素7-O-葡萄糖基转移酶为主；KEGG鉴定7条显著通路，玉米素生物合成和次生代谢物生物合成最显著，并含苯丙烷和硫代葡萄糖苷生物合成。
3.6 多项目转录组景观：212个样本的主成分分析前两主成分分别解释13.1%和9.6%方差；共检测31,896个表达基因，高表达基因数为12,699–19,228；组织特异性类别最高，果实发育次之；全基因组热图聚类对应主要生物分组。
3.7 LcUGT表达图谱：102个LcUGT层次聚类显示多数呈条件或组织偏好表达而非组成型，覆盖果实发育、MV氧化胁迫、暗果皮病、八组织及4 °C冷储等模块，显示家族功能分化。
3.8 候选LcUGT时空表达：6个基因在‘Huaizhi’（HZ）和‘Nuomici’（NM）四组织四阶段呈四类模式：果肉优势（LITCHI005027，HZ早峰、NM中峰）、早期发育（LITCHI015428）、种子优势中/晚峰（LITCHI021578、LITCHI009263）及多组织双功能（LITCHI010540、LITCHI007213）；品种间存在定量或时序差异。
3.9 共表达网络：WGCNA软阈值功率为14；blue模块2951个基因与aril正相关（r=0.87，p=2.7×10
?12），yellow模块1888个基因与seed正相关（r=0.94，p=6.4×10
?18）；枢纽基因为LcUGT16/LITCHI015428（blue早期aril；候选转录因子（transcription factor, TF）包括LOB、MYB-related、B3、bZIP、SWI/SNF-BAF60b、AUX/IAA）、LcUGT97/LITCHI005027（blue中晚aril；zf-HD、Tify、bHLH）、LcUGT102/LITCHI021578（yellow中期seed；HB-BELL、MADS-MIKC）、LcUGT90/LITCHI009263（yellow晚期seed；Trihelix、C2H2、ZF-HD、bZIP、HB-HD-ZIP、C2C2-Dof）。
3.10 qRT-PCR验证：对‘Feizixiao’假种皮四阶段检测6个候选基因，总体与TPM趋势一致；LITCHI007213和LITCHI021578渐增，LITCHI005027和LITCHI009263中峰，LITCHI010540早/晚较高而中间降低；LITCHI015428的qRT-PCR中峰与转录组早期高后降存在差异，原文归因于品种或季节效应；两数据均支持其早期活跃表达，相对表达采用2
?ΔΔCt法。
讨论部分总结：与葡萄181、百香果149和蓝莓361相比，荔枝102个UGT较少；10个组也少于拟南芥14–15组等；串联重复占34.3%，而蓝莓以WGD/片段重复为主；Ka/Ks均小于1，支持整基因水平纯化选择。GO涉及槲皮素3-O/7-O糖基转移酶、UGT78、脱落酸（abscisic acid, ABA）糖基转移酶UGT71B、水稻UGT85E1及蓝莓VcUGT160；启动子含G-box 210个、GT1 89个、ABA/茉莉酸甲酯（methyl jasmonate, MeJA）、低温响应元件（low-temperature-responsive element, LTR）和热激元件；11个‘Feizixiao’组织中5个基因τ≥0.80、18个0.60≤τ<0.80，LITCHI021578在seed的τ=0.84（seed 170.6 TPM、aril 129.7 TPM，其他组织≤13 TPM），LITCHI005027在aril达541.8 TPM；与柚84.1%广泛表达不同；HZ离区早峰、核桃约39%硬核期峰值、哈密瓜UGT73C2在20–30天转变；拟南芥UGT79B2/B3受CBF1调控，苹果MdUGT受NaCl、甘露醇、4 °C和ABA诱导；本数据中4 °C响应较MV明显，原文强调公共数据设计异质、结果为描述性，需正式差异统计；WGCNA提供组织模块和候选TF，为可检验假设。
研究结论翻译：本研究首次对荔枝UGT基因家族进行全基因组系统鉴定与表征，鉴定102个LcUGT并分为10个亚家族；近端和串联重复是家族扩张的主要驱动，重复对在整基因水平主要受纯化选择；通过整合212个公开RNA-seq样本构建泛转录组表达图谱，揭示LcUGT家族广泛的时空分化；四个枢纽基因为优先功能候选：LcUGT16（LITCHI015428）在早期aril优势表达，LcUGT97（LITCHI005027）呈aril优势表达并达541.8 TPM，LcUGT102（LITCHI021578）在seed强富集（τ=0.84）且aril高表达，LcUGT90（LITCHI009263）在晚期seed高表达且aril积累；共表达网络将这些枢纽置于aril和seed相关模块并鉴定候选上游TF；qRT-PCR在另一品种验证6个候选中的5个表达趋势；该研究为LcUGT功能解析提供基础资源，并为荔枝果实品质和耐逆性分子育种提供候选靶标。
