畜牧兽医学报  2022, Vol. 53 Issue (8): 2417-2434. DOI: 10.11843/j.issn.0366-6964.2022.08.002    PDF    
绵羊经济性状相关基因研究进展及其应用
李亭亭1, 刘秋月2, 李向臣1, 王海涛2     
1. 浙江农林大学 动物科技学院·动物医学院, 杭州 311300;
2. 中国科学院遗传与发育生物学研究所, 北京 100101
摘要:绵羊是重要的农业经济动物之一,不同绵羊品种具有各自独特的性状。利用性状存在变异的绵羊个体,对其进行基因组水平的分析,可以定位获得影响某种性状的基因区段或主效突变位点,鉴定得到的这些重要基因可为绵羊遗传育种提供分子标记,加快遗传选育进程。文章综述了在绵羊中鉴定得到的影响繁殖、生长、抗病及抗逆性状的主效基因发现过程及相应的调控机理,同时还综述了这些基因在绵羊分子育种方面的研究进展,有助于后续探索这些重要基因的具体分子机制,并为绵羊分子育种提供新思路。
关键词绵羊    经济性状    主效基因    分子标记    
Research Progress and Applications of Genes Associated with Economic Traits in Sheep
LI Tingting1, LIU Qiuyue2, LI Xiangchen1, WANG Haitao2     
1. College of Animal Science and Technology · College of Veterinary Medicine, Zhejiang A&F University, Hangzhou 311300, China;
2. Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing 100101, China
Abstract: Sheep is one of the important agricultural economic animals and their characteristics vary in different sheep breeds. Major putative regions and gene loci affecting economic traits can be mapped by genome-wide analysis in sheep breeds with phenotypic variation. These identified gene loci may be utilized as molecular markers to benefit sheep breeding and speed up the process of genetic selection. In this paper, genes related to reproduction, growth, disease resistance and stress tolerance identified from different sheep breeds were extensively reviewed. Furtherly, discovery process, regulatory mechanism and application progress in molecular breeding in sheep were also reviewed. These studies may contribute to further molecular mechanisms exploration of these major genes and provide new insights into molecular breeding of sheep.
Key words: sheep    economic traits    major gene    molecular marker    

绵羊(Ovis aries)起源于现亚洲中西部的新月沃地区,是世界上最早被人类驯化的家养动物之一[1]。根据联合国粮农组织(Food and Agriculture Organization of United Nations,FAO)的统计数据,我国是世界上最大的羊肉生产国,其产量几乎占世界总产量的三分之一[2]。随着人们对羊肉需求的增加,针对地方绵羊品种繁殖、生长以及抗病性状等方面的选育越来越受到重视,育种方式也越来越多样化。

绵羊在上万年的进化过程中经历了长期的人工和自然选择。在不同地区和环境中的绵羊群体为了适应自然环境或人类喜好,会对自身体型、毛色、繁殖、耐受等生理表型进行调整,还会在基因组上留下不同印记,通过对绵羊品种中不同印记进行检测,就能够对动物的环境适应性机制有更好的理解,同时如果能够找到影响这些表型的遗传变异,就能极大提升育种的效率。在绵羊上,目前研究者主要针对繁殖、生长、脂尾、抗病及抗逆等几个与经济效益关联紧密的性状,通过家系遗传定位、比较基因组学、候选基因关联分析以及全基因组关联分析等手段,先后鉴定获得绵羊中多个包括增加排卵数、提高生长速度、增加产肉量以及增强抗病能力的主效基因以及致因位点,如图 1所示。其中部分基因已经大范围运用于绵羊育种实践,大部分正处于动物个体水平验证的阶段。

图 1 绵羊重要经济性状及其相关基因 Fig. 1 Important economic traits and related genes of sheep

随着遗传修饰技术的进一步发展,一套新的育种策略和方法已经逐渐形成。研究者可以通过转基因或基因编辑等方式直接对获得的候选基因进行体内功能验证,极大的缩短了大动物中候选主效基因的验证周期。针对定位获得的候选基因,最普遍的做法就是通过CRISPR/Cas9系统将基因进行直接敲除,观察动物表型以确认基因功能[3-4];通过CRISPR/Cas9结合同源重组或供体DNA,可以在不完全破坏功能基因的前提下针对特定基因型进行定点编辑[5-6];更进一步,随着单碱基编辑系统的开发,目前研究者在无需依赖基因组产生双链断裂的情况下,可以进行高效定点的基因编辑[7-8],或者利用先导编辑(prime editor,PE)技术对基因进行插入、缺失等突变体的编辑以确认致因突变的价值[9]。本文将目前已研究的相对清楚并影响绵羊重要经济性状的相关基因和应用情况进行了总结,为以后的绵羊分子育种工作提供参考。

1 绵羊重要经济性状及相关基因 1.1 繁殖性状

绵羊的繁殖性状是影响养羊业最重要的经济性状,可体现在产羔和季节性发情两个方面。目前在绵羊上鉴定出的重要多羔主效基因的多个突变主要集中在变异生长因子β(transforming growth factor-β,TGF-β)通路的相关基因上,主要有FecBBMP15、GDF9,另外还包括近年来鉴定出的FecL基因,它是唯一一个不属于TGF-β信号通路的调控绵羊排卵数的基因。而对于影响绵羊季节性发情的研究,鉴定出的基因主要有TSHRMTNR1A,可通过绵羊对光照的反应来调控体内相关生殖激素的合成,从而影响发情。

1.1.1 产羔数性状   产羔数性状是养羊业经济性状中最重要的生产性状,受母羊排卵数、子宫容量等因素影响,其中排卵数可直接决定母羊的产羔数和繁殖力。母绵羊产羔数与排卵数的相关性达到0.96,但产羔数性状的遗传力较低,约为0.11[10],而排卵数性状遗传力相对较高,自然群体中约为0.19左右,而在选择群体中可高达0.30[11]。雌性绵羊在出生前,卵巢上便形成了大量原始卵泡。初情期前,卵泡虽能发育,但不能成熟排卵,当发育到一定程度时,便闭锁退化。初情期后,在每个发情周期中可发育的卵泡多达几十个,但单羔品种最终能发育到成熟排卵的卵泡一般只有1~2个,其余卵泡中途闭锁而死亡,而多羔品种则有多个卵泡能发育成熟并排卵。

目前发现的多羔主效基因大部分都属于TGF-β信号通路成员。FecB基因(fecundity booroola)是第一个被发现与绵羊产羔数有关的主效基因,1994年,Montgomery等[12]利用微卫星标记和遗传连锁的方法将与Booroola Merino绵羊多产相关的FecB基因定位到绵羊6号染色体上,之后发现该基因实际为绵羊骨形态发生蛋白受体1B型基因(bone morphogenetic protein receptor type 1B,BMPR1B) 编码区发生了A746G突变,造成其编码的氨基酸由谷氨酰胺变为精氨酸(Q249R),该位点的突变减弱了BMPR1B在BMP信号通路中的配体功能,导致卵母细胞的加速成熟并排卵,使母羊排卵数增加[13]。经过大量研究证实,FecB突变在Booroola Merino绵羊、Javanese Indonesia、小尾寒羊、湖羊、Garole绵羊和Kalehkoohi等多产绵羊品种中均有发现[14]。TGF-β家族中分别位于绵羊X和5号染色体上的骨形态发生蛋白15(bone morphogenetic protein 15,BMP15)和生长分化因子9(growth differentiation factor 9,GDF9)基因,是由欧洲科学家利用自然突变的高产羔数绵羊群体,构建遗传连锁图谱,利用微卫星扫描或者全基因组关联分析(genome-wide association study,GWAS)定位获得的影响绵羊产羔数的主效基因,二者有高度的氨基酸序列同源性,通常共同形成异源二聚体作用于BMPRⅡ型受体,并激活Smad2/3(sma- and mad-related,Smad)和Smad 1/5/8信号通路,促进颗粒细胞的增殖及分化、卵泡发育以及增加排卵数[15]。迄今为止,在不同绵羊品种中已识别出BMP15基因的10个与排卵率、繁殖力相关的点突变,分别具有不同的功能,即FecXHFecXG(终止密码子提前使蛋白合成受损),FecXRFecXBar(通过改变阅读框),FecXL(导致BMP15蛋白功能丧失),FecXIFecXB(与GDF9蛋白相互作用受损),FecXGrFecXO(细胞信号活性改变)以及FecXN(未知)[16]。与其他突变不同的是,FecXNFecXGrFecXO的纯合突变不仅不会导致母羊不育,反而使其具有更高的繁殖能力。在GDF9中发现了11个突变,分别是G1~G8,以及FecGTFecGEFecGA[17],这些突变均可对排卵和产羔造成影响。目前,基因编辑技术已经运用到大动物模型上,研究者们利用CRISPR/Cas9技术获得了具有g.A746G点突变的FecB绵羊[18-19],但有关这些突变羔羊繁殖能力的研究还未见报道。吴艳芳[20]不仅获得了突变的滩羊,而且还证明具有FecB突变的滩羊比野生型滩羊具有更高的产羔数,这些研究为获得FecB突变的绵羊从而提高绵羊产羔数提供了一个新的方向。而对于BMP15和GDF9,可能是由于其效应位点较多难以选择或一些绵羊群体中存在自然突变,很少有关于其在绵羊上的报道,但是有研究人员在猪上引入BMP15无义突变后,观察到纯合突变导致母猪不育,而杂合突变则可增加母猪排卵数[21],这一发现与存在自然突变的绵羊群体中观察到的现象一致。

上述3个与多羔相关的主效基因中,多个突变位点的效应已经在世界范围内多个绵羊品种中得到证实。近年来研究者还利用其它绵羊群体陆续在这3个基因中发现了新的突变,可为绵羊产羔性状研究提供新的参考[22-27]。另外,研究者还通过转录组及代谢组等分析手段,发现在这3个基因发生突变的个体中,与胰岛素分泌、cAMP、MAPK、Smads、Hippo等相关的信号通路被显著富集[28-29],且突变个体卵泡液中的葡萄糖6-磷酸、天冬氨酸、半胱氨酸、谷胱甘肽等重要代谢物与排卵率的增加显著相关,这些研究结果将有助于深入阐明绵羊多羔性状的发生机制[30]

