甜瓜幼果果皮颜色基因GR的精细定位

2021-08-20许昕阳沈佳张跃建李国景牛晓伟寿伟松

中国农业科学 2021年15期

许昕阳，沈佳，张跃建，李国景，牛晓伟，寿伟松

甜瓜幼果果皮颜色基因的精细定位

许昕阳，沈佳，张跃建，李国景，牛晓伟，寿伟松

浙江省农业科学院蔬菜研究所，杭州 310021

【】探究甜瓜幼果果皮颜色性状的遗传规律，精细定位目标性状基因，加深对甜瓜发育过程中果皮颜色转变的认知，为开展甜瓜果皮颜色的分子设计育种奠定基础。以幼果深绿皮的薄皮甜瓜纯系‘MR-1’和幼果浅绿皮的厚皮甜瓜纯系‘LGR’为亲本，构建F1正反交群体；以及利用F1与浅绿皮亲本‘LGR’杂交构建BC1F1回交群体，对甜瓜幼果果皮颜色基因（）进行遗传分析。选取BC1F1群体中深绿皮和浅绿皮单株各20株，混池其DNA进行BSA-seq以获取初定位区间。基于‘MR-1’和‘LGR’两亲本的重测序数据，开发初定位区段内特异性较好的分子标记，鉴定筛选扩大群体（BC1F1和F2）中的重组交换单株，验证和缩小定位区间，实现精细定位。将两亲本定位区段内注释基因的编码区进行测序以确定候选基因和关键变异位点。通过调查BC1F1回交群体中幼果果皮颜色和成熟果果皮颜色，利用相关性分析探究果皮颜色转变在甜瓜发育过程中的内在联系。通过分析F1群体果皮颜色发现所有F1单株幼果都表现为深绿皮。另外，BC1F1群体单株幼果果皮颜色会发生分离，其中深绿皮单株数﹕浅绿皮单株数约等于1﹕1，以及F2群体中深绿皮植株与浅绿皮植株的分离比为3﹕1。这些分离比都符合孟德尔遗传定律，表明幼果果皮颜色是受单个核基因控制的质量性状，并且深绿对浅绿为显性。通过BSA-seq分析将基因初步定位于4号染色体长臂，物理距离为1.8 Mb的范围内。利用开发的分子标记在扩大的定位群体中共筛选到24个重组交换单株。经过后代基因型和表型验证，最终将精细定位在标记4-102和4-81之间约17.7 kb的范围内，区段内共包含4个注释基因。经测序分析发现一个编码GLKs类转录因子CmAPRR2的基因在亲本‘MR-1’和‘LGR’中存在多处变异，其中有3处发生了同义突变，1处错义突变和1处无义突变。无义突变出现在的编码区第856位碱基处（由G变成T），导致亲本‘LGR’中蛋白翻译提前终止，其Myb-DNA结合结构域大部分缺失，推测基因（）即为影响甜瓜幼果果皮颜色的基因，而第856位的单碱基替换造成的无义突变即为关键变异位点。此外，BC1F1回交群体单株的表型调查结果显示幼果与成熟果的果皮颜色之间存在显著相关性。甜瓜幼果果皮颜色（深绿/浅绿）性状为质量性状，受单个核基因控制。通过遗传定位手段推断（）为最有可能影响甜瓜幼果果皮颜色的候选基因。

甜瓜；幼果果皮颜色；基因定位；；（）

0 引言

