花生荚果发育过程中的microRNA鉴定与表达分析

胡冬秀 刘 浩 洪彦彬 梁炫强 陈小平

花生荚果发育过程中的microRNA鉴定与表达分析

胡冬秀**刘 浩**洪彦彬 梁炫强 陈小平*

广东省农作物遗传改良重点实验室 / 国家油料作物改良中心南方分中心 / 广东省农业科学院作物研究所, 广东广州 510640

花生(L.)是典型“地上开花、地下结果”的作物, 为从转录后调控水平解析此独特的果实发育现象, 本文应用small RNA测序技术研究荚果发育11个时期果壳及种子中的microRNA及其靶基因。通过测序分别获得212个已知的microRNA和112个新microRNA, 其中, 已知microRNA包括197个保守microRNA和15个花生特异microRNA, 新microRNA来自62个新的microRNA前体序列。表达分析发现, 67个microRNA及其靶基因在荚果发育的11个时期存在时空特异性表达, 部分microRNA的表达量积累阶段性调节果壳与种子的发育, 表明microRNA参与了花生荚果暗发育的整个过程。此外, 对28个microRNA与30个靶基因进行荧光定量PCR验证发现, microRNA和靶基因的表达趋势与测序结果基本一致。本研究通过对花生荚果不同发育时期的果壳和种子进行small RNA测序, 鉴定参与调控花生荚果膨大相关的microRNA, 为研究黑暗条件下植物果实发育的调控机制与花生遗传改良奠定理论基础。

花生; 荚果发育; small RNA测序; microRNA; 靶基因

花生作为“地上开花、地下结果”的作物, 其生殖发育过程与其他作物存在差异。花生开花受精后, 子房柄伸长并向地生长形成果针, 位于其顶端的受精胚珠随着果针的伸长被推入土壤中, 果针停止伸长, 胚胎发育在黑暗条件下重新启动, 最终膨大形成荚果。花生荚果发育是一个复杂的生物学过程, 不同水平的调控对花生荚果种子发育均至关重要。研究microRNA如何调控黑暗条件下花生种子发育, 对理解植物果实暗发育调控与作物遗传改良具有重要意义。

microRNA是一类内源性非编码小分子RNA, 在植物中通常由18~24个核苷酸组成, 通过靶向降解或翻译抑制在转录后水平调控基因的表达[1]。植物中非编码的microRNA基因, 在NOT2和CDC5等转录因子的作用下通过RNA聚合酶II转录生成具有茎环结构的初级转录物Pri-microRNA (Primary microRNA)。随后Pri-microRNA经DCL1、HY1和SE组成的切割复合体剪切产生前体pre-microRNA (microRNA前体)[2], pre-microRNA经DCL1与DCL4进一步剪切生成成熟的microRNA[3], 在甲基转移酶HEN1的作用下使其3'端发生2'-O-甲基化来增强其稳定性[4-5]。microRNA在转运蛋白HST作用下从细胞核进入细胞质中, 随后与AGO效应蛋白(argonaute proteins, AGOs)结合形成RNA介导的沉默复合体(RNA-induced silencing complexes, RISCs)中, 结合互补的靶向基因信使RNA后使其被剪切降解[6]。microRNA作为一种重要的转录后调控因子, 参与植物的生长发育、信号转导、胁迫响应、器官分化以及病原菌调控等多种生物学过程[7-14]。拟南芥中, FUL蛋白与ARF (auxin response factor)相互作用并结合到miR172前体序列的启动子区域, 通过抑制miR172下游靶基因与3表达进而促进果实发育[15]; 人为干预miR319的表达活性则导致拟南芥花瓣变短、雄蕊败育、种子缺失等现象[16]。单子叶植物中, 水稻OsIDD2与SLR1相互作用, 通过赤霉素途径调控miR396介导的细胞增殖过程[17]; 而OsmiR408剪切OsUCL8编码序列会影响细胞中的铜离子稳态、质体蓝素丰度与水稻光合效率, 从而增加穗部分枝和粒数[18]。花生受到黄曲霉侵染时, 种子通过降低ayh-miRNA160、ayh-miRNA164对抗病基因TIR-NBS-LRR表达量产生的负调控效应, 从而提高抗病蛋白的表达量以抵御黄曲霉污染[19]。花生子房柄中的ayh-miR160a、ayh-miR171n与ayh-miR156e通过调控其靶基因参与激素信号转导途径, 在花生荚果、种子膨大阶段起关键作用[20]。

