肝细胞癌患者血清外泌体microRNA表达鉴定

目的:探讨肝癌细胞外泌体中差异表达的microRNAs(miRNAs)在肝细胞癌(HCC)诊断中的应用价值。方法:通过高通量测序筛选肝癌细胞外泌体中差异表达的miRNAs。实时定量PCR验证差异表达分子;检测差异表达的miRNAs在健康人(Health)、慢性乙型肝炎患者(CHB)、肝硬化患者(LC)及乙型肝炎病毒阳性的肝细胞癌患者(HCC)血清外泌体中的表达。结果:高通量测序筛选到肝癌细胞外泌体中差异表达的miRNA共88种,其中58种表达上调,30种表达下调。选择其中8种差异表达的miRNAs进行q RT-PCR验证,结果显示,此8种miRNAs在细胞上清外泌体、细胞内、癌与癌旁组织中的表达趋势与测序结果一致。miR-221-3p和miR-224-5p在HCC组外泌体中的表达水平显著高于Health组、CHB组和LC组(P0.01),miR-124-3p和let-7a-5p在HCC组外泌体中的表达水平显著低于其他各组(P0.05)。四个组中,miR-21-5p、miR-191-5p、miR-34a-5p和miR-122-5p的表达水平不存在显著性差异(P0.05)。结论:血清外泌体中的miR-221-3p、miR-224-5p、miR-124-3p和let-7a-5p可能成为肝细胞癌的候选标志物。     

猪流行性腹泻病毒感染PK-15细胞miRNA的差异表达分析

为研究猪流行性腹泻病毒(Porcine epidemic diahorrea virus,PEDV)对猪肾传代细胞PK-15miRNA表达谱的影响,将PEDV感染PK-15细胞后,提取总RNA,进行高通量测序,构建感染与未感染PEDV的PK-15细胞的miRNA表达谱,并进行差异性分析,对表达差异显著的miRNA进行heatmap聚类分析及GO(Gene ontology,GO)分子功能(Molecular function,MF)分析,选取10个显著性差异表达的miRNA通过RT-qPCR进行了验证。结果发现,显著性差异表达的miRNA有214个,其中175个miRNA表达上调,39个miRNA表达下调。Heatmap聚类分析表明,病毒组绝大部分差异表达的miRNA较对照组表达量上调,少部分下调。GO分析表明,miRNA广泛参与结合、蛋白结合、蛋白激酶活性、转移酶活性、含磷基团转移、磷酸转移酶活性等生物作用。RT-qPCR验证的各miRNA的表达趋势与高通量测序结果一致。结果表明,猪流行性腹泻病毒感染对PK-15细胞编码的miRNA表达水平有显著影响,从而对进一步研究治疗PEDV的miRNA制剂提供新思路。     

基于数据挖掘分析microRNAs在冠状动脉支架内再狭窄中的调控机制

目的:通过寻找支架内再狭窄相关的循环microRNAs(miRNAs)及其作用靶基因,阐释miRNAs在经皮冠状动脉介入治疗的同时带来的支架内再狭窄(ISR)中的网络调控作用及机制。方法:由美国NCBI的GEO公共数据平台下载GSE60959中的miRNAs芯片数据,通过GeneSpringGX芯片分析软件分析得到ISR患者和对照组之间差异表达的循环miRNAs;采用miRNAs靶基因预测算法TargetScan和miRanda预测差异表达miRNAs的靶基因,通过DAVID平台和KEGG数据库分析得到靶基因所参与的信号通路;用GSE46560的mRNA芯片数据分析得到ISR患者较非ISR患者循环mRNAs差异表达谱,用来验证miRNAs调控ISR所作用的靶基因;采用Cytoscape软件构建基于miRNAs-信号通路、miRNAs-靶基因的共表达网络。结果:与对照组相比,ISR组存在131个差异表达循环miRNAs,其中上调41个。前5个上调miRNAs分别是miR-1183、miR-512-5p、miR-187-5p、miR-144-3p和miR-1225-5p。5个miRNAs的预测靶基因共6587个,信号通路富集分析显示这些基因主要参与的信号通路包括MAPK signaling pathway、ErbB signaling pathway、Wnt signaling pathway、Focal adhesion、Neurotrophin signaling pathway和Regulation of actin cytoskeleton。进一步对mRNAs芯片进行差异分析显示,与对照组相比,ISR患者差异循环mRNAs共171个,其中最可能受上述5个miRNAs调控的基因共43个。构建的共表达网络显示,miR-512-5p是作用靶基因和参与信号通路最多的miRNAs。结论:ISR患者循环miRNAs和mRNAs表达谱均存在改变,多个差异表达miRNAs通过作用于多个靶mRNAs,进而影响多个信号通路的活性,最终对ISR的发生、发展形成网络调控作用。     