【研究意义】甜瓜（L.）是我国重要的园艺和经济作物。果皮颜色作为甜瓜最直观的外观品质性状，直接影响其商品价值和消费者的选择，开展甜瓜果皮颜色性状的遗传研究对甜瓜育种有重要意义。【前人研究进展】瓜类植物的果皮颜色调节开始于果实发育的早期，在幼果期表现为绿色，反映的是叶绿素含量[1]。随着果实发育的进行，甜瓜外果皮颜色在成熟过程中会发生较大转变[2]。成熟甜瓜果皮颜色丰富多样，主要有白色、黄色、橙色、绿色以及多种颜色混杂或者斑纹类型。早期对甜瓜果皮颜色是由单基因或者多基因控制存在较大的争议，但该争议多基于不同品种间组合的差异[3-4]。KUBICKI[5]认为幼果果皮白色对绿色为显性，并且是受单基因遗传控制。PEREIRA等[6]通过构建‘Védrantais’和‘Piel de Sapo’的高密度GBS遗传图谱发现浅皮为显性性状，并且将该基因定位在7号染色体，该位置与相近，推测为同一个基因。BURGER等[7]利用两个杂交群体：深皮‘Dulce’和浅皮‘TAM-Dew’，深绿‘Krymka’和浅绿‘Eshkolit Ha’Amaqim’探究甜瓜幼果果皮颜色遗传模式时发现浅皮亲本受隐性单基因控制，并且该基因与之前报道的浅皮显性基因不等位。OREN等[8]利用‘Tam Dew’和‘Dulce’的重组自交系以及‘Dulce’和‘Noy Amid’的F3:4家系将一个影响甜瓜幼果果皮颜色的基因定位在4号染色体上290 kb的区段内，后通过差异性分析确定基因是影响该幼果果皮颜色的基因。在成熟果中，HUGHES[9]发现‘Honeydew’品种的白色相对于‘Smiths’ Perfect cantaloupe’品种的绿色为隐性。TADMOR等[1]在黄皮‘Noy Amid’和绿皮‘TVT’杂交的F2分离群体中发现柚皮素查尔酮含量的分离比符合孟德尔3﹕1的遗传定律，证实了该色素受单基因调控。FEDER等[10]进一步利用该F2群体进行了基因定位，并结合遗传转化鉴定出影响成熟果果皮颜色的基因。另外，利用‘Piel de Sapo’和‘PI 161375’的杂交组合，研究者们也鉴定了果皮颜色相关的QTLs[11-12]。【本研究切入点】随着分子生物学技术的发展，许多甜瓜果皮颜色相关的QTLs和基因被分离和鉴定。但是，多数基因的精细定位仍存在定位区间较大、候选基因众多等情况。【拟解决的关键问题】本研究构建薄皮甜瓜两个幼果果皮颜色差异的纯系甜瓜‘MR-1’和厚皮甜瓜‘LGR’杂交组合，对其进行表型观察、群体构建、遗传模式分析、基因图位克隆等，精细定位幼果果皮颜色基因，为深入揭示甜瓜幼果果皮颜色的形成机制奠定基础。

1 材料与方法

1.1 试验材料

本研究所使用的材料‘MR-1’由国外引进，属于薄皮甜瓜亚种泡瓜变种（var.），果实为近圆形；‘LGR’由浙江省农业科学院蔬菜所甜瓜研究室提供，属于厚皮甜瓜亚种冬甜瓜变种（var.），果实为长椭圆型（图1），成熟果有网纹[13-14]。利用两个材料杂交构建了正反交F1群体、F1与‘LGR’回交BC1F1群体，主要用于遗传分析；BC1F1和F2群体主要用于BSA-seq和遗传定位。杂交和回交构建于2017年春季，定位群体种植于2017年秋季、2018年春秋、2019年春秋以及2020年春季。所有材料均种植在浙江海宁甜瓜实验基地，采用常规大田生产管理。

1.2 果皮颜色性状调查

两亲本甜瓜‘LGR’和‘MR-1’的幼果果皮颜色调查以花后10 d左右的果实为对象，利用目测对比颜色深浅，分别定为“浅绿色”和“深绿色”。遗传分析和定位群体中单株的表型调查也以花后10 d左右的果实为对象，对照两亲本的幼果果皮颜色，经3次不同时间（花后8、10和12 d）的目测定义表型。成熟果颜色调查和测色仪定量测定取自花后35 d左右的果实。利用美能达CR-400/410色彩色差计（Konica Minolta Camera Co., Ltd, Osaka, Japan）对果皮颜色定量，其中L*值代表亮度从黑（0）至白（100），a*值代表从红（+a）至绿（-a），b*值代表从黄（+b）至蓝（-b）。利用SPSS Statistics V21.0进行相关性分析。

1.3 BSA-seq分析

从BC1F1群体中挑选深绿皮和浅绿皮植株各20株，利用CTAB法提取DNA。将20个深绿皮和20个浅绿皮单株DNA分别等量混合，形成两个混池。2个混池和2个亲本DNA按照样品检测、文库构建、质量检测以及上机测序等流程送由公司操作。参考基因组为甜瓜基因组Harukei3 genome reference ver1.41（http://melonet-db.dna.affrc.go.jp/ap/jbh）。

1.4 分子标记开发

基于亲本重测序信息，在初定位区段的双侧和内部挑选有多态性的Indel标记，利用凝胶电泳鉴定分子标记的结合能力和特异性。后利用合适的引物对F2和BC1F1群体中的隐性单株个体进行基因型鉴定，验证初步连锁结果以及缩小定位区间。

1.5 DNA提取、PCR扩增与聚丙烯酰胺凝胶电泳

