畜牧兽医学报  2021, Vol. 52 Issue (11): 3126-3136. DOI: 10.11843/j.issn.0366-6964.2021.011.014    PDF    
犏牛和牦牛ESCO2基因的序列特征及其在睾丸中的表达对比分析
闵星宇1, 杨丽雪2, 杨满珍1, 胡宇磊3, 于海玲1, 杨璐瑜3, 李键1,2, 熊显荣1,2     
1. 西南民族大学畜牧兽医学院, 成都 610041;
2. 青藏高原动物遗传育种资源保护与利用国家教育部重点实验室, 成都 610041;
3. 动物科学国家民委重点实验室, 成都 610041
摘要:旨在克隆犏牛和牦牛姐妹染色单体内聚建立蛋白2(establishment of sister chromatid cohesion N-acetyltransferase 2,ESCO2)基因,并分析其在不同发育阶段睾丸中的表达与定位,为进一步解析ESCO2在减数分裂过程中的作用机制提供理论依据。本研究以健康雄性犏牛及牦牛为试验动物,根据年龄分为胎牛组(5~6月龄)、幼年组(1~2岁)和成年组(3~4岁),每组各3头。通过RT-PCR技术克隆犏牛和牦牛ESCO2基因并进行生物信息学分析,采用实时荧光定量PCR(quantitative real-time PCR,qRT-PCR)检测ESCO2基因在犏牛不同组织中的表达谱,比较分析ESCO2在犏牛和牦牛不同时期睾丸中的表达规律,利用免疫组织化学(immunohistochemistry,IHC)技术检测ESCO2蛋白的细胞定位和表达差异。结果显示,犏牛ESCO2基因(GenBank登录号:MW198470)CDS区为1 833 bp,编码610个氨基酸,与牦牛相比,犏牛ESCO2序列第301~319位多19个氨基酸,另有3个氨基酸突变;犏牛ESCO2蛋白序列与黄牛的同源性高于其他哺乳动物;ESCO2可能与SMC3、SMC1A、PDS5A、PDS5B、STAG2等蛋白相互作用,互作蛋白功能与姐妹染色单体凝聚、减数分裂细胞周期、DNA修复、细胞分裂和染色体重构等生物学过程相关。ESCO2在犏牛各组织中均有表达,但在睾丸中的相对表达水平显著高于其它组织(P < 0.05);在犏牛睾丸中的表达随年龄增长呈上升趋势,幼年和成年时期犏牛睾丸中ESCO2的表达显著低于同时期牦牛(P < 0.05);IHC染色结果发现,雄性犏牛减数分裂阻滞于初级精母细胞,ESCO2蛋白在犏牛初级精母细胞中无表达并与牦牛存在差异。本研究结果表明,犏牛与牦牛的ESCO2基因、蛋白序列差异较大,且在睾丸发育过程中的表达模式差异显著,这可能是引起雄性犏牛减数分裂阻滞及不育的原因之一,其具体作用机制有待进一步研究。
关键词犏牛    ESCO2    睾丸    减数分裂    基因表达    
Sequence Characteristics and Comparative Expression Analysis in Testis of ESCO2 Gene in Cattle-yak and Yak
MIN Xingyu1, YANG Lixue2, YANG Manzhen1, HU Yulei3, YU Hailing1, YANG Luyu3, LI Jian1,2, XIONG Xianrong1,2     
1. College of Animal Science and Veterinary Medicine, Southwest Minzu University, Chengdu 610041, China;
2. Key Laboratory of Qinghai-Tibetan Plateau Animal Genetic Resource Reservation and Utilization of Ministry of Education, Chengdu 610041, China;
3. Key Laboratory of Animal Science of State Ethnic Affairs Commission, Chengdu 610041, China
Abstract: The aim of this study was to clone the establishment of sister chromatid cohesion N-acetyltransferase 2(ESCO2) gene, analyze its expression and localization in the testis at different developmental stages in cattle-yak and yak, and provide a theoretical basis for further study the role mechanism of ESCO2 in the male meiosis process. Healthy male cattle-yaks and yaks were used as experimental animals in this study, they were divided into fetal cattle group (5-6 months old), juvenile group (1-2 years old) and adult group (3-4 years old), with 3 animals in each group. The CDS region of ESCO2 gene of cattle-yak and yak were cloned by RT-PCR and its sequence characteristics were analyzed by bioinformatics softwares. qRT-PCR was used to detect the expression patterns of ESCO2 gene in different tissues of cattle-yak. The expression of ESCO2 mRNA in testes at different developmental stages of cattle-yak and yak were detected by qRT-PCR, the differences of cell location and expression of ESCO2 protein were detected by immunohistochemistry (IHC) staining. The results showed that the CDS region of ESCO2 genes in cattle-yak was 1 833 bp(MW198470), it encoded 610 amino acids. Compared with yak, the ESCO2 protein sequence of cattle-yak had more 19 amino acids at position 301-319 aa and there were 3 amino acid mutations. The cattle-yak and cattle had higher homology than other mammals based on the protein sequence of ESCO2. Interaction network showed that ESCO2 mainly interacted with SMC3, SMC1A, PDS5A, PDS5B, STAG2 and these proteins were related to sister chromatid condensation, meiotic cell cycle, DNA repair, cell division and chromosome remodeling. ESCO2 was expressed in all tissues of the cattle-yak, but its relative expression level in the testis was significantly higher than those in other tissues (P < 0.05). In addition, the spatial expression of ESCO2 in testes of cattle-yak showed an upward trend with age, and the relative expression abundance of ESCO2 in cattle-yak at juvenile and adult stages was significantly lower than that of yak(P < 0.05). IHC staining showed that meiosis of male cattle-yak was blocked in primary spermatocytes and ESCO2 protein was differentially expressed in primary spermatocytes between cattle-yak and yak. The research suggested that there was obvious difference in the sequence and expression pattern of ESCO2 between cattle-yak and yak, which may be one of reasons caused male sterility in cattle-yak, but the specific mechanism need further study.
Key words: cattle-yak    ESCO2    testis    meiosis    gene expression    

