蜜蜂球囊菌中长链非编码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间接影响靶基因的表达,从而参与调节东方蜜蜂微孢子虫孢子中的物质和能量代谢等基本生命活动。     

苏博美利奴羊毛囊发育相关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可能的靶基因。     

东方蜜蜂微孢子虫孢子中长链非编码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间接影响靶基因的表达,从而参与调节东方蜜蜂微孢子虫孢子中的物质和能量代谢等基本生命活动。     

基于芯片筛选螺原体感染中华绒螯蟹血细胞的免疫microRNAs

为筛选分析螺原体(Spiroplasma eriocheiris)感染中华绒螯蟹(Eriocheir sinensis)血细胞免疫相关microRNAs(miRNA)及其差异表达情况,利用微流体芯片技术,对正常和感染螺原体的中华绒螯蟹血细胞miRNA进行系统的鉴定与表达检测;结合生物信息学方法,分析其靶基因和免疫信号通路。利用前期研究和已有报道897条成熟miRNA作为探针订制单色芯片,使用Cy3染料标记,通过与中华绒螯蟹血细胞总RNA样本杂交,鉴定miRNA,检测其表达情况;对差异表达的miRNA进行靶基因预测、功能聚类和通路分析;随机挑选部分差异表达的miRNA进行茎环RT-qPCR验证。结果显示,总共鉴定出中华绒螯蟹血细胞miRNA共303个;健康组和试验组共表达miRNA有48个;差异表达miRNA为27个,其靶基因共注释到116个GO和63个KEGG通路;通过茎环RT-qPCR检测,结果与芯片基本一致。研究发现,差异表达miRNA靶基因,许多均存在于重要免疫信号通路之中,研究结果可为进一步分析中华绒螯蟹血细胞螺原体感染中的miRNA与其靶基因间的免疫调控机制提供借鉴。     

杨树miRNA的靶基因预测及低氮胁迫表达分析

目的鉴定杨树受低氮胁迫后miRNA的靶基因,分析靶基因在氮胁迫后的差异表达并探讨其功能,为揭示杨树低氮胁迫下miRNA的调控功能提供参考,并为树木低氮营养高效利用育种提供重要的候选基因。方法根据miRNA的保守性及与靶基因的严谨互补配对关系,以杨树miRNA为探针利用靶基因预测软件psRNATarget,通过与毛白杨转录组的基因序列进行比对鉴定靶基因,进一步开展毛白杨受低氮胁迫后靶基因的差异表达分析及功能注释。结果获得了131个miRNA家族的242个miRNA成员对应的3 024个靶基因,分别参与了植物激素信号转导、次生代谢产物的生物合成、氨基酸合成代谢、碳代谢和RNA运输等通路。57个靶基因在低氮胁迫处理后发生显著变化,其中受到诱导（29个）和抑制（28个）的基因数目相当。14个低氮胁迫响应的miRNA,其对应的11个靶基因也发生了显著的差异表达变化,其中miRNA和靶基因表达量发生相反变化的有8个miRNA。本研究发现参与植物激素信号转导的靶基因（2个）及参与代谢途径的靶基因（6个）发生了差异表达。miR162的靶基因编码ABC转运蛋白,miR393运用于靶基因KAT2调节Na+和K+动态平衡,miR399的靶基因PIF3编码光敏色素互作因子PIFs蛋白,这些miRNA及靶基因可能在杨树响应低氮胁迫中发挥重要作用。结论本文鉴定到了毛白杨中一批低氮胁迫响应miRNA的靶基因,可调控杨树对氮逆境胁迫信号的反应。这些miRNA及靶基因为进一步揭示miRNA及靶基因在低氮胁迫下的调控功能提供了研究线索,为树木氮营养的高效利用改良提供了重要候选基因。      

捻转血矛线虫伊维菌素耐药相关长链非编码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靶基因的预测与验证技术研究进展

