雌雄蕊原基分化期低温胁迫下2个小麦品种幼穗miRNA表达谱分析

为探讨春季低温(倒春寒)胁迫下抗倒春寒能力不同的小麦品种相关miRNA及其响应规律,提供miRNA的表达谱和差异表达miRNA的信息,揭示差异表达miRNA在小麦抗倒春寒中的作用,采用高通量测序(Small RNA-seq)技术对雌雄蕊原基分化期0℃低温胁迫72 h及未低温处理(对照)的矮抗58(AK58)和郑麦366(ZM366)幼穗进行测序,通过生物信息学分析,以padj0.05且|log_2(Fold_change)|≥1为条件,筛选低温胁迫下差异表达显著的miRNA,利用psRobot(v1.2)对差异表达miRNA靶基因进行预测,然后对差异表达显著miRNA的靶基因进行GO富集分析和KEGG Pathway分析。结果表明,小麦幼穗测序共鉴定到112类miRNA,分属于38个不同的MIRNA家族。AK58、ZM366分别鉴定出85,88个差异表达miRNA,AK58有27个miRNA表达有显著差异,其中23个显著上调表达,4个显著下调表达,差异表达极显著的miRNA为2个,分别是tae-miR9672b、tae-miR9672a-3p; ZM366有48个差异表达显著的miRNA,13个显著上调表达,35个显著下调表达,差异表达极显著的miRNA为4个,分别为tae-miR9672b、tae-miR9672a-3p、tae-miR6201、tae-miR9674b-5p。对低温胁迫下AK58与ZM366差异表达的miRNA进行靶基因预测,AK58 85个差异表达miRNA对应848个靶基因;ZM366 88个差异表达miRNA对应6 537个靶基因。GO功能富集分析结果显示,AK58差异表达miRNA靶基因有20个生物学过程、6个细胞成分、6个分子功能类别极显著富集(FDR0.01);ZM366有10个生物学过程、14个细胞成分、19个分子功能类别极显著富集(FDR0.01)。对2个小麦品种幼穗中差异表达miRNA靶基因进行KEGG代谢通路(Pathway)富集分析,AK58差异表达miRNA靶基因显著富集(P0.05)到RNA聚合酶,嘌呤代谢,嘧啶代谢,光合生物固碳途径,自噬,托烷、哌啶和吡啶生物碱的生物合成,核苷酸切除修复,错配修复,DNA复制,核糖体,烟酸和烟酰胺代谢11个代谢通路;ZM366差异表达miRNA靶基因显著富集(P0.05)在植物-病原互作,植物昼夜节律,缬氨酸、亮氨酸和异亮氨酸降解,细胞周期,苯丙氨酸代谢,苯丙烷生物合成,RNA聚合酶,植物激素信号转导,氰胺酸代谢,次生代谢产物的生物合成,光合生物中碳固定,丙酮酸代谢,类胡萝卜素生物合成13个代谢通路。miRNA调控mRNA参与小麦应答低温胁迫中,小麦可能通过miRNA171a、miRNA156、miRNA398等多个miRNA的介导调控来抵御倒春寒胁迫,植物昼夜节律代谢途径可能与低温胁迫密切相关。2个小麦品种抗倒春寒能力不同的原因可能与miRNA调控光形态建成转导途径和开花过程的靶基因表达模式不同有关。