花生(L.)是我国重要的油料与经济作物, 荚果暗发育是花生最显著的特征之一, 已有相关研究报道非编码小RNA参与了花生胚胎与荚果早期的发育[21-22], 但对其荚果种子发育的整个过程缺乏全面的了解, 尤其对microRNA方面的相关研究仍然较为欠缺。本研究以‘粤油7号’为试验材料, 对其荚果发育不同时期的材料进行small RNA测序, 并利用qRT-PCR技术验证测序的结果, 鉴定和筛选出花生荚果种子发育过程中相关的microRNA及其在花生荚果种子发育过程中的表达模式。

1 材料与方法

1.1 试验材料

花生品种粤油7号种植于广东省农业科学院白云基地。对粤油7号植株的自交花进行鉴定并用彩色线标记, 在花后第8天用彩色标签系好伸长的果针。花生果针入土前, 被标记为第1个时期, 命名为P0 (Pod 0); 荚果开始膨大后通过测量荚果与种子的直径大小, 将整个荚果发育划分为10个时期, 即P1~P10 (Pod1~Pod10), 其中果壳与种子分开取样。果壳样品以缩写SH (Shell)表示, 种子样本以缩写SD (Seed)表示。本研究共对11个时期, 包括1个地上(果针即将入土)和10个地下发育阶段(图1-A), 以开花后天数为参考, 荚果与种子直径为主要指标来划分生育期, 最后获得20份样品。

1.2 RNA的提取及small RNA测序

利用TRIzol (Invitrogen)法提取20个样品的总RNA, 使用NanoDrop (Thermo Scientific)和安捷伦2100生物分析仪(安捷伦)对每个样品的RNA质量进行评估分析; 利用Illumina small RNA-seq文库构建试剂盒制备small RNA测序文库, 利用Illumina HiSeq 2000平台产生原始测序数据。

1.3 small RNA测序结果分析

利用Fastx-Tool kit, 对small RNA序列的原始数据进行过滤筛选, 去除接头序列、低质量序列, 包括冗余碱基或未被识别的碱基冗余序列, 以及小于18 nt的序列。将过滤后的数据利用Rfam 11.0数据库和GenBank非编码RNA数据库去除t/r/sn/snoRNA以及带有ployA尾巴的RNA, 对其他的small RNA序列利用microRNA19.0数据库进行BLASTn搜索, 以检索花生中已知的microRNA。此外, 利用Bowtie将small RNA序列映射到GenBank中的AHGI v2.0和花生GSS中进行比对, 完全匹配的序列将用于进一步分析。如果接近完全匹配的序列小于检索的small RNA或已知small RNA的长度, 则人工挑选检查与比较不匹配部分, 以确定匹配核苷酸的数量。为了有效识别花生中潜在的microRNA前体序列和新的microRNA, 用miREAP预测软件分析microRNA前体序列, 利用MiPred对由此产生的茎环结构进行过滤, 剩余的pre-microRNA的二级发夹结构使用CentroldFold软件进行评估。

1.4 花生荚果发育过程中microRNA趋势分析

microRNA测序数据以RPKM (reads per kilobase per million mapped reads)值表示, 将表达量数据归一化分析。差异表达的microRNA以RPKM > 10为筛选标准, 以此数据为基础利用TBtools软件绘制表达热图, 利用STEM软件对花生荚果发育过程中的microRNA和靶基因进行趋势分析, 并获得显著性和非显著性表达模式, 其中≤0.05的定义为显著表达模式。

1.5 茎环RT-PCR荧光定量检测microRNA

取5 µg总RNA加入DNase I去除基因组DNA污染, 利用TaKaRa SMART Scribea Reverse Transcriptase反转录酶对microRNA进行反转录。取2 µg的总RNA样品中加入终浓度为2 µmol L-1的microRNA的茎环引物, 将混合物在72℃预热3 min后置于冰上冷却使RNA变性, 打开二级结构。随后加5×First Strand Buffer、dNTP Mix、20 mmol L-1DTT、0.15 μL RRI (recombinant RNase inhibitor, 40 U μL-1)。加入1 mL (100 U μL-1) SMART Scribe RT再次混匀, 42℃孵育60 min, 反应结束后70℃加热15 min或加入4 µL 50 mmol L-1EDTA终止反应。荧光定量PCR方法参照试剂AceQ qPCR SYBR Green Master Mix说明书, 上述反转的cDNA样品浓度稀释为150 ng µL-1。荧光定量PCR的反应体系为20 µL, 包含10 µL SYBR Green Master Mix、2 µL模板、特异性正向引物和通用反向引物各0.4 µL、0.4 µL ROX Reference Dye I、6.8 µL的ddH2O。ABI Step OnePlus仪器检测, 循环反应为95℃ 10 s, 60℃ 30 s, 共40个循环;为内参基因, 用2-ΔΔCT法计算microRNA表达水平的相对定量, 每个样品设置3次生物学重复。

