蜜蜂球囊菌菌丝和孢子中微小RNA及其靶mRNA的比较分析

【目的】蜜蜂球囊菌(Ascosphaera apis,球囊菌)专性侵染蜜蜂幼虫而导致白垩病。本研究旨在通过small RNA-seq（sRNA-seq）技术和生物信息学方法对球囊菌纯化菌丝（AaM）和纯化孢子（AaS）进行深度测序和比较分析,明确球囊菌菌丝miRNA和孢子miRNA的数量、结构和表达谱差异,并揭示菌丝和孢子共有miRNA、特有miRNA和差异表达miRNA（differentially expressed miRNA,DEmiRNA）及其靶mRNA与球囊菌菌丝和孢子生长、发育和病原致病性的潜在关系。【方法】实验室条件下获得纯培养的球囊菌,利用sRNA-seq技术对AaM和AaS分别进行测序,通过对原始读段（raw reads）进行过滤和质控获得有效标签序列（clean tags）。通过Venn分析筛选菌丝和孢子共有miRNA和特有miRNA。根据P≤0.05且|log₂ fold change|≥1的标准筛选AaM vs AaS的DEmiRNA。对上述共有miRNA、特有miRNA和DEmiRNA的靶mRNA进行预测,并对靶mRNA进行GO及KEGG数据库注释。根据靶向结合关系构建DEmiRNA和靶mRNA的调控网络。利用RT-qPCR验证测序数据的可靠性。【结果】AaM和AaS中分别得到12 982 320和12 708 832条raw reads,经过滤和质控分别得到10 800 101和9 888 848条clean tags。AaM中miRNA的长度介于18—26 nt,AaS中miRNA的长度介于18—24 nt,分布miRNA数量最多的长度均为18 nt,AaM和AaS中首位碱基为U的miRNA数量最多。AaM和AaS中表达量最高的miRNA均为miR6478-x、miR10516-x和miR482-x。菌丝和孢子共有miRNA靶向结合5 946个mRNA,二者特有miRNA分别靶向结合6 141和6 346个mRNA。共有miRNA的靶mRNA主要参与代谢进程、细胞进程和催化活性等42个功能条目,以及翻译、碳水化合物代谢和能量代谢等120条通路。AaM vs AaS比较组包含93个DEmiRNA,可靶向结合6 090个mRNA,这些靶mRNA可注释到38个功能条目和120条通路。DEmiRNA与靶mRNA之间形成较为复杂的调控网络,miR-4968-y位于调控网络的中心且能够靶向结合多达118个mRNA。RT-qPCR结果显示5个DEmiRNA的表达趋势与测序数据一致,证实了本研究中测序数据的可靠性。【结论】球囊菌菌丝和孢子中的miRNA具有类似的结构特征,但表达谱表现出明显差异;菌丝和孢子可能通过特异性表达和差异表达部分miRNA对其生长、发育和生殖进行调控。     

蜜蜂球囊菌菌丝和孢子中微小RNA及其靶mRNA的比较分析

【目的】蜜蜂球囊菌(Ascosphaera apis,球囊菌)专性侵染蜜蜂幼虫而导致白垩病。本研究旨在通过small RNA-seq(sRNA-seq)技术和生物信息学方法对球囊菌纯化菌丝(AaM)和纯化孢子(AaS)进行深度测序和比较分析,明确球囊菌菌丝miRNA和孢子miRNA的数量、结构和表达谱差异,并揭示菌丝和孢子共有miRNA、特有miRNA和差异表达miRNA(differentially expressed miRNA,DEmiRNA)及其靶mRNA与球囊菌菌丝和孢子生长、发育和病原致病性的潜在关系。【方法】实验室条件下获得纯培养的球囊菌,利用sRNA-seq技术对AaM和AaS分别进行测序,通过对原始读段(raw reads)进行过滤和质控获得有效标签序列(clean tags)。通过Venn分析筛选菌丝和孢子共有miRNA和特有miRNA。根据P≤0.05且|log2 fold change|≥1的标准筛选AaM vs AaS的DEmiRNA。对上述共有miRNA、特有miRNA和DEmiRNA的靶mRNA进行预测,并对靶mRNA进行GO及KEGG数据库注释。根据靶向结合关系构建DEmiRNA和靶mRNA的调控网络。利用RT-qPCR验证测序数据的可靠性。【结果】AaM和AaS中分别得到12 982 320和12 708 832条raw reads,经过滤和质控分别得到10 800 101和9 888 848条clean tags。AaM中miRNA的长度介于18—26 nt,AaS中miRNA的长度介于18—24 nt,分布miRNA数量最多的长度均为18 nt,AaM和AaS中首位碱基为U的miRNA数量最多。AaM和AaS中表达量最高的miRNA均为miR6478-x、miR10516-x和miR482-x。菌丝和孢子共有miRNA靶向结合5 946个mRNA,二者特有miRNA分别靶向结合6 141和6 346个mRNA。共有miRNA的靶mRNA主要参与代谢进程、细胞进程和催化活性等42个功能条目,以及翻译、碳水化合物代谢和能量代谢等120条通路。AaM vs AaS比较组包含93个DEmiRNA,可靶向结合6 090个mRNA,这些靶mRNA可注释到38个功能条目和120条通路。DEmiRNA与靶mRNA之间形成较为复杂的调控网络,miR-4968-y位于调控网络的中心且能够靶向结合多达118个mRNA。RT-qPCR结果显示5个DEmiRNA的表达趋势与测序数据一致,证实了本研究中测序数据的可靠性。【结论】球囊菌菌丝和孢子中的miRNA具有类似的结构特征,但表达谱表现出明显差异;菌丝和孢子可能通过特异性表达和差异表达部分miRNA对其生长、发育和生殖进行调控。     