意大利蜜蜂幼虫肠道的miRNAs的生物信息学预测及分析

【目的】蜜蜂肠道是食物消化和营养吸收的主要部位,同时也是抵御病原侵染的重要场所。本研究旨在对意大利蜜蜂Apismelliferaligustica(简称意蜂)幼虫肠道的microRNAs(miRNAs)及其靶基因进行深入分析,进而解析miRNAs在幼虫肠道发育和生长过程中的作用。【方法】利用small RNA-seq(sRNA-seq)技术对意蜂4、5和6日龄幼虫肠道进行测序,通过相关生物信息学软件对意蜂幼虫肠道的miRNAs进行预测及分析。利用TargetFinder软件预测miRNAs的靶基因,然后将靶基因通过BLAST软件进行GO和KEGG数据库的功能注释。利用Cytoscape软件构建miRNAs与其靶向结合的mRNAs的调控网络。【结果】意蜂幼虫肠道样品的测序共获得96 329 456条有效标签序列(Clean tags),预测出560个miRNAs,包括45个novelmiRNAs。上述miRNAs的长度主要分布在17-27nt之间,且不同长度的miRNAs的首位碱基偏向性具有明显差异。表达量聚类分析结果显示,331个miRNAs为4、5和6日龄幼虫肠道所共有且表达量较为稳定,随着发育时间延长,特有miRNAs的表达水平呈总体增加的趋势。利用TargetFinder软件共预测出16479个意蜂幼虫肠道的靶基因,其中有8132和3361个可分别注释到GO和KEGG数据库,进一步分析发现有224个靶基因注释在Wnt信号通路,分别有2个和27个靶基因与保幼激素和蜕皮激素的调控密切相关。【结论】本研究在全基因组水平对意蜂幼虫肠道的miRNAs进行预测和分析,研究结果不仅丰富了意蜂的miRNAs信息,也为深入研究意蜂幼虫肠道的生长发育机理打下了初步基础。     

猪圆环病毒2型通过外泌体miR-125a-5p靶向Bcl-2诱导淋巴细胞凋亡

为研究miR-125a-5p在猪圆环病毒2型(porcine circovirus type 2,PCV2)诱导淋巴细胞凋亡中的作用及其作用机制,以PCV2感染PK-15细胞外泌体孵育的淋巴细胞为研究对象,采用流式细胞术、蛋白质免疫印迹试验(Western blotting)和实时荧光定量PCR,检测淋巴细胞凋亡率及凋亡相关miRNA表达;合成miR-125a-5p模拟物和抑制物转染PK-15细胞,检测miR-125a-5p过表达或抑制表达后细胞凋亡率;采用生物信息学方法预测miR-125a-5p的靶基因,双荧光素酶报告基因检测miR-125a-5p对靶基因的调控;Western blotting检测外泌体孵育淋巴细胞的线粒体凋亡信号通路相关蛋白Bcl-2、Bax、细胞色素C和caspase-3的表达。结果显示,感染PCV2的PK-15细胞分泌的外泌体极显著提高淋巴细胞凋亡率,在一定浓度范围内呈剂量依赖性;与PCV2诱导细胞凋亡相关的miRNA中,miR-125a-5p表达量极显著升高,miR-125a-5p模拟物转染细胞后极显著提高细胞凋亡率;利用TargetScan预测发现,miR-125a-5p与Bcl-2 3''UTR区有结合位点,miR-125a-5p模拟物极显著抑制pmir-Bcl-2 3''UTR-WT荧光素酶活性,对pmir-Bcl-2 3''UTR-MuT的荧光素酶活性无明显改变;外泌体孵育的淋巴细胞Bcl-2表达量显著降低,Bax、细胞色素C的释放和caspase-3表达量显著升高,Bcl-2/Bax的比值极显著降低。这表明,PCV2通过外泌体诱导淋巴细胞上调miR-125a-5p的表达,进而抑制Bcl-2 mRNA和蛋白表达,激活淋巴细胞线粒体凋亡通路诱导细胞凋亡。     