1.6 Quantative real time-PCR检测靶基因表达量

利用Prime Script RT试剂盒(TaKaRa, 中国大连)将1 µg总RNA反向转录为cDNA, 反应体系为20 µL, 采用ABI Step One Plus系统, 使用SYBR Premix Ex(TaKaRa, 中国大连)进行PCR反应, 用2-ΔΔCT法计算目标基因各组织的相对表达量, 以P0作为对照样品。作为内参基因, 每个样品设置3个重复, 数据用平均值±标准差(Mean± SD)表示。

2 结果与分析

2.1 花生荚果种子发育的表型特征

花生开花受精后, 子房柄不断伸长形成果针(P0)。随着果针的伸长, 位于其顶端的受精胚珠进入土壤中, 当果针刚进入地下时(P1), 其直径并未发生变化(图1-B)。在黑暗条件下, 果针停止伸长, 顶端迅速膨大形成荚果(P2~P6)。通过测量直径, 果壳在P6时期的直径大小是P2时期的5~6倍, 但是P6阶段过后, 果壳基本不再继续膨大, 因此进入土壤后的第6个时期是果壳发育逐渐停止的时期。种子发育在P2~P6阶段生长缓慢, 在P6期果壳不再发育之后, 种子则以较快的速度开始充实, 其直径从P6到最后的成熟期, 至少增大了3倍以上。表明, 果壳早于种子开始发育, 果壳的提前膨大为后续种子生长提供了空间和能量, 两者间发育的时间差有利于能量物质有效地供应种子生长。

A: 试验样品信息; B: 花生荚果种子发育的表型特征。P0: 地上果针; P1: 地下果针; P2~P5: 果壳膨大期; P6~P9: 果实充实期; P9~P10: 成熟期。

A: summary of experimental samples; B: phenotype characteristics of peanut pod (shell and seed) development. P0: aerial peg; P1: subterranean peg; P2–P5: pod expansion stage; P6–P9: seed filling stage; P9–P10: mature stage.

2.2 花生荚果发育过程中small RNA分析

对花生果壳和种子发育的20个样品进行small RNA测序, 共获得282,960,633个reads, 每个文库的reads数量从7,049,655到17,761,513。对接头序列和短序列进行修剪后, 共得到206,940,419个序列(图2-B), 长度大于17 nt。其中长度在18~30 nt的有177,581,859条, 有38,699,203条为花生中的特异small RNA。

microRNA属于微小的非编码RNA, 在转录后水平的基因调控中具有重要作用。为全面了解荚果发育过程中microRNA的全基因组信息, 从花生荚果种子发育过程中鉴定了长度为18~30 nt的small RNA序列(图2-A), 在21个和24个核苷酸处观察到了花生地下荚果种子发育(P1~P10)的2个主要高峰。果壳文库(P2SH~P10SH)有2个接近的高峰, 而种子文库(P2SD~P10SD)则以24个核苷酸small RNA为主。值得注意的是, 气生荚果P0中19个核苷酸处出现1个额外的高峰, 进一步分析发现, 这19个核苷酸序列中的大多数来自于几个高度表达的microRNA。此外, 在最初的地下荚果P1中, 总的small RNA在27~30个核苷酸之间高度表达, 但24个核苷酸small RNA在所有样本中均占优势, 在种子和果壳的总small RNA中分别占74%和72%, 表明花生的small RNA是类型复杂且多样的非编码RNA。

A: 全部以及特异性small RNA的长度分布和频率; B: small RNA统计分析; C: small RNA的分类。P0: 地上果针; P1: 地下果针; P2~P5: 果壳膨大期; P6~P9: 果实充实期; P9~P10: 成熟期。

A: length distribution and frequency of total and unique small RNAs; B: summary statistics of small RNAs; C: the categories of peanut small RNAs. P0: aerial peg; P1: subterranean peg; P2–P5: pod expansion stage; P6–P9: seed filling stage; P9–P10: mature stage.

2.3 花生荚果发育过程中small RNA的鉴定

