中棉所70分离群体农艺性状相关QTL定位

【目的】定位棉花产量相关性状的数量性状基因座(Quantitative trait locus,QTL)。【方法】以中棉所70的F_2分离群体为遗传作图群体,利用从14 820对简单序列重复(Simple sequence repeat,SSR)引物中筛选出的267对两亲本间的多态性引物检测F_2群体250个单株的标记基因型,利用Joinmap 4.0进行连锁分析,并通过WinQTLCart 2.5复合区间作图法对F_(2:3)群体的株高、单株结铃数和单株果枝数性状进行QTL定位。【结果】在F_2群体中共获得342个SSR标记位点,并构建了包括312个标记、35个连锁群,总长1 929.9 cM的遗传连锁图谱(标记间平均距离为9.2 cM,覆盖棉花基因组的43.4%)。经QTL定位,共检测到19个QTL,其中涉及株高的7个、单株果枝数4个、单株结铃数8个,这些QTL分布在8条染色体上,解释0.25%~11.28%的表型变异。【结论】这些与农艺性状相关的QTL有助于棉花产量分子标记辅助选择。     

陆地棉产量及农艺性状的SSR标记关联分析

【目的】挖掘与陆地棉产量和重要农艺性状密切关联的分子标记,为棉花分子标记辅助育种提供依据。【方法】本研究以147份陆地棉材料为供试群体,分别调查了7个环境下该群体的铃重和衣分数据,3个环境下的单株果枝数、单株铃数和株高数据。利用237对SSR标记对上述材料进行基因型检测,采用Tassel 2.1中的一般线性模型程序对上述5个性状进行关联分析。【结果】铃重、衣分、单株果枝数、单株铃数和株高的平均变异系数的变幅为7.28%～21.06%,平均值为13.47%。这表明供试群体具有较为丰富的表型多样性。237个标记在147份陆地棉中共检测到281个多态性位点,各位点上等位基因数的范围为2～6个,涉及690个等位基因,平均每个位点2.455 0个;多态性信息含量平均值为0.216 5;Structure2.2群体结构分析将147份供试材料划分为7个亚群;在3个及以上环境检测到的与铃重和衣分显著关联的标记位点分别有45个和1个,解释表型变异率在不同环境中最高分别达19.30%和11.58%(P0.01);在2个及以上环境检测到与单株果枝数、单株铃数和株高显著关联的标记位点分别有1、4和4个,表型变异解释率在不同环境中分别最高达9.96%、7.00%和5.75%,在BNL2448附近同时检测到与铃重、单株果枝数和株高显著关联的位点。【结论】通过关联分析挖掘了多个可重复检测到的与陆地棉产量及重要农艺性状关联的分子标记位点,此研究结果可为棉花高产育种分子辅助选择提供有益参考。     

陆地棉株型及生育期相关性状QTL定位

【目的】为了深入分析棉花株型及生育期相关性状的分子遗传机制,加速适机采棉花新品种分子标记辅助育种进程。【方法】通过构建包含413个单株的F2群体,结合高密度SNP(Single nucleotide polymorphism,单核苷酸多态性)遗传图谱,开展株高(Plant height,PH)、第一果枝节位(Node of the first fruiting branch,NFFB)及其高度(Height of NFFB,HNFFB)、第四果枝第一果节长度(First node length of the forth fruiting branch,FNLFFB)、第七果枝第一果节长度(First node length of the seventh fruiting branch,FNLSFB)等5个株型性状和开花期(Flowering time,FT)、花铃期(Flowering to boll-opening period,FBP)、全生育期(Whole growth period,WGP)等3个生育期性状的QTL(Quantitative trait loci,数量性状位点)定位研究。【结果】8个性状均呈现连续的双向超亲分布,性状之间存在广泛的正向相关性,PH与其他株型性状均为极显著正相关,NFFB与3个生育期相关性状均为显著或极显著正相关。共定位到36个加性QTL(Additive QTL,a QTL)位点,包括14个PH相关a QTL、6个NFFB相关a QTL、3个HNFFB相关a QTL、5个FNLFFB相关a QTL、3个FNLSFB相关a QTL、2个FT相关aQTL、2个FBP相关aQTL、1个WGP相关aQTL,单个aQTL贡献率为1.70%～10.38%。这些a QTL分布于20条染色体,每条染色体有1～4个a QTL。发现3个a QTL重叠区段,分别为A11染色体的qNFFB-A11-1与qWGP-A11-1、D3染色体的qHNFFB-D3-1与qPH-D3-1、D8染色体的q NFFB-D8-1与qFNLSFB-D8-1。检测到263个上位性QTL(Epistatic QTL,e QTL),单个e QTL的贡献率为1.17%～6.19%;19个a QTL与21个e QTL位置重叠;At基因组分布有17个a QTL和202个e QTL,Dt基因组分布有19个a QTL和61个e QTL。【结论】本研究为探索棉花株型的分子遗传机制奠定了研究基础,为机采棉分子标记辅助育种提供理论指导。     

