畜牧兽医学报  2021, Vol. 52 Issue (3): 596-609. DOI: 10.11843/j.issn.0366-6964.2021.03.004    PDF    
RXRB基因SNPs检测及其与生长育肥和繁殖性状的关联分析
彭雅鑫1, 刘军1, 赵诗瑜1, 徐在言2, 左波1     
1. 华中农业大学 农业部猪遗传育种重点实验室, 武汉 430070;
2. 华中农业大学动物医学院 基础兽医系, 武汉 430070
摘要:旨在研究猪RXRB基因多态性位点,并筛选与猪生长育肥和繁殖性状显著相关的SNPs,为种猪的遗传改良提供新的遗传标记位点。本研究通过采集962头健康大白猪(美系大白459头,法系大白503头)的血液DNA,利用直接测序法和PCR-RFLP技术分别检测两个品系大白猪RXRB基因SNP位点,并进行遗传多态性分析。利用SAS 8.0软件中的混合线性模型(Mixed)对RXRB基因SNP位点与表型性状值进行关联分析。利用Haploview 4.2软件对RXRB基因SNP位点进行连锁不平衡分析。结果显示,在两个品系大白猪群中共检测到10个SNPs,其中,美系大白有5个SNPs,法系大白有5个SNPs。经X2适合性检验,该10个SNPs在大白猪群体中均处于Hardy-Weinberg平衡状态(P>0.05)。关联分析结果表明,rs340542491、rs330162688和rs326226767与达100 kg体重日龄达到显著相关(P < 0.05);rs340542491、rs330162688、rs336609453和rs332169586与活体背膘厚达到显著相关(P < 0.05);rs340542491、rs325538588、rs691889709和rs323107853与初生重达到显著相关(P < 0.05);在体长中,仅rs324141460的GA基因型个体显著大于GG基因型(P < 0.05);rs325538588和rs80789331与眼肌面积达到显著相关(P < 0.05);在乳头数中,rs340542491和rs330162688与左乳头数达到极显著相关(P < 0.01),rs324141460与左、右乳头数达到极显著相关(P < 0.01),rs336609453和rs332169586与右乳头数达到显著相关(P < 0.05)。连锁不平衡分析结果显示,在美系大白中,rs326226767-rs325538588处于强连锁不平衡状态(D’=0.96,r2=0.91>0.33),rs330162688-rs324141460处于完全连锁状态(D’=1),rs324141460-rs340542491处于完全连锁状态(D’=1),单倍型组合关联分析结果与单个SNP关联分析结果一致,单倍型组合与达100 kg体重日龄、活体背膘厚、初生重和乳头数等性状存在极显著(P < 0.01)或显著(P < 0.05)相关;在法系大白中,rs691889709-rs323107853处于完美连锁状态(D’=1,r2=1),单倍型组合关联分析结果显示,单倍型组合与体长和乳头数达到极显著(P < 0.01)或显著(P < 0.05)相关。以上结果提示,RXRB基因内的10个SNPs可以作为猪育种改良的分子标记,同时对深入研究该基因功能具有一定的参考意义。
关键词RXRB    关联分析    连锁不平衡    单倍型    生长育肥    繁殖性状        
Detection of SNPs in Porcine RXRB Gene and Their Association Analysis with Growth, Fattening and Reproduction Traits
PENG Yaxin1, LIU Jun1, ZHAO Shiyu1, XU Zaiyan2, ZUO Bo1     
1. Key Laboratory of Pig Genetics and Breeding of Ministry of Agriculture, Huazhong Agricultural University, Wuhan 430070, China;
2. Department of Basic Veterinary Medicine, College of Veterinary Medicine, Huazhong Agricultural University, Wuhan 430070, China
Abstract: The study aimed to identify the polymorphism of RXRB gene, screen SNPs significantly associated with growth, fattening and reproductive traits, and provide novel genetic markers for the genetic improvement of pigs. In this study, the blood genomic DNA of 962 healthy Large White pigs (459 American Large White pigs and 503 French Large White pigs) were sampled. The SNPs of RXRB gene were detected by direct sequencing, and PCR-RFLP methods and used to perform the genetic polymorphism analysis. The mixed linear model (Mixed) in SAS 8.0 software was used to analyze the association between the SNPs of RXRB gene and phenotypic traits. The Haploview 4.2 software was used to analyze the linkage disequilibrium of the SNPs. The results showed that a total of 10 SNPs were detected in the two strains of Large White pigs, of which there were 5 SNPs in the American Large White pigs and 5 SNPs in the French Large White pigs. According to the result of χ2 test, the 10 SNPs were all under the Hardy-Weinberg equilibrium status (P>0.05). The results of association analysis showed that rs340542491, rs330162688 and rs326226767 were significantly associated with days at 100 kg body weight (P < 0.05). The rs340542491, rs330162688, rs336609453 and rs332169586 were significantly associated with backfat thickness (P < 0.05). The rs340542491, rs325538588, rs691889709 and rs323107853 were significantly associated with birth weight (P < 0.05). For body length, only the individuals with genotype GA at rs324141460 were significantly larger than those with GG genotype (P < 0.05). The rs325538588 and rs80789331 were significantly associated with the loin muscle area (P < 0.05). The rs340542491 and rs330162688 were extremely significantly associated with number of left nipples (P < 0.01). The rs324141460 was extremely significantly associated with number of left and right nipples (P < 0.01). The rs336609453 and rs332169586 were significantly associated with the number of right nipples (P < 0.05). The results of linkage disequilibrium analysis showed that rs326226767-rs325538588 was in a strong linkage disequilibrium status (D'=0.96, r2=0.91>0.33), rs330162688-rs324141460 was in a fully linkage status (D'=1), and rs324141460-rs340542491 was in a fully linked status (D'=1) in American Large White pigs. The results of association analysis of haplotype combinations were consistent with the results of the association analysis of single SNP; they were extremely significant (P < 0.01) or significant (P < 0.05) associated with days at 100 kg body weight, backfat thickness, birth weight and the number of nipples. In French Large White pigs, rs691889709-rs323107853 was in a fully linkage status (D'=1, r2=1). The associations between haplotype combinations and the body length, number of nipples reached extremely significant (P < 0.01) or significant (P < 0.05) levels. The results suggest that the 10 SNPs of RXRB gene can be used as molecular markers for pig breeding improvement, and provide the basis for the in-depth study of the function of RXRB.
Key words: RXRB    association analysis    linkage disequilibrium    haplotype    growth and fattening    reproduction traits    pig    