Bmo-miR-375-3p对家蚕丝素蛋白轻链基因BmFib-L表达的调控作用

【目的】探究家蚕Bombyx mori miRNAs对丝素蛋白轻链基因(fibroin light chain gene,Bm FibL)和P25蛋白基因Bm P25表达的调控作用。【方法】分别以Bm Fib-L和Bm P25 mRNA 3'UTR为靶序列,应用在线预测软件RNAhybrid从前期家蚕4和5龄幼虫后部丝腺(posterior silk gland)sRNA高通量测序获得的29个差异表达bmo-miRNAs中,筛选对Bm Fib-L和Bm P25具有调控作用的候选bmo-miRNAs;采用半定量RT-PCR方法在mRNA水平分析候选bmo-miRNAs及其靶基因Bm Fib-L和Bm P25在4和5龄幼虫期以及5龄第3天幼虫不同组织的表达水平;利用双荧光报告基因检测系统在家蚕细胞Bm N中验证候选bmo-miRNAs对Bm Fib-L和Bm P25表达的调控作用。【结果】筛选到一个在Bm Fib-L和Bm P25 mRNA 3'UTR上都有靶位点的bmo-miR-375-3p(曾称bmo-miR-375*)。在后部丝腺中bmo-miR-375-3p和其预测靶基因Bm Fib-L和Bm P25都有较高的表达水平,提示bmo-miR-375-3p具备调控Bm Fib-L和Bm P25表达的时空条件。miR-375-3p重组表达载体与融合了Bm Fib-L 3'UTR的萤火虫荧光素酶(luciferase)报告基因(luc)表达载体共转染Bm N细胞,报告基因luc的表达水平显著下降;而在miR-375-3p重组表达载体与融合了Bm P25 3'UTR的luc表达载体共转染Bm N细胞中报告基因luc的表达水平下降不显著。【结论】miR-375-3p显著下调Bm Fib-L基因的表达,而对Bm P25基因的表达无明显调控作用。     

miR-15a靶基因的预测及生物信息学分析

目的:对目前研究较为广泛的miR-15a的靶基因进行预测及相关生物信息学分析,以期为miR-15a靶基因的实验验证提供数据支持,并为深入研究miR-15a的调控机制及生物学功能奠定基础和提供理论指导。方法:选择TargetScan5.1与PicTar两种计算方法预测miR-15a的靶基因的交集作为分析的基因集合,分别进行GO注释描述、GO富集分析和生物通路富集分析。结果与结论:预测靶基因集合分别富集在转录调控、蛋白质修饰、细胞周期等生物学过程和蛋白激酶活性等分子功能上(P0.01);经典miR-15a预测靶基因集合显著富集于KEGG通路数据库中的Wnt信号通路、细胞周期和p53信号通路等5个信号转导通路及前列腺癌、慢性髓细胞性白血病、黑素瘤等7个疾病通路中(P0.05)。     

家蚕miR-301对化学感受蛋白基因csp9表达的调控

【目的】探索家蚕Bombyx mori miRNAs对化学感受蛋白基因表达的调控作用,以进一步研究miRNAs及其靶基因在昆虫化学识别中的作用。【方法】利用生物信息学方法预测和筛选可能作用于家蚕化学感受蛋白CSPs基因家族成员的miRNAs;实时荧光定量PCR分析预测获得的候选miR-301和其作用的靶基因在家蚕成虫不同组织中的表达变化;构建miR-301预测靶基因3′-UTR的双荧光素酶报告载体,与合成的miR-301 mimics或阴性对照转染人胚肾细胞HEK293,通过双荧光素酶报告基因检测系统中荧光素酶活性变化检测miR-301对其靶基因表达的调控作用。【结果】生物信息学分析结果发现,家蚕化学感受蛋白基因csp9是miR-301的预测靶基因,二者的结合位点位于csp9的3′-UTR区。实时荧光定量PCR检测结果表明,miR-301在交配后家蚕雌雄成虫触角和雄成虫头部都显著上调表达,靶基因csp9在对应组织中表达则显著下调。二者共转染HEK293细胞后,双荧光素酶检测结果表明,miR-301可以通过与csp9 3′-UTR的互作,显著抑制上游荧光素酶报告基因的表达。【结论】家蚕化学感受蛋白基因csp9是miR-301的靶基因,miR-301通过与靶基因3′-UTR的结合,在翻译水平上抑制csp9的表达。     

miR-17-92基因簇microRNAs对哺乳动物器官发育及肿瘤发生的调控