蜜蜂球囊菌中长链非编码RNA的调控作用

【目的】长链非编码RNA（long non-coding RNA,lncRNA）是一类二、三级结构高度保守,长度>200 nt且不具蛋白编码能力的RNA,在转录和转录后水平广泛参与调控剂量补偿、细胞分化和生长发育等生命活动。本研究基于前期获得的蜜蜂球囊菌（Ascosphaera apis,简称球囊菌）菌丝和孢子混合样品的高质量lncRNA组学数据进行球囊菌lncRNA的顺式（cis）作用、反义lncRNA（antisense lncRNA）作用和竞争性内源RNA（competing endogenous RNA,ceRNA）作用的分析和探讨,以期揭示lncRNA在球囊菌中的潜在功能。【方法】基于lncRNA基因在染色体上的位置,预测lncRNA上下游100 kb以内的蛋白编码基因;使用Blast软件将上下游基因比对到GO和KEGG数据库,以获得功能和通路注释。利用LncTar软件对反义lncRNA的靶mRNA进行预测,并使用Blast软件将上述靶mRNA比对到KEGG和eggNOG数据库。利用TargetFinder软件预测lncRNA靶向结合的miRNA及miRNA靶向结合的mRNA,根据靶向结合关系建立lncRNA-miRNA和lncRNA-miRNA-mRNA调控网络,进而通过Cytoscape v3.7.1软件进行调控网络的可视化。利用RT-PCR对调控网络中的lncRNA、靶miRNA和靶mRNA进行表达验证。【结果】共预测出371个lncRNA的5 852个上下游基因,可注释到细胞进程、代谢进程和应激反应等48个功能条目,以及新陈代谢途径、次生代谢产物的生物合成和抗生素的生物合成等121条通路,表明部分lncRNA可通过顺式作用调节上下游基因的表达,从而参与调控球囊菌的生长发育及物质能量代谢等基础生命活动。球囊菌的7个lncRNA与7个靶mRNA存在序列互补关系,其中5个mRNA在eggNOG数据库中仅注释为假定蛋白,仅gene3444在KEGG数据库注释为核孔复合体蛋白An-Nup120和假定蛋白,表明上述1个反义lncRNA可能参与调控核孔复合体蛋白的生物合成等生物学过程。此外,共预测出227个lncRNA与73个miRNA之间存在靶向结合关系,其中多数lncRNA(79.02%)仅能结合1—2个miRNA,部分miRNA可被多个lncRNA靶向结合;进一步构建氧化磷酸化通路和MAPK信号通路相关的lncRNA-miRNA-mRNA调控网络,分析结果显示氧化磷酸化通路相关的222个lncRNA靶向78个miRNA及50个mRNA,MAPK信号通路相关的222个lncRNA靶向76个miRNA及46个mRNA,表明部分lncRNA通过ceRNA作用调控此两条通路,从而影响球囊菌的能量合成、环境适应以及生长发育等过程。【结论】部分lncRNA可能通过顺式作用和ceRNA作用调节球囊菌的生长发育、物质能量代谢以及环境适应等生物学过程;MSTRG.5393.1可能作为反义lncRNA调控球囊菌的核孔复合体的蛋白合成。     

东方蜜蜂微孢子虫孢子中长链非编码RNA的竞争性内源RNA调控网络及潜在功能

