四川动物  2018, Vol. 37 Issue (5): 497-506

扩展功能

文章信息

朱未, 朱磊, 何凡, 付磊, 帅素蓉, 江建平
ZHU Wei, ZHU Lei, HE Fan, FU Lei, SHUAI Surong, JIANG Jianping
阿坝藏猪蒲江种群SLA-DRB基因片段SNP检测及其与生长发育性状关联分析
Single Nucleotide Polymorphisms of the Partial Sequence of SLA-DRB Gene in Pujiang Population of Tibetan Pigs and Correlation with Growth Traits
四川动物, 2018, 37(5): 497-506
Sichuan Journal of Zoology, 2018, 37(5): 497-506
10.11984/j.issn.1000-7083.20180102

文章历史

收稿日期: 2018-03-27
接受日期: 2018-05-25
阿坝藏猪蒲江种群SLA-DRB基因片段SNP检测及其与生长发育性状关联分析
朱未1 , 朱磊2 , 何凡3 , 付磊1 , 帅素蓉3 , 江建平1*     
1. 中国科学院成都生物研究所, 成都 610041
2. 四川中测环境技术有限公司, 成都 611731
3. 四川农业大学动物科技学院, 成都 611130
摘要:猪Sus scrofa domesticus的白细胞抗原复合体(SLA)基因的单核苷酸多态性(SNP)与猪生长发育等性状有关。为了开发藏猪品系选育的遗传标记,本文测定了阿坝藏猪蒲江种群89头个体的SLA-DRB基因部分序列,共鉴定出14个SNP位点。根据位点间的连锁强度分为3个组合(组合1:g.11470、g.11945、g.11997、g.11387、g.11369与g.11829;组合2:g.11532与g.11838;组合3:g.11911与g.11708)和4个独立多态位点(g.11468、g.11440、g.11478和g.11944)。将基因型与体质量、体长等14项生长发育性状进行关联性分析,结果表明,组合3的基因型与体高、体长、胸围等11项体尺指标间的相关性均显著(FDR < 0.05),基因型为TT/CG的个体较小,表明SLA-DRB基因可用于藏猪小型化选育的遗传标记。管围和12项管围率(除体质量以外)与SLA-DRB基因多态性的相关性最高,受组合1、组合2和位点g.11468、g.11478基因型的极显著影响(FDR < 0.001),其中,基因型为组合1=GGGAAC/GGGAAC、组合2=GT/GT、位点g.11468=C/C和g.11478=C/C的个体拥有更大的管围和管围率。管骨骨髓是SLA基因发挥免疫功能的主要场所之一。本研究为基于管围的猪种筛选提供了理论依据和可用的SNP位点,也为藏猪抗病遗传筛选提供了基础数据和参考资料。
关键词藏猪     SLA-DRB基因     管围     免疫功能     遗传筛选    
Single Nucleotide Polymorphisms of the Partial Sequence of SLA-DRB Gene in Pujiang Population of Tibetan Pigs and Correlation with Growth Traits
ZHU Wei1 , ZHU Lei2 , HE Fan3 , FU Lei1 , SHUAI Surong3 , JIANG Jianping1*     
1. Chinese Academy of Sciences, Chengdu Institute of Biology, Chengdu 610041, China;
2. Sichuan Zhongce Environmental Technology Co., Ltd., Chengdu 611731, China;
3. College of Animal Science, Sichuan Agricultural University, Chengdu 611130, China
Abstract: The single nucleotide polymorphism (SNP) of swine leukocyte antigen (SLA) gene is associated with the growth and development traits of swine. In this study, the partial sequence of SLA-DRB gene of 89 Aba Tibetan pigs (Pujiang county, Sichuan province) was sequenced to screen the stable markers for genetic breeding, and 14 loci of SNPs were identified. Based on the linkage relationships, these 14 loci of SNPs were allocated into 3 groups (loci group 1:g.11470, g.11945, g.11997, g.11387, g.11369 and g.11829; loci group 2:g.11532 and g.11838; loci group 3:g.11911 and g.11708) and 4 single polymorphic loci (g.11468, g.11440, g.11478 and g.11944). Correlation analyses were conducted between the genotypes of these loci or loci groups and 14 growth traits (i.e. body mass and body length). The polymorphism of loci group 3 was significantly associated with 11 growth traits (FDR < 0.05), including withers height, body length and chest circumference, and the individuals with genotype TT/CG had smaller somatotype than that of others. This suggested that SLA-DRB gene might be a stable genetic marker for the miniaturization selection of Tibetan pigs. The cannon bone circumference (CBC) and ratios of CBC to other growth traits (except body mass), were tightly associated with the polymorphism of loci group 1, loci group 2, g.11468 and g.11478 (FDR < 0.001). The individuals with genotypes GGGAAC/GGGAAC at loci group 1, GT/GT at loci group 2, C/C at g.11468 and C/C at g.11478 tend to have large CBC and CBC ratios. Given that the cannon bone marrow is a major site for the immune function of SLA genes, our results here provide theoretical basis and stable SNP loci for CBC based swine selective breeding and disease resistance oriented genetic breeding.
Keywords: Tibetan pig     SLA-DRB gene     cannon bone circumference     immune function     genetic breeding    