除TGF-β信号通路之外的产羔数主效基因也已有报道。FecL基因是法国农科院Bodin等[11]从1975到1996年间通过5个世代的选育构建Lacaune绵羊的回交家系,又经过基因组扫描、高通量测序以及分子生物学验证,于2013年定位到了B4GALNT2(glycosylation enzyme beta-1, 4-N-acetyl galactosaminyltransferase 2)为FecL的主效基因[31-32],可能的机制是突变导致了B4GALNT2基因在卵泡中的异位表达,此蛋白在颗粒细胞中的高表达会使卵泡中特定靶蛋白的糖基化状态改变,进而引起了排卵数的改变。这一发现揭开了除TGF-β信号通路之外调控绵羊排卵数的全新机制。除了FecL基因外,最新的研究也鉴定到一些新的与绵羊产羔相关的重要基因,并发现了这些基因上存在的影响产羔数的多个突变位点,例如FSHRLEPRHIRACOIL等基因[33],这些功能基因突变之后,会对繁殖相关的生物学过程产生影响,如FSHR可通过影响排卵过程而导致产羔数的变化[34]LEP及其受体LEPR可参与卵泡的成熟和发育,影响胚胎附着和排卵过程[35]HIRA可能通过参与卵母细胞的发育而影响产羔数[36]COIL则可影响绵羊卵巢和卵泡发育[37]。另外,研究者还利用产羔数显著差异的绵羊个体进行比较生物学研究,从转录组和蛋白层面鉴定到了多个相关基因及信号通路,如DLG4、GSTK1、RRM2B、STMN1等可分别参与昼夜节律、GSH代谢和胰岛素分泌合成等相关激素调节过程[38]PRLPOMCGNRH1能影响GnRH活性和生殖激素的释放[39]FOXO4可调节卵泡发育和卵母细胞成熟[40]FLT3与KILTG能协同促进卵泡成熟和排卵;NLRP5可与卵泡发育相关的蛋白质相互作用,参与卵泡生成[41],这些结果使研究者对影响绵羊繁殖性状的机制多了新的认识。

1.1.2 发情性状   季节性发情是限制绵羊繁殖能力的瓶颈因素之一,经过研究者们几十年的探索,现在基本确定了哺乳动物中光照通过下丘脑-垂体-性腺轴(hypothalamic-pituitary-gonadal axis,HPGA)系统引起季节性发情,其机制为:眼球视网膜接受光照刺激产生视神经冲动,并由下视丘传至下丘脑视交叉上核(suprachiasmatic nucleus,SCN),SCN作为动物体内最主要的生物钟,将信号传至松果体,引起松果体中褪黑素(melatonin,MEL)分泌的改变,之后褪黑素通过HPGA轴进一步引起下丘脑中促性腺激素释放激素(gonadotropin releasing hormone,GnRH)以及相关生殖激素的合成、释放[42],最终调控动物的季节性发情活动。近年来,促甲状腺激素受体(thyroid-stimulating hormone receptor,TSHR)基因和褪黑素受体1A(melatonin receptor 1A,MTNR1A)基因被认为是影响绵羊季节性发情的重要调控基因,可作为改良绵羊季节性发情的重要分子标记。

TSHR最初在狗的甲状腺中被发现,之后通过对生活在不同环境下的野生和家养品种进行选择分析,发现TSHR在鸡和绵羊的驯化过程中均受到强烈选择[43-44],并且可通过HPGA轴通路反馈调节季节性繁殖活动[45]。目前认为的调控方式是下丘脑中的促甲状腺激素β亚基(thyrotrophin β subunit,TSHβ)可与TSHR结合,通过cAMP通路调控四碘甲状腺原氨酸(tetraiodothyronine,T4)与三碘甲状腺原氨酸(triodothyronine,T3)的转换[46],进而调节哺乳动物对光照周期的应答与季节性发情。研究发现, 在Soay绵羊中光照长短能通过调节TSHR的表达从而影响绵羊的季节性繁殖活动[47]。通过对国内多个地方绵羊品种的重测序以及分群验证,证实该基因的多态性与中国本土绵羊的季节性发情行为和繁殖活动显著相关[48-50]MTNR1A基因于1996年首先在西伯利亚仓鼠中被发现[51],参与介导褪黑素对哺乳动物生殖活动的影响,随后在Merinos d’Arles、Aragonesa、Sarda、小尾寒羊等绵羊品种中陆续发现MTNR1A基因与绵羊卵巢的排卵活动和产羔频率相关[52-53]。目前,研究者发现绵羊MTNR1A基因的不同突变不仅与绵羊的季节性发情和产羔有关[53-55],同时某些突变还对雄性绵羊的繁殖活动(如外阴嗅探、爬跨等)产生促进作用[56],被公认是调节哺乳动物繁殖活动的重要基因。

1.2 生长性状

畜禽生长性能包括生长速度、活体背膘厚及饲料转化率等,而其中生长速度可体现在肌细胞的增殖发育以及骨骼的生长速度两个方面。一方面,肌肉的形成和再生需通过一个高度精密的调控过程,使成肌细胞分化和融合成多核合胞体。在胚胎发育过程中,中胚层前体细胞受到周围组织信号的控制分化为肌源细胞,其中需要包括Pax7、MyoD和Myf5等因子的上调[57],增殖的前体细胞/成肌细胞停止细胞周期并诱导肌肉特异性基因表达[58],之后通过一系列高度有序的程序先发生肌细胞之间的融合形成新生的肌管,随后募集和融合额外的成肌细胞产生多核肌管,最终形成成熟的肌肉纤维。另一方面,骨骼系统包括硬骨、软骨、肌腱和韧带等,共同起源于中胚层和神经嵴产生的间充质干细胞[59],进一步分化为构成特定骨骼类型的细胞,这个过程受多种信号通路和转录因子的调节,如FGF、IGF1、MAPK、PI3K/Akt通路都可参与软骨细胞增殖和分化的调控[60],软骨形成后一部分作为永久性软骨,另一部分可作为生长板形成的模板,这些模板可形成钙化骨,成为纵向生长的骨骼,BMP信号通路在此过程中起着关键作用。在软骨中,肥大的软骨细胞产生Ihh(indian hedgehog),作用于生长板增厚带附近的软骨膜细胞,诱导成骨细胞分化[61],成骨细胞成熟后,产生骨基质蛋白,包括I型胶原蛋白和骨钙素,使骨骼有一定的弹性和硬度。

生长激素(growth hormone,GH)基因是研究较多的促进生长发育的基因,也是线性骨生长的主要调节剂,可通过全身和局部胰岛素样生长因子1(insulin-like growth factor 1,IGF1)的产生来刺激骨骼生长,使机体的生长速度更快。分子育种技术运用之初,GH就被通过转基因的方式导入绵羊体内获得生长更快的个体[62-64],但是近年来研究较少。

肌肉生长抑制素(myostatin,MSTN),也称为生长和分化因子8(growth and differentiation factor 8,GDF8),位于绵羊2号染色体,主要表达于骨骼肌中,在肌肉分化和生长过程中起负调控作用,通过部分抑制成肌细胞的增殖来调节肌肉质量[65],不同畜种中的“双肌现象”已被证明均是MSTN基因突变的结果。研究表明,MSTN可降低肌肉发育相关基因Pax7和MyoD的表达[66-67],阻断成肌细胞的细胞周期,使其停滞在G1期,抑制成肌细胞分化[68]。因此,当MSTN基因敲除之后,蛋白功能受损,则会产生相反的效果,促进肌细胞的增殖和分化,导致双肌性状。迄今为止,已经在多个物种中发现MSTN的突变与肌肉质量增加相关,在不同绵羊品种中也发现了多个突变位点[69-71],其中大多数位于该基因的非编码区,这表明位于编码区之外的遗传变异在调节肌肉发育方面有重要作用。研究表明,绵羊MSTN基因3′UTR区的SNP可将其表达水平降低约三分之一,且在显著增加肌肉占比的同时不影响肉品质[72]。近年来,随着CRISPR/Cas9技术的广泛应用,研究者已利用原核注射或体细胞核移植的方式获得MSTN成功敲除的绵羊个体[73-74],这些基因编辑个体与野生型绵羊相比,体重显著增加,骨骼肌纤维肥大,具有双肌性状。利用基因编辑的方法能够快速验证候选基因在绵羊中的效应并评价其运用效果,对绵羊产业的未来发展具有重要经济价值。

细胞因子信号转导抑制因子2(suppressors of cytokine signaling,SOCS2)是SOCS家族的一员,作为GH信号通路和骨骼生长的负调节因子,可通过SH2结构域与生长激素受体(growth hormone receptor,GHR)结合,对JAK/STAT5b通路起到抑制作用,影响机体的生长、骨骼发育以及新陈代谢等,表型上可对动物的骨骼线性增长和体重增加产生抑制作用[75-76]。当SOCS2发生突变时,导致蛋白结构发生改变,功能丧失,与受体的结合也被破坏,从而调节下游的GH/IGF1通路,对动物的生长产生促进作用。研究表明,SOCS2-/-小鼠体型可比野生型大30%~40%,而且这种体重增加是由于内脏器官、骨骼长度以及身体长度的显著增加引起的[77-78]。Rupp等[79]于2015年首次在Lacaune绵羊中定位到与绵羊体型大小、体重及产奶量显著相关的QTL,区段中包含SOCS2基因,当编码SOCS2蛋白SH2结构域的区段发生p.R96C点突变时,其编码的蛋白结构也发生改变,导致其与GHR的结合能力丧失,与野生型绵羊相比,p.R96C纯合突变的绵羊的体型、体重分别增加了24%和18%。之后,Zhou等[80]对绵羊SOCS2基因进行单碱基编辑,将BE3 mRNA和sgRNA共同注射到绵羊受精卵中,实现了SOCS2中p.R96C的单碱基替换,获得具有点突变的绵羊,证明该基因突变对于绵羊体重和体型有促进作用。

除以上提到的基因以及Callipyge[81]VRTN[82]Carwell[83]基因都可以影响绵羊的生长性状,Callipyge基因与MSTN有类似的效应,也是目前研究人员公认的与绵羊臀部肌肉发育有关的基因,VRTN是影响绵羊肋骨数的重要分子标记,而Carwell能够影响绵羊的眼肌宽度、深度及面积,这些分子标记都可对以后的绵羊育种提供重要的指导。

