甜瓜短蔓基因Cmdm1 的精细定位及候选基因分析

奥鹏网院作业

甜瓜短蔓基因Cmdm1 的精细定位及候选基因分析
马建1，李丛丛2，王建设1

（1 北京市农林科学院蔬菜研究中心/农业农村部华北地区园艺作物生物学与种质创制重点实验室，北京 100097；2 北京市农林科学院北京农业生物技术研究中心/农业基因资源与生物技术北京市重点实验室，北京 100097）

摘要：【目的】对甜瓜短蔓突变体Z8 进行短蔓基因的精细定位并确定候选基因，为甜瓜株型的分子改良奠定基础。【方法】考察短蔓突变体Z8 和野生型B15 的主蔓节数、主蔓长度、主蔓节间长度以及侧枝长度等农艺性状。配制Z8/B15 杂交组合并进行遗传分析，利用F2群体中的短蔓单株进行基因精细定位。通过对定位区间内注释基因编码区进行测序以确定候选基因。【结果】与野生型B15 相比，突变体Z8 节间显著变短导致植株矮化，顶端花序紧凑簇生，遗传分析表明其短蔓性状由一对隐性核基因Cmdm1 控制。采用基因图位克隆策略，利用780 个F2短蔓单株最终将该基因精细定位于第7 染色体短臂标记c7-112 和s2 之间约56 kb 的区间内，并与标记dm-1 共分离，区间内共包含4 个注释基因。经测序鉴定，发现Z8 中与拟南芥ERECTA 同源的MELO3C016916 ATG 下游第1 995 位碱基由T 突变为G 而产生终止密码子，导致蛋白翻译提前终止，致使后面激酶结构域完全缺失，推测MELO3C016916即为控制蔓长的Cmdm1。【结论】Z8 短蔓性状受隐性核基因Cmdm1 控制，利用分子标记最终将该基因定位于7 号染色体短臂标记c7-112 和s2 之间约56 kb 区间内，推测MELO3C016916 为最有可能的候选基因。

关键词：甜瓜；短蔓；基因定位；Cmdm1；ERECTA