Sus scrofa domesticus在解剖学、生理学、疾病发生机理方面和人Homo sapiens极其相似(任丽华,2006),是人类疾病研究的潜力模型。小型猪的一般生物学特征与普通家猪基本相同,但其具有养殖成本低、便于实验操作、节约试剂与受试物等突出优势,目前已被广泛用于人类疾病模型构建、药理毒理学研究以及新药与化妆品评价等方面(Bronner et al., 2002Garthoff et al., 2002Ryska et al., 2004Sebert et al., 2005Turk et al., 2005Imai et al., 2006;袁进,顾为望,2011)。目前,我国小型猪资源和品种主要有版纳微型猪、五指山猪、广西巴马小型猪、贵州香猪、滇南小耳猪、藏猪等(李群,2008)。其中,藏猪是我国唯一的高原高寒型放牧猪种,环境适应力强、抗病性高、饲养成本低,6月龄体质量约20 kg(何凡等,2013),适用于药物代谢实验,因而具有重要的经济价值和科研价值。在此基础上,进一步选育抗病能力强、生长发育稳定、小型化的品系具有良好的应用前景,而筛选可以作为标记辅助选择(marker-assisted selection,MAS)的遗传标记无疑是实现这一目标的重要手段。

研究表明,猪的白细胞抗原复合体(swine leukocyte antigen,SLA)基因多态性不仅影响猪对疾病的抵抗能力,还与个体生长发育、繁殖性能、胴体品质以及肉质等性状有关(Gautschi & Gaillard,1990Mallard et al., 1991彭勇波等,2005)。SLA基因主要分为Ⅰ、Ⅱ和Ⅲ三大类,其中,SLA Ⅱ类基因位于SLA-D区,控制机体的免疫应答与调控,在调节机体抗病能力方面起到非常重要的作用。SLA Ⅱ类基因包括DPADPBDMADMBDOADOBDQADQBDRADRB等,但仅SLA-DQSLA-DR基因能编码蛋白质(Shia et al., 1991, 1995)。SLA-DRB基因包含6个外显子和5个内含子,其产物由237个氨基酸残基组成(Gustafsson et al., 1990a, 1990b)。吴圣龙等(2007)对苏太猪SLA-DQBSLA-DRB基因外显子多态性及其与繁殖性能的关联分析发现,组合基因型BBDF的第3胎繁殖性能(总产仔数、产活仔数、初生窝重和断奶仔猪数)显著高于组合基因型AADD、AADF、ABDD和CCDF。张冬杰(2010)对黑龙江野猪与北京黑猪杂交的F1代个体的研究发现,SLA-DRB基因经BcnⅠ和MboⅠ酶切后产生BB型和AB型的个体,其仔猪的初生质量更大。这些研究结果表明,SLA-DRB基因的多态性在猪种的选育方面具有重要参考价值。

本文检测了近亲选配的阿坝藏猪SLA-DRB基因部分内含子3、外显子4、内含子4、外显子5共1 220 bp序列,分析了其单核苷酸多态性(single nucleotide polymorphism,SNP)位点多态性及其与生长发育性状的相关性,以筛选可以用于藏猪早期筛选的MAS遗传标记。这将为藏猪的选育、资源评价、利用及保护提供理论依据,促进藏猪实验动物化的研发工作。

1 材料和方法 1.1 实验动物

本实验的藏猪是引自四川省阿坝藏族羌族自治州红原县的1雄2雌全同胞藏猪的后代(F4~F6),饲养于四川省蒲江县复兴镇。雌雄藏猪分圈饲养,每圈2头,按日龄供料,每日09:00和17:00各饲喂1次,自动饮水器喂水。