MicroRNA(miRNA)是真核生物中具有重要调控功能的一类小分子RNA。植物miRNA通过碱基互补配对方式引导RNA沉默复合体识别靶基因并介导靶基因mRNA被剪切或蛋白质翻译被抑制,所以对植物miRNA靶基因的分析和验证是研究miRNA功能必不可少的一部分。植物miRNA靶基因可通过生物信息学方法预测,主要根据miRNA与mRNA的互补程度进行,然而理论上预测的靶基因尚有待于进一步的试验验证。简要介绍了植物miRNA产生过程的最新研究进展以及靶基因的预测方法,重点综述了植物miRNA靶基因的验证方法,包括RLM-5'RACE、降解组测序、瞬时共表达和转基因系统,为更深入地研究植物miRNA介导的调控机制提供参考。     

苏尼特羊不同生长时期肌肉组织lncRNA的差异表达分析

为鉴定和分析苏尼特羊肌肉不同发育时期lncRNA的差异表达,利用RNA-seq技术对120d苏尼特羊胎儿和5月龄羔羊肌肉组织进行了lncRNA比较分析,筛选出差异表达的lncRNA及其靶基因,并对预测的mRNA进行了GO和KEGG分析,同时对随机筛选的差异表达lncRNA进行qRT-PCR验证。结果显示,120d和5月龄羔羊肌肉之间差异表达的lncRNA数共为626个,共靶向1 273个mRNA编码基因。GO和KEGG分析表明,这些差异基因主要涉及代谢过程的调控,生物合成,基因表达,蛋白结合以及mTOR信号通路、FoxO信号通路、胰岛素信号、肌动蛋白细胞骨架的调节、Wnt信号通路和间隙连接等信号通路。qRT-PCR验证证明,随机选取的8个lncRNA在肌肉组织中的表达量与RNA-seq结果一致。综上,本研究利用RNA-Seq技术分析了苏尼特羊120d胎儿和5月龄羔羊的差异表达lncRNA及其靶基因,发现其参与了不同生长时期苏尼特羊的肌肉发育和生长过程,为从lncRNA的角度更好地理解苏尼特羊肌肉生长发育的遗传调控机制提供了理论基础,也为绵羊分子辅助育种提供参考。     

亚洲棉短纤维发育相关长链非编码RNA的鉴定及表达

【目的】长链非编码RNA(long non-coding RNAs,lncRNAs)是一类无蛋白质编码能力，但参与许多重要生命活动调控过程的长度大于200 nt的RNA。通过对亚洲棉无短纤维突变体（GA0149）和野生型（GA0146）纤维发育早期的转录组数据进行分析，挖掘调控短纤维发育的lncRNA，并明确其调控网络，为进一步解析棉花纤维发育机制奠定基础。【方法】选择GA0146和GA0149 2个材料在开花后当天（0 DPA）及花后3 d(3 DPA)、5 d(5 DPA)和8 d(8 DPA)的胚珠和纤维为材料进行转录组测序。鉴定lncRNA并预测其调控的靶基因；通过mRNA和lncRNA的差异表达分析，比较2个材料在不同纤维发育时期的差异。进一步利用KOBAS软件预测对差异lncRNA的靶基因进行富集分析并预测其参与的生物过程；最后通过实时荧光定量（RT-qPCR）技术对25个差异表达的lncRNA转录组数据进行验证。【结果】共鉴定获得15 339个lncRNA，其中11 595个lncRNA位于基因间区，包括2 428个反义lncRNA、350个内含子lncRNA及966个正...     

深度测序鉴定玉米雄穗花器官分化期响应干旱胁迫的miRNA和其靶基因

【目的】 运用现代分子生物学及高通量测序技术,研究不同品种玉米雄穗花器官分化期响应干旱胁迫的差异miRNA和其靶基因,挖掘和鉴定与玉米发育相关的基因,分析其性状差异的信号通路及分子调控网络。【方法】 以耐旱自交系"PHBA6"和干旱敏感自交系“吉63”为研究对象,应用深度测序技术进行miRNA文库构建,鉴定其差异表达的 miRNA,对差异表达mi RNA及其靶基因进行功能注释、聚类分析及通路富集。【结果】 鉴定了总共337种前体miRNA,其中包含289种已知的miRNA和48种新的miRNA。在3个文库,两个分组中共有155种差异表达miRNA。基于GO功能分类及遗传和基因组(KEGG)的功能富集显示,这些miRNA可能通过靶向一系列与胁迫相关的基因而在干旱胁迫中发挥作用。至少55个预测的靶基因进一步被60个miRNA调控。NAC,MYB和MAPK基因家族在干旱胁迫下评分最高,在植物抗旱中重要作用。一系列miRNA与这些排名靠前的基因相关,包括miR164、miR172、miR1520、miR6158、ghr-n24、ghr-n56等。【结论】 miRNAs可能在玉米雄穗花器官分化期耐旱中发挥重要作用,miRNAs的筛选为分子辅助育种和转基因育种将提供新的靶标。     