microRNAs(miRNAs)是近年发现的一种高度保守的非编码小RNA, 它们通过抑制靶基因mRNA的翻译或将其降解, 在转录后水平调控基因的表达, 参与调控哺乳动物多个器官的发育过程和人类疾病的发生。miR-17-92基因簇是一个高度保守的基因簇, 编码miR-17-5p、miR-17-3p、miR-18a、miR-19a、miR-20a、miR-19b-1和miR-92-1等7个miRNAs。大量证据表明, miR-17-92基因簇miRNAs参与了心、肺、免疫系统的发育、血管生长及前脂肪细胞的分化等过程。此外, miR-17-92基因簇miRNAs在多种肿瘤中高表达, 能作为致癌基因诱发淋巴瘤和血管化肿瘤的发生, 但它也可以作为抑癌基因抑制乳腺癌细胞的增殖。文章对miR-17-92基因簇miRNAs在哺乳动物器官发育及肿瘤发生中的作用进行综述     

草地贪夜蛾取食Cry1Ab和Cry1Fa蛋白后中肠转录组及ABC基因表达分析

【目的】为揭示草地贪夜蛾Spodoptera furgiperda幼虫取食Bt蛋白后与中肠上相关ATP结合盒转运子(ATP-binding cassette transporter, ABC)蛋白基因的表达量变化的关系。【方法】分别使用含活化晶体蛋白Cry1Ab (LC₇₀=240.2 μg/g)和Cry1Fa (LC₇₀=270.0 μg/g)蛋白的人工饲料饲喂草地贪夜蛾4龄幼虫48 h,利用高通量测序对中肠进行转录组测序并进行生物信息学分析,筛选处理后差异表达基因;利用RT-qPCR验证差异表达ABC基因的表达量。【结果】与饲喂正常人工饲料的对照相比,饲喂含240.2 μg/g Cry1Ab和270.0 μg/g Cry1Fa的人工饲料后草地贪夜蛾4龄幼虫中肠转录组中分别检测到1 305和1 202个差异表达基因。Cry1Ab和Cry1Fa处理组与对照组之间分别有994和912个差异表达基因被GO功能注释到生物学过程、分子功能和细胞组分三大类。在最终筛选到的9个差异表达的ABC家族基因中,Cry1Ab处理组与对照组之间有4个差异表达ABC基因,3个上调,1个下调; Cry1Fa处理组与对照组之间有5个差异表达ABC基因,2个上调,3个下调;Cry1Ab和Cry1Fa处理组与对照组之间有2个ABC基因(LOC118267200和LOC118267201)表达量均显著上调。RT-qPCR验证结果表明,与对照组相比,Cry1Ab处理组有3个ABC基因表达量极显著上调,2个ABC基因表达量下调;Cry1Fa处理组有5个ABC基因表达量上调,1个ABC基因表达量下调。【结论】Cry1Ab和Cry1Fa蛋白的摄入可以影响草地贪夜蛾幼虫中肠一些ABC家族基因的表达量变化,这些基因的表达量变化与昆虫抗性产生有关。经比对后发现,ABCC家族与ABCG8基因表达量变化显著。本研究为下一步明确草地贪夜蛾体内ABC转运蛋白在Bt蛋白杀虫机制中的作用,以及合理使用Bt蛋白防治草地贪夜蛾及延缓抗性提供了理论依据。     

Identification and bioinformatics analysis of microRNAs associated with stress and immune response in serum of heat-stressed and normal Holstein cows

MicroRNAs (miRNAs) are small single-stranded non-coding RNAs that have an important regulatory function in animal growth and developmental processes. However, the differential expression of miRNA and the role of these miRNAs in heat-stressed Holstein cows are still unknown. In this study, the profile of differentially expressed miRNAs and the target genes analysis in the serum of heat-stressed and normal Holstein cows were investigated by a Solexa deep-sequencing approach and bioinformatics. The data identified 52 differentially expressed miRNAs in 486 known miRNAs which were changed significantly between heat-stressed and normal Holstein cows (fold change >2, P < 0.001). Target genes analysis showed that at least 7 miRNAs (miR-19a, miR-19b, miR-146a, miR-30a-5p, miR-345-3p, miR-199a-3p, and miR-1246) were involved in the response to stress, oxidative stress, development of the immune system, and immune response among the identified 52 differentially expressed miRNAs. Five miRNAs (miR-27b, miR-181a, miR-181b, miR-26a, and miR-146b) were involved in stress and immune responses and the expression of five miRNAs was striking (P < 0.001). In addition, RT-qPCR and deep-sequencing methods showed that 8 miRNAs among the 12 selected miRNAs (miR-19a, miR-19b, miR-27b, miR-30a-5p, miR-181a, miR-181b, miR-345-3p, and miR-1246) were highly expressed in the serum of heat-stressed Holstein cows. GO and KEGG pathway analysis showed that these differentially expressed miRNAs were involved in a pathway that may differentially regulate the expression of stress response and immune response genes. Our study provides an overview of miRNAs expression profile and the interaction between miRNAs and their target genes, which will lead to further understanding of the important roles of miRNAs in heat-stressed Holstein cows.     