随着高通量测序技术的飞速发展和测序成本的降低,国内外许多研究人员利用猪的高密度芯片对猪重要经济性状进行了GWAS研究,鉴定出大量的SNPs标记、QTLs和候选基因,为猪的遗传育种提供了必要的分子基础[1]。其中,SNP分子标记多态性高、密度大、遗传稳定性高、易于实现自动化高通量基因分型,在表型与基因型关联性研究中起重要作用[2],因而被广泛应用于人类和动植物各项遗传学研究,包括遗传图谱构建、连锁分析、关联分析、基因组预测等[3]。分子遗传标记与动物分子育种技术加快了品种改良和新品系培育速度[4]。目前,SNP的应用已成为了研究猪遗传育种领域不可或缺的工具之一。

视黄酸X受体β(retinoid X receptor beta, RXRB)基因定位在猪7号染色体上,由10个外显子和9个内含子组成。视黄素X受体(retinoid X receptor, RXR)是一种能使视黄酸受体(retinoic acid receptor, RAR)发挥正常作用的辅助蛋白,它与RAR可形成异源二聚体,使RAR与其配体更容易结合,从而增强RAR的转录和与配体的亲和力[5-7],也可以通过与其他核激素受体形成二聚体来激活多种功能,从而促进不同信号通路的活动[8]。RXR作为一种转录因子,属于核转录调节因子超家族一员,该家族还包含了维生素D受体、甲状腺激素受体、类固醇激素受体以及所谓的孤儿受体。自1990年RXR被分离并鉴定为核受体超家族成员以来,人们对核受体的生理调节产生了新的认识,从生物学上讲,RXR通过与过氧化物酶体增殖物激活受体(PPAR)和肝X受体(LXR)等其他核受体结合形成异二聚体,发挥关键作用[9-10]。研究发现,RXR既可以与其他核受体形成异源二聚体启动转录,也可以以同源二聚体形式启动转录[11]。通过不同的结合体在细胞质和核中定位不同,从而在调控细胞增殖、分化和凋亡中发挥重要作用[12]。目前,已知RAR受体家族是由全反式视黄酸(all-trans retinoic acid, at-RA)和9-顺式视黄酸(9-cis-RA)激活的,而RXR受体家族仅由9-cis-RA激活,表明9-cis-RA对RAR/RXR的形成是必需的[13]

RXRB作为视黄酸的一个编码基因,可以介导视黄酸的生理活性。Kastner等[14]研究发现,在RXRB基因突变和缺乏RXRB的小鼠中约有50%在出生前或出生时死亡,而幸存的雄性小鼠因其精子发生异常而不育。另外,RXRB无效突变导致Sertoli细胞脂质代谢异常和精子发生异常[15]。Mascrez等[16]进一步研究发现,Sertoli细胞中的胆固醇稳态取决于TIF2共激活因子和RXRB/LXRB异二聚体。另有研究报道,在Sertoli细胞中RXRB基因的体细胞靶向性失活可导致精子释放失败,胆固醇蓄积并随后使睾丸变性[17]。在表达模式上,有研究发现,RXRB基因mRNA在不同日龄大白猪的心、肾、脾、大肠、子宫、卵巢中持续表达,其中在大肠、子宫和卵巢中持续高表达[18]。黄锐等[19]以长白山野猪和东北民猪为材料,检测了RXRB基因在不同骨骼肌组织中mRNA表达量,结果发现,该基因在两个品种猪肌肉组织中的表达量普遍高于非肌肉组织,表明RXRB可能是影响猪肉品质的候选基因之一。Pena等[20]利用高密度SNP阵列对382头伊比利亚猪进行GWAS分析,筛选到了与肌内脂肪中饱和脂肪酸和不饱和脂肪酸含量密切相关的候选基因RXRB。以上这些研究都表明了RXRB基因可能与动物的生长繁殖性能相关。本研究将RXRB基因,作为候选基因研究其多态性与猪重要经济性状的相关性,以期筛选到与猪生长、繁殖性状显著相关的SNPs,为今后种猪遗传改良提供新的遗传标记位点。

1 材料与方法 1.1 试验猪群及性状测定