为进一步鉴定花生中保守的和新的特异性microRNA, 我们通过图3所示的策略分析small RNA。对接头和短序列进行过滤后, 共得到38,699,203条花生的特异small RNA, 随后去除t/r/sn/snoRNA共15,326,397条, 其中rRNA的数量最多, 其次是tRNA、snRNA和snoRNA数量相对较少(图2-C)。得到38,326,363条候选的small RNA, 鉴定到212个已知的microRNA, 包括39个microRNA家族的197个保守的microRNA和13个microRNA家族的15个特异microRNA (附表1)。经过逐步筛选, 本研究还获得62个新的microRNA二级前体结构以及其对应生成的112个新的成熟microRNA。将前体命名为ahy-MIR-s (1~62)。利用RNAfold在线软件(http://rna.tbi.univie.ac.at//cgi-bin/RNAWeb- Suite/RNAfold.cgi)预测其二级结构(图4), 所有前体都能形成茎环结构, 多数前体形成4~6个茎环, 此外, ahy-MIR-s23、ahy-MIR-s31和ahy-MIR-s57等3个前体形成了3个茎环, ahy-MIR-s2、ahy-MIR-s36和ahy-MIR-s56 3个前体有7个茎环, ahy-MIR-s45和ahy-MIR-s20分别形成1个和8个茎环。进一步研究发现, 112个新的成熟microRNA (表1)定位于62个的新的Pre-miRNA的3¢臂或5¢臂, 并将其命名为ahy-miR-s(1-62)-3/5P, 大小为19~24 nt, 其中ahy-miR-s8-3p为19 nt, ahy-miR-s3-5p、ahy-miR-s7- 5p和ahy-miR-s29-5p等7个microRNA为20 nt, ahy-miR-s1-3p、ahy-miR-s2-3p和ahy-miR-s3-3p等40个microRNA为21 nt, 序列长度为22 nt和23 nt的microRNA分别为12个和6个, ahy-miR-s6-5p、ahy-miR-s15-5p和ahy-miR-s16-3p等46个microRNA长度为24 nt。

2.4 花生荚果发育过程中microRNA的表达分析

对RPKM大于10的microRNA进行表达分析发现, 在花生荚果11个发育阶段共有67个microRNA差异表达。利用TBtools软件[23]进行聚类分析和绘制热图(图5)发现, microRNA在荚果发育过程动态表达, 且呈时期特异性。在种子中microRNA的表达可分为K1、K2和K3三个簇, K1簇中microRNA在种子发育中表达相对较低。在K2簇中microRNA在种子成熟期(P7~P10)高表达。例如, ahy-miR166e、ahy-miR-s21-5p和ahy-miR-s26-5p在P8具有较高的表达量, 而ahy-miR-s1-3p、ahy-miR319c、ahy- miR-s29-5p和ahy-miR-s53-5p等4个microRNA在P10表达显著高于其他时期; ahy-miR166h-3p、ahy-miR168a和ahy-miR-s12-5p等microRNA在P8~P10都呈现了较高的表达。K3簇中microRNA主要在P5~P7中表达, 其中ahy-miR319i、ahy- miR390a-5p和ahy-miR-s31-5p等3个microRNA在种子初期(P2)有较高表达。ahy-miR167-3p、ahy-miR- s36-5p和ahy-miR-s10-5p等3个microRNA在P4~P6均有较高表达。在种子发育过程中P3时期microRNA的表达都相对较低。microRNA在果壳中的表达, 聚类将其分为了5个簇, C1簇中microRNA在果壳P8~P10高度表达, 例如, ahy-miR319c、ahy-miR156b-5p和ahy-miR167e分别在P10、P9和P8时期高丰度表达。C2簇主要在果壳发育初期(P2)高丰度表达, C3簇中microRNA呈现阶段性表达, 主要在P5~P9期有较高的表达量, 例如ahy- miR167c在P7~P9均表达较高。C4簇中microRNA在果壳中呈现低表达, C5簇主要在P3呈现高表达。表达分析结果进一步表明, microRNA在果壳和种子中都以高度阶段性或时期特异性的方式表达。

花生荚果发育过程中microRNA的趋势分析结果显示, 花生荚果发育中的microRNA共分为10种表达模式(图6-A), 大约84%的microRNA聚类到两种模式, 其中有43个microRNA聚类为profile.9表达模式(图6-B), 13个microRNA聚类为profile.8表达模式, 说明花生荚果发育过程中多数microRNA归入上述2个显著表达模式。microRNA在花生荚果发育过程中的表达是动态的, 在profile.9表达模式中, 花生荚果发育初期(P2~P5期)中microRNA的表达在种子中呈现下调, 在果壳中表达也较低, 随后microRNA表达上调, 种子和果壳发育中期(P6~P7期) microRNA表达达到最高, 花生荚果发育后期microRNA表达下调。在profile.8表达模式中, 花生荚果发育过程中microRNA的表达均处于上调, 在种子发育P4和P7期表达最高, 而在果壳中P9期表达最高。在其他几个表达模式中, microRNA在花生荚果发育的不同时期的积累水平也存在一定的差异。

P0: 地上果针; P1: 地下果针; P2~P5: 果壳膨大期; P6~P9: 果实充实期; P9~P10: 成熟期。

P0: aerial peg; P1: subterranean peg; P2–P5: pod expansion stage; P6–P9: seed filling stage; P9–P10: mature stage.

2.5 花生荚果发育过程中靶基因的表达分析

microRNA在转录后水平调控其下游靶基因来发挥生物学功能, 对67个microRNA的靶基因进行预测分析, 按照错配数≤5进行筛选, 共得到552个靶基因, 其中ahy-miR482e的靶基因数目最多(32个), ahy-miR-s51-3p和ahy-miR-s58-5p的靶基因数目最少(仅1个)。大多数保守的microRNA的靶基因是MYB、ARF和BEL1-like等转录因子, 其他microRNA的靶基因包括NB-LRR型抗病蛋白、F-box/FBD/LRR重复蛋白和锌指蛋白等, 此外还有部分未知功能的靶基因。为了进一步了解所预测靶基因的可能功能, 对其进行表达模式分析, 如图7所示, 其靶基因在花生果壳和种子发育过程中也呈现了高度阶段性和特异性的表达方式。例如,为BEL1-like转录因子, 在种子的P6期高丰度表达, 而和分别在果壳的P4和P3期高丰度表达。编码了NB-LRR抗病蛋白, 在种子中表达量较高, 而在果壳中几乎不表达。介导丝氨酸羟甲基转移酶(SHMT)途径, 在荚果发育的后期, 即种子和果壳的P7~P10期均有较高的表达。值得注意的是,在荚果发育的整个过程中均高丰度表达, 但对其功能有待研究。

P0: 地上果针; P1: 地下果针; P2~P5: 果壳膨大期; P6~P9: 果实充实期; P9~P10: 成熟期。

P0: aerial peg; P1: subterranean peg; P2–P5: pod expansion stage; P6–P9: seed filling stage; P9–P10: mature stage.

对花生荚果发育过程中靶基因表达趋势进行分析表明, 花生荚果发育过程中靶基因共有20个表达模式(图8)。其中有profile.7、profile.11、profile.18和profile.19等4个显著表达模式, 在profile.7表达模式中, 基因在种子P6~P9中高丰度表达, 在种子成熟期基因表达逐渐下调, 而在果壳中表达相对较低。在profile.11表达模式中, 基因在种子和果壳的发育早期表达最高。在profile.18中基因在荚果发育的过程中均处于上调表达, 在种子P3和P6以及果壳的P9时期表达量最高。在profile.19中, 种子发育早期基因表达量较低, 在P6时期基因上调表达, 而在果壳发育早期基因表达量较高, P8期开始下调表达。

2.6 RT-PCR实时定量验证

为了验证测序结果数据的准确性, 选取28个microRNA和30个靶基因分别进行茎环RT-PCR和实时荧光定量PCR检测发现, microRNA和靶基因的表达趋势与测序分析结果基本一致(图9)。例如, 在microRNA中(图9-A), ahy-miR-s29-5p和ahy- miR166d分别在种子P10和果壳P3高丰度表达; ahy-miR2118a在果壳中高丰度表达, 而在种子中表达相对较低; ahy-miR156e在种子发育过程中从P3时期开始有表达, P7时期表达量最高, 随后逐渐降低, 而在果壳中几乎不表达。在靶基因表达中(图9-B),编码的AUX1蛋白基因在果壳P10时期的表达显著高于其他时期,为GATA转录因子, 在种子的P8时期高表达,为乙醇酸氧化酶基因, 该基因在种子P7~P9时期均高丰度表达, 呈现阶段性的表达方式。总体而言, 荧光定量检测结果表明, microRNA及其靶基因的表达量与测序结果基本一致, 该结果为后续深入分析、克隆microRNA及其靶基因调控花生荚果发育奠定了基础。

P0: 地上果针; P1: 地下果针; P2~P5: 果壳膨大期; P6~P9: 果实充实期; P9~P10: 成熟期。

P0: aerial peg; P1: subterranean peg; P2–P5: pod expansion stage; P6–P9: seed filling stage; P9–P10: mature stage.

3 讨论

为了解花生荚果发育过程的表型特征, 对花生荚果的1个地上发育时期和10个地下发育时期的研究发现, 果针向地生长过程中(即P0~P1时期)不断伸长但其直径几乎不变, 果针入土后停止生长, 胚胎恢复发育, 荚果迅速膨大, 即P2~P6时期是荚果体积的快速增长阶段, 而种子的发育相对滞后, 即在这一时期生长相对缓慢。在P6时期后荚果大小基本不变, 种子则进入快速生长阶段。因此, P6时期既是荚果停止膨胀的分界点, 也是种子进入快速发育阶段的分界点。

本研究获得了62个新的Pre-microRNA和112个新的成熟的microRNA, 对Pre-microRNA进行二级结构预测发现, 所有前体均具有典型的茎环结构, 且多数序列(A+U)含量高于(G+C)含量, 说明鉴定得到的microRNA前体结构较稳定。成熟的microRNA在茎环结构的3¢臂或5¢臂上, 并能检测到其表达。此外, 对microRNA进行表达模式分析表明, microRNA在荚果发育过程中的表达是动态的, 且呈现阶段性和特异性表达。例如, ahy-miR390a-5p在种子P2时期表达较高, 在甘蓝型油菜中, miR390通过调控靶基因的表达介导早期胚胎发育, 从而影响其种子的含油量, 推测ahy-miR390a-5p可能参与花生荚果发育早期胚胎的发育[24]。ahy-miR156d和ahy-miR156e在花生种子成熟期(P7~P9)的表达显著高于其他时期, 相关研究表明, miR156通过靶向和基因, 能够导致不正常的细胞分裂[25], 从而控制种子的发育[26-27]。ahy-miR397a在P0期高度表达, 在水稻中, miR397的过表达通过靶向L-抗坏血酸氧化酶(AO)[28], 抑制抗坏血酸(AA)的氧化, 从而使脱氢抗坏血酸(DHA)维持在较低水平, 进而促进细胞分裂[29]。而种子在早期发育阶段需要具备很高的细胞分裂能力, 推测ahy-miR397a可能通过抗氧化途径而调控花生种子的发育过程。ahy-miR167e在种子P4~P7时期表达显著高于其他时期, ahy-miR167-5p在种子发育过程中均高表达, 在拟南芥和水稻中, miR167通过调控靶基因和的表达参与种子的发育[30-31], 因此ahy- miR167可能在花生荚果发育过程中发挥重要作用。microRNA319家族(ahy-miR319a、ahy-miR319b、ahy- miR319e和ahy-miR319f)在P0时期高度表达, miR319在花生荚果发育过程中参与生长素和GA信号传导途径, 能够促进花生胚胎和荚果早期发育[22]。

P0: 地上果针; P1: 地下果针; P2~P5: 果壳膨大期; P6~P9: 果实充实期; P9~P10: 成熟期。

P0: aerial peg; P1: subterranean peg; P2–P5: pod expansion stage; P6–P9: seed filling stage; P9–P10: mature stage.

4 结论

花生荚果种子发育过程中, microRNA以时空特异的方式表达, 表明不同的microRNA在花生荚果、种子发育的不同时期发挥作用, 该结果为阐释花生荚果暗发育的分子机制提供了转录后调控的理论参考依据。

附表 请见网络版: 1) 本刊网站http://zwxb.chinacrops.org/; 2) 中国知网http://www.cnki.net/; 3) 万方数据http://c.wanfangdata.com.cn/Periodical-zuo wxb.aspx。

