脊椎动物胚胎发育从一个单细胞到一个独立的个体,是通过高度复杂而又协调的机制在时间和空间上共同发挥作用的,其中一个关键阶段就是脊椎轴向骨骼的形成。在胚胎发育过程中,体节是轴向骨骼形成的前提,最初的体节是由中胚层分化形成,并在形态上是相似的,但它们具有沿体轴方向的解剖学特征。最终形成具有复杂特征的枕骨、颈椎、胸椎、腰椎、荐椎和尾椎[1-3]。在物种内每个域中的椎骨数量是固定的,但在物种之间每个域中的椎骨数量各不相同,从而产生了特定于物种的轴向公式[4]。而猪是为数不多的胸、腰椎数可变的经济动物之一[5]。猪胸、腰椎总数与胴体性状相关性较强,研究表明随着胸、腰椎数的增加,胴体长[6]和胴体重[7]也随之增加。
体节作为脊椎形成最重要的前体结构,在形成过程中受到多个基因及基因家族的调节,包括时钟基因Hes7[8]和Lfng[9],以及Hox基因家族[10]等。其中Homeobox(Hox)基因在组织特性规范中起关键作用,其表达受到严格的时间调控,赋予了轴向特性,还可以控制肢芽的位置,形成侧板中胚层,因此,Hox基因的表达指定了解剖学上不同的脊柱亚型[2]。Hox基因被认为在脊椎动物胚胎发生过程中参与前后轴(a-p)模式[11],在生物体的轴向发育模式中起着重要作用[12]。哺乳动物有将近40个Hox基因家族成员,分成4个簇:HoxA、HoxB、HoxC和HoxD,这些基因的形成主要是通过祖先簇的重复以及随后的基因丢失或重复引起[13]。Hox基因具有相似的序列和相对位置被划分为13个组,13组基因在簇内的位置具有功能相关性,反映了基因在前后轴上的时间和空间表达顺序,这一特征称为共线性[14]。脊椎动物的Hox簇在组中是紧凑的,具有统一的转录方向,共线Hox基因的激活被翻译成沿轴的空间表达模式[15]。其中HoxB簇基因在不同发育阶段和多个组织中的表达存在空间共线性[16]。HoxB簇基因共包含10个基因(HoxB1-9和HoxB13)。前期本团队对北京黑猪进行GWAS和GO富集分析将HoxB1-7、HoxB9和HoxB13作为影响脊椎数的候选基因[17]。
北京黑猪是1960年由巴克夏、大白和中国地方猪种杂交形成的一个培育猪种[18],在北方猪肉市场中占有重要地位,但是北京黑猪在Hox基因多态性与性状的关联分析方面未见报道。本研究对HoxB簇中9个基因所有外显子区进行PCR及测序检测突变位点并将其与表型进行关联分析,以期探索影响性状变异候选基因功能位点,为性状变异遗传机理的阐释奠定基础。
1 材料与方法 1.1 试验群体北京黑猪群体均饲养在北京黑六牧业科技有限公司,共251头个体,所有猪健康状况良好。管理方式和饲喂方式均一致。251头北京黑猪于220日龄在北京二商集团有限责任公司进行屠宰,在屠宰中收集其组织样本并记录脊椎数和胴体性状(胴体长、胴体斜长、胴体直长和胴体重)。
1.2 DNA提取及检测使用QIAamp DNA Mini Kit试剂盒进行组织样品的DNA提取,使用IMPLEN超微量分光光度计进行DNA质量检测(DNA浓度、核酸与蛋白及酚类物质最高吸收峰的吸收波长比值(OD260 nm/OD280 nm)、核酸与碳水化合物最高吸收峰的波长比值(OD260 nm/OD230 nm)的检测),并用2%琼脂糖凝胶进行DNA质检,检测合格的DNA样品用于后期试验。
1.3 引物设计及合成本试验的引物均使用Primer Premier 5.0软件设计,引物序列均由赛默飞(ThermoFisher)合成。本研究设计9个基因(HoxB1 (ENSSSCT0000 0019086.4)、HoxB2 (ENSSSCT00000019087.4)、HoxB3 (ENSSSCT00000019088.4)、HoxB4 (EN SSSCT00000019089.5)、HoxB5 (ENSSSCT00000 019094.5)、HoxB6 (ENSSSCT00000019093.5)、HoxB7 (ENSSSCT00000019092.4)、HoxB9 (ENS SSCT00000101474.1)、HoxB13 (ENSSSCT00-000019095.5))共45对引物,其中检测到SNPs的引物有6对(表 1)。
使用诺唯赞(Code:P505-d1,Vazyme)的高保真酶进行PCR,本试验对251头北京黑猪个体采用25 μL反应体系(2×Phanta Max Buffer 12.5 μL,0.5 μL聚合酶,0.5 μL dNTP,1 μL上游引物,1 μL下游引物,8.5 μL无菌水,1 μL模板)进行PCR。程序如下(ABI, Singapore):95 ℃预变性5 min;95 ℃变性30 s,57~62 ℃退火30 s,72 ℃延伸30~90 s,35个循环;72 ℃延伸2 min。PCR产物通过切胶回收单一条带在北京六合华大基因科技有限公司平台进行双向Sanger测序。
1.5 基因分型及数据统计分析将测序结果使用SeqMan进行基因分型,使用Microsoft Excel 2016进行基因型频率和等位基因频率的统计。采用SAS软件的GLM程序,以基因型为固定效应,采用单因素方差分析(ANOVA)分析基因型效应。动物模型如下: Y=μ+基因型+e,其中Y为表型值,μ为共同均值,e为随机误差。采用Duncan’s多重检验估计不同基因型间最小二乘均值(LSM)的差异。数据使用“平均值±标准误”表示,P < 0.05判定为差异显著。
2 结果 2.1 北京黑猪脊椎数和胴体性状的表型统计251头北京黑猪个体的表型值(脊椎数和胴体性状)见表 2,脊椎数的均值为20.95节,变异系数为2.67%;胴体长的均值为98 cm,变异系数为5.19%;胴体直长的均值为89.74 cm,变异系数为4.83%;胴体斜长的均值为76.06 cm,变异系数为4.38%;胴体重的均值为64.37 kg,变异系数为12.47%。
通过对HoxB簇中9个基因的mRNA序列进行变异位点检测共发现12个SNPs(表 3)均在UTR区,其中1个5′UTR突变(HoxB2),11个3′UTR突变分别在HoxB5、HoxB9、HoxB13中(表 3)。这12个突变位点的等位基因频率变化范围为4%~96%,其基因型频率变化范围为1%~93%。
本研究将HoxB簇中的9个基因通过PCR扩增及变异检测筛选到12个(1个HoxB2基因5′UTR,2个HoxB5基因3′UTR,2个HoxB9基因3′UTR,7个HoxB13基因3′UTR)SNPs与脊椎数和胴体性状进行关联分析,结果见表 4。HoxB2上的1个5′UTR突变c.40 C>T、HoxB5上的2个3′UTR突变c.1766 A>G、c.1767 A>G与性状关联不显著;HoxB9上的2个3′UTR突变中c.2743 C> T与胴体性状关联显著,其中CT基因型对应的胴体性状有优势,而c.2254 A>G与性状关联不显著;HoxB13上的7个3′UTR突变中c.1863 G>A、c.1440 A>C与脊椎数性状显著关联,两个位点中AA基因型对应的脊椎数更有优势,而c.1736 G> A、c.1444 G>C、c.1433 A>C、c.1310 A>G及c.1221 G>A与性状关联不显著。
猪的脊椎数作为一个高遗传力的有重要经济价值的优良性状,一直以来都受到研究者的高度青睐。研究者利用QTL定位和GWAS技术在猪的多个经济性状进行检测,发现了大量的QTLs[19]。猪QTL数据库共记录34 333个QTLs(http://www.animalgenome.org/cgi-bin/QTLdb/index,2021年12月27日),其中对于脊椎数的研究显示,QTLs主要集中在SSC1和SSC7上,而对于胴体长和胴体重每条染色体上都存在QTL。Mikawa等[20]通过对3个F2群体进行最小二乘区间映射分析,发现SSC1和SSC7上存在与脊椎数相关的QTL。大量的研究表明,影响脊椎数变异的QTL主要集中在SSC7上[21]。对于胴体性状也进行了大量的QTL定位和GWAS研究,发现QTL在SSC7上也有定位[22-24]。但在北京黑猪的GWAS结果中除了SSC7上出现了经典的QTL以外,在SSC12上出现了一个新的QTL,并且在QTL中定位到了HoxB簇的9个基因为影响脊椎数变异的候选基因[17]。
每簇的3′端Hox基因首先表达,而更多的5′端Hox基因在稍后时间表达,这种现象被称为“时间共线性”[25]。Hox基因在动物胚胎发育过程中调节轴向区域特征,最初在原肠胚形成期间被激活。HoxB簇包含10个基因,HoxB1-9和HoxB13。HoxB簇基因中缺失90 kb的小鼠表现出一系列有序的沿颈椎和胸椎柱的单节前部同源异位转化以及胸骨形态发生缺陷[26]。在鸡胚中,对HoxB簇基因的研究表明,HoxB1-3基因均在原始条纹附近表现出更高水平的表达,而在周围细胞中表达水平逐渐降低[16];HoxB5基因表达模式与其他HoxB簇基因显著不同,它始于后原始条纹,随后沿条纹向前延伸;HoxB8基因沿原始条纹和Hensen节点表达,偶尔以不对称方式表达;且HoxB9与HoxB8基因表达非常相似。研究表明,HoxB9基因在Alpha5亚基缺陷小鼠的胚胎中表达减少,Alpha5beta1纤连蛋白的互作对于中胚层衍生物的维持以及某些神经嵴细胞的存活至关重要,影响了中胚层的发育[27]。并且HoxB9基因还对体节发育中胸腔的形成以及前肢的发育有重要作用[28-30]。HoxB13基因在主动终止沿前后轴的后延伸中至关重要[31-34]。此外,HoxB13基因敲除对于小鼠体节的增加有显著影响[35]。同时将其启动子进行过表达,发现小鼠胚胎缺失大部分甚至全部尾巴[31]。因此,先前大量的研究均表明Hox家族在发育中发挥着重要作用,从原肠胚形成阶段开始,到胚胎前后轴的形成、肢体发育以及器官形成都起到关键作用[15, 30, 36-37]。因此,本研究对北京黑猪SSC12上HoxB簇的9个基因进行多态性分析,发现在mRNA的5′和3′UTRs存在SNPs。其中,有个通常高度保守的从几个到几百个核苷酸不等的3′UTRs[38]。储存在3′UTRs中的遗传信息可以通过3′UTRs介导的蛋白质互作(PPI)的形成传递给蛋白质[39]。高度保守的序列特征表明3′UTRs具有重要的调控作用,如mRNA的定位、mRNA的稳定性和翻译等占有重要地位;另一方面,3′UTRs存在异构体,可以用来区分高度相似的蛋白质功能[40]。
先前对于HoxB簇基因进行了多态性研究,对东亚人HoxB簇基因多态性与牙齿和咬合特征进行关联性分析[41]。对HoxB13基因的多态性进行探索也发现了与前列腺癌关联的SNPs[42]。对肺癌的研究中,在HoxB2基因中也发现了与其关联显著的SNPs[43]。在本研究中,发现3′UTRs存在与脊椎数和胴体性状显著相关的SNPs。本研究在HoxB9基因的3′UTR中通过关联分析发现1个SNP(c.2743 C>T)与胴体性状关联显著。因此,HoxB13基因对前后轴的完成有至关重要的作用,主要是在主动终止沿前AP轴的后延伸[31-34]。但在本研究中对其两个外显子进行PCR扩增及关联分析发现,3′UTR中2个SNPs(c.1863 G>A、c.1440 A>C)与脊椎数性状关联显著,与胴体性状关联不显著。
4 结论在北京黑猪群体中,HoxB9基因上的c.2743 C> T位点与胴体性状关联显著,HoxB13基因中c.1863 G>A、c.1440 A>C与脊椎数性状关联显著。这3个位点的发现为猪脊椎数和胴体性状的分子标记辅助选育工作提供了基础。
[1] |
SCAAL M. Early development of the vertebral column[J]. Semin Cell Dev Biol, 2016, 49: 83-91. DOI:10.1016/j.semcdb.2015.11.003 |
[2] |
WELDON S A, MVNSTERBERG A E. Somite development and regionalisation of the vertebral axial skeleton[J]. Semin Cell Dev Biol, 2022, 127: 10-16. DOI:10.1016/j.semcdb.2021.10.003 |
[3] |
WELDON S A, MVNSTERBERG A E. Somite development and regionalisation of the vertebral axial skeleton[J]. Semin Cell Dev Biol, 2022, 127: 10-16. DOI:10.1016/j.semcdb.2021.10.003 |
[4] |
IIMURA T, DENANS N, POURQUIÉ O. Establishment of Hox vertebral identities in the embryonic spine precursors[J]. Curr Top Dev Biol, 2009, 88: 201-234. |
[5] |
NARITA Y, KURATANI S. Evolution of the vertebral formulae in mammals: a perspective on developmental constraints[J]. J Exp Zool B Mol Dev Evol, 2005, 304(2): 91-106. |
[6] |
KING J W B, ROBERTS R C. Carcass length in the bacon pig; its association with vertebrae numbers and prediction from radiographs of the young pig[J]. Anim Prod, 1960, 2(1): 59-65. DOI:10.1017/S0003356100033493 |
[7] |
BURGOS C, LATORRE P, ALTARRIBA J, et al. Allelic frequencies of NR6A1 and VRTN, two genes that affect vertebrae number in diverse pig breeds: a study of the effects of the VRTN insertion on phenotypic traits of a Duroc×Landrace-Large White cross[J]. Meat Sci, 2015, 100: 150-155. DOI:10.1016/j.meatsci.2014.09.143 |
[8] |
ANDERSON M J, MAGIDSON V, KAGEYAMA R, et al. Fgf4 maintains Hes7 levels critical for normal somite segmentation clock function[J]. Elife, 2020, 9: e55608. DOI:10.7554/eLife.55608 |
[9] |
OKUBO Y, SUGAWARA T, ABE-KODUKA N, et al. Lfng regulates the synchronized oscillation of the mouse segmentation clock via trans-repression of Notch signalling[J]. Nat Commun, 2012, 3(1): 1141. DOI:10.1038/ncomms2133 |
[10] |
BÖHMER C. Correlation between Hox code and vertebral morphology in the mouse: towards a universal model for Synapsida[J]. Zoological Lett, 2017, 3: 8. DOI:10.1186/s40851-017-0069-4 |
[11] |
GODSAVE S, DEKKER E J, HOLLING T, et al. Expression patterns of Hoxb genes in the Xenopus embryo suggest roles in anteroposterior specification of the hindbrain and in dorsoventral patterning of the mesoderm[J]. Dev Biol, 1994, 166(2): 465-476. DOI:10.1006/dbio.1994.1330 |
[12] |
李桃桃, 金美林, 费晓娟, 等. Hox基因家族及其对动物脊椎形成的影响[J]. 畜牧兽医学报, 2022, 53(4): 999-1009. LI T T, JIN M L, FEI X J, et al. The Hox gene family and its effects on spine formation in animals[J]. Acta Veterinaria et Zootechnica Sinica, 2022, 53(4): 999-1009. (in Chinese) |
[13] |
KRUMLAUF R. Hox genes, clusters and collinearity[J]. Int J Dev Biol, 2018, 62(11-12): 659-663. DOI:10.1387/ijdb.180330rr |
[14] |
DURSTON A J. Some questions and answers about the role of Hox temporal collinearity in vertebrate axial patterning[J]. Front Cell Dev Biol, 2019, 7: 257. DOI:10.3389/fcell.2019.00257 |
[15] |
MALLO M, WELLIK D M, DESCHAMPS J. Hox genes and regional patterning of the vertebrate body plan[J]. Dev Biol, 2010, 344(1): 7-15. DOI:10.1016/j.ydbio.2010.04.024 |
[16] |
GOUVEIA A, MARCELINO H M, GONÇALVES L, et al. Patterning in time and space: HoxB cluster gene expression in the developing chick embryo[J]. Cell Cycle, 2015, 14(1): 135-145. DOI:10.4161/15384101.2014.972868 |
[17] |
NIU N, WANG H, SHI G, et al. Genome scanning reveals novel candidate genes for vertebral and teat number in the Beijing Black Pig[J]. Anim Genet, 2021, 52(5): 734-738. DOI:10.1111/age.13111 |
[18] |
ZHANG L C, WANG L G, LI Y, et al. A substitution within erythropoietin receptor gene D1 domain associated with litter size in Beijing Black pig, Sus scrofa[J]. Anim Sci J, 2011, 82(5): 627-632. DOI:10.1111/j.1740-0929.2011.00901.x |
[19] |
ERNST C W, STEIBEL J P. Molecular advances in QTL discovery and application in pig breeding[J]. Trends Genet, 2013, 29(4): 215-224. DOI:10.1016/j.tig.2013.02.002 |
[20] |
MIKAWA S, HAYASHI T, NII M, et al. Two quantitative trait loci on Sus scrofa chromosomes 1 and 7 affecting the number of vertebrae[J]. J Anim Sci, 2005, 83(10): 2247-2254. DOI:10.2527/2005.83102247x |
[21] |
FAN Y, XING Y Y, ZHANG Z Y, et al. A further look at porcine chromosome 7 reveals VRTN variants associated with vertebral number in Chinese and western pigs[J]. PLoS One, 2013, 8(4): e62534. DOI:10.1371/journal.pone.0062534 |
[22] |
QIAO R M, GAO J, ZHANG Z Y, et al. Genome-wide association analyses reveal significant loci and strong candidate genes for growth and fatness traits in two pig populations[J]. Genet Sel Evol, 2015, 47(1): 17. DOI:10.1186/s12711-015-0089-5 |
[23] |
GUO Y M, HUANG Y X, HOU L J, et al. Genome-wide detection of genetic markers associated with growth and fatness in four pig populations using four approaches[J]. Genet Sel Evol, 2017, 49(1): 21. DOI:10.1186/s12711-017-0295-4 |
[24] |
LI L Y, XIAO S J, TU J M, et al. A further survey of the quantitative trait loci affecting swine body size and carcass traits in five related pig populations[J]. Anim Genet, 2021, 52(5): 621-632. DOI:10.1111/age.13112 |
[25] |
DESCHAMPS J, VAN NES J. Developmental regulation of the Hox genes during axial morphogenesis in the mouse[J]. Development, 2005, 132(13): 2931-2942. DOI:10.1242/dev.01897 |
[26] |
MEDINA-MARTÍNEZ O, BRADLEY A, RAMÍREZ-SOLIS R. A large targeted deletion of Hoxb1-Hoxb9 produces a series single-segment anterior homeotic transformations[J]. Dev Biol, 2000, 222(1): 71-83. DOI:10.1006/dbio.2000.9683 |
[27] |
GOH K L, YANG J T, HYNES R O. Mesodermal defects and cranial neural crest apoptosis in alpha5 integrin-null embryos[J]. Development, 1997, 124(21): 4309-4319. DOI:10.1242/dev.124.21.4309 |
[28] |
CHEN F, CAPECCHI M R. Targeted mutations in Hoxa-9 and Hoxb-9 reveal synergistic interactions[J]. Dev Biol, 1997, 181(2): 186-196. DOI:10.1006/dbio.1996.8440 |
[29] |
MCINTYRE D C, RAKSHIT S, YALLOWITZ A R, et al. Hox patterning of the vertebrate rib cage[J]. Development, 2007, 134(16): 2981-2989. DOI:10.1242/dev.007567 |
[30] |
XU B, WELLIK D M. Axial Hox9 activity establishes the posterior field in the developing forelimb[J]. Proc Natl Acad Sci U S A, 2011, 108(12): 4888-4891. DOI:10.1073/pnas.1018161108 |
[31] |
AIRES R, DE LEMOS L, NÓVOA A, et al. Tail bud progenitor activity relies on a network comprising Gdf11, Lin28, and Hox13 genes[J]. Dev Cell, 2019, 48(3): 383-395. DOI:10.1016/j.devcel.2018.12.004 |
[32] |
AMIN S, NEIJTS R, SIMMINI S, et al. Cdx and T brachyury Co-activate growth signaling in the embryonic axial progenitor niche[J]. Cell Rep, 2016, 17(12): 3165-3177. DOI:10.1016/j.celrep.2016.11.069 |
[33] |
MALLO M. Reassessing the role of Hox genes during vertebrate development and evolution[J]. Trends Genet, 2018, 34(3): 209-217. DOI:10.1016/j.tig.2017.11.007 |
[34] |
STEVENTON B, MARTINEZ ARIAS A. Evo-engineering and the cellular and molecular origins of the vertebrate spinal cord[J]. Dev Biol, 2017, 432(1): 3-13. DOI:10.1016/j.ydbio.2017.01.021 |
[35] |
ECONOMIDES K D, ZELTSER L, CAPECCHI M R. Hoxb13 mutations cause overgrowth of caudal spinal cordand tail vertebrae[J]. Dev Biol, 2003, 256(2): 317-330. DOI:10.1016/S0012-1606(02)00137-9 |
[36] |
LAPPIN T R J, GRIER D G, THOMPSON A, et al. HOX GENES: seductive science, mysterious mechanisms[J]. Ulster Med J, 2006, 75(1): 23-31. |
[37] |
WELLIK D M. Hox genes and kidney development[J]. Pediatr Nephrol, 2011, 26(9): 1559-1565. DOI:10.1007/s00467-011-1902-1 |
[38] |
XIE X H, LU J, KULBOKAS E J, et al. Systematic discovery of regulatory motifs in human promoters and 3'UTRs by comparison of several mammals[J]. Nature, 2005, 434(7031): 338-345. DOI:10.1038/nature03441 |
[39] |
BERKOVITS B D, MAYR C. Alternative 3'UTRs act as scaffolds to regulate membrane protein localization[J]. Nature, 2015, 522(7556): 363-367. DOI:10.1038/nature14321 |
[40] |
MAYR C. What are 3' UTRs doing?[J]. Cold Spring Harb Perspect Biol, 2019, 11(10): a034728. DOI:10.1101/cshperspect.a034728 |
[41] |
YAMAGUCHI T, KAWAGUCHI A, KIM Y I, et al. The role of polymorphisms associated with early tooth eruption in dental and occlusal traits in East Asian populations[J]. Korean J Orthod, 2014, 44(2): 96-102. DOI:10.4041/kjod.2014.44.2.96 |
[42] |
SAUNDERS E J, DADAEV T, LEONGAMORNLERT D A, et al. Fine-mapping the HOXB region detects common variants tagging a rare coding allele: evidence for synthetic association in prostate cancer[J]. PLoS Genet, 2014, 10(2): e1004129. DOI:10.1371/journal.pgen.1004129 |
[43] |
CLEMENCEAU A, BOUCHERAT O, LANDRY-TRUCHON K, et al. Lung cancer susceptibility genetic variants modulate HOXB2 expression in the lung[J]. Int J Dev Biol, 2018, 62(11-12): 857-864. DOI:10.1387/ijdb.180210yb |
(编辑 郭云雁)