取新鲜的甜瓜叶片于冷冻干燥机中放置1 d，采用CTAB法提取DNA：将叶片磨成粉末状，加入600 μL预热的CTAB提取液，摇匀，65℃水浴30 min；加入600 μL氯仿﹕异戊醇（24﹕1）抽提液，用力摇晃2 min后，12 000 r/min离心5 min；取上清液，加入360 μL（0.6倍体积上清液）-20℃预冷的异丙醇，-20℃中放置30 min；12 000 r/min离心5 min，弃上清液，加入700 μL 70%无水乙醇洗涤2次，晾干后溶于去离子水中，-20℃保存。

PCR扩增体系为10 μL，其中包括：2×buffer mix 5 μL，上下游引物各0.3 μL，DNA 2 μL，加去离子水2.4 μL。PCR反应程序：98℃预变性2 min，然后扩增33个循环（包括：94℃ 30 s；55℃ 30 s；72℃ 40 min），接着72℃ 5 min，4℃ 10 min。扩增好的样品于4℃保存。

2 结果

2.1 亲本果皮颜色表型

通过对花后10 d的果皮颜色观察，发现薄皮甜瓜‘MR-1’幼果果皮表现为深绿色，花后35 d的成熟果果皮表现为亮黄色；而厚皮甜瓜‘LGR’幼果果皮表现为浅绿色，成熟果果皮表现为灰色有网纹（图1）。

2.2 幼果果皮颜色遗传模式分析

为探究幼果果皮颜色的遗传规律，以‘MR-1’和‘LGR’为亲本进行正反交组配，所得F1代的幼果果皮颜色均与‘MR-1’表型相同，表明幼果果皮颜色受核基因控制，且深绿对浅绿为显性（图1）。F2群体植株的幼果果皮颜色发生分离，其中132个果实为深绿色果皮，44个果实为浅绿色果皮，符合3﹕1的孟德尔遗传分离比，BC1F1代群体表现为1﹕1的分离比（表1）。这些结果表明幼果果皮颜色是受单基因控制的质量性状。因该基因主要控制果皮绿色的深浅，因此将其命名为（）。

图1 亲本以及正反交F1果皮颜色表型

表1 幼果果皮颜色表型的遗传学分析

2.3 GR的初步定位

为了分离和克隆，首先通过从BC1F1群体中挑选浅绿色和深绿色果皮植株各20株进行DNA混池测序（BSA-seq）。参照两亲本的重测序信息，全基因组连锁分析显示4号染色体长臂上存在一个目标性状区域（图2）。该基因座的物理位置为203 744—1 974 725 bp，置信区间跨越1.8 Mb。

图2 GR的初定位（红圈代表目的区段）

2.4 GR的精细定位

基于亲本重测序得到的覆盖了13条染色体，共1 320 497个多态性分子标记，在4号染色体上开发位于初定位区间两侧的标记4-6和4-58，对BC1F1群体中的隐性单株（浅绿皮）进行基因型鉴定。结果显示浅绿皮单株与标记4-6有十分明显的共分离趋势（图3-A），证实了BSA初定位的结果。为了精细定位该基因，开发了初定位结果内特异性好的引物（图3-B，表2），并对扩大的BC1F1和F2群体共984株（BC1F1：877株，F2：107株）进行基因型鉴定，结果共筛选出24个重组交换单株（图4-B）。通过对重组交换单株进行后代表型和基因型鉴定，最终将精细定位于标记4-102和4-81之间，物理距离为17.7 kb（Chr04: 602468-620197）（图4-A、B）。

由于重组交换单株的缺乏，对17.7 kb区间内包含的基因进行序列比对分析。根据最新的甜瓜参考基因组的基因注释发现该区间包含了4个基因：2个新注释基因（和）、1个完整注释基因（）以及的一部分（图4-C）。其中，编码未知蛋白；预测编码一个类生长素蛋白；编码一个双组分反应调节蛋白，该蛋白属于Golden2-like类转录因子，故被命名为CmAPRR2；编码一个烟草花叶病毒复制蛋白。通过序列比对发现和等3个基因的编码区序列在两亲本中存在差异（表3）。其中，仅存在由单碱基替换造成的错义变异，除了3处同义突变和1处错义突变，还存在1处单碱基替换造成的无义变异（表3）。通过对亲本‘MR-1’和‘LGR’测序，证明‘LGR’中编码区第856位碱基（G变为T）产生了终止密码子，致使MELO3C003375蛋白翻译提前终止，导致了Myb-DNA结合结构域大部分缺失（图5）。在黄瓜、番茄以及西瓜中，的同源基因、以及被报道参与叶绿素和类胡萝卜素的合成以及质体发育等过程[8,15-16]。由此推测，（）为控制甜瓜幼果果皮颜色的基因，而单碱基替换造成的无义突变为造成两亲本幼果果皮颜色深浅的关键差异位点。