本研究样本采集及数据收集于浙江省某企业。选取来自不同血缘世代且生长健康的后备大白猪作为研究对象,将猪群按品系分为美系大白猪和法系大白猪,其中美系大白猪共459头(公猪85头,母猪374头),法系大白猪共503头(公猪18头,母猪485头),上述试验猪群的饲养环境和管理方式一致。通过对试验猪群进行前腔静脉采血,并利用血液基因组DNA中量提取试剂盒(离心柱型)提取基因组DNA,将样品DNA保存于-80 ℃超低温冰箱中备用。测定的表型性状主要包括:达100 kg体重日龄(d)、活体背膘厚(mm)、初生重(kg)、体长(cm)、眼肌面积(cm2)、左乳头数和右乳头数。其中,当猪的体重达到(100±5) kg时,利用B型超声设备(ESAOTE,Mylob Touch VET)测定猪活体背膘厚和眼肌面积,活体背膘厚和眼肌面积测定位置在猪倒数第3、4根肋骨距离背中线4~5 cm处。所有测量数据由育种软件GBS 5.0校正后得到。

1.2 引物设计及PCR扩增

根据NCBI中GenBank上提供的猪RXRB基因序列(登录号为NC_010449.5),应用Primer Premier 5.0软件进行引物设计,并由上海生工生物工程技术有限公司合成。所用引物及目的片段长度如表 1所示。PCR扩增反应总体积为20.0 μL:DNA模板1.0 μL,2×Taq PCR Mix 10.0 μL,ddH2O 8.0 μL,正、反向引物(10 μmol·L-1)各0.5 μL。PCR反应程序:95 ℃预变性5 min;95 ℃变性40 s,退火30 s(退火温度见表 1),72 ℃延伸,共36个循环;最后72 ℃延伸5 min。PCR产物于1.5% 琼脂糖凝胶电泳检测,以5 μL的DNA Marker DL2000作为参照。电泳结束后利用凝胶成像系统观察扩增结果。

表 1 RXRB引物序列 Table 1 RXRB primer sequences
1.3 PCR产物酶切反应体系及条件

在本研究中检测到的10个SNPs中,仅多态性位点rs326226767存在合适的酶切位点,可利用PCR-BglⅠ-RFLP方法对该位点进行基因分型。其余9个SNPs则利用直接测序法进行基因分型。测序结果用Chromas软件查看。酶切反应体系及条件:PCR产物酶切反应体积为20.0 μL,其中10× FastDigest Green Buffer 2.0 μL,限制性内切酶0.5 μL(5 U),ddH2O 7.5 μL,PCR产物10.0 μL,将样品混匀后离心,37 ℃酶切1 h。酶切产物于1.5%琼脂糖凝胶电泳检测,以5 μL的DNA Marker DL2000作为参照。电泳结束后利用凝胶成像系统观察基因型。

1.4 数据处理与统计分析

1.4.1 RXRB基因SNPs位点的多态性信息含量(PIC)计算   根据试验结果计算各品系群体中不同等位基因的基因型频率与基因频率分布情况,多态信息含量用于对标记基因多态性的估计。PIC>0.5,则该位点为高度多态性,PIC < 0.25,则该位点为低度多态性,0.25 < PIC < 0.5,则该位点为中度多态性。公式:$ PIC = 1 - \left( {\sum\nolimits_{i = 1}^n {P_i^2} } \right) - \left( {\sum\nolimits_{i = 1}^{n - 1} {\sum\nolimits_{j = i + 1}^n {2P_i^2P_j^2} } } \right) $。其中,PiPj为第i个和第j个等位基因频率,n为复等位基因数。

1.4.2 RXRB基因各SNPs位点的连锁不平衡分析   本研究利用Haploview 4.2软件分析RXRB基因SNPs位点之间的连锁不平衡水平并计算出各单倍型频率。

1.4.3 分子标记的基因效应分析   本研究利用SAS 8.0软件中的混合线性模型(Mixed)进行多态性位点不同基因型和单倍型组合与猪生长繁殖性状表型值之间的关联分析,结果以“平均值±标准误”表示,P < 0.05判定为差异显著,P < 0.01判定为差异极显著。应用Excel软件分别计算基因型频率和等位基因频率,并进行χ2检验,分析其是否偏离Hardy-Weinberg平衡状态。分析模型:Tijkl= μ + Gi + Fj+ Sk + Bl+ eijklm。式中,Tijkl为性状表型值;μ为平均值;Gi为基因型效应;FjSkBl为固定效应,分别为家系、性别、批次效应,eijklm为残差效应。

2 结果 2.1 PCR扩增与酶切产物的电泳检测

将引物P1、P2和P3分别在美系大白猪和法系大白猪DNA混合池中进行PCR扩增,并将PCR产物送至上海生工生物工程技术有限公司进行测序。结果发现,在美系大白猪中,引物P1扩增产物检测到4个SNPs,分别为第一外显子区域的rs340542491、rs324141460和第一内含子区域的rs330162688、rs325538588;引物P2扩增产物中检测到1个SNP(rs326226767),位于第七外显子区域,且该位点由碱基G突变成碱基A,进而导致限制性内切酶BglⅠ酶切位点的改变,可采用PCR-BglⅠ-RFLP方法来进行基因型分型,酶切电泳结果见图 1。在法系大白猪中,引物P3扩增产物检测到5个SNPs,分别为第八内含子区域的rs691889709、rs336609453,第九外显子区域的rs80789331,第九内含子区域的rs332169586,第十外显子区域的rs323107853。上述10个SNPs的测序峰图如图 2所示。