[1] Rogers K, Chen X. Biogenesis, turnover, and mode of action of plant microRNAs., 2013, 25: 2383–2399.

[2] Chen X. A microRNA as a translational repressor of APETALA2 inflower development., 2004, 303: 2022–2025.

[3] Dong Z, Han M H, Fedoroff N. The RNA-binding proteins HYL1 and SE promote accurateprocessing of pri-miRNA by DCL1., 2008, 105: 9970–9975.

[4] Yu B, Yang Z Y, Li J J, Minakhina S, Yang M C, Padgett R W, Steward R, Chen X M. Methylation as a crucial step in plant microRNA biogenesis., 2005, 307: 932–935.

[5] Park W, Li J J, Song R T, Messing J, Chen X M. CARPEL FACTORY, a dicer homolog, and HEN1, a novel protein, act in microRNA metabolism in., 2002, 12: 1484–1495.

[6] Xie M, Zhang S X, Yu B. microRNA biogenesis, degradation and activity in plants., 2015, 72: 87–99.

[7] 邢利娟, 刘悦萍, 王磊, 徐妙云. miRNA参与植物胚和胚乳发育调控的研究进展. 生物技术进展, 2020, 10(2): 109–116. Xing L J, Liu Y P, Wang L, Xu M Y. Research progress of miRNA involved in regulation of plant embryo and endosperm development., 2020, 10(2): 109–116 (in Chinese with English abstract).