【背景】 长链非编码RNA（long non-coding RNA,lncRNA）是一类长度>200 nt且不具有编码能力的转录本,能够在真核生物的转录水平和转录后水平通过顺式（cis）、反式（trans）或竞争性内源RNA（competing endogenous RNA,ceRNA）机制调控基因的表达,已被证明在生物体的生长和发育等生物学过程中发挥重要的调控作用。【目的】 结合测序得到的东方蜜蜂微孢子虫（Nosema ceranae）纯化孢子的small RNA（sRNA）组学数据和前期获得的lncRNA组学数据对lncRNA的调控方式进行深入细致的分析和探讨,揭示lncRNA在东方蜜蜂微孢子虫孢子中的潜在功能。【方法】 利用small RNA-seq（sRNA-seq）技术对东方蜜蜂微孢子虫纯化孢子进行测序,通过相关的生物信息学软件对测序数据进行质控。根据lncRNA与mRNA的位置关系预测lncRNA的上下游基因;使用Blast软件比对GO和KEGG数据库,对上下游基因进行功能和代谢通路注释;利用Target Finder软件预测lncRNA-miRNA及lncRNA-miRNA-mRNA调控网络,进而通过Cytoscape v3.7.2软件对调控网络进行可视化;利用RT-PCR对调控网络中的部分lncRNA、miRNA和mRNA的表达进行验证。【结果】 东方蜜蜂微孢子虫纯化孢子样品的sRNA-seq分别得到16 597 883、15 451 791和12 248 316条raw reads,经质控分别得到15 608 370、14 249 255和11 440 684条clean reads,各样品的Q30≥98.58%。共预测出lncRNA的上下游基因310个,它们可注释到代谢进程、细胞进程、催化活性、结合和细胞等相关的35个功能条目;此外还可注释到56条代谢通路,涉及嘌呤代谢、碳代谢和丙酮酸代谢等物质代谢通路,以及甲烷代谢、氧化磷酸化和糖酵解/糖异生等能量代谢通路。上述结果表明相应的lncRNA可能通过cis作用调节上下游基因的表达,从而参与调控东方蜜蜂微孢子虫孢子中的物质和能量代谢以及细胞生命活动。进一步构建lncRNA-miRNA调控网络,分析结果显示MSTRG.3636.1、MSTRG.4498.1和MSTRG.4883.1可靶向结合4个miRNA（nce-miR-7729、nce-miR-7502、nce-miR-8639和nce-miR-8565）,说明此3个lncRNA可作为ceRNA在东方蜜蜂微孢子虫孢子中发挥潜在作用。LncRNA-miRNA-mRNA调控网络分析结果显示,三者之间形成较为复杂的调控网络关系,miRNA处于调控网络的中心并联系lncRNA和mRNA,nce-miR-7502和nce-miR-8639结合的mRNA最多,达到28个;与MSTRG.3636.1、MSTRG.4498.1和MSTRG.4883.1存在靶向关系的miRNA结合多个mRNA,表明上述3个lncRNA可能通过ceRNA机制发挥作用,从而影响东方蜜蜂微孢子虫孢子中的新陈代谢和生命活动。【结论】 lncRNA可能通过cis作用调控上下游基因的表达,以及作为ceRNA通过吸附miRNA间接影响靶基因的表达,从而参与调节东方蜜蜂微孢子虫孢子中的物质和能量代谢等基本生命活动。     

东方蜜蜂微孢子虫孢子中长链非编码RNA的竞争性内源RNA调控网络及潜在功能

【背景】长链非编码RNA(long non-coding RNA,lncRNA)是一类长度200 nt且不具有编码能力的转录本,能够在真核生物的转录水平和转录后水平通过顺式(cis)、反式(trans)或竞争性内源RNA(competing endogenous RNA,ceRNA)机制调控基因的表达,已被证明在生物体的生长和发育等生物学过程中发挥重要的调控作用。【目的】结合测序得到的东方蜜蜂微孢子虫(Nosema ceranae)纯化孢子的small RNA(sRNA)组学数据和前期获得的lncRNA组学数据对lncRNA的调控方式进行深入细致的分析和探讨,揭示lncRNA在东方蜜蜂微孢子虫孢子中的潜在功能。【方法】利用small RNA-seq(sRNA-seq)技术对东方蜜蜂微孢子虫纯化孢子进行测序,通过相关的生物信息学软件对测序数据进行质控。根据lncRNA与mRNA的位置关系预测lncRNA的上下游基因;使用Blast软件比对GO和KEGG数据库,对上下游基因进行功能和代谢通路注释;利用Target Finder软件预测lncRNA-miRNA及lncRNA-miRNA-mRNA调控网络,进而通过Cytoscape v3.7.2软件对调控网络进行可视化;利用RT-PCR对调控网络中的部分lncRNA、miRNA和mRNA的表达进行验证。【结果】东方蜜蜂微孢子虫纯化孢子样品的sRNA-seq分别得到16 597 883、15 451 791和12 248 316条raw reads,经质控分别得到15 608 370、14 249 255和11 440 684条clean reads,各样品的Q30≥98.58%。共预测出lncRNA的上下游基因310个,它们可注释到代谢进程、细胞进程、催化活性、结合和细胞等相关的35个功能条目;此外还可注释到56条代谢通路,涉及嘌呤代谢、碳代谢和丙酮酸代谢等物质代谢通路,以及甲烷代谢、氧化磷酸化和糖酵解/糖异生等能量代谢通路。上述结果表明相应的lncRNA可能通过cis作用调节上下游基因的表达,从而参与调控东方蜜蜂微孢子虫孢子中的物质和能量代谢以及细胞生命活动。进一步构建lncRNA-miRNA调控网络,分析结果显示MSTRG.3636.1、MSTRG.4498.1和MSTRG.4883.1可靶向结合4个miRNA(nce-miR-7729、nce-miR-7502、nce-miR-8639和nce-miR-8565),说明此3个lncRNA可作为ceRNA在东方蜜蜂微孢子虫孢子中发挥潜在作用。LncRNA-miRNA-mRNA调控网络分析结果显示,三者之间形成较为复杂的调控网络关系,miRNA处于调控网络的中心并联系lncRNA和mRNA,nce-miR-7502和nce-miR-8639结合的mRNA最多,达到28个;与MSTRG.3636.1、MSTRG.4498.1和MSTRG.4883.1存在靶向关系的miRNA结合多个mRNA,表明上述3个lncRNA可能通过ceRNA机制发挥作用,从而影响东方蜜蜂微孢子虫孢子中的新陈代谢和生命活动。【结论】lncRNA可能通过cis作用调控上下游基因的表达,以及作为ceRNA通过吸附miRNA间接影响靶基因的表达,从而参与调节东方蜜蜂微孢子虫孢子中的物质和能量代谢等基本生命活动。     