0 引言
【研究意义】甜瓜（Cucumis melo L.）为葫芦科（Cucurbitaceae）甜瓜属（Cucumis）甜瓜种一年生草本植物，是一种世界范围内广泛种植的重要园艺经济作物。株型对于甜瓜的生产栽培具有重要意义，可影响甜瓜的种植密度、产量、品质、坐果性等方面，而节间长度是甜瓜株型的重要构成因素，短蔓（节间）品种往往具有株型紧凑的特点，生产上可通过增加种植密度来提高单位面积产量而增加经济效益。因此，短蔓性状作为株型的关键性状，逐步受到育种家的关注。深入挖掘并利用甜瓜短蔓资源，鉴定短蔓控制基因，从而利用分子手段培育短蔓甜瓜新品种，进而解决甜瓜栽培上的密植问题，对促进甜瓜产业化发展具有重要意义。【前人研究进展】研究表明，植物矮化基因在抗倒伏、产量、早熟、分蘖能力上都有潜在作用而具有重要的经济价值[1]。20 世纪五六十年代，由于半矮杆基因semidwarf 1（sd-1）的应用，在世界范围内掀起了一场“绿色革命”，使作物产量大幅度提高[2]。葫芦科作物是蔬菜作物的主要大类之一，常见的有黄瓜、西瓜、南瓜、甜瓜等10 余种，而短蔓现象似乎在葫芦科作物中普遍存在，比如已报道的黄瓜上的cp[3]、cp-2[4]、scp[5]、si[6]和csdm1[7]等矮化突变体；西瓜上的dw-1[8]和DMSW（dwarf male sterile watermelon）[9]等矮化材料；南瓜上利用无蔓（矮生）突变体育成的无蔓1 号和无蔓4 号无蔓南瓜的新品种以及无蔓南瓜 Bu[10]、南瓜矮化突变体 dm1[11]和SQ026[12]等突变体。甜瓜的短蔓性状研究起步较早，1968 年，美国学者DENNA 报道了一份甜瓜短蔓甜瓜材料UC Topmark bush，后来表明其携带基因si-1，si-1表达受1 个至多个修饰基因影响，其短蔓植株株型紧凑、节间细、脆、短、蔓弯曲、果小，属于矮生型[13]；PARIS 等[13]报道了Persia202 携带基因si-2，其短蔓植株株型紧凑直立，座瓜靠近茎基部，属于丛生型或鸟巢型；KNAVEL[14]报道了品种Maindwarf 携带基因si-3，报道的3 个基因彼此间非等位且为独立遗传的隐性单基因。HWANG 等[15]对一份短节间甜瓜PNU-D1的遗传模式进行了研究，表明短节间表型受一对隐性基因mdw1 控制。王建设等[16]以长蔓品种状元和短蔓资源1A533、1A440 为材料，分析了甜瓜短蔓性状的遗传规律，研究结果表明两份短蔓甜瓜各携带l 对非等位的隐性短蔓基因。随着瓜类作物基因组的测序完成，近几年通过图位克隆等方法在瓜类作物上也陆续鉴定出一些植株矮化基因，LI 等[3]将源于黄瓜PI 308915 的cp 精细定位到220 kb 区间内，鉴定到其中的一个细胞分裂素氧化酶基因（cytokinin oxidase，CKX），该基因由于第一个外显子上发生3 bp 缺失而造成矮化表型。LIN 等[6]对一份EMS 诱变来源的短节间矮化黄瓜突变体si 进行了研究，发现编码F-box 蛋白的CsaVFB1（Csa4P641640）与表型最相关。XU等[7]对黄瓜矮化突变体Csdw 进行研究，鉴定出一个编码CLAVATA1 类型的受体激酶基因Csa3G872760（CsCLAVATA1）发生突变，可能是导致Csdw 产生矮化表型的原因。ZHANG 等[12]对一份矮化突变体SQ026 进行研究，利用genotyping-by-sequencing 方法鉴定到一个编码赤霉素20 氧化酶基因（gibberellin 20-oxidase，qCmB2），该基因的启动子区在突变体中存在1 249 bp 缺失，可能与矮化表型相关。HWANG等[15]对短蔓甜瓜PNU-D1 中的隐性控制基因mdw1 进行了研究，最终将该基因定位于1.8 cM 距离内。【本研究切入点】茎蔓长度是甜瓜株型重要的性状之一，直接影响单位面积的种植密度。尽管前人已对甜瓜的短蔓性状进行了大量的研究，但尚未提升到基因水平。王建设等[16]研究表明短蔓自然突变体材料1A440 的短蔓性状受l 对隐性基因控制，本研究在此基础上通过对1A440 再次进行自交提纯获得了稳定遗传的甜瓜短蔓突变体Z8，在此基础上拟对Z8 中的短蔓基因进行精细定位并分析可能的目的基因。【拟解决的关键问题】本研究以甜瓜短蔓突变体Z8 为材料，通过表型和遗传分析、基因定位和候选基因分析等研究手段，拟确定造成短蔓表型的关键目的基因，为将来该基因的转基因功能验证以及利用分子标记辅助手段 培育短蔓甜瓜新品种，提供资源、基因和方法。

1 材料与方法
试验于2018 年在北京市农林科学院蔬菜研究中心进行。

1.1 实验材料、群体构建及表型调查
父本材料：厚皮甜瓜短蔓材料1A440 自交5 代获得的短蔓突变体Z8，母本材料：源于甜瓜地方品种‘一窝猴’的高代自交系B15。两者杂交、自交获得BC1F1（31 株）和F2（780 株）代分离群体，所有材料于春季育苗后移栽于北京市农林科学院蔬菜研究中心四季青基地大棚内。选取父、母本各10 个单株自然生长，不整枝不掐头，于果实成熟期分别调查各单株主蔓节数，用软尺测量各单株主蔓以及基部最长的5 个侧枝长度，并进行数据统计分析。

1.2 DNA 提取及PCR 扩增、检测
取甜瓜新鲜叶片，利用CTAB 法[17]提取基因组总DNA。PCR 扩增反应体系为10 μL，包括0.8 μL Buffer、0.8 μL dNTPs（2.5 mmol·L-1 each）、1 μL Primers（10μmol·L-1，F+R）、1 μL 基因组DNA（50 ng·μL-1）、0.1 μL EasyTaq® DNA Polymerase for PAGE（购自北京全式金生物技术有限公司）和6.3 μL ddH2O。反应程序：94℃预变性4 min；94℃变性30 s，59℃退火30 s，72℃延伸30 s，运行35 个循环；72℃延伸5 min。扩增产物采用8%聚丙烯酰胺凝胶电泳检测，银染后在观片机上拍照、记录。