两品种鸭下丘脑miRNA的鉴定及其表达验证分析

为了鉴定与鸭性成熟启动相关的miRNA,根据鸭卵巢组织的发育情况,利用Solexa技术分别对105日龄连城白鸭和樱桃谷鸭的下丘脑组织进行小RNA测序分析。两个小RNA文库共检测到939个成熟miRNA,共有miRNA数量806个。以樱桃谷鸭为对照,两个文库之间总共检测到290个差异表达miRNA,其中上调miRNA 85个,下调miRNA 205个。靶基因预测分析结果表明,较多的靶基因被注释到神经活性配体-受体相互作用、黏着斑、MAPK信号通路、胞吞作用等通路中。定量PCR结果表明,2个差异表达miRNA靶基因的mRNA表达量在两品种鸭下丘脑组织中存在着极显著差异。     

1月龄五指山猪与长白猪骨骼肌miRNA转录组比较

为研究microRNA(miRNA)在五指山猪和长白猪肌肉生长发育中的差异和探究差异miRNA对骨骼肌发育转录后调控的影响,以1月龄的五指山猪和长白猪背最长肌为材料,通过转录组测序和生物信息学分析筛选与骨骼肌发育相关的差异miRNA。结果显示,从2个文库中共鉴定获得318种己知的猪miRNA,五指山猪和长白猪分别鉴定出312个和306个已知miRNA,同时分别发现了83个和88个新的miRNA。获得17个差异显著的miRNA,其中16个miRNA表达上调,1个miRNA表达下调。靶基因预测和生物学功能预测发现17个差异miRNA,共预测到535个靶基因,靶向作用于mTOR信号通路和肌动蛋白细胞骨架调控的ssc-miR-486,ssc-miR-204,ssc-miR-338和ssc-miR-885-3p 4个miRNA可能参与骨骼肌生长发育过程。     

黔北黑猪骨骼肌中肌纤维类型转化相关circRNAs的筛选

选取9头6月龄健康黔北黑猪,随机分成3组,采集背最长肌(LDM)和腰大肌(PMM)肌肉组织,通过RNA–seq构建环状RNA(circRNAs)的表达谱,筛选差异表达的circRNAs,在黔北黑猪2种肌肉中共鉴定出63个差异表达circRNAs,其中31个在LDM中表达上调,32个表达下调;对差异表达circRNAs的源基因进行GO、KEGG富集分析,结果表明,差异表达的circRNAs的源基因主要富集在磷酸酶活性的调控、肌节组织活动、肌球蛋白–轻链–磷酸酶活性的正调控等条目;富集的KEGG通路中与肌纤维类型转化相关的通路包括肌动蛋白细胞骨架的调控、ECM–受体互作、PI3K–Akt信号通路。预测差异表达circRNA的靶miRNAs,并筛选对应的靶mRNAs,构建的circRNA–miRNA–mRNA调控网络显示,circ_0003749–miR–27a–MSTN、circ_0003749–miR–27a–PPARγ、circ_0012316–miR–23a–Myh1/2/4、circ_0007093–miR–214–3p–MEF2C、circ_0010118–miR–132–FoxO1和circ_0010118 –miR–374a–5p–FoxO1可能通过ceRNA网络调控肌纤维类型的转化。     

低致病性A/Anhui/1/2013(H_7N_9)流感病毒感染小鼠肺组织的miRNA表达谱分析