A：隐性单株在分子标记4-6的基因型。泳道1为‘MR-1’，2为‘LGR’，8、19、58为重组交换单株；B：引物多态性。上方标记为引物名，每个引物左边泳道为‘MR-1’，右边为‘LGR’。M为DNA marker

表2 定位所用引物序列

A：GR初定位于4号染色体4-6和4-58标记之间；B：利用984个单株精细定位GR于4-102和4-81标记约17.7 kb范围内。标记下面的数字代表交换单株数；C：定位区间内的候选基因

2.5 幼果果皮颜色与成熟果皮颜色存在显著相关性

为探究甜瓜幼果果皮颜色与成熟果皮颜色之间的相关性，对BC1F1群体共153个单株的幼果进行目测，并对照亲本将其划分深绿（1）和浅绿（2）两类。因成熟果果皮颜色变化差异较大，使用测色仪定量分析成熟果颜色变异（表4）。通过统计分析和相关性检测，发现幼果果皮颜色与成熟果皮颜色的亮度指标和绿色强度指标的相关系数分别为0.849和0.857，呈极显著相关。这说明幼果果皮颜色与成熟果果皮颜色存在显著相关性。

表3 定位区段内基因编码区变异

参考和比对分别代表‘MR-1’和‘LGR’中碱基和翻译氨基酸REF and ALT represent the base and amino acid in ‘MR-1’ and ‘LGR’, respectively

A：MELO3C003375的基因结构。数字代表到转录起始位点的距离。前面的碱基代表‘MR-1’中的碱基，后面代表‘LGR’中的碱基；B：引起LGR蛋白截短的变异位点；C：‘MR-1’和‘LGR’中MELO3C003375的蛋白结构域

3 讨论

随着甜瓜基因组研究的深入，功能基因组学取得了快速发展[17]。作为重要的农艺性状，果皮颜色的遗传模式一直是研究者关注的热点。4号染色体长臂被广泛报道含有一个调控甜瓜果皮颜色的基因[8,18-20]。OREN等[8]利用3个亲本‘Tam Dew’‘Dulce’和‘Noy Amid’将幼果果皮颜色基因定组位在4号染色体290 kb区段内；ZHAO等[18]利用全基因组关联分析（GWAS）将果皮颜色基因定位于4号染色体500 kb的范围内。因该定位范围内含有一个编码双组分反应调节蛋白CmAPRR2的基因，且该基因的同源基因曾在黄瓜、番茄和胡椒中被报道可调控果皮的叶绿素代谢和色素积累，因此（）被前人认定为候选基因[8,18]。有别于前人使用厚皮甜瓜进行亚种内组配，本研究利用厚皮与薄皮甜瓜开展亚种间杂交；同时，本研究的亲本在果形、网纹以及成熟果果皮颜色等方面与前人选材存在较大差异，这不仅给该性状的遗传定位选材增添了新类型，也暗示了该基因功能在亚种间的保守性[8,13-14]。另外，本研究采用BSA-seq确定初定位区间、重测序开发新标记和图位克隆缩小区段这一思路清晰且完整的研究方案，通过多代表型和基因型的验证，将精细定位在17.7 kb的区间内。但重组交换单株的缺乏使该区间较难进行进一步的缩小。通过对区间内注释基因的测序分析发现（）在‘LGR’中存在1处无义变异，致使该转录因子调控下游基因的重要结构域（Myb-DNA结合结构域）缺失，因此，（）被认为是最佳候选基因。相比以往报道，本研究极大地缩短了幼果果皮颜色性状基因的定位区间，提高了（）作为候选基因的可信度。此外，本研究还发现‘MR-1’亲本中第二个外显子处的单碱基变异与已报道的材料都不相同，这也为甜瓜中的变异类型增添了新认知。

表4 BC1F1群体植株果皮颜色表型

续表4 Continued table 4

-a*/b*代表绿色的强度 -a*/b* represents the intensity of green color