微小RNA及其介导的竞争性内源RNA调控网络在意大利蜜蜂工蜂中肠发育过程中的潜在作用

【目的】微小RNA（microRNA,miRNA）能够通过靶向结合导致mRNA的抑制或降解,从而对基因表达进行负调控。MiRNA在昆虫的生长、发育、免疫和细胞生命活动调控方面扮演关键角色。本研究旨在对意大利蜜蜂（Apis mellifera ligustica,简称意蜂）工蜂中肠的差异表达miRNA（differentially expressed miRNA,DEmiRNA)及其调控网络进行深入分析,并结合前期在mRNA、长链非编码RNA(long non-coding RNA,lncRNA)和环状RNA(circular RNA,circRNA)组学层面的相关研究结果,系统解析中肠发育过程的miRNA差异表达谱及竞争性内源RNA（competing endogenous RNA,ceRNA）调控网络介导的中肠发育机理。【方法】利用small RNA-seq技术对正常的意蜂7日龄和10日龄工蜂中肠（Am7和Am10）进行测序。将质控后的测序数据比对西方蜜蜂基因组,再将比对上的序列标签（tags）比对到miRBase数据库以鉴定已知miRNA。利用每百万标签序列（tags per million,TPM）算法对miRNA表达量进行计算和归一化处理。根据|log₂fold change|≥1且P≤0.05的标准筛选得到显著性DEmiRNA;利用相关软件进行靶mRNA的预测及GO和KEGG数据库注释。根据靶向结合关系及前期研究结果,利用Cytoscape软件对AMPK、PI3K-Akt、Wnt、cAMP、Hippo、mTOR、Toll/Imd、TGF-beta和MAPK 9条信号通路相关的DElncRNA/DEcircRNA-DEmiRNA-DEmRNA网络进行可视化;构建以miR-342-y为核心的ceRNA调控网络,进一步对网络中的靶mRNA进行代谢通路注释。利用茎环实时荧光定量PCR（Stem-loop RT-qPCR）验证测序数据及miRNA差异表达的可靠性。【结果】Am7 vs Am10比较组中共有112个显著性DEmiRNA,包括38个显著上调miRNA和74个显著下调miRNA,分别靶向结合7 434和9 559个mRNA,它们可注释到生物学进程相关的21和23个功能条目,细胞组分相关的16和17个功能条目,以及分子功能相关的10和11个功能条目;此外还能注释到果糖和甘露糖代谢、嘌呤代谢、甘氨酸、丝氨酸和苏氨酸代谢等物质代谢相关的83和86条通路;内吞作用、泛素蛋白水解、黑色素生成等细胞免疫相关的10和10条通路;MAPK、Jak-STAT、NF-κB等体液免疫相关的5和5条通路;以及Hippo、FoxO、Notch等涉及生长发育的13和11条信号通路。进一步构建上述9条信号通路相关的ceRNA调控网络,分析发现意蜂中肠发育过程中DEmiRNA与DElncRNA、DEcircRNA和DEmRNA之间存在复杂的调控关系;DEmiRNA位于网络中心,而DElncRNA、DEcircRNA和DEmRNA位于网络的外周。miR-342-y在意蜂工蜂中肠和幼虫肠道发育过程中均为显著性下调表达,可靶向结合3个DEcircRNA、4个DElncRNA和327个mRNA。Stem-loop RT-qPCR结果显示4个DEmiRNA的差异表达趋势与测序结果一致,说明本研究中测序数据及miRNA差异表达趋势真实可靠。【结论】DEmiRNA可能通过参与调控物质和能量代谢通路、Hippo和Wnt等信号通路以及细胞和体液免疫通路相关基因的表达影响意蜂中肠的生长和发育;miR-182-x、miR-291-y、miR-342-y、ame-miR-6001-3p等关键DEmiRNA介导的ceRNA调控网络可能在意蜂中肠发育过程发挥重要的调控作用。     