1.3 RNA 提取与基因测序、表达分析
两周左右的甜瓜幼苗经液氮研磨后，利用天根生化科技（北京）有限公司的RNAprep pure Plant Kit 试剂盒进行总RNA 提取，第一链cDNA 采用北京全式金生物技术有限公司的TransScript® One-Step gDNA Removal and cDNA Synthesis SuperMix 试剂盒合成，操作步骤参照说明书。扩增候选基因长片段PCR 反应体系如下：2×PCR Buffer for KOD FX 25 μL，10 μL dNTPs（2 mmol·L-1 each），2 μL 引物（10 μmol·L-1，F+R），1 μL KOD FX 酶（1 U·μL-1，购自东洋纺（上海）生物科技有限公司），2 μL DNA 模板（100 ng·μL-1），加ddH2O 至50 μL。扩增反应条件：94℃预变性2 min；98℃变性10 s，59℃退火30s，68℃依片段大小延伸～1 kb/min，35 个循环；68℃延伸5 min。扩增产物经1%琼脂糖胶检测后送至北京博迈德科技发展有限公司测序，获得序列用DNAMAN 软件拼接比对，确定突变位点。

基因表达分析采用实时定量荧光PCR 方法[18]，利用宝日医生物技术（北京）有限公司（TaKaRa）的SYBR Premix Ex TaqTM II 试剂盒进行扩增，选用Applied Bio systems公司的96孔PCR板在ABI PRISM 7500 实时定量PCR 仪上运行。以甜瓜的CmActin（AY859055）作为内参[19]，反应程序：95℃预变性15 s；95℃变性5 s，60℃延伸34 s，40 个循环。每个反应重复3 次，采用2-ΔΔCT 算法分析结果[20]。引物序列见表1。

1.4 基因定位及引物的开发
参考甜瓜基因组数据库（https://www.melonomics. net/）公布的基因组及基因注释数据，利用在线引物设计程序Primer-BLAST（http://www.ncbi.nlm.nih.gov/ tools/primer-blast/index.cgi?LINK_LOC=BlastHome）设计引物，引物由北京博迈德科技发展有限公司合成，引物序列见表1。

2 结果
2.1 突变体Z8 植株的表型分析
与野生型对照B15 相比，突变体Z8 在子叶期即表现为矮化（图1-A），伸蔓期植株生长较慢，节间长度变短，植株矮化（图1-B）。此外，突变体另一个明显的特征是顶端花序紧凑簇生（图1-C），最终影响到整个植株形态建成。突变体Z8 主蔓平均长度为（181±13.9） cm，表现为短蔓，与对照差异极显著（P=0.006）；平均节间长度（3.2±0.2 cm）与对照差异极显著（P＜0.01）；而总节数和平均最长侧枝长度与对照无显著差异（表2）。

2.2 突变体Z8 短蔓性状的遗传分析
为了解突变体短蔓形成的遗传机制，以突变体Z8为父本，野生型B15 为母本配制杂交组合，所得F1代植株蔓长均表现为野生型B15 表型，表明Z8 的短蔓性状受隐性核基因控制（表2）。F2 代群体中蔓长出现分离，其中似野生型植株175 株，短蔓植株55株，表现3﹕1 的孟德尔遗传分离比（χ2=0.4＜3.84），BC1F1 代群体表现为1﹕1 的分离比（表2），这些结果表明Z8 中的短蔓性状由单隐性核基因控制，暂定该基因名称为Cmdm1（dwarf melon 1）。

2.3 Cmdm1 的定位和物理图谱构建
为分离控制蔓长性状的Cmdm1，首先，利用两亲本B15 和Z8 的DNA 对实验室保存的分布在12 对染色体上的331 对InDel 引物进行筛选，其中表现出多态性的引物共120 对。随后，选取15 个短蔓F2单株DNA等量混合建立混池，加上亲本B15 和Z8，共3 个样品对120 对引物进行筛选。结果表明，7 号染色体短臂引物c7-3 和c7-206 与混池DNA 表现出共分离。进一步利用两对引物对55 个短蔓F2 单株进行分子标记数据分析，发现引物与短蔓表型表现出连锁关系，并将控制短蔓性状的Cmdm1 初步定位于两个引物c7-3 和 c7-206 之间（图2-A）。

表1 本研究所用的PCR 引物序列
Table 1 PCR primer sequences used in this study

pagenumber_ebook=159,pagenumber_book=805
pagenumber_ebook=159,pagenumber_book=805
图1 B15（野生型）和Z8（短蔓）表型
Fig. 1 Phenotypes of wild-type B15 and Z8 mutant

