蜜蜂繁殖性状编码RNA和非编码RNA的研究进展

蜜蜂蜂王的繁殖性状是蜂群产生经济效益的基础性状，繁殖性状的分子机制研究是蜜蜂遗传育种领域的重要课题。在蜜蜂的重要繁殖活动产卵和卵巢激活过程中，编码RNA和非编码RNA发挥了重要作用。介绍了蜜蜂繁殖性状相关的信使RNA、小RNA、长链非编码RNA、环状RNA的基本特征和功能，综述了编码RNA和非编码RNA对蜜蜂卵巢激活和产卵过程的调控机制。蜜蜂繁殖性状RNA的研究对于阐明蜜蜂繁殖性状的分子遗传机制意义重大，也为挖掘繁殖性状分子标记基因、开展分子育种工作，进而提高蜜蜂繁殖性能、增加养蜂经济效益奠定基础。     

蜜蜂球囊菌中长链非编码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及其靶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数据库注释。根据靶向...     

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

【目的】微小RNA(micro RNA,miRNA)能够通过靶向结合导致m RNA的抑制或降解,从而对基因表达进行负调控。Mi RNA在昆虫的生长、发育、免疫和细胞生命活动调控方面扮演关键角色。本研究旨在对意大利蜜蜂(Apis mellifera ligustica,简称意蜂)工蜂中肠的差异表达miRNA(differentially expressed miRNA,DEmiRNA)及其调控网络进行深入分析,并结合前期在m RNA、长链非编码RNA(long non-coding RNA,lncRNA)和环状RNA(circular RNA,circ RNA)组学层面的相关研究结果,系统解析中肠发育过程的miRNA差异表达谱及竞争性内源RNA(competing endogenous RNA,ceRNA)调控网络介导的中肠发育机理。【方法】利用small RNA-seq技术对正常的意蜂7日龄和10日龄工蜂中肠(Am7和Am10)进行测序。将质控后的测序数据比对西方蜜蜂基因组,再将比对上的序列标签(tags)比对到miRBase数据库以鉴定已知miRNA。利用每百万标签序列(tags p...     

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

[目的]利用第三代纳米孔(nanopore)长读段测序技术对蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)的纯化菌丝(Aam)和孢子(Aas)进行测序,构建和注释球囊菌的高质量全长转录组.[方法]通过Oxford Nanopore PromethION平台对Aam和Aas进行测序.利用Guppy软件对原始读...     

球囊菌侵染对意大利蜜蜂成年工蜂的影响

为了探究球囊菌侵染对意大利蜜蜂成年工蜂存活率、免疫基因、发育基因及肠道菌群的影响，将纯化培养的球囊菌孢子接种6日龄意大利蜜蜂成年工蜂，利用荧光定量PCR(qRT-PCR)检测球囊菌侵染对意大利蜜蜂成年工蜂免疫基因、发育基因及肠道菌群的影响。结果表明，连续3 d接种球囊菌孢子对意大利蜜蜂成年工蜂存活率无显著影响(χ²=0.133,P=0.715)。意大利蜜蜂工蜂中肠与后肠球囊菌孢子含量在接种球囊菌孢子4 d时分别显著高于接种8 d时(P<0.05)。与未接种组相比，接种球囊菌孢子4 d，意大利蜜蜂成年工蜂免疫相关基因abaecin、apidaecin、defensin-1、defensin-2、hymenoptaecin、GNBP-1和IMD均显著上调表达(P<0.05)；接种球囊菌孢子8 d，意大利蜜蜂成年工蜂免疫相关基因相对表达量无显著差异(P>0.05)；接种球囊菌孢子4 d和8 d，发育相关基因VGMC、usp、kr-h1、AMHex10869相对表达量无显著差异(P>0.05)。与未接种组相比，接种球囊菌孢子4 d，意大利蜜蜂成年工...     

蜜蜂患白垩病虫体内一株球囊菌拮抗细菌的分离与鉴定