[8] 秦耀旭, 张关元, 刘司奇, 刘洋, 许志茹. 植物重金属胁迫相关miRNA的研究进展. 分子植物育种, 2019, 17: 2855–2861.Qin Y X, Zhang G Y, Liu S Q, Liu Y, Xu Z R. Research progress of miRNA related to heavy metal stress in plants., 2019, 17: 2855–2861 (in Chinese with English abstract).

[9] 郁佳雯, 韩荣鹏, 仇婷, 李宗芸, 潘沈元, 董婷婷. microRNA在植物生长发育中的研究进展. 分子植物育种, 2020, 18: 1496–1504. Yu J W, Han R P, Qiu T, Li Z Y, Pan S Y, Dong T T. Research progress on microRNA in plant growth and development., 2020, 18: 1496–1504 (in Chinese with English abstract).

[10] Long R C, Li M N, Li X, Gao Y L, Zhang T J, Sun Y, Kang J M, Wang T H, Cong L L, Yang Q C. A novel miRNA sponge form efficiently inhibits the activity of miR393 and enhances the salt tolerance and ABA insensitivity in., 2017, 35: 409–415.

[11] 李科学, 曲德杰, 黄慧梅, 曹金山, 樊仔慧, 李国纪, 王幼宁. miRNA调控大豆根系结瘤及共生固氮的分子机制研究进展. 植物生理学报, 2019, 55: 1587–1594. Li K X, Qu D J, Huang H M, Cao J S, Fan Z H, Li G J, Wang Y N. Research progress on miRNA-mediated molecular mechanisms of nodulation and symbiotic nitrogen fixation in soybean., 2019, 55: 1587–1594 (in Chinese with English abstract).