A：子叶期表型（标尺1 cm）；B：伸蔓期表型（标尺20 cm）；C：Z8 植株顶端表型（红色箭头所示，标尺5 cm）
A: Cotyledon stage (Scale bars, 5 cm); B: Vine elongation stage (Scale bars, 20 cm); C: The top inflorescence of Z8 (Scale bars, 5 cm)

表2 B15（野生型）和Z8（突变体）4 个农艺性状比较分析
Table 2 Comparison of agronomic traits between the wild-type B15 and Z8 mutant (mean±SD, n=10)

pagenumber_ebook=160,pagenumber_book=806
**：差异极显著 Significant differences (P＜0.01)

表3 Z8 短蔓性状的遗传分析
Table 3 Genetic analysis of the Z8 mutant

pagenumber_ebook=160,pagenumber_book=806
pagenumber_ebook=160,pagenumber_book=806
图2 Cmdm1 的物理图谱
Fig. 2 Physical maps of the Cmdm1 locus

A：Cmdm1 定位于7 号染色体短臂c7-3 和c7-206 标记之间；B：利用210 个短蔓F2 单株的初定位；标记下面的数字代表重组单株数目；C：利用780个短蔓F2 单株的精细定位；D：定位区间的物理图谱，箭头代表注释的4 个基因
A: the Cmdm1 gene was located on the short arm of chromosome 7 between markers c7-3 and c7-206; B: the Cmdm1 gene was preliminarily delimited to the c7-101-c7-132 interval using 210 short internodes F2 individuals; the numerals below the corresponding markers indicate the numbers of recombinants (Rec) identified among F2 plants; C: fine genetic mapping of the Cmdm1 gene based on 780 short internodes F2 individuals; D: physical map of the Cmdm1 region based on the reference sequence and arrows indicated four annotated ORFs in this region

在此基础上，根据甜瓜参考基因组序列（https:// www.melonomics.net/）在两引物之间开发出7 对多态性分子标记，随后对210 个短蔓F2 单株进行分子标记数据分析，初步将Cmdm1 定位于标记c7-101 和c7-132之间（图2-B）。进一步将定位群体扩大到780 个短 蔓F2 单株，最终将Cmdm1 精细定位于标记c7-112 和s2 之间，物理距离约为56 kb（chr7：1127514—1183716），并与标记dm-1 共分离（图2-C）。

pagenumber_ebook=161,pagenumber_book=807
图3 MELO3C016916 结构、蛋白结构域及表达分析
Fig. 3 Gene structure, predicated function domains and gene expression levels of MELO3C016916 gene

A：MELO3C016916 基因结构及其突变位点分析；B：野生型B15 和突变体Z8 中的MELO3C016916 基因编码的蛋白结构域示意图
A: Gene structure of MELO3C016916 and T1995 to G1995 mutation site was identified in the 25th exon generating a premature termination codon; B: The predicted function domains of MELO3C016916 protein in B15 and Z8, respectively

2.4 定位区间内候选基因分析
根据甜瓜参考基因组对该区间内的基因注释，56 kb 区间内共包含 4 个预测的编码基因，分别为MELO3C016916、MELO3C016917、MELO3C016918和MELO3C016919（图2-D）。其中MELO3C016916编码一个富含亮氨酸重复序列的类受体激酶蛋白（LRR receptor-like serine/threonine-protein kinase），MELO3C016917 编码一个未知功能蛋白（Uncharacterized protein），MELO3C016918 编码一个脂肪酸结合蛋白2（Fatty-acid-binding protein 2），MELO3C016919 编码鸟嘌呤核苷酸结合蛋白α-1 亚基（Guanine nucleotide- binding protein alpha-1 subunit）。随后，分别对野生型B15 和突变体Z8 中4 个基因的编码区DNA 进行测序，结果显示，相比B15 和参考基因组序列，Z8 中MELO3C016917编码区发生两处未造成编码氨基酸变化的同义突变C180T 和C255T；MELO3C016918 编码区发生一处点突变A359G，造成编码氨基酸由天冬酰胺突变为丝氨酸；MELO3C016919 编码区无差异；MELO3C016916 编码区共发生4 处同义突变，分别位于ATG 下游C381T、C961T、C2112T 和T2388C 处。此外，Z8 中该基因第1995 位置处碱基位点由T 突变为G，造成密码子由TAT 突变为终止密码子TAG（图3-A）。基因组注释表明MELO3C016916 编码区包含27 个外显子，突变位点位于第25 外显子。SMART 软件（http://smart.embl-heidelberg.de/）预测该基因编码991 个氨基酸组成的类受体激酶蛋白，突变造成翻译蛋白提前终止，仅编码664 个氨基酸组成的截短蛋 白，造成后面激酶结构域完全缺失（图3-B），从而影响了该蛋白的功能。另外，荧光定量PCR 检测显示，与对照B15 相比，Z8 中该基因表达水平明显降低（图4）。MELO3C016916与已报道的拟南芥的ERECTA同源[21]，编码蛋白的氨基酸序列相似性为88.6%，推测MELO3C016916 即为控制蔓长的基因Cmdm1。