犏牛是牦牛(Bos grunniens)与普通牛(Bos taurus)通过远缘杂交得到的F1代,按生产功能可分为乳用、肉用和役用犏牛等。牦牛生产性能低下,普通牛又无法适应青藏高原恶劣的高寒低氧环境,而犏牛具有显著的杂种优势,其体格大、生长快,产奶和产肉性能与牦牛相比都具有显著的优势,且具有良好的高原适应性[1-3],深受藏区牧民喜爱。但犏牛雄性不育的缺陷,使得无法通过横交固定其杂种优势,这严重阻碍了杂种优势的利用[4-5]。因此,解析雄性犏牛不育对提高青藏高原畜牧生产水平、改善牧民生活品质和发展青藏高原特色畜牧产业都具备巨大的潜在价值。

ESCO2是一种黏连蛋白乙酰基转移酶,其主要功能是在S期介导黏连蛋白复合体中染色体结构维持蛋白2(SMC3)亚基的乙酰化,而黏连蛋白复合体在Sororin的作用下将姐妹染色单体束缚在一起,可确保它们在细胞分裂过程中准确对齐和分离[6-9],是细胞分裂过程中黏连蛋白复合体乙酰化和着丝粒结合的必要条件[10]。人ESCO2基因的突变会导致一种常染色体隐形遗传疾病——罗伯茨综合征(Roberts syndrome, RBS),其患者表现为出生前和出生后发育迟缓,颅面、小头、肢体畸形和精神缺陷[11-12]。Hogarth等[13]研究发现,ESCO2是哺乳动物有丝分裂到减数分裂过渡的潜在调节因子,在特异性敲除不孕小鼠卵巢[14]和低生殖力鲶鱼性腺[15]中表现为下调。小鼠卵母细胞减数分裂成熟过程中,ESCO2通过与组蛋白H4结合并介导H4K16乙酰化来调节纺锤体组装检查点(SAC),耗竭ESCO2会导致极体的提前排出,并造成非整倍体卵的产生[16-18]。ESCO2能够修复双链染色体的断裂,在雄性减数分裂中发挥重要的作用,其异常表达可引起精子发生障碍和雄性不育[16]。Mcnicoll等[19]的研究揭示,ESCO2是雄性配子发生过程中必不可少的调节因子,特异性敲除ESCO2导致雄性小鼠不育。然而,关于ESCO2基因的研究主要集中在小鼠和人,该基因近期在犏牛和牦牛上的相关研究尚未见报道。

本研究以犏牛为研究对象,利用RT-PCR技术克隆获得犏牛ESCO2编码区完整序列,并对该序列进行生物信息学分析,同时采用qRT-PCR和IHC染色检测其在犏牛各组织及不同发育阶段睾丸组织中的表达,以期为揭示ESCO2在减数分裂中的作用及调控精子发生机制提供参考,为进一步从分子水平解析犏牛雄性不育机制奠定基础。

1 材料与方法 1.1 试验动物与样品采集

动物样品采自成都市周边牦牛屠宰场,健康犏牛和牦牛宰杀后立即使用无菌剪采集心、肝、脾、肺、肾、大肠、小肠、胃、肌肉、脂肪和附睾组织,随后采集睾丸组织,根据年龄分为胎牛时期(5~6月龄)、幼年时期(1~2岁)和成年时期(3~4岁)。以上所有样本每组各取3头为生物学重复,采集的样本使用高压灭菌处理的生理盐水冲洗剪成0.5 cm×0.5 cm× 0.5 cm的组织块,部分置于液氮罐中带回实验室,于-80 ℃保存备用;部分使用4%多聚甲醛固定,用于免疫组织化学染色。

1.2 组织总RNA的提取与反转录

犏牛和牦牛的组织经液氮快速研磨后,使用Trizol Reagent(Invitrogen,美国)进行总RNA的提取,使用紫外分光光度计(Biospec-nano,日本)检测RNA的浓度和纯度,选取OD260 nm/OD280 nm介于1.8~2.0的RNA,经1%琼脂糖凝胶电泳检测其完整性,备用。根据RevertAid First Strand cDNA Synthesis Kit反转录试剂盒(Thermo Scientific,美国)使用说明合成第一链cDNA,置于-20 ℃保存备用。

