反刍动物Toll样受体多基因家族的分子进化及表达模式分析

doi:10.11843/j.issn.0366-6964.2023.09.013

摘要/Abstract

摘要： 为了了解Toll样受体(Toll-like receptors,TLRs)多基因家族在反刍动物中的分子进化关系及表达分析。本研究基于6个科53种反刍动物基因组,利用同源基因对反刍动物TLR1-10基因进行系统鉴定,分析其染色体定位和基因结构、进化关系、选择压力(ω),并结合转录组数据明晰TLR基因在各组织的表达模式。研究结果表明:1)反刍动物的TLR1-10均为单拷贝基因,按进化关系分为TLR1 (TLR1/6/10)、TLR2、TLR4、TLR3 (TLR3/5)、TLR7 (TLR7/8/9)亚家族。2)反刍动物10个TLR基因整体上经受了纯化选择作用(ω<1),但共有35个氨基酸位点受到强烈的正选择作用;牛科比鹿科具有更高的ω值和更高比例的正选择位点,非病毒性TLRs比病毒性TLRs具有更高的ω和更多的正选择位点,提示鹿科TLR以及病毒性的TLR受到较强的进化约束力。3)在绵羊中,10个TLRs基因在免疫组织中均有所表达,其中TLR2和TLR9在外周血单个核细胞(PBMC)中存在高度表达。结果提示,反刍动物TLRs基因家族在进化水平上是相对保守的,均受到强烈的纯化选择作用,TLR2和TLR9在PBMC组织起主要的免疫作用。

关键词: Toll样受体, 反刍动物, 系统进化, 选择压力, 表达模式

Abstract: The study aimed to understand the molecular evolutionary relationships and expression analysis of the multigene family of Toll-like receptors (TLRs) in ruminants. The TLR1-10 genes was systematically identified in ruminants based on the genomes of 53 ruminant species from 6 families, their chromosomal localization and gene structure, evolutionary relationships, selection pressure (ω) was analyzed, and the expression patterns of TLR genes was clarified combined with transcriptomic data in various tissues. The results showed that:1) TLR1-10 in ruminants were all single-copy genes, which were classified into TLR1 (TLR1/6/10), TLR2, TLR4, TLR3 (TLR3/5) and TLR7 (TLR7/8/9) subfamilies according to evolutionary relationships. 2) The 10 TLR genes in ruminants as a whole underwent purifying selection (ω<1), but a total of 35 amino acid loci were subject to strong positive selection; bovidae had higher ω values and a higher proportion of positive selection loci than cervidae, and non-viral TLRs had higher ω and more positive selection loci than viral TLRs, suggesting that cervidae TLRs as well as viral TLRs were subject to stronger evolutionary binding. 3) In sheep, all 10 TLR genes were expressed in immune tissues, with TLR2 and TLR9 expressed at relatively high levels in PBMC. The results suggest that the gene family of ruminant TLRs is relatively conserved at the evolutionary level and both are subject to strong purifying selection, with TLR2 and TLR9 playing a major immune role in PBMC tissues.

Key words: Toll-like receptor, ruminant, phylogeny, selection pressure, expression pattern

中图分类号:

S813.2

马淑娟, 徐祎洁, 何珂, 马瑞丰, 朱英. 反刍动物Toll样受体多基因家族的分子进化及表达模式分析[J]. 畜牧兽医学报, 2023, 54(9): 3722-3734.

MA Shujuan, XU Yijie, HE Ke, MA Ruifeng, ZHU Ying. Molecular Evolution and Expression Patterns of a Multigene Family of Toll-like Receptors in Ruminants[J]. Acta Veterinaria et Zootechnica Sinica, 2023, 54(9): 3722-3734.

参考文献