【目的】寻找抑制蜜蜂球囊菌的拮抗菌,用于白垩病的生物防治。【方法】利用细菌纯化技术从患白垩病的蜜蜂幼虫体内分离纯化了一株对蜜蜂球囊菌(Ascosphaera apis)有抑制作用的细菌,并结合形态学、革兰氏染色以及16S rDNA序列分析技术进行鉴定,同时通过蜜蜂体内、体外接种试验研究其对蜜蜂球囊菌的抑制作用。【结果】该细菌被初步鉴定为蜡样芽胞杆菌(Bacillus cereus),它在培养基中能够完全抑制蜜蜂球囊菌的生长,但其发酵液对蜜蜂球囊菌的生长没有影响。将含菌饲料喂蜜蜂幼虫,蜜蜂幼虫的死亡率和生长发育没有受到影响,说明该细菌对蜜蜂幼虫没有致病性。蜜蜂幼虫体内同时接种细菌和球囊菌孢子,发现细菌对球囊菌在蜜蜂体内的萌发和初期菌丝的生长没有影响,蜜蜂依然能够患病,但生长后期该细菌能够降解球囊菌菌丝,对球囊菌子代孢子的产生表现出一定的抑制作用。【结论】 该研究结果为蜜蜂白垩病的生物防治提供了理论依据。     

RNA干扰技术在害虫防治中的应用及其安全性

RNA干扰（RNA interference, RNAi）是一种在生物体内高度保守的、由20～30个核苷酸（nucleotide, nt）组成的小分子非编码RNA引发的基因沉默现象。在昆虫等真核生物中,干扰小RNA(small interfering RNA, siRNA)、微RNA(microRNA, miRNA)和PIWI蛋白互作RNA(PIWI-interacting RNA, piRNA)等小分子非编码RNA在RNA干扰中发挥着重要作用。RNA干扰在害虫防治中的应用主要体现在基因功能研究、转基因抗虫植物及新型核酸农药研发等方面。RNA干扰技术及其产品在应用之前,需严格评价其在非靶标效应、靶标害虫的抗性及环境持久性方面的安全风险。本文就RNA干扰技术的作用机制、RNA干扰在农林害虫防治中的应用及其安全性等方面进行了综述,以期为RNA干扰技术在害虫防治中获得更广泛的应用提供一定的理论依据。     

绵羊繁殖相关编码RNA和非编码RNA的研究进展

随着分子生物学和生物信息学的发展及应用，已有诸多研究从影响绵羊卵巢发育、排卵、产羔、睾丸发育和精子发生等主要繁殖活动相关编码RNA和非编码RNA方面着手，对绵羊繁殖相关的编码和非编码RNA进行筛选和鉴定。本文对影响绵羊繁殖的编码RNA和非编码RNA研究成果进行总结梳理，共整理了61个关键mRNAs、9个miRNAs、17个lncRNAs和4个circRNAs,并进一步阐述RNA的发展历史和生物学功能，以及编码RNA与非编码RNA构成调控绵羊繁殖的复杂ceRNA网络。结果表明：1)不同年龄公羊睾丸和母羊卵巢在繁殖过程中mRNA呈动态变化；2)非编码RNA与DNA、RNA或蛋白质之间相互作用，调控转录、翻译以及翻译后修饰等过程，实现不同品种、年龄和饲喂水平绵羊对精子发生和卵泡发育等过程的调控；3)编码RNA和非编码RNA形成稳定的ceRNA调控机制使繁殖活动趋于动态平衡。本综述为挖掘绵羊繁殖性状分子标记基因、开展分子育种工作、改善绵羊繁殖性能及提高羊产业经济效益打下基础。     

长链非编码RNA分子作用机制研究进展

长链非编码RNA(Long non-coding RNA,lncRNA)是缺乏开放阅读框架,且长度大于200 nt不编码蛋白质的转录本,其绝对数量大、种类多,在表达模式上具有明显的细胞特异性.lncRNA通过碱基配对与DNA或RNA,或通过RNA高级结构与蛋白质结合,发挥多种生物学功能,共同构成了一个复杂而精细的分子调...     

东方蜜蜂微孢子虫孢子中长链非编码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(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的筛选和功能研究提供了必要的数据支持。     

不同温度下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和KE...     

意大利蜜蜂幼虫肠道在球囊菌侵染前期的 差异表达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,进行转录组高通量测序,然后利用生物信息学分析,以甜橙基因组为参考,以|log2 fold change|≥1,q-value≤0.01为阈值选取感病和健康红橘叶片转录组的差异表达基因,应用GO数据库对差异表达基因(differentially expressed gene,DEG)进行功能分类,KEGG分析代谢途径,MapMan软件分析生物胁迫信号通路相关基因的表达变化。采用qRT-PCR方法对测序结果进行验证。【结果】红橘接种链格孢菌橘致病型28 h后产生大量与胁迫相关的差异表达基因,获得上调差异基因5 173个,下调差异基因6 555个。GO功能分类显示差异基因主要与蛋白结合、膜、氧化还原过程等相关。通过KE...