中国农业科学 ›› 2021, Vol. 54 ›› Issue (3): 662-674.doi: 10.3864/j.issn.0578-1752.2021.03.019

• 畜牧·兽医·资源昆虫 • 上一篇    

BLV-miRNA跨界调控人类靶基因预测及生物信息学分析

王雍(),李思妍,何思锐,张迪,连帅,王建发,武瑞()   

  1. 黑龙江八一农垦大学动物科技学院/黑龙江省牛病防制重点实验室,黑龙江大庆 163319
  • 收稿日期:2020-02-23 接受日期:2020-07-29 出版日期:2021-02-01 发布日期:2021-02-01
  • 通讯作者: 武瑞
  • 作者简介:王雍,Tel:13251599676;E-mail: bywy0209@126.com
  • 基金资助:
    国家自然科学基金(2041340046);黑龙江省自然科学基金(YQ2019C014);研究生创新科研项目(YJSCX2019-Y39)

Prediction and Bioinformatics Analysis of BLV-miRNA Transboundary Regulation of Human Target Genes

WANG Yong(),LI SiYan,HE SiRui,ZHANG Di,LIAN Shuai,WANG JianFa,WU Rui()   

  1. College of Animal Science and Veterinary Medicine, Heilongjiang Bayi Agricultural University/Heilongjiang Provincial Key Laboratory of Prevention and Control of Bovine Diseases, Daqing 163319, Heilongjiang
  • Received:2020-02-23 Accepted:2020-07-29 Online:2021-02-01 Published:2021-02-01
  • Contact: Rui WU

摘要:

【目的】评估牛白血病病毒(BLV)来源的miRNAs跨界调控人源基因的风险。对BLV-miRNA可能带来的食品安全问题及对人体健康可能造成何种影响进行前瞻性研究,为未来实际生产中地方流行性白血病防控措施执行的必要性研究奠定基础,对BLV与人类疾病间关联性的研究提供理论指导。【方法】首先使用mirbase网站对BLV miRNA的成熟序列进行查询,通过miRanda软件对BLV编码的10种miRNA(BLV-miR-B1-3P,5P、BLV-miR-B2-3P,5P、BLV-miR-B3-3P,5P、BLV-miR-B4-3P,5P、BLV-miR-B5-3P,5P)进行靶基因预测,并选取每个BLV-miRNA评分前10的候选靶基因(去除重复基因后共88个)进行功能分析,对受到多个BLV miRNA共同调控的候选靶基因使用RNAhybrid软件进行二次预测验证,并对其功能进行分析。【结果】 BLV编码的10种miRNA经预测后分别获得1 630—16 383个靶基因不等。对评分前十的共计88个候选靶基因进行功能分析后发现,其中18个基因无相关功能报道;36个候选靶基因与肿瘤性疾病的发生发展存在相关性。2个候选靶基因可以对细胞周期起调控作用;16个候选靶基因参与细胞信号转导的调控;14个候选靶基因在细胞结构/骨架蛋白的形成中发挥作用;细胞的增殖与凋亡的功能表现成拮抗关系,往往促进增殖的基因同时也可以抑制细胞凋亡,共有13个基因对细胞的增殖和凋亡起调节作用,有趣的是,这13的候选靶基因对细胞增殖凋亡功能的调节是双向性的,但不能明确BLV miRNA对细胞的调节到底是更趋向于增殖还是凋亡,因此仍需要后续研究深入探讨;2个候选靶基因对细胞分化起调节作用;16个候选靶基因对细胞的迁移/侵袭功能起调节作用,再次提示BLV miRNA与肿瘤性疾病可能存在更重要的关联性。7个候选靶基因可能在乳腺细胞的分化、迁移、侵袭过程中发挥重要作用,提示BLV与人乳腺癌相关性的研究中,可以从BLV miRNA的角度深入探讨;BLV-B4-3P的2个候选靶基因Ⅰ型胶原α1链基因(COL1A1)、断裂点簇集区(BCR)对人急性淋巴细胞白血病(ALL)具有调节作用。此外,可以被多个BLV miRNA共同靶向的候选靶基因均属于黏蛋白家族(MUC5B、MUC12和 MUC16),且均可以在结肠中表达,对结肠黏膜的形成产生影响。【结论】外源性BLV miRNA可能跨界调控细胞周期、信号转导、结构/骨架、增殖、凋亡、分化、迁移/侵袭相关等细胞功能相关基因,破坏细胞结构;BLV miRNA与人乳腺癌的相关性可能表现在人乳腺癌细胞的分化、迁移、和侵袭过程中;而BLV-miR-B4-3p本身与白血病相关miR 29a共享同一种子区域,可能对人急性淋巴细胞白血病的发生发展造成影响;外源性BLV miRNA具有靶向抑制黏蛋白基因(MUC5B、MUC12、MUC16)表达,通过破坏肠黏膜形成这一途径,跨界调控人源基因的风险。