小麦矮秆突变体je0098的遗传分析与其矮秆基因定位

倒伏易引发小麦严重减产,发掘和利用优异矮秆基因是培育高产抗倒伏小麦新品种的关键。本研究以京411(WT)及其经EMS诱变获得的产量相关性状优良的矮秆突变体je0098为试验材料,对其株高进行遗传分析,结合外显子捕获测序和遗传连锁分析定位矮秆基因。3年田间株高数据统计分析表明,je0098与WT相比株高降低15cm,组织细胞学观察结果显示,je0098与WT相比节间细胞长度缩短18%,暗示je0098的矮化是由于节间细胞长度变短所致;赤霉素敏感性分析表明, je0098为赤霉素敏感型矮秆突变体。利用WT和je0098杂交构建的由344个单株组成的F2分离群体,结合F2:3家系表型数据,选取矮秆纯合和高秆单株构建混池,对两亲本和子代混池分别进行外显子捕获测序,在2D染色体上定位到一个具有降秆效应的数量性状位点(QTL)。结合全基因组重测序所得SNP位点,在2D染色体开发了6个KASP分子标记,对F2单株进行基因分型。利用QTL IciMapping作图软件构建遗传连锁图谱,结合3年田间表型数据,将矮秆基因定位在20.77～28.84 Mb区间内,遗传距离为11.48 cM。本研究结果为突变...     

利用陆地棉MAGIC群体定位产量、生育期和株高性状的QTL

棉花的产量、生育期和株高是重要的农艺性状, 决定着棉花的经济效益和生产方式。为了研究这些性状的遗传基础, 利用8个亲本构建的包含960个株系的陆地棉MAGIC (multi-parent advanced generation inter-cross)群体和SSR标记对这些性状进行关联分析。对产量、生育期和株高共8个相关性状进行了3年3个地点共5个环境的田间试验。经过最佳线性无偏预测(BLUP)综合多环境的表型数据, 发现MAGIC群体比亲本有更丰富的表型变异。8个性状的广义遗传力(H ²)变化范围为0.17~0.71。结合284个SSR标记基因型数据, 利用混合线性模型对产量、生育期和株高相关性状进行关联分析, 分别检测到51、27和9个显著关联的位点, 这些位点都表现出微效性, 表明该陆地棉MAGIC群体在性状位点的挖掘方面具有高效性。检测到20个标记位点或区间控制多个性状, 还发现单株有效铃数、单铃皮棉重、第一果枝节位和株高存在染色体热点区域, 对多性状综合研究或单性状深入挖掘具有重要价值。本研究为后续深入利用MAGIC群体进行遗传研究提供参考, 一些表型优良的材料和关联到的位点为育种改良奠定了基础。     

机采棉主要农艺性状与密度相关性分析

为研究棉花适宜机采农艺性状与密度的相关性,明确通过密度塑造适宜机采株型的方法,采用大田试验,以K836为试验材料,设计3×10 ⁴、6×10 ⁴和9×10 ⁴株/hm ² 3个密度处理,研究密度对机采棉子棉产量的影响及其与机采棉主要农艺性状的相关性分析。结果表明,在本研究密度范围内,棉花株高、果枝长度、果枝第一节位长度、果枝节数、单株果枝数和单株干物质量随着密度的增加而下降,密度对第一果枝节位和果枝夹角没有显著影响。与低密度（3×10 ⁴株/hm ²）相比,高密度处理（9×10 ⁴株/hm ²）棉花单铃重降低,衣分提高,单株铃数降低,总铃数增加,密度对子棉产量没有显著影响。除第一果枝高度外,株高、第一果枝节位、果枝长度、果枝第一节位长度、果枝节数、果枝夹角、单株果枝数与单株干物质量和单株铃数呈正相关,与单铃重和子棉产量呈负相关,其中果枝长度和单株果枝数呈极显著正相关。因此,适当增加种植密度使棉花株高降低,果枝变短,株型更为紧凑,可以通过密度塑造适合机械采收的株型,冀南地区高密度（9×10 ⁴株/hm ²）处理棉花株型符合机采要求。     

