水稻明恢826光身基因GL826的精细定位与候选基因分析

核农学报2024，38（3）：0413~0423
Journal of Nuclear Agricultural Sciences
水稻明恢826光身基因GL826的精细定位与
候选基因分析
韦新宇1， 2张文静1， 2曾跃辉1， 2肖长春1， 2张锐1， 2
黄建鸿1， 2许旭明2， 3， *
（1三明市农业科学研究院生物技术研究所，福建三明365500；
2福建省（山区）作物遗传改良与创新利用重点实验室，福建三明365500；3三明市农业科学研究院水稻研究所，福建三明365500）摘要：光身稻是一种重要的种质资源，在杂交水稻育种中具有重要的利用价值，其中明恢826（MH826）是自主选育的优良光身稻恢复系。

为探究MH826光身基因功能和分子机制，本研究对MH826光身表型进行鉴定、对光身基因进行精细定位和候选基因分析。

扫描电镜分析结果发现，其叶片和籽粒颖壳表现为表面光滑，未观察到明显的尖刺状表皮毛特征。

利用MH826分别与具有明显表皮毛特性的水稻品种华占、明恢2155和9311杂交构建的F1和F2代分离群体对MH826光身基因进行遗传学分析，结果表明，MH826的光身性状由单个细胞核隐性基因控制，并将其命名为GL826（Glabrous Leaf 826）。

进一步应用MH826与华占杂交获得的F2分离群体，结合分离群体分析法（BSA）和隐性群体分析法（RCA），并通过Mapmaker/ Exp3.0和MapDraw v2.1软件分析，将光身基因GL826精细定位在水稻第5号染色体短臂端，且位于分子标记InDel-70和InDel-11之间，相对遗传距离均为0.2 cM。

区间物理距离为141.03 kb，共预测有23个开放阅读框，其中LOC_Os05g02730功能注释为已报道的和水稻表皮毛发育相关基因。

测序分析结果发现，MH826中候选基因LOC_Os05g02730在编码区和启动子区均未存在DNA水平上的碱基变异。

实时荧光定量PCR（qRT-PCR）分析结果表明，LOC_Os05g02730在MH826不同生长发育阶段的叶片、茎和幼穗中的表达受到了严重抑制，表现为不同程度的下调。

推测候选基因LOC_Os05g02730即为GL826。

本研究结果为光身稻恢复系MH826及其光身基因GL826在杂交水稻育种上的进一步应用提供了理论基础。

关键词：水稻；明恢826；光身基因；精细定位；候选基因分析
DOI：10.11869/j.issn.1000‑8551.2024.03.0413
光身稻（Oryza glaberrima）因其茎秆、叶片和谷壳等组织表面光滑无毛而得名，又叫光叶稻或光壳稻。

光身稻具有广亲和性、对强光适应能力强、抗倒伏性强、适合机械化操作和节省储运空间等优良特性［1-4］，已作为杂交水稻育种重要的种质资源被加以研究和利用［5］。

植物表皮毛起源于表皮细胞，是一种表皮组织向外延伸的毛状突起，在植物体中普遍存在，可由单细胞或多细胞构成，形态结构多样［6-8］。

水稻的表皮毛为双细胞毛状体，呈圆柱状分布［9］。

目前，对水稻表皮毛性状的研究主要集中在叶片和颖壳［5］。

光身稻的无表皮毛性状主要表现为叶片和颖壳表面毛状体的减少或缺失。

已报道的与水稻茸毛生长发育相关的基因
HL6［10］、HL1［11］和GL6［12］，均已被定位在水稻第6号染色体上，且叶片长毛表型均受细胞核单基因控制，属显性遗传。

已报道的水稻光身基因wmg2［13］位于第
3号染色体上，wmg1［13］位于第4号染色体上，GL5［10］、NUDA/GL-1［14］、dep［9］、gl-l［15］、GLAG5［16］位于第5号染色体上，GLR3［17］、OsSPL10［18］、gl6［19］、GLL［20］位于第6号染色体上。

其中，GL5基因被精细定位在水稻第5号染色体33 kb区间内，GL5基因的表达水平降低导致表皮毛的缺失［10］；NUDA/GL-1基因座被精细定位在水稻第5号染色体28.5 kb区间内，通过RNAi和互补分析，确认一个编码具有WUSCHEL结构域蛋白的同源基因（OsWOX3B）作为候选基因［14］；dep基因被精细定位在
文章编号：1000‑8551（2024）03‑0413‑11
收稿日期：2023‑10‑09 接受日期：2023‑11‑08
基金项目：国家现代农业产业技术体系专项（CARS-01），福建省自然科学基金项目（2021J01535），三明市科技计划项目（2022-N-7）
作者简介：韦新宇，男，副研究员，主要从事水稻遗传育种研究。

E-mail：wxy1209@
*通讯作者：许旭明，男，研究员，主要从事水稻遗传育种研究。