Analysis of microRNA expression profiles in porcine PBMCs after LPS stimulation

In the present study, we used microRNA (miRNA) sequencing to discover and explore the expression profiles of known and novel miRNAs in 1000 ng/ml LPS stimulated for 8 h vis-à-vis non-stimulated (i.e. control) PBMCs isolated from the blood of healthy pigs. A total of 291 known miRNAs were bio-computationally identified in porcine PBMCs, and 228 novel miRNAs (not enlisted in the swine mirBase) were identified. Among these miRNAs, ssc-miR-148a-3p, ssc-let-7g, ssc-let-7f, 3_8760, ssc-miR-26a, ssc-miR-451, ssc-miR-21, ssc-miR-30d, ssc-miR-99a and ssc-miR-103 were the top 10 most abundant miRNAs in porcine PBMCs. Through miRNA differential analysis combined with quantitative PCR, we found the expressions of ssc-miR-122, ssc-miR-129b, ssc-miR-17-5p and ssc-miR-152 were significantly changed in porcine PBMCs after LPS stimulation. Furthermore, targets prediction and function analysis indicated a significant enrichment in gene ontology functional categories related to diseases, immunity and inflammation. In conclusion, this study on profiling of miRNAs expressed in LPS-stimulated PBMCs provides an important reference point for future studies on regulatory roles of miRNAs in porcine immune system.     

MicroRNA Transcriptome Profile Analysis in Porcine Muscle and the Effect of miR-143 on the MYH7 Gene and Protein

Porcine skeletal muscle fibres are classified based on their different physiological and biochemical properties. Muscle fibre phenotype is regulated by several independent signalling pathways, including the mitogen-activated protein kinase (MAPK), nuclear factor of activated T cells (NFAT), myocyte enhancer factor 2 (MEF2) and peroxisome proliferator-activated receptor (PPAR) signalling pathways. MicroRNAs are non-coding small RNAs that regulate many biological processes. However, their function in muscle fibre type regulation remains unclear. The aim of our study was to identify miRNAs that regulate muscle fibre type during porcine growth to help understand the miRNA regulation mechanism of fibre differentiation. We performed Solexa/Illumina deep sequencing for the microRNAome during 3 muscle growth stages (63, 98 and 161 d). In this study, 271 mature miRNAs and 243 pre-miRNAs were identified. We detected 472 novel miRNAs in the muscle samples. Among the mature miRNAs, there are 23 highest expression miRNAs (over 10000 RPM), account for 85.3% of the total counts of mature miRNAs., including 10 (43.5%) muscle-related miRNAs (ssc-miR-133a-3p, ssc-miR-486, ssc-miR-1, ssc-miR-143-3p, ssc-miR-30a-5p, ssc-miR-181a, ssc-miR-148a-3p, ssc-miR-92a, ssc-miR-21, ssc-miR-126-5p). Particularly, both ssc-miR-1 and ssc-miR-133 belong to the MyomiRs, which control muscle myosin content, myofibre identity and muscle performance. The involvement of these miRNAs in muscle fibre phenotype provides new insight into the mechanism of muscle fibre regulation underlying muscle development. Furthermore, we performed cell transfection experiment. Overexpression/inhibition of ssc-miR-143-3p in porcine skeletal muscle satellite cell induced an/a increase/reduction of the slow muscle fibre gene and protein (MYH7), indicating that miR-143 activity regulated muscle fibre differentiate in skeletal muscle. And it regulate MYH7 through the HDAC4-MEF2 pathway.     

MicroRNA expression profile in RAW264.7 cells in response to Brucella melitensis infection