1.2 性状测定

共测定89头藏猪(45雄44雌)不同月龄的体高、体长、头长、最大额宽、头深、尻高(荐部最高点至地面的垂直距离)、尻长(腰角前缘到臀端后缘的直线距离)、腰角宽、胸深、胸宽、胸围、腹围、管围(左前肢前臂骨上三分之一处的水平周径)和体质量共14项指标,具体测量标准参考何凡(2013)。体质量精确到0.05 kg,其他指标精确到0.2 cm。所有被测猪体尺均在喂食2 h后实地测定,体质量在早上饲喂前称量。管围率为管围除以其他体尺指标后所得的衍生指标。

1.3 SLA-DRB基因序列获取

取约25 mg耳组织样品,切碎后置于无菌的1.5 mL离心管中,细胞组织DNA提取试剂盒(TransGen Biotech,北京)提取DNA,通用型DNA纯化回收试剂盒(天根生化科技,北京)纯化组织DNA。

参照SLA-DRB基因序列(GenBank登录号:AY303991),利用Primer Premier 5.0设计PCR扩增引物,主要扩增区域为该基因的部分内含子3、外显子4、内含子4、外显子5,扩增片段长度为1 220 bp。上游引物序列为SLA-DRB F:5'-ACCCACTGAGCG-AGGGAAGGGATTG-3',下游引物序列为SLA-DRB R:5'-GGAGGACGCAGAGCATAGCAGGAGC-3'。

PCR反应体系为25 μL,包括2×Mix 12 μL,上、下游引物各0.7 μL(10 pmol·μL-1),DNA模板1.6 μL(50 ng·μL-1),ddH2O 10 μL。PCR扩增条件为:94 ℃预变性4 min;94 ℃变性35 s,68 ℃退火30 s,72 ℃延伸40 s,35个循环;72 ℃延伸10 min,最后4 ℃保存。PCR产物用1.5%琼脂糖凝胶电泳检测,切胶后用试剂盒纯化,纯化产物送生工生物工程(上海)股份有限公司和北京六合华大基因科技有限公司测序。

1.4 数据分析