E-mail：fj63xxm@
413
核农学报38 卷
水稻第5号染色体22 kb区间内，互补试验证实dep基因对无毛表型具有调节作用［9］。

位于第6号染色体上的OsSPL10启动子区发生变异导致籼稻品种WD-17993出现无毛性状［18］，也有研究表明OsSPL10是水稻茸毛发育的正向调控基因，其与OsWOX3B、OsPLT2基因间存在遗传互作，共同调控水稻表皮毛的起始发生和伸长［21］。

明恢826（MH826）是三明市农业科学研究院利用亲本IR249和NO3S2杂交配组，经多世代自交选择、抗病鉴定以及与大量不育系测配筛选育成的优良光身稻恢复系。

该恢复系具有株型好、结实率高、米质优、配合力强等优点，所配杂交稻组合已陆续进入生产应用，其育种利用和研究价值较好。

但目前对MH826的光身性状及控制光身性状的基因尚不了解。

因此，本研究对
MH826的光身性状进行显微观察和遗传学分析，利用分离群体分析法（bulked segregation analysis，BSA）和隐性群体分析法（recessive class analysis，RCA），对其光身基因进行精细定位和候选基因分析，以期为进一步明确其基因功能和分子机制奠定基础，同时为光身稻品种及其光身基因在水稻育种上的应用提供理论依据。

1　材料与方法
1.1　试验材料
利用光身稻恢复系MH826与具有明显表皮毛特性的水稻品种华占、明恢2155和9311杂交，获得F1代。

F1代自交后获得F2代分离群体，用于光身基因的遗传学分析和精细定位。

2020年晚季在福建沙县三明市农业科学研究院水稻试验基地种植亲本、F1和F2代群体，常规田间管理。

1.2　表型鉴定与遗传分析
分别在苗期、分蘖期和抽穗期对亲本、F1和F2代分离群体单株进行表型鉴定和统计学分析。

以手指自叶尖沿叶面和边缘向叶基部触摸无粗糙感作为光身表型鉴定依据。

通过统计F2代分离群体中光身单株和有毛单株的分离比进行遗传学分析，并利用χ2进行适合性检验。

1.3　扫描电镜观察
选取水稻分蘖期的叶片和灌浆期的种子用2.5%戊二醛固定液于4 ℃冰箱内固定1 h以上，然后用0.2 mol·L-1的磷酸缓冲液（pH值7.2）漂洗3次，每次20 min，再用1%锇酸固定1 h，然后用0.2 mol·L-1的磷酸缓冲液（pH值7.2）漂洗3次，每次20 min，再用梯度乙醇溶液进行脱水，每次20 min。

之后在K850 CO2临界点干燥仪（Quorum，英国）内进行干燥处理，再将样品粘在样品台上用MC1000离子溅射仪（HITACHI，日本）进行离子溅射镀膜处理，最后利用SU8100扫描电子显微镜（HITACHI，日本）观察并拍照［22］。

1.4　基因的精细定位
采用改良的十六烷基三甲基溴化铵（cetyltrimethyl ammonium bromide，CTAB）法［23］提取水稻苗期叶片全基因组DNA。

在F2分离群体中分别随机选择15个光身单株和15个有毛单株提取基因组DNA，浓度统一稀释为50 ng·μL-1，等量混合从而构建光身基因池和有毛基因池［24］。

通过水稻生物信息学网站（http：//www. /）搜索PCR扩增产物大小在300 bp以内，且均匀分布于水稻12条染色体上的356对简单重复序列（simple sequence repeat，SSR）引物对亲本和基因池进行多态性检测并对F2隐性单株进行基因型分析，筛选与目标基因紧密连锁的分子标记。

PCR扩增反应体系为10 μL：2×Taq PCR MasterMix（艾德莱，北京）5.0 μL，50 ng·μL-1 DNA模板0.8 μL，10 μmol·L-1正、反向引物各0.5 μL，ddH2O 3.2 μL。

PCR扩增反应程序：94 ℃预变性3 min；94 ℃变性10 s，55～60 ℃退火10 s，72 ℃延伸10 s，35个循环；72 ℃终延伸5 min；4 ℃保存。

PCR扩增反应结束后，将PCR扩增产物用8.0%的非变性聚丙烯酰胺凝胶进行电泳，通过GelDoc XR+凝胶成像系统（Bio-Rad，美国）进行拍照观察。

另外，根据粳稻品种日本晴和籼稻品种9311之间的基因组序列差异，利用Primer Premier 5.0软件进一步设计和开发分布密度更高的插入缺失（insertion-deletion，InDel）标记，用于光身基因GL826的精细定位。

PCR反应及电泳程序方法和步骤同上。

利用Mapmaker/Exp3.0软件进行连锁分析，采用Kosambi运算函数将重组率转换为遗传距离，结合MapDraw v2.1软件绘制遗传连锁图谱。

根据水稻生物信息学网站（http：//www.gramene. org/）提供的水稻基因组序列信息并通过染色体步移的方法构建物理图谱。