M. DNA相对分子质量标准;1、2、5. GG基因型;3. AA基因型;4、6、7. GA基因型 M. DL2000 marker; 1, 2, 5. GG genotype; 3. AA genotype; 4, 6, 7. GA genotype 图 1 RXRB基因引物P2扩增片段的PCR-Bgl Ⅰ-RFLP酶切分型结果 Fig. 1 The PCR-Bgl Ⅰ-RFLP results of amplified fragments by primer P2 in RXRB gene
图 2 RXRB基因SNPs测序结果 Fig. 2 SNPs sequencing results in RXRB gene
2.2 大白猪RXRB基因10个SNPs的遗传多态性分析

利用直接测序法分别检测了459头美系大白猪中引物P1扩增区域内的4个SNPs和503头法系大白猪中引物P3扩增区域内的5个SNPs,同时利用PCR-Bgl Ⅰ-RFLP技术对418头美系大白猪中引物P2扩增区域内的1个SNP(rs326226767)进行了基因分型。上述10个SNPs均检测到3种基因型,基因型与等位基因频率分布见表 2X2适合性检验结果表明,上述10个SNPs在大白猪群体中均处于Hardy-Weinberg平衡状态(P>0.05)。

表 2 RXRB基因10个SNPs位点的基因型及等位基因频率 Table 2 The genotype and allele frequencies of 10 SNPs of RXRB gene

从群体遗传学角度分析大白猪RXRB基因SNPs的基因型和等位基因,由表 2可知,除位点rs80789331的多态性含量(PIC)为0.133,属于低度多态(PIC < 0.25),其余9个SNPs均属于中度多态(0.25 < PIC < 0.50),遗传变异处于中等水平,因此选择潜力较大。

2.3 大白猪RXRB基因多态性与生长育肥和繁殖性状的关联分析

从大白猪RXRB基因多态性与经济性状的关联分析(表 3)可以看出,在达100 kg体重日龄中,位点rs340542491、rs330162688和rs326226767均存在显著的关联性(P < 0.05);在活体背膘厚中,位点rs340542491和rs330162688的GG基因型个体的背膘厚度均显著大于GA和AA基因型(P < 0.05),GA和AA基因型个体间无显著差异(P>0.05),表明了等位基因G会降低瘦肉率,为不利等位基因,同时,位点rs336609453和rs332169586也存在显著的关联性(P < 0.05);在初生重中,位点rs340542491和rs325538588存在显著的关联性(P < 0.05),且GA基因型个体相比于其他基因型的初生重更大,在实际生产过程中初生重越大的仔猪通常具有更高的成活率,即GA基因型为有利基因型,同时,位点rs691889709和rs323107853也存在显著的关联性(P < 0.05);在体长中,位点rs324141460的GA基因型个体显著大于GG基因型(P < 0.05);在眼肌面积中,位点rs325538588的AA基因型个体显著大于GG和GA基因型(P < 0.05),位点rs80789331的AA基因型个体显著大于CC和CA基因型(P < 0.05);在乳头数中,位点rs340542491和rs330162688的AA基因型个体的左乳头数均极显著高于GA基因型(P < 0.01),位点rs324141460的AA基因型个体的左、右乳头数均极显著高于GG和GA基因型(P < 0.01),位点rs336609453的AA基因型个体的右乳头数显著高于AG基因型(P < 0.05),位点rs332169586的GG基因型个体的右乳头数显著高于GA基因型(P < 0.05),通常情况下,母猪的乳头数越多,其繁殖性能越强。另外,通过对两个品系猪表型值数据的统计分析发现,法系大白猪的平均体长和平均乳头数均要高于美系大白猪,这可能是不同的选育方向造成了两个品系猪群的差异。

表 3 RXRB基因10个SNPs与生长育肥和繁殖性状间的关联分析 Table 3 Association analysis between 10 SNPs of RXRB gene and growth, fattening and reproduction traits
2.4 RXRB基因10个SNPs连锁不平衡及单倍型分析

利用Haploview 4.2软件分别对美系大白猪和法系大白猪RXRB基因内处于Hardy-Weinberg平衡状态的10个SNPs进行连锁不平衡及单倍型分析。如图 3所示,在美系大白猪中的5个SNPs可组成一个单倍型模块,其中, rs326226767-rs325538588处于强连锁不平衡状态(D’=0.96,r2=0.91>0.33),rs330162688-rs324141460处于完全连锁状态(D’=1),rs324141460-rs340542491处于完全连锁状态(D’=1);在法系大白猪中的5个SNPs也可组成一个单倍型模块,其中, rs691889709-rs323107853处于完美连锁状态(D’=1,r2=1),这也解释了该2个SNPs基因型频率和关联分析结果一致性问题。如表 4表 5所示,美系大白猪和法系大白猪中的SNPs各能形成4种单倍型,分别将4种单倍型进行组合,剔除频率小于0.03的组合后,在美系大白猪中可产生8种双倍型,法系大白猪中可产生6种双倍型,并分别进行下一步的关联分析。