1.3 脂尾性状

绵羊尾部脂肪是肥尾绵羊储存能量的重要组织,反映了绵羊对恶劣和极端环境的适应性,尾部堆积的大量脂肪不仅可以增强绵羊在极端环境下的耐寒能力,也可作为迁徙和食物缺乏期间维持机体所需的重要能量来源,对于人类来说,它还可以在干旱和饥荒期间为人类提供高能量的食物[84-85]。脂肪细胞来源于间充质干细胞,尾部脂肪的沉积过程包括脂肪细胞的增殖、分化和成熟阶段。间充质干细胞首先分化为前体脂肪细胞,再转化为具有脂质运输和合成、分泌脂肪细胞特异性蛋白质等功能的成熟脂肪细胞[86-87],之后沉积在内脏器官、皮下或尾部,可为机体维持正常的生命活动提供能量。这些过程由复杂的基因网络之间相互协调,涉及PPAR、MAPK、Wnt/β-catenin、BMP/SMADs、胰岛素/IGF1等相关通路[88-89],多与脂质形成和代谢相关[90]。然而,近几十年来,随生活水平的提高和饮食习惯等的改变,绵羊大脂肪尾巴作为能量储备的优势逐渐减少,生产者和消费者对低脂肉类更感兴趣,形成相同质量的尾部脂肪比生成等量的肌肉组织需要付出更大的饲养成本,且相同质量的肥尾的价格远低于瘦肉价格,肥尾还会对母羊繁殖能力产生负面影响,过多的尾部脂肪沉积极大的降低了养羊业的经济效益。因此,探索调节绵羊尾部脂肪沉积和发育的候选基因及遗传分子机制对培育尾脂较少的绵羊品种具有重要意义。

通过基因组学方法,已经定位到数个与绵羊尾部脂肪沉积相关的候选基因,研究较多的为血小板衍生生长因子D(platelet derived growth factor D,PDGFD)和骨形态发生蛋白2(bone morphogenetic protein 2,BMP2)基因,在不同的绵羊脂尾相关研究中被多次提到[91-97],证明BMP2和PDGFD基因是影响绵羊肥尾表型的主要因素。其中PDGFD在全世界不同尾型的绵羊中选择信号最强,PDGFD通过激活PDGFRβ来调节PPARγ2和C/EBPα的表达,抑制白色脂肪细胞分化[97],也可以通过促进成熟脂肪细胞的形成或增殖来影响脂质代谢(如脂质合成),造成绵羊尾部的脂肪沉积[96]BMP2是脂肪发生的重要调节因子,不仅能诱导干细胞分化为脂肪细胞,还可能有助于将能量储存分配到内脏和皮下脂肪组织中[98]。Zhu等[99]的研究表明, BMP2的错义突变(Ala>Val)与绵羊肥尾表型直接相关,并参与尾部脂肪沉积。最新研究也表明, BMP2和PDGFD基因的SNP可能与阿勒泰绵羊的尾脂沉积相关,但这些基因的具体作用机制还不清晰,后期可作为重要的候选基因继续研究和验证。

近年来有多项全基因组关联研究通过利用地方绵羊品种定位到了多个影响脂尾性状的候选区域,其中,GLIS1、PPARGPPP1CC、LOC101117953等基因在绵羊尾型性状中受到强烈选择[95, 100-102]GLIS1作为促脂肪因子,可能在胎儿发育过程中对中胚层细胞的分化起关键作用,并影响绵羊尾部的脂肪沉积[102]PPARG(也称为PPARγ)参与脂肪的发育和功能,是控制脂肪生成、分化和发育的关键转录因子,在脂肪细胞分化的早期阶段就已开始发挥作用[89]PPP1CC位于BMP2附近,与BMP2有很强的相关性[93]LOC101117953为PPP1CC的反转录本(r-PPP1CC),但由于大多数反转录拷贝被认为是没有功能的,所以猜测其可能与PPP1CC相似,通过影响上游的BMP2来发挥作用[94]。目前, 对绵羊脂尾性状的研究仍主要集中在BMP2和PDGFD两个基因上,研究者们也逐渐发现了其他有效应的突变位点,但还需进一步对其进行功能验证和作用机制的探索。除此之外,一些新的相关基因和通路也逐渐被发现,也可作为研究其作用机制的参考。

1.4 抗病性状

在畜牧业养殖中,疫病常为养殖业带来巨大的经济损失。而在绵羊中,感染源通常包括细菌、病毒、寄生虫等,防治办法通常有添加抗生素、疫苗预防或治疗,但都会存在相应的问题,比如产生耐药株及疫苗失效等。不同品种的绵羊在对疾病的抵抗力方面存在差异,一旦通过现代基因组学方法找到与某些疾病抗性相关的主效基因以及确定相关突变,就可以更深入地了解与抗性相关的分子基础和遗传变异,结合现代基因操作技术就有可能获得稳定的绵羊抗病新品系。动物的抗病性状的是病原体在与宿主免疫系统对抗过程中被打败的结果,其抗性基因的定位大多也集中在病原体识别、免疫因子的分泌、抗感染及相关通路方面。

线虫是寄生虫的一种,而胃肠道线虫(gastrointestinal nematode,GIN)是对绵羊产业危害最大的一种线虫,胃肠道线虫病是影响全球畜牧生产系统的主要健康问题,动物感染之后可在死亡率、生产力和治疗成本方面给全世界造成重大经济损失,澳大利亚绵羊产业中每年约有4.3亿美元用于与线虫感染相关的生产损失和治疗成本[103],是世界各地绵羊和其他放牧反刍动物的最重要的健康疾病之一。在过去的几十年中,养羊业越来越依赖驱虫药物作为一种控制寄生虫病的方法,然而,这种方式不仅花费巨大又通常不是很有效,还会导致耐药株的出现。该病对反刍动物的影响主要反映在体重减轻、腹泻、厌食、贫血和死亡等方面,而幼龄羔羊又特别容易受到线虫感染,这是影响绵羊生产的主要经济负担[104],因此,大多数关于GIN抗性的QTL逐渐被定位出来[105]。Al Kalaldeh等[106]使用GWAS发现位于绵羊2号染色体上的PALLDGALNTL6基因与寄生虫抗性之间存在很强的关联,MASP2、DIS3L2也被发现与内部胃肠道线虫感染的抗性相关[107-108]。另外,绵羊对线虫的抗性及恢复能力也可能与介导病原体识别(TLRCLRNLR基因家族)、先天免疫反应(MASPNLRTLRIL20R)以及对感染的适应性反应(CLRIL2R、OLA-DRA)方面的基因有关[109],这些结果都可为研究绵羊抗胃肠道线虫感染提供参考。

布鲁氏菌病由布鲁氏菌(Brucella)感染所致,它是对反刍动物养殖和人类健康危害极大的人畜共患病,该病主要对动物的生殖系统造成危害,不仅能造成母羊流产、不孕不育,也会造成雄性动物生殖障碍,严重影响牲畜生产性能。科研人员对于布鲁氏菌病的抗性基因也进行了很多的挖掘,目前也有一些收获。Li等[110]通过构建褪黑素合成酶ASMT(acetylserotonin O-methyltransferase)基因过表达载体,成功获得具有抗布病的ASMT转基因绵羊,并证明ASMT基因过表达绵羊可能通过影响免疫相关的信号通路和肠道微生物群来增强对布鲁氏菌病的抵抗力。另一团队通过全基因组重测序分析,确定了4种与绵羊布鲁氏菌病易感性相关的重要途径以及相关的候选基因,即粘附小带(CTNNA3、PARD3和PTPRM)、细胞粘附分子(NLGN1、CNTNAPNCAM1和PTPRM)、唾液分泌(LOC101102109、PRKG1和ADCY2)和Hippo信号通路(CTNNA3、YAP1、PARD3)[111],这些基因与细胞粘附过程密切相关,可能有助于机体发挥免疫功能的第一道屏障。

目前,病原菌疾病相关的研究很少是在建立动物模型的基础上来开展的,大多数仍然是在群体中通过抗病个体来寻找可能的相关基因,或直接在发病个体或细胞水平研究其致病机制,从而来研制治疗方法,原因可能是目前大家对这些疾病深入的致病机制仍不清晰,导致即使想要构建动物疾病模型却无从下手。

1.5 羊毛性状

绵羊毛作为一种天然的蛋白纤维在绵羊养殖业和纺织业中都占有重要的经济地位,也是毛用绵羊的重要经济性状,其品质和质量可体现在羊毛的长度、细度、强度等方面。毛发生长发生在毛囊发育的生长期,始于真皮乳头部周围的基底膜干细胞,随着细胞系分化的开始而停止分裂,之后细胞的形状和位置停止变化,纤维细胞系开始合成角蛋白(keratin,K),皮层逐渐角质化,非角蛋白细胞成分逐渐分解,在上层的毛囊中,被皮脂腺和汗腺包裹的富含油脂的纤维从体表长出[112],形成毛发。羊毛纤维主要是由角蛋白组成的角质层,羊毛角蛋白包括含量较稳定的角蛋白中间丝蛋白(intermediate filament proteins,IFPs)和影响羊毛特性的角蛋白联合蛋白(keratin associate proteins,KAPs)两种[113-114],而角蛋白基因的表达可影响羊毛的弯曲、粗细、重量等重要经济特征。

绵羊的羊毛性状是受遗传、环境和营养等因素影响的复杂生理生化特性,许多信号通路和相关因素都参与到这些经济特征的调节中。成纤维细胞生长因子是一系列在促进和调节成纤维细胞生长中起重要作用的因子,其中成纤维细胞生长因子5(fibroblast growth factor 5, FGF5)基因已被证实是毛发生长的抑制剂,沉默FGF5基因的表达会延长毛发周期的生长期VI期,导致极长毛发的表型[115],而通过靶向使FGF5功能丧失可抑制真皮细胞的活化和毛发纤维的合成,促进羊毛的生长和增加羊毛产量[116-117]。除了FGF5基因,影响绵羊毛色的刺鼠信号蛋白(agouti signaling protein,ASIP)基因也是影响羊毛质量的重要基因,ASIP在哺乳动物中可抑制黑色素皮质素受体1(melanocortin 1 receptor,MC1R)信号转导和黑色素的产生,调控毛发的颜色,Zhang等[3]利用基因编辑技术破坏其功能,使绵羊的毛色发生改变。另外,角蛋白相关蛋白基因的表达及其多态性对绵羊羊毛直径、韧性、弯曲程度的影响也已在多项研究中被报道[118-119]