1.5　候选基因分析
通过水稻基因组注释数据库（Rice Genome Annotation Project，RGAP，http：///）对精细定位区间内候选基因的功能进行预测，并结合实时荧光定量PCR（quantitative real-time PCR，qRT-PCR）以及克隆和测序技术对候选基因进行表达和序列分析，具体方法参照文献［25］。

2　结果与分析
2.1　MH826表型鉴定和扫描电镜分析
本研究以手指触摸的方式对研究材料的表型进行
414
3 期水稻明恢826光身基因GL826的精细定位与候选基因分析初步鉴定。

结果表明，具有表皮毛特性的水稻品种华占、明恢2155和9311的叶片和籽粒颖壳具有明显的粗糙感，手指在叶片和颖壳表面滑动时存在明显的阻滞感；而光身稻恢复系MH826的叶片和籽粒颖壳触摸光滑，无粗糙感。

进一步采用扫描电镜分析了具有表皮毛特性的水稻品种华占和光身稻恢复系MH826分别在分蘖期叶片以及灌浆期籽粒的表皮毛发育情况。

结
果表明，华占和MH826在叶片和籽粒表皮存在显著的表皮毛差异。

在华占的叶片和籽粒表皮可观察到明显的尖刺状表皮毛，包括大刺毛和小刺毛，且在叶片表面分布相对密集（图1-A~D ）。

而在MH826的叶片和籽粒表皮仅存在少量的腺毛，未观察到明显的尖刺状表皮毛，呈现表面光滑特性（图1-E~H ），其与手指触感鉴定结果表现一致。

2.2　MH826光身性状的遗传学分析
为了对MH826的光身性状进行遗传特性分析，本研究利用MH826分别与具有表皮毛特性的水稻品种华占、明恢2155和9311配制杂交组合，获得F 1代并构建了3个独立的F 2代分离群体。

在分蘖盛期，分别对华占×MH826，明恢2155×MH826和MH826×9311获得的F 1代植株以及F 2代群体进行叶片表型鉴定和统计学分析。

结果表明，杂交获得的所有F 1代
植株叶片均表现为明显的表皮毛性状，与具有表皮毛特性的亲本叶片表型表现一致，说明MH826的光身性状属于隐性遗传。

另外，在3个独立的杂交F 2代群体中，显性单株（有毛单株）和隐性单株（光身单株）表现出明显的分离，且均符合孟德尔单基因控制的3∶1分离比（χ2<χ20.05）（表1）。

以上结果表明，MH826的光身性状由单个细胞核隐性基因控制。

本研究将其命
名为GL826 （Glabrous Leaf 826）。

注：A~C ：分蘖期华占叶片上表皮，B 、C 分别为A 、B 红色矩形区域放大视野；D ：灌浆期华占籽粒表皮表型；E~G ：分蘖期MH826叶片上表皮，F 、G 分别
为E 、F 红色矩形区域放大视野；H ：灌浆期MH826籽粒表皮表型。

Ma ：大刺毛；Mi ：小刺毛；Gla ：腺毛；Sto ：气孔。

Note ： A -C ： Morphology of upper epidermis of leaves in Huazhan at the tilling stage ， B and C show magnified map of the area of red rectangle in A and B ，
respectively. D ： Morphology of grain epidermis of Huazhan at grain filling stage. E -G ： Morphology of upper epidermis of leaves in MH826 at the tilling stage ，
F and
G show magnified map of the area of red rectangle in E and F ， respectively.
H ： Morphology of grain epidermis of MH826 at grain filling stage.
Ma ： Macro hairs. Mi ： Micro hairs. Gla ： Glandular hairs. Sto ： Stomata.
图1 华占和MH826的叶片、籽粒扫描电镜观察
Fig.1　Observation of leaf and grain phenotypes of two parental lines （Huazhan and MH826） by scanning electron microscope
表1　MH826杂交组合的分离分析Table 1　Segregation analysis of MH826 crosses
杂交组合
Cross combinations 华占×MH826Huazhan×MH826明恢2155×MH826
Minghui2155×MH826
MH826×9311
F 2表型 F 2 Phenotype
分离群体
Segregation population
1 785590844
显性单株Dominant plants
1 356432621
隐性单株Recessive plants
429158224
显性∶隐性Dominant ∶Recessive
452∶143216∶79155∶56
χ20.891.101.07
χ20.053.843.843.54
415
核农学报38 卷
2.3　GL826基因的精细定位
由表1可知，华占×MH826、明恢2155×MH826和
MH826×9311构建的3个F 2代分离群体植株数量分别是1 785、590和844个。

综合考虑基因精细定位对大群体数量的要求，本研究优先选用群体植株数量较大的华占×MH826获得的429个F 2隐性单株对光身基因GL826进行精细定位。

首先利用平均分布于水稻12条染色体上的356对SSR 引物对亲本以及两个极端基因池进行多态性筛选，通过对F 2隐性单株进行基因型鉴定，并结合Mapmaker/Exp 3.0进行连锁分析。