意大利蜜蜂工蜂中肠发育过程中长链非编码RNA的差异表达分析

【目的】长链非编码RNA(lncRNA)在真核生物的基因表达、表观遗传和细胞周期调控等方面发挥重要功能。本研究旨在探究意大利蜜蜂(Apis mellifera ligustica,简称意蜂)工蜂中肠发育过程中lncRNA的表达谱及其作用。【方法】利用RNA-seq技术和链特异性建库方法对意蜂7和10日龄工蜂中肠(Am7、Am10)进行深度测序,下机的原始数据经过Perl脚本过滤得到高质量有效读段。利用bowtie工具将有效读段比对核糖体数据库,进一步利用Top Hat2软件将未比对到核糖体数据库上的数据比对到参考基因组。利用CPC和CNCI软件对转录本的编码能力进行预测。通过RT-PCR对部分lncRNA进行鉴定。利用edgeR软件进行差异表达lncRNA(DElncRNA)分析,进而预测lncRNA的上下游基因,并对上下游基因进行GO及KEGG代谢通路富集分析。联用RNAhybrid、Miranda和Target Scan软件预测DElncRNA靶向结合的mi RNA及mi RNA靶向结合的靶基因,并通过Cytoscape软件构建DElncRNAs-mi RNAs-m RNAs的调控网络。最后,通过RT-qPCR验证测序数据的可靠性。【结果】Am7和Am10的深度测序分别获得134 802 058和147 051 470条原始读段,经严格过滤分别得到134 166 157和146 293 288条有效读段;共得到3 890个DElncRNA,包括2 005个上调lncRNA与1 885个下调lncRNA。RT-PCR验证结果显示共有8个lncRNA能扩增出符合预期的目的片段,表明预测出的lncRNA真实存在。DElncRNA的上下游基因数为1 793个,它们涉及42个GO条目,包括代谢进程、发育进程、细胞进程、应激反应和免疫系统进程等;这些上下游基因还涉及251条代谢通路,包括碳代谢、嘌呤代谢和脂肪酸的生物合成等物质代谢通路,硫代谢、甲烷代谢和氧化磷酸化等能量代谢通路,Hippo信号通路、Wnt信号通路和Notch信号通路等信号通路,溶酶体、内吞作用和泛素介导的蛋白水解等细胞免疫通路,以及MAPK信号通路、Jak-STAT信号通路和NF-kappa B信号通路等体液免疫通路,上述结果表明DElncRNA在意蜂中肠发育过程中参与物质和能量代谢、细胞生命活动和免疫调控。进一步分析发现TCONS_00020918可通过调控西方蜜蜂王浆主蛋白1编码基因在意蜂工蜂中肠的营养吸收、级型分化中发挥功能。DElncRNA的调控网络分析结果显示DElncRNA与mi RNA、m RNA间存在复杂的调控关系,部分DElncRNA处于调控网络的中心位置且能靶向结合较多的mi RNA,也有部分mi RNA可被多个DElncRNA共同靶向,表明这些DElncRNA可能在中肠发育中发挥重要作用。随机挑取5个DElncRNA进行RT-qPCR验证,结果显示它们的表达量变化趋势与测序结果一致,证实了本研究测序数据的可靠性。【结论】差异表达长链非编码RNA(DElncRNA)广泛参与意蜂工蜂中肠的新陈代谢、细胞活动和免疫调控并作为竞争性内源RNA(ce RNA)发挥作用,研究结果为关键lncRNA的筛选和功能研究提供了必要的数据支持。     

利用第三代纳米孔长读段测序技术构建和注释蜜蜂球囊菌的全长转录组