[目的]本试验旨在利用基因芯片测序技术分析宿主miRNA在调控流感病毒感染小鼠肺组织过程中的作用。[方法]应用低致病性H_7N_9亚型A/Anhui/1/2013流感病毒及PBS感染BALB/c小鼠,感染后3 d取小鼠肺组织分别利用转录组测序技术和miRNA芯片测序技术筛选差异表达mRNA和miRNA,并利用RT-qPCR验证差异miRNA。之后分别采用Targetscan、PITA及microRNAorg软件预测靶基因并结合转录组mRNA信息筛选候选miRNA的假定靶基因,最后利用Gene Ontology(GO)和Kyoto Encyclopedia of Gene and Genomes(KEGG)分析差异miRNA的生物学功能及其可能调控的信号通路。[结果]与PBS组相比,低致病性H_7N_9病毒感染组共筛选到265个差异表达miRNA,其中143个miRNA显著上调,122个miRNA显著下调。差异表达miRNA经RT-qPCR验证,10条候选miRNA的RT-qPCR结果与芯片结果有非常好的一致性。进一步KEGG分析表明,这些差异表达miRNA主要富集在Rap1、PI3K-Akt、Hippo、MAPK、Wnt、黏着斑、自噬等免疫相关信号通路中。结合差异mRNA信息进行miRNA-mRNA调控网络分析,显示主要有15条差异表达miRNA和31个差异靶基因富集到这些信号通路中。[结论]成功筛选到了低致病性H_7N_9感染小鼠肺组织后引起的差异表达miRNA,且RT-qPCR证明其与基因芯片测序结果表达趋势具有很好的一致性,为深入研究宿主miRNA在调控低致病性H_7N_9流感病毒与宿主相互作用的分子机制奠定了基础。     

不同发育时期亚麻茎秆中木质素积累相关miRNA及其靶基因的挖掘分析

[目的]深入挖掘分析亚麻茎秆不同发育时期木质素积累相关的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负调控,并参与亚麻茎秆中木质素的积累.     

梅山猪和杜洛克猪卵泡中差异表达lncRNA克隆鉴定及与miRNAs相关性分析

【目的】通过克隆鉴定梅山猪和杜洛克猪卵泡期第4天M2卵泡中差异表达的lncRNA-ALDBSSCT0000005583,分析其在猪颗粒细胞与miRNAs表达量的相关性,为探究lncRNA调控miRNA在母猪卵泡发育过程中的作用提供理论依据。【方法】对课题组前期在梅山和杜洛克M2卵泡中筛选出的差异表达lncRNA-ALDBSSCT0000005583,进行RT-qPCR验证,并利用RACE克隆其全长序列;根据编码潜力评估工具CPAT、CPC对该lncRNA的编码能力进行预测,原核表达试验进一步鉴定其编码能力;核质分离试验对lncRNA-ALDBSSCT0000005583进行亚细胞定位,RT-qPCR检测其在各组织中的表达水平;利用miRBase网站查找猪的miRNA数据库,结合RNAhybrid、miRanda在线软件预测与lncRNA-ALDBSSCT0000005583有相互作用的物种间保守miRNA,使用TargetScan和miRanda预测与lncRNA-ALDBSSCT0000005583具有相互作用miRNA的靶基因,并对其靶基因进行GO富集和KEGG信号通路分析;通过过...     

基于WGCNA鉴定茶树响应草甘膦相关的基因共表达模块