根据D’值和r2判断位点之间是否存在连锁,其中D’=1称之为完全连锁,D’=0为无连锁或连锁平衡,r2=1称之为完全连锁,r2>0.33提示为强连锁,r2=0为无连锁或连锁平衡 The linkage among SNPs was evaluated by the value of D' and r2, where D' = 1 indicated full linkage, D' = 0 indicated no linkage or linkage equilibrium, r2 = 1 indicated full linkage, r2>0.33 indicated strong linkage, r2=0 indicated no linkage or linkage equilibrium 图 3 RXRB基因10个SNPs位点连锁不平衡分析 Fig. 3 The linkage disequilibrium analysis among 10 SNPs of RXRB gene
表 4 RXRB基因单倍型及其频率 Table 4 The haplotypes and their frequencies of RXRB gene
表 5 RXRB基因单倍型组合及其频率 Table 5 The haplotype combinations and their frequencies of RXRB gene

表 6可知,美系大白猪RXRB基因产生的8种双倍型在达100 kg体重日龄、活体背膘厚、初生重和左乳头数等性状上存在极显著(P < 0.01)或显著差异(P < 0.05),对于达100 kg体重日龄、活体背膘厚和初生重性状,在生产实践中H1-H3为有利双倍型;对于左乳头数性状,H2-H2为有利双倍型。而法系大白猪RXRB基因产生的6种双倍型在体长和右乳头数性状上存在极显著(P < 0.01)或显著差异(P < 0.05),其中D1-D1个体能极显著提高右乳头数(P < 0.01),但极显著降低了体长(P < 0.01),因此需要根据实际情况进行选留。

表 6 RXRB基因单倍型组合与生长育肥和繁殖性状的关联分析 Table 6 Association analysis of RXRB haplotype combinations with growth, fattening and reproduction traits
3 讨论 3.1 RXRB基因群体遗传特征分析

在自然选择下,人工选择干预、遗传漂变等因素都会导致动物的一些核苷酸发生改变,从而使生物的遗传性状变得更加丰富多样[21]。群体遗传变异度的大小代表着群体遗传丰富度的高低,一般来说,群体遗传变异度越高,则群体遗传丰富度也越高,多态信息含量(PIC)常用来指示这一指标[22]。本研究中,除rs80789331外,其余9个SNPs均属于中度多态(0.25 < PIC < 0.50),说明遗传丰富度较高,具有一定的遗传变异潜力,能提供的选择效果更加丰富[23]。在两个品系大白猪群体中检测到的SNPs均处于Hardy-Weinberg平衡状态,说明在这两个群体中这些SNPs的选择压力不强,可能是未受到突变、选择、遗传漂变等影响[24],也可能是由于等位基因和基因型在适应性方面具有遗传优势,经过长期人工选择及大规模繁育后又重新达到平衡[25]

3.2 RXRB基因多态性与生长繁殖性状的关联分析

相关研究表明,RXRB基因作为核转录调节因子超家族一员,在生物体的胚胎发育、器官形成、细胞分化、细胞凋亡和维持生殖功能等方面具有重要作用[26]。李华振等[27]总结前人研究结果得出,RXRs参与组成的二聚体对动物繁殖发挥着多种重要的调控作用,但目前RXRs组成的二聚体在发情周期和季节性发情中具体的分子机制不够明晰。在支持细胞中,RXRB基因的表达能够维持胆固醇的稳态与精子形成[28]。成年哺乳动物体内都存在白色脂肪(WAT)和棕色脂肪(BAT),WAT的主要作用是将体内多余能量以脂肪的形式储存起来,过多的WAT会形成肥胖,而BAT能够将脂肪转化为热量,加强新陈代谢,促进白色脂肪消耗[29]。有研究表明,在大鼠脂肪组织中,WAT和BAT表达RARRXR亚型有着不同的方式,并且在不同的哺乳动物中都有相似的表达差别[30-31]。本研究发现,在美系大白猪群体中检测到的RXRB多态性位点rs340542491和rs330162688均与活体背膘厚性状显著相关,且GA和AA基因型均为有利基因型;在法系大白猪群体检测到的RXRB多态性位点rs336609453和rs332169586均与活体背膘厚性状显著相关,且杂合子基因型均为不利基因型,以上4个SNPs位点均可以作为降低大白猪活体背膘厚,提高瘦肉率的潜在分子标记。除此之外,我们发现在美系大白猪中RXRB基因rs325538588位点的AA基因型个体的眼肌面积显著大于GG和GA基因型个体,但AA基因型个体的初生重显著小于GA基因型个体,虽然AA基因型能显著增加眼肌面积,但是会降低初生重,可能会影响仔猪成活率,因此在生产实践中,需要根据实际情况进行选留。

有研究发现,PPAR-γ/RXR异二聚体在鼠胎盘发育中起着关键的调节作用[32]RXR基因在哺乳动物下丘脑-垂体-卵巢轴中表达,能作为连接动物繁殖发育和生理发情过程的关键转录因子[33]。另有研究发现,RXR基因表达与绵羊的繁殖有关,可能参与绵羊季节性发情上游基因调控[34]。另有研究报道,产仔数与总乳头数呈正相关[35]。通过与繁殖性状的关联分析发现,在美系大白猪中RXRB基因多态性位点rs340542491、rs324141460和rs330162688与左乳头数性状均达到极显著相关,AA基因型均为有利基因型;在法系大白猪中RXRB基因多态性位点rs336609453和rs332169586与右乳头数性状均达到显著相关,以上5个SNPs可作为提高母猪乳头数性状的候选位点,为母猪繁殖力提供新的分子标记和理论依据。

