浙江大学学报（农业与生命科学版）

浙江大学学报（农业与生命科学版）, 2021, 47(6): 719-728 doi: 10.3785/j.issn.1008-9209.2021.04.121

研究论文

柑橘果实大小与质量的遗传分析和数量性状位点定位

罗艾,,, 龚桂芝, 彭祝春, 杨程, 常珍珍, 洪棋斌,,

西南大学柑桔研究所/中国农业科学院柑桔研究所，重庆 400712

Genetic analysis and quantitative trait locus mapping of citrus fruit size and mass

LUO Ai,,, GONG Guizhi, PENG Zhuchun, YANG Cheng, CHANG Zhenzhen, HONG Qibin,,

Citrus Research Institute, Southwest University/Chinese Academy of Agricultural Sciences, Chongqing 400712, China

通讯作者: 洪棋斌（https://orcid.org/0000-0002-0039-4025），E-mail：hongqb@sina.com

收稿日期: 2021-04-12   接受日期: 2021-06-13  

基金资助: 国家重点研发计划.  2019YFD1001402
国家科技支撑计划.  2013BAD02B02

Received: 2021-04-12   Accepted: 2021-06-13  

作者简介 About authors

罗艾（https://orcid.org/0000-0001-6100-6279），E-mail：1012502243@qq.com , E-mail：1012502243@qq.com

摘要

为探究柑橘果实大小与质量的遗传调控机制，开展了果实大小和质量的数量性状位点（quantitative trait locus, QTL）定位分析，以挖掘调控基因。以‘晚蜜2号’和‘梨橙2号’的94株F1代材料为分离群体，采用保守同源序列（conserved ortholog sequence, COS）和微卫星标记即简单重复序列（simple sequence repeat, SSR）构建遗传连锁图谱，测定单株成熟果实大小和质量等性状，并进行遗传传递力分析和QTL定位分析。结果表明：本研究构建了一张包含201个标记、10个连锁群，总长为1 194.5 cM的柑橘遗传连锁图谱；定位到与果实质量相关的4个QTLs，与横径相关的3个QTLs，与纵径相关的4个QTLs，都分别位于WL3和WL8连锁群上。通过对这些区域内的编码基因进行挖掘，发现了一些与果实发育相关的候选基因：GATA转录因子、生长素响应因子（auxin response factor, ARF）和GDSL酯酶。本研究鉴定到的QTLs和候选基因可为柑橘果实大小和质量的遗传研究、育种和物种改良提供理论参考和分子标记工具。

关键词： 柑橘 ; 果实大小 ; 果实质量 ; 遗传图谱 ; 数量性状位点

Abstract

In order to explore the genetic regulation mechanism of citrus fruit size and mass, mapping analysis of quantitative trait loci (QTLs) of fruit size and mass was carried out to find regulatory genes. The segregating population in this study was the F1 generation with a total of 94 plants, from the cross of ‘Wanmi 2’ [Citrus unshiu (Mark.) Marc.×C. sinensis (L.) Osb.] and a local sweet orange cultivar ‘Licheng 2’ [C. sinensis (L.) Osb.]. Conserved ortholog sequence (COS) and simple sequence repeat (SSR) markers were used to construct a genetic linkage map of the segregating population. Fruit size and mass were measured at the full maturing stage. Genetic transmitting ability analysis and QTL mapping were carried out. The results showed that this study constructed a genetic linkage map, which contained 201 markers, 10 linkage groups, with a total length of 1 194.5 cM. Four QTLs related to citrus fruit mass, three QTLs related to citrus transverse diameter, and four QTLs related to citrus longitudinal diameter were detected. All the QTLs were distributed on WL3 and WL8 linkage groups. A batch of functional genes (GATA transcription factor, auxin response factor and GDSL esterase) related to the development of citrus fruits were identified according to the functional annotations of the genes in the genome defined by the markers. The QTLs located and the candidate genes identified in this study will be helpful for citrus breeding and theoretical research.

Keywords： citrus ; fruit size ; fruit mass ; genetic map ; quantitative trait locus (QTL)