【目的】利用第三代纳米孔（nanopore）长读段测序技术对蜜蜂球囊菌（Ascosphaera apis,简称球囊菌）的纯化菌丝（Aam）和孢子（Aas）进行测序,构建和注释球囊菌的高质量全长转录组。【方法】通过Oxford Nanopore PromethION平台对Aam和Aas进行测序。利用Guppy软件对原始读段（raw reads）进行碱基识别（base calling）,通过过滤短片段和低质量原始读段得到有效读段（clean reads）。通过识别两端引物鉴定全长转录本序列。通过比对Nr、Swissprot、KOG、eggNOG、Pfam、GO和KEGG数据库获得全长转录本的注释信息。分别利用CPC、CNCI、CPAT、Pfam 4种方法对长链非编码RNA（long non-coding RNA,lncRNA）进行预测,取四者的交集作为高可信度的lncRNA。【结果】Aam和Aas的纳米孔测序分别测得6 321 704和6 259 727条原始读段,经质控得到5 669 436和6 233 159条有效读段,其中包含的全长有效读段分别为4 497 102（79.32%）和4 963 101（79.62%）条。共鉴定到9 859和16 795条非冗余全长转录本,N50分别为1 482和1 658 bp,平均长度分别为1 187和1 303 bp,最大长度分别为6 472和6 815 bp。Venn分析结果显示有6 512条非冗余全长转录本为菌丝和孢子所共有,分别有3 347和10 283个非冗余全长转录本为二者特有。此外,在球囊菌菌丝和孢子中共鉴定到20 142条全长转录本,其中分别有20 809、11 151、17 723、12 164、11 340和9 833条全长转录本可注释到Nr、KOG、eggNOG、Pfam、GO和KEGG数据库。注释全长转录本数量最多的物种是球囊菌、Polytolypa hystricis和荚膜组织胞浆菌（Histoplasma capsulatum）。GO数据库注释结果显示,上述全长转录本可注释到45个功能条目,涉及细胞组件、细胞和细胞器等细胞组分相关条目;催化活性、结合和转运器活性等分子功能相关条目;以及细胞进程、代谢进程和单一组织进程等生物学进程相关条目。KEGG数据库注释结果显示,上述全长转录本还可注释到抗生素的生物合成、核糖体、氨基酸的生物合成、碳代谢和剪接体等49条通路。此外,鉴定到648条高可信度的lncRNA,包含480条基因间区lncRNA、119条反义链lncRNA和49条正义链lncRNA。【结论】构建和注释了球囊菌的首个高质量全长转录组,为探究球囊菌转录组的复杂性,完善参考基因组的序列和功能注释信息以及深入开展球囊菌可变剪接体的功能研究提供了关键依据。     

不同温度下TMV侵染枯斑三生烟的LncRNA差异表达

【目的】 筛选不同温度下烟草花叶病毒（Tobacco mosaic virus,TMV）侵染后枯斑三生烟（Nicotiana tabacum var. Samsun NN）差异表达的长链非编码RNA（long non-coding RNA,lncRNA）,研究lncRNA在枯斑三生烟抗性反应中的作用。【方法】 N基因的温度敏感性使枯斑三生烟在25℃时具备对TMV的抗性、在31℃抗性丧失,在这两个温度条件下对枯斑三生烟接种TMV和磷酸盐缓冲盐水（phosphate buffered saline,PBS）,48 h后提取系统叶总RNA,构建链特异性文库后进行深度测序。对测序结果进行过滤后利用HTSeq将有效数据与近缘品种TN90（N. tabacum var. TN90）基因组比对,筛选得到lncRNA后利用FPKM法估计lncRNA的表达水平。通过edgeR筛选差异表达lncRNA（differentially expressed lncRNA,DElncRNA）,并利用qRT-PCR技术对这一结果进行验证。通过共定位及共表达分析预测DElncRNA的靶基因,通过参考基因组注释、GO和KEGG富集分析研究靶基因的功能。【结果】 4个处理共12个样本经lncRNA-seq各测得约8 000万条clean reads,共获得4 737条已知lncRNA、40 169条新lncRNA。其中64个lncRNA在不同温度条件下TMV侵染后存在差异表达,qRT-PCR测定结果显示这些lncRNA的测序正确率在80%左右,表明本研究所得测序数据具备较高的可信度。对DElncRNA进行靶基因预测,发现一些基因同时被25℃下调和31℃上调的DElncRNA靶向。靶基因注释功能丰富,主要参与植物抗病、激素和代谢等生理过程。部分可能与激素通路相关的lncRNA,在25℃下TMV侵染时呈现下调趋势,而在31℃下TMV侵染则呈现上调趋势。GO富集分析显示靶基因主要参与构成膜、囊泡等组分,具备钙、钾离子通道抑制剂活性等分子功能,使相应离子得以转运引发随后的反应,同时也参与发病、抗原加工和呈现、细胞分裂素代谢等生理过程。KEGG分析发现靶基因显著富集在植物激素信号转导通路,25℃下调和31℃上调的DElncRNA靶基因同时富集在激素信号传导、ABC运输蛋白、苯丙烷类生物合成等通路。【结论】 不同温度（25℃和31℃）条件下TMV侵染枯斑三生烟后,长链非编码RNA差异表达,DElncRNA通过作用于激素信号传导、物质转运等过程参与寄主系统获得性抗性反应。研究结果可为揭示植物系统获得性抗性中lncRNA的调控功能以及新型抗病毒技术开发提供依据。     

捻转血矛线虫伊维菌素耐药相关长链非编码RNA及其调控功能分析

为探究捻转血矛线虫伊维菌素耐药分子标记和关键调控元件,通过高质量的长链非编码RNA(Long non-coding RNA,lncRNA)表达谱对捻转血矛线虫敏感与耐药虫株lncRNA的顺式调控(Cis-regulation, cis)和竞争性内源RNA(Competing endogenous RNA,ceRNA)调控进行分析,筛选影响虫体伊维菌素耐药机制的微小调控分子。本研究利用RNA-seq技术对捻转血矛线虫伊维菌素敏感虫株和耐药虫株进行测序,根据生物学软件预测结果,筛选出位于蛋白编码基因上下游10 kb以内的lncRNA作为顺式调控作用的元件。对差异lncRNA进行GO和KEGG富集分析,并将显著差异的lncRNA和miRNA通过Target finder和Cytoscape v3.7.1软件建立lncRNA-miRNA调控网络和可视化分析。根据miRNA调控元件对lncRNA-circRNA-miRNA-mRNA调控网络进行连通性分析,结合耐药相关信号通路的富集结果构建lncRNA-miRNA-mRNA可视化网络,同时对随机选取的10个差异lncRNA进行qRT-PCR验证。...     