甜瓜果皮颜色丰富多样，FALLIK等[21]对不同成熟度的‘Galia’甜瓜进行了果皮颜色变化调查并将其分成6个阶段：（1）深绿色；（2）绿色；（3）淡黄色夹杂绿色；（4）淡黄色；（5）黄色；（6）深黄色至橙色。REID等[22]通过评估‘Crenshaw’‘Persian’和‘PMR45’3个甜瓜品种发现果实的成熟会伴随着叶绿素含量的下降和类胡萝卜素的上升。而FLÜGEL和GROSS[23]在‘Galia’甜瓜中观察到类胡萝卜素含量在发育期间并没有变化，外果皮变黄主要是因为在成熟过程中叶绿素的降解增加。因此，果皮颜色变化被认为主要由叶绿素和类胡萝卜素（主要是-胡萝卜素）引起[24]。在本研究中，（）作为调控甜瓜幼果果皮颜色的基因，其是否参与了单一色素或者组合色素的合成或者分解过程？在黄瓜中，被发现影响质体发育和叶绿素的合成[15]。PAN等[16]发现番茄和辣椒果实中过表达APRR类转录因子会增加幼果的叶绿素和成熟果的类胡萝卜素。OREN等[8]发现的表达与-类胡萝卜素的含量存在紧密的相关性。另外，与APRR2相关的GLK2转录因子也曾被报道影响质体的发育以及调节植物的叶绿素和类胡萝卜素的含量[25-27]。这些报道都为下一步比较分析亲本色素含量提供了依据和方向。

作为双组分反应调节蛋白（ARR）的一类，APRR2属于B型响应调节蛋白，该类型的蛋白被广泛报道参与细胞分裂素的信号转导[28]。在的三突变体中，ARGYROS等[29]发现叶绿素的含量显著减少并且大部分的细胞分裂素调节基因显著下调。在拟南芥中，CHOI等[30]发现ARR2可通过Myb-DNA结合结构域与TGA3互作，进而通过结合植物防御相关基因启动子的ACGTCATAGA基序来调控的表达。这些展现了ARRs类转录因子对下游基因的调控作用以及Myb-DNA结合结构域对ARRs类转录因子功能的重要影响。在本研究中，‘LGR’与‘MR-1’的（）共有5处SNP差异，其中仅有无义突变在蛋白功能结构域中发生，造成‘LGR’中MELO3C003375（CmAPRR2）的Myb-DNA结合结构域大部分缺失。类似地，OREN等[8]在甜瓜‘Tam Dew’中也发现了相同的无义突变，并且通过等位性测验证实了该变异是造成‘Tam Dew’幼果浅色的原因。基于Myb-DNA结合结构域的功能特点，猜测该无义突变会导致MELO3C003375（CmAPRR2）蛋白无法识别和结合下游基因，进而无法调控色素代谢，以致产生浅绿色果皮。对于APRR2或者ARRs类转录因子调控色素相关基因的方式（直接或者有其他基因介导），暂时未有相关深入的研究。鉴于拟南芥中转座酶类的转录因子FHY3和FAR1，以及ABRE类转录因子ABF2、ABF3和ABF4常通过直接调节叶绿素代谢途径中的结构基因来影响叶绿素水平[31-32]，推测（）通过直接调控色素代谢途径中的结构基因来影响幼果果皮颜色，而缺失了Myb-DNA结合结构域的MELO3C003375（CmAPRR2）则极有可能因为无法直接识别与结合结构基因的启动子，进而无法调节结构基因表达，造成‘LGR’无法积累色素，产生浅绿色幼果。

作为对消费者有吸引力的性状，果皮颜色一直是甜瓜育种的重要目标。相较于成熟果，幼嫩果实颜色相对单一，主要由叶绿素影响。但成熟果与幼嫩果实之间是否存在关联性，研究却很少涉及。本研究中，两亲本幼嫩果果皮颜色差异明显，定位群体（BC1F1）单株表型也分离显著，随着发育进程，亲本以及单株果皮颜色发生变化，但有趣的是，浅绿皮植株和深绿皮成熟果果皮颜色的绿色强度却形成了显著的差异分化，这表明即使果实在发育过程中会发生转变，但成熟果果皮颜色仍然与幼嫩果果皮颜色存在紧密的相关性，该结论为未来进行甜瓜外观品质育种提供了思路。另外，本研究发现在BC1F1群体幼果为深绿皮的植株中，部分单株的成熟果会表现出黄色，而部分单株仍然表现为深绿色。这表明虽然成熟果颜色与幼果颜色有显著相关性，但两者之间的内在联系以及它们与果实发育进程之间的关系还需进一步的探索和解析。这为后续进行成熟果的相关研究提供了思路，也暗示选择合适时期的材料对果实皮色的研究十分重要。

4 结论

甜瓜幼果果皮颜色（深绿和浅绿）是受单个核基因控制的质量性状，且深绿对浅绿为显性。通过遗传定位将该定位于4号染色体分子标记4-102和4-81之间，物理距离约为17.7 kb。此区间内（）的第8个外显子在‘LGR’中发生单碱基替换（G856T），造成蛋白翻译提前终止，推测（）即为调控甜瓜幼果果皮颜色的基因。