【目的】分析茶树响应草甘膦相关的基因表达规律和调控途径,在转录水平上探究草甘膦对茶树的作用,确定茶树响应草甘膦的关键基因。【方法】以茶树品种‘金观音’为试验材料,将推荐使用浓度的草甘膦施于茶树土壤基质表面,经0、0.25、1、3和7 d后,取叶片进行转录组测序,并测定莽草酸含量。利用WGCNA方法联合分析转录组和莽草酸含量数据,鉴定与草甘膦响应相关的共表达基因模块,筛选关键调控基因。【结果】茶树叶片中的莽草酸含量在草甘膦处理后0.25、1和3 d降低,而在7 d时大量积累（为未处理的6.99倍）。从表达谱数据中筛选得到12 568个差异表达基因（DEGs）,草甘膦处理不同时间点与未处理数据比对的DEGs均显著富集在苯丙烷、类黄酮生物合成及植物激素信号转导途径;此外,草甘膦处理分别诱导茶树莽草酸代谢及其下游苯丙烷、类黄酮生物合成和激素信号转导途径相关的24、52、31和69个基因差异表达。通过加权基因共表达网络（WGCNA）方法鉴定得到19个基因模块,将转录组与莽草酸含量数据相关联,筛选到两个与草甘膦响应高度相关的关键基因模块,分别包含2 024和2 305个基因。选取关键模块中连通度最高的前50个基因进行共表达分析,获得6个关键调控基因,包括2个抗性基因（SHMT和RPM）、1个耐药性基因（PDR）、1个离子转运基因（At）、1个膜转运基因（GPT）和1个转录因子（ERF）。【结论】草甘膦通过干扰茶树叶片中莽草酸代谢,影响其下游代谢途径苯丙烷、类黄酮生物合成及激素信号转导途径的基因转录。此外,研究还鉴定到2个草甘膦响应密切相关的共表达模块,发现SHMT、RPM、At、PDR、ERF和GPT等多个潜在候选基因和转录因子在茶树抵御草甘膦逆境中发挥重要作用。     

刺参补体AjC3活性相关miRNA的筛选与初步研究

为了从miRNA调控方面探索刺参Apostichopus japonicus的免疫机制,用脂多糖刺激刺参,在第0、6、12 h取样(记为样本A00、A06、A12),以海胆Strongylocentrotus purpuratus基因组为参照基因,运用高通量测序和相关生物信息学方法筛选与刺参补体(AjC3)相关的miRNAs,并通过实时定量PCR验证相关miRNA和AjC3 mRNA的表达变化。结果表明:经比对筛选出刺参体腔细胞的miRNA共41个;通过靶基因预测和差异表达分析,共得到靶基因为补体AjC3的miRNA有3个,分别为spu-miR-133、spu-miR-137、spu-miR-2006;3种miRNA和AjC3 mRNA实时定量PCR显示,spu-miR-133表达量呈先升高后降低的趋势,12 h达到最高值,spu-miR-137表达量呈先降低后升高的趋势,9 h达到最低值,12 h达到最高值,spu-miR-2006表达量呈先升高后降低的趋势,9 h达到最高值,AjC3 mRNA表达量呈先升高后降低的趋势,9 h达到最高值。研究表明,spu-miR-133、spu-miR-137可能共同抑制调控补体AjC3基因表达,spu-miR-2006在免疫反应中对其他靶基因进行调控。     

南阳黑猪肌内脂质沉积相关lncRNA鉴定和功能预测

为了探讨长链非编码RNA(Long non-coding RNA,lncRNA)在南阳黑猪肌内脂肪沉积中的重要作用，通过分析已发表的南阳黑猪高脂组和低脂组背最长肌转录组数据，通过编码潜能预测获得可能的lncRNA,结合DESeq2软件筛选出在南阳黑猪高脂组和低脂组背最长肌差异表达基因和lncRNA,并通过GO和KEGG通路富集分析，找出与南阳黑猪肌内脂肪沉积相关的差异表达lncRNA。结果表明，通过编码潜能预测获得61个新的lncRNAs,新鉴定lncRNAs和已知lncRNAs与编码基因相比较表现出一些典型特征，如转录本长度较短、外显子数较少。差异表达分析发现，814个编码基因和98个lncRNA在低脂和高脂组背最长肌间差异表达。差异表达的基因参与了多种与脂质代谢相关的信号通路，如甘油三脂代谢、PPAR信号通路、脂肪酸生物合成、脂肪酸降解等。联合差异表达lncRNAs的靶基因预测和差异表达基因结果分析表明，一些lncRNA可能作用于潜在的靶基因，参与脂质代谢过程，调控背最长肌脂质沉积。研究结果为进一步深入研究与猪脂质沉积相关的lncRNA生物学功能及调控机制奠定了基础，在未来育种中改...