感染产肠毒素大肠杆菌的猪小肠上皮细胞microRNA表达谱分析

doi:10.16431/j.cnki.1671-7236.2021.11.034

摘要/Abstract

摘要： 试验旨在探索产肠毒素大肠杆菌（enterotoxigenic Escherichia coli，ETEC）感染猪小肠上皮细胞（IPEC-J2）诱导的microRNA （miRNA）表达谱变化，为解析宿主miRNA在ETEC感染过程中的调控作用提供理论基础。利用Illumina 6000 Novoseq SE50测序平台分别对ETEC感染前后的IPEC-J2进行高通量测序，用Bowtie与参考基因组比对，用DESeq R Package进行miRNA差异性分析。通过miRanda和RNAhybid共同预测差异表达miRNA的靶基因，对差异表达miRNA靶基因进行GO功能和KEGG通路分析。随机选取5个miRNAs，对测序结果进行实时荧光定量PCR验证。结果显示，IPEC-J2在感染前后的sRNA文库经过滤分别得到12 889 260和11 203 056条clean reads。感染前后文库中，miRNA所占比例最高，分别为73.16%和54.10%；分别有97.98%和69.83%长度为18~40 nt的sRNA可比对到参考基因组，表明测序质控良好。长度在22~24 nt的序列大部分首位碱基偏向U，2~8位点出现频率最高的碱基分别为AGCUUAU。共发现311个已知miRNAs，128个新miRNAs。在2个文库中，长度为23 nt的miRNA序列占比最高，分别为41.42%和23.56%。感染后共筛选到140个差异表达miRNAs，其中74个表达上调，66个表达下调。GO分析表明，miRNA靶基因显著富集于代谢过程、正向调节代谢过程、细胞成分或生物合成、免疫系统、细胞内部分和细胞器等功能。KEGG分析表明，差异表达miRNA靶基因显著富集于赖氨酸降解、生产IgA的肠道免疫网络、NF-κB信号通路和T细胞受体信号通路等。实时荧光定量PCR验证结果表明，随机选取的5个miRNAs表达趋势与测序结果一致，表明测序准确可靠。综上所述，IPEC-J2的miRNAs参与了ETEC感染过程，为进一步揭示调控ETEC感染的关键miRNA及其作用机制提供科学依据。

关键词: 猪小肠上皮细胞(IPEC-J2); 产肠毒素大肠杆菌(ETEC); miRNA; 表达谱

Abstract: This study was aimed to explore the changes of microRNA (miRNA) expression profile in porcine small intestine epithelial cells (IPEC-J2) infected with enterotoxigenic Escherichia coli (ETEC), and provide a theoretical basis for the analysis of the regulatory role of host miRNA in the process of ETEC infection. IPEC-J2 pre- and post-infection of ETEC were sequenced by Illumina 6000 Novoseq SE50 sequencing platform, Bowtie was used to compare it with reference genome, and DESeq R Package was used to analysis the miRNA differential expression. The target genes of the differentially expressed miRNA were predicted jointly by miRanda and RNAhybid, and the GO and KEGG analysis of the target genes were performed. Five miRNAs were randomly selected to verify the sequencing results by Real-time quantitative PCR. The results showed that after filtration, 12 889 260 and 11 203 056 clean reads were obtained in the sRNA libraries of IPEC-J2 before and after ETEC infection, respectively. In the pre-infection and post-infection libraries, miRNA accounted for the highest proportion, which were 73.16% and 54.10%. A total of 97.98% and 69.83% of sRNA with length of 18-40 nt matched to the reference genome, indicating the sequencing quality was well controlled. Most of the first bases of the 22-24 nt sequences were inclined to U, and AGCUUAU was the most frequent base at 2-8 sites. A total of 311 known miRNAs and 128 novel miRNAs were discovered. The miRNA sequences of 23 nt length accounted for the highest proportion, taking up 41.42% and 23.56% in the two libraries, respectively. A total of 140 differentially expressed miRNAs were revealed in the comparison between the two libraries, among which 74 were up-regulated and 66 were down-regulated after ETEC infection. GO analysis showed that miRNA target genes were significantly enriched in metabolic processes, positive regulation of metabolic processes, cell composition or biosynthesis, immune system, intracellular parts and organelles. KEGG analysis showed that the differentially expressed miRNA target genes were significantly enriched in lysine degradation, IgA-producing intestinal immune network, NF-κB signaling pathway and T cell receptor signaling pathway. Verification by Real-time quantitative PCR showed that the expression trends of 5 randomly selected miRNAs were consistent with the sequencing results, indicating that the sequencing was accurate and reliable. In conclusion, the miRNAs in IPEC-J2 were involved in the process of ETEC infection, which provided a scientific basis for further revealing the key miRNA regulating ETEC infection and its mechanism of action.