基于两个陆地棉低世代群体定位纤维品质相关QTL

【目的】定位棉花纤维品质性状相关的数量性状位点(Quantitative trait locus,QTL)。【方法】以陆地棉高强纤维品系中棉所679和纤维品质一般的农垦5号为亲本构建包含200个单株的F2群体及对应的F2:3家系群体,对2个群体的纤维长度、断裂比强度等5个纤维品质性状进行检测。用6 688对简单重复序列(Simple sequence repeat, SSR)引物在双亲间筛选,得到149对多态性引物,以F2为作图群体,使用QTL IciMapping软件进行连锁图谱构建,并对F2及F2:3群体进行QTL定位。【结果】根据F2群体基因型信息构建了1张包含119个标记、28个连锁群、总长为1 173.5 cM(centiMorgan)的遗传连锁图谱。分别在F_2、F2:3群体中检测到9个和11个与纤维品质性状相关的QTLs,这些QTLs分布在11个连锁群上。其中F2群体的qFL-D11-1、q BT-D11-1与F2:3群体的qFL-D11-1、q MIC-D11-1均定位在标记DPL0062与HAU0423之间,推测这些位点可能是控制纤维品质性状的重要QTL。【结论】利用多个群体进行QTL定位有益于发现稳定的QTL位点,控制纤维品质性状的基因可能成簇存在,为挖掘纤维品质性状相关基因及分子标记辅助育种奠定基础。     

陆海高代回交群体抗黄萎病QTL定位

【目的】定位棉花抗黄萎病数量性状位点(Quantitative trait loci, QTL)。【方法】以海7124和TM-1配制抗感组合F1,再以鲁棉研28为轮回亲本构建的137个BC4F1家系为作图群体,筛选出多态性重复序列(Simple sequence repeat, SSR)标记,并与已发表的整合高密度遗传连锁图谱相比对,构建遗传图谱。采用复合区间作图法(Composite interval mapping,CIM)进行大田和病圃两个环境下抗黄萎病QTL定位。【结果】216个多态性SSR位点分布在26条染色体上,可覆盖棉花基因组3 380 cM(centi Morgan),标记间平均距离15.77 cM。定位到6个QTLs,分布在6条染色体上,可解释表型变异8.56%~20.26%,其中5个QTLs与前人研究结果相一致,在第1染色体上新定位到一个QTL。本研究可为分子标记辅助选择抗病育种提供帮助。【结论】定位到6个黄萎病相关QTLs,其中1个是在第1染色体上新发现的QTL。     

基于两个陆地棉低世代群体定位纤维品质相关QTL

【目的】定位棉花纤维品质性状相关的数量性状位点(Quantitative trait locus,QTL)。【方法】以陆地棉高强纤维品系中棉所679和纤维品质一般的农垦5号为亲本构建包含200个单株的F2群体及对应的F2:3家系群体,对2个群体的纤维长度、断裂比强度等5个纤维品质性状进行检测。用6 688对简单重复序列(Simple sequence repeat, SSR)引物在双亲间筛选,得到149对多态性引物,以F2为作图群体,使用QTL IciMapping软件进行连锁图谱构建,并对F2及F2:3群体进行QTL定位。【结果】根据F2群体基因型信息构建了1张包含119个标记、28个连锁群、总长为1 173.5 cM(centiMorgan)的遗传连锁图谱。分别在F_2、F2:3群体中检测到9个和11个与纤维品质性状相关的QTLs,这些QTLs分布在11个连锁群上。其中F2群体的qFL-D11-1、q BT-D11-1与F2:3群体的qFL-D11-1、q MIC-D11-1均定位在标记DPL0062与HAU0423之间,推测这些位点可能是控制纤维品质性状的重要QTL。【结论】利用多个群体进行QTL定位有益于发现稳定的QTL位点,控制纤维品质性状的基因可能成簇存在,为挖掘纤维品质性状相关基因及分子标记辅助育种奠定基础。     

基于BSA-seq的甜瓜抗枯萎病相关基因的定位与候选基因的注释