获得一个位于水稻第5号染色体上的SSR 标记Rm13与GL826基因连锁。

进一步在Rm13附近筛选分布密度更高的SSR 标记对两亲本、基因池以及F 2隐性单株进行多态性检测和基因型分析，获得与GL826连锁的SSR 标记共7对：Rm17731、Rm6300、Rm17787、Rm17812、Rm17819、Rm5579和Rm405。

并将GL826初步定位在Rm17787和Rm17812之间，相对遗传距离分别为
1.6和
2.4 cM （图2-A ）。

通过比对粳稻品种日本晴和
籼稻品种9311基因组序列，在Rm17787和Rm17812标记之间开发获得插入缺失InDel 标记95对，进一步对亲本、基因池以及F 2隐性单株进行多态性检测和基因型分析，获得与GL826紧密连锁的特异性InDel 标记19对：InDel-02、InDel-11、InDel-13、InDel-14、InDel-17、InDel-24、InDel-26、InDel-28、InDel-31、InDel-40、InDel-44、InDel-48、InDel-61、InDel-64、InDel-68、InDel-70、InDel-71、InDel-72和InDel-76。

最终将光身基因GL826精细定位在分子标记InDel-70和InDel-11之间，且相对遗传距离均为0.2 cM （图2-A ）。

结合
RGAP 数据库，该精细定位区间对应粳稻品种日本晴的物理距离为141.03 kb （图2-B ）。

参考日本晴各细菌人工染色体（bacterial artificial chromosome ，BAC ）克隆序列，并通过染色体步移的方法将每个与光身基因GL826连锁的分子标记在染色体上的位置锚定在相应的BAC 克隆上，从而构建了覆盖GL826基因位点的物理图谱（图2-B ）。

本研究用于GL826精细定位的InDel 标记序列及扩增片段大小如表2
所示。

注：A ：光身基因GL826精细遗传连锁图谱；B ：光身基因GL826物理图谱。

括号中的数字表示相应连锁标记在429个F 2隐性单株中发生交换重组的单
株数；垂直虚线对应的是相应连锁标记在水稻第5号染色体上的实际物理位置；n 表示定位群体的植株数量。

Note ： A ： A fine genetic linkage map of the glabrous leafgene GL826 on chromosome 1. B ： A physical contig map of the region spanning the GL826 locus was constructed using 429 F 2 recessive individuals. The numbers in parentheses below each of the molecular markers linked to the GL826 represent the number of
recombinants among the 429 recessive F 2 individuals. The vertical and dashed lines designate the actual physical location of the corresponding linkage
molecular markers on chromosome 5.n is the number of mapping populations.
图2 光身基因GL826精细定位
Fig.2　Fine mapping of glabrous leaf gene GL826
416
3 期水稻明恢826光身基因GL826的精细定位与候选基因分析
2.4　GL826候选基因分析
参考RGAP水稻基因组注释数据库，发现在该精细定位区间共预测有23个开放阅读框（ORFs）（表3）。

其中候选基因LOC_Os05g02730为已报道的与水稻表皮毛发育相关的基因OsWOX3B，该基因编码一个
WOX3B转录因子家族蛋白，参与植物特异性横向器官的发育调节，OsWOX3B基因突变将导致水稻叶片和籽粒颖壳表皮毛缺失，表现为植株整体光滑无毛的表型［9，14］。

因此，本研究将LOC_Os05g02730作为GL826的重点候选基因并进行进一步的序列和表达分析。

为了鉴定候选基因LOC_Os05g02730是否为GL826，本研究通过PCR扩增和测序技术首先克隆和分析了光身稻恢复系MH826和具有表皮毛特性的水稻品种华占、明恢2155和9311中LOC_Os05g02730基因的编码序列。

结果表明，华占、明恢2155、9311和MH826在LOC_Os05g02730基因的编码区未检测到碱基差异，且序列和已报道的OsWOX3B［9，14］基因序列完全一致（图3-A）。

即光身稻恢复系MH826中含有LOC_ Os05g02730（OsWOX3B）基因完整的编码序列（coding sequence，CDS）。

另外，进一步选择LOC_Os05g02730基因起始密码子上游1.8 kb启动子区域进行PCR扩
增和序列比对分析。

结果表明，华占、明恢2155、9311和MH826在该启动子区域同样未检测到碱基差异并且和OsWOX3B基因启动子序列表现一致（图3-B）。

因此，LOC_Os05g02730基因在具有表皮毛特性的水稻品种华占、明恢2155、9311和光身稻恢复系MH826中未鉴定到DNA水平差异，包括编码区和基因上游的启动子区域。

本研究用于LOC_Os05g02730基因CDS和启动子序列测序的特异性引物如表4所示。