[1] TADMOR Y, BURGER J, YAAKOV I, FEDER A, LIBHABER S E, PORTNOY V, MEIR A, TZURI G, SA’AR U, ROGACHEV I, AHARONI A, ABELIOVICH H, SCHAFFER A A, LEWINSOHN E, KATZIR N. Genetics of flavonoid, carotenoid, and chlorophyll pigments in melon fruit rinds. Journal of Agricultural and Food Chemistry, 2010, 58(19): 10722-10728.

[2] GUSMINI G, WEHNER T C. Genes determining rind pattern inheritance in watermelon: A review. HortScience, 2005, 40(6): 1928-1930.

[3] WHITAKER T W, DAVIS G N. Cucurbits: Botany, cultivation and utilization. Journal of Association of Official Agricultural Chemists, 1962, 45(4): 1052.

[4] PARTHASARATHY V A, SAMBANDAM C N. Inheritance in Indian melons. Indian Journal of Genetics and Plant Breeding, 1981, 41(1): 114-117.

[5] KUBICKI B. Inheritance of some characters in muskmelons. Genetiea Polonlca, 1962, 3: 265-274.

[6] PEREIRA L, RUGGIERI V, PÉREZ S, ALEXIOU K G, FERNÁNDEZ M, JAHRMANN T, PUJOL M, GARCIA-MAS J. QTL mapping of melon fruit quality traits using a high-density GBS-based genetic map. BMC Plant Biology, 2018, 18(1): 324.

[7] BURGER Y, BHASTEKER D, SA’AR U, KATZIR N, PARIS H S. A recessive gene for light immature exterior color of melon. Cucurbit Genetics Cooperative Report, 2006, 28/29: 17-18.

[8] OREN E, TZURI G, VEXLER L, DAFNA A, MEIR A, FAIGENBOIM A, KENIGSWALD M, PORTNOY V, SCHAFFER A A, LEVI A, BUCKLER E S, KATZIR N, BURGER J, TADMOR Y, GUR A. The multi-allelic APRR2 gene is associated with fruit pigment accumulation in melon and watermelon. Journal of Experimental Botany, 2019, 70(15): 3781-3794.

[9] HUGHES M. The inheritance of two characters ofand their interrelationship. Journal of the American Aociety for Horticultural Science, 1948, 52: 399-402.

[10] FEDER A, BURGER J, GAO S, LEWINSOHN E, KATZIR N, SCHAFFER A A, MEIR A, DAVIDOVICH-RIKANATI R, PORTNOY V, GAL-ON A, FEI Z J, KASHI Y, TADMOR Y. A kelch domain- containing F-box coding gene negatively regulates flavonoid accumulation in muskmelon. Plant Physiology, 2015, 169(3): 1714-1726.

[11] EDUARDO I, ARÚS P, MONFORTE A J, OBANDO J, FERNANDEZ-TRUJILLO J P, MARTINEZ J A, ALARCÓN A L, ALVAREZ J M, VAN DER KNAAP E. Estimating the genetic architecture of fruit quality traits in melon (L.) using a genomic library of near-isogenic lines. Journal of the American Society for Horticultural Science, 2007, 132(1): 80-89.

[12] OBANDO J, FERNÁNDEZ-TRUJILLO J P, MARTÍNEZ J A, ALARCÓN A L, EDUARDO I, ARÚS P, MONFORTE A J. Identification of melon fruit quality quantitative trait loci using near-isogenic lines. Journal of the American Society for Horticultural Science, 2008, 133(1): 139-151.

[13] PITRAT M, HANELT P, HAMMER K. Some comments on infraspecific classification of cultivars of melon. Acta Horticulturae, 2000, 510: 29-36.

[14] 林德佩. 中国栽培甜瓜植物的起源、分类及进化. 中国瓜菜, 2010(4): 34-36.

LIN D P. Origin classification and evolution for cultivated plants of Chinese melon. China Cucurbits and Vegetables, 2010(4): 34-36. (in Chinese)

[15] LIU H Q, JIAO J Q, LIANG X J, LIU J, MENG H W, CHEN S X, LI Y H, CHENG Z H. Map-based cloning, identification and characterization of the w gene controlling white immature fruit color in cucumber (L.). Theoretical and Applied Genetics, 2016, 129(7): 1247-1256.

