0 引言
【研究意义】中国是世界油菜生产大国,种植面积7×106 hm2,总产1.1×107 t,均占全世界的三分之一（https://www.chinaagrisci.com/article/2018/0578-1752/www.fao.org）。油菜生产在中国农业中有着举足轻重的作用。然而,随着中国农村劳动力成本增加,油菜种植效益和种植面积逐年下降[1],急需筛选和培育具有早熟、矮秆、株型紧凑、分枝短、熟期整齐度好等特性的品种,以适应油菜机械化生产的需要。由于油菜驯化历史短,变异类型有限,很大程度限制了油菜的遗传改良。因此,发现和创制新的油菜变异类型对拓宽现有油菜遗传资源、加快选育适合机械化作业的油菜品种有重要意义。【前人研究进展】高等植物的花序分限花序和无限花序[2]。自然生长的（野生型）油菜属于无限花序类型,即随着花序轴的伸长,不断产生花芽、形成花蕾。具有无限生长习性的油菜在生产上存在一些缺陷,如：植株较高、易倒伏;结角层不集中,整齐度不好;熟期不一致,无效角果多[3]。实现油菜从无限花序向有限花序转变,克服无限花序生长习性的不足,是选育适合机械化收获的油菜品种的新思路。高等植物开花现象涉及一系列精确的发育进程[4,5,6,7,8,9,10]。如果植株花序分生组织表现持续激活,就会源源不断产生花分生组织,表现为无限花序类型。目前,在一些物种中已发现无限花序植株的有限花序变异。SHANNON等[11]通过EMS诱变产生了有限花序类型的拟南芥;BRADLEY等[12]发现了一个自然变异的金鱼草有限花序突变体;DHANASEKAR等[13]通过γ射线诱导豌豆产生了有限花序突变体;KAUR等[14]在人工合成芥菜型油菜创制过程中产生了有限花序变异类型的芥菜型油菜;ZHANG等[15]通过EMS诱变产生了芝麻有限花序突变体;刘毅[16]发现了自然突变的有限花序芝麻突变体。在不同物种的有限花序突变体中,花序形态基因突变被认为是控制无限花序向有限花序变异的关键因素。例如,金鱼草的基因CENTRORADIALI（CEN）是最早克隆到的导致无限花序向有限花序变异的花序形态基因[12]。在随后的研究中发现,不同物种的CEN直系同源基因也与植株从无限花序向有限花序变异相关,包括拟南芥TERMINAL FLOWER 1（TFL1）[17]、番茄SELF-PRUNING（SP）[18]、烟草CEN-LIKE GENES FROM TOBACCO（CET）[19]、黑麦草LpTFL1[20]和豌豆PsTFL1-a[21],近年来报道的柑橘[22]、水稻[23]、苹果[24]等物种中的CEN/TFL1同源基因均被认为是控制花序变异的关键基因。在一些物种的有限花序突变体中,已经定位了基因所在的关键区间：如在芝麻中,有限花序调控基因定位到连锁群LG8上18.0—19.2 cM的遗传区间[15];芸薹属人工合成芥菜型油菜中的有限花序突变体[25,26],其突变基因定位到B5染色体上;甘蓝型油菜中,通过小孢子培养得到的一个有限花序双单倍体（DH）株系,经遗传分析发现有限花序的遗传模式为隐性单基因调控,其调控基因定位到A10染色体68 kb的物理区段,预测BnTFL1为关键候选基因[27]。【本研究切入点】尽管目前已经在若干无限花序物种中发现了有限花序突变体,并且相关基因已经被定位甚至克隆,但在油菜上关于有限花序突变体及相关遗传机制的研究鲜见报道[27]。【拟解决的关键问题】本研究发现了一个与前人报道不同的自然产生的油菜有限花序突变体（denterminate inflorescence 1,di1）,基于高通量测序平台对di1突变体和野生型以及F2代的2个极端混合池开展全基因组重测序,定位与油菜有限花序性状显著关联的区域,并通过基因克隆测序筛选关键候选基因,为适合机械化收获的油菜新品种培育提供新思路和新材料,同时也为揭示油菜有限花序遗传机制奠定基础。1 材料与方法
1.1 材料
2008年在育种材料中发现了一株油菜有限花序自然突变体di1,连续5年套袋自交获得稳定遗传的突变株系FM8。2015年以突变株系FM8与正常自交系FM7正反交获得F1和F2。从FM8为母本与FM7杂交构建F2分离群体中挑选20个有限花序性状明显的F2单株和20个正常F2单株构建成突变类型混合池和野生类型混合池,用于有限花序性状基因定位。1.2 田间试验和性状调查
经过多年田间观察,确定油菜有限花序性状的评价方法如下：花期花序生长轴顶端发育形成顶生花,花序继续生长和后续花蕾形成均受阻,结角期表现出主枝和侧枝花序顶端没有出现无效角果,有封顶的特征。2016年花期,F1单株套袋,收获自交种子获得F2。2016年10月5日采用穴盘育苗方式播种F2以及2个亲本种子,11月10日移栽到浙江省农业科学院试验大田。田间种植方式为3行区,行距为0.4 m,株距为0.3 m。田间管理按照当地常规方式进行,水平一致。初花期观察每个F2单株花序生长情况,并挂牌记录,终花期再次观察确认F2单株的花序状态。
1.3 极端混合池构建与基因分型
F2群体中挑选20株有限花序性状明显的突变单株和20株野生类型单株,连同2个亲本FM8和FM7,取幼嫩植株叶片0.15 g,利用OMEGA HP Plant DNA提取试剂盒提取DNA。将20个有限花序F2单株和20个野生型F2单株的DNA分别等量混合,构建成突变性状混合池和野生类型混合池。2个混合池DNA和2个亲本DNA按照标准流程构建文库,并通过illumina HiSeqTM PE150对2个混合池和两亲本开展20×和10×覆盖度的全基因组重测序。测序得到的原始序列去除接头和低质量序列,再通过BWA软件[28]（参数：mem-t4-k32-M）比对到油菜参考基因组（http://www.genoscope.cns.fr/brassicanapus/data/Brassica_napus_v4.1.chromosomes.fa.gz）。比对结果经SAMTOOLS的rmdup命令去除重复[29]。采用GATK3.3软件的UnifiedGenotyper模块进行多个样本单核苷酸多态性（SNPs）标记和插入缺失（InDels）标记的检测,使用VariantFiltration进行过滤[30],并利用ANNOVAR软件对SNP和InDel进行注释[31]。1.4 基于混合池的有限花序性状基因定位与候选基因预测
基于基因分型的结果,筛选亲本间纯合差异的多态性位点。选择亲本FM8作为参照亲本,参考TAKAGI等[32]的方法计算2个子代混合池在每个多态性位点上的SNP频率（SNP-index）。为减少测序错误和比对错误造成的影响,对计算出SNP-index后的亲本多态性位点进行过滤,过滤标准如下：（1）2个子代中SNP-index都小于0.3,并且SNP深度都小于7的位点,过滤掉;（2）1个子代SNPindex缺失的位点,过滤掉。随后,计算2个子代SNP-index作差（△SNP-index）。选择1 Mb为窗口、1 kb为步长对△SNP-index在各个染色体上的分布进行作图,选取95%置信水平作为筛选的阈值,置信水平以上的窗口作为候选区间。为了预测关联区间内的候选基因,我们收集了拟南芥中与花序和花发育相关的基因（数据来源于BRAD数据库http://brassicadb.org/brad/）,利用blast2go[33]中的blastall将这些基因的CDS序列与101 040个油菜基因的CDS序列进行比对,以比对覆盖率（>70%）和E值（E<1e-50）为标准寻找对应到油菜上的同源基因。如果花序和花发育相关基因对应到油菜上的同源基因位于显著关联的候选区间,则认为该基因为候选基因。1.5 候选基因分子克隆和序列分析
参照油菜基因组数据库（http://brassicadb.org/brad/）公布的目的基因序列设计引物（电子附表1）。以有限花序和无限花序亲本,以及有限花序和无限花序的F2株系的基因组DNA为模板进行扩增,电泳后回收纯化,连接、转化,PCR检测,挑选出与预测片段分子量大小一致的阳性克隆送到公司测序。Table 1
表1
表1有限花序突变株系与野生型杂交组合后代分离比例
Table 1Segregation of the crosses between determinate inflorescence mutant and wildtype in progeny populations
| 群体 Population | 总株数 Total plants | 有限花序 Determinate inflorescence | 无限花序 Indeterminate inflorescence | 期望比 Expected ration | χ2 | P |
|---|---|---|---|---|---|---|
| FM8 | 40 | 40 | 0 | - | - | - |
| FM7 | 40 | 0 | 40 | - | - | - |
| F1 (FM8×FM7) | 40 | 0 | 40 | - | - | - |
| F1 (FM7×FM8) | 40 | 0 | 40 | - | - | - |
| F2 (FM8×FM7) | 260 | 33 | 227 | 3﹕1 | 12.88 | <0.01 |
| F2 (FM8×FM7) | 260 | 33 | 227 | 13﹕3 | 3.71 | >0.05 |
新窗口打开
2 结果
2.1 油菜有限花序性状的形态观察与鉴定
油菜有限花序突变体di1在初花期表现出主花序和侧枝花序顶部形成1个或者若干个柱头外露的顶生花（图1-A）,此类型顶生花不仅自身无法正常开花,同时也限制了花序生长和后续花蕾形成。结角期,突变体表现出主枝和侧枝有封顶的特征（图1-C）。由于存在营养竞争,无限花序顶端角果结实率往往显著降低,产生无效角果（图1-D）。而有限花序顶部通常不存在无效角果（图1-C）。由于有限花序的形成,突变体表现出株高降低、成熟期更一致并且熟期提早等有利特征。除此之外,有限花序突变体与野生型油菜的其他农艺性状并无明显差别,表现为植株粗壮,叶色较深,分裂较浅,叶面光滑,初花期为3月10日。除顶端花蕾柱头外露,其他部位的花器发育良好,花瓣平滑。
显示原图|下载原图ZIP|生成PPT图1油菜有限花序突变体与野生型的花序形态比较 A：花期的油菜有限花序突变体;B：花期的野生型油菜;C：结角期的油菜有限花序突变体;D：结角期的野生型油菜
-->Fig. 1The inflorescence morphology comparison of rapeseed between determinate growth habit mutant and wildtypeA: Morphology of determinate inflorescence oilseed rape mutant at flower stage; B: Morphology of wildtype oilseed rape at flower stage; C: Morphology of determinate inflorescence oilseed rape mutant at pod development stage; D: Morphology of wildtype oilseed rape at pod development stage
-->
2.2 油菜有限花序性状的遗传分析
以稳定遗传的油菜有限花序株系FM8作为亲本,与野生型油菜自交系FM7正反交,F1均表现为野生型无限花序,说明该有限花序突变为隐性突变。突变株系FM8与正常自交系FM7构建的F2群体中出现了有限花序与无限花序的分离,卡方测验表明分离比符合13﹕3分离（表1）,说明油菜有限花序性状受2对隐性基因和1对隐性上位抑制基因互作控制。2.3 极端群体混合池构建与测序数据分析
对20个有限花序F2单株和20个无限花序F2单株构建的混合池以及两亲本开展全基因组重测序,总共得到53.17 G的原始数据。过滤后4个样本的有效序列数据量在8 722.79—17 717.90 M,总数据量为53.026 G。测序数据Q20>96.9%,Q30>95.29%,GC含量在37.35%—38.73%,测序数据有95.36%—98.56%的序列可以成功比对到参考基因组。由此可知,所有样本的数据量足够,测序质量合格,GC分布正常,测序数据与油菜参考基因组比对结果正常,可用于后续的变异检测及性状的基因定位。将突变体亲本FM8和野生型亲本FM7测序结果比较,共获得1 186 750个SNPs和764 900个InDels。亲本间检测到的SNPs明显多于InDels。根据变异发生的位置进行比较发现,无论SNPs变异还是InDels变异,位于基因非编码区的多态性位点显著多于基因编码区。基于基因分型的结果,筛选2个亲本间纯合差异的标记,最终挑选出137 759个多态性位点,包括30 123个SNPs和107 636个InDels,用于后续基因定位。这些多态性位点在染色体间呈现非均匀分布,其中A07染色体拥有最丰富的多态性位点,而A04染色体拥有最少的多态性位点。
2.4 油菜有限花序性状基因定位
分析2个混合池△SNP-index,用于甘蓝型油菜有限花序性状的基因定位。定位结果显示,分布于甘蓝型油菜5条染色体的7个区域出现了超过临界值水平的极显著峰,说明这些区域可能包含调控油菜有限花序变异的基因（图2）。这7个显著关联区间在油菜基因组上的分布为Chr.A08：12.36—12.51 Mb、Chr.A09：26.01—26.05 Mb、Chr.A09：30.34—30.68 Mb、Chr.A10：14.43—17.37 Mb、Chr.C08：27.38—27.75 Mb、Chr.C09：45.80—46.45 Mb和Chr.C09：46.54—46.60 Mb。这些显著关联区间覆盖的染色体长度为0.04—3.04 Mb。其中,染色体A10关联区间顶点峰值最高,暗示该区间可能具有引起油菜有限花序变异的主效基因。
显示原图|下载原图ZIP|生成PPT图2与油菜有限花序性状显著关联的7个基因组区间及其包含的候选基因红色曲线代表利用滑动窗口计算得到的整个染色体上ΔSNP-index的变化趋势。绿色虚线代表显著性阈值,如果一个区域的滑动窗口曲线出现明显的峰谷并超过临界值,则认为这个区域与性状显著关联,每个点应对一个SNP,横坐标代表每个SNP在油菜基因组上的物理位置,纵坐标代表ΔSNP-index 值
-->Fig. 2The seven significantly associated genomic regions for the trait of determinate inflorescence in rapeseed and the containing candidate genesThe figure showed the graph of ΔSNP-index spots, red curve was the result of sliding window analysis which was applied to visualize the tendence of ΔSNP-index plots, green curve was threshold value. The significantly associated regions can be identified as peaks or valleys of the ΔSNP-index plot, each spot corresponded to a ΔSNP-index, the x-axis corresponded to the chromosomal position
-->
2.5 染色体A10主效关联区域与C09关联区域的微共线性比较研究
从表2可以看出,A10染色体关联区间对应到拟南芥基因组上的共线性区段（Chr.5：0—4.02 Mb）完全覆盖并延伸了C09染色体2个关联区间对应到拟南芥上的共线性区段（Chr.5：3.10—3.57 Mb和Chr.5：2.98—3.04 Mb）。进一步分析A10上的关联区间发现,该区间内的14.63—15.07 Mb区段与C09的2个关联区间对应到拟南芥基因组的相同位置（图3）,这说明A10染色体14.36—15.07 Mb的关联区域与C09染色体的2个候选区域具有高度同源性。油菜A和C亚基因组上的同源区域同时被检测出与油菜有限花序性状显著连锁不仅证实了这三个候选区间的可信性,并且可以推测这三个区域对油菜有限花序的调控可能由同源遗传因子决定。
显示原图|下载原图ZIP|生成PPT图3与油菜有限花序关联的A10候选区间部分区段与和C09候选区间的微共线性比较左右两边绿色线条代表A10关联区间部分区段和C09关联区间所所对应的物理位置,中间红色和绿色线条为两者对应到拟南芥基因组上的物理位置,其中红色部分代表A10关联区间部分区段和C09关联区间对应到拟南芥基因组相同位置
-->Fig. 3Synteny blocks of candidate regions for the trait of determinate inflorescence in rapeseed between A10 and C09The green bars on the left and right sides represented the physical location of the partial associated region on A10 and associated regions on C09 in rapeseed, respectively. The bars filled with green and red colors represented the physical location on A. thaliana of the beside rapeseed regions, where red regions were the overlapped regions of the beside rapeseed regions
-->
Table 2
表2
表2有限花序性状7个关联区间对应到拟南芥基因组上的共线性区段
Table 2Synteny blocks of seven associated regions of rapeseed determinate growth habit on genomes of A. thaliana
| 关联区间 Associated regions (Mb) | 拟南芥共线性区段 Synteny blocks on genomes of A. thaliana (Mb) |
|---|---|
| Chr.A08: 12.36-12.51 | Chr.4: 12.52-12.60; Chr4: 16.91-17.03 |
| Chr.A09: 26.01-26.05 | Chr.3: 20.87-20.91 |
| Chr.A09: 30.34-30.68 | Chr.2: 8.44-8.47; Chr1: 6.24-6.33; Chr1: 6.49-6.74 |
| Chr.A10: 14.43-17.37 | Chr.5: 0-4.02 |
| Chr.C08: 27.38-27.75 | Chr.3: 20.30-20.57 |
| Chr.C09: 45.80-46.45 | Chr.5: 3.10-3.57 |
| Chr.C09: 46.54-46.60 | Chr.5: 2.98-3.04 |
新窗口打开
2.6 油菜有限花序性状候选基因预测与克隆
通过数据库收集拟南芥与花序和花发育相关的基因（数据来源于BRAD数据库http://brassicadb.org/ brad/）,并比对得到这些基因对应到油菜基因组上的同源基因（电子附表2）。8个基因被鉴定到位于A08、A09、A10、C08和C09的5个关联区间（表3）：GSBRNA2T00016119001（FLOWERING LOCUS C,FLC）、GSBRNA2T00016124001（FLC）、GSBRNA2T00027297001（ATBZIP14,FD）、GSBRNA2T00135921001（FLC）、GSBRNA2T00136426001（TFL1）、GSBRNA2T00016144001（MULTICOPY SUPPRESSOR OF IRA1 4,FVE）、GSBRNA2T00016146001（FVE）和GSBRNA2T00108914001（SCHLAFMUTZE,SMZ）。其中,TFL1是花序形态关键基因,被锚定到具有最高峰值的A10染色体主效关联区间,因此,推测TFL1为调控有限花序性状的关键基因。Table 3
表3
表3关联区间内的候选基因预测
Table 3The predicted candidate genes in the mapped associated regions
| 候选基因 Candidate genes | 关联区间 Associated regions | |||||||
|---|---|---|---|---|---|---|---|---|
| 油菜基因编号 B. napus ID | 名称 Name | 拟南芥基因编号 A. thaliana ID | 染色体 Chr. | 起点Start (Mb) | 终点End (Mb) | 染色体 Chr. | 起点Start (Mb) | 终点End (Mb) |
| GSBRNA2T00016119001 | FLC | AT5G10140 | Chr.C09 | 46.37 | 46.37 | Chr.C09 | 46.54 | 46.6 |
| GSBRNA2T00016124001 | FLC | AT5G10140 | Chr.C09 | 46.35 | 46.35 | Chr.C09 | 46.54 | 46.6 |
| GSBRNA2T00027297001 | FD | AT4G35900 | Chr.A08 | 12.45 | 12.45 | Chr.A08 | 12.35 | 12.51 |
| GSBRNA2T00135921001 | FLC | AT5G10140 | Chr.A10 | 15 | 15 | Chr.A10 | 14.33 | 17.37 |
| GSBRNA2T00136426001 | TFL1 | AT5G03840 | Chr.A10 | 16.77 | 16.77 | Chr.A10 | 14.33 | 17.37 |
| GSBRNA2T00016144001 | FVE | AT2G19520 | Chr.A09 | 30.34 | 30.34 | Chr.A09 | 30.34 | 30.68 |
| GSBRNA2T00016146001 | FVE | AT2G19520 | Chr.A09 | 30.35 | 30.35 | Chr.A09 | 30.34 | 30.68 |
| GSBRNA2T00108914001 | SMZ | AT3G54990 | Chr.C08 | 27.41 | 27.42 | Chr.C08 | 27.38 | 27.75 |
新窗口打开
利用有限花序亲本FM8和有限花序F2株系,野生型亲本FM7和野生型F2株系的基因组DNA对有限花序性状关联区间内的8个潜在候选基因进行基因克隆和测序。结果显示,有3个潜在候选基因发生了编码区的碱基序列变异,并最终导致氨基酸序列变异：位于A10染色体上的主效关联区间内的TFL1,有限花序亲本在其编码区存在2个SNPs变异,并最终导致TFL1蛋白序列发生了2个氨基酸的变异（图4）;位于A09染色体关联区间内的FVE,有限花序亲本在其编码区存在6个碱基的插入和22个SNPs变异,并最终导致FVE蛋白序列发生了2个氨基酸的插入和8个氨基酸序列变异（图4）;位于C08染色体关联区间内的SMZ,有限花序亲本在其编码区发生了31个碱基片段的丢失,导致第146—174位置的氨基酸序列发生变异,并且第175位置发生氨基酸终止编码（图4）。
显示原图|下载原图ZIP|生成PPT图4有限花序突变体与野生型候选基因蛋白质序列的比较箭头和虚线区域代表氨基酸序列发生变异的位置
-->Fig. 4Sequence polymorphisms of the protein sequence of candidate genes between lines with determinate and indeterminate growth habit The arrows and dotted area represented protein polymorphisms of candidate genes
-->
3 讨论
3.1 di1突变体是一个新型的油菜有限花序突变体
在芸薹属作物中已有关于有限花序变异的报道。KAUR等[26]在远缘杂交获得的人工合成芥菜型油菜后代中发现了受隐性单基因控制的芥菜型油菜有限花序突变体。LI等[27]在甘蓝型油菜不同生态型杂交的小孢子培养后代中也发现了受隐性单基因控制的甘蓝型油菜有限花序突变体。与前人报道的芸薹属中的有限花序突变体进行比较发现：（1）从遗传模式上看,本研究中的油菜有限花序突变受2对隐性基因和1对隐性上位抑制基因互作控制。这个遗传模式与前人报道的芸薹属中的有限花序突变体并不完全相同。（2）从突变体产生方式看,上述得到的芥菜型油菜和甘蓝型油菜突变体为远缘杂交或者不同生态型间杂交诱导产生,而本研究中的油菜有限花序突变体为自然突变产生。（3）从基因定位结果看,芥菜型油菜有限花序调控基因定位在B5染色体1.9 cM的区段,但其调控基因尚未报道。LI等[27]报道的甘蓝型油菜有限花序调控基因定位于A10染色体,预测BnTFL1为候选基因,该基因与本研究预测得到的位于主效区域的候选基因相同。然而,除此之外,本研究还定位得到其他调控位点,这与该性状受2对隐性基因和1对隐性上位抑制基因互作控制的遗传模式吻合。（4）尽管BnTFL1在本研究和LI等[27]的报道中均为关键候选基因,但变异位点存在差异。LI等[27]报道中BnTFL1蛋白序列的第46和第47位的氨基酸产生变异,而本研究BnTFL1蛋白变异位点位于第37位和第139位。综上所述,本研究得到的油菜有限花序突变体与前人报道的有限花序突变体并不相同,该新型油菜有限花序突变体的发现对于拓宽甘蓝型油菜的变异类型具有一定意义。针对预测得到的候选基因开展更加深入的功能验证将在后续研究中继续开展。3.2 di1突变体有限花序候选基因的预测
利用本研究定位到的7个关联区间总共预测得到8个潜在候选基因,包括1个TFL1、3个FLC、2个FVE、1个FD和1个SMZ。TFL1、FLC、FVE、FD和SMZ在拟南芥花序和花形态建成的过程中具有非常重要的作用[17,34-41]（图5）。基因克隆揭示有限花序突变体中的TFL1、FVE和SMZ发生了基因编码区的碱基序列变异,并最终导致氨基酸序列变异。根据前人研究报道,TFL1属于磷脂酰乙醇胺结合蛋白家族,是维持花序分生组织特性所必需的,有抑制开花的作用[17],因此,TFL1在植物从营养生长到生殖生长这一成花转化过程中扮演非常重要的角色[42,43,44]。在拟南芥中,AtTFL1的突变将导致花序分生组织转化为花分生组织,形成单个的顶生花,最终使得花序由无限花序变为有限花序[11,17,45]。在豌豆和豇豆中TFL1的作用均被揭示为维持花序分生组织的无限生长习性,该基因的突变将产生有限花序类型的豌豆和豇豆[13,21]。此外,金鱼草[12]、番茄[18]、烟草[19]、黑麦草[20]、柑橘[22]、水稻[23]、苹果[24]等多个物种中均报道了由TFL1同源基因变异导致的有限花序突变。本研究中,TFL1位于A10染色体主效关联区域内,并且该基因在野生型和突变体中存在编码区和蛋白序列的变异,因此,推断TFL1对控制油菜有限花序变异具有关键作用。FVE是自主开花途径中的一个关键基因,编码了一个MSI1同源基因,跟组蛋白的脱乙酞化复合体的形成有关。它的存在对FLC染色体组蛋白去乙酞化是必要的,FLC染色体组蛋白去乙酞化后,FLC由活化状态转变为非活化状态,从而启动植物开花[46]。刘新庆等[47]将拟南芥fve突变体与野生型相比较发现在野生型植株正常开花的阶段突变体植株尚未发育出花序,因此认为fve突变体影响了拟南芥花序的发育,并最终导致晚花表型。SMZ位于光周期途径,编码AP2/EREBP转录因子,但其作用却独立于CONSTANS,是通过与miR172的结合来抑制FT表达,最终抑制成花[48]。尽管这些基因位于不同的途径,但均可以作用于下游的花序分生组织特异基因,最终参与到花序和花的形态建成（图5）。
显示原图|下载原图ZIP|生成PPT图5位于关联区间内的候选基因所涉及的花序和花发育调控网络箭头表示促进作用,线段表示抑制作用;基因右上角星号代表候选基因
-->Fig. 5The involved regulation network of inflorescence and floral development of candidate genes in associated regions The arrows and segments represented positive and negative regulation, respectively; the stars in the upper right corner of genes represented candidate genes
-->
4 结论
油菜有限花序性状受2对隐性基因和1对隐性上位抑制基因互作控制。有7个区间与有限花序性状显著关联。TFL1、FVE和SMZ被推断为有限花序性状的候选基因。（责任编辑 李莉）
The authors have declared that no competing interests exist.
