Correlation between Single Nucleotide Polymorphism in 5' Regulation Region of IL-8 and Coccidiosis-Resistance of Jinghai Yellow Chicken
-
摘要:目的 研究白介素8（IL-8）基因5′调控区单核苷酸突变对鸡柔嫩艾美耳球虫（E. tenella）抗性指标的影响。方法 本试验利用DNA直接测序技术，检测京海黄鸡IL-8基因5′调控区单核苷酸多态性（SNPs），并对5′调控区的SNPs突变前后的转录因子进行预测，然后分析SNPs与柔嫩艾美耳球虫抗性指标的关联性。结果 测序结果表明：IL-8基因5′调控区共检测到3个突变位点（T-550C、G-398T和T-360C），均形成了3种基因型，杂合度在0.436~0.471，PIC值在0.25~0.50，均属于中度多态，且3个突变位点均处于哈代-温伯格平衡状态。生物信息学分析表明：3个突变位点均改变了其原有的转录因子结合位点。关联分析显示：T-550C突变位点的TC型个体的IL-8表达量与TT型差异显著，TC型GSH-PX、CAT、IL-2、IL-6、IFN-γ指标均高于其他2种基因型，但差异不显著。G-398T突变位点的TT型个体SOD活性和GT型个体CAT活性均显著高于GG型；TT型个体NO含量与GG型差异极显著，与GT型差异显著；TT型IL-2表达量显著高于GT型。T-360C突变位点的TT型和TC型个体SOD活性与CC型差异极显著或差异显著；TT型和TC型个体的NO含量均显著高于CC型；TT型个体的IL-2、IL-8表达量均显著高于CC型。结论 G-398T和T-360C突变位点的TT型个体球虫抗性优于其他基因型个体，T-550C突变位点的杂合型个体球虫抗性优于纯合型个体，表明IL-8基因5′调控区的突变对球虫抗性指标有显著的调控作用，可作为球虫病抗性鸡新品种或新品系的育种参考依据。Abstract:Objective Correlation between the single nucleotide polymorphism (SNPs) in the 5' regulation region of IL-8 and the disease resistance to Eimeria tenella of chickens was studied.Method Direct DNA sequencing was performed to determine SNPs in IL-8 gene of Jinghai yellow chicken. Predicted transcription factors before and after mutation were analyzed to correlate between SNPs and coccidiosis-resistance of the birds.Result There were 3 SNPs at the mutation sites of T-550C, G-398T and T-360C detected in the target region that formed CT-, TT- and GG-genotypes with heterozygous degrees between 0.436 and 0.471 and PIC between 0.25 and 0.5. In a Hardy-Weinberg equilibrium state, SNPs were moderately polymorphic. A bioinformatics analysis confirmed that the original transcription factor binding sites in the genes had all been altered after the mutation. The expression of IL-8 of TC-genotype mutated at T-550C was significantly higher than that of TT-genotype. The coccidiosis-resistance indicators including GSH-PX, CAT, IL-2, IL-6 and IFN-γ in TC-genotype were higher than that in the other two genotypes, though not significantly. In genotypes mutated at G-398T, the SOD activity of TT-genotype and the CAT activity of GT-genotype were significantly higher than those of GG-genotype; the NO content of TT-genotype extremely significantly different from that of GG-genotype and significantly different from that of GT-genotype; and, the IL-2 expression of TT-genotype significantly higher than that of GT-genotype. Whereas, in the genotypes mutated at T-360C, the SOD activities of TT- and TC-genotypes were extremely significantly different or significantly different from that of CC-genotype; the NO contents of TT- and TC-genotypes significantly higher than that of the CC genotype; and, the IL-2
and IL-8 expressions of the TT genotype significantly higher than those of CC-genotypes. Conclusion It appeared that the TT-genotype with mutated G-398T or T-360C site would be more resistant to the parasitic attack by E. tenella than the other genotypes, and that hybrids of the genotypes carrying T-550C mutation would be more resistant than the homozygous counterparts. Thus, the significant regulation function of the polymorphism in the 5' regulation region of IL-8 gene could conceivably be targeted for breeding coccidiosis-resistant chickens. -
0. 引言
【研究意义】球虫病是由原生动物寄生虫引起的一种重要的肠道疾病，可导致家畜体重下降、营养不良、失血、脱水和对其他疾病因子的易感性增加，降低动物生产性能，给养殖业带来巨大的经济损失[1-2]。鸡对艾美耳球虫属易感，寄生于盲肠的柔嫩艾美耳球虫感染尤为严重，病理表现为肠道黏膜损伤、盲肠肿胀、出血及排出血便[3-4]。迄今为止，药物防治和活疫苗是两种主要的球虫病控制策略[5]，但抗球虫药和接种疫苗会引起鸡只的药物残留、抗药性及虫株变异等风险[6]。因此在分子水平上研究球虫病抗性基因，通过遗传选育、培育抗病品种是解决这一问题的有效途径。【前人研究进展】白介素8（IL-8）是一种趋化性细胞因子，一方面能激活嗜中性粒细胞，诱导中性粒细胞释放溶酶体酶，清除病原菌[7]，另一方面参与免疫反应，实现免疫细胞的抵抗功能，促进伤口愈合[8]。Cornelissen等[9]研究表明IL-8能够有效地招募Th1 CD4+、巨噬细胞和单核细胞，进而诱导IFN-γ的产生。Swaggerty等[10]通过连续3个世代选择外周血白细胞促炎性细胞因子白介素6（IL-6）、IL-8和趋化因子CCLi2的高水平表达（非球虫感染）的肉鸡个体进行配种，发现其后代柔嫩艾美耳球虫的抵抗能力极显著高于低表达水平群体。陈仁金等[11]采用混合动物模型将IL-8基因遗传多态性与荷斯坦奶牛7个性状进行关联分析，发现IL8基因连锁突变对荷斯坦牛泌乳性状以及抗病性状有较大的遗传效应，可用于中国荷斯坦牛的分子标记辅助选择。林雨鑫等[12]采用RNA-seq技术对E. tenella感染鸡和非感染对照鸡的盲肠组织进行转录组分析，筛选到许多显著富集的通路和差异表达基因，主要有IL-6、IL-8、IL-12β、IL-15、IL-17和TGFB2等。辛世杰等[8]研究发现，IL-8基因在感染柔嫩艾美耳球虫的京海黄鸡脾脏和盲肠组织中表达量显著升高，表明IL-8基因在鸡球虫感染中发挥一定的作用。【本研究切入点】虽然已有研究发现IL-8参与免疫应答、具有诱导中性粒细胞释放溶酶体酶和清除病原体等作用，但关于IL-8基因5′调控区单核苷酸多态性对球虫抗性指标的影响的研究鲜见报道。【拟解决的关键问题】本试验采用DNA直接测序技术，检测京海黄鸡IL-8基因5′调控区单核苷酸多态性（SNPs），并对5′调控区的SNPs突变前后转录因子进行预测，同时还分析这些SNPs与球虫抗性指标的关联性，旨在筛选与京海黄鸡柔嫩艾美耳球抗性选育有利的突变基因型，为培育球虫病抗性新品种或新品系提供参考依据。
1. 材料与方法
1.1 试验动物
选购江苏省海门集团1日龄京海黄鸡共92只，公母各半，饲养于无病原的鸡舍中，所食饲料中不含任何球虫及抗球虫性药物，试验鸡预试期内经粪检无球虫。E. tenella孢子化卵囊为扬州大学兽医学院寄生虫教研室保存虫种，经非球虫免疫健康鸡体内传代一次后收集获得。
1.2 样本处理及采集
30日龄时，每只鸡灌饲1.5万个柔嫩艾美耳球虫孢子化卵囊。记录感染期间0～8 d体增重。在感染后的第8 d翅静脉采血，肝素钠抗凝，立即4 000 r·min−1离心5 min，分离血浆和血细胞后于−20℃冰箱中保存。
1.3 基因组DNA提取及检测
根据文献[8]的方法提取所有试验鸡的全血基因组DNA，用微量紫外可见光分光光度计对基因组DNA进行OD值以及浓度测定，经测定所提取的DNA样品OD260 nm/OD280 nm值在1.8～2.0，质量浓度在500～1 000 μg·μL−1，均符合试验要求。将合格样品稀释至100 ng·μL−1，4℃或−20℃（长期）保存备用。
1.4 引物设计及PCR扩增
根据GenBank中已公布的鸡IL-8基因序列（Gene ID: 396495），使用Primer Premier 5.0对IL-8基因上游−2 000 bp至−1 bp的核苷酸序列共设计4对引物。引物由上海生工生物工程有限公司合成，引物设计详细信息见表1。
表 1 引物序列信息Table 1. Information on primer sequence引物 Primer 引物序列（5′→3′）Primers sequence(5′→3′) 退火温度 Annealing temperature/℃ 长度 Length/bp P1 F: TTCCATTCGCATAAGTCATC 51 638 R: AAAGTTGATTTGGGGATACC P2 F: TGTAATTGGGAATTCAAGGGGGA 58 708 R: CCCATTTGGTGTGTGATAAGATGA P3 F: AGTCCACAGACCACAAAGCA 58 693 R: TCGCAATATAAGTTTCTGATGGCTT P4 F: AAACCAGCAACACAAAGTC 60 574 R: CATCTCAGCAAGTGCCAAG PCR反应体系：10 μL的2×Taq Master Mix，上下游引物各1 μL，DNA模板1 μL，最后加入水补足至20 μL体系。PCR扩增程序：95℃预变性5 min；95℃变性30 s，引物最佳退火温度退火30 s，72℃ 延伸1 min，35个循环；72℃延伸10 min，4℃保存。
1.5 基因测序
PCR产物送至上海生物工程技术服务有限公司进行测序。使用仪器为ABI测序仪（3730xlDNAAnalyzer）。DNA序列拼接后，用MEGA6.06和DNAMAN5.2软件对测序结果进行比对分析，并确定京海黄鸡IL-8基因的SNPs。
1.6 转录因子预测
使用在线软件AliBaba2.1（http://gene-regulation.com/pub/programs/alibaba2/index.html）进行预测IL-8基因5′调控区核苷酸序列突变后可能存在的转录因子结合位点和可能结合的转录因子。
1.7 抗性指标测定
血浆抗氧化指标按试剂盒（购自南京建成生物工程研究所）说明书采用生物素双抗体夹心酶联免疫吸附法（ELISA）法测定[13-14]，测定指标包括：超氧化物歧化酶（SOD）、丙二醛（MDA）、谷胱甘肽过氧化物酶（GSH-PX）、过氧化氢酶（CAT）、一氧化氮（NO）、白介素1-β（IL-1β）、白介素2（IL-2）、白介素6（IL-6）、白介素8（IL-8）、干扰素-γ（IFN-γ）。
1.8 数据统计分析
使用Popgene1.32软件对IL-8基因5′调控区SNPs数据进行群体遗传多样性分析，包括杂合度（H）、有效等位基因数（Ne）和多态信息含量（PIC），统计整理各突变位点的基因型频率、等位基因频率，用χ2检验方法对群体Hardy-Weinberg进行平衡检测。运用SPSS25.0对各抗性指标数据进行关联分析。
2. 结果与分析
2.1 IL-8基因PCR测序结果
以提取的92只京海黄鸡基因组DNA为模板，依据NCBI数据库中鸡（Gallus gallus）IL-8基因上游−2 000 bp至−1 bp的核苷酸序列设计的引物进行PCR扩增测序及拼接。测序共发现了3个SNPs突变位点（见表2），1号突变位点在−550碱基发生T→C的突变，命名为T-550C。2号突变位点在−398碱基发生G→T的突变，命名为G-398T。3号突变位点在−360碱基发生T→C的突变，命名为T-360C。详细信息见图1。
![]()
图 1 IL-8基因位点位点峰图及序列比对注：A：T-550C位点峰图比对；B：G-398T位点峰图比对；C：T-360C位点峰图比对；D：T-550C位点序列比对图；E：G-398T位点序列比对图；F：T-360C位点序列比对图Figure 1. Peak and sequence alignments of IL-8Note： A: T-550C peak alignment; B: G-398T peak alignment; C: T-360C peak alignment; D: T-550C sequence alignment; E: G-398T sequence alignment; F: T-360C sequence alignment.表 2 IL-8基因5′调控区单核苷酸突变位点信息Table 2. Information on SNPs in 5′ regulation region of IL-8序号 Site number 染色体位置 Chromosome position 序列号 Serial number 单核苷酸多态性 SNP 1 51282560 rs740065165 T-550C 2 51282712 rs731947764 G-398T 3 51282750 rs16409254 T-360C 注：A：T-550C位点峰图比对；B：G-398T位点峰图比对；C：T-360C位点峰图比对；D：T-550C位点序列比对图；E：G-398T位点序列比对图；F：T-360C位点序列比对图
Note： A: T-550C peak alignment; B: G-398T peak alignment; C: T-360C peak alignment; D: T-550C sequence alignment; E: G-398T sequence alignment; F: T-360C sequence alignment.2.2 IL-8基因5′调控区SNPs突变引起的转录因子改变预测结果
由表3可见，IL-8基因5′调控区在−550 bp发生了T→C的突变，此处突变导致原有的Oct-1转录因子结合位点消失；在−398 bp发生了G→T的突变，导致此处原先C/EBPalp转录因子的结合位点发生位置改变，且新增了1个Pit-1a转录因子的结合位点；在−360 bp的突变（T→C）导致原有的NF-1转录因子结合位点消失。
表 3 IL-8基因5′调控区SNPs突变前后转录因子变化预测结果Table 3. Predicted transcription factors before and after mutation of IL-8 in 5' regulation region突变位点
Mutation site碱基
Base转录因子
Transcription factor转录因子结合位点碱基序列
Transcription factor binding site base sequence转录因子位置
Transcription factor position−550 bp T Oct-1 GTTGCATTTG −551~−542 bp C −398 bp G C/EBPalp GAAATAAATA −398~−389 bp T Pit-1a TAAATAAATA −398~−389 bp C/EBPalp ACATAAATAA −401~−392 bp −360 bp T NF-1 AGCCAGTTAT −362~−353 bp C 注：图中下划线标注为突变碱基。
Note: The underline in the figure is the mutant base.2.3 IL-8基因5′调控区多态性分析
对IL-8基因5′调控区发现的3个SNPs位点进行多态性及遗传多样性分析（表4），结果表明：3个突变位点均形成了3种基因型，杂合度均在0.436～0.471，PIC值在0.25～0.50，均属于中度多态，遗传差异性较大。卡方检验结果显示，3个突变位点均处于哈代-温伯格平衡状态。其中，除了G-398T突变位点外，原等位基因相对于突变后的等位基因均为优势等位基因，而G-398T突变位点突变后的等位基因为优势等位基因。
表 4 IL-8基因5′调控区多态性Table 4. SNPs in 5′ regulation region of IL-8突变位点
Mutation site基因型
Genotype数量
Number基因型频率
Genotype frequency等位基因
Allelic gene等位基因频率
Allele frequencyχ2值
χ2 valueP值
P value杂合度 H 有效等位
基因数 Ne多态信息含量
PICT-550C TT 40 0.435 T 0.679 1.390 0.500 0.436 1.772 0.341 CC 7 0.076 C 0.321 TC 45 0.489 G-398T GG 13 0.141 G 0.380 0.023 0.990 0.471 1.892 0.360 TT 35 0.381 T 0.720 GT 44 0.478 T-360C TT 40 0.435 T 0.674 0.725 0.699 0.440 1.784 0.343 CC 8 0.087 C 0.326 TC 44 0.478 注：PIC＞0.50为高度多态；0.25＜PIC＜0.50为中度多态；PIC＜0.25为低度多态[15]。
Note: PIC＞0.50 is highly polymorphic; 0.25＜PIC＜0.50 is moderately polymorphic; PIC＜0.25 means low polymorphism[15].2.4 IL-8基因5′调控区多态性与血浆抗性指标的关联分析
2.4.1 T-550C突变位点各基因型与鸡球虫抗性指标关联分析
T-550C突变位点各基因型与鸡球虫抗性指标关联分析显示（表5），除了SOD、MDA、NO和IL-1β指标外，TC型个体的其他指标均高于其他2种基因型，但差异不显著（P＞0.05）；TC型的IL-8表达量显著高于TT型（P＜0.05），其加性效应值和显性效应值分别为9.10和19.15 ng·L−1；TT型个体的NO含量极显著高于CC型（P＜0.01），其加性效应值和显性效应值分别为10.12和7.39 μmol·L−1。
表 5 IL-8基因T-550C突变位点各基因型与鸡球虫抗性指标关联分析Table 5. Correlation between genotypes mutated at T-550C in IL-8 and coccidiosis-resistance indicators指标
Indices基因型 Genotype TT（40） TC（45） CC（7） 超氧化物歧化酶（SOD）/（U·mL−1） 750.63±88.37 720.62±78.91 723.95±67.57 丙二醛（MDA）/（nmol·mL−1） 3.66±0.40 3.74±1.26 4.16±1.03 谷胱甘肽过氧化物酶（GSH-PX）/（U·mL−1） 151.77±16.94 152.22±17.53 150.83±22.52 过氧化氢酶（CAT）/（U·mL−1） 0.87±0.17 1.01±0.16 0.95±0.19 一氧化氮（NO）/（μmol·L−1） 52.35±5.25 A 49.62±3.79 AB 32.12±3.69 B 白介素1-β（IL-1β）/（ng·L−1） 20.54±2.00 19.75±3.04 20.22±1.19 白介素2（IL-2）/（ng·L−1） 14.78±2.66 15.82±2.04 15.18±2.69 白介素6（IL-6）/（ng·L−1） 69.93±3.55 72.40±8.64 67.81±4.29 白介素8（IL-8）/（ng·L−1） 90.87±7.63 b 119.11±10.42 a 109.06±11.76 ab 干扰素-γ（IFN-γ）/（ng·mL−1） 189.47±15.25 201.21±13.67 195.58±16.06 注：同行数据后不同大写字母表示差异极显著（P＜0.01），不同小写字母表示差异显著（P＜0.05），相同字母或无字母表示差异不显著（P＞0.05），表6～7同。
Note: Different uppercase letters on a same row indicate significant differences (P<0.01); different lowercase letters, significant differences (P<0.05); and, same or no letter, no significant differences (P>0.05). Same for Tables 6–7.2.4.2 G-398T突变位点各基因型与鸡球虫抗性指标关联分析
由表6可知，TT型个体SOD活性和GT型个体CAT活性均显著高于GG型（P＜0.05），SOD与CAT活性的加性效应值分别为46.77和0.035 U·mL−1，显性效应值分别为10.99和0.175 U·mL−1；TT型个体NO含量均高于其他两种基因型，且与GG型差异极显著（P＜0.01），与GT型差异显著（P＜0.05），其加性效应值和显性效应值分别为10.86和2.365 μmol·L−1；TT型IL-2表达量显著高于GT型（P＜0.05），其加性效应值和显性效应值分别为0.67和2.18 ng·L−1；TT型IL-1β、IL-8表达量均高于GG型个体，但未达到显著水平（P＞0.05）。
表 6 IL-8基因G-398T突变位点各基因型与鸡球虫抗性指标关联分析Table 6. Correlation between genotypes mutated at G-398T in IL-8 and coccidiosis-resistance indicators指标
Indices基因型 Genotype GG（13） GT（44） TT（35） 超氧化物歧化酶（SOD）/（U·mL−1） 631.43±50.73 b 667.21±27.41 ab 724.97±87.91 a 丙二醛（MDA）/（nmol·mL−1） 3.66±0.40 3.63±0.39 3.19±0.73 谷胱甘肽过氧化物酶（GSH-PX）/（U·mL−1） 164.04±21.41 155.20±20.63 158.58±19.50 过氧化氢酶（CAT）/（U·mL−1） 0.84±0.10 b 1.05±0.15 a 0.91±0.18 一氧化氮（NO）/（μmol·L−1） 24.18±4.26 Bab 32.67±9.88 ABb 45.89±11.63 Aa 白介素1-β（IL-1β）/（ng·L−1） 20.45±1.70 22.51±1.30 23.17±1.65 白介素2（IL-2）/（ng·L−1） 15.17±2.54 ab 13.66±1.38 b 16.51±3.10 a 白介素6（IL-6）/（ng·L−1） 65.66±6.05 71.01±7.81 68.40±7.30 白介素8（IL-8）/（ng·L−1） 94.31±9.86 104.24±7.70 105.89±13.30 干扰素-γ（IFN-γ）/（ng·mL−1） 183.25±17.99 179.94±15.83 175.99±17.26 2.4.3 T-360C突变位点各基因型与鸡球虫抗性指标关联分析
由T-360C突变位点各基因型与鸡球虫抗性指标关联分析（表7）可知，除MDA、CAT、IL-1β、IFN-γ指标以外，TT型个体的其他指标均高于其他2种基因型；对于SOD活性，TT型个体与CC型差异极显著（P＜0.01），TC型个体与CC型个体差异显著（P＜0.05），其加性效应值和显性效应值分别为64.43和10.55 U·mL −1；对于NO含量，TT型和TC型个体均显著高于CC型个体（P＜0.05），其加性效应值和显性效应值分别为5.94和4.53 μmol·L−1；对于IL-2、IL-8表达量，TT型个体均显著高于CC型个体（P＜0.05），IL-2与IL-8活性的加性效应值分别为1.68和10.02 ng·L−1，显性效应值分别为0.32和1.87 ng·L−1；对于IFN-γ含量，CC型个体显著高于其他2种基因型个体（P＜0.05），其加性效应值和显性效应值分别为15.015和9.55 ng·L−1；TT型GSH-PX、IL-6指标均高于其他2种基因型，但差异不显著（P＞0.05）。
表 7 IL-8基因T-360C突变位点各基因型与鸡球虫抗性指标关联分析Table 7. Correlation between genotypes mutated at T-360C in IL-8 and coccidiosis-resistance indicators指标 Indices 基因型 Genotype TT（40） TC（44） CC（8） 超氧化物歧化酶（SOD）/（U·mL−1） 781.29±76.94 Aab 706.31±72.85 ABa 652.43±44.92 Bb 丙二醛（MDA）/（nmol·mL−1） 4.01±1.20 3.91±0.49 4.23±0.79 谷胱甘肽过氧化物酶（GSH-PX）/（U·mL−1） 124.34±11.73 117.17±13.46 116.68±14.53 过氧化氢酶（CAT）/（U·mL−1） 0.77±0.16 0.86±0.30 0.74±0.14 一氧化氮（NO）/（μmol·L−1） 46.91±11.24 a 45.50±9.28 a 35.03±8.11 b 白介素1-β（IL-1β）/（ng·L−1） 20.11±1.03 20.65±1.84 21.55±2.73 白介素2（IL-2）/（ng·L−1） 15.82±2.04 a 13.82±1.84 ab 12.46±1.55 b 白介素6（IL-6）/（ng·L−1） 62.74±10.21 60.80±9.47 57.92±9.38 白介素8（IL-8）/（ng·L−1） 103.19±14.94 a 91.30±13.82 ab 83.15±13.01 b 干扰素-γ（IFN-γ）/（ng·mL−1） 163.28±13.49 b 168.75±15.29 b 193.31±16.32 a 3. 讨论与结论
3.1 IL-8基因5′调控区的多态性检测
5′调控区是基因转录调控的核心区域，该区域突变形成的SNPs可改变原有的转录因子结合位点和转录因子，从而影响与该基因有关性状的表达。因此研究基因5′调控区SNPs及其与相关性状的关联分析是发掘新的分子遗传标记的重要途径之一[4, 15]。宁娟[16]研究表明结构蛋白N和非结构蛋白Nsp2在Marc-145细胞中表达能激活转录因子NF-κB，而活化的NF-κB进入核内与IL-8启动子上的NF-κB位点结合后，能诱导IL-8的表达。Wu等[17]体外试验结果表明，转录激活因子3（ATF3）的抑制可显著提高人支气管上皮细胞中促炎细胞因子IL-6和IL-8的表达。本研究在5′调控区共检测到3个SNP位点：T-550C、G-398T和T-360C。这3个突变位点均形成了3种基因型，杂合度均在0.436～0.471，PIC值在0.25～0.50，均属于中度多态，遗传差异性较大，有利于鸡育种中的群体遗传分化。其中，T-550C和T-360C突变位点均是杂合型所占比例较高，TC为优势基因型，T为优势基因，而G-398T位点突变后的等位基因为优势等位基因，杂合型为优势基因型。卡方检验结果显示，3个突变位点均处于哈代-温伯格平衡状态，这可能是由于群体规模较大，突变位点没有受到选种选配的影响或是在选育过程中达到了新的平衡。生物信息学分析发现：在−550 bp发生了T→C的突变，导致原有的Oct-1转录因子结合位点消失；在−398 bp发生了G→T的突变，导致此处原先C/EBPalp转录因子结合位点发生位置改变，且新增了1个Pit-1a转录因子的结合位点；在−360 bp的突变（T→C）导致原有的NF-1转录因子结合位点消失。本研究发现的这3个突变位点均改变了其原有的转录因子结合位点，转录因子结合位点的改变可能使IL-8基因表达发生变化，进而影响血浆中IL-8的浓度。IL-8参与免疫反应，在炎症和免疫应答中起着重要作用。有研究表明鸡感染球虫后，体内IL-8表达量越高，其抗球虫能力越强[10]。因此推测转录因子结合位点的改变可能会影响鸡抗球虫能力，但具体机制有待进一步研究。
3.2 IL-8基因5′调控区的多态性与鸡球虫抗性指标的关系
目前有关鸡IL-8基因的SNPs与球虫抗性指标关联分析的研究较少。本研究将5′调控区发现的3个SNPs与鸡球虫抗性指标关联分析，结果表明：在T-550C位点突变后形成的3种基因型中，大多数指标均是TC型高于其他2种基因型，且TC型个体IL-8表达量显著高于TT型。对于G-398T突变位点，TT型个体显著或极显著高于GG型的SOD、NO指标，显著高于GT型的IL-2表达量，并且TT型IL-1β、IL-8表达量也都高于GG型个体，因此突变后的TT型个体可能抵抗球虫病能力较强。对于T-360C突变位点，除了对MDA、GSH-Px、CAT、IL-1β、IL-6指标无显著影响外，TT型和TC型个体的SOD活性极显著或显著高于CC型，TT型和TC型个体的NO含量均显著高于CC型个体，TT型个体的IL-2、IL-8表达量均显著高于CC型。IL-8参与免疫应答，能诱导中性粒细胞释放溶酶体酶和清除病原体，在炎症和免疫应答中发挥关键调节作用。本研究在IL-8基因5′调控区检测到的3个SNPs均改变了其原有的转录因子结合位点，进而影响IL-8基因的表达，导致血浆中IL-8的表达量升高，有效地提高鸡体内的免疫水平，增强抵抗E. tenella的能力。
-
![]()
图 1 IL-8基因位点位点峰图及序列比对
注：A：T-550C位点峰图比对；B：G-398T位点峰图比对；C：T-360C位点峰图比对；D：T-550C位点序列比对图；E：G-398T位点序列比对图；F：T-360C位点序列比对图
Figure 1. Peak and sequence alignments of IL-8
Note： A: T-550C peak alignment; B: G-398T peak alignment; C: T-360C peak alignment; D: T-550C sequence alignment; E: G-398T sequence alignment; F: T-360C sequence alignment.
表 1 引物序列信息
Table 1 Information on primer sequence
引物 Primer 引物序列（5′→3′）Primers sequence(5′→3′) 退火温度 Annealing temperature/℃ 长度 Length/bp P1 F: TTCCATTCGCATAAGTCATC 51 638 R: AAAGTTGATTTGGGGATACC P2 F: TGTAATTGGGAATTCAAGGGGGA 58 708 R: CCCATTTGGTGTGTGATAAGATGA P3 F: AGTCCACAGACCACAAAGCA 58 693 R: TCGCAATATAAGTTTCTGATGGCTT P4 F: AAACCAGCAACACAAAGTC 60 574 R: CATCTCAGCAAGTGCCAAG 
下载: 导出CSV
表 2 IL-8基因5′调控区单核苷酸突变位点信息
Table 2 Information on SNPs in 5′ regulation region of IL-8
序号 Site number 染色体位置 Chromosome position 序列号 Serial number 单核苷酸多态性 SNP 1 51282560 rs740065165 T-550C 2 51282712 rs731947764 G-398T 3 51282750 rs16409254 T-360C 注：A：T-550C位点峰图比对；B：G-398T位点峰图比对；C：T-360C位点峰图比对；D：T-550C位点序列比对图；E：G-398T位点序列比对图；F：T-360C位点序列比对图
Note： A: T-550C peak alignment; B: G-398T peak alignment; C: T-360C peak alignment; D: T-550C sequence alignment; E: G-398T sequence alignment; F: T-360C sequence alignment.
下载: 导出CSV
表 3 IL-8基因5′调控区SNPs突变前后转录因子变化预测结果
Table 3 Predicted transcription factors before and after mutation of IL-8 in 5' regulation region
突变位点
Mutation site碱基
Base转录因子
Transcription factor转录因子结合位点碱基序列
Transcription factor binding site base sequence转录因子位置
Transcription factor position−550 bp T Oct-1 GTTGCATTTG −551~−542 bp C −398 bp G C/EBPalp GAAATAAATA −398~−389 bp T Pit-1a TAAATAAATA −398~−389 bp C/EBPalp ACATAAATAA −401~−392 bp −360 bp T NF-1 AGCCAGTTAT −362~−353 bp C 注：图中下划线标注为突变碱基。
Note: The underline in the figure is the mutant base.
下载: 导出CSV
表 4 IL-8基因5′调控区多态性
Table 4 SNPs in 5′ regulation region of IL-8
突变位点
Mutation site基因型
Genotype数量
Number基因型频率
Genotype frequency等位基因
Allelic gene等位基因频率
Allele frequencyχ2值
χ2 valueP值
P value杂合度 H 有效等位
基因数 Ne多态信息含量
PICT-550C TT 40 0.435 T 0.679 1.390 0.500 0.436 1.772 0.341 CC 7 0.076 C 0.321 TC 45 0.489 G-398T GG 13 0.141 G 0.380 0.023 0.990 0.471 1.892 0.360 TT 35 0.381 T 0.720 GT 44 0.478 T-360C TT 40 0.435 T 0.674 0.725 0.699 0.440 1.784 0.343 CC 8 0.087 C 0.326 TC 44 0.478 注：PIC＞0.50为高度多态；0.25＜PIC＜0.50为中度多态；PIC＜0.25为低度多态[15]。
Note: PIC＞0.50 is highly polymorphic; 0.25＜PIC＜0.50 is moderately polymorphic; PIC＜0.25 means low polymorphism[15].
下载: 导出CSV
表 5 IL-8基因T-550C突变位点各基因型与鸡球虫抗性指标关联分析
Table 5 Correlation between genotypes mutated at T-550C in IL-8 and coccidiosis-resistance indicators
指标
Indices基因型 Genotype TT（40） TC（45） CC（7） 超氧化物歧化酶（SOD）/（U·mL−1） 750.63±88.37 720.62±78.91 723.95±67.57 丙二醛（MDA）/（nmol·mL−1） 3.66±0.40 3.74±1.26 4.16±1.03 谷胱甘肽过氧化物酶（GSH-PX）/（U·mL−1） 151.77±16.94 152.22±17.53 150.83±22.52 过氧化氢酶（CAT）/（U·mL−1） 0.87±0.17 1.01±0.16 0.95±0.19 一氧化氮（NO）/（μmol·L−1） 52.35±5.25 A 49.62±3.79 AB 32.12±3.69 B 白介素1-β（IL-1β）/（ng·L−1） 20.54±2.00 19.75±3.04 20.22±1.19 白介素2（IL-2）/（ng·L−1） 14.78±2.66 15.82±2.04 15.18±2.69 白介素6（IL-6）/（ng·L−1） 69.93±3.55 72.40±8.64 67.81±4.29 白介素8（IL-8）/（ng·L−1） 90.87±7.63 b 119.11±10.42 a 109.06±11.76 ab 干扰素-γ（IFN-γ）/（ng·mL−1） 189.47±15.25 201.21±13.67 195.58±16.06 注：同行数据后不同大写字母表示差异极显著（P＜0.01），不同小写字母表示差异显著（P＜0.05），相同字母或无字母表示差异不显著（P＞0.05），表6～7同。
Note: Different uppercase letters on a same row indicate significant differences (P<0.01); different lowercase letters, significant differences (P<0.05); and, same or no letter, no significant differences (P>0.05). Same for Tables 6–7.
下载: 导出CSV
表 6 IL-8基因G-398T突变位点各基因型与鸡球虫抗性指标关联分析
Table 6 Correlation between genotypes mutated at G-398T in IL-8 and coccidiosis-resistance indicators
指标
Indices基因型 Genotype GG（13） GT（44） TT（35） 超氧化物歧化酶（SOD）/（U·mL−1） 631.43±50.73 b 667.21±27.41 ab 724.97±87.91 a 丙二醛（MDA）/（nmol·mL−1） 3.66±0.40 3.63±0.39 3.19±0.73 谷胱甘肽过氧化物酶（GSH-PX）/（U·mL−1） 164.04±21.41 155.20±20.63 158.58±19.50 过氧化氢酶（CAT）/（U·mL−1） 0.84±0.10 b 1.05±0.15 a 0.91±0.18 一氧化氮（NO）/（μmol·L−1） 24.18±4.26 Bab 32.67±9.88 ABb 45.89±11.63 Aa 白介素1-β（IL-1β）/（ng·L−1） 20.45±1.70 22.51±1.30 23.17±1.65 白介素2（IL-2）/（ng·L−1） 15.17±2.54 ab 13.66±1.38 b 16.51±3.10 a 白介素6（IL-6）/（ng·L−1） 65.66±6.05 71.01±7.81 68.40±7.30 白介素8（IL-8）/（ng·L−1） 94.31±9.86 104.24±7.70 105.89±13.30 干扰素-γ（IFN-γ）/（ng·mL−1） 183.25±17.99 179.94±15.83 175.99±17.26
下载: 导出CSV
表 7 IL-8基因T-360C突变位点各基因型与鸡球虫抗性指标关联分析
Table 7 Correlation between genotypes mutated at T-360C in IL-8 and coccidiosis-resistance indicators
指标 Indices 基因型 Genotype TT（40） TC（44） CC（8） 超氧化物歧化酶（SOD）/（U·mL−1） 781.29±76.94 Aab 706.31±72.85 ABa 652.43±44.92 Bb 丙二醛（MDA）/（nmol·mL−1） 4.01±1.20 3.91±0.49 4.23±0.79 谷胱甘肽过氧化物酶（GSH-PX）/（U·mL−1） 124.34±11.73 117.17±13.46 116.68±14.53 过氧化氢酶（CAT）/（U·mL−1） 0.77±0.16 0.86±0.30 0.74±0.14 一氧化氮（NO）/（μmol·L−1） 46.91±11.24 a 45.50±9.28 a 35.03±8.11 b 白介素1-β（IL-1β）/（ng·L−1） 20.11±1.03 20.65±1.84 21.55±2.73 白介素2（IL-2）/（ng·L−1） 15.82±2.04 a 13.82±1.84 ab 12.46±1.55 b 白介素6（IL-6）/（ng·L−1） 62.74±10.21 60.80±9.47 57.92±9.38 白介素8（IL-8）/（ng·L−1） 103.19±14.94 a 91.30±13.82 ab 83.15±13.01 b 干扰素-γ（IFN-γ）/（ng·mL−1） 163.28±13.49 b 168.75±15.29 b 193.31±16.32 a
下载: 导出CSV
-
[1] HUANG G P, TANG X L, BI F F, et al. Eimeria tenella Infection Perturbs the Chicken Gut Microbiota from the Onset of Oocyst Shedding [J]. Veterinary Parasitology, 2018, 258: 30−37. DOI: 10.1016/j.vetpar.2018.06.005
[2] BLAKE D P, TOMLEY F M. Securing poultry production from the ever-present Eimeria challenge [J]. Trends in Parasitology, 2014, 30(1): 12−19. DOI: 10.1016/j.pt.2013.10.003
[3] 林雨鑫, 张菁菁, 戴国俊, 等. 京海黄鸡杂交配套系亲本对柔嫩艾美耳球虫的敏感性 [J]. 中国兽医学报, 2015, 35(7):1074−1078. LIN Y X, ZHANG J J, DAI G J, et al. Comparation of the sensitivity of E.tenella among Jinghai yellow chicken hybrid matching system [J]. Chinese Journal of Veterinary Science, 2015, 35(7): 1074−1078.(in Chinese)
[4] 王诗琴, 辛世杰, 王晓慧, 等. IL-6基因启动子区单核苷酸多态对京海黄鸡柔嫩艾美尔球虫抗性指标的影响 [J]. 扬州大学学报(农业与生命科学版), 2018, 39(1):36−41. WANG S Q, XIN S J, WANG X H, et al. The effects of IL-6 gene promoter region single nucleotide polymorphism to Eimeria tenella resistance index in Jinghai yellow chickens [J]. Journal of Yangzhou University(Agricultural and Life Science Edition), 2018, 39(1): 36−41.(in Chinese)
[5] JANG S I, LILLEHOJ H S, LEE S H, et al. Immunoenhancing effects of MontanideTM Isa oil-based adjuvants on recombinant Coccidia antigen vaccination against Eimeria acervulina infection [J]. Veterinary Parasitology, 2010, 172(3/4): 221−228.
[6] ABBAS R Z, IQBAL Z, BLAKE D, et al. Anticoccidial drug resistance in fowl Coccidia: the state of play revisited [J]. World's Poultry Science Journal, 2011, 67(2): 337−350. DOI: 10.1017/S004393391100033X
[7] COUSSENS L M, WERB Z. Inflammation and cancer [J]. Nature, 2002, 420(6917): 860−867. DOI: 10.1038/nature01322
[8] 辛世杰, 王晓慧, 戴国俊, 等. 京海黄鸡柔嫩艾美尔球虫感染对脾脏和盲肠IL-6、IL-8、CCLi2基因表达量的影响及其相关性 [J]. 浙江农业学报, 2019, 31(1):39−46. DOI: 10.3969/j.issn.1004-1524.2019.01.05 XIN S J, WANG X H, DAI G J, et al. Effect and correlation analysis of Eimeria tenella infection on IL-6, IL-8 and CCLi2 genes expression in spleen and Caecum of Jinghai Yellow Chicken(Gallus gallus) [J]. Acta Agriculturae Zhejiangensis, 2019, 31(1): 39−46.(in Chinese) DOI: 10.3969/j.issn.1004-1524.2019.01.05
[9] CORNELISSEN J B W J, SWINKELS W J C, BOERSMA W A, et al. Host response to simultaneous infections with Eimeria acervulina, maxima and tenella: A cumulation of single responses [J]. Veterinary Parasitology, 2009, 162(1/2): 58−66.
[10] SWAGGERTY C L, PEVZNER I Y, KOGUT M H. Selection for pro-inflammatory mediators produces chickens more resistant to Eimeria tenella [J]. Poultry Science, 2015, 94(1): 37−42. DOI: 10.3382/ps/peu053
[11] 陈仁金, 杨章平, 毛永江, 等. 中国荷斯坦牛IL-8基因遗传多态性与泌乳性状以及体细胞评分的关联 [J]. 遗传, 2010, 32(12):1256−1262. CHEN R J, YANG Z P, MAO Y J, et al. Genetic polymorphism of the IL-8 gene and its associations with milk traits and SCS in Chinese Holstein [J]. Hereditas, 2010, 32(12): 1256−1262.(in Chinese)
[12] 林雨鑫. 京海黄鸡柔嫩艾美耳球虫感染后盲肠转录组分析及其杂交配套系亲本的抗性评价[D]. 扬州: 扬州大学, 2015. LIN Y X. RNA Sequencing Analysis of Chicken Cecum Tissues Following E. tenella Infection and Coccidiosis Evaluation of Jinghai Yellow Chicken Cross-breeding System Parents[D]. Yangzhou: Yangzhou University, 2015.(in Chinese)
[13] 杨军. 益生菌对柔嫩艾美尔球虫感染雏鸡免疫水平的影响及免疫保护效果评价[D]. 长春: 吉林农业大学, 2015. YANG J. Studies on effect of probiotics on chicken immunity levels and the evaluation of immune protective effect on Eimeria tenella[D]. Changchun: Jilin Agricultural University, 2015. (in Chinese)
[14] 李玲秀. 茶多酚与茶皂素对隐性感染球虫的肉鸡生长、免疫力及抗氧化性能的影响[D]. 合肥: 安徽农业大学, 2014. LI L X. The influence of tea polyphenols and tea saponin on growth, immunity and antioxidant ability of broilers with Coccidium latent infection[D]. Hefei: Anhui Agricultural University, 2014. (in Chinese)
[15] 李秀秀, 王文文, 郝园丽, 等. 寿光鸡AANAT基因5′调控区的单核苷酸多态性及其与产蛋性状的关联分析 [J]. 农业生物技术学报, 2017, 25(6):921−929. LI X X, WANG W W, HAO Y L, et al. SNPs in 5′ regulatory region of AANAT gene and the association analysis with laying traits in Shouguang chickens(Gallus gallus) [J]. Journal of Agricultural Biotechnology, 2017, 25(6): 921−929.(in Chinese)
[16] 宁娟. 猪繁殖与呼吸综合征病毒感染Marc-145细胞激活白细胞介素8的表达的研究[D]. 武汉: 华中农业大学, 2009. NING J. Study on the Induction of Interleukin-8 by Porcine Reproductive and Respiratory Syndrome Virus in Marc-145 Cells[D]. Wuhan: Huazhong Agricultural University, 2009. (in Chinese)
[17] WU Y P, CAO C, WU Y F, et al. Activating transcription factor 3 represses cigarette smoke-induced IL6 and IL8 expression via suppressing NF-κB activation [J]. Toxicology Letters, 2017, 270: 17−24. DOI: 10.1016/j.toxlet.2017.02.002
-
期刊类型引用(7)
1. 符乐毅，陈智雄. 茶梗栽培食用菌研究进展. 现代农业科技. 2024(05): 53-56 . 
百度学术
2. 马婧嘉，李雅茹，雷萍，张文隽，吴亚召. 基于人工神经网络和遗传算法优化用于栽培皱环球盖菇的农业废弃物培养基. 菌物学报. 2024(11): 146-154 . 百度学术
3. 李杨，黄艺鸿，林银娜，王雅婷. 茶梗的综合利用研究进展. 农产品加工. 2023(02): 87-90 . 百度学术
4. 王爱仙，王怡暄，刘福阳，刘金仙，巫仁高. 茶枝屑代料栽培黑木耳配方. 北方园艺. 2021(12): 121-126 . 百度学术
5. 王伟伟，张建勇，陈琳，江和源. 茶梗的综合利用研究进展. 茶叶通讯. 2020(01): 20-24 . 百度学术
6. 刘金仙，邓家耀，傅仙玉，王爱仙. 茶梗、茶枝屑代料对黑木耳生长及品质的影响. 宜春学院学报. 2020(03): 108-113 . 百度学术
7. 王丽丽，杨军国，林清霞，项丽慧，陈林. 基于化学组分的铁观音茶梗和叶茶模式识别研究. 茶叶学报. 2019(01): 1-8 . 百度学术
其他类型引用(3)
计量
- 文章访问数: 1015
- HTML全文浏览量: 400
- PDF下载量: 17
- 被引次数: 10
