中国农业科学 ›› 2021, Vol. 54 ›› Issue (2): 435-448.doi: 10.3864/j.issn.0578-1752.2021.02.018

• 畜牧·兽医·资源昆虫 • 上一篇    下一篇

蜜蜂球囊菌菌丝和孢子中长链非编码RNA的比较及潜在功能分析

陈华枝1(),王杰1(),祝智威1,蒋海宾1,范元婵1,范小雪1,万洁琦1,卢家轩1,郑燕珍1,付中民1,2,3,徐国钧1,陈大福1,2,3,郭睿1,2,3()   

  1. 1福建农林大学动物科学学院(蜂学学院),福州350002
    2福建农林大学蜂疗研究所,福州 350002
    3福建农林大学蜂产品加工与应用教育部工程研究中心,福州 350002
  • 收稿日期:2020-02-10 接受日期:2020-02-25 出版日期:2021-01-16 发布日期:2021-02-03
  • 通讯作者: 郭睿
  • 作者简介:陈华枝,E-mail: CHZ0720@outlook.com。|王杰,E-mail: wanglegejie@163.com
  • 基金资助:
    国家自然科学基金(31702190);国家现代农业产业技术体系建设专项资金(CARS-44-KXJ7);福建省自然科学基金(2018J05042);福建农林大学硕士生导师团队(郭睿);福建农林大学科技创新专项基金(CXZX2017342);福建农林大学科技创新专项基金(CXZX2017343);福建省大学生创新创业项目(202010389016);福建省大学生创新创业项目(202010389162)

Comparison and Potential Functional Analysis of Long Non-Coding RNAs Between Ascosphaera apis Mycelium and Spore

CHEN HuaZhi1(),WANG Jie1(),ZHU ZhiWei1,JIANG HaiBin1,FAN YuanChan1,FAN XiaoXue1,WAN JieQi1,LU JiaXuan1,ZHENG YanZhen1,FU ZhongMin1,2,3,XU GuoJun1,CHEN DaFu1,2,3,GUO Rui1,2,3()   

  1. 1College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, Fuzhou 350002
    2Apitherapy Research Institution, Fujian Agriculture and Forestry University, Fuzhou 350002
    3Engineering Research Center of Processing and Application of Bee Products of Ministry of Education, Fujian Agriculture and Forestry University, Fuzhou 350002
  • Received:2020-02-10 Accepted:2020-02-25 Online:2021-01-16 Published:2021-02-03
  • Contact: Rui GUO

摘要:

【背景】蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)是一种专性侵染蜜蜂幼虫的真菌病原,可引起成年蜜蜂数量和蜂群群势的急剧下降。长链非编码RNA(long non-coding RNA,lncRNA)是一类新近发现的非编码RNA,在表观遗传、细胞周期、剂量补偿等众多生命活动中发挥重要生物学功能。【目的】明确球囊菌菌丝和孢子中lncRNA的数量、种类和表达谱差异,并探究共有lncRNA、特有lncRNA和差异表达lncRNA(differentially expressed lncRNA,DElncRNA)在菌丝与孢子中的潜在功能。【方法】利用基于链特异性建库的lncRNA-seq技术对球囊菌的纯化菌丝(AaM)和纯化孢子(AaS)分别进行测序。根据FPKM(Fragment Per Kilobase of per Million mapped reads)法计算lncRNA在AaM和AaS中的表达水平。通过Venn分析筛选AaM与AaS的共有lncRNA和特有lncRNA。按照P≤0.05且|log2 fold change|≥1的标准筛选AaM vs AaS比较组中的DElncRNA。通过Blast工具将共有lncRNA、特有lncRNA和DElncRNA的上下游基因比对GO和KEGG数据库,以进行功能及通路注释。根据靶向结合关系构建共有lncRNA、特有lncRNA和DElncRNA的竞争性内源RNA(competing endogenous RNA,ceRNA)调控网络并利用Cytoscape软件进行可视化。利用RT-qPCR验证测序数据的可靠性。【结果】AaM和AaS分别测得108 614 646和105 675 408条原始读段(raw reads),经严格过滤得到107 780 032和104 621 402条有效读段(clean reads),Q20分别为98.76%和98.72%,Q30分别达到95.84%和95.78%。共鉴定到850个lncRNA。AaM和AaS的共有lncRNA为701个,二者的特有lncRNA分别为39和110个。上述共有lncRNA通过顺式作用调控3 992个上下游基因,它们涉及细胞进程、代谢进程和催化活性等42个功能条目,以及代谢途径、次生代谢物的生物合成和抗生素的生物合成等117条通路;AaM的特有lncRNA和AaS的特有lncRNA通过顺式作用分别调控243和672个上下游基因。AaM vs AaS比较组包含的255个DElncRNA通过顺式作用调控1 479个上下游基因,它们涉及代谢进程、细胞进程和催化活性等41个功能条目,以及代谢途径、次生代谢产物的生物合成和抗生素的生物合成等107条通路。从共有lncRNA、孢子特有lncRNA和DElncRNA中分别预测到41、5和13个微小RNA(microRNA,miRNA)的前体序列。调控网络分析结果显示,菌丝lncRNA、孢子lncRNA形成较为复杂的ceRNA调控网络;菌丝lncRNA可靶向结合8个miRNA,进而调控77个mRNA;孢子lncRNA可靶向结合7个miRNA,进而调控87个mRNA;2个DElncRNA(TCONS_00008630与TCONS_00009302)可靶向结合miR-4968-y,进而调控10个mRNA。RT-qPCR验证结果显示4个DElncRNA的差异表达趋势与测序结果一致,表明测序数据真实可靠。【结论】共有lncRNA、特有lncRNA和DElncRNA可能通过调控上下游基因的表达,作为miRNA的前体,以及充当ceRNA影响菌丝和孢子的物质和能量代谢、自噬、转录、MAPK信号通路、泛素介导的蛋白水解、蛋白酶体以及次生代谢产物的生物合成等生物学过程,从而调节球囊菌的生长、发育、生殖和致病性。