1.6 抗逆性状

极端气候地区的动物品种通过长期的自然选择及人工选择,在生理、生化、形态等方面获得了稳定的遗传特征,以适应高海拔、低氧、寒冷、干旱或高温等恶劣环境。因此,极端环境中具有特殊环境适应性的物种是挖掘抗逆机制以及候选基因的资源宝库,如果能在全基因组水平上了解其遗传基础并阐明这些分子标记,解析动物不同适应性表型的遗传基础,可为未来气候环境变化背景下制定适当的育种计划提供方向。

1.6.1 高海拔适应性状   生活在高海拔地区的动物为了适应低氧带来的挑战,在生理、细胞水平都发生了适应性变化,比如机体通过降低代谢活性、增加呼吸速率和红细胞总量来维持氧气运输[120],以满足动物维持其正常生长、发育和繁殖所需的组织氧气水平。此外,高原动物还会启动线粒体氧化代谢和信号转导通路,引起一系列的生理变化,使动物能够在缺氧条件下生存,而缺氧诱导因子(hypoxia-inducible factor,HIF)通路及相关因子在这其中发挥了重要的作用[121]。藏羊是一种耐低氧物种,生活在极其荒凉的高海拔环境中,与低海拔地区相比,该地区紫外线辐射强、温度低、氧分压低。藏羊除了提供优质的羊毛和肉类外,还具有抗粗饲和适应高海拔缺氧等优点,可作为一种对极端环境具有较强适应性的模型,从中发掘一些有助于指导育种实践的关键基因。目前已有大量关于高海拔适应性的研究报道,例如藏獒(EPAS1、EGLN1)、牦牛(ADAM17、ARG2和MMP3)、藏羚羊(PKLRATP12ANOS3)和藏野猪(EPAS1、CYP4F2和THSD7A)等[122-123],这些基因功能涉及缺氧、能量代谢、血管生成、Ca2+代谢、促进红细胞生成等方面。

利用高原上的低氧来完成正常的生理功能,这是机体在低氧的恶劣环境中生存的一种有效方式。近些年的研究报道了许多与低氧适应性相关的候选基因,其中研究最为透彻的基因是EPAS1和EGLN1,这两个候选基因都是在HIF通路上游起作用的关键基因[124-125]EPAS1(endothelial PAS domain protein 1,EPAS1),又称为HIF2α,通过对不同海拔地区的山羊进行外显子组测序和基因分型分析,发现EPAS1多态性在藏绒山羊中显示出强烈的选择信号[126-127],而Q579L突变对藏绒山羊适应高海拔环境起了关键作用。在绵羊上,Wei等[128]发现,EPAS1基因3′UTR的一个SNP与藏绵羊适应高海拔、促进红细胞生成有关,另外,该研究也发现了一些与代谢相关的基因,如DPP4能维持葡萄糖稳定性,参与糖酵解过程;PPARG与磷脂代谢相关;SOD1可激活ATP生成并减少线粒体ROS的产生。绵羊HIF的靶基因ZEB1以及位于7号染色体的FGF7基因在高海拔地区的绵羊群体中受到显著的选择[129-130],可能参与高海拔条件下细胞环境中可用氧的调节或其他调节。对于EGLN1基因,在绵羊中还未见到有关其与高海拔适应性相关的报道。

高海拔地区的生物除了要面对高海拔带来的缺氧和寒冷问题之外,还要应对紫外线强度对身体带来的损伤,因此紫外线信号通路扮演着重要的角色。在哺乳动物中,高强度的紫外线辐射可诱导大量ROS的产生,导致DNA损伤、细胞凋亡、组织损伤和癌症,破坏细胞内环境[131]。Zhang等[132]对藏鸡的研究发现,与炎症和免疫反应以及与紫外辐射导致的DNA损伤刺激和DNA修复过程相关的基因受到选择。在对绵羊的研究中,发现藏羊的HBBMITF基因[133]、印度肥尾羊的ERCC3和TGM3基因受到紫外线照射方面的选择[134]。从以上研究可以看出,目前参与高强度紫外线辐射适应性的研究较少,发现的相关基因也有限,并且相互之间的相似性很低。

除了以上提到的相关基因外,高海拔适应相关的FGF5、MTORTRHDETUBB3基因也受到强烈选择[135],但还需要进一步的研究和功能验证。

1.6.2 抗寒冷性状   由于中国北方高海拔地区有着漫长而严酷的高寒冬季,绵羊经常出现体温过低的情况,因此,抵御寒冷的能力对其生存至关重要。为了在处于低温环境时维持体温,绵羊必须增加机体的产热(生热),这包括颤抖生热(shivering thermogenesis,ST)和非颤抖生热(non-shivering thermogenesis,NST),一是通过骨骼肌的快速收缩水解ATP释放热量[136],二是通过中枢神经系统释放去甲肾上腺素,激活棕色脂肪细胞表面的受体或募集富含解偶联蛋白1(uncoupling protein,UCP1)的产热脂肪细胞分解甘油三酯产生热量[137]。据报道,同一品种内不同品系的绵羊抗寒能力存在差异,而绵羊抗寒能力增强的遗传力约为0.27[138]

在与抗寒相关的研究中,UCP家族通常受到强烈选择。该家族主要包括UCP1、UCP2和UCP3,三者之间有高度同源性,UCP1主要存在于线粒体内膜中,是棕色脂肪的标志物[139]UCP2在全身各个组织中均有表达,UCP3主要分布于骨骼肌和心组织中[140]。先前的报道称,该家族基因与不同家畜的寒冷适应性有关,存在于绵羊15号染色体上的UCP1基因参与寒冷诱导的非颤抖产热[141]UCP2和UCP3基因在绵羊骨骼肌中的表达对于其适应寒冷发挥着重要作用[140, 142]。通过将印度的一种高海拔绵羊品种Changthangi与分布在平原地区的绵羊品种相比较,15号染色体UCP基因家族存在最显著信号[143],说明该基因可能是参与绵羊适应寒冷气候环境的重要分子标记。

1.6.3 耐高温性状   从春季到夏季,气温和环境湿度均升高的综合影响会导致农业动物的热应激。热应激对动物的饲料消耗、生长速度和繁殖效率影响比较明显。据统计,热应激对全球畜牧生产的经济影响超过12亿美元[144]。研究认为,热应激会诱导ROS的产生和脂质过氧化,可造成DNA和蛋白质的损伤,而ROS的产生被认为是诱导热休克蛋白基因(heat shock protein,HSP)表达的关键因素[145]。热应激对动物健康有严重影响,它可能导致牲畜的多器官衰竭,降低其生产力和繁殖效率,严重时可造成死亡[146-147],但是在沙漠地区的炎热气候下,仍有保持正常生长的绵羊品种,这表明在绵羊上存在耐热性的遗传基础。

目前发现,与绵羊热应激相关的基因大多都是属于HSP家族或相关成员。热休克蛋白,也被称为分子伴侣,可由各种应激诱导表达,被认为是细胞应激的保护剂,可以介导降解受损细胞内的蛋白折叠,这种保护机制使细胞能够在高温条件下生存。此外,HSPs也能促进细胞的生长和分化,调节细胞凋亡,而且还和免疫功能相关。HSPs的激活主要是通过热休克转录因子(heat shock factor 1,HSF1)实现的[148],HSF1是一种响应热应激而被激活的转录因子。白绵羊中发现,HSP90AA1(HSP90伴侣)基因的多态性可有效克服热应激并保持正常的生产性能[149]HSPA8基因中的T177P多态性也与Arabi绵羊对炎热气候的高度耐热适应性有关[150]。另外,多浪羊与中国本土其他绵羊品种遗传距离较大,具有良好的抗逆性和耐热性[91],Wang等[151]发现DNAJB5(dnaJ heat shock protein family (Hsp40) member B5)基因多态性可能参与多浪羊的热应激,DNAJB5为DNAJ热休克蛋白40家族伴侣蛋白的一个成员,在免疫细胞的热应激耐受性中起着至关重要的作用[152]

2 总结与展望

很久以前人们就观察到某一性状在不同个体中的差别,但是由于当时技术的局限性,只能通过传统的杂交育种方式获得所想要的性状,而对其背后隐藏的生物学机制发掘进展缓慢。随着测序技术的飞速发展,基于基因组重测序以及GWAS等遗传定位研究方法的普遍运用,针对绵羊经济性状主效基因的研究也较之前有较大进展。目前,鉴定出的部分绵羊经济性状主效基因及效应影响如表 1所示,其中部分基因的作用机制已比较明晰,也成功构建了动物模型对其功能进行验证。另外通过遗传定位的方式获得了一些候选基因,但还需要进一步阐明机制和验证其功能。再加上近年来基因编辑技术的迅猛发展,这为快速建立基于基因编辑的动物模型,验证候选基因在个体水平的效应提供了解决方案,也为未来建立具有优良性状的绵羊新品系提供参考。