[12] Wang T Y, Ping X K, Cao Y R, Jian H J, Gao Y M, Wang J, Tan Y C, Xu X F, Lu K, Li J N, Liu L Z. Genome-wide exploration and characterization of miR172/euAP2 genes inL. for likely role in flower organ development., 2019, 19: 336.

[13] Jiao X M, Wang H C, Yan J J, Kong X Y, Liu Y W, Chu J F, Chen X Y, Fang R X, Yan Y S. Promotion of BR biosynthesis by miR444 is required for ammonium-triggered inhibition of root growth., 2020, 182: 1454–1466.

[14] Guo G H, Liu X Y, Sun F L, Cao J, Huo N, Wuda B, Xin M M, Hu Z R, Du J K, Xia R, Rossi V, Peng H R,Ni Z F, Sun Q X, Yao Y Y. Wheat miR9678 affects seed germination by generating phased siRNAs and modulating abscisic acid/gibberellin signaling., 2018, 30: 796–814.

[15] José Ripoll J, Bailey L J, Mai Q A, Wu S L, Hon C T, Chapman E J, Ditta G S, Estelle M, Yanofsky M F. microRNA regulation of fruit growth., 2015, 1: 15036.

[16] Todesco M, Rubio-Somoza I, Paz-Ares J, Weigel D. A collection of target mimics for comprehensive analysis of microRNA function in., 2010, 6: e1001031.

[17] Lu Y Z, Feng Z, Meng Y L, Bian L Y, Xie H, Mysore K S, Liang J S. SLENDER RICE1 andINDETERMINATE DOMAIN2 regulating osmiR396 are involved in stem elongation., 2020, 182: 2213–2227.

[18] Zhang J P, Yu Y, Feng Y Z, Zhou Y F, Zhang F, Yang Y W, Lei M Q, Zhang Y C, Chen Y Q. MiR408 regulates grain yield and photosynthesis via a phytocyanin protein., 2017, 175: 1175–1185.

[19] Zhao C, Li T, Zhao Y, Zhang B, Li A, Zhao S, Hou L, Xia H, Fan S, Qiu J, Li P, Zhang Y, Guo B, Wang X. Integrated small RNA and mRNA expression profiles reveal miRNAs and their target genes in response togrowth in peanut seeds., 2020, 20: 215.

[20] Ma X L, Zhang X G, Zhao K K, Li F P, Li K, Ning L L, He J L, Xin Z Y, Yin D M. Small RNA and degradome deep sequencing reveals the roles of microRNAs in seed expansion in peanut (L.)., 2018, 9: 349.

[21] Chen H, Yang Q, Chen K, Zhao S S, Zhang C, Pan R L, Cai T C, Deng Y, Wang X J, Chen Y T, Chu W T, Xie W P, Zhuang W J. Integrated microRNA and transcriptome profiling reveals a miRNA-mediated regulatory network of embryo abortion under calcium deficiency in peanut (L.)., 2019, 20: 392.

