《Genome》：The genome sequence of the a?aí berry (Euterpe oleracea Mart.) and RNA-Seq analysis of the fruit ripening
编辑推荐：
Euterpe oleracea Mart.（阿萨伊棕榈，夹竹桃科Arecaceae）是原生分布于巴西亚马逊东部湿地的多年生湿地棕榈，其产出的浆果状果实称为a?aí（阿萨伊果），是制备在当地具重大社会经济意义的饮品的原料。研究人员利用Oxford Nanopo
Euterpe oleracea Mart.（阿萨伊棕榈，夹竹桃科Arecaceae）是原生分布于巴西亚马逊东部湿地的多年生湿地棕榈，其产出的浆果状果实称为a?aí（阿萨伊果），是制备在当地具重大社会经济意义的饮品的原料。研究人员利用Oxford Nanopore长读长（long-read）测序平台完成了阿萨伊果基因组测序，并通过RNA-Seq分析了紫色与白色两个品种果实在成熟过程中的基因表达。组装所得基因组大小为3?066?969?163 bp，最长contig为430?834 bp，N50为54?646 bp，GC含量为49.97%。转录组分析鉴定到调控紫色阿萨伊果花青素（anthocyanin）合成的关键基因——例如黄酮类3′5′-羟化酶（flavonoid 3′5′-hydroxylase，F3′5′H）在开花后70~130天上调表达，促进飞燕草素-3-葡萄糖苷（delphinidin 3-glucoside）等紫色花青素的合成；与之相反，白色品种中花青素合成初期阶段相关酶基因普遍下调。本研究结果为理解阿萨伊果成熟过程的分子调控机制提供了重要依据，并报道了首个E. oleracea基因组。
论文解读：《The genome sequence of the a?aí berry (Euterpe oleracea Mart.) and RNA-Seq analysis of the fruit ripening》发表于《Genome》
研究背景与意义
Euterpe oleracea Mart.（阿萨伊棕榈，棕榈科Arecaceae）是亚马逊东部原生湿地多年生棕榈，其果实阿萨伊果（a?aí berry）是巴西北部具有重要社会经济价值的生物经济产品，可加工为饮品，富含α-生育酚（α-tocopherol）、膳食纤维、脂质、多酚（主要为花青素 anthocyanins）及矿物质，具抗炎、镇痛及抗氧化活性。除常见紫色果实外，部分单株产低花青素含量的白色阿萨伊果（white a?aí）。同属尚有E. precatoria（西部亚马逊及中美单干种）与E. edulis（大西洋森林种），二者果浆产量与感官品质低于E. oleracea。既往流式细胞术预估E. oleracea基因组约4.13 Gbp，前人仅完成叶绿体基因组及部分转录组研究，缺乏核基因组完整组装及果实花色差异的分子机制解析。为此，研究人员首次使用Oxford Nanopore长读长测序获得E. oleracea核基因组草图，并对紫、白两品种果实三个成熟期进行比较RNA-Seq，旨在阐明果实成熟尤其是花青素积累/缺失的转录调控基础。
主要关键技术方法
研究人员选取巴西Embrapa亚马孙东部分离保存的E. oleracea株系，根组织提取基因组DNA；果实取样为紫色品种（编号33.1）与白色品种（编号294-7）各一株，分别采集开花后70天（未熟）、130天（中转）、160天（成熟）三个时期的果实，各设3次生物学重复（共18个RNA样本）。基因组DNA经Oxford Nanopore PromethION 2 Solo平台测序，去除叶绿体污染读段后用Shasta进行de novo组装，QUAST评估组装指标，BUSCO检验基因完整性，Helixer预测基因模型，eggNOG-mapper与KEGG注释功能并与椰枣（Phoenix dactylifera）、椰子（Cocos nucifera）、油棕（Elaeis guineensis、E. oleifera）比较。RNA经mRNA富集建cDNA文库，Ion GeneStudio S5 Plus测序，Trimmomatic质控后用Trinity de novo组装转录本，TransDecoder预测ORF，DESeq2进行差异表达分析（DEGs，阈值|log2FC|≥1且padj≈<0.001），COG及KEGG对DEGs分类，重点分析花青素生物合成途径基因的表达模式。
研究结果
Whole genome sequencing（全基因组测序）
经Nanopore测序获12,202,098条读段，去除匹配E. oleracea叶绿体（NC_057603.1）的64,802条读段后用于组装。最终组装基因组大小为3,066,969,163 bp（约为流式细胞术预估值的75%），最长contig 430,834 bp，N50=54,646 bp，GC含量49.97%。Helixer注释得到41,703个基因。BUSCO评估表明基因集完整性良好。COG功能注释显示"L—复制与修复"（主要对应转座元件transposable elements, TEs）与"S—功能未知"占比最高；与其他四种棕榈比较发现E. oleracea与非洲油棕E. guineensis基因功能聚类最近，椰枣P. dactylifera因基因组较小且TE少而功能谱偏异。KEGG富集显示大基因组棕榈（E. oleracea、E. guineensis、E. oleifera）共有萜类骨架生物合成（terpenoid backbone biosynthesis）及硫循环相关基因，区别于椰子与椰枣。
Differential expression during fruit ripening（果实成熟过程差异表达分析）
样品相关性热图与PCA显示同品种同期重复聚为一类，白色阿萨伊果70 DAF与130 DAF转录谱相似度高于紫色品种同期。差异表达统计如下：白果70 vs 130 DAF检出76,541个转录本异构体（3,972上调、8,938下调），130 vs 160 DAF检出54,198个（3,188上调、4,766下调）；紫果70 vs 130 DAF检出62,577个（4,330上调、4,812下调），130 vs 160 DAF检出68,926个（1,208上调、4,457下调）。COG分类显示紫果70–130 DAF富集DNA复制修复（O）与蛋白翻译修饰（K），130–160 DAF富集碳水化合物代谢运输（G）与次生代谢产物合成（Q，含花青素）；白果70–130 DAF富集DNA复制修复（O）与氨基酸代谢（E），130–160 DAF富集细胞壁代谢（X）与核苷酸代谢（F），无次生代谢物合成相关明显变化。
花青素合成通路分析表明：白色阿萨伊果中CHS（chalcone synthase）、CHI（chalcone isomerase）、F3H（flavanone 3-hydroxylase）、DFR（dihydroflavonol 4-reductase）、ANS（anthocyanidin synthase）等早期及中期步骤酶编码基因均下调，UFGT（UDP-glucose:flavonoid 3-O-glucosyltransferase）表达稳定但因底物不足无法活化，故白色果不积累花青素。紫色阿萨伊果中F3′5′H（flavonoid 3′,5′-hydroxylase，F3′5′H）在70–130 DAF显著上调，促使二氢杨梅素（dihydromyricetin）生成增加，进而合成飞燕草素-3-葡萄糖苷（delphinidin 3-glucoside）等紫色花青素，是紫色果皮呈色关键；CHS/CHI在后期有下调但前期产物及可能转录后调控仍维持通路流量。
讨论与结论总结
研究人员指出此为E. oleracea首个核基因组版本，虽仅覆盖预估大小约75%，但已揭示棕榈科基因组普遍富含TE致COG中"L"类偏高，E. oleracea与E. guineensis功能相近。RNA-Seq比较表明紫、白品种成熟调控根本差异在于花青素途径基因的差异化表达——紫色品种F3′5′H特异性上调驱动花青素积累，白色品种该途径上游结构基因普遍抑制导致无色果实现象。此差异为阿萨伊果品质育种与分子标记开发提供靶基因。
结论（翻译）
研究人员报道了首版Euterpe oleracea基因组，大小约3.1 Gbp（约为预期基因组的75%）。通过该海量遗传信息鉴定到与41,703个注释基因关联的主要COG类别，最丰富的是"L—复制与修复"（主要由转座遗传元件代表）和"S—功能未知"。由于Euterpe属缺乏基因组数据，研究人员将E. oleracea与其他棕榈科物种比较，发现其与Elaeis属具高度功能相似性。注释基因数与转录组组装异构体数之差反映转录后调控的高度复杂性。在白、紫果不同成熟阶段鉴定到数千差异表达基因（DEGs）。黄酮类3′5′-羟化酶（F3′5′H）的上调对紫色阿萨伊果花青素合成增加似乎至关重要，反之白色品种花青素合成途径初始阶段编码基因普遍下调。本研究增进了对调控阿萨伊果成熟分子机制的理解，并提供重要亚马逊生物多样性象征物种E. oleracea的首个基因组序列。
