新闻媒体-太阳成tyc7111cc

发稿时间：2019-10-17来源：天昊生物

当你手中有好的实验材料时，从全转录组（mrna、lncrna、mirna和circrna）角度入手，不失为一种高效的研究策略。刚刚在上月发表的一篇文章，就利用该策略系统阐述了油菜细胞质雄性不育系花药发育中的circrna表达模式及cerna和mirna-mrna间的调控网络。

发表期刊：int. j. mol. sci.  发表时间：2019-9-27  影响因子：4.183

植物的雄性不育，是指它们不能产生正常功能的花粉粒。细胞质雄性不育是一种重要的植物繁殖特性，可作为开发作物杂种优势的有用工具。之前研究发现了一些mrna和mirnas可以作为细胞质雄性不育系统的重要调节因子。环状rna(circrna)是一类新型共价闭合的单链内源性非编码rna，也可以通过竞争性内源rna(cernas)调控机制，参与转录后调控过程。

在本文中，研究者在油菜上建立了细胞质雄性不育系（sterile line）bcpol97-05a和可育系（fertile line）bcah97-01b，通过全转录组测序对不育系和可育系共6个样本（2类样本各3次重复）的花蕾进行了rna表达谱比较，共鉴定出31个差异表达的circrnas、47个差异表达的mirnas和4779个差异表达的mrnas (注：本文未提到lncrna相关信息，小编估计是作者发掘到很多感兴趣的候选lncrnas，留到后续研究)。利用上述结果，构建了mirna介导的调控网络和cerna网络，推测circrna a02:23507399|23531438是在转录后水平调控花药发育的重要circrnas。go分析表明，mirnas和circrnas可以调节纤维素、孢粉素、果胶和胰蛋白酶的有序分泌和沉积；脂质的及时降解；绒毡层细胞的程序性细胞死亡在花药发育中起着关键作用。本研究揭示了一个新的circrna-mirna-mrna网络，它参与了油菜花药的发育，丰富了对开花植物细胞质雄性不育的理解，为进一步研究circrna和mirna在花药发育中的作用奠定了基础。

实验材料及方法

本文材料为油菜雄性不育系bcpol97-05a，其保持系为bcajh97-01b。开花后，鉴定可育和不育植物，将单株植物花序进行取样，液氮速冻后-75℃保存。利用trizol提取花序的总rna。全转录组分析的总rna来自每个株系30棵单株植物的混合，共三个生物重复。在验证了rna样品的质量后，使用epicentre ribo-zero试剂盒去除rrna，之后分别构建各类rna文库。构建好文库并质控后，利用illumina hi-seq平台进行测序。使用find_circ等软件预测circrna。通过与mirbase中成熟mirnas的序列进行比对来鉴定已知的mirnas。使用mirdeep2软件预测了新的mirnas。后续还进行了目标基因预测、go分析和目标基因的kegg分析等，并进行了qrt-pcr验证，最终利用cytoscape构建mirna-mrna和circrna-mirna网络。

实验结果

油菜细胞质雄性不育系的建立

本实验以可育植株bcajh97-01b为轮回亲本，回交八代选育不育系，最终建立了雄性不育系bcpol97-05a，其保持系为bcajh97-01b。与可育系相比，bcpol97-05a的雄蕊似乎已经耗尽了色素，呈现出白化的样子。它们也枯萎了，没有成熟的花粉粒。此外，不育系bcpol97-05a的雄蕊似乎较短，并显示出相连的花药(图1)，其花粉的活力比可育系低，约85%的花粉粒是败育的(图2a，b)。dapi染色，也发现不育系的花粉形态异常，核发育异常(图2c-f)。实验还观察了两个品系花粉发育的表型特征。不育系的花药囊比可育系的小，最终成熟花粉粒比可育系中少得多(图3)。

图1、油菜可育系和不育系花的形态

图2、油菜可育系和不育系的花粉特性

图3、油菜可育系和不育系的花药横切面

不育系和可育系circrna的鉴定

通过全基因组测序，研究者从雄性不育系和可育系的花蕾中鉴定并注释了1443个circrnas。共有121个circrnas (图4a)。其转录本广泛分布在所有10条染色体中(图4b)。大约12.68%的circrnas来自chr3，10.46%来自chr6，14.00%来自chr9，而来自其他染色体的circrnas的百分比都小于10%(图4c)。circrnas分为三类，其中外显子环状911个(63.13%)，内含子circrnas 182个(12.61%)，外显子内含子circrnas 350个(24.26%；图5a)，其中大部分长度为200-400 bp(图5b)。

图4、油菜可育系和不育系中鉴定的circrna

图5、油菜可育系和不育系中circrna特征