1.3 引物的设计与基因克隆

从NCBI数据库中选取野牦牛(Bos mutus)ESCO2基因预测序列(登录号:XM_005900074)和β-actin基因序列(登录号:DQ838049.1),使用Primer Premier 5软件设计引物(表 1),并送至南京金斯瑞生物科技有限公司进行合成。以睾丸组织cDNA为模板,利用RT-PCR技术在PCR仪(Applied Biosystems,美国)扩增犏牛ESCO2基因CDS区序列,PCR扩增体系为25 μL:2×Phanta Max Master Mix(Dye Plus)(南京诺唯赞生物科技股份有限公司)12.5 μL,模板1 μL,上、下游引物各1 μL,ddH2O 9.5 μL。PCR扩增程序:95 ℃预变性3 min;95 ℃变性15 s,59 ℃退火15 s,72 ℃延伸60 s,30个循环;72 ℃彻底延伸5 min,4 ℃保存。使用浓度为1.5%琼脂糖凝胶电泳分离后,按照胶回收试剂盒(康宁生命科学有限公司)说明对目的扩增产物进行纯化回收,交由生工生物工程(上海)股份有限公司成都测序实验室测序。

表 1 PCR引物序列及产物长度 Table 1 PCR primer sequences and products length
1.4 ESCO2生物信息学分析