PDF (4098KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

罗艾, 龚桂芝, 彭祝春, 杨程, 常珍珍, 洪棋斌. 柑橘果实大小与质量的遗传分析和数量性状位点定位. 浙江大学学报（农业与生命科学版）[J]. 2021, 47(6): 719-728 doi:10.3785/j.issn.1008-9209.2021.04.121

LUO Ai, GONG Guizhi, PENG Zhuchun, YANG Cheng, CHANG Zhenzhen, HONG Qibin. Genetic analysis and quantitative trait locus mapping of citrus fruit size and mass. Journal of Zhejiang University（Agriculture & Life Sciences）[J]. 2021, 47(6): 719-728 doi:10.3785/j.issn.1008-9209.2021.04.121

柑橘（Citrus L.）属于芸香科（Rutaceae）柑橘亚科（Aurantiodeae），在我国的栽种历史久远，栽培种类众多，同时也是国内栽培面积和产量最大的水果，具有重要的经济价值[1-2]。柑橘因富含维生素C、类胡萝卜素和黄酮类化合物等功能成分而具有较高的健康属性，且储藏运输方便、加工特性突出，是全世界重要的水果商品，其生产和消费量一直在不断增加[3]。近年来，随着社会经济的快速增长和消费水平的提高，人们对柑橘产品的要求也在不断改变和提高，对柑橘育种提出了更高要求。但柑橘是木本植物，具有树体大、童期长、遗传高度杂合等生物学特性，再加上一些品种的无融合生殖特性、雌性败育和/或雄性败育的存在，影响了研究者对其性状，尤其是数量性状的遗传研究，极大阻碍了柑橘育种产业发展[4-5]。分子标记辅助选择（marker-assisted selection, MAS）以及转基因育种（genetically modified breeding, GMB）技术在植物研究中的广泛应用，为加快柑橘等长周期的果树育种带来了希望[6]。针对柑橘的品质和产量等重要农艺性状开展研究，发掘其相关的调控基因或相关标记，成为利用新技术提高柑橘育种效率的重要基础。

果实大小是柑橘的重要农艺性状，不仅影响柑橘产量和外观品质，也影响生产成本和消费者的选择，过小的果实采摘成本高，过大的果实消费者一次性食用不完，影响消费者的购买欲望，因此，日本柑橘育种者在宽皮柑橘育种中，提出了果实大小介于150～300 g的育种目标[7]。果实大小相关性状是由多基因控制的数量性状，受到内外多种因素的影响，研究难度不小，目前关于其遗传调控的机制尚不清楚，更缺乏可供育种利用的调控基因或分子标记。因此，本研究基于培育多年的‘晚蜜2号’和‘梨橙2号’杂交的F1代分离群体，借助果实大小和质量的数量性状位点（quantitative trait locus, QTL）定位，发掘相关的调控基因或相关标记，进而为柑橘果实大小和质量的遗传改良提供分子辅助育种工具和理论参考。

1 材料与方法

1.1 试验材料

分离群体是以‘晚蜜2号’[Citrus unshiu (Mark.) Marc.×C. sinensis (L.) Osb.]为母本、‘梨橙2号’[C. sinensis (L.) Osb.]为父本杂交所得到的94株F1代材料。‘晚蜜2号’为尾张温州蜜柑和细叶薄皮甜橙杂交所得，较抗寒，雄性部分败育，果实外观呈扁圆形，果皮橙红，种子数目较少，单胚。‘梨橙2号’为甜橙品种，果实外观呈梨形，果皮橙色，种子少，单胚。

分离群体材料和亲本均采用飞龙枳作砧木，区组重复，2012年定植在中国农业科学院柑桔研究所育种圃砧木比较试验园，到2021年完成测试时已是11年生树，但仍有部分材料未正常结果，本研究选择结果较多的材料，总计61株（2株亲本和59株杂交后代）进行果实性状的观测。

1.2 试验方法

1.2.1 分子标记分析和遗传图谱的构建

提取分离群体和亲本的DNA以开展分子标记检测。DNA提取、保守同源序列（conserved ortholog sequence, COS）和微卫星标记即简单重复序列（simple sequence repeat, SSR）分析参考王炯[8]的方法。采用JoinMap® 4.0软件[9]进行分离群体的遗传连锁分析。

1.2.2 柑橘果实大小和质量测定及遗传分析

于2020年柑橘果实成熟期，每株选取6个果实，分别使用电子天平和游标卡尺对单果质量、纵横径进行测量。每个果实重复测定3次，最后取平均值作为后续试验分析。亲本单果质量和纵横径的遗传传递力（Ta）计算公式[10]为：Ta=杂交后代性状均值/双亲性状均值×100%。

1.2.3 数量性状位点（QTL）分析

利用SPSS 22.0软件对果实纵横径和质量数据进行统计分析，绘制频率分布直方图。基于构建的连锁图谱，采用MapQTL® 6.0[11]软件对果实的单果质量和纵横径进行QTL定位分析。

2 结果与分析

2.1 柑橘遗传连锁图谱的构建结果

本试验中，杂交亲本均为遗传学上高度杂合的材料，符合JoinMap® 4.0软件的“Cross Pollinator, CP”模式，CP模式杂交后代会表现出5种分离类型：ab×cd、ef×eg、hk×hk、lm×ll、nn×np。对本试验中杂交分离群体的电泳结果进行解读分析，可以比较准确地辨识出后4种分离类型，但未发现“ab×cd”类型。不同分离类型的代表性样品电泳结果如图1所示。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1   不同分离类型的代表性扩增检测结果

W为母本，L为父本，数字1～94分别为W×L的F1代。

Fig. 1   Representative amplification testing results of different segregation types

W is the female parent, and L is the male parent, and numbers of 1-94 are the F1 generation of W×L, respectively. CSSRC2-14: lm×ll; COSC4-10: nn×np; COSC8-10: hk×hk; CSSRC2-13: ef×eg.


研究发现，在所有标记中共有108个（约54%）标记表现出明显的偏分离现象，而且大多偏向母本 （P＞0.05，表1）。标记偏分离比率大于CURTOLO等[12]和王炯[8]构建的柑橘遗传连锁图谱的偏分离比率（43%和45%，P=0.05），小于DE OLIVEIRA等[13]构建的柑橘连锁图谱中的偏分离比率（61%，P=0.05）。

表1   标记的分离类型

Table 1  Segregation types of all markers

分离类型

Segregation type

标记数目

No. of markers

期望比率

Expected ratio

偏分离位点数目

No. of skewed loci*

偏分离位点方向（数目）

Direction of skewed loci (No.)

合计 Total201108108
ef×eg601∶1∶1∶149偏母本（40）＋偏父本（9）
hk×hk571∶2∶145偏杂合体（45）
lm×ll201∶116偏父本（16）
nn×np641∶128偏母本（13）＋偏父本（15）

*表示在P＜0.05水平差异有统计学意义的偏分离位点数目。

Single asterisk (*) indicates significant differences of No. of skewed loci at the 0.05 probability level.

新窗口打开| 下载CSV



采用比较严格的作图标准似然函数比值的对数[logarithm (base) of odds, LOD]构建遗传连锁图谱（LOD=6），并将少数未入群的标记按“最强关联组（strong cross linkage group, SCL-Group）”进行归群后，得到本研究的柑橘遗传连锁图谱。该图谱包含201个COS和SSR标记、10个连锁群，总长度为 1 194.5 cM，标记间平均距离约为5.9 cM。分析发现，连锁群的主位标记与测序C. clementina（克里迈丁）的Scaffold顺序号具有一致性。以亲本‘晚蜜2号’和‘梨橙2号’的拼音首字母缩写及测序C. clementina 的Scaffold顺序号将连锁群命名为WL1～WL9；而COSC9系列标记在标准作图连锁过程中不能自然整合到同一连锁群，而是分成2个连锁群，本文将其分别命名为WL9-1和WL9-2（图2）。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2   基于COS和SSR标记构建的连锁图谱

Fig. 2   Linkage maps constructed by COS and SSR markers


通过与其他已发表的柑橘遗传连锁图谱[3,7-8,12]进行对比，发现分离群体、标记类型及数量对构建的图谱影响很大。RUIZ等[14]利用随机扩增多态性DNA（random amplified polymorphic DNA, RAPD）、限制性片段长度多态性（restriction fragment length polymorphism, RFLP）和其他一些分子标记技术构建了柑橘遗传图谱，但其所构建的连锁群数与柑橘染色体数目不一致，同时还存在图距大、基因组覆盖率低等问题。近年来，研究者希望采用单核苷酸多态性（single nucleotide polymorphism, SNP）标记构建高密度遗传连锁图，也同样出现了连锁群数与染色体数不一致的情况，而且还存在分型单一、后续分析验证较难等问题。王炯[8]在2017年构建了含240个COS和SSR标记、9个连锁群，总长为1 186.4 cM的遗传连锁图谱，但其连锁群1的密度很低。本研究采用其反交材料构建图谱，构建了一张包含201个标记、10个连锁群，总长为1 194.5 cM的柑橘遗传连锁图谱，增加了连锁群1的密度，其中多数相同的标记仍然定位到同一连锁群上，少数标记在连锁群上的排列顺序有差异或者定位到了不同的连锁群上（图2）。

2.2 柑橘果实大小和质量的相关性及遗传分析

对分离后代的果实质量、横纵径的均值进行比较，发现其均低于亲中值和亲本值，表现为趋小变异。亲本果实大小和质量的遗传传递力由大到小排列为横径＞纵径＞果实质量，其中果实质量的遗传传递力为74%，横径的遗传传递力为93%（表2）。这3个性状的变异大，均出现明显高于或低于亲本个体的情况，表明控制柑橘果实大小和质量的增减效基因在双亲中分散分布，因此，可以通过基因重组等方式培育超亲后代。

表2   亲本及分离群体的果实质量和横纵径分析

Table 2  Analysis of the fruit mass, transverse and longitudinal diameters of parents and the segregation populations

性状

Trait

‘晚蜜2号’

‘Wanmi 2’

‘梨橙2号’

‘Licheng 2’

中值

Mid-value

平均值

Mean

标准差

Standard

deviation

范围

Range

遗传传递力

Ta/%

偏离度系数

Skewness

coefficient

峰度系数

Kurtosis

coefficient

果实质量 Fruit mass/g194.05238.60216.33161.005.2782.90～266.02741.040.46

横径

Transverse diameter/cm

8.087.577.837.250.105.41～9.52930.271.68

纵径

Longitudinal diameter/cm

6.538.507.516.280.104.56～8.56841.040.65

新窗口打开| 下载CSV


对测试的各性状进行相关性分析发现，果实质量与横纵径均呈极显著正相关，相关系数分别为0.824和0.835，说明果实质量受纵径影响更大。横径和纵径也呈极显著正相关，相关系数为0.756（表3）。该结果与前人研究报道[12]一致。

表3   分离群体果实大小和质量的相关性分析

Table 3  Correlation analysis of fruit size and mass in the segregation populations

性状

Trait

横径

Transverse

diameter

纵径

Longitudinal

diameter

果实质量 Fruit mass0.824**0.835**
横径 Transverse diameter0.756**

**表示在P＜0.01水平极显著相关。

Double asterisks (**) indicate extremely significant correlations at the 0.01 probability level.

新窗口打开| 下载CSV


2.3 柑橘果实大小和质量的数量性状位点（QTL）定位分析

对分离群体的果实质量、横径、纵径数据进行分析，发现三者均为连续正态分布（图3），果实质量、横径和纵径的偏离度系数分别为1.04、0.27和1.04，峰度系数分别为0.46、1.68和0.65（表2），绝对值均低于2，因此，可以对其进行QTL定位。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3   果实质量和横纵径在分离群体中的频率分布

Fig. 3   Frequency distributions of fruit mass, transverse and longitudinal diameters in the segregation populations


基于构建的连锁图谱，采用MapQTL® 6.0软件的区间作图法对柑橘果实的质量、横径、纵径进行QTL定位，检测到与单果质量相关的4个QTL（LOD≥4.39），与果实横径相关的3个QTL（LOD≥3.08），与纵径相关的4个QTL（LOD≥3.70）（表4）。由图4可见：这些QTLs均位于WL3和WL8连锁群上，表现出集聚性；从标记界定的QTL区间来看，果实质量与纵径在连锁群WL8上有重叠的基因组区间（CSSRC9-18～COSC8-6），也表现出一致性。这可能与果实质量、横径、纵径三者之间具有极显著正相关关系（表3）有关。

表4   柑橘果实大小和质量的QTLs分析

Table 4  QTLs analysis of citrus fruit size and mass

性状

Trait

数量性状位点

QTL

连锁群

Linkage group

位置

Position/cM

标记区间

Marked interval

LOD

贡献率

Contribution rate/%

果实质量

Fruit mass (FM)

qtl3.1 (FM)WL339.905COSC3-9～LcSSR1305.9434.4
LcSSR130 (FM)WL340.176LcSSR1305.8834.1
qtl3.2 (FM)WL342.117COSC3-1～COSC3-335.0430.0
qtl8.1 (FM)WL835.563CSSRC9-18～COSC8-64.3926.7

横径

Transverse diameter (TD)

LcSSR120 (TD)WL80.000LcSSR1203.9124.2
qtl8.1 (TD)WL829.405COSC8-8～CSSRC9-183.5322.1
COSC8-8 (TD)WL828.405COSC8-83.0819.6

纵径

Longitudinal diameter (LD)

qtl8.1 (LD)WL834.563CSSRC9-18～COSC8-63.8824.0
qtl8.2 (LD)WL837.447COSC8-6～COSC9-123.8223.7
qtl8.3 (LD)WL835.563CSSRC9-18～COSC8-63.8023.6
qtl8.4 (LD)WL829.405COSC8-8～CSSRC9-183.7023.0

新窗口打开| 下载CSV


图4

新窗口打开| 下载原图ZIP| 生成PPT

图4   柑橘果实大小和质量QTLs在连锁图谱上的定位

Fig. 4   QTLs locations of citrus fruit size and mass in the linkage map


2.4 柑橘果实大小和质量QTL区域的功能基因挖掘

以QTL所在的分子标记或区间两侧的分子标记序列选择C. clementina为参考基因组，在Phytozome（https://phytozome-next.jgi.doe.gov/pz/portal.html）网站进行Blast分析，对标记所界定的基因组区间内的编码基因进行功能注释查询，初步挖掘到一批可能与柑橘果实发育有关的功能基因，部分结果如表5所示。进一步选择LOD值和贡献率较高，或在果实质量和纵横径QTL定位中均出现的候选区域开展进一步的基因挖掘分析。qtl8.1（FM）的LOD值为4.39，qtl8.1（LD）的LOD值为3.88，均定位于CSSRC9-18～COSC8-6界定的区域，对该区域的序列进行分析发现：在参考基因组中临近COSC8-6引物序列的上游，有注释功能为GATA转录因子的功能基因；在临近COSC8-6引物序列的下游，有注释功能为生长素响应因子（auxin response factor, ARF）的功能基因。GATA为一类植物转录因子，能与靶基因启动子上的WGATAR区域结合，激活或抑制下游靶基因的转录，从而调控植物的发育[15]。ARF是由ULMASOV等[16]新发现的一类调控生长素响应基因表达的转录因子家族，能够特异地与生长素响应基因启动子区域的生长素响应元件结合。梅丽华[17]研究发现，番茄ARF家族成员之一的SlARF10通过控制淀粉合成相关酶的表达来促进糖类物质的积累，在番茄果实发育过程中发挥着重要作用。说明ARF可能具有一因多效的性质，对果实质量、横纵径有共同调控作用。本试验中，LcSSR130（FM）的LOD值为5.88，LcSSR130是基于表达序列标签（expressed sequence tag, EST）开发的引物，该EST所在的基因功能注释为GDSL酯酶（GDSL esterase）。有研究报道，来自水稻脆叶鞘1突变体BS1的GDSL酯酶成员是一种次级壁乙酰-木聚糖酯酶，能调控木聚糖的脱乙酰化，进而导致次生壁发育受损和异常的多变性表型[18]。柑橘果实中的GDSL酯酶成员也很可能通过动态微调多糖的乙酰化结构来发挥结构调控作用，从而调节细胞壁结构，进而对果实的大小和质量产生影响。

表5   与柑橘果实发育可能相关的基因

Table 5  Genes putatively associated with citrus fruit development

基因序列号

GenBank No.

注释

Description

克里迈丁位置

Scaffold on C. clementina

Ciclev10018887m过渡性内质网腺苷三磷酸酶（ATPase），是一种保守的ATPase，参与多种细胞过程[19]Scaffold_3:3367083..3371823
Ciclev10018599m富含亮氨酸的重复序列，能够为蛋白质-蛋白质相互作用的形成提供一个通用的结构框架[20]Scaffold_3:3372883..3377278
Ciclev10023159m阿拉伯半乳聚糖肽，是一种可溶性膳食纤维，由与15个残基肽相连的阿拉伯半乳聚糖多糖组成[21]Scaffold_3:1160992..1161777
Ciclev10029679m线粒体导入内膜转位亚基TIM9，在线粒体膜蛋白的导入过程中起着重要的伴侣蛋白样作用[22]Scaffold_8:23446258..23448342
Ciclev10028134m重复序列蛋白，参与植物发育、细胞器生物发生、线粒体和叶绿体中的RNA加工和编辑、质体到核逆行信号传导以及对环境压力的反应[23]Scaffold_8:23439132..23442638
Ciclev10029015mGATA转录因子，能与靶基因启动子上的WGATAR区域结合，激活或抑制下游靶基因的转录，从而调控植物的发育[15]Scaffold_8:21875917..21880999
Ciclev10027901m生长素响应因子，能够特异地与生长素响应基因启动子区域的生长素响应元件结合[16]Scaffold_8:21887372..21893128
Ciclev10028912m锌指蛋白，在许多植物蛋白中高度保守，与类固醇激素受体和GATA1显著相关[24]Scaffold_8:22487087..22489505
Ciclev10020515mGDSL酯酶，能调控木聚糖的脱乙酰化，进而导致次生壁发育受损和异常的多变性表型[18]Scaffold_3:40426486..40429683

新窗口打开| 下载CSV


3 讨论与结论

3.1 偏分离标记的作图应用

孟德尔分离定律是经典遗传学的基本原则，能帮助遗传学家预测简单遗传性状的表达情况。人们从20世纪初就开始研究偏离孟德尔分离预期的现象，即偏分离现象，主要是配子或合子选择、种间或种内杂交引起染色体重排而导致的。在物种的大多数遗传作图研究中都观察到了偏分离标记[8,10]，其可能会使标记之间遗传距离的估计出现偏差，并可能影响标记之间的顺序，从而影响表型性状QTL定位的准确性。因此，在遗传作图中常有将偏分离标记舍弃后再进行图谱构建的报道[3,7]。DE OLIVEIRA等[25]通过对仅包含正常分离标记图谱和包含偏分离标记的图谱进行比较发现，偏分离标记一般位于染色体的末端部位，并未显著改变标记的距离。同时，随着对偏分离的深入研究，人们发现偏分离标记对QTL定位中的加性效应影响较小，并且偏分离标记可能与含生物关键遗传特性的染色体区域相关[26]，如果将其排除在遗传图谱之外，则无法检测到一些与生长和成熟相关的生物学功能。偏分离标记大部分是聚集的，并且每个聚集内的标记具有相同的偏斜方向，比如多数位点偏向于同一亲本或杂合体，这表明偏分离标记是由生物因素引起的[27]。本研究的偏分离标记大多偏向于母本，说明来自母本基因型的配子在形成合子时占主导地位。同时，本研究所定位到的QTL相关的分子标记，比如LcSSR130（FM）、LcSSR120（TD）和COSC8-8（TD）等大多是偏分离标记，因此将其用于图谱的构建必不可少。目前，已有关于偏分离分析模型与软件的报道，说明大部分的偏分离标记确实可以用于遗传图谱构建，只是需根据实际情况进行评估和调整[27]。

3.2 柑橘果实大小和质量的QTL比较

本研究采用的COS标记是根据测序克里迈丁从头（de novo）组装结果（https://phytozome.jgi.doe.gov/pz/portal.html）开发的。因此，得到的连锁群与克里迈丁参考基因组的染色体顺序具有一致性，可以比较全面地与其他研究者基于一致的连锁群定义而开展的柑橘果实大小和质量的QTL定位进行比较分析。以克里迈丁为参考基因组，GARCÍA等[28]检测到与平均果实大小相关的3个QTLs，分别分布在V2、V7、V8连锁群上；检测到与平均单果质量相关的2个QTLs，分布在V2和V7连锁群上。但其报道中并未表明这些位点的LOD值和贡献率。YU等[3]结合‘Fortune’בMurcott’（Citrus reticulata Blanco）的F1代检测到3个与柑橘横径相关的QTLs，3个与单果质量相关的QTLs，这些位点分布在MUR4.2、MUR8、FOR5.1和FOR9.3连锁群上，其中只有位于连锁群FOR5.1上的FM5.1和位于MUR8上的FM8在2012年和2013年均有被检测到。IMAI等[7]定位到4个与柑橘单果质量相关的QTLs，分布在LG3、LG4和LG7连锁群上，其中只有位于LG4上的位点FMq3在连续2年里被检测到。本研究定位到的与果实大小相关的QTLs分布在WL3和WL8连锁群上。综上所述，在柑橘第3、4和8染色体上可能存在控制果实大小的QTLs，而在第8染色体上可能存在控制柑橘果实大小的相关基因，在第3染色体上可能主要存在与果实质量相关的功能基因。

YU等[3]和IMAI等[7]的研究中，与果实质量和横径相关的QTLs的LOD值在2.69～3.68之间，贡献率在14.9%～33.74%之间；本研究中与果实大小和质量相关的QTLs的LOD值在3.08～5.94之间，贡献率在19.6%～34.4%之间，研究结果与前人相近。CURTOLO等[12]研究发现，与柑橘果实横纵径相关的QTLs的LOD值虽然都大于5.94，但是其表型变异率最高只有9.7%；GARCÍA等[28]却没有报道其定位到的QTLs的LOD值和解释表型变异的决定系数，因此，其定位到的QTLs还有待考证。BUDAHN等[29]研究表明，超过45%的表型变异是由主要的QTLs效应造成的，而这在本文和以上研究中均没有观察到。但也有报道指出贡献率大于30%即可认为是主效QTLs[30]。因此，对于主效QTLs的判断标准还有待进一步探讨，但这也说明果实大小和质量是由多基因共同控制的。

基因型×环境互作效应对柑橘果实品质性状起着重要作用，基因型×取样日期对所有柑橘果实品质性状也有显著影响，特别是柑橘的外部特征比内部特征波动更大。YU等[3]使用2年的表型数据检测到与果实直径和果实质量相关的QTLs有6个，在该研究中，2年里取样时间和表型材料并不完全相同，所以只有关于果实质量的2个QTLs在2年里被检测到，而与果实直径相关的QTLs在2年里均不相同。IMAI等[7]采用3年的表型数据定位到与果实质量相关的QTLs位点有4个，且只有1个位点在2年里被连续检测到。CURTOLO等[12]结合1年的表型数据，只检测到3个与果实纵横径相关的QTLs。与以上研究比较，我们对分离群体的果实发育全程进行了分析测试，发现到发育后期，即开花后246—274 d，果实的生长量变化很小，果实大小总体已进入比较稳定的状态；此外，本研究仅选用果实发育停滞时段的测试数据进行分析，减少了数据波动对QTL定位的影响。本研究定位到了与果实大小相关的10个QTLs，数量较多，这可能与我们所研究的分离群体果实大小和质量分离比较大有关，如在本研究中单果质量的变化范围为82.90～266.02 g，变异范围较大。分离群体中尚有部分未结果，连锁图谱存在标记数目不足，连锁群间和群内分布不均匀等，都对本研究的QTL数量确定和精准鉴定等带来了较大影响。后续可通过开发和应用更多分子标记优化遗传图谱，对更多结果个体进行分析以得到更好的QTL鉴定结果。在对LOD值和贡献率较高，或在果实质量和纵横径QTL定位中均出现的候选区域进行基因挖掘时，发现本研究所定位到的QTL所在的基因组区间确实存在对果实发育和大小产生影响的基因，这些基因对柑橘果实大小和质量的遗传和育种的深入研究和应用探索具有一定的参考价值。

参考文献

周常勇.中国果树科学与实践:柑橘.西安:陕西科学技术出版社,2020:34-38.

[本文引用: 1]

ZHOU C Y. Chinese Fruit Science and Practice: Citrus. Xi’an: Shaanxi Science and Technology Press,2020:34-38. (in Chinese)

[本文引用: 1]

杨晓明.柑橘亚科植物系统发育基因组学及野生枸橼、宜昌橙谱系地理学研究.武汉:华中农业大学,2017. DOI:10.1007/s11295-017-1113-4

[本文引用: 1]

YANG X M. The phylogenomic of aurantioide and phylogeography of wild citron(Citrus medica) and Ichang papeda (Citrus ichangensis). Wuhan: Huazhong Agricultural University, 2017. (in Chinese with English abstract)

DOI:10.1007/s11295-017-1113-4      [本文引用: 1]

YU Y, CHEN C X, GMITTER F G, Jr, et al.

QTL mapping of mandarin (Citrus reticulata) fruit characters using high-throughput SNP markers

Tree Genetics & Genomes, 2016,12(4):77. DOI:10.1007/s11295-016-1034-7

[本文引用: 6]

邓秀新.

柑橘产业发展趋势与桂林柑橘品种结构调整

.南方园艺,2020,31(6):1-4. DOI:10.3969/j.issn.1674-5868.2020.06.001

[本文引用: 1]

DENG X X.

The development trend of citrus industry and the adjustment of citrus variety structure in Guilin

Southern Horticulture, 2020,31(6):1-4. (in Chinese with English abstract)

DOI:10.3969/j.issn.1674-5868.2020.06.001      [本文引用: 1]

RUIZ C, BRETO M P, ASÍNS M J.

A quick methodology to identify sexual seedlings in citrus breeding programs using SSR markers

Euphytica, 2000,112(1):89-94. DOI:10.1023/A:1003992719598

[本文引用: 1]

QAIM M.

Role of new plant breeding technologies for food security and sustainable agricultural development

Applied Economic Perspectives and Policy, 2020,42(2):129-150. DOI:10.1002/aepp.13044

[本文引用: 1]

IMAI A, YOSHIOKA T, HAYASHI T.

Quantitative trait locus (QTL) analysis of fruit-quality traits for mandarin breeding in Japan

Tree Genetics & Genomes, 2017,13(4):79. DOI:10.1007/s11295-017-1162-8

[本文引用: 6]

王炯.基于COS Marker构建柑橘连锁图谱及作图群体的光合特性研究.重庆:西南大学,2017. DOI:10.25165/j.ijabe.20181102.3189

[本文引用: 5]

WANG J. Citrus linkage map construction based on COS marker and photosynthetic characteristics research of the mapping population. Chongqing: Southwest University, 2017. (in Chinese with English abstract)

DOI:10.25165/j.ijabe.20181102.3189      [本文引用: 5]

OOIJEN J W VAN.

JoinMap® 4

.0: software for the calculation of genetic linkage maps in experimental populations. Wageningen, the Netherlands: Kyazma B V, 2006.

[本文引用: 1]

陈克玲,陈力耕.

柑橘果形遗传的研究

.西南农业大学学报,1994,16(2):120-123. DOI:10.3390/plants9020196

[本文引用: 2]

CHEN K L, CHEN L G.

The heredity of fruit shape in citrus

Journal of Southwest Agricultural University, 1994,16(2):120-123. (in Chinese with English abstract)

DOI:10.3390/plants9020196      [本文引用: 2]

OOIJEN J W VAN.

MapQTL® 6

.0: software for the mapping of quantitative trait loci in experimental populations of diploid species. Wageningen, the Netherlands: Kyazma B V, 2009.

[本文引用: 1]

CURTOLO M, CRISTOFANI-YALY M, GAZAFFI R, et al.

QTL mapping for fruit quality in Citrus using DArTseq markers

BMC Genomics, 2017,18:289. DOI:10.1186/s12864-017-3629-2

[本文引用: 5]

DE OLIVEIRA A C, BASTIANEL M, CRISTOFANI-YALY M, et al.

Development of genetic maps of the citrus varieties ‘Murcott’ tangor and ‘Pêra’ sweet orange by using fluorescent AFLP markers

Journal of Applied Genetics, 2007,48:219-231. DOI:10.1007/BF03195216

[本文引用: 1]

RUIZ C, ASINS M J.

Comparison between Poncirus and Citrus genetic linkage maps

Theoretical and Applied Genetics, 2003,106:826-836. DOI:10.1007/s00122-002-1095-x

[本文引用: 1]

王娟,兰海燕.

GATA转录因子对植物发育和胁迫响应调控的研究进展

.植物生理学报,2016,52(12):1785-1794. DOI:10.13592/j.cnki.ppj.2016.0259

[本文引用: 2]

WANG J, LAN H Y.

Advances in regulation of GATA transcription factor to plant development and stress responses

Plant Physiology Journal, 2016,52(12):1785-1794. (in Chinese with English abstract)

DOI:10.13592/j.cnki.ppj.2016.0259      [本文引用: 2]

ULMASOV T, HAGEN G, GUILFOYLE T J.

ARF1, a transcription factor that binds to auxin response elements

Science, 1997,276(5320):1865-1868. DOI:10.1126/science.276.5320.1865

[本文引用: 2]

梅丽华.生长素响应因子基因SlARF10在番茄果实发育过程中的功能研究.重庆:重庆大学,2017. DOI:10.1101/253237

[本文引用: 1]

MEI L H. Functional studies of auxin response factor gene SlARF10 in tomato fruit development. Chongqing: Chongqing University, 2017. (in Chinese with English abstract)

DOI:10.1101/253237      [本文引用: 1]

ZHANG B C, ZHANG L J, FENG L, et al.

Control of secondary cell wall patterning involves xylan deacetylation by a GDSL esterase

Nature Plants, 2017,3(3):17017. DOI:10.1038/nplants.2017.17

[本文引用: 2]

FANG C Y, ZHANG X Q, ZHANG L, et al.

Identification of palmitoylated transitional endoplasmic reticulum ATPase by proteomic technique and pan antipalmitoylation antibody

Journal of Proteome Research, 2016,15(3):956-962. DOI:10.1021/acs.jproteome.5b00979

[本文引用: 1]

KOBE B, KAJAVA A V.

The leucine-rich repeat as a protein recognition motif

Current Opinion in Structural Biology, 2001,11(6):725-732. DOI:10.1016/S0959-440X(01)00266-4

[本文引用: 1]

HARRIS S, POWERS S, MONTEAGUDO-MERA A, et al.

Determination of the prebiotic activity of wheat arabino-galactan peptide (AGP) using batch culture fermentation

European Journal of Nutrition, 2020,59(1):297-307. DOI:10.1007/s00394-019-01908-7

[本文引用: 1]

SPILLER M P, GUO L, WANG Q, et al.

Mitochondrial Tim9 protects Tim10 from degradation by the protease Yme1

Bioscience Reports, 2015,35(3):e00193. DOI:10.1042/BSR20150038

[本文引用: 1]

范敏,金黎平,刘庆昌,等.

马铃薯PPR蛋白家族基因SoDIPPR的克隆及其在干旱条件下的表达特征分析

.中国农业科学,2008,41(8):2249-2257. DOI:10.3864/j.issn.0578-1752.2008.08.005

[本文引用: 1]

FAN M, JIN L P, LIU Q C, et al.

Cloning of SoDIPPR gene of pentatricopeptide repeat (PPR) protein family in potato and analysis of expression characteristics under drought conditions

Scientia Agricultura Sinica, 2008,41(8):2249-2257. (in Chinese with English abstract)

DOI:10.3864/j.issn.0578-1752.2008.08.005      [本文引用: 1]

SANDLER L, HIRAIZUMI Y, SANDLER I.

Meiotic drive in natural populations of Drosophila melanogaster

Ⅰ. The cytogenetic basis of segregation-distortion. Genetics, 1959,44(2):233-241. DOI:10.1093/genetics/44.2.233

[本文引用: 1]

DE OLIVEIRA R P, CRISTOFANI M, MACHADO M A, et al.

Genetic linkage maps of ‘Pêra’ sweet orange and ‘Cravo’ mandarin with RAPD markers

Pesquisa Agropecuária Brasileira, 2004,39(2):159-165. DOI:10.1590/S0100-204X2004000200009

[本文引用: 1]

ZHAO J L, HAN D D, SHI K T, et al.

Influence of epistatic segregation distortion loci on genetic marker linkages in Japanese flounder

Genomics, 2018,110(1):59-66. DOI:10.1016/j.ygeno.2017.08.006

[本文引用: 1]

SONG X L, SUN X Z, ZHANG T Z.

Segregation distortion and its effect on genetic mapping in plants

Chinese Journal of Agricultural Biotechnology, 2006,3(3):163-169. DOI:10.1079/CJB2006110

[本文引用: 2]

GARCÍA M, ASÍNS M J, CARBONELL E A.

QTL analysis of yield and seed number in Citrus

Theoretical & Applied Genetics, 2000,101(3):487-493. DOI:10.1007/s001220051507

[本文引用: 2]

BUDAHN H, PETERKA H, MOUSA M A A, et al.

Molecular mapping in oil radish (Raphanus sativus L.) and QTL analysis of resistance against beet cyst nematode (Heterodera schachtii)

Theoretical & Applied Genetics, 2009,118(4):775-782. DOI:10.1007/s00122-008-0937-6

[本文引用: 1]

阮成江,何祯祥,钦佩.

我国农作物QTL定位研究的现状和进展

.植物学通报,2003,20(1):10-22. DOI:10.3969/j.issn.1674-3466.2003.01.002

[本文引用: 1]

RUAN C J, HE Z X, QIN P.

Research advancements on crop QTL mapping in China

Chinese Bulletin of Botany, 2003,20(1):10-22. (in Chinese with English abstract)

DOI:10.3969/j.issn.1674-3466.2003.01.002      [本文引用: 1]

/

〈 〉