关键词: 蜜蜂球囊菌, 长链非编码RNA, 菌丝, 孢子, 竞争性内源RNA

Abstract:

【Background】Ascosphaera apis is a fungal pathogen that exclusively infects honeybee larvae, causing a sharp decrease of the population of adult honeybees and colony population. Long non-coding RNA (lncRNA), a recently discovered non-coding RNA (ncRNA), plays a vital biological role in various activities such as epigenetics, cell cycle and dose compensation.【Objective】This study aimed to clarify the differences of number, type and expression profile of lncRNAs between A. apis mycelium and spore, and investigate the potential role of the common lncRNAs, specific lncRNAs and differentially expressed lncRNAs (DElncRNAs). 【Method】The purified mycelia (AaM) and purified spores (AaS) of A. apis were respectively sequenced using strand specific library-based lncRNA-seq technology. The expression levels of lncRNAs in AaM and AaS were calculated using FPKM (Fragment Per Kilobase of per Million mapped reads) method. Common lncRNAs and specific lncRNAs were filtered out following Venn analysis. DElncRNAs within AaM vs AaS comparison group were screened out following the standard of P≤0.05 and |log2 fold change|≥1. Upstream and downstream genes of common lncRNAs, specific lncRNAs and DElncRNAs were aligned against GO and KEGG databases to obtain function and pathway annotations. The competing endogenous RNA (ceRNA) regulation networks of common lncRNAs, specific lncRNAs and DElncRNAs were constructed following target binding relationships, followed by visualization using Cytoscape software. RT-qPCR was performed to verify the reliability of the sequencing data.【Result】In total, 108 614 646 and 105 675 408 raw reads were gained from AaM and AaS, and after strict filtering, 107 780 032 and 104 621 402 clean reads were obtained, respectively, with Q20 of 98.76% and 98.72%, and Q30 of 95.84% and 95.78%. A total of 850 lncRNAs were identified. Seven hundred and one lncRNAs were shared by AaM and AaS, and there were 39 and 110 specific lncRNAs. Via cis function, these shared lncRNAs could regulate 3 992 upstream and downstream genes involving in 42 functional terms such as cellular process, metabolism process and catalytic activity; and 117 pathways such as metabolism pathway, biosynthesis of secondary metabolites and biosynthesis of antibiotics. Specific lncRNAs for AaM and AaS could respectively regulate 243 and 672 upstream and downstream genes. In AaM vs AaS comparison group, 255 DElncRNAs were identified and found to regulate 1 479 upstream and downstream genes, which were associated with 41 functional terms including cellular process, metabolism process and catalytic activity; and 107 pathways including metabolism pathway, biosynthesis of secondary metabolites and biosynthesis of antibiotics. Forty one, five and 13 miRNA precursors were predicted from common lncRNAs, specific lncRNAs and DElncRNAs. The result of regulatory network analysis showed the formation of ceRNA networks of mycelium lncRNAs and spore lncRNAs; ten lncRNAs in mycelium could bind to eight miRNAs, further targeting 77 mRNAs; while eight lncRNAs in spore could link to seven miRNAs, further targeting 87 mRNAs; two DElncRNAs including TCONS_00008630 and TCONS_00009302 could simultaneously target miR-4968-y, further regulating ten mRNAs. The result of RT-qPCR suggested the differential expression trend of four DElncRNAs were in accordance of that in sequencing result, indicating the reliability of our sequencing data.【Conclusion】Common lncRNAs, specific lncRNAs and DElncRNAs are likely to affect material and energy metabolisms, autophagy, transcription, MAPK signaling pathway, ubiquitin-mediated proteolysis, proteasome and biosynthesis of secondary metabolites, by regulating the expression of upstream and downstream genes, serving as miRNA precursors or ceRNAs, thus regulating the growth, development, reproduction and pathogenicity of A. apis.

Key words: Ascosphaera apis, long non-coding RNA (lncRNA), mycelium, spore, competing endogenous RNA (ceRNA)