为了进一步确定候选基因LOC_Os05g02730即为光身基因GL826，本研究利用qRT-PCR技术分析了LOC_Os05g02730基因在光身稻恢复系MH826和具有表皮毛特性的水稻品种华占在不同时期以及不同组织部位中的表达差异，其中分别包括LOC_Os05g02730在苗期、分蘖期和抽穗期的根、茎、叶以及抽穗期幼穗组织中的表达。

结果表明，LOC_Os05g02730基因在MH826苗期、分蘖期和抽穗期的茎、叶片以及抽穗期的幼穗中的表达水平均极显著低于具有表皮毛特性的水稻品种华占，尤其在叶片中的表达下调最为明显（图4）。

说明LOC_Os05g02730基因在光身稻恢复系MH826中的表达受到了严重的抑制，其表达量在不同生长发育阶段不同组织中表现为不同程度的下调。

因
表2　用于GL826基因精细定位的InDel标记
Table 2　New developed InDel markers used for fine mapping of the GL826 gene
标记名称Marker name
InDel-02 InDel-04 InDel-11 InDel-13 InDel-17 InDel-24 InDel-26 InDel-28 InDel-31 InDel-40 InDel-44 InDel-48 InDel-61 InDel-64 InDel-68 InDel-70 InDel-71 InDel-72 InDel-76
引物序列
Primer sequence
上游引物（5′→3′）
Forward primer （5′→3′）
GATCGGTTGGAGTTGCTGAT TCCCACAACCTCTATTTTACCGA
GCTTTGCCAACCAGTGAACA
CACTGATTGAACACGCCGT
GGTACACGTGATGGCTTGAG TCCATATAGTCCTGAGAAGCAATC CCTAATGTCACCAGATCAAAGCA
ACCGAAACTGCCATCTCTGA TGGAGGGAGTATTTCGAGAAGA
AGCACGTCCTGGGTAGAATG
TTGCGCAACAGGTGAATACT
CGCACGCCATTACAAACTGA CGGAATTGTGAACCTTGGTGT
ACGCGTCATCATTGTCATGTT
GGGAGGGAGGCCATATGATA
CGCGTAACCCTCCTTACACA
TCCGTAACCTTGACACAAGT
ATAACTCGTTACGGGCTTGC
CCACGTATCAAAAGCCTCCG
下游引物（5′→3′）
Reverse primer （5′→3′）
ACAACCAGTGCAAGTCTCCT
GCAAGCCCGTAACGAGTTAT
TGGCTGCTAACCTCGTTTCT
AGCTTTCCCACGGTCATG
TAGACCGTTTTCCACCCGTT
CAGTGCACACCGTCAATGTC
TGAAGAACCAGTCGTAGCGA
AGGGTTGCTGGTAATTGGGA
GGTGCCAGATATCATACACACA
AACAATTCCAAAGAGCGCGG
TGCCGGATTCCTAATTGCAG
TGTTGTTAAGGTACCCATGTCA
CAAACACAGGGAAGCGTTCA
TGATGTTCGTGGCCTTTTGT
TTCCGCCGGTAAACTGTGA
AGGGCAAAGAAAGGGTCAGA
AGGCCTTCTAGAACCTGGAT
AGAACATGGTCCAGCTATCACT
GCAGCATCCTTTTCCTCTCG
预测PCR产物大小
Expected PCR product size/bp
日本晴
Nipponbare
163
161
257
213
143
223
278
127
134
296
117
130
285
276
269
159
278
141
169
9311
146
174
266
194
134
213
296
117
121
305
128
122
279
270
289
153
292
128
184
417
核农学报38 卷
此，推测LOC_Os05g02730即为GL826基因。

本研究用于LOC_Os05g02730基因表达分析的特异性引物如表4所示。

3　讨论
MH826是三明市农业科学研究院经多世代自交选择、抗病鉴定以及大量不育系测配筛选育成的优良光身稻恢复系。

本研究通过手指触感和扫描电镜观察进行光身表型鉴定，结果显示MH826的叶片和颖壳表面手指触摸光滑，无粗糙感。

扫描电镜同样未观察到明显的尖刺状表皮毛，呈现表面光滑表型，二者鉴定结果一致。

在明确了MH826光身表型的基础上，利用MH826与具有表皮毛特性的水稻品种华占、明恢2155和9311杂交，构建多个F2代分离群体对其光身性状进行遗传学分析。

结果表明MH826的光身性状由细胞核隐性单基因控制。

目前已报道的调控水稻光身性状的基因NUDA/GL-1［14］、OsWOX3B［26］、dep［9］、GLR1［13］、GL5［27］、gl1［28］、OsSPL10［29］、GLR3［17］、glr2［30］、GLL［20］等也均受细胞核隐性单基因控制。

本研究利用MH826与华占杂交构建的429个F2代隐性（光身性状）单株对GL826基因进行精细定位。

将GL826基因定位在水稻第5号染色体短臂端，位于分子标记InDel-70和InDel-11之间，相对遗传距离均为
0.2 cM，区间物理距离为141.03 kb。