Key words: porcine small intestine epithelial cells (IPEC-J2); enterotoxigenic Escherichia coli (ETEC); miRNA; expression profile

中图分类号:

S852.23

刘迎迎, 张迪, 刘国梅, 李培, 李妍, 左玉柱. 感染产肠毒素大肠杆菌的猪小肠上皮细胞microRNA表达谱分析[J]. 中国畜牧兽医, 2021, 48(11): 4242-4253.

LIU Yingying, ZHANG Di, LIU Guomei, LI Pei, LI Yan, ZUO Yuzhu. Analysis of microRNA Expression Profile in Porcine Small Intestinal Epithelial Cells Infected with Enterotoxigenic Escherichia coli[J]. China Animal Husbandry and Veterinary Medicine, 2021, 48(11): 4242-4253.

参考文献

[1] 王红宁, 代敏.仔猪腹泻的病因及综合防治进展[J]. 肉品卫生, 2004, 8:23-25. WANG H N, DAI M.Progress in etiology and comprehensive control of piglet diarrhea[J]. Meat Hygiene, 2004, 8:23-25.(in Chinese)
[2] KWON D, CHOI C, JUNG T, et al. Genotypic prevalence of the fimbrial adhesins (F4, F5, F6, F41 and F18) and toxins (LT, STa, STb and STx2e) in Escherichia coli isolated from postweaning pigs with diarrhoea or oedema disease in Korea[J]. Veterinary Record, 2002, 150(2):35-37.
[3] LUPPI A, GIBELLINI M, GIN T, et al. Prevalence of virulence factors in enterotoxigenic Escherichia coli isolated from pigs with post-weaning diarrhoea in Europe[J]. Porcine Health Management, 2016, 2:20.
[4] TOLEDO A, GÓMEZ D, CRUZ C, et al. Prevalence of virulence genes in Escherichia coli strains isolated from piglets in the suckling and weaning period in Mexico[J]. Journal of Medical Microbiology, 2012, 61(Pt 1):148-156.
[5] BARIGYE R, GAUTAM A, PICHE L M, et al. Prevalence and antimicrobial susceptibility of virulent and avirulent multidrug-resistant Escherichia coli isolated from diarrheic neonatal calves[J]. American Journal of Veterinary Research, 2012, 73(12):1944-1950.
[6] 葛晨玲, 石大丽, 胡文, 等.广西猪源产肠毒素大肠杆菌毒力因子及多重耐药性分析[J]. 中国畜牧兽医, 2019, 46(3):913-923. GE C L, SHI D L, HU W, et al. Analysis of virulence factors and multidrug resistance of enterotoxin-producing Escherichia coli in Guangxi[J]. China Animal Husbandry & Veterinary Medicine, 2019, 46(3):913-923.(in Chinese)
[7] YAN X, HUANG X, REN J, et al. Adhesion phenotypes of pigs of Chinese and Western breeds and a White Duroc-Erhualian crossbreed with regard to susceptibility to enterotoxigenic Escherichia coli with fimbrial adhesins K99, 987P, and F41[J]. American Journal of Veterinary Research, 2011, 72(1):80-84.
[8] SAYED D, ABDELLATIF M.microRNAs in development and disease[J]. Physiological Reviews, 2011, 91(3):827-887.
[9] BARTEL D P.microRNAs:Genomics, biogenesis, mechanism, and function[J]. Cell, 2004, 116(2):281-297.
[10] LIU X, ZHU L, LIAO S, et al. The porcine microRNA transcriptome response to transmissible gastroenteritis virus infection[J]. PLoS One, 2015, 10(3):e0120377.
[11] YE L, SU X, WU Z, et al. Analysis of differential miRNA expression in the duodenum of Escherichia coli F18-sensitive and -resistant weaned piglets[J]. PLoS One, 2012, 7(8):e43741.
[12] 刘曼迪.猪源产肠毒素大肠杆菌总RNA的免疫保护作用研究[D].保定:河北农业大学, 2019. LIU M D.The immune protection of porcine enterotoxigenic Escherichia coli total RNA[D].Baoding:Hebei Agricultural University, 2019.(in Chinese)
[13] 苗志国, 张金洲, 张晓建, 等.microRNA在猪中的研究进展[J]. 中国畜牧兽医, 2018, 45(6):1626-1632. MIAO Z G, ZHANG J Z, ZHANG X J, et al. Research progress on microRNA in pigs[J]. China Animal Husbandry & Veterinary Medicine, 2018, 45(6):1626-1632.(in Chinese)
[14] WU Z, QIN W, WU S, et al. Identification of microRNAs regulating Escherichia coli F18 infection in Meishan weaned piglets[J]. Biology Direct, 2016, 11(1):59.
[15] MI S, CAI T, HU Y, et al. Sorting of small RNAs into Arabidopsis argonaute complexes is directed by the 5' terminal nucleotide[J]. Cell, 2008, 133(1):116-127.
[16] OHTSUKA M, ARASE H, TAKEUCHI A, et al. NFAM1, an immunoreceptor tyrosine-based activation motif-bearing molecule that regulates B cell development and signaling[J]. Proceedings of The National Academy of Sciences of The United States of America, 2004, 101(21):8126-8131.
[17] ZHOU C, LIU Z, JIANG J, et al. Differential gene expression profiling of porcine epithelial cells infected with three enterotoxigenic Escherichia coli strains[J]. BMC Genomics, 2012, 13:330.
[18] 傅仅, 王芳, 李仕国, 等.FADD、IPO-38与肿瘤的研究进展[J]. 世界最新医学信息文摘, 2015, 15(86):52-53. FU J, WANG F, LI S G, et al. Research progress of FADD, IP-38 and tumor[J]. World Latest Medical Information Abstracts, 2015, 15(86):52-53.(in Chinese)
[19] 乔新安, 陈丽颖, 王艳玲, 等.CCL25/CCR9和CCL28/CCR10在仔猪胃肠道组织中的表达[J]. 中国农学通报, 2007, 23(5):1-5. QIAO X A, CHEN L Y, WANG Y L, et al. Expression of CCL25/CCR9 and CCL28/CCR10 in the gastrointestinal tract of piglets[J]. Chinese Agricultural Science Bulletin, 2007, 23(5):1-5.(in Chinese)
[20] VOGLER M, HAMALI H A, SUN X M, et al. BCL2/BCL-X(L) inhibition induces apoptosis, disrupts cellular calcium homeostasis, and prevents platelet activation[J]. Blood, 2011, 117(26):7145-7154.
[21] SUN L, WU S, DAI C H, et al. Insight into the molecular mechanism of miR-192 regulating Escherichia coli resistance in piglets[J]. Bioscience Reports, 2018, 38(1):BSR20171160.
[22] DAI C H, WANG F, WANG S Q, et al. miR-215 Targeting novel genes EREG, NIPAL1 and PTPRU regulates the resistance to E.coli F18 in piglets[J]. Genes, 2020, 11(9):1053
[23] SAILO L, KUMAR A, SAH V, et al. Genome-wide integrated analysis of miRNA and mRNA expression profiles to identify differentially expressed miR-22-5p and miR-27b-5p in response to classical Swine fever vaccine virus[J]. Functional & Integrative Genomics, 2019, 19(6):901-918.
[24] WANG J, CHEN W D, WANG Y D.The relationship between gut microbiota and inflammatory diseases:The role of macrophages[J]. Frontiers in Microbiology, 2020, 11:1065.