为挖掘甜瓜抗枯萎病相关基因,本研究以高抗枯萎病自交系YN为母本,高感枯萎病自交系CG为父本构建F₂群体,基于F₂群体分别构建极端高抗混池和高感混池,对子代混池和双亲池开展30×覆盖深度的全基因组重测序,通过计算2个子代池的SNP-index,比较SNP-index在高抗池与高感池染色体各区段的差异,定位与抗病性的关联区域,根据基因注释信息,预测候选基因。结果表明4个样本获得的SNP位点和InDel分别在91万～250万个之间和17万～45万个之间。基于SNP-index关联分析,在99%的置信区间内将候选区域定位到甜瓜染色体Chr.9的1个区间,总长度为1.37 Mb,位于1～1 370 000 bp区间,包含197个基因,SNP-inDel位点数7 442个。通过对候选区域内的基因进行KEGG、GO数据库分析发现,候选基因分别被注释到18个生物学进程、9个细胞组分和6个分子功能中;主要参与氨基酸生物合成、氨基酸代谢、酯类代谢、糖类代谢等,经分析在候选区域内与抗病性相关的基因共有22个,3个与F-box基因相关,3个水解蛋白相关基因,4个乙烯...     

牛大力漆酶基因CsLAC17的克隆与载体构建

为研究漆酶在牛大力生长发育过程中的生物学功能,本实验以牛大力叶片为材料,从反转录PCR获得的cDNA中扩增出漆酶基因CsLAC17全长。序列分析表明,CsLAC17 c DNA全长为2 010 bp,开放阅读框大小为1 755 bp,编码一个由584个氨基酸组成的蛋白质。结构域分析表明,Cs LAC17蛋白的保守结构域具有漆酶典型结构域的特征——铜离子结合域(Cu-oxidase和Cu-oxidase-2)。同源序列分析表明,Cs LAC17蛋白序列与绿豆(Vigna radiata var. radiata)和野生大豆(Glycine soja)同源性为88%、藜豆(Mucuna pruriens)87%、蒺藜苜蓿(Medicago truncatula) 85%。组织特异性表达分析显示,CsLAC17在茎中表达量最高,叶中表达量次之,根中较少。此外,进一步构建了pBI121-CsLAC17过表达载体并转入农杆菌。本研究为日后CsLAC17基因的功能验证以及为牛大力开展分子生物学研究提供帮助。     

基于高通量测序的棉花InDel标记开发及其应用

【目的】获得可以用于棉花品种鉴定和纯度检测的二态性InDel标记,提高棉花种子的检验精确度和效率,为棉花的分子育种发挥作用。【方法】基于来源不同的121份棉花全基因组信息,根据高多态性信息含量筛选多态性高的InDel位点,开发二态性InDel标记,并在我国66个棉花品种的遗传距离分析和聚类分析中进行应用。【结果】基于121份棉花的二代测序数据获得10967个InDel位点,合成了85对InDel引物,选择其中有效引物64对,其中At亚组的特异引物35对、Dt亚组特异引物29对,At、Dt亚组染色体的平均最小等位频率分别为0.45、0.32,多态信息含量分别为0.49、0.40;所用的66个棉花品种的遗传距离范围是0.04~0.65cM(厘摩),平均为0.39cM,遗传距离最大的2个品种是泗棉3号和中棉所36,遗传距离最小的是徐棉18和徐杂3号。【结论】开发的64个棉花二态性InDel标记能有效地通过揭示品种间的遗传距离反映品种之间的亲缘关系,区分来源不同的棉花品种,具有一定的理论意义和应用价值。     

植物生长促进剂对赣北移栽棉生长发育及产量的影响

【目的】明确叶面喷施6种植物生长促进剂对赣北移栽棉生长发育的影响,筛选出对棉花增产效果好的类型在赣北植棉区棉花生产上推广应用。【方法】在赣北育苗移栽植棉模式下,利用萘乙酸、赤霉素、细胞分裂素、芸苔素、复硝酚钠和胺鲜酯于棉花苗期(三叶期)、现蕾期、开花期和盛花期喷施,研究其对棉花生长发育的影响。【结果】细胞分裂素和复硝酚钠在棉花前中期对生长的促进作用明显,且复硝酚钠表现盛花结铃期成铃速度快;胺鲜酯对棉花生长的促进作用前后期保持较好的态势,单株叶面积系数、单株鲜物质质量和干物质质量均最大;芸苔素在棉花受到不良环境的伤害后可增强棉株的抗逆能力,明显提高铃重和产量;萘乙酸虽然对棉花前中期生长发育影响小,但最终增产效果还较好;赤霉素对棉花后期生长不利,棉花表现个体小、铃重略降低;芸苔素+胺鲜酯不及各自单用效果好。【结论】6种植物生长促进剂在赣北育苗移栽植棉模式下,于棉花苗期、现蕾期、开花期和盛花期叶面喷施,对棉花生长发育均有促进作用,可提高籽棉产量2.43%～10.04%,建议棉花苗期和蕾期叶面喷施细胞分裂素或复硝酚钠,花铃期叶面喷施芸苔素或胺鲜酯。     