参考RGAP数据库，发现在该精细定位区间共注释有23个完整的ORF。

其中LOC_Os05g02730基因［9，14］为前人已报道的与水稻表皮毛性状发育相关的基因。

因此，本研究重点将LOC_Os05g02730作为GL826的候选基因，并作进一步的序列比对和表达模式分析。

结果发现LOC_ Os05g02730基因在光身稻恢复系MH826苗期、分蘖期叶片、茎以及抽穗期叶片、茎和幼穗组织中的表达表现出不同程度的下调，说明LOC_Os05g02730基因在光身稻恢复系MH826中的表达受到了严重抑制。

Zhang 等［14］认为LOC_Os05g02730基因与WUSCHEL-like基因（OsWOX3B）同源。

OsWOX3B属于WOX3家族中的一个亚族，其主要调控水稻叶片形态的发生以及表皮毛的生长和发育，同时影响水稻侧生器官的发育和形成。

序列分析结果表明，LOC_Os05g02730与水稻的OsWOX3以及拟南芥的PRS、玉米的NS1和NS2具有高
表3　23个候选基因功能注释
Table 3　Functional annotations of twenty-three candidate genes of GL826 based on the RGAP database
ORF ORF1 ORF2 ORF3 ORF4 ORF5 ORF6 ORF7 ORF8 ORF9 ORF10 ORF11 ORF12 ORF13 ORF14 ORF15 ORF16 ORF17 ORF18 ORF19 ORF20 ORF21 ORF22 ORF23
基因名称
Gene name
LOC_Os05g02650
LOC_Os05g02670
LOC_Os05g02680
LOC_Os05g02690
LOC_Os05g02700
LOC_Os05g02710
LOC_Os05g02720
LOC_Os05g02730
LOC_Os05g02740
LOC_Os05g02750
LOC_Os05g02754
LOC_Os05g02760
LOC_Os05g02770
LOC_Os05g02780
LOC_Os05g02790
LOC_Os05g02810
LOC_Os05g02820
LOC_Os05g02840
LOC_Os05g02850
LOC_Os05g02860
LOC_Os05g02870
LOC_Os05g02880
LOC_Os05g02890
基因功能注释
Gene function annotation
表达蛋白
驱动蛋白
核因子相关的kappa-B结合蛋白
红光敏感性蛋白1
表达蛋白
反转座子蛋白
表达蛋白
同源盒结构域蛋白
表达蛋白
ABC转运蛋白
表达蛋白
表达蛋白
富含甘氨酸的细胞壁结构蛋白
富含甘氨酸蛋白A3
细胞周期蛋白
鸟苷酸结合蛋白
蛋白质相关转运蛋白
假定蛋白
反转录转座子蛋白
反转录转座子蛋白
白棕色复合体同源蛋白
富含精氨酸/丝氨酸剪接因子
白棕色复合体同源蛋白
418
3 期水稻明恢826光身基因GL826的精细定位与候选基因分析度相似性［9，30-31］
，其功能对于表皮细胞的分化至关重
要。

因此，综合以上表型鉴定、构建3个独立群体的遗传分析、精细定位、LOC_Os05g02730基因的表达模式分析及其功能注释等分析结果，本研究认为候选基因LOC_Os05g02730（OsWOX3B ）即为GL826基因。

然而本研究在进行基因序列比对时发现，LOC_Os05g02730基因在光身稻恢复系MH826和具有表皮毛特性的水稻品种华占、明恢2155和9311之间均未检测到DNA
水平上的差异，包括CDS 编码区和基因上游的启动子区域（图3）。

这与NUDA/GL-1［14］、OsWOX3B ［26］、dep ［9］报道的研究结果一致，即LOC_Os05g02730基因在有毛水稻和光身水稻品种间不存在DNA 水平上的差异，因而推测其光身表型可能由表观遗传修饰引起。

其中关于DNA 甲基化的研究较多，如水稻甲基化转移酶SDG714的功能缺失，会导致基因组甲基化水平明显降
低，进而导致植株表现出光叶和光壳表型［32］。

无毛非
注：A ： LOC_Os05g02730 CDS 编码区序列比对分析；B ：LOC_Os05g02730启动子区序列比对分析。

Note ： A ： Comparisons of coding sequences of the candidate gene LOC_Os05g02730. B ： Comparisons of promoter sequences of the candidate gene LOC_Os05g02730.
图3 候选基因LOC_Os05g02730在华占、明恢2155、9311和MH826中的序列比对分析
Fig.3　Sequences alignment of the candidate gene LOC_Os05g02730in Huazhan （HZ ）， Minghui2155 （MH2155）， 9311 and MH826
419
核农学报38 卷
洲栽培稻的颖壳无毛性状是由于hAT 转座子插入到控制毛状体发育基因GLAG5（GLABROUS GLUME 5）上，引发DNA 甲基化，进而抑制该基因的表达而产生无毛性状［16］。