[16] PAN Y, BRADLEY G, PYKE K, BALL G, LU C G, FRAY R, MARSHALL A, JAYASUTA S, BAXTER C, VAN WIJK R, BOYDEN L, CADE R, CHAPMAN N H, FRASER P D, HODGMAN C, SEYMOUR G B. Network inference analysis identifies an APRR2-like gene linked to pigment accumulation in tomato and pepper fruits. Plant Physiology, 2013, 161(3): 1476-1485.

[17] 贾继增, 高丽锋, 赵光耀, 周文斌, 张卫健. 作物基因组学与作物科学革命. 中国农业科学, 2015, 48(17): 3316-3332.

JIA J Z, GAO L F, ZHAO G Y, ZHOU W B, ZHANG W J. Crop genomics and crop science revolutions. Scientia Agricultura Sinica, 2015, 48(17): 3316-3332. (in Chinese)

[18] ZHAO G W, LIAN Q, ZHANG Z H, FU Q S, HE Y H, MA S W, RUGGIERI V, MONFORTE A J, WANG P Y, JULCA I, WANG H S, LIU J P, XU Y, WANG R Z, JI J B, XU Z H, KONG W H, ZHONG Y, SHANG J L, PEREIRA L, ARGYRIS J, ZHANG J A, MAYOBRE C, PUJOL M, OREN E, OU D D, WANG J M, SUN D X, ZHAO S J, ZHU Y C, LI N, KATZIR N, GUR A, DOGIMONT C, SCHAEFER H, FAN W, BENDAHMANE A, FEI Z J, PITRAT M, GABALDÓN T, LIN T, GARCIA-MAS J, XU Y Y, HUANG S W. A comprehensive genome variation map of melon identifies multiple domestication events and loci influencing agronomic traits. Nature Genetics, 2019, 51(11): 1607-1615.

[19] 杨光华, 范荣, 杨小锋, 侯军亮, 袁士臣, 曹明, 王学林, 李劲松. 甜瓜果实颜色3个质量性状基因的定位. 园艺学报, 2014, 41(5): 898-906.

YANG G H, FAN R, YANG X F, HOU J L, YUAN S C, CAO M, WANG X L, LI J S. Construction of a highly dense genetic map using SNP and mapping of three qualitative traits in. Acta Horticulturae Sinica, 2014, 41(5): 898-906. (in Chinese)

[20] 欧点点, 赵光伟, 贺玉花, 王平勇, 徐志红, 孔维虎, 张健, 徐永阳. 甜瓜果皮颜色遗传分析及基因定位. 中国农学通报, 2019, 35(13): 64-69.

OU D D, ZHAO G W, HE Y H, WANG P Y, XU Z H, KONG W H, ZHANG J A, XU Y Y. Genetic analysis and gene mapping for melon rind color. Chinese Agricultural Science Bulletin, 2019, 35(13): 64-69. (in Chinese)

[21] FALLIK E, ALKALI-TUVIA S, HOREV B, COPEL A, RODOV V, AHARONI Y, ULRICH D, SCHULZ H. Characterization of ‘Galia’ melon aroma by GC and mass spectrometric sensor measurements after prolonged storage. Postharvest Biology and Technology, 2001, 22(1): 85-91.

[22] REID M S, LEE T H, PRATT H K, CHICHESTER C O. Chlorophyll and carotenoid changes in developing muskmelons. Journal of the American Society for Horticultural Science, 1970, 95(6): 814-815.

[23] FLÜGEL M, GROSS J. Pigment and plastid changes in mesocarp and exocarp of ripening muskmeloncv. Galia. Angewandte Botanik, 1982, 56(5/6): 393-406.

[24] LIGHTBOURN G J, GRIESBACH R J, NOVOTNY J A, CLEVIDENCE B A, RAO D D, STOMMEL J R. Effects of anthocyanin and carotenoid combinations on foliage and immature fruit color ofL. Journal of Heredity, 2008, 99(2): 105-111.

[25] CHEN M, JI M L, WEN B B, LIU L, LI S X, CHEN X D, GAO D S, LI L. GOLDEN 2-LIKE transcription factors of plants. Frontiers in Plant Science, 2016, 7: 1509.

[26] WATERS M T, MOYLAN E C, LANGDALE J A. GLK transcription factors regulate chloroplast development in a cell-autonomous manner. The Plant Journal, 2008, 56(3): 432-444.

[27] WATERS M T, WANG P, KORKARIC M, CAPPER R G, SAUNDERS N J, LANGDALE J A. GLK transcription factors coordinate expression of the photosynthetic apparatus in. The Plant Cell, 2009, 21(4): 1109-1128.

[28] MASON M G, MATHEWS D E, ARGYROS D A, MAXWELL B B, KIEBER J J, ALONSO J M, ECKER J R, SCHALLER G E. Multiple type-B response regulators mediate cytokinin signal transduction in. The Plant Cell, 2005, 17(11): 3007-3018.