棉花CLE多肽家族的全基因组鉴定与生物信息学分析

【目的】CLE(CLAVATA3/Embryo surrounding region-related)多肽是一类小分子分泌蛋白,在植物生长发育过程中,对于维持分生组织细胞增殖与分化间的平衡起到重要的信号转导和调控作用。CLE家族是迄今为止报道的最大的植物多肽分子家族,也是近十年来研究最为热门的植物多肽激素。尽管CLE基因家族在拟南芥和水稻等模式植物中取得较多的研究进展,但在棉花中有关CLE多肽基因家族的研究尚属空白。通过对CLE多肽家族的鉴定及演化研究,可以为进一步挖掘棉花中的功能多肽及其信号转导研究提供一定的理论基础。【方法】本文以亚洲棉(Gossypium arboreum L.)、雷蒙德氏棉(G. raimondii Ulbrich)、陆地棉(G. hirsutum L.)和海岛棉(G. barbadense L.)为研究对象,利用生物信息学方法和已鉴定的拟南芥CLE基因序列,在棉花基因组中进行同源比对,并对4个棉种中的CLE家族成员进行全基因组鉴定和分析。【结果】共得到148个棉花CLE候选基因。陆地棉和海岛棉中CLE基因数量都分别是亚洲棉和雷蒙德氏棉的2倍。聚类分析将所有CLE基因分为5个亚组。Ka/Ks分析表明大部分基因都经历了负选择作用,有11对基因经历了显著的正选择作用。转录组数据分析发现,除海岛棉中的GbCLE39、GbCLE13、GbCLE43,陆地棉中GhCLE34、GhCLE9以及亚组棉中的GaCLE4外,大部分CLE基因在棉花组织中的表达量较低,这与多肽在植物体内含量较低相一致。CLE核心保守基序分析发现了12个棉花特有的多肽,其中有2个多肽来源于受到强烈正选择作用的基因,因此在后续研究中可以对其进行进一步深入分析。【结论】本文利用生物信息学、比较基因组学等方法,首次对棉花的CLE基因家族进行了鉴定和分析,对于进一步挖掘棉花中的功能多肽及其信号调控网络研究提供了一定的理论基础。     

菠萝AcSAP转录因子对非生物胁迫和生物胁迫的应答响应

STERILE APETALA (SAP)是调节花序、花和胚珠发育,调控分生组织细胞的增殖和器官大小的多功能基因。本研究利用生物信息学的方法对菠萝SAP转录因子进行序列分析,并通过qRT-PCR技术研究了非生物胁迫和生物胁迫对菠萝SAP基因表达的影响。结果表明,菠萝SAP蛋白为稳定疏水酸性蛋白且含有3个WD40重复区域,蛋白质二级结构显示,菠萝SAP蛋白主要结构元件为延伸链和无规则卷曲;qRT-PCR分析显示,逆境胁迫下,AcSAP的表达量与对照差异显著。在H2O2、NaCl、SA、ABA、Eth、低温(4℃)和病菌侵染胁迫下,AcSAP的表达量显著升高。研究表明,逆境胁迫能使AcSAP基因的表达受到影响,进而直接或间接影响植物生长发育,为今后菠萝的抗逆研究和分子育种提供了理论依据。     

复合性状转基因棉花对土壤线虫多样性的影响