[1] FITZGERALD K A,KAGAN J C.Toll-like receptors and the control of immunity[J].Cell,2020,180(6):1044-1066.
[2] ZHANG E J,MA Z Y,LU M J.Contribution of T-and B-cell intrinsic toll-like receptors to the adaptive immune response in viral infectious diseases[J].Cell Mol Life Sci,2022,79(11):547.
[3] 周明.鸟类TLRs基因的适应性进化研究[D].雅安:四川农业大学,2019.
ZHOU M.Adaptive evolution of TLRs in birds[D].Yaan:Sichuan Agricultural University,2019.(in Chinese)
[4] BEUTLER B,REHLI M.Evolution of the TIR,Tolls and TLRs:functional inferences from computational biology[J]. Curr Top Microbiol Immunol,2002,270:1-21.
[5] 张磊.两栖类TLR基因家族分子进化研究[D].曲阜:曲阜师范大学,2019.
ZHANG L.Molecular evolution of amphibian TLR gene family[D].Qufu:Qufu Normal University,2019.(in Chinese)
[6] ZHANG L,LIU G S,XIA T,et al.Evolution of toll-like receptor gene family in amphibians[J].Int J Biol Macromol, 2022, 208:463-474.
[7] 吴昌华,陈伟东,杨文政,等.Toll样受体7/8在公猪生殖系统中的表达及其配体对猪X/Y精子分选效果评价[J].畜牧兽医学报, 2023,54(4):1490-1499.
WU C H,CHEN W D,YANG W Z,et al.Expression of Toll-like receptors 7/8 in boar reproductive system and evaluation of effect of their ligands on boar X/Y sperm sorting[J].Acta Veterinaria et Zootechnica Sinica,2023,54(4):1490-1499.(in Chinese)
[8] 甘心梦,陈凯文,刘贺贺,等.TLRs基因家族在鸭法氏囊中的表达模式及分子进化分析[J].畜牧兽医学报,2018,49(9):1851-1860.
GAN X M,CHEN K W,LIU H H,et al.The expression patterns of TLRs gene family during duck bursa of fabricius development and their molecular evolutionary analysis[J].Acta Veterinaria et Zootechnica Sinica,2018,49(9):1851-1860.(in Chinese)
[9] ROACH J C,GLUSMAN G,ROWEN L,et al.The evolution of vertebrate Toll-like receptors[J].Proc Natl Acad Sci U S A,2005, 102(27):9577-9582.
[10] LIU G S,ZHANG H X,ZHAO C,et al.Evolutionary history of the Toll-Like receptor gene family across vertebrates[J].Genome Biol Evol,2020,12(1):3615-3634.
[11] 徐燕.山羊TOLL样受体基因家族的克隆表达及群体进化分析研究[D].武汉:华中农业大学,2009.
XU Y.Molecular characterization,expression and phylogenetic analysis of toll-like receptor gene family in goats[D].Wuhan:Huazhong Agricultural University,2009.(in Chinese)
[12] RAJA A,VIGNESH A R,MARY B A,et al.Sequence analysis of Toll-like receptor genes 1-10 of goat (Capra hircus)[J].Vet Immunol Immunopathol,2011,140(3-4):252-258.
[13] YU D D,WU Y,XU L,et al.Identification and characterization of toll-like receptors (TLRs) in the Chinese tree shrew (Tupaia belangeri chinensis)[J].Dev Comp Immunol,2016,60:127-138.
[14] CHEN L,QIU Q,JIANG Y,et al.Large-scale ruminant genome sequencing provides insights into their evolution and distinct traits[J].Science,2019,364(6446):eaav6202.
[15] CHEN S Q,LIU S L,SHI S L,et al.Comparative epigenomics reveals the impact of ruminant-specific regulatory elements on complex traits[J].BMC Biol,2022,20(1):273.
[16] HERNÁNDEZ FERNÁNDEZ M,VRBA E S.A complete estimate of the phylogenetic relationships in ruminantia:a dated species-level supertree of the extant ruminants[J].Biol Rev Camb Philos Soc,2005,80(2):269-302.
[17] 马志杰,王世康,张卫忠,等.柴达木牛和蒙古牛Y染色体基因组遗传多样性与父系起源研究[J].畜牧兽医学报,2023,54(3):1026-1033.
MA Z J,WANG S K,ZHANG W Z,et al.Y-chromosome genomic diversity and paternal origins of qaidam cattle and Mongolian cattle[J].Acta Veterinaria et Zootechnica Sinica,2023,54(3):1026-1033.(in Chinese)
[18] LAN D L,LIN B S,XIONG X R,et al.Identification and characteristics analysis of toll-like receptors family genes in yak[J]. Genes Genom,2016,38(5):429-438.
[19] DUBEY P K,GOYAL S,KATHIRAVAN P,et al.Sequence characterization of river buffalo Toll-like receptor genes 1-10 reveals distinct relationship with cattle and sheep[J].Int J Immunogenet,2013,40(2):140-148.
[20] YOU Z Y,LU B G,DU B B,et al.Spatio-temporal niche of sympatric tufted deer (Elaphodus cephalophus) and Sambar (Rusa unicolor) based on camera traps in the Gongga Mountain National Nature Reserve,China[J].Animals (Basel),2022,12(19):2694.
[21] MENZIES M,INGHAM A.Identification and expression of Toll-like receptors 1-10 in selected bovine and ovine tissues[J].Vet Immunol Immunopathol,2006,109(1-2):23-30.
[22] OGOREVC J,SIM AČG I AČG M,ZORC M,et al.TLR2 polymorphism (rs650082970) is associated with somatic cell count in goat milk[J]. Peerj,2019,7:7340.
[23] WANG M Y,QI Y,CAO Y T,et al.Domain fusion TLR2-4 enhances the autophagy-dependent clearance of Staphylococcus aureus in the genetic engineering goat[J].eLife,2022,11:e78044.
[24] 顾棚.偶蹄目TLR1-4基因的分子进化研究[D].雅安:四川农业大学,2014.
GU P.The study of molecular evolution on TLR1-4 gene of artiodactyla[D].Yaan:Sichuan Agricultural University,2014.(in Chinese)
[25] CHEN C J,CHEN H,ZHANG Y,et al.TBtools:an integrative toolkit developed for interactive analyses of big biological data[J].Mol Plant,2020,13(8):1194-1202.
[26] KUMAR S,STECHER G,LI M,et al.MEGA X:molecular evolutionary genetics analysis across computing platforms[J].Mol Biol Evol,2018,35(6):1547-1549.
[27] MACHADO D J,DE LUNA MARQUES F P,JIMÉNEZ-FERBANS L,et al.An empirical test of the relationship between the bootstrap and likelihood ratio support in maximum likelihood phylogenetic analysis[J].Cladistics,2022,38(3):392-401.
[28] NASER-KHDOUR S,MINH B Q,LANFEAR R.Assessing confidence in root placement on phylogenies:an empirical study using nonreversible models for mammals[J].Syst Biol,2022,71(4):959-972.
[29] R Core Team.R:a language and environment for statistical computing.R foundation for statistical computing,Vienna[Z].2021.
[30] ROYSTON J P.The W test for normality[J].Appl Stat,1982,31(2):176-180.
[31] WANG K,LI Y,ZHAI X W,et al.A method for evaluating the coal spontaneous combustion index by the coefficient of variation and Kruskal-Wallis test:a case study[J].Environ Sci Pollut Res Int,2023,30(20):58956-58966.
[32] GOYAL S,DUBEY P K,KUMARI N,et al.Caprine Toll-like receptor 8 gene sequence characterization reveals close relationships among ruminant species[J].Int J Immunogenet,2014,41(1):81-89.
[33] MCGUIRE K,JONES M,WERLING D,et al.Radiation hybrid mapping of all 10 characterized bovine Toll-like receptors[J].Anim Genet,2006,37(1):47-50.
[34] SARTORIUS R,TROVATO M,MANCO R,et al.Exploiting viral sensing mediated by Toll-like receptors to design innovative vaccines[J].npj Vaccines,2021,6(1):127.
[35] WERLING D,JANN O C,OFFORD V,et al.Variation matters:TLR structure and species-specific pathogen recognition[J].Trends Immunol,2009,30(3):124-130.
[36] 梁春鸿,李丹蕾,朱英,等.鸟类TLR基因进化分析[J].西南民族大学学报(自然科学版),2022,48(3):260-266.
LIANG C H,LI D L,ZHU Y,et al.Evolution of toll-like receptors in avian[J].Journal of Southwest University for Nationalities (Natural Science Edition),2022,48(3):260-266.(in Chinese)
[37] 周建波.大弹涂鱼Toll样受体基因的分子克隆、免疫应答及分子进化分析[D].舟山:浙江海洋大学,2019.
ZHOU J B.Molecular cloning,molecular evolution and expression analysis of 5 TLR genes in great blue-spotted mudskipper(Boleophthalmus pectinirostris)[D].Zhoushan:Zhejiang Ocean University,2019.(in Chinese)
[38] KATO S,OKAMURA E,MATSUNAGA T M,et al.Cyanidioschyzon merolae aurora kinase phosphorylates evolutionarily conserved sites on its target to regulate mitochondrial division[J].Commun Biol,2019,2:477.
[39] HUGHES J F,SKALETSKY H,PAGE D C.Sequencing of rhesus macaque Y chromosome clarifies origins and evolution of the DAZ (Deleted in AZoospermia) genes[J].Bioessays,2012,34(12):1035-1044.
[40] DENG W W,SU Z,LIANG P P,et al.Single-cell immune checkpoint landscape of PBMCs stimulated with Candida albicans[J].Emerg Microbes Infect,2021,10(1):1272-1283.