利用ORFfinder(https://www.ncbi.nlm.nih.gov/orffinder/)查找犏牛ESCO2克隆序列开放阅读框并翻译成氨基酸序列,采用EMBOSS Needle(https://www.ebi.ac.uk/Tools/psa/emboss_needle/)在线软件比对犏牛和牦牛ESCO2基因CDS区核苷酸序列。采用T-COFFEE(http://tcoffee.crg.cat/apps/tcoffee/index.html)在线软件并结合Jalview软件多重比较犏牛和牦牛与其他物种的氨基酸序列,通过MEME(https://meme-suite.org/meme/)在线工具分析查找具有高度保守性的基序,并用MEGA7软件构建ESCO2蛋白的NJ系统发生树。利用ExPASy(https://web.expasy.org/protparam/)、ProtParam(https://web.expasy.org/protparam/)、InterPro(http://www.ebi.ac.uk/interpro/)和SignalP 4.0(http://www.cbs.dtu.dk/services/SignalP/)在线软件进行ESCO2蛋白基本理化性质、功能结构域和信号肽分析。通过SOPMA(https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html)、SWISS-MODEL(https://swissmodel.expasy.org/)和STRING 11.0(https://string-db.org/)在线工具预测ESCO2蛋白二级结构、三级结构和互作蛋白。

1.5 犏牛ESCO2基因的组织表达谱

β-actin为内参基因,通过实时荧光定量仪(Bio-Rad,美国)检测ESCO2基因在犏牛各组织(心、肝、脾、肺、肾、大肠、小肠、胃、肌肉、脂肪、附睾和睾丸组织)中的表达量并进行表达谱分析,及其在不同时期(胎牛时期、幼年时期和成年时期)犏牛和牦牛睾丸组织中的表达量。20 μL反应体系如下:2×NovaStart SYBR qPCR Super Mix plus 10 μL,cDNA 1 μL, 上、下游引物各0.5 μL,ddH2O 8 μL。PCR扩增条件如下:95 ℃预变性1 min;95 ℃变性20 s,60 ℃退火1 min,39个循环;熔解曲线为65~95 ℃每10 s增加0.5 ℃。

1.6 免疫组织化学染色检测ESCO2在犏牛及牦牛睾丸中的定位

将犏牛和牦牛胎牛时期、幼年时期和成年时期的睾丸组织从4%多聚甲醛组织固定液中取出,制作不同发育时期睾丸石蜡切片。依次将切片放入二甲苯中脱蜡,使用梯度酒精洗涤后,置于EDTA抗原修复缓冲液(pH 9.0)的修复盒中,在微波炉中加热进行抗原修复;切片自然冷却后在PBS中晃动洗涤后放入3%双氧水溶液中阻断内源性过氧化物酶;3%BSA均匀覆盖组织室温封闭30 min后,滴加一抗(多克隆兔抗ESCO2,1∶100稀释,博奥森)4 ℃孵育过夜;使用PBS洗涤玻片3次后,滴加相应种属的二抗(辣根过氧化物酶标记)覆盖组织,室温孵育50 min;PBS洗涤后滴加新鲜配置的DAB显色液显色,自来水冲洗切片终止显色;苏木精复染3 min左右,复染细胞核;依次将切片放入梯度酒精和二甲苯脱水透明后,使用中性树胶封片。使用激光共聚焦显微镜(Zeiss,德国)观察并拍照。

1.7 数据统计分析

采用2-ΔΔCt法对qRT-PCR结果进行计算,所有试验数据使用“平均值±标准误(Mean±SEM)”表示。使用SPSS 19.0软件进行单因素方差分析,P < 0.05表示差异显著,P < 0.01表示差异极显著。每组试验至少重复3次。

2 结果 2.1 犏牛和牦牛ESCO2基因CDS区序列特征分析

本研究以犏牛和牦牛睾丸组织cDNA为模板,F1、R1为上下游引物,对ESCO2基因的CDS区进行扩增,经1.5%琼脂糖凝胶电泳检测,目标条带与预期结果一致(图 1)。经测序获得犏牛2 742 bp的核苷酸序列,其犏牛ESCO2基因CDS长度为1 833 bp,可编码610个氨基酸;牦牛2 690 bp的核苷酸序列,其牦牛ESCO2基因CDS长度为1 776 bp,可编码591个氨基酸。犏牛与牦牛ESCO2基因CDS区核苷酸序列的比对结果显示共有7处变化,同源性为96.6%,其中在第4处犏牛比牦牛多57个碱基,导致氨基酸序列N端第301~319位多编码19个氨基酸,第2处C→T突变,编码的亮氨酸变为苯丙氨酸,第3处G→A突变,编码的甘氨酸变为天冬氨酸,第6处C→T突变,编码的苏氨酸变为甲硫氨酸,另外3处均为同义突变。犏牛ESCO2蛋白的分子量、等电点和分子式分别为69 049.02、9.36和C3038H4890N852O927S27;牦牛ESCO2蛋白的分子量、等电点和分子式分别为66 781.62、9.48和C2947H4755N827O890S24。犏牛和牦牛ESCO2蛋白均不含信号肽且无跨膜结构域。

M. DNA相对分子质量标准; 1. 犏牛ESCO2基因扩增产物; 2. 牦牛ESCO2基因扩增产物; 3. 空白对照; 4. 犏牛β-actin扩增产物 M. DL5000 marker; 1.PCR product of cattle-yak ESCO2 gene; 2. PCR product of yak ESCO2 gene; 3. Black control; 4. PCR product of cattle-yak β-actin gene 图 1 犏牛和牦牛ESCO2基因的PCR扩增结果 Fig. 1 PCR amplification results of cattle-yak and yak ESCO2 gene
2.2 犏牛ESCO2蛋白的同源性分析

多重比较犏牛和其它14种哺乳动物ESCO2蛋白的氨基酸序列,结果发现,犏牛与黄牛的氨基酸序列同源性最高(100%),依次为水牛(98.69%)、牦牛(96.39%)、山羊(96.23%)、抹香鲸(88.07%)、猪(84.80%)、单峰驼(84.31%)、家马(84.20%)、灰熊(81.43%)、家猫(80.46%)、家犬(76.80%)、家兔(76.62)、智人(76.26%)和小鼠(67.59%)。通过MEME软件在线分析获得ESCO2蛋白最保守的3个基序,分别位于犏牛ESCO2氨基酸序列N端第383~432、434~483和526~565位点,对比发现不同物种ESCO2蛋白对应基序氨基酸序列具有高度保守性。使用MEGA7软件分析15种哺乳动物的氨基酸序列并构建NJ系统发生树,发现犏牛和黄牛聚为一支,在进化上与水牛、牦牛和山羊亲缘关系较近,与其他物种较远,这与同源性结果一致(图 2)。

图 2 ESCO2蛋白氨基酸序列系统进化树 Fig. 2 Phylogenetic tree for ESCO2 amino acid sequences
2.3 ESCO2蛋白结构功能预测

采用SOPMA预测ESCO2蛋白的二级结构,犏牛340个氨基酸(55.74%)形成无规则卷曲、125个氨基酸(20.49%)形成α螺旋、119个氨基酸(19.51%)形成延伸链和26个氨基酸(4.26%)形成β转角(图 3A);通过SWISS-MODEL预测到犏牛ESCO2蛋白分子三级结构(图 3B)。采用STRING数据库分析犏牛ESCO2潜在相互作用蛋白显示,ESCO2可能与SMC3、SMC1A、PDS5A、PDS5B、STAG2等蛋白存在相互作用(图 3C)。GO分析显示,ESCO2蛋白互作网络中蛋白质能够参与姐妹染色单体凝聚、减数分裂细胞周期、DNA修复、细胞分裂和染色体重构等生物学过程。

A. 犏牛ESCO2蛋白二级结构;B. 犏牛ESCO2蛋白三级结构;C. 犏牛ESCO2的蛋白互作网络 A. Secondary structure of cattle-yak ESCO2 protein; B. Tertiary structure of cattle-yak ESCO2 protein; C. Interaction network for cattle-yak ESCO2 protein 图 3 犏牛ESCO2蛋白的空间结构及互作蛋白 Fig. 3 Spatial structure and interaction proteins of cattle-yak ESCO2 protein
2.4 犏牛ESCO2基因的组织表达谱

β-actin作为内参基因,利用qRT-PCR技术检测ESCO2基因在犏牛肌肉、脂肪、睾丸、附睾、心、肝、脾、肺、肾、大肠、小肠、胃组织中的相对表达水平(图 4)。结果显示,犏牛ESCO2 mRNA在各个组织中均有表达,其中睾丸组织表达最高,显著高于其他组织(P < 0.05)。

不同的字母表示差异显著(P < 0.05)。下同 Different letters represent the significant differences(P < 0.05). The same as below 图 4 ESCO2 mRNA在犏牛各组织的表达谱 Fig. 4 Tissue expression profile of ESCO2 mRNA in cattle-yak
2.5 不同发育时期睾丸中ESCO2基因的表达

β-actin作为内参基因,检测犏牛和牦牛胎牛时期、幼年时期和成年时期睾丸中ESCO2基因的相对表达水平(图 5)。结果显示,ESCO2 mRNA的表达在犏牛和牦牛不同发育时期睾丸中均呈上升趋势;5~6月龄犏胎牛睾丸中的表达量低于牦胎牛,差异不显著(P>0.05),但1~2岁和3~4岁犏牛睾丸中ESCO2 mRNA的表达量显著低于同时期牦牛(P < 0.05)。

图 5 ESCO2 mRNA在不同发育阶段牦牛和犏牛睾丸中的表达规律 Fig. 5 Expression patterns of ESCO2 mRNA in testis at different developmental stages of cattle-yak and yak
2.6 犏牛及牦牛ESCO2蛋白的细胞定位

免疫组织化学技术检测ESCO2在犏牛和牦牛不同年龄阶段睾丸组织中的定位和表达(图 6),其中细胞着色为棕黄色或者黄色判定为阳性信号,颜色越深说明蛋白表达水平越高。结果显示,ESCO2蛋白在犏牛各年龄段睾丸组织中均有表达。其中胎牛睾丸组织中,ESCO2蛋白主要定位于支持细胞(SC),在精原干细胞(SSC)几乎不表达,与同时期牦牛一致;幼年和成年时期犏牛睾丸组织中,生精小管发育不良,管径小且管壁皱缩,无次级精母细胞(SS)和精子细胞,在支持细胞、精原细胞(SP)和间质细胞(IC)检测到ESCO2蛋白阳性信号,与同时期牦牛相比,犏牛初级精母细胞(PS)数量少且细胞核皱缩,无ESCO2蛋白阳性信号。

SP. 精原细胞; PS. 初级精母细胞; SS. 次级精母细胞; SC. 支持细胞; IC. 间质细胞; MC. 肌样细胞; RS. 圆形精子细胞; ES. 长形精子细胞 SP. Spermatogonium; PS. Primary spermatocyte; SS. Secondary spermatocyte; SC. Sertoli cell; IC. Leydig cell; MC. Myoid cell; RS. Round spermatid; ES. Elongating spermatid 图 6 不同阶段雄性犏牛和牦牛生殖细胞发育及ESCO2蛋白定位 Fig. 6 The growth of male cattle-yak and yak germ cells and localization of ESCO2 protein at different stages
3 讨论

在哺乳动物精子发生过程中,蛋白质的乙酰化能调控减数分裂,其异常会导致雄性生殖力下降[20]。其中ESCO2是一种非组蛋白乙酰基转移酶,能调节减数分裂过程中姐妹染色单体的分离,是哺乳动物配子发生所必须的[21-23]。因此,探索ESCO2在睾丸中的表达规律对于解析犏牛雄性不育机制具有重要的研究价值。

本研究克隆获得犏牛和牦牛ESCO2完整CDS区序列,并对其编码蛋白的结构和功能进行了预测和分析。对犏牛和牦牛ESCO2编码区序列进行对比,发现犏牛增加了57个碱基并有3处错义突变。基因的错义突变能够改变多肽链的氨基酸种类从而使蛋白质结构和功能发生改变并造成疾病,如白化病、镰刀型贫血病和早衰症等[24-26]ESCO2基因的突变在人当中能够造成罗伯茨综合征[27],而犏牛与牦牛ESCO2蛋白的差异是否改变蛋白功能有待进一步研究。鉴于ESCO2蛋白主要通过与相应蛋白结合发生作用,本研究进一步分析犏牛ESCO2的互作蛋白,发现与SMC3、SMC1A、PDS5A等10个蛋白存在相互作用。其中,SMC1A、SMC1B、SMC3、RAD21和STAG2是黏连蛋白复合体的组成亚基,有研究表明,ESCO2能够介导SMC3的去乙酰化[6],黏连蛋白复合体具有细胞分裂过程中凝聚姐妹染色单体的作用,姐妹染色单体内聚力的过早丧失会导致染色体的错误分离[28-29]。Al-Jomah等[30]的研究表明,PDS5A和PDS5B是一种富含HEAT重复序列的蛋白质,可与黏附素结合并介导姐妹染色单体黏附,PDS5与WAPL直接相互作用,在有丝分裂过程中去除黏附素,它们的丢失会导致DNA损伤。精母细胞和卵母细胞中NIPBL与黏连蛋白复合体有相同的定位,人NIPBL基因的突变能够导致德朗热综合征[31-32]。本研究发现,犏牛与牦牛ESCO2基因和编码蛋白存在差异,并对其结构和互作蛋白进行了预测,这些与犏牛ESCO2互作的潜在蛋白均在细胞有丝分裂和减数分裂中扮演不同的角色,这对于下一步解析ESCO2在雄性犏牛不育中的分子机制提供了新的见解和思路。

组织表达谱分析显示,ESCO2 mRNA在各组织中均有表达,其中睾丸、脾、大肠、小肠和胃中表达较高。该结果与现有研究结果一致,已报道ESCO2在人睾丸、骨髓、各级消化道组织中表达较高,在褐家鼠睾丸、脾和胸腺中具有较高的表达量[33-34]。Chen等[35]研究揭示,ESCO2能够抑制细胞增殖并诱导胃癌细胞凋亡,在大肠癌中亦可抑制肿瘤转移[36],而ESCO2在犏牛和人各消化道组织中高表达,推测能够维持消化道细胞的正常功能,其具体作用机制有待进一步研究。ESCO2在人、褐家鼠、犏牛等物种睾丸中高表达,因此推测其在雄性生殖中发挥了重要作用。

本试验发现,ESCO2在犏牛睾丸中的表达量随着生殖发育呈上升趋势,在各年龄阶段犏牛睾丸中ESCO2的表达量均低于牦牛。有研究报道,在小鼠卵巢和睾丸中特异性敲除ESCO2能够造成生殖力下降[14, 19],在低生殖力动物的性腺中ESCO2也表现出低表达[15],因此推测,雄性犏牛ESCO2的低表达也能造成生殖力的下降。为进一步探讨ESCO2在睾丸中发挥的功能,利用IHC染色对比分析ESCO2在犏牛和牦牛睾丸中的细胞定位和表达水平,观察到雄性犏牛减数分裂阻滞于初级精母细胞,生精小管内只能观察到支持细胞、精原细胞和少量初级精母细胞,这与前人报道一致[2, 37-39]。Evans等[16]的研究表明,ESCO2蛋白定位于减数分裂Ⅰ期精母细胞染色体上,本试验在牦牛初级精母细胞检测到阳性信号,但犏牛初级精母细胞上未检测到表达和定位,可能ESCO2在犏牛初级精母细胞上的缺失是造成减数分裂阻滞的原因之一。此外,在精原细胞、次级精母细胞、支持细胞和间质细胞中也检测到ESCO2蛋白定位,推测ESCO2在雄性生殖中还发挥了多种生物学功能。山羊精原干细胞中过表达ESCO2基因能够促进细胞进入减数分裂[40],而犏牛精原细胞上ESCO2蛋白阳性信号也弱于牦牛,提示ESCO2蛋白在犏牛精原细胞上的低表达可能导致初级精母细胞减少。综上推测,ESCO2是精子发生中重要的功能基因,而犏牛睾丸中的低表达可能是导致减数分裂阻滞的原因之一。

4 结论

本试验成功获得了犏牛和牦牛的ESCO2基因序列,发现ESCO2基因在睾丸组织中表达最高,其在成年期睾丸中的表达量高于幼年期和胎牛期,犏牛睾丸组织中ESCO2 mRNA表达水平显著低于同时期牦牛,该蛋白在犏牛和牦牛睾丸组织中存在差异定位,表明ESCO2基因与雄性生殖相关。本研究为进一步探讨犏牛雄性不育的机制和ESCO2基因在雄性生殖系统中的调控机制提供了试验依据。

参考文献
[1]
WU Y, ZHANG W X, ZUO F, et al. Comparison of mRNA expression from Y-chromosome X-degenerate region genes in taurine cattle, yaks and interspecific hybrid bulls[J]. Anim Genet, 2019, 50(6): 740-743. DOI:10.1111/age.12841
[2]
XU C F, SHAH M A, MIPAM T, et al. Bovid microRNAs involved in the process of spermatogonia differentiation into spermatocytes[J]. Int J Biol Sci, 2020, 16(2): 239-250. DOI:10.7150/ijbs.38232
[3]
王泳杰, 王之盛, 胡瑞, 等. 不同品种肉牛产肉性能、牛肉营养品质及风味物质的比较[J]. 动物营养学报, 2019, 31(8): 3621-3631.
WANG Y J, WANG Z S, HU R, et al. Comparison of meat performance, nutritional quality and flavor substance in beef of different breeds cattle[J]. Chinese Journal of Animal Nutrition, 2019, 31(8): 3621-3631. (in Chinese)
[4]
QIU Q, ZHANG G J, MA T, et al. The yak genome and adaptation to life at high altitude[J]. Nat Genet, 2012, 44(8): 946-949. DOI:10.1038/ng.2343
[5]
阿果约达, 熊显荣, 王艳, 等. 犏牛雄性不育相关lncRNA的鉴定与分析[J]. 畜牧兽医学报, 2019, 50(3): 551-561.
AGUO Y D, XIONG X R, WANG Y, et al. Identification and analysis of long non-coding RNA associated with cattle-yak male infertility[J]. Acta Veterinaria et Zootechnica Sinica, 2019, 50(3): 551-561. (in Chinese)
[6]
ALOMER R M, DA SILVA E M L, CHEN J R, et al. Esco1 and Esco2 regulate distinct cohesin functions during cell cycle progression[J]. Proc Natl Acad Sci U S A, 2017, 114(37): 9906-9911. DOI:10.1073/pnas.1708291114
[7]
BENDER D, DA SILVA E M L, CHEN J R, et al. Multivalent interaction of ESCO2 with the replication machinery is required for sister chromatid cohesion in vertebrates[J]. Proc Natl Acad Sci U S A, 2020, 117(2): 1081-1089. DOI:10.1073/pnas.1911936117
[8]
LADURNER R, KREIDL E, IVANOV M P, et al. Sororin actively maintains sister chromatid cohesion[J]. EMBO J, 2016, 35(6): 635-653. DOI:10.15252/embj.201592532
[9]
MFAREJ M G, SKIBBENS R V. An ever-changing landscape in Roberts syndrome biology: implications for macromolecular damage[J]. PLoS Genet, 2020, 16(12): e1009219. DOI:10.1371/journal.pgen.1009219
[10]
MCKAY M J, CRAIG J, KALITSIS P, et al. A Roberts syndrome individual with differential genotoxin sensitivity and a DNA damage response defect[J]. Int J Radiat Oncol Biol Phys, 2019, 103(5): 1194-1202. DOI:10.1016/j.ijrobp.2018.11.047
[11]
GORDILLO M, VEGA H, TRAINER A H, et al. The molecular mechanism underlying Roberts syndrome involves loss of ESCO2 acetyltransferase activity[J]. Hum Mol Genet, 2008, 17(14): 2172-2180. DOI:10.1093/hmg/ddn116
[12]
KOUZNETSOVA E, KANNO T, KARLBERG T, et al. Sister chromatid cohesion establishment factor ESCO1 operates by substrate-assisted catalysis[J]. Structure, 2016, 24(5): 789-796. DOI:10.1016/j.str.2016.03.021
[13]
HOGARTH C A, MITCHELL D, EVANOFF R, et al. Identification and expression of potential regulators of the mammalian mitotic-to-meiotic transition[J]. Biol Reprod, 2011, 84(1): 34-42. DOI:10.1095/biolreprod.110.086215
[14]
LI H, YOU L J, TIAN Y F, et al. DPAGT1-mediated protein N-Glycosylation is indispensable for oocyte and follicle development in mice[J]. Adv Sci (Weinh), 2020, 7(14): 2000531. DOI:10.1002/advs.202000531
[15]
LALDINSANGI C, VIJAYAPRASADARAO K, RAJAKUMAR A, et al. Two-dimensional proteomic analysis of gonads of air-breathing catfish, Clarias batrachus after the exposure of endosulfan and malathion[J]. Environ Toxicol Pharmacol, 2014, 37(3): 1006-1014. DOI:10.1016/j.etap.2014.03.007
[16]
EVANS E B, HOGARTH C, EVANOFF R M, et al. Localization and regulation of murine Esco2 during male and female meiosis[J]. Biol Reprod, 2012, 87(3): 61.
[17]
LU Y J, CHEN Y, CUI Z K, et al. Distinct roles of cohesin acetyltransferases Esco1 and Esco2 in porcine oocyte meiosis I[J]. Cell Cycle, 2019, 18(19): 2481-2494. DOI:10.1080/15384101.2019.1651162
[18]
LU Y J, DAI X X, ZHANG M Q, et al. Cohesin acetyltransferase Esco2 regulates SAC and kinetochore functions via maintaining H4K16 acetylation during mouse oocyte meiosis[J]. Nucleic Acids Res, 2017, 45(16): 9388-9397. DOI:10.1093/nar/gkx563
[19]
MCNICOLL F, KVHNEL A, BISWAS U, et al. Meiotic sex chromosome cohesion and autosomal synapsis are supported by Esco2[J]. Life Sci Alliance, 2020, 3(3): e201900564. DOI:10.26508/lsa.201900564
[20]
LUENSE L J, DONAHUE G, LIN-SHIAO E, et al. Gcn5-mediated histone acetylation governs nucleosome dynamics in spermiogenesis[J]. Dev Cell, 2019, 51(6): 745-758. DOI:10.1016/j.devcel.2019.10.024
[21]
AKSNES H, REE R, ARNESEN T. Co-translational, post-translational, and non-catalytic roles of N-terminal acetyltransferases[J]. Mol Cell, 2019, 73(6): 1097-1114. DOI:10.1016/j.molcel.2019.02.007
[22]
CHEN Y, ZHOU Y L, YIN H S. Recent advances in biosensor for histone acetyltransferase detection[J]. Biosens Bioelectron, 2021, 175: 112880. DOI:10.1016/j.bios.2020.112880
[23]
RIVERA-COLÓN Y, MAGUIRE A, LISZCZAK G P, et al. Molecular basis for cohesin acetylation by establishment of sister chromatid cohesion N-acetyltransferase ESCO1[J]. J Biol Chem, 2016, 291(51): 26468-26477. DOI:10.1074/jbc.M116.752220
[24]
ALEXEEV V, IGOUCHEVA O, DOMASHENKO A, et al. Localized in vivo genotypic and phenotypic correction of the albino mutation in skin by RNA-DNA oligonucleotide[J]. Nat Biotechnol, 2000, 18(1): 43-47. DOI:10.1038/71901
[25]
KOBLAN L W, ERDOS M R, WILSON C, et al. In vivo base editing rescues Hutchinson-Gilford progeria syndrome in mice[J]. Nature, 2021, 589(7843): 608-614. DOI:10.1038/s41586-020-03086-7
[26]
PECKER L H, LANZKRON S. Sickle cell disease[J]. Ann Intern Med, 2021, 174(1): ITC1-ITC16. DOI:10.7326/AITC202101190
[27]
VEGA H, WAISFISZ Q, GORDILLO M, et al. Roberts syndrome is caused by mutations in ESCO2, a human homolog of yeast ECO1 that is essential for the establishment of sister chromatid cohesion[J]. Nat Genet, 2005, 37(5): 468-470. DOI:10.1038/ng1548
[28]
CHENG H Z, ZHANG N G, PATI D. Cohesin subunit RAD21:from biology to disease[J]. Gene, 2020, 758: 144966. DOI:10.1016/j.gene.2020.144966
[29]
ZHANG J, FENG C, SU H D, et al. The cohesin complex subunit ZmSMC3 participates in meiotic centromere pairing in maize[J]. Plant Cell, 2020, 32(4): 1323-1336. DOI:10.1105/tpc.19.00834
[30]
AL-JOMAH N, MUKOLOLO L, ANJUM A, et al. Pds5A and Pds5B display non-redundant functions in mitosis and their loss triggers Chk1 activation[J]. Front Cell Dev Biol, 2020, 8: 531. DOI:10.3389/fcell.2020.00531
[31]
DEARDORFF M A, BANDO M, NAKATO R, et al. HDAC8 mutations in Cornelia de Lange syndrome affect the cohesin acetylation cycle[J]. Nature, 2012, 489(7415): 313-317. DOI:10.1038/nature11316
[32]
KULESZEWICZ K, FU X W, KUDO N R. Cohesin loading factor Nipbl localizes to chromosome axes during mammalian meiotic prophase[J]. Cell Div, 2013, 8(1): 12. DOI:10.1186/1747-1028-8-12
[33]
FAGERBERG L, HALLSTRÖM B M, OKSVOLD P, et al. Analysis of the human tissue-specific expression by genome-wide integration of transcriptomics and antibody-based proteomics[J]. Mol Cell Proteomics, 2014, 13(2): 397-406. DOI:10.1074/mcp.M113.035600
[34]
YU Y, FUSCOE J C, ZHAO C, et al. A rat RNA-Seq transcriptomic BodyMap across 11 organs and 4 developmental stages[J]. Nat Commun, 2014, 5: 3230. DOI:10.1038/ncomms4230
[35]
CHEN H M, ZHANG L, HE W T, et al. ESCO2 knockdown inhibits cell proliferation and induces apoptosis in human gastric cancer cells[J]. Biochem Biophys Res Commun, 2018, 496(2): 475-481. DOI:10.1016/j.bbrc.2018.01.048
[36]
GUO X B, HUANG B, PAN Y H, et al. ESCO2 inhibits tumor metastasis via transcriptionally repressing MMP2 in colorectal cancer[J]. Cancer Manag Res, 2018, 10: 6157-6166. DOI:10.2147/CMAR.S181265
[37]
LOU Y N, LIU W J, WANG C L, et al. Histological evaluation and Prdm9 expression level in the testis of sterile male cattle-yaks[J]. Livest Sci, 2014, 160: 208-213. DOI:10.1016/j.livsci.2013.12.017
[38]
SATO Y, KURIWAKI R, HAGINO S, et al. Abnormal functions of Leydig cells in crossbred cattle-yak showing infertility[J]. Reprod Domest Anim, 2020, 55(2): 209-216. DOI:10.1111/rda.13609
[39]
YAN P, XIANG L, GUO X, et al. The low expression of Dmrt7 is associated with spermatogenic arrest in cattle-yak[J]. Mol Biol Rep, 2014, 41(11): 7255-7263. DOI:10.1007/s11033-014-3611-x
[40]
李波, 庄梦茹, 李娜, 等. 奶山羊ESCO2基因克隆及其功能的初步研究[J]. 农业生物技术学报, 2016, 24(2): 176-185.
LI B, ZHUANG M R, LI N, et al. Molecular cloning of ESCO2 gene and the preliminary study of its function in dairy goat (Capra hircus)[J]. Journal of Agricultural Biotechnology, 2016, 24(2): 176-185. (in Chinese)

(编辑   郭云雁)