在以往的研究中,一般认为基因突变发生在内含子不会影响生物性状,发生在其他区域则可能影响生物性状,但随着研究者们的深入研究发现,真核生物中内含子参与基因剪切,内含子突变可能导致剪切准确性或效率发生变化,从而导致氨基酸编码改变,最终改变真核生物表型,如IGF2基因第三内含子3072位点的G/A突变可影响其与ZBED6转录因子的结合,进而改善巴马猪的肉产量[36]。即使内含子突变没有造成对应氨基酸的改变,但其可能会改变调控剪接的外显子基序,导致剪接位点失活;另外,有些内含子含有增强子序列,其可促进基因的专一性转录,从而对动物的表型产生影响[37]。本研究在两个品系猪群中共鉴定出5个位于RXRB内含子区域的SNPs(rs330162688、rs325538588、rs691889709、rs336609453和rs332169586),并且这5个SNPs分别与大白猪性状表型值有显著或极显著的关联。这5个SNPs是因果突变位点,还是仅与因果突变位点紧密连锁的分子标记,还需进行后续功能验证来阐明其遗传机理。

3.3 RXRB基因单倍型分析

动物基因对表型的影响可能受到多个突变位点的共同作用[38]。研究发现,由于强关联SNPs直接可以互相随机替代,单个SNP分析作用较小,而单倍型可以提供比单个SNP更为丰富的信息量,单倍型分析对于单染色体上杂合SNP的定位、疾病基因的挖掘和疾病治疗新方法的寻找有重要作用[39]。构建单倍型前需要进行LD分析,连锁不平衡常用D’和r2度量,当D’=0,r2=0时说明处于完全连锁平衡状态,当D’=1,r2=1时说明处于完全连锁不平衡状态。先前有研究表明,等位基因效应在不同品系猪中存在差异,推测是由于不同品系猪的连锁区间和随机突变不同所引起的[40-41]。在本研究的美系大白猪和法系大白猪中,分别对符合Hardy-Weinberg平衡的10个SNPs进行了连锁不平衡分析发现,美系大白猪中单倍型组合关联分析结果与单个SNP关联分析结果一致,均与达100 kg体重日龄、活体背膘厚、初生重和乳头数性状达到显著相关;但法系大白猪中单倍型组合关联分析结果显示出不同双倍型个体间的体长具有极显著差异,这与单个SNP关联分析结果不一致,考虑到该研究的样本量,有些单倍型所占比例太小,产生的试验结果具有一定的误差,后期还需进一步扩大样本量来验证。

4 结论

本研究检测到RXRB基因内的10个SNPs,对大白猪达100 kg体重日龄、活体背膘厚、初生重和乳头数等性状均具有显著的遗传效应,表明RXRB基因可以作为影响猪生长繁殖性状的候选基因,上述10个SNPs可以作为猪育种改良的分子标记,同时为深入研究该基因功能具有一定的参考意义。