pagenumber_ebook=162,pagenumber_book=808
图4 MELO3C016916 在B15 和Z8 中的相对表达分析
Fig. 4 Relative expression analysis of MELO3C016916 in B15 and Z8

3 讨论
“高产、优质、高效”始终是甜瓜生产追求的目标，合理密植是实现这一目标的关键措施之一。甜瓜的种植密度依赖于品种类型、熟性、种植方式、土壤肥力等条件。目前，西北甜瓜主产区以露地匍匐栽培为主，现有的推广品种瓜蔓长达3 m 以上，栽培上株距、行距较大；东部地区保护地栽培一般采用直立吊蔓栽培，与匍匐栽培相比，利用了空间资源，但生育前期需要不间断绑蔓、吊蔓，生育后期基部叶片相互遮蔽，光合效率有所降低[16]。因此，在现有栽培方式下，通过改良甜瓜品种株型，有效利用土地及光热资源，在改善品质的同时进一步提高产量是甜瓜株型育种的新方向。主蔓长度是甜瓜株型的重要构成因素，利用短蔓基因，培育株型紧凑、适应密植、能充分利用土地与光热资源而且早熟的短蔓品种，可显著降低绑蔓、吊秧次数，提高单位面积产量，从而节约劳动力、降低成本、提高经济效益。之前的报道中，美国科学家就曾选育出‘Mininesota 266’‘U.C. Grenshaw Bush’‘Honey Dew Bush’‘Top Mark Bush’‘Perlita Bush’‘G515’‘KyP7’‘Main Dwarf’等早熟、丰产、优质和短蔓的甜瓜新品种[13-14,22-23]，也先后报道了3 个独立遗传的短蔓基因si-1、si-2 和si-3[13-14]。本研究的短蔓突变体Z8 源于短蔓材料1A440 的多代自交系，短蔓性状受一对隐性基因控制，该突变体幼苗期即表现矮化，成株期株型紧凑、主蔓较短、顶端花序紧凑簇生，对该突变体的研究极大的丰富了甜瓜的短蔓资源。此外，OH[24]研究数据表明，短蔓材料PNU-D1 的短蔓性状与短节间表现出正相关性（r2=0.69，P＜0.005），即PNU-D1 的侧枝长度也变短，但本研究结果表明Z8 中5 个最长侧枝长度与对照无明显的差异，这与OH[24]报道不一致，可能是由于材料不同，也可能是统计方法不同所造成。

基因组学研究是近年来发展最为迅速的一个研究领域[25]。随着甜瓜全基因组测序的完成，甜瓜功能基因组学取得了快速发展，顺利定位和克隆了一批控制甜瓜重要园艺性状的基因，如甜瓜枯萎病抗性基因Fom1[26]、橙色果肉基因CmOr[27]、卷须控制基因CmTCP1[28]、白粉病抗性基因Pm-Edisto47-1[29]等。2014 年，HWANG 等[15]对一份韩国的短蔓突变体PNU-D1 进行研究，通过构建遗传连锁图谱最终将控制蔓长的mdwl 定位于7 号染色体1.8 cM 内，两侧标记mERE 和mUBI 距离分别为0.6 cM 和1.7 cM，但定位区间较大，作者没有确定最终的候选基因。本研究通过对短蔓突变体Z8 中的短蔓基因进行精细定位，最终将控制蔓长的Cmdm1 定位于7 号染色体短臂56 kb 区间内，该位置与mdwl 所在位置距离较近，不排除两者互为等位基因的可能性。通过对定位区间内候选基因编码区序列分析，发现定位区间内包含一个拟南芥基因 ERECTA（AT2G26330）的同源基因MELO3C016916，两者编码蛋白序列相似性为88.6%，ERECTA 是从拟南芥中分离到的一个编码富含亮氨酸重复序列的类受体激酶基因（leucine-rich repeats receptor-like serine/threonine-protein Kinase，LRR- RLK）[18]。拟南芥ERECTA 在调控顶端分生组织的发育中起重要作用，该基因发生突变会导致节间缩短、花梗变短，具有紧凑的花序和粗而直立的花序茎，造成花簇集中在茎的顶端，最终影响植株形态建成[18]，这些表型与Z8 的表型极为相似。突变体Z8 中该基因第25外显子发生T1995G 单碱基突变，造成编码蛋白的激酶结构域完全缺失，表达水平显著降低，而拟南芥的研究表明ERECTA 参与了生长素信号途径[30]，C 端激酶结构域对下游信号转导具有重要作用[31]。因此，本研究推测MELO3C016916为Cmdm1最终的候选基因。将来利用转基因手段对该基因功能进行验证，将为深入阐明Cmdm1 控制节间长度的分子机制奠定基础，而且为培育短蔓甜瓜新品种提供支持。