意大利蜜蜂幼虫肠道在球囊菌侵染前期的 差异表达microRNA及其调控网络

【目的】 微小RNA（microRNA,miRNA）是一类重要的基因表达调控因子,可影响宿主与病原间的互作过程。蜜蜂球囊菌（Ascosphaera apis）是一种特异性侵染蜜蜂幼虫的致死性真菌病原。本研究旨在对意大利蜜蜂（Apis mellifera ligustica,简称意蜂）幼虫肠道在球囊菌胁迫前期的差异表达miRNA（differentially expressed miRNA,DEmiRNA）及其靶基因进行深入分析,在miRNA组学水平探究意蜂幼虫在球囊菌侵染前期的胁迫应答,并通过构建显著DEmiRNA的调控网络筛选出与宿主应答相关的关键miRNA。【方法】 利用small RNA-seq（sRNA-seq）技术对正常及球囊菌侵染的意蜂4日龄幼虫肠道（AmCK和AmT）进行高通量测序,首先对原始数据进行质控和评估,随后将过滤后的数据与西方蜜蜂（Apis mellifera）参考基因组进行比对;将比对上的序列标签（tags）注释到miRBase数据库,得出已知miRNA的表达量;通过TPM（tags per million）算法对各样本中miRNA的表达量进行归一化处理,以|log₂ fold change|≥1且P≤0.05作为标准筛选得到显著DEmiRNA;利用TargetFinder 软件预测显著DEmiRNA的靶基因,并对其进行GO和KEGG代谢通路（pathway）富集分析。利用Cytoscape软件对miRNA-mRNA调控网络进行可视化。最后,利用茎环反转录PCR（Stem-loop RT-PCR）和荧光定量PCR（qPCR）验证测序数据的可靠性。【结果】 AmCK和AmT样品的测序分别得到13 553 302和10 777 534条原始读段（raw reads）,经严格过滤后得到的有效读段（clean reads）数分别为13 186 921和10 480 913条。各样品的生物学重复间的Pearson相关性系数分别在0.9822和0.9508以上。共有10个显著DEmiRNA,包括4个上调miRNA和6个下调miRNA。显著DEmiRNA在AmT的整体表达水平低于AmCK。10个显著DEmiRNA可靶向结合3 788个靶基因,其中上调miRNA的1 240个靶基因可注释到GO数据库中的39个GO条目,主要富集在结合、细胞进程、代谢进程和应激反应等;下调miRNA的749个靶基因可注释到34个GO条目,主要富集在细胞进程、结合、代谢进程和应激反应等。KEGG数据库注释结果显示,上调miRNA和下调miRNA的靶基因分别注释到95和66条代谢通路,富集基因数最多的分别是Wnt信号通路、Hippo信号通路、光传导以及内吞作用、磷脂酰肌醇信号系统、嘌呤代谢。对于上调和下调miRNA,分别有31和52个靶基因注释到内吞作用,15和7个靶基因注释到泛素介导的蛋白水解,11和5个靶基因注释到Jak-STAT信号通路,1和3个靶基因注释到MAPK信号通路。显著DEmiRNA与靶mRNA之间形成复杂的调控网络,7个显著DEmiRNA靶向结合96个与Wnt信号通路相关的mRNA,8个显著DEmiRNA靶向结合55个与内吞作用相关的mRNA。Stem-loop RT-PCR和qPCR结果验证了测序数据的可靠性。【结论】 对意蜂幼虫肠道在球囊菌侵染前期的DEmiRNA及其靶基因进行预测和分析,并构建和分析了DEmiRNA-mRNA调控网络,研究结果提供了宿主miRNA的表达谱和差异表达信息,揭示了DEmiRNA通过调控细胞生命活动、新陈代谢以及部分细胞和体液免疫等生物学过程参与宿主的胁迫应答。miR-4331-y、miR-4968-y、miR-8440-y、novel-m0023-5p和novel-m0025-3p共同参与了宿主的Wnt信号通路和内吞作用的调控,可作为白垩病治疗的潜在分子靶标。     

苏博美利奴羊毛囊发育相关lncRNA与mRNA共表达网络的构建