[29] ARGYROS R D, MATHEWS D E, CHIANG Y H, PALMER C M, THIBAULT D M, ETHERIDGE N, ARGYROS D A, MASON M G, KIEBER J J, SCHALLER G E. Type B response regulators ofplay key roles in cytokinin signaling and plant development. The Plant Cell, 2008, 20(8): 2102-2116.

[30] CHOI J, HUH S U, KOJIMA M, SAKAKIBARA H, PAEK K H, HWANG I. The cytokinin-activated transcription factor ARR2 promotes plant immunity via TGA3/NPR1-dependent salicylic acid signaling in. Developmental Cell, 2010, 19: 284-295.

[31] TANG W J, WANG W Q, CHEN D Q, JI Q A, JING Y J, WANG H Y, LIN R C. Transposase-derived proteins FHY3/FAR1 interact with PHYTOCHROME-INTERACTING FACTOR1 to regulate chlorophyll biosynthesis by modulating HEMB1 during deetiolation in. The Plant Cell, 2012, 24(5): 1984-2000.

[32] GAO S, GAO J, ZHU X Y, SONG Y, LI Z P, REN G D, ZHOU X, KUAI B K. ABF2, ABF3, and ABF4promote ABA-mediated chlorophyll degradation and leaf senescence by transcriptional activation of chlorophyll catabolic genes and senescence-associated genes in. Molecular Plant, 2016, 9(9): 1272-1285.

Fine Mapping of an Immature Rind Color Genein Melon

XU XinYang, SHEN Jia, ZHANG YueJian, LI GuoJing, NIU XiaoWei, SHOU WeiSong

Institute of Vegetables, Zhejiang Academy of Agricultural Sciences, Hangzhou 310021

【】The aim of this study was to explore the inheritance pattern of immature rind color of melon(L.), to fine map target gene, and to deepen the understanding of rind color change during fruit growth, so as to provide a guidance for the improvement of melon color with molecular design breeding.【】The dark-green rind line ‘MR-1’ (.ssp.) and the light-green rind line ‘LGR’ (C.ssp.) were used as parents to construct the F1hybrid population, and the BC1F1backcross population was constructed from crossing by F1and ‘LGR’ for genetic analysis of immature rind color. By selecting 20 dark-green and light-green plants each in BC1F1population, the DNA was mixed, respectively, and the BSA-seq was operated for initial mapping. Based on the resequencing of two parents, the molecular markers with better specificity were developed in the initial region and recombinant individuals were identified and selected to verify and narrow the interval for fine mapping thegene. By sequencing of coding regions between two parents according to the gene annotation, the candidate gene and key variant site were determined. Moreover, the phenotype of immature and mature fruit in the BC1F1backcross population was investigated and assessed by correlation analysis to explore the underlying mechanism of rind color transition in the fruit development.【】According to the investigation, the phenotype of all F1individual plants exhibited dark-green rind color. In addition, the immature rind color of BC1F1 backcross population was found to be separated, and the ratio of the number of dark-green to light-green was approximately 1﹕1. As well as, the ratio of the number of dark-green to light-green in the F2population was 3﹕1. These ratios corresponded to Mendel’s law of inheritance, indicating that the immature rind color of melon was a quality trait, controlled by a single nuclear gene, and the dark-green was dominant to light-green. Through BSA-seq, the gene was initially mapped to a 1.8 Mb interval on the long arm of chromosome 4. With developed molecular markers, 24 recombinant individuals were selected in expanded mapping population. The gene was further narrowed down to a region between markers 4-102 and 4-81 with a physical distance of 17.7 kb by genotype and phenotype verification of offspring, where were four predicted genes with latest annotation. Sequencing analysis revealed that a geneencoding GLKs transcription factor CmAPRR2 had several variations in ‘MR-1’ and ‘LGR’. Among them, there were three synonymous, a missense mutation, and a nonsense mutation. The nonsense mutation (G to T) appeared in the 856th base of the coding region led to premature translation termination and most of the Myb-DNA binding domain to be lost in ‘LGR’. Thus, the() was speculated to be thegene and the nonsense mutation was the key variation that affected the immature rind color of melon. It was found that there was a significant correlation between the rind color of immature fruit and mature fruit by investigating the phenotype of fruits in the BC1F1backcross population.【】Immature rind color (dark green/light green) of melon was a quality trait and controlled by a single nuclear gene. By mapping,() was presumed to be the candidate gene that affected immature rind color.

melon; immature rind color; gene mapping;;()