4 结论
甜瓜短蔓突变体Z8 平均节间缩短，顶端花序簇生。遗传分析表明该突变性状受一对隐性核基因Cmdm1 控制，利用分子标记最终将该基因定位于7号染色体短臂标记c7-112 和s2 之间，物理距离约为 56 kb，并与标记dm-1 共分离。此定位区间内MELO3C016916 与拟南芥的ERECTA 同源，Z8 中该基因发生T1995G 碱基突变，造成编码蛋白提前终止，表达水平显著降低，推测该基因为控制蔓长的Cmdm1。

References

[1] EVANS L T. Crop Evolution, Adaptation and Yield. Cambridge: Cambridge University Press, 1993.

[2] SASAKI A, ASHIKARI M, UEGUCHI-TANAKA M, ITOH H, NISHIMURA A, SWAPAN D, ISHIYAMA K, SAITO T, KOBAYASHI M, KHUSH G S, KITANO H, MATSUOKA M. Green revolution: A mutant gibberellin-synthesis gene in rice. Nature, 2002, 416(6882): 701-702.

[3] LI Y H, YANG L M, PATHAK M, LI D W, HE X M, WENG Y Q. Fine genetic mapping of cp: A recessive gene for compact (dwarf) plant architecture in cucumber, Cucumis sativus L. Theoretical and Applied Genetics, 2011, 123(6): 973-983.

[4] KUBICKI B, SOLTYSIAK U. KORZENIEWSKA A. Induced mutations in cucumber (Cucumis sativus L.) V. Compact type of growth. Genetica Polonica, 1986, 27(3/4): 289-298.

[5] NIEMIROWICZ-SZCZYTT K, RUCINSKA M, KORZENIEWSIA A. An induced mutation in cucumber (Cucumis sativus L.): Super compact. Cucurbit Genetics Cooperative Report, 1996, 19: 1-3.

[6] LIN T, WANG S H, ZHONG Y, GAO D L, CUI Q Z, CHEN H M, ZHANG Z H, SHEN H L, WENG Y Q, HUANG S W. A truncated F-Box protein confers the dwarfism in cucumber. Journal of Genetics and Genomics, 2016, 43(4): 223-226.

[7] XU L L, WANG C, CAO W, ZHOU S M, WU T. CLAVATA1-type receptor-like kinase CsCLAVATA1 is a putative candidate gene for dwarf mutation in cucumber. Molecular Genetics and Genomics, 2018, 293(6): 1393-1405.

[8] GUNER N, WEHNER T C. Gene list for watermelon. Cucurbit Genetics Cooperative Report, 2003, 26: 76-92.

[9] HUANG H X, ZHANG X Q, WEI Z C, LI Q H, LI X. Inheritance of male-sterility and dwarfism in watermelon [Citrullus lanatus (Thunb.) Matsum and Nakai]. Scientica Horticulturae, 1998, 74(3): 175-181.

[10] 王深浩, 李海真, 张忠华, 贺俊, 贾长才, 张帆, 黄三文. 南瓜矮生基因Bu 的比较定位. 园艺学报, 2011, 38(1): 95-100. WANG S H, LI H Z, ZHANG Z H, HE J，JIA C C，ZHANG F，HUANG S W. Comparative mapping of the dwarf gene Bu from tropical pumpkin (Cucurbita moschata Duchesne). Acta Horticulturae Sinica, 2011, 38(1): 95-100. (in Chinese)