[22] Gao C, Wang P F, Zhao S Z, Zhao C Z, Xia H, Hou L, Ju Z, Zhang Y, Li C S, Wang X J. Small RNA profiling and degradome analysis reveal regulation of microRNA in peanut embryogenesis and early pod development., 2017, 18: 220.

[23] Chen C J, Chen H, Zhang Y, Thomas H R, Frank M H, He Y H, Xia R. Tbtools: an integrative toolkit developed for interactive analyses of big biological data., 2020, S1674-2052: 30187–30188.

[24] Zhao Y T, Wang M, Fu S X, Yang W C, Qi C K, Wang, X J. Small RNA profiling in twocultivars identifies microRNAs with oil production-and development-correlated expression and new small RNA classes., 2012, 158: 813–823.

[25] Nodine M D, Bartel D P. MicroRNAs prevent precocious gene expression and enable pattern formation during plant embryogenesis., 2010, 24: 2678–2692.

[26] Palatnik J F, Allen E, Wu X, Schommer C, Schwab R, Carrington J C, Weigel D. Control of leaf morphogenesis by microRNAs., 2003, 425: 257–263.

[27] Wang S K, Wu K, Yuan Q B, Liu X Y, Liu Z B, Lin X Y, Zeng R Z, Zhu H T, Dong G J, Qian Q, Zhang G Q, Fu X D. Control of grain size, shape and quality by OsSPL16 in rice., 2012, 44: 950–954.

[28] Pignocchi C, Kiddle G, Hernández I, Foster S J, Asensi A, Taybi T, Barnes J, Foyer C H. Ascorbate oxidase-dependent changes in the redox state of the apoplast modulate gene transcript accumulation leading to modified hormone signaling and orchestration of defense processes in tobacco., 2006, 141: 423–435.

[29] Potters G, Horemans N, Caubergs R J, Asard H. Ascorbate and dehydroascorbate influence cell cycle progression in a tobacco cell suspension., 2000, 124: 17–20.

[30] Yang J H, Han S J, Yoon E K, Lee W S. Evidence of an auxin signal pathway, microRNA167-ARF8-GH3, and its response to exogenous auxin in cultured rice cells., 2006, 34: 1892–1899.

[31] Jones-Rhoades M W, Bartel D P. Computational identification of plant microRNAs and their targets, including a stress induced miRNA., 2004, 14: 787–799.

Identification and expression analysis of microRNA during peanut (L.) pod development

HU Dong-Xiu**, LIU Hao**, HONG Yan-Bin, LIANG Xuan-Qiang, and CHEN Xiao-Ping*

Guangdong Provincial Key Laboratory of Crop Genetic Improvement / South China Peanut Sub-Center of National Center of Oilseed Crops Improvement / Crops Research Institute, Guangdong Academy of Agricultural Sciences, Guangzhou 510640, Guangdong, China

“Aerial flower and subterranean fruit” is a distinct feature in peanut (L.). To dissect this character at post-transcription level, small RNA sequencing was performed to identify microRNAs in peanut pod shell and seed during eleven developmental stages. Sequencing analysis identified 212 known microRNAs, including 197 conserved and 15 specific microRNA. In addition, 112 novel microRNAs from 62 novel microRNA precursors were identified. Among the known and new microRNAs, 67 microRNAs and their target genes showed differentially expressed patterns during peanut pod development. Expression trend analysis revealed stage-specific and tissue-specific expression of microRNA and their target genes during pod shell and seed development, implying that microRNAs probably played a role in peanut pod development. To validate expression profiles from small RNA sequencing, quantitative real-time RT-PCR were performed using 28 microRNAs and 30 target genes, revealing consistent expression profiles with sequencing results. The data regarding microRNA and their target genes generated in this study would contribute to understanding the molecular mechanism of plant fruit development under darkness and to crop improvement.

peanut; pod development; small RNA sequencing; microRNA; target gene

10.3724/SP.J.1006.2021.04144

本研究由广东省重点研发计划项目(2020B020219003), 国家自然科学基金项目(31771841)和广东省基础与应用基础研究基金(2020A1515010021)资助。

This study was supported by the Guangdong Provincial Major Research and Development Project (2020B020219003), the National Natural Science Foundation of China (31771841), and the Guangdong Basic and Applied Basic Research Foundation (2020A1515010021).

陈小平, E-mail: chenxiaoping@gdaas.cn

**同等贡献(Contributed equally to this work)

胡冬秀, E-mail: HuDX1017@163.com; 刘浩, E-mail: liuhao@gdaas.cn

2020-07-02;

2020-10-14;

2020-11-18.

URL: https://kns.cnki.net/kcms/detail/11.1809.S.20201118.1557.006.html