MicroRNA (miRNA) is small non-coding RNA with approximate 22 nt in length. Recent studies indicate that miRNAs play significant roles in pathogen-host interactions. Brucella organisms are Gram-negative facultative intracellular bacteria that cause Brucellosis. Brucella strains infect macrophages and establish chronic infection by altering host life activities including apoptosis and autophagy. Here, we report a comprehensive analysis of miRNA expression profiles in mock- and Brucella-infected RAW264.7 cells using high-throughput sequencing approach. In total, 344 unique miRNAs were co-expressed in the two libraries, in which 57 miRNAs were differentially expressed. Eight differentially expressed miRNAs with high abundance were subjected to further analysis. The GO enrichment analysis suggests that the putative target genes of these differentially expressed miRNAs are involved in apoptosis, autophagy and immune response. In particular, a total of 25 target genes are involved in regulating apoptosis and autophagy, indicating that these miRNAs may play important regulatory roles in the Brucella-host interactions. Furthermore, the interactions of miR-1981 and its target genes, Bcl-2 and Bid, were validated by luciferase assay. The results show that miR-1981 mimic up-regulated the luciferase activity of psiCHECK-2 Bcl-2 3' UTR, but the luciferase activity of psiCHECK-2 Bid 3' UTR was not changed significantly. Taken together, these data provide valuable framework on Brucella induced miRNA expression in RAW264.7 cells, and suggest that Brucella may establish chronic infection by regulating miRNA expression profile.     

感染GCRV的草鱼miRNA测序与比较分析

为探索基于高通量筛选抗性主效miRNA功能分子的方法应用于抗性草鱼(Ctenopharyngodon idellus)的选育,研究借助高通量测序鉴定分析感染GCRV前后的草鱼头肾组织中miRNAs.测序结果共鉴定出821个成熟miRNAs,其中118个在GCRV攻毒组特异表达,82个为未攻毒组特有;差异分析结果显示,G...     

意大利蜜蜂幼虫肠道在球囊菌胁迫后期的差异表达微小RNA及其靶基因分析

【目的】蜜蜂球囊菌Ascosphaera apis(简称球囊菌)是一种特异性侵染蜜蜂幼虫肠道的致死性真菌病原。微小RNA(microRNA,mi RNA)是一种在转录后水平对m RNA进行负调控的关键调控因子。本研究旨在对意大利蜜蜂Apis mellifera ligustica(简称意蜂)幼虫肠道在球囊菌胁迫后期的差异表达mi RNAs(DEmiRNAs)及其靶基因进行深入分析,为揭示DEmiRNAs在胁迫应答中的作用提供重要信息。【方法】利用small RNA-seq(s RNA-seq)技术对正常及球囊菌侵染的意蜂6日龄幼虫肠道(分别表示为Am CK和Am T)进行测序,通过相关生物信息学软件对DEmiRNAs进行预测和分析。利用TargetFinder软件预测DEmiRNAs的靶基因,然后利用Blast软件对靶基因进行GO和KEGG数据库的功能注释。通过Cytoscape软件构建DEmiRNAs与其靶m RNAs的调控网络。通过茎环实时荧光定量PCR(stem-loop RT-qPCR)验证测序数据的可靠性。【结果】Am CK和Am T的测序分别得到9 230 496和10 823 667条有效序列标签; Am CK vs Am T比较组中包含15个上调和6个下调mi RNAs,分别结合3 503和3 252个靶基因,它们可分别注释到40和39个GO terms以及104和99条代谢通路;进一步分析发现调控网络中的17个DEmiRNAs靶向结合116个与丝氨酸蛋白酶相关的m RNAs,14个DEmiRNAs靶向结合54个与泛素介导的蛋白水解相关的m RNAs。Stem-loop RT-qPCR验证结果显示随机选择的4个DEmiRNAs (mi R-251-x,mi R-9277-y,mi R-1672-x和mi R-4968-y)的表达量变化趋势与测序结果一致,表明测序数据真实可信。【结论】本研究率先对意蜂幼虫肠道在球囊菌胁迫后期的DEmiRNAs及其靶基因进行预测与分析,提供了mi RNAs的表达谱和差异表达信息。ame-miR-927b,mi R-429-y和mi R-8440-y等DEmiRNAs可能参与对丝氨酸蛋白酶的调控和对泛素介导的蛋白水解的调控,DEmiRNAs与靶基因之间存在复杂的调控网络关系。DEmiRNAs可参与到意蜂幼虫肠道与球囊菌间的互作。