[11] WANG R, HUANG H X, LIN Y E, CHEN Q H, LIANG Z J, WU T Q. Genetic and gene expression analysis of dm1, a dwarf mutant from Cucurbita maxima Duch. ex Lam, based on the AFLP method. Canadian Journal of Plant Science, 2014, 94(2): 293-302.

[12] ZHANG G Y, REN Y, SUN H H, GUO S G, ZHANG F, ZHANG J, ZHANG H Y, JIA Z C, FEI Z J, XU Y, LI H Z. A high-density genetic map for anchoring genome sequences and identifying QTLs associated with dwarf vine in pumpkin (Cucurbita maxima Duch.). BMC Genomics, 2015, 16: 1101.

[13] PARIS H S, NERSON H, KARCHI Z. Genetics of internode length in melons. Journal of Heredity, 1984, 75(5): 403-407.

[14] KNAVEL D E. Inheritance of a short-internode mutant of ‘mainsteam’ muskmelon. Hortscience, 1990, 25(10): 1274-1275.

[15] HWANG J, OH J, KIM Z, STAUB J E, CHUNG S M, PARK Y. Fine genetic mapping of a locus controlling short internode length in melon (Cucumis melo L.). Molecular Breeding, 2014, 34(3): 949-961.

[16] 王建设, 张立杰, 唐晓伟, 宋曙辉. 甜瓜短蔓性状的遗传分析. 华北农学报, 2003, 18(4): 58-60. WANG J S, ZHANG L J, TANG X W, SONG S H. Inheritance of short stem character on Cucumis melo L. Acta Agriculturae Boreali- Sinica, 2003, 18(4): 58-60. (in Chinese)

[17] MURRAY M G, THOMPSON W F. Rapid isolation of high molecular weight plant DNA. Nucleic Acids Research, 1980, 8(19): 4321-4325.

[18] MA J, LEI C L, XU X T, HAO K, WANG J L, CHENG Z J, MA X D, MA J, ZHOU K N, ZHANG X, GUO X P, WU F Q, LIN Q B, WANG C M, ZHAI H Q, WANG H Y, WAN J M. Pi64, Encoding a novel CC-NBS-LRR protein, confers resistance to leaf and neck blast in rice.

Molecular Plant-Microbe Interactions, 2015, 28(5): 558-568.

[19] CHENG H, KUN W P, LIU D S, SU Y Q, HE Q W. Molecular cloning and expression analysis of CmMlo1 in melon. Molecular Biology Reports, 2012, 39(2): 1903-1907.

[20] SCHMITTGEN T D, LIVAK K J. Analyzing real-time PCR data by the comparative C(T) method. Nature Protocols, 2008, 3(6): 1101-1108.

[21] TORII K U, MITSUKAWA N, OOSUMI T, MATSUURA Y, YOKOYAMA R, WHITTIER R F, KOMEDA Y. The Arabidopsis ERECTA gene encodes a putative receptor protein kinase with extracellular leucine-rich repeats. The Plant Cell, 1996, 8(4): 735-746.

[22] ZINK F W. U. C. SR-91 Bush, U. C. Top Mark Bush and U. C. Perlita Bush muskmelon breeding lines. Hortscience, 1979, 13(4): 486.

[23] ZINK F W. U. C. Crenshaw bush breeding line. Hortscience, 1980, 15(3): 318-319.

[24] OH J Y. Genetic analysis and development of molecular markers of dwarf trait in melon (Cucumis melon L.) [D]. South Korea: Gyeongsang National University, 2013.

[25] 贾继增, 高丽锋, 赵光耀, 周文斌, 张卫健. 作物基因组学与作物科学革命. 中国农业科学, 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)

[26] BROTMAN Y, NORMANTOVICH M, GOLDENBERG Z, ZVIRIN Z, KOVALSKI I, STOVBUN N, DONIGER T, BOLGER A M, TROADEC C, BENDAHMANE A, COHEN R, KATZIR N, PITRAT M, DOGIMONT C, PERL-TREVES R. Dual Resistance of melon to Fusarium oxysporum races 0 and 2 and to Papaya ring-spot virus is controlled by a pair of head-to-head-oriented NB-LRR genes of unusual architecture. Molecular Plant, 2013, 6(1): 235-238.