关键词: 生物信息学, 牛白血病病毒, 跨界调控, 人源基因, miRNA

Abstract:

【Objective】To assess risk of regulation of human-derived genes by miRNAs derived from bovine leukemia virus (BLV), the prospective research on the possible food safety problems and the possible impact on human health caused by BLV-miRNA were carried out, which would lay the foundation for the necessary research on the implementation of Enzootic Bovine Leukosis (EBL) prevention and control measures in actual production in the future, and provide theoretical guidance for the study of the relationship between BLV and human diseases.【Method】 In this study, the mature sequence of BLV miRNAs was first queried using mirbase website, and the miRanda software was used to predict target genes. The predictive 10 miRNAs (BLV-miR-B1-3P,5P, BLV-miR-B2-3P, 5P, BLV-miR-B3-3P,5P, BLV-miR-B4-3P,5P, and BLV-miR-B5-3P,5P) were encoded by BLV. The top 10 candidate target genes of each BLV-miRNA score were selected for functional analysis, including a total of 88 duplicated genes. The candidate target genes co-regulated by multiple BLV miRNAs were verified by secondary prediction using RNAhybrid software, and their functions were analyzed. 【Result】The ten miRNAs encoded by BLV were predicted to obtain 1 630-16 383 target genes, respectively. After functional analysis of eighty-eight candidate target genes in the top ten, it was found that eighteen of them had no relevant functional reports. Thirty-six candidate target genes were related to the occurrence and development of neoplastic diseases. Two candidate target genes could regulate cell cycle. Sixteen candidate target genes were involved in the regulation of cell signal transduction. Fourteen candidate target genes played a role in the formation of structure/cytoskeleton proteins. The function of cell proliferation and apoptosis showed an antagonistic relationship and the genes that often promoting proliferation could also suppress apoptosis. A total of thirteen genes played a regulatory role in cell proliferation and apoptosis. Interestingly, the regulation of the thirteen candidate target genes on cell proliferation and apoptosis was bidirectional. However, it was not clear whether the regulation of BLV miRNA towards cells was more prone to proliferation or apoptosis, so further studies were still needed to discuss in depth. Two candidate target genes could regulate cell differentiation. The sixteen candidate target genes played a role in regulating cell migration/invasion function, again suggesting that BLV miRNA might have a more important correlation with neoplastic diseases. The seven candidate target genes might play an important role in the differentiation, migration and invasion of breast cells, suggesting that the study on the correlation between BLV and human breast cancer could be further discussed from the perspective of BLV miRNA. Two candidate target genotypes of BLV-B4-3P, Collagen 1 chain gene (COL1A1), had a regulatory effect on human acute lymphoblastic leukemia (ALL). In addition, candidate target genes that could be co-targeted by multiple BLV miRNAs belong to the mucin family (MUC5B, MUC12 and MUC16), and be expressed in the colon, influencing the formation of colon mucosa.【Conclusion】Exogenous BLV miRNA might transboundary regulate cell cycle signal transduction structure/cytoskeleton proliferation apoptosis differentiation migration/invasion related cell function related genes and destroy cell structure. The correlation between BLV miRNA and human breast cancer might be shown in the process of differentiation, migration and invasion of human breast cancer cells. BLV-miR-B4-3p shared a seed sequence with miR 29a, which might affect the occurrence and development of human acute lymphoblastic leukemia. Exogenous BLV miRNA had the target of inhibiting the expression of mucin genes, such as MUC5B, MUC12, and MUC16, through the destruction of intestinal mucosa formation to achieve transboundary regulation of human gene risk.

Key words: bioinformatics, bovine leukemia virus, transboundary regulation, human-derived genes, miRNA