另外，组蛋白乙酰化也被发现参与调控植物表皮毛的发育，如GCN5（GENERAL CONTROL NON DEREPRESS-IBLE 5）通过组蛋白乙酰化作用于GL1、GL2和GL3的启动子，参与抑制拟南芥的表皮毛发育
起始［33］。

miRNA 作为一类进化上保守的小分子RNA ，已有研究表明其参与了植物表皮毛的发育，如在拟南芥中，SPL （SQUAMOSA PROMOTER BINDING PROTEIN LIKE ）是调控表皮毛发育的一类转录因子，miRNA156通过调控SPLs 的表达来调控表皮毛发育，miRNA156的参与会导致SPLs 的表达量增加，拟南芥主茎上表皮毛的数量相应减少［34］。

另有研究报道，
水稻光身突变
注：**表示在P <0.01水平上差异极显著；误差线表示标准差（SD ）。

Note ： ** indicates significant differences at 0.01 level. Error bars represent the standard deviation （SD ）.
图4 LOC_Os05g02730在华占和MH826不同时期不同组织中的表达模式分析
Fig.4　Expression patterns of LOC_Os05g02730 in various tissues of huazhan and MH826 during different growth stages ，
determined by qRT-PCR analysis
表4　用于LOC_Os05g02730基因测序和表达分析的引物
Table 4　Primers used for sequencing and expression analysis of LOC_Os05g02730 in this study
引物名称 Primer name
GL730-S1-F GL730-S1-R GL730-S2-F GL730-S2-R GL730-S3-F GL730-S3-R
GL730-Pro1-F GL730-Pro1-R GL730-Pro2-F GL730-Pro2-R GL730-Q-F GL730-Q-R OsActin -F
OsActin -R
引物序列 Primer sequence （5′→3′）TTCTCGTATATCCCGGGGTG
AGCTAGCTAGGAGTAACTCTGA AGTTGAGTGTTGCCCTCCAA AAGCAAATTAAACGGCCGGA TCTGTCTCCTCACCGGCAAT TCGACCCATCCATCCATCAT TATGCTCACCTGCCCAAAGA
CGCATCCGTAGCTCAGAAAC TGTGGCAGCTCTCAAGTACA ACAAGCAACCTTCTCCTCCA
CGCAGATACAGCAGATCACG GTGGTGGTGGTAGTAGTGGG ACCTTCAACACCCCTGCTAT
CACCATCACCAGAGTCCAAC 用途 Usage
LOC_Os05g02730基因CDS 测序LOC_Os05g02730基因CDS 测序LOC_Os05g02730基因CDS 测序LOC_Os05g02730基因启动子测序LOC_Os05g02730基因启动子测序LOC_Os05g02730基因qRT-PCR 分析内参基因OsActin qRT-PCR 分析
注：F ： 正向引物；R ：反向引物。

Note ： F ： Forward primer. R ： Reverse primer.
420
3 期水稻明恢826光身基因GL826的精细定位与候选基因分析
体的光叶或光壳表型的形成可能由RNA的折叠等二级结构变化导致［5，28］。

由此可知，关于光身稻恢复系
MH826光身基因GL826的具体调控机理有待进一步研究。

前人对水稻光身基因的研究多采用光身突变体或远缘种质作为研究材料，而本研究则采用新育成的已在生产上应用的光身稻恢复系作为研究材料，这使得定位和克隆的基因更能与生产实际结合，能够直接在育种和生产实践中加以利用。

近年来，光身稻品种的选育和种质资源的应用越来越多。

张世辉等［35］利用光壳遗传背景材料AF6-82/R80与三合占/农香16杂交，经多代自交、回交，育成香型优质稻品种农香32。

陈锦文等［36］以光身稻澳大利亚-112与明恢77杂交，经多世代选育，育成光身稻恢复系泉恢404，并配组育成杂交稻组合福农优404。

本课题组利用光身稻恢复系MH826配制的杂交稻新组合泸优826、广优826、明优826均表现产量高、米质好、千粒重大、抗倒伏性强等优点，应用前景较好。

4　结论
光身稻恢复系MH826的光身表型由单个细胞核隐性基因控制。

利用MH826与华占杂交构建的F2分离群体将GL826基因精细定位在水稻第5号染色体短臂端，位于InDel-70和InDel-11之间约141.03 kb区域内。

通过序列分析和表达模式分析，确定候选基因LOC_Os05g02730即为GL826基因。