表 1 绵羊重要经济性状主效基因及致因突变(Oar_v4.0(GCF_000298735.2)) Table 1 Major genes and causative mutations of important economic traits in sheep(Oar_v4.0 (GCF_000298735.2))
参考文献
[1]
ZEDER M A. Domestication and early agriculture in the Mediterranean Basin: origins, diffusion, and impact[J]. Proc Natl Acad Sci USA, 2008, 105(33): 11597-11604. DOI:10.1073/pnas.0801317105
[2]
WANG H H, ZHANG L, CAO J X, et al. Genome-wide specific selection in three domestic sheep breeds[J]. PLos One, 2015, 10(6): e0128688. DOI:10.1371/journal.pone.0128688
[3]
ZHANG X M, LI W R, LIU C X, et al. Alteration of sheep coat color pattern by disruption of ASIP gene via CRISPR Cas9[J]. Sci Rep, 2017, 7(1): 8149. DOI:10.1038/s41598-017-08636-0
[4]
NIU Y, JIN M, LI Y, et al. Biallelic β-carotene oxygenase 2 knockout results in yellow fat in sheep via CRISPR/Cas9[J]. Anim Genet, 2017, 48(2): 242-244. DOI:10.1111/age.12515
[5]
CHANG Y Y, SHAO J J, GAO Y, et al. Reporter gene knock-in into Marc-145 cells using CRISPR/Cas9-mediated homologous recombination[J]. Biotechnol Lett, 2020, 42(8): 1317-1325. DOI:10.1007/s10529-020-02860-x
[6]
MENCHACA A, DOS SANTOS-NETO P C, SOUZA-NEVES M, et al. Otoferlin gene editing in sheep via CRISPR-assisted ssODN-mediated homology directed repair[J]. Sci Rep, 2020, 10(1): 5995. DOI:10.1038/s41598-020-62879-y
[7]
KOMOR A C, KIM Y B, PACKER M S, et al. Programmable editing of a target base in genomic DNA without double-stranded DNA cleavage[J]. Nature, 2016, 533(7603): 420-424. DOI:10.1038/nature17946
[8]
GAUDELLI N M, KOMOR A C, REES H A, et al. Publisher correction: programmable base editing of A·T to G·C in genomic DNA without DNA cleavage[J]. Nature, 2018, 559(7714): E8. DOI:10.1038/s41586-018-0070-x
[9]
ANZALONE A V, RANDOLPH P B, DAVIS J R, et al. Search-and-replace genome editing without double-strand breaks or donor DNA[J]. Nature, 2019, 576(7785): 149-157. DOI:10.1038/s41586-019-1711-4
[10]
WALDRON D F, THOMAS D L. Increased litter size in Rambouillet sheep: I.Estimation of genetic parameters[J]. J Anim Sci, 1992, 70(11): 3333-3344. DOI:10.2527/1992.70113333x
[11]
BODIN L, SANCRISTOBAL M, LECERF F, et al. Segregation of a major gene influencing ovulation in progeny of Lacaune meat sheep[J]. Genet Sel Evol, 2002, 34(4): 447-464. DOI:10.1186/1297-9686-34-4-447
[12]
MONTGOMERY G W, LORD E A, PENTY J M, et al. The booroola fecundity (FecB) gene maps to sheep chromosome 6[J]. Genomics, 1994, 22(1): 148-153. DOI:10.1006/geno.1994.1355
[13]
CHEN S, GUO X F, HE X Y, et al. Transcriptome analysis reveals differentially expressed genes and long non-coding RNAs associated with fecundity in sheep hypothalamus with different FecB genotypes[J]. Front Cell Dev Biol, 2021, 9: 633747. DOI:10.3389/fcell.2021.633747
[14]
KUMAR S, RAJPUT P K, BAHIRE S V, et al. Differential expression of BMP/SMAD signaling and ovarian-associated genes in the granulosa cells of FecB introgressed GMM sheep[J]. Syst Biol Reprod Med, 2020, 66(3): 185-201. DOI:10.1080/19396368.2019.1695977
[15]
MOTTERSHEAD D G, SUGIMURA S, AL-MUSAWI S L, et al. Cumulin, an oocyte-secreted heterodimer of the transforming growth factor-β family, is a potent activator of granulosa cells and improves oocyte quality[J]. J Biol Chem, 2015, 290(39): 24007-24020. DOI:10.1074/jbc.M115.671487
[16]
CHANTEPIE L, BODIN L, SARRY J, et al. Genome-wide identification of a regulatory mutation in BMP15 controlling prolificacy in sheep[J]. Front Genet, 2020, 11: 585. DOI:10.3389/fgene.2020.00585
[17]
BRAVO S, LARAMA G, PAZ E, et al. Polymorphism of the GDF9 gene associated with litter size in Araucana creole sheep[J]. Anim Genet, 2016, 47(3): 390-391. DOI:10.1111/age.12404
[18]
ZHOU S W, DING Y G, LIU J, et al. Highly efficient generation of sheep with a defined FecBB mutation via adenine base editing[J]. Genet Sel Evol, 2020, 52(1): 35. DOI:10.1186/s12711-020-00554-6
[19]
ZHOU S W, YU H H, ZHAO X E, et al. Generation of gene-edited sheep with a defined Booroola fecundity gene (FecBB) mutation in bone morphogenetic protein receptor type 1B (BMPR1B) via clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated (Cas) 9[J]. Reprod Fertil Dev, 2018, 30(12): 1616-1621. DOI:10.1071/RD18086
[20]
吴艳芳. MSTN基因敲除和FecB基因突变滩羊扩繁试验[D]. 杨凌: 西北农林科技大学, 2021.
WU Y F. Breeding experiment of MSTN gene knockout tan sheep and FecB gene mutation tan sheep[D]. Yangling: Northwest A&F University, 2021. (in Chinese)
[21]
FLOSSMANN G, WURMSER C, PAUSCH H, et al. A nonsense mutation of bone morphogenetic protein-15 (BMP15) causes both infertility and increased litter size in pigs[J]. Bmc Genomics, 2021, 22(1): 38. DOI:10.1186/s12864-020-07343-x
[22]
ÇELIKELOǦLU K, TEKERLI M, ERDOǦAN M, et al. An investigation of the effects of BMPR1B, BMP15, and GDF9 genes on litter size in Ramlιç and Daǧlιç sheep[J]. Arch Anim Breed, 2021, 64(1): 223-230. DOI:10.5194/aab-64-223-2021
[23]
NAJAFABADI H A, KHANSEFID M, MAHMOUD G G, et al. Identification of sequence variation in the oocyte-derived bone morphogenetic protein 15 (BMP15) gene (BMP15) associated with litter size in New Zealand sheep (Ovisaries) breeds[J]. Mol Biol Rep, 2021, 48(9): 6335-6342. DOI:10.1007/s11033-021-06627-z
[24]
NAJAFABADI H A, KHANSEFID M, MAHMOUD G G, et al. Identification of polymorphisms in the oocyte-derived growth differentiation growth factor 9 (GDF9) gene associated with litter size in New Zealand sheep (Ovisaries) breeds[J]. Reprod Domest Anim, 2020, 55(11): 1585-1591. DOI:10.1111/rda.13813
[25]
WANG F Y, CHU M X, PAN L X, et al. Polymorphism detection of GDF9 gene and its association with litter size in Luzhong mutton sheep (Ovisaries)[J]. Animals (Basel), 2021, 11(2): 571.
[26]
WEN Y L, GUO X F, MA L, et al. The expression and mutation of BMPR1B and its association with litter size in small-tail Han sheep (Ovisaries)[J]. Arch Anim Breed, 2021, 64(1): 211-221. DOI:10.5194/aab-64-211-2021
[27]
LI H X, XU H W, AKHATAYEVA Z, et al. Novel indel variations of the sheep FecB gene and their effects on litter size[J]. Gene, 2021, 767: 145176. DOI:10.1016/j.gene.2020.145176
[28]
LI Z F, HE X Y, ZHANG X S, et al. Analysis of expression profiles of circRNA and miRNA in oviduct during the follicular and luteal phases of sheep with two fecundity (FecB gene) genotypes[J]. Animals (Basel), 2021, 11(10): 2826.
[29]
JIA J L, JIN J P, CHEN Q, et al. Eukaryotic expression, Co-IP and MS identify BMPR-1B protein-protein interaction network[J]. Biol Res, 2020, 53(1): 24. DOI:10.1186/s40659-020-00290-7
[30]
GUO X F, WANG X Y, DI R, et al. Metabolic effects of FecB gene on follicular fluid and ovarian vein serum in sheep (Ovis aries)[J]. Int J Mol Sci, 2018, 19(2): 539. DOI:10.3390/ijms19020539
[31]
DROUILHET L, TARAGNAT C, FONTAINE J, et al. Endocrine characterization of the reproductive axis in highly prolific Lacaune sheep homozygous for the FecLL mutation[J]. Biol Reprod, 2010, 82(5): 815-824. DOI:10.1095/biolreprod.109.082065
[32]
DROUILHET L, MANSANET C, SARRY J, et al. The highly prolific phenotype of lacaune sheep is associated with an ectopic expression of the B4GALNT2 gene within the ovary[J]. PLoS Genet, 2013, 9(9): e1003809. DOI:10.1371/journal.pgen.1003809
[33]
GOOTWINE E. Invited review: opportunities for genetic improvement toward higher prolificacy in sheep[J]. Small Ruminant Res, 2020, 186: 106090. DOI:10.1016/j.smallrumres.2020.106090
[34]
WANG W M, LIU S J, LI F D, et al. Polymorphisms of the ovine BMPR-IB, BMP-15 and FSHR and their associations with litter size in two Chinese indigenous sheep breeds[J]. Int J Mol Sci, 2015, 16(5): 11385-11397.
[35]
JUENGEL J L, FRENCH M C, O'CONNELL A R, et al. Mutations in the leptin receptor gene associated with delayed onset of puberty are also associated with decreased ovulation and lambing rates in prolific Davisdale sheep[J]. Reprod Fertil Dev, 2016, 28(9): 1318-1325. DOI:10.1071/RD14382
[36]
ZHOU M, PAN Z Y, CAO X H, et al. Single nucleotide polymorphisms in the HIRA gene affect litter size in small tail han sheep[J]. Animals (Basel), 2018, 8(5): 71.
[37]
MA H Y, FANG C, LIU L L, et al. Identification of novel genes associated with litter size of indigenous sheep population in Xinjiang, China using specific-locus amplified fragment sequencing technology[J]. Peer J, 2019, 7: e8079. DOI:10.7717/peerj.8079
[38]
LI C Y, HE X Y, ZHANG Z J, et al. Pineal gland transcriptomic profiling reveals the differential regulation of lncRNA and mRNA related to prolificacy in STH sheep with two FecB genotypes[J]. BMC Genomic Data, 2021, 22(1): 9. DOI:10.1186/s12863-020-00957-w
[39]
ZHANG Z B, TANG J S, DI R, et al. Integrated hypothalamic transcriptome profiling reveals the reproductive roles of mRNAs and miRNAs in sheep[J]. Front Genet, 2020, 10: 1296. DOI:10.3389/fgene.2019.01296
[40]
XU S S, GAO L, SHEN M, et al. Whole-genome selective scans detect genes associated with important phenotypic traits in sheep (Ovis aries)[J]. Front Genet, 2021, 12: 738879. DOI:10.3389/fgene.2021.738879
[41]
CHEN S, TAO L, HE X Y, et al. Single-nucleotide polymorphisms in FLT3, NLRP5, and TGIF1 are associated with litter size in Small-tailed Han sheep[J]. Arch Anim Breed, 2021, 64(2): 475-486. DOI:10.5194/aab-64-475-2021
[42]
EBLING F J P. Photoperiodic regulation of puberty in seasonal species[J]. Mol Cell Endocrinol, 2010, 324(1-2): 95-101. DOI:10.1016/j.mce.2010.03.018
[43]
RUBIN C J, ZODY M C, ERIKSSON J, et al. Whole-genome resequencing reveals loci under selection during chicken domestication[J]. Nature, 2010, 464(7288): 587-591. DOI:10.1038/nature08832
[44]
FARIELLO M I, SERVIN B, TOSSER-KLOPP G, et al. Selection signatures in worldwide sheep populations[J]. PLoS One, 2014, 9(8): e103813. DOI:10.1371/journal.pone.0103813
[45]
DARDENTE H, WYSE C A, BIRNIE M J, et al. A molecular switch for photoperiod responsiveness in mammals[J]. Curr Biol, 2010, 20(24): 2193-2198. DOI:10.1016/j.cub.2010.10.048
[46]
KARLSSON AC, FALLAHSHAHROUDI A, JOHNSEN H, et al. A domestication related mutation in the thyroid stimulating hormone receptor gene (TSHR) modulates photoperiodic response and reproduction in chickens[J]. Gen Comp Endocrinol, 2016, 228: 69-78. DOI:10.1016/j.ygcen.2016.02.010
[47]
HANON E A, LINCOLN G A, FUSTIN J M, et al. Ancestral TSH mechanism signals summer in a photoperiodic mammal[J]. Curr Biol, 2008, 18(15): 1147-1152. DOI:10.1016/j.cub.2008.06.076
[48]
窦立静, 赵赓, 高磊, 等. TSHR基因T315A位点在不同繁殖特性绵羊群体中的多态性研究[J]. 家畜生态学报, 2016, 37(4): 9-15.
DOU L J, ZHAO G, GAO L, et al. Study on distribution of T315A locus of TSHR gene in sheep populations with different reproductive characteristics[J]. Acta Ecologae Animalis Domastici, 2016, 37(4): 9-15. DOI:10.3969/j.issn.1673-1182.2016.04.003 (in Chinese)
[49]
夏青, 狄冉, 刘艳琴, 等. 绵羊TSHR基因表达及其多态性与季节性发情之间的关系[J]. 家畜生态学报, 2020, 41(3): 15-20.
XIA Q, DI R, LIU Y Q, et al. Expression and polymorphism of TSHR gene and its association with seasonal estrus in sheep[J]. Acta Ecologae Animalis Domastici, 2020, 41(3): 15-20. DOI:10.3969/j.issn.1673-1182.2020.03.003 (in Chinese)
[50]
轩俊丽, 马晓萌, 王慧华, 等. 绵羊季节性繁殖相关基因TSHR外显子多态性研究[J]. 畜牧兽医学报, 2016, 47(7): 1342-1353.
XUAN J L, MA X M, WANG H H, et al. Study on exon polymorphism of seasonal breeding related gene TSHR in sheep[J]. Acta Veterinaria et Zootechnica Sinica, 2016, 47(7): 1342-1353. (in Chinese)
[51]
WEAVERD R, LIU C, REPPERTS M. Nature's knockout: the Mel1b receptor is not necessary for reproductive and circadian responses to melatonin in Siberian hamsters[J]. Mol Endocrinol, 1996, 10(11): 1478-1487.
[52]
HE X Y, ZHANG Z B, LIU Q Y, et al. Polymorphisms of the melatonin receptor 1A gene that affects the reproductive seasonality and litter size in Small Tail Han sheep[J]. Reprod Domest Anim, 2019, 54(10): 1400-1410. DOI:10.1111/rda.13538
[53]
LURIDIANA S, COSSO G, PULINAS L, et al. New polymorphisms at MTNR1A gene and their association with reproductive resumption in sarda breed sheep[J]. Theriogenology, 2020, 158: 438-444. DOI:10.1016/j.theriogenology.2020.10.006
[54]
MURA M C, LURIDIANA S, PULINAS L, et al. Reproductive response to male joining with ewes with different allelic variants of the MTNR1A gene[J]. Anim Reprod Sci, 2019, 200: 67-74. DOI:10.1016/j.anireprosci.2018.11.012
[55]
CALVO J H, SERRANO M, MARTINEZ-ROYO A, et al. SNP rs403212791 in exon 2 of the MTNR1A gene is associated with reproductive seasonality in the Rasa aragonesa sheep breed[J]. Theriogenology, 2018, 113: 63-72. DOI:10.1016/j.theriogenology.2018.02.013
[56]
ABECIA J A, MURA M C, CARVAJAL-SERNA M, et al. Polymorphisms of the melatonin receptor 1A (MTNR1A) gene influence the age at first mating in autumn-born ram-lambs and sexual activity of adult rams in spring[J]. Theriogenology, 2020, 157: 42-47. DOI:10.1016/j.theriogenology.2020.07.030
[57]
KABLAR B, ASAKURA A, KRASTEL K, et al. MyoD and Myf-5 define the specification of musculature of distinct embryonic origin[J]. Biochem Cell Biol, 1998, 76(6): 1079-1091. DOI:10.1139/o98-107
[58]
HASTY P, BRADLEY A, MORRIS J H, et al. Muscle deficiency and neonatal death in mice with a targeted mutation in the myogeningene[J]. Nature, 1993, 364(6437): 501-506. DOI:10.1038/364501a0
[59]
KOBAYASHI T, KRONENBERG H M. Overview of skeletal development[M]//HILTON M J. Skeletal Development and Repair. New York: Humana, 2021: 3-16.
[60]
KAWAMOTO T, KAWAMOTO K. Preparation of thin frozen sections from nonfixed and undecalcified hard tissues using kawamot's film method (2012)[M]//HILTON M J. Skeletal Development and Repair: Methods and Protocols. Totowa: Humana Press, 2014.
[61]
LONG F X, CHUNG U I, OHBA S, et al. Ihh signaling is directly required for the osteoblast lineage in the endochondral skeleton[J]. Development, 2004, 131(6): 1309-1318. DOI:10.1242/dev.01006
[62]
HAMMER RE, PURSEL V G, REXROADC E Jr, et al. Production of transgenic rabbits, sheep and pigs by microinjection[J]. Nature, 1985, 315(6021): 680-683. DOI:10.1038/315680a0
[63]
MURRAYJ D, NANCARROWC D, MARSHALLJ T, et al. Production of transgenic merino sheep by microinjection of ovine metallothionein-ovine growth hormone fusion genes[J]. Reprod Fertil Dev, 1989, 1(2): 147-155. DOI:10.1071/RD9890147
[64]
ADAMS N R, BRIEGEL J R, WARD K A. The impact of a transgene for ovine growth hormone on the performance of two breeds of sheep[J]. J Anim Sci, 2002, 80(9): 2325-2333.
[65]
MCPHERRON AC, GUO T, WANG Q, et al. Soluble activin receptor type ⅡB treatment does not cause fat loss in mice with diet-induced obesity[J]. Diabet Obes Metab, 2012, 14(3): 279-282. DOI:10.1111/j.1463-1326.2011.01520.x
[66]
HAN D S, HUANG H P, WANG T G, et al. Transcription activation of myostatin by trichostatin A in differentiated C2C12myocytes via ASK1-MKK3/4/6-JNK and p38 mitogen-activated protein kinase pathways[J]. J Cell Biochem, 2010, 111(3): 564-573. DOI:10.1002/jcb.22740
[67]
LANGLEY B, THOMAS M, BISHOP A, et al. Myostatin inhibits myoblast differentiation by down-regulating MyoD expression[J]. J Biol Chem, 2002, 277(51): 49831-49840. DOI:10.1074/jbc.M204291200
[68]
CHEN M M, ZHAO Y P, ZHAO Y, et al. Regulation of myostatin on the growth and development of skeletal muscle[J]. Front Cell Dev Biol, 2021, 9: 785712. DOI:10.3389/fcell.2021.785712
[69]
GROCHOWSKA E, BORYS B, LISIAK D, et al. Genotypic and allelic effects of the myostatin gene (MSTN) on carcass, meat quality, and biometric traits in Colored Polish Merino sheep[J]. Meat Sci, 2019, 151: 4-17. DOI:10.1016/j.meatsci.2018.12.010
[70]
GROCHOWSKA E, BORYS B, MROCZKOWSKI S. Effects of intronic SNPs in the myostatin gene on growth and carcass traits in colored polish merino sheep[J]. Genes (Basel), 2019, 11(1): 2. DOI:10.3390/genes11010002
[71]
OSMAN N M, SHAFEY H I, ABDELHAFEZ M A, et al. Genetic variations in the Myostatin gene affecting growth traits in sheep[J]. Vet World, 2021, 14(2): 475-482. DOI:10.14202/vetworld.2021.475-482
[72]
KIJAS J W, MCCULLOCH R, EDWARDS J E H, et al. Evidence for multiple alleles effecting muscling and fatness at the OvineGDF8 locus[J]. BMC Genet, 2007, 8: 80.
[73]
LI H H, WANG G, HAO Z Q, et al. Generation of biallelic knock-out sheep via gene-editing and somatic cell nuclear transfer[J]. Sci Rep, 2016, 6: 33675. DOI:10.1038/srep33675
[74]
CRISPO M, MULET A P, TESSON L, et al. Efficient generation of myostatin knock-out sheep using CRISPR/Cas9 technology and microinjection into zygotes[J]. PLoS One, 2015, 10(8): e0136690. DOI:10.1371/journal.pone.0136690
[75]
DOBIE R, MACRAE V E, HUESA C, et al. Direct stimulation of bone mass by increased GH signalling in the osteoblasts of Socs2-/- mice[J]. J Endocrinol, 2014, 223(1): 93-106. DOI:10.1530/JOE-14-0292
[76]
GREENHALGH C J, BERTOLINO P, ASA S L, et al. Growth enhancement in suppressor of cytokine signaling 2 (SOCS-2)-deficient mice is dependent on signal transducer and activator of transcription 5b (STAT5b)[J]. Mol Endocrinol, 2002, 16(6): 1394-1406. DOI:10.1210/mend.16.6.0845
[77]
METCALF D, GREENHALGH C J, VINEY E, et al. Gigantism in mice lacking suppressor of cytokine signalling-2[J]. Nature, 2000, 405(6790): 1069-1073. DOI:10.1038/35016611
[78]
PASS C, MACRAE V E, HUESA C, et al. SOCS2 is the critical regulator of GH action in murine growth plate chondrogenesis[J]. J Bone Miner Res, 2012, 27(5): 1055-1066. DOI:10.1002/jbmr.1544
[79]
RUPP R, SENIN P, SARRY J, et al. A point mutation in suppressor of cytokine signaling 2 (Socs2) increases the susceptibility to inflammation of the mammary gland while associated with higher body weight and size and higher milk production in a sheep model[J]. PLoS Genet, 2015, 11(12): e1005629. DOI:10.1371/journal.pgen.1005629
[80]
ZHOU S W, CAI B, HE C, et al. Programmable base editing of the sheep genome revealed no genome-wide off-target mutations[J]. Front Genet, 2019, 10: 215. DOI:10.3389/fgene.2019.00215
[81]
MA D Y, YU Q Q, HEDRICK V E, et al. Proteomic and metabolomic profiling reveals the involvement of apoptosis in meat quality characteristics of ovine M.longissimus from different callipyge genotypes[J]. Meat Sci, 2020, 166: 108140. DOI:10.1016/j.meatsci.2020.108140
[82]
LI C Y, LI M, LI X Y, et al. Whole-genome resequencing reveals loci associated with thoracic vertebrae number in sheep[J]. Front Genet, 2019, 10: 674. DOI:10.3389/fgene.2019.00674
[83]
COCKETT N E, SMIT M A, BIDWELL C A, et al. The callipyge mutation and other genes that affect muscle hypertrophy in sheep[J]. Genet Sel Evol, 2005, 37(Suppl 1): S65-S81. DOI:10.1186/1297-9686-37-S1-S65
[84]
BAKHTIARIZADEH M R, SALAMI S A. Identification and expression analysis of long noncoding RNAs in fat-tail of sheep breeds[J]. G3 (Bethesda), 2019, 9(4): 1263-1276. DOI:10.1534/g3.118.201014
[85]
MORADI M H, NEJATI-JAVAREMI A, MORADI-SHAHRBABAK M, et al. Genomic scan of selective sweeps in thin and fat tail sheep breeds for identifying of candidate regions associated with fat deposition[J]. BMC Genet, 2012, 13: 10.
[86]
ROSEN E D, MACDOUGALD O A. Adipocyte differentiation from the inside out[J]. Nat Rev Mol Cell Biol, 2006, 7(12): 885-896. DOI:10.1038/nrm2066
[87]
DAVIS G H, MONTGOMERY G W, ALLISON A J, et al. Segregation of a major gene influencing fecundity in progeny of Booroola sheep[J]. New Zeal J Agr Res, 1982, 25(4): 525-529. DOI:10.1080/00288233.1982.10425216
[88]
BAKHTIARIZADEH M R, SALEHI A, ALAMOUTI A A, et al. Deep transcriptome analysis using RNA-Seq suggests novel insights into molecular aspects of fat-tail metabolism in sheep[J]. SciRep, 2019, 9(1): 9203.
[89]
HU G, WANG S Z, WANG Z P, et al. Genetic epistasis analysis of 10 peroxisome proliferator-activated receptorγ-correlated genes in broiler lines divergently selected for abdominal fat content[J]. Poult Sci, 2010, 89(11): 2341-2350. DOI:10.3382/ps.2010-00857
[90]
ZHANG W, XU M S, WANG J J, et al. Comparative transcriptome analysis of key genes and pathways activated in response to fat deposition in two sheep breeds with distinct tail phenotype[J]. Front Genet, 2021, 12: 639030. DOI:10.3389/fgene.2021.639030
[91]
WEI C H, WANG H H, LIU G, et al. Genome-wide analysis reveals population structure and selection in Chinese indigenous sheep breeds[J]. BMC Genomics, 2015, 16(1): 194. DOI:10.1186/s12864-015-1384-9
[92]
MOIOLI B, PILLA F, CIANI E. Signatures of selection identify loci associated with fat tail in sheep[J]. J Anim Sci, 2015, 93(10): 4660-4669. DOI:10.2527/jas.2015-9389
[93]
YUAN Z, LIU E, LIU Z, et al. Selection signature analysis reveals genes associated with tail type in Chinese indigenous sheep[J]. Anim Genet, 2017, 48(1): 55-66. DOI:10.1111/age.12477
[94]
PAN Z Y, LI S D, LIU Q Y, et al. Rapid evolution of a retro-transposable hotspot of ovine genome underlies the alteration of BMP2 expression and development of fat tails[J]. BMC Genomics, 2019, 20(1): 261. DOI:10.1186/s12864-019-5620-6
[95]
ZHAO F P, DENG T Y, SHI L Y, et al. Genomic scan for selection signature reveals fat deposition in Chinese indigenous sheep with extreme tail types[J]. Animals (Basel), 2020, 10(5): 773.
[96]
DONG K Z, YANG M, HAN J G, et al. Genomic analysis of worldwide sheep breeds reveals PDGFD as a major target of fat-tail selection in sheep[J]. BMC Genomics, 2020, 21(1): 800. DOI:10.1186/s12864-020-07210-9
[97]
LI X, YANG J, SHEN M, et al. Whole-genome resequencing of wild and domestic sheep identifies genes associated with morphological and agronomic traits[J]. Nat Commun, 2020, 11(1): 2815. DOI:10.1038/s41467-020-16485-1
[98]
GUIU-JURADO E, UNTHAN M, BÖHLER N, et al. Bone morphogenetic protein 2 (BMP2) may contribute to partition of energy storage into visceral and subcutaneous fat depots[J]. Obesity, 2016, 24(10): 2092-2100. DOI:10.1002/oby.21571
[99]
ZHU C Y, LI N, CHENG H P, et al. Genome wide association study for the identification of genes associated with tail fat deposition in Chinese sheep breeds[J]. Biol Open, 2021, 10(5): bio054932. DOI:10.1242/bio.054932
[100]
FAN H Y, HOU Y L, SAHANA G, et al. A transcriptomic study of the tail fat deposition in two types of Hulun Buir Sheep according to tail size and sex[J]. Animals (Basel), 2019, 9(9): 655.
[101]
ZHU C Y, FAN H Y, YUAN Z H, et al. Genome-wide detection of CNVs in Chinese indigenous sheep with different types of tails using ovine high-density 600K SNP arrays[J]. Sci Rep, 2016, 6: 27822. DOI:10.1038/srep27822
[102]
LUO R S, ZHANG X R, WANG L K, et al. GLIS1, a potential candidate gene affect fat deposition in sheep tail[J]. Mol Biol Rep, 2021, 48(5): 4925-4931. DOI:10.1007/s11033-021-06468-w
[103]
LANE J, SHEPHARD R, WEBB-WARE J, et al. Priority list of endemic diseases for the red meat industries[R]. Final Report, North Sydney: Meat&Livestock Australia Limited, 2015.
[104]
TESFAYE T. Prevalence, species composition, and associated risk factors of small ruminant gastrointestinal nematodes in South Omo zone, South-western Ethiopia[J]. J Adv Vet Anim Res, 2021, 8(4): 597-605.
[105]
ATLIJA M, ARRANZ J J, MARTINEZ-VALLADARES M, et al. Detection and replication of QTL underlying resistance to gastrointestinal nematodes in adult sheep using the ovine 50K SNP array[J]. Genet Sel Evol, 2016, 48: 4. DOI:10.1186/s12711-016-0182-4
[106]
AL KALALDEH M, GIBSON J, LEE S H, et al. Detection of genomic regions underlying resistance to gastrointestinal parasites in Australian sheep[J]. Genet Sel Evol, 2019, 51(1): 37. DOI:10.1186/s12711-019-0479-1
[107]
SAYRE B L, HARRIS G C. Systems genetics approach reveals candidate genes for parasite resistance from quantitative trait loci studies in agricultural species[J]. Anim Genet, 2012, 43(2): 190-198. DOI:10.1111/j.1365-2052.2011.02231.x
[108]
BECKER G M, DAVENPORT K M, BURKE J M, et al. Genome-wide association study to identify genetic loci associated with gastrointestinal nematode resistance in Katahdin sheep[J]. Anim Genet, 2020, 51(2): 330-335. DOI:10.1111/age.12895
[109]
RASCHIA M A, DONZELLI M V, MEDUS P D, et al. Single nucleotide polymorphisms from candidate genes associated with nematode resistance and resilience in Corriedale and Pampinta sheep in Argentina[J]. Gene, 2021, 770: 145345. DOI:10.1016/j.gene.2020.145345
[110]
LI G D, LV D Y, YAO Y J, et al. Overexpression of ASMT likely enhances the resistance of transgenic sheep to brucellosis by influencing immune-related signaling pathways and gut microbiota[J]. FASEB J, 2021, 35(9): e21783.
[111]
LI X L, WU Q M, ZHANG X X, et al. Whole-genome resequencing to study brucellosis susceptibility in sheep[J]. Front Genet, 2021, 12: 653927. DOI:10.3389/fgene.2021.653927
[112]
HARLAND D P, PLOWMAN J E. Development of hair fibres[M]//PLOWMAN J E, HARLAND D P, DEB-CHOUDHURY S. The Hair Fibre: Proteins, Structure and Development. Singapore: Springer, 2018: 109-154.
[113]
MATSUNAGA R, ABE R, ISHII D, et al. Bidirectional binding property of high glycine-tyrosine keratin-associated protein contributes to the mechanical strength and shape of hair[J]. J Struct Biol, 2013, 183(3): 484-494. DOI:10.1016/j.jsb.2013.06.005
[114]
HE D Q, CHEN L Y, LUO F, et al. Differentially phosphorylated proteins in the crimped and straight wool of Chinese Tan sheep[J]. J Proteomics, 2021, 235: 104115. DOI:10.1016/j.jprot.2021.104115
[115]
HÉBERT J M, ROSENQUIST T, GÖTZ J, et al. FGF5 as a regulator of the hairgrowth cycle: evidence from targeted and spontaneous mutations[J]. Cell, 1994, 78(6): 1017-1025. DOI:10.1016/0092-8674(94)90276-3
[116]
LI W R, LIU C X, ZHANG X M, et al. CRISPR/Cas9-mediated loss of FGF5 function increases wool staple length in sheep[J]. FEBS J, 2017, 284(17): 2764-2773. DOI:10.1111/febs.14144
[117]
ZHANG R, LI Y, JIA K, et al. Crosstalk between androgen and Wnt/β-catenin leads to changes of wool density in FGF5-knockout sheep[J]. Cell Death Dis, 2020, 11(5): 407. DOI:10.1038/s41419-020-2622-x
[118]
胡馨予, 姚逸安, 胡情情, 等. 羊毛角蛋白基因家族及其启动子调控作用研究进展[J]. 中国细胞生物学学报, 2021, 43(8): 1705-1713.
HU X Y, YAO Y A, HU Q Q, et al. Research progress of wool keratin gene family and its promoter regulation role[J]. Chinese Journal of Cell Biology, 2021, 43(8): 1705-1713. (in Chinese)
[119]
ULLAH F, JAMAL S M, ZHOU H T, et al. Variation in the KRTAP6-3 gene and its association with wool characteristics in Pakistani sheep breeds and breed-crosses[J]. Trop Anim Health Prod, 2020, 52(6): 3035-3043. DOI:10.1007/s11250-020-02322-6
[120]
RAMIREZ J M, FOLKOW L P, BLIX A S. Hypoxia tolerance in mammals and birds: from the wilderness to the clinic[J]. Annu Rev Physiol, 2007, 69: 113-143. DOI:10.1146/annurev.physiol.69.031905.163111
[121]
BIGHAM A W, LEE F S. Human high-altitude adaptation: forward genetics meets the HIF pathway[J]. Genes Dev, 2014, 28(20): 2189-2204. DOI:10.1101/gad.250167.114
[122]
LIU X X, ZHANG Y L, LI Y F, et al. EPAS1 gain-of-function mutation contributes to high-altitude adaptation in Tibetan horses[J]. Mol Biol Evol, 2019, 36(11): 2591-2603. DOI:10.1093/molbev/msz158
[123]
MA Y F, HAN X M, HUANG C P, et al. Population genomics analysis revealed origin and high-altitude adaptation of Tibetan pigs[J]. Sci Rep, 2019, 9(1): 11463. DOI:10.1038/s41598-019-47711-6
[124]
XU X H, HUANG X W, QUN L, et al. Two functional loci in the promoter of EPAS1 gene involved in high-altitude adaptation of Tibetans[J]. Sci Rep, 2014, 4: 7465. DOI:10.1038/srep07465
[125]
LORENZO F R, HUFF C, MYLLYMÄKI M, et al. A genetic mechanism for Tibetan high-altitude adaptation[J]. Nat Genet, 2014, 46(9): 951-956. DOI:10.1038/ng.3067
[126]
SONG S, YAO N, YANG M, et al. Exome sequencing reveals genetic differentiation due to high-altitude adaptation in the Tibetan cashmere goat (Caprahircus)[J]. BMC Genomics, 2016, 17: 122. DOI:10.1186/s12864-016-2449-0
[127]
SASAZAKI S, TOMITA K, NOMURA Y, et al. FGF5 and EPAS1 gene polymorphisms are associated with high-altitude adaptation in Nepalese goat breeds[J]. Anim Sci J, 2021, 92(1): e13640.
[128]
WEI C H, WANG H H, LIU G, et al. Genome-wide analysis reveals adaptation to high altitudes in Tibetan sheep[J]. Sci Rep, 2016, 6: 26770. DOI:10.1038/srep26770
[129]
LIU J B, YUAN C, GUO T T, et al. Genetic signatures of high-altitude adaptation and geographic distribution in Tibetan sheep[J]. Sci Rep, 2020, 10(1): 18332. DOI:10.1038/s41598-020-75428-4
[130]
GORKHALI N A, DONG K Z, YANG M, et al. Genomic analysis identified a potential novel molecular mechanism for high-altitude adaptation in sheep at the Himalayas[J]. Sci Rep, 2016, 6: 29963. DOI:10.1038/srep29963
[131]
WU P P, ZHANG B, HAN X Y, et al. HucMSC exosome-delivered 14-3-3ζ alleviates ultraviolet radiation-induced photodamage via SIRT1 pathway modulation[J]. Aging, 2021, 13(8): 11542-11563. DOI:10.18632/aging.202851
[132]
ZHANG Q, GOU W Y, WANG X T, et al. Genome resequencing identifies unique adaptations of Tibetan chickens to hypoxia and high-Dose ultraviolet radiation in high-altitude environments[J]. Genome Biol Evol, 2016, 8(3): 765-776. DOI:10.1093/gbe/evw032
[133]
HU X J, YANG J, XIE X L, et al. The genome landscape of Tibetan sheep reveals adaptive introgression from argali and the history of early human settlements on the Qinghai-Tibetan plateau[J]. Mol Biol Evol, 2019, 36(2): 283-303. DOI:10.1093/molbev/msy208
[134]
MWACHARO J M, KIM E S, ELBELTAGY A R, et al. Genomic footprints of dryland stress adaptation in Egyptian fat-tail sheep and their divergence from East African and Western Asia cohorts[J]. Sci Rep, 2017, 7(1): 17647. DOI:10.1038/s41598-017-17775-3
[135]
EDEA Z, DADI H, DESSIE T, et al. Genomic signatures of high-altitude adaptation in Ethiopian sheep populations[J]. Genes Genom, 2019, 41(8): 973-981. DOI:10.1007/s13258-019-00820-y
[136]
NOWACK J, GIROUD S, ARNOLD W, et al. Muscle non-shivering thermogenesis and its role in the evolution of endothermy[J]. Front Physiol, 2017, 8: 889. DOI:10.3389/fphys.2017.00889
[137]
GAUDRY M J, JASTROCH M. Comparative functional analyses of UCP1 to unravel evolution, ecophysiology and mechanisms of mammalian thermogenesis[J]. Comp Biochem Physiol PartB: Biochem Mol Biol, 2021, 255: 110613. DOI:10.1016/j.cbpb.2021.110613
[138]
SLEE J, STOTT A W. Genetic selection for cold resistance in Scottish Blackface lambs[J]. Anim Sci, 1986, 43(3): 397-404. DOI:10.1017/S0003356100002622
[139]
WANG Z H, YU X F, CHEN Y. Recruitment of thermogenic fat: trigger of fat burning[J]. Front Endocrinol, 2021, 12: 696505. DOI:10.3389/fendo.2021.696505
[140]
POHL E E, RUPPRECHT A, MACHER G, et al. Important trends in UCP3 investigation[J]. Front Physiol, 2019, 10: 470. DOI:10.3389/fphys.2019.00470
[141]
CARON A, LABBÉ S M, CARTER S, et al. Loss of UCP2 impairs cold-induced non-shivering thermogenesis by promoting a shift toward glucose utilization in brown adipose tissue[J]. Biochimie, 2017, 134: 118-126. DOI:10.1016/j.biochi.2017.01.006
[142]
HENRY B A, ANDREWS Z B, RAO A, et al. Central leptin activates mitochondrial function and increases heat production in skeletal muscle[J]. Endocrinology, 2011, 152(7): 2609-2618. DOI:10.1210/en.2011-0143
[143]
AHMAD S F, MEHROTRA A, CHARLES S, et al. Analysis of selection signatures reveals important insights into the adaptability of high-altitude Indian sheep breed Changthangi[J]. Gene, 2021, 799: 145809. DOI:10.1016/j.gene.2021.145809
[144]
BAUMGARD L H, RHOADS R P Jr. Effects of heat stress on postabsorptive metabolism and energetics[J]. Annu Rev Anim Biosci, 2013, 1: 311-337. DOI:10.1146/annurev-animal-031412-103644
[145]
ABDELNOUR S A, EL-HACK M E A, KHAFAGA A F, et al. Stress biomarkers and proteomics alteration to thermal stress in ruminants: a review[J]. J Therm Biol, 2019, 79: 120-134. DOI:10.1016/j.jtherbio.2018.12.013
[146]
QUINN C M, DURAN R M, AUDET G N, et al. Cardiovascular and thermoregulatory biomarkers of heat stroke severity in a conscious rat model[J]. J Appl Physiol, 2014, 117(9): 971-978. DOI:10.1152/japplphysiol.00365.2014
[147]
ROMO-BARRON C B, DIAZ D, PORTILLO-LOERA J J, et al. Impact of heat stress on the reproductive performance and physiology of ewes: a systematic review and meta-analyses[J]. Int J Biometeorol, 2019, 63(7): 949-962. DOI:10.1007/s00484-019-01707-z
[148]
PIRKKALA L, NYKÄNEN P, SISTONEN L. Roles of the heat shock transcription factors in regulation of the heat shock response and beyond[J]. FASEB J, 2001, 15(7): 1118-1131. DOI:10.1096/fj00-0294rev
[149]
SALCES-ORTIZ J, GONZÁLEZ C, BOLADO-CARRANCIO A, et al. Ovine HSP90AA1 gene promoter: functional study and epigenetic modifications[J]. Cell Stress Chaperon, 2015, 20(6): 1001-1012. DOI:10.1007/s12192-015-0629-5
[150]
AL-THUWAINI T M, AL-SHUHAIB M B S, HUSSEIN Z M. A novel T177P missense variant in the HSPA8 gene associated with the low tolerance of Awassi sheep to heat stress[J]. Trop Anim Health Prod, 2020, 52(5): 2405-2416. DOI:10.1007/s11250-020-02267-w
[151]
WANG W M, ZHANG X X, ZHOU X, et al. Deep genome resequencing reveals artificial and natural selection for visual deterioration, plateau adaptability and high prolificacy in Chinese domestic sheep[J]. Front Genet, 2019, 10: 300. DOI:10.3389/fgene.2019.00300
[152]
LAMPIS A, CAROTENUTO P, VLACHOGIANNIS G, et al. MIR21 drives resistance to Heat Shock Protein 90 inhibition in cholangiocarcinoma[J]. Gastroenterology, 2018, 154(4): 1066-1079.e5. DOI:10.1053/j.gastro.2017.10.043

(编辑   郭云雁)