【目的】 随着高通量测序技术的发展和完善,海量的转录组测序数据随之涌现,基因网络的方法也越来越多被用于研究中;细毛羊毛囊发育受多个基因及lncRNA共同调控,单独研究某一分子并不足以发现其调控机制,本研究旨在构建毛囊发育相关lncRNA和mRNA的共表达网络,挖掘细毛羊毛囊发育的潜在候选基因。【方法】 以苏博美利奴羊胎龄65 d（D65）、85 d（D85）、105 d（D105）、135 d（D135）的胎儿以及出生后7 d（A7）、30 d（A30）的羔羊肩胛部皮肤组织,每个时期有3个生物学重复,共18个样本,进行转录组测序以获得6个不同发育时期的lncRNA与mRNA表达谱数据,筛选出相邻时期的差异表达lncRNA与mRNA,运用加权基因共表达网络分析（weighted gene co-expression network analysis,WGCNA）方法构建共表达模块,使用DAVID在线工具进行GO（gene ontology）和 KEGG pathway富集分析以找到与毛囊发育相关的模块,最后从目标模块筛选高连通度的lncRNA与mRNA使用Cytoscape软件进行网络可视化。【结果】 从表达谱数据中筛选得到9 070个差异表达的lncRNA与mRNA,运用WGCNA方法得到11个模块,使用DAVID富集分析发现honeydew1、paleturquoise、skyblue2模块中的基因富集在皮肤发育、毛囊发育、毛囊形态发生、Wnt信号通路负调控、细胞粘附等生物过程,以及Wnt信号通路、TGF-β信号通路、Hedgehog信号通路、紧密连接、MAPK信号通路、ECM-受体相互作用等毛囊发育相关的信号通路,并筛选这3个模块中连通度高的lncRNA和mRNA构建了子网络,得到包括TGFB2、CTSB、SFN、SPINT1、FAM83G、GSDMA、MPZL3、VIM、CRABP1在内的毛囊发育相关基因,预测出ENSOART00000029117、TCONS_00489976、TCONS_00376759等24个lncRNA的可能靶基因。【结论】 首次运用WGCNA方法构建了苏博美利奴羊毛囊发育相关lncRNA与mRNA的共表达网络,鉴别出3个毛囊发育相关的共表达模块,发现多个与毛囊发育相关的潜在候选基因,并预测出24个lncRNA可能的靶基因。     

红橘响应褐斑病菌侵染的转录组学分析

【目的】 明确红橘（Citrus reticulata Blanco, cv. Hongjv）接种褐斑病致病菌——链格孢菌橘致病型（Alternaria alternata tangerine pathotype）后基因种类和表达量在转录水平的变化规律,确定红橘响应该致病型侵染的关键基因。【方法】 采用链格孢菌橘致病型接种红橘离体叶片,28 h后选取感病叶片和未接种叶片提取RNA,进行转录组高通量测序,然后利用生物信息学分析,以甜橙基因组为参考,以|log₂ fold change|≥1,q-value≤0.01为阈值选取感病和健康红橘叶片转录组的差异表达基因,应用GO数据库对差异表达基因（differentially expressed gene,DEG）进行功能分类,KEGG分析代谢途径,MapMan软件分析生物胁迫信号通路相关基因的表达变化。采用qRT-PCR方法对测序结果进行验证。【结果】 红橘接种链格孢菌橘致病型28 h后产生大量与胁迫相关的差异表达基因,获得上调差异基因5 173个,下调差异基因6 555个。GO功能分类显示差异基因主要与蛋白结合、膜、氧化还原过程等相关。通过KEGG富集和MapMan软件分析发现,红橘在受链格孢菌橘致病型胁迫的过程中基础代谢被严重破坏。乙烯、水杨酸和生长素等寄主防御相关的植物激素信号转导途径多个基因表达出现差异,其中乙烯起主导作用,乙烯受体ETR 3个成员被不同程度激活,下游激酶和乙烯响应因子均上调,生长素大部分关键信号基因、绝大部分生长素响应因子ARF和水杨酸合成途径的基因均下调表达。同时,黄酮醇、花青素、萜类化合物和生物碱合成相关基因受该菌诱导显著变化,萜类合成中大部分基因下调,而黄酮类合成相关上调基因数量和表达趋势均强于表达下调基因,有抗虫和抑菌作用的硫代葡萄糖苷基因呈上调趋势。进一步研究发现,大量参与抗逆过程的转录因子如WRKY、bZIP、ERF、MYB、NAC被诱导激活,其中大部分WRKY和bZIP转录因子受该菌正向调控,超过50%的ERF家族基因表达上调;在转录因子调控下,PTI及ETI响应基因如受体激酶、R蛋白、NBS抗病蛋白等大量表达,多个PR家族抗菌蛋白基因上调表达,22个抗氧化保护酶系统POD成员基因受到活性氧信号激发大量表达。以上结果表明链格孢菌橘致病型侵染对寄主内部生理状态产生显著影响。选取了19个与植物抗病相关基因进行qRT-PCR分析,其基因表达趋势与测序数据一致。【结论】 获得了红橘响应链格孢菌橘致病型侵染的差异表达基因及显著上调表达基因,其主要富集于代谢过程、应激反应及转录调控等条目中,这些基因的相互协同调控是红橘对该致病型产生防御反应的重要机制。