参考文献：
［1］胡培松，唐绍清，罗炬，黄发松.美国光身稻品种的利用与超高产品种的选育［J］.作物学报， 1999， 25（1）： 32-38
［2］毛恒.4个新选育水稻光身光温敏核不育系的应用基础研究［D］.
长沙：湖南农业大学， 2020
［3］尹忠良，毛恒，熊跃东，易嘉杰，叶红兵，唐文帮，邓化冰.4个新选育水稻光身光温敏核不育系的主要农艺性状及其配合力研
究［J］.杂交水稻， 2022， 37（4）： 21-29
［4］刘文炳，张建新，罗维禄，林琼，池晓雯，蔡为铭.光身杂交稻选育研究进展［J］.杂交水稻， 2005， 20（1）：6-10
［5］柳蕾，孙健，吴殿星，舒小丽.水稻光叶性状及表皮毛发育研究进展［J］.核农学报， 2015， 29（11）：2110-2116
［6］Song S S， Liu B， Song J Q， Pang S H， Song T X， Gao S， Zhang Y，Huang H，Qi T C. A molecular framework for signaling crosstalk
between jasmonate and ethylene in anthocyanin biosynthesis， trichome
development， and defenses against insect herbivores in Arabidopsis［J］.
Journal of Integrative Plant Biology， 2022， 64（9）： 1770-1788［7］Hulskamp M. Trichomes［J］. Current Biology， 2019， 29（8）： R273-R274［8］Han G L， Li Y X， Yang Z R， Wang C F， Zhang Y Y， Wang B S.
Molecular mechanisms of plant trichome development［J］. Frontiers
in Plant Science， 2022， 13： 910228
［9］Angeles-Shim R B，Asano K，Takashi T，Shim J，Kuroha T，Ayano M，Ashikari M. A WUSCHEL-related homeobox 3B gene，depilous（dep）， confers glabrousness of rice leaves and glumes ［J］.
Rice， 2012， 5： 28
［10］李晨光. 水稻光壳基因GLR3、GL5的克隆及功能研究和叶片茸毛基因HL6的精细定位［D］. 北京：中国农业大学， 2016
［11］牟文斌. 水稻多蘖矮秆基因htd（t）和茸毛基因HL1的精细定位及候选基因分析［D］. 南京：南京农业大学， 2013
［12］Zeng Y H，Zhu Y S，Lian L，Xie H G，Zhang J F，Xie H A.
Genetic analysis and fine mapping of the pubescence gene GL6
inrice （Oryza sativa L.）［J］. Chinese Science Bulletin，2013，58（24）： 2992-2999
［13］柳蕾. 水稻白脉光叶性状的表征及基因定位［D］. 杭州：浙江大学， 2015
［14］Zhang H L， Wu K， Wang Y F， Peng Y， Hu F Y， Wen L， Han B，Qian Q，Teng S. A WUSCHEL-like homeobox gene，OsWOX3B
responses to NUDA/GL-1 locus in rice［J］. Rice， 2012， 5： 30［15］王东. 水稻光滑叶基因精细定位和品种间SSR多态性分析［D］.
成都：四川师范大学， 2007
［16］Chang L Q， Hu M， Ning J， He W， Gao J Y， Ndjiondjop M N， Fu Y C，Liu F X，Sun H Y，Gu P，Sun C Q，Zhu Z F. The genetic
control of glabrous glume during African rice domestication［J］.
Journal of Genetics and Genomics， 2022， 49（5）： 427-436
［17］宋海冰，汪斌，陈壬杰，郑小雅，于世波，兰涛.水稻“光身”突变体glr3的遗传分析及基因定位［J］. 遗传， 2016， 38（11）： 1011-1018［18］Ishida T，Kurata T，Okada K，Wada T. A genetic regulatory network in the development of trichomes and root hairs［J］. Annual
Review of Plant Biology， 2008， 59： 365-386
［19］张洪凯. 一份水稻光滑叶突变体gl6的遗传分析与基因鉴定［D］.
成都：四川农业大学， 2018
［20］董陈文华，张小玲，朱骞，熊海波，魏振飞，吕永刚，张利东，伍腾飞，李伟，吴超，陈丽娟，李东宣. 水稻光叶突变新基因的克隆和
亚细胞定位［J］. 分子植物育种， 2015， 13（4）： 716-726
［21］郑小雅. 水稻表皮毛发育特征及相关基因的功能分析［D］. 福州：福建农林大学， 2018
［22］王穆穆，何艳芳，郑永胜，王晖，王丽媛，王东建，张晗，李汝玉. 水稻落粒基因SH8的精细定位与克隆［J］. 作物学报， 2022，48（8）： 1948-1956
［23］Murray M G， Thompson W F.Rapid isolation of high molecular weight plant DNA［J］. Nucleic Acids Research， 1980， 8（19）：4321-4325［24］Michelmore R W，Paran I，Kesseli R V.Identification of markers linked to disease-resistance genes by bulked segregation analysis：A
rapid method to detect markers in specific genomic regions by using
segregation population［J］. Proceedings of the National Academy of
Sciences of the United States of America，1991，88（21）：9828-
9832
［25］韦新宇，曾跃辉，杨旺兴，肖长春，候新坡，黄建鸿，邹文广，许旭明. 利用CRISPR-Cas9技术编辑Badh2基因创制优质香型籼
稻三系不育系［J］. 作物学报， 2023， 49（8）： 2144-2159
421。