可育系与不育系中差异rnas的鉴定

本研究在不育系中鉴定出31个decircrnas (9个上调和22个下调)和47个demirna (6个上调和41个下调)。我们分别在表1和表2中总结了基于log2fc的前10个脱氧核糖核酸和核糖核酸。此外，与可育系相比，研究者在不育系中检测到4779个demrnas(1021个上调和3758个下调)。

表1、与油菜可育系和不育系中差异表达的前10位circrnas

表2、油菜可育系和不育系差异表达的前10位mirnas

对demrnas的kegg途径富集分析表明，“淀粉和蔗糖代谢”、“苯丙素类生物合成”和“戊糖和葡糖醛酸相互转化”是最富集的代谢途径。研究还表明，大多数脱氧核糖核酸都被注释为与植物激素信号转导有关(图6)。

图6、油菜可育系和不育系差异表达基因的kegg途径富集分析

此外还对demrnas进行了go分析，分析了差异基因和所有基因二级功能的富集，反映了每个二级功能的状态(图7)。此外，我们分析了每一个go术语的富集情况，并总结了在“生物过程”类别、“细胞成分”类别和“分子功能”类别中最显著富集的go术语。在bp类别中，最重要的三个节点是植物型细胞壁修饰、花粉管生长和肌动蛋白丝基运动。大多数差异基因参与花粉壁发育，包括花粉外壁和胞质形成(表3)。

图7、油菜可育系和不育系差异基因go分析

表3、油菜可育系和不育系差异表达基因的bp类别前10位

此外，研究者使用相同的不育系和可育系样品，通过qrt-pcr证实了核糖核酸序列的结果。随机选择五个circrnas、八个mirnas和八个mrnas。定量的结果与测序结果一致，这证明了rna图谱的高度可靠性(图8)。

图8、rna测序数据的qrt-pcr验证

2.4 .半核糖核酸-半核糖核酸网络的构建

为了更好地理解雄性不育花药发育过程中的基因调控网络，本研究使用targetfinder软件鉴定了推测的demirna–demrna相互作用。结果总共获得了170个demirna–demrna相互作用。我们选择了18个mirnas和37个mrnas的demirna–demrna对，其中mirnas和相应的mrnas具有相反的表达，并利用cytoscape构建了mirna介导的调控网络(图9)。为了研究mirnas的潜在功能，对该网络中假定的靶基因进行了go分析(图10)。

图9、油菜可育系和不育系花药发育中的demirna-demrna网络图

图10、油菜可育系和不育系demirna-demrna网络中mirnas潜在靶标的go分类

decircrna–demirna–demrna网络的构建

通过targetfinder软件，在不育系中获得了两个由一个上调的circrna(a02:23507399|23531438)和两个下调的mirna(非保守_a02_5092和非保守_a07_27586)靶基因组成的decircrna–demirna相互作用。为了研究cerna在花药发育中的调控，并鉴定与花药发育相关的circrna，构建了decircrna–demirna–demrna网络。最终，通过分析，构建了一个cerna网络，由一个circrna、两个mirnas、五个mrnas(bra006799、bra013352、bra022668、bra029377和bra038700)和10个交互对组成(图11)。

图11、油菜可育系和不育系花药发育中的决定decircrna-demirna-demrna三重网络

结论

本研究通过全转录组测序，鉴定了几种与花药发育相关的rnas，包括mrna、circrna和mirna，还重建了涉及这些rna分子的假定调节网络，这将加深和丰富人们对花药发育和雄性生育控制的理解。此外，对cerna网络的解析，为花药发育过程中rna的调控作用提供新的视角。

近几年植物cerna相关文章：

全转录组测序技术优势：

• rrna去除建库，保留了完整的rna种类信息；
•  链特异性文库，可以保留转录本的链信息，更准确地检测反义rna；
• 使用高通量测序，能够获得更加全面的rna信息，包括低丰度的rna；
• 通过测定的序列信息精确地分析不同类型rna的表达丰度变化及其生物学功能；
• 分析同一样本中的mrna、lncrna、mirna和circrna四种类型rna，明确这些rna之间的共表达和调控关系，揭示cerna调控机制。

关于天昊：

天昊生物具有丰富的全转录组测序和qrt-pcr实验经验，我们致力于为研究者提供高质量的科研策略咨询、实验技术服务和遗传数据分析服务，期待成为大家科研工作中的“昊”助手与“昊”伙伴。欢迎联系太阳成tyc7111cc具体咨询！邮箱：techsupport@geneskies.com 电话：400-065-6886

欢迎联系太阳成tyc7111cc！邮箱：techsupport@geneskies.com电话：400-065-6886