使用DNAStar中的MegAlign对测序获得的所有SLA-DRB基因序列进行比对,检测SNP位点。基于所有SLA-DRB基因序列计算其基因频率、基因型频率,并进行Hardy-Weinberg平衡卡方适合性检验(王铁岗,2009),计算位点纯合度(Ho)、位点杂合度(He)、有效等位基因(Ne)和多态信息含量(PIC)(Nei,1978)。参照Shi和He(2005)的方法,利用SHEsis(http://analysis2.bio-x.cn/myAnalysis.php)对SLA-DRB基因的SNP位点进行连锁不平衡分析,并根据连锁强度将各SNP位点分组,每组含1到多个SNP位点。对于仅含1个SNP的位点组(即独立SNP),利用Matlab的多因素方差分析评价该位点的单体型组合和日龄与各体尺指标的关联;对于含有2个及以上的SNP位点组,参照黄萌和蔡琳(2009)的方法,利用PHASE 2.1分析其在该藏猪群体中存在的单体型、单体型组合及各自频率,然后利用Matlab的多因素方差分析评价该位点组的单体型组合和日龄与各体尺指标的关联;合并两部分的分析结果,并用BHFDR对P值进行矫正。

2 结果 2.1 SNP位点的识别

共获得89头个体的SLA-DRB基因序列,通过序列比对发现存在14个SNP位点(图 1),其中,6个突变位于内含子3,分别为g.11369 C>G、g.11387 A>G、g.11440 C>T、g.11468 C>T、g.11470 A>G、g.11478 C>T;1个突变位于外显子4 (g.11532 A>G);7个突变位于内含子4,分别为g.11708 C>T、g.11829 A>G、g.11838 C>T、g.11911 T>G、g.11944 C>G、g.11945 A>G、g.11997 C>T;外显子5无突变。

图 1 藏猪SLA-DRB基因14个单核苷酸多态性位点(箭头所指) Fig. 1 14 loci of single nucleotide polymorphisms of SLA-DRB gene in Tibetan pigs (arrows)
2.2 SLA-DRB基因遗传结构分析

遗传结构分析结果(表 1)表明,位点g.11369、g.11387、g.11470、g.11478、g.11532、g.11829、g.11838、g.11944、g.11945、g.11997处于中度多态(0.25<PIC<0.50),多态信息含量较高,遗传变异较大,可以获得更多的选择结果,暗示这些位点的整齐度还存在进一步提高的潜能;位点g.11440、g.11468、g.11708、g.11911处于低度多态(PIC<0.25),多态信息含量较低,说明这些位点的遗传变异较小,选择潜力不大。位于外显子4的位点g.11532 A>G使原始密码子GGA突变为GGG,均编码甘氨酸,属同义突变。卡方适合性检验显示,除位点g.11468处于Hardy-Weinberg平衡状态(P>0.05)外,其余13个SNP位点均处于不平衡状态,即均受到选择压的影响。

表 1 SLA-DRB基因突变位点的遗传结构 Table 1 Genetic structure of mutation sites in SLA-DRB gene
位点Site 基因型Genotype
(频度Frequentness)
等位基因Allelic
(频率Frequency)
群体遗传结构指标Population genetic structure index
纯合度
(Ho)
位点杂
合度(He)
有效等位
基因(Ne)
多态信息
含量(PIC)
χ2 P
g.11369 CC(48) CG(9) GG(51) C(0.486) G(0.514) 0.500 0.500 2.000 0.374 8 74.98 5.00E-18
g.11387 AA(47) AG(10) GG(51) A(0.481) G(0.519) 0.500 0.500 2.000 0.374 6 71.66 3.00E-17
g.11440 CC(89) CT(11) TT(8) C(0.875) T(0.125) 0.781 0.219 1.280 0.194 8 30.84 3.00E-08
g.11468 CC(81) CT(27) TT(0) C(0.875) T(0.125) 0.781 0.219 1.280 0.194 8 2.20 0.138
g.11470 AA(39) AG(20) GG(49) A(0.454) G(0.546) 0.504 0.496 1.984 0.372 9 42.38 8.00E-11
g.11478 CC(78) CT(12) TT(18) C(0.778) T(0.222) 0.655 0.345 1.527 0.285 8 49.73 2.00E-12
g.11532 AA(33) AG(13) GG(62) A(0.366) G(0.634) 0.536 0.464 1.866 0.356 4 59.23 1.00E-14
g.11708 CC(12) CT(8) TT(88) C(0.148) T(0.852) 0.748 0.252 1.337 0.220 4 68.42 1.00E-16
g.11829 AA(51) AG(18) GG(39) A(0.556) G(0.444) 0.506 0.494 1.976 0.371 8 47.40 6.00E-12
g.11838 CC(29) CT(17) TT(62) C(0.347) T(0.653) 0.547 0.453 1.828 0.350 5 46.02 1.00E-11
g.11911 TT(88) TG(7) GG(13) T(0.847) G(0.153) 0.741 0.259 1.350 0.225 6 60.69 7.00E-15
g.11944 CC(78) CG(10) GG(20) C(0.769) G(0.231) 0.645 0.355 1.550 0.292 2 59.10 1.00E-14
g.11945 AA(48) AG(21) GG(39) A(0.542) G(0.458) 0.504 0.496 1.984 0.373 2 39.98 3.00E-10
g.11997 CC(50) CT(16) TT(42) C(0.537) T(0.463) 0.503 0.497 1.988 0.373 6 52.23 3.00E-13
2.3 SLA-DRB基因位点连锁分析及单体型分析

根据连锁不平衡分析的定义,2个位点间的r2≥0.33,即2个位点处于强连锁状态。对SLA-DRB基因的14个SNP位点进行连锁不平衡分析,发现所有位点均存在直接或间接的强连锁关系(图 2)。以14个SNP位点识别单体型,发现89头猪共有27种单体型,30种单体型组合。

图 2 单核苷酸多态性位点间的相关性网络(基于r2值的连锁不平衡分析) Fig. 2 Network between the loci of single nucleotide polymorphisms (linkage disequilibrium analysis based on r2 value)

通过提高强连锁的标准(r2),发现g.11468、g.11440、g.11478与g.11944各自与其他位点的关联性相对较低;g.11911与g.11708间的关联性高,而与其他位点的关联性低;g.11470、g.11945、g.11997、g.11387、g.11369与g.11829两两间的关联性高;g.11532与g.11838两者间的关联性高(图 2)。因此,将g.11470、g.11945、g.11997、g.11387、g.11369与g.11829归为紧密相关的SNP位点,定义为组合1;g.11532与g.11838归为紧密相关的SNP位点,定义为组合2;g.11911与g.11708归为紧密相关的SNP位点,定义为组合3。在此基础上,将所有样本的序列重新划分单体型(表 2):组合1共有11种单体型(H11~H111)和12种单体型组合(H11H11=35、H111H111=37、H11H111=4、H11H51=3、H11H41=3、H31H51=1、H91H111=1、H71H111=1、H61H101=1,H71H81=1、H11H61=1、H21H111=1);组合2共有4种单体型(H12~H42)和5种单体型组合(H42H42=48、H12H12=26、H12H42=9、H12H22=5、H12H32=1);组合3共有2种单体型(H13、H23)和3种单体型组合(H13H13= 73、H23H23=12、H13H23=4)。

表 2 单体型序列及频次信息 Table 2 Haplotype forecast and frequency
编号
ID
单体型
Haplotype
位点Locus 频次
Frequency
g.11369 g.11387 g.11470 g.11532 g.11708 g.11829 g.11838 g.11911 g.11945 g.11997
H11 CAAGGT C A A G G T 81
H21 CAAAGT C A A A G T 1
H31 CAGGGT C A G G G T 1
H41 CAGAAT C A G A A T 3
H51 CAGAAC C A G A A C 4
H61 CGGAAC C G G A A C 2
H71 GAAGGT G A A G G T 2
H81 GAGAAC G A G A A C 1
H91 GGAAGT G G A A G T 1
H101 GGAAGC G G A A G C 1
H111 GGGAAC G G G A A C 81
H12 AC A C 67
H22 AT A T 5
H32 GC G C 1
H42 GT G T 105
H13 TT T T 150
H23 CG C G 28
2.4 SLA-DRB基因单体型与生长发育性状的关联分析

性别、各位点频次大于3的单体型组合与生长发育性状的关系分析结果(图 3)表明:性别与腹围、体质量、胸宽、胸深、胸围和腰角宽的相关性显著,雌性较雄性大。g.11944和g.11440的位点多态性与所有体尺指标的相关性均不显著;g.11478的位点多态性与管围、尻高、体长和头长的相关性显著,单体型组合为C/T的个体具有较小的管围、尻高、体长和头长;g.11468的位点多态性与管围显著相关。组合3的位点多态性与除尻高、体质量和头长以外的所有体尺指标显著相关,这些体尺指标在单体型组合为TT/CG(H13H23)的个体中均偏小。组合2的位点多态性与管围、尻长和胸围的相关性显著,其中,单体型组合为AC/AT(H12H22)的个体有较小的尻长和胸围。组合1的位点多态性与管围、尻长、体高和胸深显著相关,其中,单体型组合为CAAGGT/CAGAAT(H11H41)的个体有较小的管围、尻长、体高和胸深。

图 3 单体型组合与各体尺指标的关联分析 Fig. 3 Association analysis between haplotype combination and growth traits 每一方格的颜色代表-log10转换的P值,颜色越深代表该值越大,即FDR值越小,该组合与该体尺指标的关联越显著;* * * FDR<0.001,* * FDR<0.01,* FDR<0.05;方格中标注的单体型组合的含义为在这一多态性位点上拥有该单体型组合的个体,其对应的体尺指标最小;下同 The color of tiles stands for -log10 transformed P value, the darker, the bigger, then FDR value is smaller; * * * FDR < 0.001, * * FDR < 0.01, * FDR < 0.05; individuals with haplotype combinations labelled on a certain tile tend to have smaller corresponding growth trait; the same below

管围受SLA-DRB基因单体型的影响最突出,共与5个位点存在关联(图 3),且g.11478、g.11468、组合1和组合2与管围存在极显著的关联(FDR<0.001)。g.11478的单体型组合为C/C、g.11468的单体型组合为C/C、组合1的单体型组合为GGGAAC/ GGGAAC(H111H111)和组合2的单体型组合为GT/GT(H42H42)的个体具有更大的管围(图 4)。

图 4 单体型组合与管围的关联 Fig. 4 Association between haplotype combination and cannon bone circumference
2.5 SLA-DRB基因单体型与管围率的关联分析

单体型组合与管围率间关系的分析结果(图 5)表明:g.11478、g.11468、组合1和组合2的位点多态性与基于额宽、腹围、尻高、体长、体高、头长、头深、胸宽、胸深、胸围和腰角宽的管围率均存在极显著的关联(FDR<0.001),提示这些位点多态性与藏猪的管围密切相关。与这些位点多态性关联最紧密的是基于体高的管围率,其中,g.11478的单体型组合为C/C、g.11468的单体型组合为C/C、组合1的单体型组合为GGGAAC/GGGAAC(H111H111)和组合2的单体型组合为GT/GT(H42H42)的个体有更大的管围率(图 6),与管围的结果一致,提示其可以有效反映个体的体型、健康状况。基于其他长度指标的管围率大致具有相同的趋势。与基于长度指标的管围率不同,基于体质量的管围率与g.11478、g.11468、组合1和组合2的位点多态性无显著相关性。

图 5 单体型组合与管围率的关联分析 Fig. 5 Association analysis between haplotype combination and the ratio of cannon bone circumference (CBC) to other growth traits 方格中标注了2种单体型组合,具有上排单体型组合的个体的管围率偏大,具有下排单体型组合的个体的管围率偏小;M.雄性,F.雌性 Two types of haplotype combinations were labelled in tiles, individuals with the upper haplotype combinations tend to have larger ratio of CBC to corresponding growth trait, while those with lower haplotype combinations tend to have smaller ratio of CBC to corresponding growth trait; M. male, F. female

图 6 单体型组合与基于体高的管围率的关联 Fig. 6 Association between haplotype combination and the ratio of cannon bone circumference to the body height
3 讨论

为了推进藏猪实验动物化的培育工作,本课题组从四川省阿坝藏族羌族自治州引入全同胞藏猪,在四川省蒲江县进行近亲繁育,早期的研究工作探讨了垂体特异性转录因子(pituitary specific transcription factor-1,POU1F1)第4内含子SNP位点多态性对该藏猪种群生长状况的影响,发现POU1F1基因的g.13952A、g.14192T、g.14384G、g.14464T位点分别为A/A、T/T、G/G和T/T的个体在所测个体中生长最缓慢,提示POU1F1基因对该藏猪群体的小型化筛选具有指导价值(何凡等,2013)。

本文进一步研究了该群体中SLA-DRB基因部分序列(部分内含子3、外显子4、内含子4、外显子5)的多态性,共发现了14个SNP位点。根据这些位点的连锁相关性,共识别出3个组合和4个单一位点。分析表明,不同位点组合/单一位点与藏猪个体生长发育性状间相关的程度不同。

3.1 SLA-DRB基因SNP位点多态性与个体生长发育的关系

已有研究表明,SLA-DRB基因SNP位点多态性与猪的繁殖性能(总产仔数、产活仔数、初生窝重、断奶仔猪数等)密切相关(吴圣龙等,2007张冬杰,2010)。本研究结果表明,SLA-DRB基因与藏猪个体的生长发育速率亦存在关联。其中,组合3与藏猪的发育速率最相关,影响个体的额宽、腹围、管围、尻长、体长、体高、头深、胸宽、胸深、胸围、腰角宽等体尺指标。值得注意的是,具有H13H23(TT/CG)单体型组合个体的上述所有指标均最小,意味着该基因型的个体生长较为缓慢,整体体型偏小。组合H13H13与H23H23个体间的体尺指标差异不大,对小型化藏猪个体的筛选意义有限。除组合3以外,组合1和g.11478分别与体高和体长相关联,具有单体型组合H11H41(CAAGGT/CAGAAT)的个体具有较小的体高,而g.11478的单体型组合为C/T的个体具有较小的体长。表明SLA-DRB基因SNP位点多态性与藏猪个体生长趋势存在显著关联,提示其可作为潜在的MAS遗传标记候选基因应用于特定生长趋势藏猪品系的筛选。

3.2 SLA-DRB基因SNP位点与管围的关系

本研究发现,与SLA-DRB基因SNP位点多态性最为相关的藏猪体尺指标为管围。管围与管骨的粗细密切相关,而管骨为长骨,其内具有丰富的骨髓。脊椎动物骨髓内含有丰富的间充质干细胞,属于多能干细胞,可以分化为成骨细胞、软骨细胞、心肌细胞、骨髓脂肪细胞和beta-胰岛细胞等(Ankrum et al., 2014)。管围和管围率较大的藏猪个体理应具有较多的管骨骨髓,因而可能在骨发育、心肌功能等方面更具优越性。除此以外,骨髓更是淋巴系统的重要组成部分,是造血组织细胞分裂分化为淋巴细胞的重要场所,与胸腺共同完成淋巴细胞的产生和早期选择(Yoffey,1975)。SLA基因作为猪的主要组织相容性复合体基因,在淋巴细胞的产生和选择过程中发挥着重要作用(Nakayama,2015)。研究结果表明,管围与SLA-DRB基因多态性之间存在极显著的关联,这与SLA-DRB基因在免疫系统中所扮演的角色相符合。因此,下一步的研究应着力于证实管围与管骨骨髓含量、淋巴血细胞造血能力以及个体免疫能力之间的关联性,为基于SLA-DRB基因的MAS遗传标记开发提供完善的理论依据。

除潜在的免疫力指示作用以外,管围及管围率直接影响藏猪的体型。研究表明,马的管围率(管围/体高)与其速度、步态特征具有显著关联,管围率大的个体奔跑速率较低(李敏等, 2014a, 2014b),藏猪可能存在类似的关联。另外,猪的管围是肉猪市场的重要选择指标,因为管围可能与生长速率、性成熟时间等因素相关(Myers et al., 2012)。因此,管围在培育实验动物猪、宠物猪、娱乐猪和食用猪等方面均具有重要参考价值。而SLA-DRB基因与管围之间的显著关联则为早期遗传辅助筛选提供了潜在途径。

根据本研究结果,g.11468、g.11478、组合1和组合2的单体型组合分别为C/C、C/C、GGGAAC/GGGAAC和GT/GT时,个体趋向于有较大的管围率,而这些单体型组合均为纯合性。这意味着一旦证实了管围与免疫能力之间的关联,这些单体型组合可以作为实验用藏猪的筛选标准。

综上所述,本研究测定了89头人工藏猪种群的部分SLA-DRB基因序列,并从中鉴定出了14个SNP位点,根据这些SNP位点的连锁紧密程度识别出3个组合和4个位点。经过与体尺指标的关联性分析,发现组合1的多态性与藏猪个体的体型存在关联,而g.11468、g.11478、组合1和组合2的位点多态性则与个体的管围、管围率显著相关。本研究结果为SLA-DRB基因作为MAS遗传标记候选基因提供了理论依据,为藏猪的定向遗传筛选提供了基础数据和参考资料。

参考文献
何凡, 帅素容, 江建平. 2013. 藏猪人工培育种群POU1F1第4内含子SNP检测及其与生长发育性状关联分析[J]. 四川动物, 32(5): 689–694.
何凡. 2013.人工培育藏猪种群SLA-DRBPOU1F1基因SNP与生长发育性状的关联分析[D].雅安: 四川农业大学. http://cdmd.cnki.com.cn/Article/CDMD-10626-1014102228.htm
黄萌, 蔡琳. 2009. Haploview与PHASE在单体型研究中的应用[J]. 福建医科大学学报, 43(4): 310–313. DOI:10.3969/j.issn.1672-4194.2009.04.008
李敏, 孟军, 刘运飞, 等. 2014a. 伊犁马速步赛成绩与体尺性状的相关性研究[J]. 中国畜牧兽医, 41(10): 209–212.
李敏, 孟军, 王建文, 等. 2014b. 伊犁马速度赛体尺与速度的相关性分析研究[J]. 新疆农业科学, 51(11): 2117–2122.
李群. 2008. 几个小型猪种[J]. 农家参谋, 9: 7.
彭勇波, 李奎, 樊斌, 等. 2005. 猪SLA-DQA一个新SNP的发现及其遗传效应的研究[J]. 中国农业科学, 38(12): 2526–2530. DOI:10.3321/j.issn:0578-1752.2005.12.025
任丽华. 2006.西藏小型猪实验动物化的初步研究[D].广州: 第一军医大学. http://cdmd.cnki.com.cn/Article/CDMD-90023-2006176902.htm
王铁岗. 2009. 动物遗传育种基础[M]. 北京: 化学工业出版社.
吴圣龙, 鞠慧萍, 孙鹏翔, 等. 2007. 苏太猪SLA-DQBSLA-DRB基因第2外显子多态性及其与繁殖性能的关联分析[J]. 农业生物技术学报, 15(4): 606–611. DOI:10.3969/j.issn.1674-7968.2007.04.012
袁进, 顾为望. 2010.小型猪作为人类疾病动物模型在生物医学研究中的应用[C].乌鲁木齐: 第九届中国实验动物科学年会论文集. http://www.wanfangdata.com.cn/details/detail.do?_type=conference&id=7659321
张冬杰. 2010. SLA-DQA基因的多态性与仔猪初生重和30日龄体质量的相关性分析[J]. 黑龙江畜牧兽医(1): 1–2. DOI:10.3969/j.issn.1004-7034.2010.01.001
Ankrum JA, Ong JF, Karp JM. 2014. Mesenchymal stem cells:immune evasive, not immune privileged[J]. Nature Biotechnology, 32(3): 252–260. DOI:10.1038/nbt.2816
Bronner S, Murbach V, Peter JD, et al. 2002. Ex vivo pharmacodynamics of amoxicillin-clavulanate against β-lactamase-producing Escherichia coli in a yucatan miniature pig model that mimics human pharmacokinetics[J]. Antimicrobial Agents & Chemotherapy, 46(12): 3782–3789.
Garthoff LH, Henderson GR, Sager AO, et al. 2002. The Autosow raised miniature swine as a model for assessing the effects of dietary soy trypsin inhibitor[J]. Food & Chemical Toxicology, 40(4): 487–500.
Gautschi C, Gaillard C. 1990. Influence of major histocompatibility complex on reproduction and production traits in swine[J]. Animal Genetics, 21(2): 161–170. DOI:10.1111/age.1990.21.issue-2
Gustafsson K, Germana S, Hirsch F, et al. 1990a. Structure of miniature swine class Ⅱ DRB genes:conservation of hypervariable amino acid residues between distantly related mammalian species[J]. Proceedings of the National Academy of Sciences, 87(24): 9798–9802. DOI:10.1073/pnas.87.24.9798
Gustafsson K, Leguern C, Hirsch F, et al. 1990b. Class Ⅱ genes of miniature swine. Ⅳ. Characterization and expression of two allelic class Ⅱ DQB cDNA clones[J]. Journal of Immunology, 145(6): 1946–1951.
Imai H, Konno K, Nakamura M, et al. 2006. A new model of focal cerebral ischemia in the miniature pig[J]. Journal of Neurosurgery, 104(2 Suppl): 123–132.
Mallard BA, Kennedy BW, Wilkie BN. 1991. The effect of swine leukocyte antigen haplotype on birth and weaning weights in miniature pigs and the role of statistical analysis in this estimation[J]. Journal of Animal Science, 69(2): 559–564. DOI:10.2527/1991.692559x
Myers N, Marshall TT, Walker WR, et al. 2012. Selecting pigs for youth swine shows[EB/OL]. (2012-04)[2017-12-11]. https://edis.ifas.ufl.edu/4h210#FOOTNOTE_2.
Nakayama M. 2015. Antigen presentation by MHC-dressed cells[J]. Frontiers in Immunology, 5(672): 1–8.
Nei M. 1978. Estimation of average heterozygosity and genetic distance from a small number of individuals[J]. Genetics, 89(3): 583–590.
Ryska M, Kieslichová E, Pantoflícek T, et al. 2004. A model of acute hepatic failure in a minipig from the surgical and anesthesia point of view[J]. Rozhl Chir, 83(9): 436–442.
Sebert SP, Lecannu GS, Hucteau S, et al. 2005. Obesity induced during sexual maturation is linked to LDL-triacylglycerols in yucatan miniature swine[J]. British Journal of Nutrition, 94(2): 282–289. DOI:10.1079/BJN20051478
Shi YY, He L. 2005. SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci[J]. Cell Research, 15(2): 97–98. DOI:10.1038/sj.cr.7290272
Shia YC, Bradshaw M, Rutherford MS, et al. 1995. Polymerase chain reaction based genotyping for characterization of SLA-DQB and SLA-DRB alleles in domestic pigs[J]. Animal Genetics, 26(2): 91–100.
Shia YC, Gautschi C, Ling MS, et al. 1991. RFLP analysis of SLA haplotypes in Swiss large white and American hampshire pigs using SLA class Ⅰ and class Ⅱ probes[J]. Animal Biotechnology, 2(1): 75–91. DOI:10.1080/10495399109525749
Turk JR, Henderson KK, Vanvickle GD, et al. 2005. Arterial endothelial function in a porcine model of early stage atherosclerotic vascular disease[J]. International Journal of Experimental Pathology, 86(5): 335–345. DOI:10.1111/j.0959-9673.2005.00446.x
Yoffey JM. 1975. Bone marrow lymphocytes[J]. Lymphology, 8(4): 154–157.