参考文献
[1] 赵谦, 浦亚斌, 关伟军, 等. 猪重要性状全基因组关联分析的研究进展[J]. 畜牧兽医学报, 2015, 46(6): 873–881.
ZHAO Q, PU Y B, GUAN W J, et al. Research progress of genome-wide association studies for important traits in pig[J]. Acta Veterinaria et Zootechnica Sinica, 2015, 46(6): 873–881. (in Chinese)
[2] 王晨阳, 王璐, 张锐虎, 等. SNP标记在动物遗传育种及人类疾病动物模型研究中的应用[J]. 中国比较医学杂志, 2019, 29(4): 120–125.
WANG C Y, WANG L, ZHANG R H, et al. Application of single nucleotide polymorphism markers in animal genetic breeding and animal models of disease[J]. Chinese Journal of Comparative Medicine, 2019, 29(4): 120–125. (in Chinese)
[3] 王琦, 朱迪, 王宇哲, 等. 全基因组SNP分型策略及基因组预测方法的研究进展[J]. 畜牧兽医学报, 2020, 51(2): 205–216.
WANG Q, ZHU D, WANG Y Z, et al. Research progress of genomic-wide SNP genotyping and genomic prediction methods[J]. Acta Veterinaria et Zootechnica Sinica, 2020, 51(2): 205–216. (in Chinese)
[4] 朱向博, 张丽. 单核苷酸多态性及其在畜牧兽医领域的研究进展[J]. 现代畜牧兽医, 2020(7): 48–51.
ZHU X B, ZHANG L. Review of single nucleotide polymorphism and research progress in animal husbandry[J]. Modern Journal of Animal Husbandry and Veterinary Medicine, 2020(7): 48–51. (in Chinese)
[5] EGEA P F, ROCHEL N, BIRCK C, et al. Effects of ligand binding on the association properties and conformation in solution of retinoic acid receptors RXR and RAR[J]. J Mol Biol, 2001, 307(2): 557–576. DOI: 10.1006/jmbi.2000.4409
[6] VIERA-VERA J, GARCÍA-ARRARÁS J E. Molecular characterization and gene expression patterns of retinoid receptors, in normal and regenerating tissues of the sea cucumber, Holothuria glaberrima[J]. Gene, 2018, 654: 23–35. DOI: 10.1016/j.gene.2018.01.102
[7] LE MAIRE A, TEYSSIER C, BALAGUER P, et al. Regulation of RXR-RAR heterodimers by RXR- and RAR-specific ligands and their combinations[J]. Cells, 2019, 8(11): 1392. DOI: 10.3390/cells8111392
[8] DE ALMEIDA N R, CONDA-SHERIDAN M. A review of the molecular design and biological activities of RXR agonists[J]. Med Res Rev, 2019, 39(4): 1372–1397. DOI: 10.1002/med.21578
[9] WAGNER C E, JURUTKA P W, MARSHALL P A, et al. Retinoid X receptor selective agonists and their synthetic methods[J]. Curr Top Med Chem, 2017, 17(6): 742–767. DOI: 10.2174/1568026616666160617091559
[10] HIEBL V, LADURNER A, LATKOLIK S, et al. Natural products as modulators of the nuclear receptors and metabolic sensors LXR, FXR and RXR[J]. Biotechnol Adv, 2018, 36(6): 1657–1698. DOI: 10.1016/j.biotechadv.2018.03.003
[11] LE MAIRE A, BOURGUET W. Retinoic acid receptors: structural basis for coregulator interaction and exchange[J]. Subcell Biochem, 2014, 70(1): 37–54.
[12] BRAZDA P, SZEKERES T, BRAVICS B, et al. Live-cell fluorescence correlation spectroscopy dissects the role of coregulator exchange and chromatin binding in retinoic acid receptor mobility[J]. J Cell Sci, 2011, 124(21): 3631–3642. DOI: 10.1242/jcs.086082
[13] DILWORTH F J, FROMENTAL-RAMAIN C, YAMAMOTO K, et al. ATP-driven chromatin remodeling activity and histone acetyltransferases act sequentially during transactivation by RAR/RXR in vitro[J]. Mol Cell, 2000, 6(5): 1049–1058. DOI: 10.1016/S1097-2765(00)00103-9
[14] KASTNER P, MARK M, LEID M, et al. Abnormal spermatogenesis in RXR beta mutant mice[J]. Genes Dev, 1996, 10(1): 80–92. DOI: 10.1101/gad.10.1.80
[15] GANSMULLER A, MARK M, KASTNER P, et al. RXRβ null mutation leads to defects in sertoli cell lipid metabolism, and abnormal spermatogenesis[J]. Biol Cell, 1996, 86(2-3): 191.
[16] MASCREZ B, GHYSELINCK N B, WATANABE M, et al. Ligand‐dependent contribution of RXRβ to cholesterol homeostasis in sertoli cells[J]. EMBO Rep, 2004, 5(3): 285–290. DOI: 10.1038/sj.embor.7400094
[17] VERNET N, DENNEFELD C, KLOPFENSTEIN M, et al. Retinoid X receptor beta (RXRB) expression in sertoli cells controls cholesterol homeostasis and spermiation[J]. Reproduction, 2008, 136(5): 619–626. DOI: 10.1530/REP-08-0235
[18] 杨国伟, 刘娣, 张冬杰, 等. 不同日龄大白猪视黄酸β基因的表达研究[J]. 黑龙江畜牧兽医, 2009(1): 14–16.
YANG G W, LIU D, ZHANG D J, et al. Research on expression of landrace RARβ gene on different days[J]. Heilongjiang Animal Science and Veterinary Medicine, 2009(1): 14–16. (in Chinese)
[19] 黄锐, 秦宁, 张英英, 等. 长白山野猪和东北民猪RXRB基因mRNA组织表达特性分析[J]. 黑龙江畜牧兽医, 2014(11): 5–8.
HUANG R, QIN N, ZHANG Y Y, et al. Characteristic analysis of the mRNA tissue expression of the RXRB gene in Changbai Mountain wild boar and Northeast Min pig[J]. Heilongjiang Animal Science and Veterinary Medicine, 2014(11): 5–8. (in Chinese)
[20] PENA R N, NOGUERA J L, GARCÍA-SANTANA M J, et al. Five genomic regions have a major impact on fat composition in Iberian pigs[J]. Sci Rep, 2019, 9(1): 2031.
[21] MÄKI-TANILA A, WEBSTER L. Heritability, SNP, inbreeding, dairy cattle, genomic selection—and other keywords[J]. J Anim Breed Genet, 2019, 136(1): 1–2.
[22] AZIMU W, MANATBAY B, LI Y, et al. Genetic diversity and population structure analysis of eight local chicken breeds of Southern Xinjiang[J]. Br Poult Sci, 2018, 59(6): 629–635.
[23] SANTOS C P, AGUIAR A F, GIOMETTI I C, et al. High final energy of gallium arsenide laser increases MyoD gene expression during the intermediate phase of muscle regeneration after cryoinjury in rats[J]. Lasers Med Sci, 2018, 33(4): 843–850. DOI: 10.1007/s10103-018-2439-3
[24] 覃媛钰, 张依裕, 罗华伦, 等. SFRP5基因遗传变异对长顺绿壳蛋鸡胸肌肌肉品质的影响[J]. 农业生物技术学报, 2019, 27(4): 703–711.
QIN Y Y, ZHANG Y Y, LUO H L, et al. Effect of genetic variation of SFRP5 gene on breast muscle quality of Changshun blue-eggshell chicken (Gallus gallus domesticus)[J]. Journal of Agricultural Biotechnology, 2019, 27(4): 703–711. (in Chinese)
[25] 黄永震, 王居强, 马桂变, 等. 郏县红牛SREBP-1c基因84 bp缺失突变的检测及其对6个生长性状的影响[J]. 畜牧兽医学报, 2010, 41(10): 1225–1231.
HUANG Y Z, WANG J Q, MA G B, et al. Deletion mutation of the bovine SREBP-1c gene of Jiaxian red cattle and its affecting on growth traits[J]. Acta Veterinaria et Zootechnica Sinica, 2010, 41(10): 1225–1231. (in Chinese)
[26] NAGATA T, KANNO Y, OZATO K, et al. The mouse Rxrb gene encoding RXRβ: genomic organization and two mRNA isoforms generated by alternative splicing of transcripts initiated from CpG island promoters[J]. Gene, 1994, 142(2): 183–189.
[27] 李华振, 狄冉, 刘秋月, 等. 类视黄醇X受体基因调控动物季节性发情和发情周期的研究进展[J]. 畜牧兽医学报, 2019, 50(8): 1525–1535.
LI H Z, DI R, LIU Q Y, et al. Research progress on the regulative role of retinoid X receptors gene in animal seasonal estrus and estrous cycle[J]. Acta Veterinaria et Zootechnica Sinica, 2019, 50(8): 1525–1535. (in Chinese)
[28] GELY-PERNOT A, RAVERDEAU M, TELETIN M, et al. Retinoic acid receptors control spermatogonia cell-fate and induce expression of the SALL4A transcription factor[J]. PLoS Genet, 2015, 11(10): e1005501.
[29] BAUTERS D, COBBAUT M, GEYS L, et al. Loss of ADAMTS5 enhances brown adipose tissue mass and promotes browning of white adipose tissue via CREB signaling[J]. Mol Metab, 2017, 6(7): 715–724.
[30] RIAZ-UL-HAQ, CHYTIL F. Expression of nuclear retinoic acid receptors in rat adipose tissue[J]. Biochem Biophys Res Commun, 1991, 176(3): 1539–1544.
[31] VILLARROYA F. Differential effects of retinoic acid on white and brown adipose tissues.An unexpected role for vitamin A derivatives on energy balance[J]. Ann N Y Acad Sci, 1998, 839: 190–195.
[32] FOURNIER T, PAVAN L, TARRADE A, et al. The role of PPAR-γ/RXR-α heterodimers in the regulation of human trophoblast invasion[J]. Ann N Y Acad Sci, 2002, 973(1): 26–30.
[33] FROMENT P, FABRE S, DUPONT J, et al. Expression and functional role of peroxisome proliferator-activated receptor-γ in ovarian folli-culogenesis in the sheep[J]. Biol Reprod, 2003, 69(5): 1665–1674.
[34] 李华振, 狄冉, 郭晓飞, 等. 绵羊繁殖轴相关组织RXRs基因的表达及RXRA基因多态性与季节性发情之间的关系[J]. 农业生物技术学报, 2019, 27(8): 1401–1409.
LI H Z, DI R, GUO X F, et al. Expression of RXRs genes in sheep (Ovis aries) reproductive axis related organizations and the relationship between RXRA gene polymorphism and seasonal estrus[J]. Journal of Agricultural Biotechnology, 2019, 27(8): 1401–1409. (in Chinese)
[35] ROHRER G A, NONNEMAN D J. Genetic analysis of teat number in pigs reveals some developmental pathways independent of vertebra number and several loci which only affect a specific side[J]. Genet Sel Evol, 2017, 49(1): 4.
[36] XIANG G H, REN J L, HAI T, et al. Editing porcine IGF2 regulatory element improved meat production in Chinese Bama pigs[J]. Cell Mol Life Sci, 2018, 75(24): 4619–4628. DOI: 10.1007/s00018-018-2917-6
[37] 谭光辉, 李杰章, 吴磊, 等. ITPR2基因内含子遗传变异与鸭蛋壳品质关联性分析[J]. 山地农业生物学报, 2020, 39(2): 40–45.
TAN G H, LI J Z, WU L, et al. The relationship between genetic variation of ITPR2 gene intron and duck eggshell quality of Arias platyrhyncha domesticca[J]. Journal of Mountain Agriculture and Biology, 2020, 39(2): 40–45. (in Chinese)
[38] 文禹粱, 王翔宇, 郭晓飞, 等. 不同产羔数小尾寒羊BMP4基因表达及其错义突变与产羔数关联分析[J]. 农业生物技术学报, 2019, 27(1): 80–88.
WEN Y L, WANG X Y, GUO X F, et al. BMP4 gene expression in small tail Han sheep (Ovis aries) with different litter size and association analysis between its missense mutation and litter size[J]. Journal of Agricultural Biotechnology, 2019, 27(1): 80–88. (in Chinese)
[39] LIU N J, ZHANG K, ZHAO H Y. Haplotype-association analysis[J]. Adv Genet, 2008, 60: 335–405.
[40] DROGEMULLER C, HAMANN H, DISTL O. Can-didate gene markers for litter size in different German pig lines[J]. J Anim Sci, 2001, 79(10): 2565–2570.
[41] 王围围, 肖玉净, 马云龙, 等. 猪3号染色体上与产仔性状和窝内仔猪初生重整齐度显著关联的标记鉴定[J]. 中国畜牧兽医, 2020, 47(11): 3546–3554.
WANG W W, XIAO Y J, MA Y L, et al. Identification of markers significantly associated with litter size traits and uniformity of the piglet's birth weight on porcine chromosome 3[J]. China Animal Husbandry & Veterinary Medicine, 2020, 47(11): 3546–3554. (in Chinese)