[27] TZURI G, ZHOU X J, CHAYUT N, YUAN H, PORTNOY V, MEIR A, SA'AR U, BAUMKOLER F, MAZOUREK M, LEWINSOHN E, FEI Z J, SCHAFFER A A, LI L, BURGER J, KATZIR N, TADMOR Y. A ‘golden’ SNP in CmOr governs the fruit flesh color of melon (Cucumis melo). The Plant Journal, 2015, 82(2): 267-279.

[28] MIZUNO S, SONODA M, TAMURA Y, NISHINO E, SUZUKI H, SATO T, OIZUMI T. Chiba Tendril-Less locus determines tendril organ identity in melon (Cucumis melo L.) and potentially encodes a tendril-specific TCP homolog. Journal of Plant Research, 2015, 128(6): 941-951.

[29] NING X F, WANG X L, GAO X W, ZHANG Z Q, ZHANG L H, YAN W L, LI G. Inheritances and location of powdery mildew resistance gene in melon Edisto47. Euphytica, 2014, 195(3): 345-353.

[30] QU X Y, ZHAO Z, TIAN Z X. ERECTA regulates cell elongation by activating auxin biosynthesis in Arabidopsis thaliana. Frontiers in Plant Science, 2017, 8: 1688.

[31] KOSENTKA P Z, ZHANG L, SIMON Y A, SATPATHY B, MARADIAGA R, MITOUBSI O, SHPAK E D. Identification of critical functional residues of receptor-like kinase ERECTA. Journal of Experimental Botany, 2017, 68(7): 1507-1518.

Fine Mapping and Candidate Gene Analysis of a Short Internodes Gene Cmdm1 in Melon (Cucumis melo L.)

MA Jian1, LI CongCong2, WANG JianShe1
(1Beijing Vegetable Research Center, Beijing Academy of Agriculture and Forestry Sciences/Key Laboratory of Biology and Genetic Improvement of Horticultural Crops (North China), Ministry of Agriculture and Rural Affairs, Beijing 100097; 2Beijing Agro-Biotechnology Research Center, Beijing Academy of Agricultural and Forestry Sciences/Beijing Key Laboratory of Agricultural Genetic Resources and Biotechnology, Beijing 100097)

Abstract: 【Objective】In this study, the Z8 was a mutant with short internode in melon, and its gene fine mapping and candidate gene analysis were conducted to provide the theoretical guidance for the improvement of melon architecture. 【Method】Agronomic traits, such as internode number, main stem length, main stem internodes length and lateral branch length were measured in the Z8 and wild-type B15. The Z8 was crossed with B15, and the F2 and BC1F1 generations were used for genetic analysis. The recessive individuals in F2 were used to locate the gene by the map-based cloning strategy, and the candidate gene in the mapping region was determined by gene sequencing. 【Result】The Z8 mutant was dwarf because the internodes were significantly shorter compared with the wild-type B15. In addition, the flower buds of Z8 mutant were clustered at the top, which showed significantly compact inflorescences. Genetic test indicated that the mutant phenotype segregated as a single recessive gene locus named Cmdm1, which had been mapped to a 56 kb interval on the short arm of chromosome 7, flanked by markers c7-112 and s2 and co-segregated with marker dm-1. The MELO3C016916 gene was homologous to Arabidopsis ERECTA gene, which was one of the four annotated open reading frames (ORFs) in the 56 kb region. Sequence analysis of MELO3C016916 gene in the Z8 mutant indicated that a single nucleotide substitution (T1995→G1995) in its 25th exon caused premature termination of the encoded protein, thus placing MELO3C016916 as a candidate gene for Cmdm1.【Conclusion】The short internode trait of Z8 was controlled by a single recessive gene named Cmdm1, which was finally mapped to a 56 kb interval between markers c7-112 and s2 on chromosome 7, and the MELO3C016916 gene was presumed to be the candidate gene.

Key words: melon (Cucumis melo L.); short internodes; gene mapping; Cmdm1; ERECTA

收稿日期：2019-08-07；

接受日期：2019-09-16

基金项目：国家自然科学基金（31701937）、北京市农林科学院院青年科研基金（QNJJ201701）、蔬菜中心科研创新基金（KYCX201705-04）、北京市农林科学院科技创新能力建设专项（KJCX20170102-12）

联系方式：马建，E-mail：majian@nercv.org。

通信作者王建设，E-mail：wangjianshe@nercv.org

开放科学（资源服务）标识码（OSID）：

pagenumber_ebook=156,pagenumber_book=802
（责任编辑 赵伶俐）