摘要: 目的深入挖掘分析亚麻茎秆不同发育时期木质素积累相关的miRNA及其靶基因,为解析亚麻纤维层木质素动态积累提供理论依据.方法以双亚4号(低木质素)和NEW(高木质素)花后20、30和40 d的茎秆为试验材料,基于miRNA测序和降解组测序结果分析miRNA转录水平动态变化及其靶基因注释功能,并采用实时荧光定量PCR检测2个亚麻品种间表达差异最显著的miRNA及其靶基因的表达模式.结果构建双亚4号和NEW花后20、30和40 d茎秆的miRNA文库,从中鉴定出305个miRNA,包括143个保守的miRNA和162个新鉴定的miRNA,且这些miRNA的靶基因不同转录本并非均受miRNA调控.将鉴定的305个miRNA与亚麻基因组序列信息进行靶基因预测,结果鉴定出286个miRNA的4807个靶基因,另外19个miRNA未预测到靶基因.通过对降解组测序数据进行整合分析,共获得21个miRNA的97个靶基因,共检测到97个降解位点,仅占鉴定出亚麻茎秆中miRNA靶基因总数(4807个)的2%,仍有大量的靶基因未被鉴定.对降解组中被降解的97个靶基因进行功能注释,发现有22个与植物生长发育有关,16个为转录因子,7个与植物次生代谢物积累有关.结合miRNA和降解组的测序结果,发现ama-miR156(Lus10003126)、bdi-miR397b-5p(Lus10017175)、lja-miR397(Lus10027782)、csi-miR160(Lus10021467)、aly-miR319c-p(Lus10008685)和gma-miR369h(Lus10011558)在双亚4号和NEW茎秆不同发育时期中最显著差异表达.实时荧光定量PCR检测结果显示,除gma-miR369h与其靶基因Lus10011558随亚麻茎秆的生长发育表达模式无明显规律外,其余5个miRNA与其靶基因表达模式恰好相反.结论亚麻茎秆木质素的积累过程与转录因子MYB、漆酶、生长素响应因子等编码基因密切相关,且这些基因可被其相应的miRNA负调控,并参与亚麻茎秆中木质素的积累.