为了评价复合性状转基因棉花种植对土壤线虫群落结构的影响,分析了大田条件下分别种植3种复合性状转基因棉花s GK321、转Cry1Ac+Cry2Ab基因棉(双Bt抗虫棉)和转Cry1Ac+Epsps基因棉(抗虫抗除草剂棉)及常规棉石远321的土壤线虫数量、营养类群组成和生态指数。结果表明,与石远321相比,种植3种复合性状转基因棉的土壤线虫数量无显著差异。土壤中共鉴定出20科、36属土壤线虫,包括10属植物寄生线虫、10属食细菌线虫、5属食真菌线虫和11属杂食/捕食性线虫。其中食细菌线虫为优势营养类群,真头叶属(Eucephalobus)、螺旋属(Helicotylenchus)、板唇属(Chiloplacus)、绕线属(Plectus)和原杆属(Protorhabditis)为优势属。3种复合性状转基因棉和石远321种植下土壤线虫群落Shannon多样性指数(H)、优势度指数(λ)、丰富度指数(SR)和成熟度指数(MI、PPI)均无显著差异,只有双Bt抗虫棉土壤线虫群落Pielou均匀度指数有显著差异。说明大田种植3种复合性状转基因棉对土壤线虫群落结构和种类组成的影响不显著。     

基于单拷贝SNP标记的棉花杂交种纯度高通量检测技术

利用有代表性的材料进行SNP位点的全基因组扫描分析与综合评估,基于KASP技术开发1套适用于我国棉花杂交种纯度高通量检测的核心SNP组合。从63K的棉花全基因组芯片中筛选获得具有单拷贝特性的SNP标记分别为5474个(中棉所TM-1基因组版本)和1850个(南京农大TM-1基因组版本)。根据芯片扫描分析结果,权衡考虑位点多态性、分型效果、纯合率与杂合率等因素,最终从每条染色体上优选1个杂交种杂合率高且分型效果理想的核心SNP位点,合计26个。采用KrakenTM软件将SNP位点转化成KASP标记,利用SNPline平台进行SNP分型检测,实现了对大量样品的高通量基因分型,尤其适用于品种纯度快速检测,为SNP标记技术在棉花品种鉴定及指纹数据库构建等方面的应用奠定基础。     

棉花磷胁迫响应基因GhWRKY6克隆及功能分析

【目的】提高棉花在低磷胁迫下的抗性,挖掘棉花磷胁迫相关功能基因。【方法】用生物信息学方法分析Gh WRKY6基因的结构特征和进化关系;构建基因超表达载体转化拟南芥,分析转基因拟南芥在低磷胁迫下的抗性;利用荧光定量PCR(Real-time quantitative polymerase chain reaction, qRT-PCR)分析拟南芥中磷信号途径Marker基因表达量的变化。【结果】以陆地棉cDNA为模板克隆得到Gh WRKY6基因,生物信息学分析表明该基因序列含有一个WRKY保守结构域,是WRKY家族转录因子。在低磷胁迫下,转基因拟南芥与野生型相比,种子萌发速度、主根长度、侧根数目、生物量均低于野生型,而在正常生长条件下两者并无差别。qRT-PCR分析表明GhWRKY6能够影响关键的磷转运蛋白基因的表达。【结论】GhWRKY6作为一个负调控因子参与到植物低磷胁迫响应信号途径中。     

沙棘果实发育阶段苹果酸代谢关键基因的表达分析

为探讨沙棘果实苹果酸积累与相关基因表达间的关系,以偏低酸品种‘橙色’和偏高酸品种‘芬兰’不同发育期(D1～D5期)的果实为试材,测定可溶性固形物和苹果酸含量,采用实时荧光定量PCR分析苹果酸代谢关键基因的表达模式。结果表明,第一,可溶性固形物含量在果实发育过程中呈上升趋势,成熟时(D5期)最高,‘橙色’比‘芬兰’高;苹果酸含量的变化与之相反,果实发育阶段呈下降趋势,成熟时最低,‘橙色’比‘芬兰’低。第二,果实发育过程中,与苹果酸合成相关的PEPC基因表达水平呈下降趋势,且在成熟期(D4～D5期)‘橙色’比‘芬兰’下降迅速;cytMDH基因在‘芬兰’果实的表达整体上呈先升后降趋势,转色期(D3期)为转折点,在‘橙色’各时期的表达水平变化较小。第三,NADP-ME基因在整个果实发育阶段的表达呈上升趋势,‘橙色’比‘芬兰’明显;PEPCK基因的表达在两个品种间存在显著差异,在‘芬兰’中的表达水平变化较小;在‘橙色’中呈现升-降-升-降的趋势,转色期后迅速表达至峰值,与NADP-ME共同作用促进苹果酸在成熟期含量迅速降低,可溶性固形物含量迅速提高。本研究结果为进一步深入探讨沙棘苹果酸代谢调控基因的功能及表达提供科学依据。