畜牧兽医学报  2022, Vol. 53 Issue (2): 481-492. DOI: 10.11843/j.issn.0366-6964.2022.02.015    PDF    
犏牛PPP1R11基因的克隆及在睾丸中的表达规律研究
闵星宇1, 杨丽雪2, 于海玲1, 胡宇磊3, 杨满珍3, 杨璐瑜3, 李键1,2, 熊显荣1,2     
1. 西南民族大学畜牧兽医学院, 成都 610041;
2. 青藏高原动物遗传育种资源保护与利用国家教育部重点实验室, 成都 610041;
3. 动物科学国家民委重点实验室, 成都 610041
摘要:旨在克隆犏牛蛋白磷酸酶1调节亚基11(protein phosphatase 1 regulatory subunit 11,PPP1R11)基因,分析其在不同发育阶段睾丸中的表达与定位,为解析其在雄性生殖中的功能机制提供理论依据。本研究采集成年犏牛睾丸、附睾、心、肝、脾、肺、肾、大肠、小肠、胃、肌肉和脂肪组织(n=3),犏牛和牦牛在胎牛时期(5~6月龄)、幼年时期(1~2岁)和成年时期(3~4岁)睾丸(n=3)作为试验材料。通过RT-PCR技术克隆犏牛PPP1R11基因并进行生物信息学分析,采用实时荧光定量PCR(quantitative real-time PCR,qRT-PCR)构建PPP1R11基因在犏牛不同组织中的表达谱,比较分析PPP1R11在犏牛和牦牛不同时期睾丸中的表达,利用免疫组织化学(immunohistochemistry,IHC)技术检测PPP1R11蛋白的细胞定位和表达。结果显示,犏牛PPP1R11基因CDS区为324 bp,编码107个氨基酸,犏牛PPP1R11蛋白与其他哺乳动物同源性较高;PPP1R11可与PPP1R2、PPP1R7、PPP1CB和UTP20等蛋白相互作用,互作蛋白功能主要与蛋白质的磷酸化相关,参与调控睾丸发育、精子发生和性成熟等生物学过程。PPP1R11在犏牛组织中广泛表达,其中在睾丸组织中的表达量最高,相对表达水平极显著高于其它组织(P < 0.01);该基因在犏牛睾丸中的表达量随年龄增长呈上升趋势,幼年和成年时期犏牛睾丸中PPP1R11的表达与同时期牦牛相比差异极显著(P < 0.01);IHC染色结果发现,PPP1R11蛋白在犏牛精原细胞和支持细胞中低表达并与牦牛存在显著差异,雄性犏牛初级精母细胞数量减少并且减数分裂阻滞于此。本研究表明,犏牛与牦牛PPP1R11在睾丸中存在时空表达差异,提示该基因可能与犏牛雄性不育相关,其具体作用机制有待进一步研究。
关键词犏牛    PPP1R11    基因克隆    基因表达    减数分裂    睾丸发育    
Cloning and Expression Analysis of Cattle-yak PPP1R11 Gene in Testis
MIN Xingyu1, YANG Lixue2, YU Hailing1, HU Yulei3, YANG Manzhen3, 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: This study was conducted to clone the PPP1R11 gene in cattle-yak, detect the expression and location in testes of cattle-yak at different developmental stages, as to provide a theoretical basis for further research of the mechanism of function in male reproduction. The testis, epididymis, heart, liver, spleen, lung, kidney, large intestine, small intestine, stomach, muscle and adipose tissue of adult cattle-yak were collected (n=3), testicles of cattle-yak and yak in the fetal period (5-6 months old), juvenile period (1-2 years old) and adult period (3-4 years old) (n=3) were collected as experimental materials. The PPP1R11 CDS region of cattle-yak was cloned by RT-PCR and analyzed by bioinformatics softwares. qRT-PCR was used to detect the expression patterns of PPP1R11 mRNA in different tissues of cattle-yak, and compared the expression difference in testes at different development stages between cattle-yak and yak. The cell location and expression of PPP1R11 protein were detected by immunohistochemistry (IHC) staining. The results showed that the CDS region of PPP1R11 gene in cattle-yak was 324 bp, which encoded 107 amino acids. The corresponding PPP1R11 protein of cattle-yak had high homology with other mammals. The protein function prediction showed that PPP1R11 could interact with PPP1R2, PPP1R7, PPP1CB, UTP20 and other proteins, and mainly related to the phosphorylation of proteins, which could regulate biological processes such as testicular development, spermatogenesis and sexual maturity. PPP1R11 gene was wildly expressed in various tissues of cattle-yak, and the relative expression level in testis was significantly higher than other tissues (P < 0.01). The expression of PPP1R11 gene in cattle-yak testis increased with age. Furthermore, there was significant difference of the expression level of PPP1R11 (P < 0.01) in testis at juvenile and adult stage of cattle-yak compared with the counterpart of yak. In addition, IHC staining result showed that PPP1R11 protein was significantly lower expressed in cattle-yak spermatogonia and sertoli cells, and there was significant difference compared with yak, and the number of primary spermatocytes in cattle-yak was reduced significantly and meiosis arrested at this stage. This study showed that there were temporal and spatial differences in the expression of PPP1R11 between cattle-yak and yak testis, suggesting that PPP1R11 may be related to to male cattle-yak infertility. However, its specific mechanism needs to be further studied.
Key words: cattle-yak    PPP1R11    gene cloning    gene expression    meiosis    testis development    

犏牛为母牦牛(Bos grunniens)与平原牛(Bos taurus)杂交得到的F1代,主要分布于我国青藏高原地区的优良家畜资源。牦牛生产性能低下,已无法满足青藏高原地区日益增长的奶肉需求,而平原牛又无法适应高原恶劣的高寒低氧环境。犏牛具有显著的杂种优势,产奶性能和产肉性能均优于牦牛,并有良好的高原适应性[1-3]。然而,雄性犏牛精子发生阻滞,使得无法通过横交固定其杂种优势,这严重制约了杂种优势的利用[4-5]。因此,开展犏牛繁殖生理的研究对提高青藏高原畜牧生产水平、改善牧民生活品质和发展青藏高原特色畜牧产业都具有重要的价值。

丝氨酸/苏氨酸蛋白磷酸酶1(serine/threonine protein phosphatase,PP1)是一种去磷酸化酶,能够通过翻译后调控参与多种生物学过程[6-7]。在哺乳动物中共发现4种亚型,包括PP1α、PP1β、PP1γ1和PP1γ2,其中PP1γ2蛋白由蛋白磷酸酶1催化亚基γ(protein phosphatase 1 catalytic subunit gamma, PPP1CC)基因编码,主要富集于睾丸组织并参与精子细胞发生和成熟[8-9],敲除PPP1CC基因造成雄性不育[10-11]。蛋白磷酸酶1调节亚基11(protein phosphatase 1 regulatory subunit 11,PPP1R11) 最早通过酵母双杂交试验鉴定到是PP1的结合蛋白,又被称作INH3、TCTEX5、HCGV、IPP3,能够抑制PP1γ2的去磷酸化作用[12-13]。前期研究表明,PPP1R11基因是t-complex的重要组成部分,而t-complex与精子活力受损和雄性不育有关[14]。Northern印迹杂交和免疫组织荧光染色显示,PPP1R11 mRNA和蛋白主要分布在小鼠睾丸组织中,其表达水平显著高于其他组织,在雄性不育小鼠中检测到PPP1R11基因的下调,推测该基因的下调导致生殖细胞的凋亡[13, 15]。可见,PPP1R11基因在调控雄性配子发生和精子成熟中有重要作用。

目前,关于PPP1R11基因的研究主要集中在小鼠上,而在雄性不育犏牛上的研究还尚未报道。本研究以犏牛为研究对象,利用RT-PCR技术克隆获得犏牛PPP1R11编码区完整序列,并对该序列进行生物信息学分析,同时采用qRT-PCR和IHC染色检测该基因和蛋白在犏牛各组织及各发育阶段睾丸组织中的表达,以期为研究PPP1R11调控精子发生机制提供参考,为进一步从分子水平解析雄性犏牛不育机制奠定基础并提供理论依据。

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

将胎牛时期(5~6月龄)、幼年时期(1~2岁)、成年时期(3~4岁)犏牛作为本研究试验对象,样品采自成都市周边屠宰场。健康的犏牛屠宰后迅速采集睾丸、附睾、心、肝、脾、肺、肾、大肠、小肠、胃、肌肉和脂肪。以上所有样品每组各取3头为生物学重复,样品使用高压灭菌处理的生理盐水冲洗剪成0.5 cm× 0.5 cm×0.5 cm的组织块,一部分放入冻存管,置于液氮罐中带回实验室,于-80 ℃保存;另一部分使用4%多聚甲醛固定,用于IHC染色检测。

1.2 犏牛和牦牛组织总RNA的提取与反转录

冷冻的犏牛和牦牛各组织使用高通量组织研磨器研磨后,按照Trizol法,经离心、沉淀、清洗和溶解等步骤进行总RNA的提取,使用紫外分光光度计(Biospec-nano,日本)检测总RNA的浓度和纯度,选取OD260 nm/OD280 nm介于1.8~2.0的RNA,并经1%琼脂糖凝胶电泳检测其完整性,置于-80 ℃保存备用。根据Revert Aid First Strand cDNA Synthesis Kit反转录试剂盒(Thermo scientific,美国)说明书合成第一链cDNA,置于-20 ℃保存。

1.3 引物的设计与基因克隆

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

表 1 PCR引物序列及产物信息 Table 1 PCR primers sequence and products length
1.4 PPP1R11基因的生物信息学分析

利用ORF Finder查找犏牛PPP1R11克隆序列开放阅读框并翻译成氨基酸序列。采用T-COFFEE在线软件并结合DNAMAN 9.0软件对不同物种PPP1R11的氨基酸序列同源性进行比较,通过MEME在线工具分析查找具有高度保守性的基序,并用MEGA 7.0软件构建系统发生树。利用ExPASy ProtParam、TMHMM、PSORTⅡ和SignalP在线软件进行犏牛PPP1R11蛋白基本理化性质、跨膜区、亚细胞定位和信号肽分析。糖基化位点和磷酸化位点预测使用NetPhos 3.1、NetOGlyc 4.0和NetNGlyc 1.0软件在线分析。该蛋白二级结构、三级结构和互作蛋白通过SOPMA、VMD 1.9.3和STRING 11.0预测。

1.5 犏牛PPP1R11基因的表达分析

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

1.6 免疫组织化学染色

将不同时期睾丸组织从固定液中取出,制作石蜡切片。切片经二苯甲中脱蜡,酒精洗涤,EDTA抗原修复,PBS洗涤后,放入3%双氧水阻断内源性过氧化物酶;5%胎牛血清室温封闭30 min,滴加一抗(多克隆兔抗PPP1R11,1∶200稀释,博奥森)后4 ℃孵育过夜,阴性对照不加一抗;滴加二抗(多聚化山羊抗兔IgG,辣根过氧化物酶标记),室温孵育50 min;滴加DAB显色液显色,自来水冲洗切片终止显色;苏木精复染3 min左右,复染细胞核;中性树胶封片。使用激光共聚焦显微镜(Zeiss,德国)观察并拍照。

1.7 数据统计分析

应用2-ΔΔCt法分析qRT-PCR检测结果,结果使用“平均值±标准误(Mean ± SEM)”表示。使用SPSS 19.0软件的T-test检验进行数据分析,P<0.05代表差异显著,P<0.01代表差异极显著。

2 结果 2.1 犏牛PPP1R11基因克隆及生物信息学分析

2.1.1 犏牛PPP1R11克隆   以犏牛睾丸组织cDNA为模板,F1、R1为上、下游引物,PCR扩增犏牛PPP1R11基因编码区序列,1.5%琼脂糖凝胶电泳检测显示,目标条带与预期结果相符(图 1A)。经测序得到893 bp的核苷酸序列,其中犏牛PPP1R11基因CDS为324 bp,共编码107个氨基酸(图 1B)。

A: M. DNA相对分子质量标准DL2000; 1. 犏牛PPP1R11基因扩增产物;B.犏牛PPP1R11基因核苷酸及其推测的氨基酸序列 A: M. DL2000 marker; 1.PCR product of cattle-yak PPP1R11 gene; B. Nucleotide sequence and the deduced amino acid sequence of cattle-yak PPP1R11 图 1 犏牛PPP1R11基因的PCR扩增结果 Fig. 1 PCR amplification results of cattle-yak PPP1R11 gene

2.1.2 犏牛PPP1R11同源性比较与系统进化树   对不同哺乳动物PPP1R11的氨基酸序列同源性分析发现,PPP1R11在不同哺乳动物中高度保守,其中犏牛与牦牛和黄牛同源性最高(100%),其次为瘤牛(99.07%)、水牛(99.07%)、绵羊(99.07%)、山羊(99.07%)、羊驼(98.06%)、猪(97.09%)、马(97.09%)、大狐蝠(96.12%)、蓝鲸(95.28%)、港海豹(95.24%)、智人(95.15%),与鸭嘴兽的同源性最低(84.88%)(图 2A)。通过MEME软件在线分析获得PPP1R11蛋白最保守的2个基序(图 2B),分别位于犏牛PPP1R11氨基酸序列N端第6~55和58~107位点。基于15种哺乳动物PPP1R11不同氨基酸序列,采用邻接算法构建系统进化树(图 3),发现犏牛和黄牛聚为一支,在进化上与牦牛、瘤牛、水牛、绵羊和山羊亲缘关系较近,与非反刍类动物亲缘关系较远,这与同源性比较结果一致。

A. 15种不同哺乳动物PPP1R11氨基酸序列的多重比较;B. PPP1R11氨基酸序列中2个最保守基序 A. Multiple alignment of PPP1R11 amino acid sequences from 15 different mammalian species; B. Two of most conserved motifs in PPP1R11 amino acid sequence 图 2 不同哺乳动物PPP1R11氨基酸序列对比 Fig. 2 Alignment of amino acid sequence of PPP1R11 among different mammals
图 3 PPP1R11蛋白氨基酸序列系统进化树 Fig. 3 Phylogenetic tree for PPP1R11 amino acid sequences

2.1.3 犏牛PPP1R11理化性质   犏牛PPP1R11蛋白的分子量、理论等电点和分子式分别为12.014 ku、8.33和C503H803N163O164S8;氨基酸组成中,脯氨酸(Pro)占比最高,其含量为13.1%;带负电荷(Asp + Glu)和带正电荷(Arg + Lys)的氨基酸残基总数分别16和18个,表明该蛋白可能带正电。糖基化位点预测表明,犏牛PPP1R11蛋白含1个N糖基化位点和18个O糖基化位点;磷酸化位点预测显示,该蛋白有19个磷酸化位点,其中丝氨酸(Ser)磷酸化位点11个,苏氨酸(Thr)磷酸化位点7个,酪氨酸(Tyr)磷酸化位点1个。该蛋白亚细胞定位主要存在于细胞核(47.8%),其次为线粒体(34.8%)和细胞质(17.4%);不含信号肽且无跨膜结构域。

2.1.4 犏牛PPP1R11结构与互作蛋白预测   采用SOPMA预测PPP1R11蛋白的二级结构,其中18个氨基酸(16.82%)形成α螺旋,16个氨基酸(14.95%)形成延伸链,9个氨基酸(8.41%)形成β转角,64个氨基酸(59.81%)形成无规则卷曲(图 4A);通过VMD 1.9.3预测到犏牛PPP1R11蛋白三级结构(图 4B)。采用STRING数据库分析犏牛PPP1R11潜在相互作用蛋白显示,PPP1R11可能与PPP1R2、PPP1R7、PPP1CA、PPP1CB、PPP1CC、PPA2、CIRH1A和UTP20等10个蛋白存在相互作用(图 4C)。PPP1R11蛋白处于互作网络核心位置,且与PPP1R7、PPP1CB、PPP1R2、UTP20蛋白间的相互作用具有更高可信度。

A. 犏牛PPP1R11蛋白二级结构;B. 犏牛PPP1R11蛋白三级结构;C. 犏牛PPP1R11的蛋白互作网络 A. Secondary structure of cattle-yak PPP1R11 protein; B. Tertiary structure of cattle-yak PPP1R11 protein; C. Interaction network for cattle-yak PPP1R11 protein 图 4 犏牛PPP1R11蛋白的空间结构及互作蛋白 Fig. 4 Spatial structure and interaction proteins of cattle-yak PPP1R11
2.2 犏牛PPP1R11基因的表达模式

2.2.1 犏牛PPP1R11组织表达谱   以β-actin作为内参基因,利用qRT-PCR技术检测PPP1R11 mRNA在犏牛睾丸、附睾、心、肝、脾、肺、肾、大肠、小肠、胃、肌肉和脂肪组织中的表达量(图 5)。结果显示,该基因在犏牛12个组织中均有表达,在睾丸、脾和胃组织中表达量较高,其在睾丸中的表达水平极显著高于其它各个组织(P<0.01)。

不同的字母表示差异极显著(P<0.01)。下同 Different letters represent the significant differences(P < 0.01). The same as below 图 5 PPP1R11 mRNA在犏牛中的组织表达谱 Fig. 5 Tissue expression profile of PPP1R11 mRNA in cattle-yak

2.2.2 犏牛睾丸中PPP1R11的表达规律   采用qRT-PCR检测PPP1R11 mRNA在犏牛不同发育时期睾丸组织中的相对表达水平(图 6)。结果显示,该基因在犏牛睾丸组织中的表达随年龄的增加呈上升趋势;5~6月龄胎牛睾丸组织中几乎不表达,表达量低于同龄的牦牛,但差异不显著(P>0.05);1~2岁和3~4岁犏牛睾丸中的表达水平极显著低于同时期牦牛(P<0.01)。

图 6 PPP1R11 mRNA在不同发育阶段犏牛睾丸中的表达规律 Fig. 6 Expression patterns of PPP1R11 mRNA in testes at different developmental stages of cattle-yak

2.2.3 犏牛PPP1R11的细胞定位   基于犏牛和牦牛PPP1R11 mRNA在睾丸中的表达差异极显著,本试验进一步通过免疫组织化学染色方法检测PPP1R11蛋白在犏牛和牦牛不同年龄段睾丸组织中的细胞定位和表达情况(图 7)。结果显示,犏牛和牦牛5~6月龄胎牛睾丸组织生精小管形态无明显差异;1~2岁和3~4岁犏牛与牦牛睾丸组织中,犏牛生精小管管径缩小,空腔增大,生精上皮皱缩,细胞组成差异显著,仅有精原细胞、支持细胞和少量初级精母细胞(PS)松散分布,无次级精母细胞(SS)和精子细胞。IHC染色发现,PPP1R11蛋白在犏牛和牦牛各年龄段睾丸组织中均有表达。其中5~6月龄犏胎牛睾丸组织中,PPP1R11蛋白主要定位于精原干细胞(SSC),支持细胞(SC)几乎不表达,与同时期牦牛具有相同的蛋白定位;1~2岁和3~4岁犏牛睾丸组织中,PPP1R11蛋白主要定位于初级精母细胞和间质细胞(IC),精原细胞和支持细胞中蛋白阳性信号较弱,而牦牛各级生殖细胞、支持细胞和间质细胞中均检测到PPP1R11阳性信号。

A. 犏胎牛睾丸(5~6月龄);B. 幼年犏牛睾丸(1~2岁);C. 成年犏牛睾丸(3~4岁);D. 牦胎牛睾丸(5~6月龄);E. 幼年牦牛睾丸(1~2岁);F. 成年牦牛睾丸(3~4岁);G. 成年犏牛睾丸(阴性对照);H. 成年牦牛睾丸(阴性对照)。SP. 精原细胞;PS. 初级精母细胞;SS. 次级精母细胞;SC. 支持细胞;IC. 间质细胞;RS. 圆形精子细胞;ES. 长形精子细胞 A. Fetal cattle-yak testis (5-6 months); B. Juvenile cattle-yak testis (1-2 years old); C. Adult cattle-yak testis (3-4 years old); D. Fetal yak testis (5-6 months); E. Juvenile yak testis (1-2 years old); F. Adult yak testis (3-4 years old); G. Adult cattle-yak testis (Negative control); H. Adult yak testis (Negative control). SP. Spermatogonium; PS. Primary spermatocyte; SS. Secondary spermatocyte; SC. Sertoli cell; IC. Leydig cell; RS. Round spermatid; ES. Elongating spermatid 图 7 不同阶段雄性犏牛生殖细胞发育及PPP1R11蛋白定位 Fig. 7 The growth of male cattle-yak germ cells and location of PPP1R11 protein at different stages
3 讨论

雄性生殖力是种畜的重要评价指标,性状优良的种公畜往往能够产生更大的经济效益,故而提高种公畜生殖力对养殖业的发展具有促进作用。犏牛是牦牛和黄牛的杂交后代,具有杂种优势显著和雄性不育的特点[16],其雄性不育犏牛可作为研究雄性生殖调控机理的天然理想模型。据报道,蛋白质磷酸化与去磷酸化参与了30%~70%的细胞信号转导过程[6],精子发生是一个高度动态变化的复杂过程,需要多种磷酸化酶的参与[17]。其中,去磷酸化酶PP1γ2能够直接调控哺乳动物的精子发生,而PPP1R11能够通过抑制PP1γ2的表达来维持精子的正常发生和成熟[15, 18]。因此,探索PPP1R11在睾丸中的表达规律对于研究生殖调控和解析犏牛雄性不育机制具有重要的价值。

本研究成功克隆了犏牛PPP1R11基因CDS区全长序列,由324 bp的核苷酸组成,编码107个氨基酸,这与NCBI中提供的野牦牛预测序列(XM_014483599.1)和已报道黄牛[19]的序列完全一致,该基因CDS区长度与瘤牛、水牛、山羊和绵羊一致。通过氨基酸序列对比发现,犏牛PPP1R11蛋白与其它14种哺乳动物同源性在95%以上(鸭嘴兽除外),并识别到2个最为保守的基序。这些结果表明,犏牛PPP1R11与其它哺乳动物相比具有高度同源性和进化上的高度保守性,PPP1R11可能在不同物种中发挥较为一致的生物学功能。鉴于PPP1R11主要通过与相应蛋白互作发挥功能,本研究进一步分析了犏牛PPP1R11的互作蛋白发现,该蛋白与PPP1R2、PPP1CB等10个蛋白质分子存在相互作用。其中PPP1CA、PPP1CB、PPP1CC和PP2A互作为不同类型蛋白磷酸酶,被广泛报道能够参与调控哺乳动物睾丸发育、精子发生和成熟[20-23]。有研究表明,在小鼠原始生殖细胞中特异性敲除PP2A可导致雄性小鼠不育[24],通过转录组测序技术检测到PP2A在雄性不育犏牛睾丸中也表现为下调[25],此外PP2A还可以通过调控减数分裂来调节卵母细胞的成熟[26]。PPP1R2、PPP1R7和PPP1R11均为PP1的调节剂,能够在小鼠附睾中共同调控精子的成熟[15],PPP1R11能够通过抑制PPP1R7与PP1结合从而确保有丝分裂过程中动粒与纺锤体的结合[27],而PPP1R11与PPP1R2能否直接互作还未见报道。有研究表明,UTP20能激活RNA聚合酶Ⅰ转录并促进细胞增殖[28],推测PPP1R11可能与UTP20互作从而参与调控细胞增殖。综上,这些潜在与犏牛PPP1R11互作的蛋白分子均在睾丸发育和精子发生中发挥了重要的作用,这为深入理解PPP1R11在哺乳动物睾丸发育及精子发生中的分子机制提供了思路和见解。

通过qRT-PCR检测了PPP1R11基因在犏牛不同组织中的表达,组织表达谱分析显示,该基因在犏牛各组织中均有表达,在睾丸组织中表达最高,这与小鼠[29]、大鼠[30]和绵羊[31]PPP1R11基因组织表达谱相一致,PPP1R11基因在多种哺乳动物睾丸组织中的高表达表明该基因与雄性生殖相关,且在犏牛睾丸中也可能发挥了相似的功能。Goswami等[15]研究表明,PPP1R11基因在小鼠附睾中能够促进精子的成熟,而犏牛附睾中该基因表达水平较低,推测是因为犏牛精子发生阻滞,附睾中无精子细胞可成熟。此外,本试验还检测到该基因在脾和胃中有较高表达,推测该基因除参与睾丸发育、精子生成和成熟外,在其他组织中也可能发挥作用,其具体功能有待进一步研究。

本试验检测到PPP1R11在犏牛睾丸组织中的表达水平随睾丸生长发育呈上升趋势,与同时期牦牛相比均较低,且在幼年和成年时期差异极显著。有研究表明,过表达PPP1R11能够抑制细胞的凋亡[32],在T-单倍型突变雄性不育小鼠模型中,PPP1R11也表现为下调[33],由此推测PPP1R11在睾丸中的表达水平与雄性生殖力成正相关,而雄性犏牛不育也可能与此有关联。为进一步探讨PPP1R11基因在犏牛睾丸中发挥的作用,利用IHC染色检测PPP1R11蛋白在犏牛3个时期睾丸中的表达和细胞定位。成年期犏牛生精小管内只观察到精原细胞、支持细胞和少量初级精母细胞,未观察到次级精母细胞和精子细胞,说明雄性犏牛精子发生阻滞于初级精母细胞,这与前人研究结果一致[34-36]。在幼年期和成年期犏牛睾丸组织中,PPP1R11蛋白主要定位于犏牛初级精母细胞和间质细胞,其蛋白阳性信号强度显著低于同时期牦牛,这与qRT-PCR的结果一致。正常生殖力小鼠睾丸中,PPP1R11蛋白主要定位于圆形期精子和长形期精子,其他各级生殖细胞和支持细胞亦有表达[13-15],正常生殖力牦牛生精小管中各级生殖细胞均检测到PPP1R11蛋白阳性信号,本试验中犏牛精原细胞几乎不表达PPP1R11蛋白,这与小鼠和牦牛的结果不一致,推测该蛋白在犏牛精原细胞中的下调导致初级精母细胞生成的减少,从而造成精子发生阻滞。综上,PPP1R11基因可能参与调控哺乳动物精子发生,其mRNA和蛋白在犏牛睾丸中的差异表达可能是造成雄性犏牛减数分裂阻滞的原因之一。

4 结论

本试验成功克隆了犏牛PPP1R11基因序列,该基因在犏牛与牦牛睾丸组织中的mRNA表达水平和蛋白定位存在差异,提示PPP1R11可能与维持雄性正常生殖力相关。本研究结果为进一步探究雄性犏牛不育的机制和PPP1R11在雄性生殖系统中的生理功能和调控机制提供了理论依据。

参考文献
[1]
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
[2]
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
[3]
SONG C C, HUANG Y Z, YANG Z X, et al. RNA-Seq analysis identifies differentially expressed genes insubcutaneous adipose tissuein Qaidamford cattle, cattle-yak, and angus cattle[J]. Animals (Basel), 2019, 9(12): 1077.
[4]
NIAYALE R, CUI Y, ADZITEY F. Male hybrid sterility in the cattle-yak and other bovines: a review[J]. Biol Reprod, 2021, 104(3): 495-507. DOI:10.1093/biolre/ioaa207
[5]
阿果约达, 熊显荣, 王艳, 等. 犏牛雄性不育相关lncRNA的鉴定与分析[J]. 畜牧兽医学报, 2019, 50(3): 551-561.
A GUO 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]
BOLLEN M, PETI W, RAGUSA M J, et al. The extended PP1 toolkit: designed to create specificity[J]. Trends Biochem Sci, 2010, 35(8): 450-458. DOI:10.1016/j.tibs.2010.03.002
[7]
UBERSAX J A, FERRELL JR J E. Mechanisms of specificity in protein phosphorylation[J]. Nat Rev Mol Cell Biol, 2007, 8(7): 530-541. DOI:10.1038/nrm2203
[8]
DUDIKI T, JOUDEH N, SINHA N, et al. The protein phosphatase isoform PP1γ1 substitutes for PP1γ2 to support spermatogenesis but not normal sperm function and fertility[J]. Biol Reprod, 2019, 100(3): 721-736. DOI:10.1093/biolre/ioy225
[9]
DEY S, BROTHAG C, VIJAYARAGHAVAN S. Signaling enzymes required for sperm maturation and fertilization in mammals[J]. Front Cell Dev Biol, 2019, 7: 341. DOI:10.3389/fcell.2019.00341
[10]
CHAKRABARTI R, CHENG L N, PURI P, et al. Protein phosphatase PP1γ2 in sperm morphogenesis and epididymal initiation of sperm motility[J]. Asian J Androl, 2007, 9(4): 445-452. DOI:10.1111/j.1745-7262.2007.00307.x
[11]
CHAKRABARTI R, KLINE D, LU J, et al. Analysis of Ppp1cc-null mice suggests a role for PP1 gamma2 in sperm morphogenesis[J]. Biol Reprod, 2007, 76(6): 992-1001. DOI:10.1095/biolreprod.106.058610
[12]
LI H H, ROKAVEC M, JIANG L C, et al. Antagonistic effects of p53 and HIF1A on microRNA-34a regulation of PPP1R11 and STAT3 and hypoxia-induced epithelial to mesenchymal transition in colorectal cancer cells[J]. Gastroenterology, 2017, 153(2): 505-520. DOI:10.1053/j.gastro.2017.04.017
[13]
CHENG L N, PILDER S, NAIRN A C, et al. PP1γ2 and PPP1R11 are parts of a multimeric complex in developing testicular germ cells in which their steady state levels are reciprocally related[J]. PLoS One, 2009, 4(3): e4861. DOI:10.1371/journal.pone.0004861
[14]
HAN Y B, FENG H L, CHEUNG C K, et al. Expression of a novel T-complex testis expressed 5 (Tctex5) in mouse testis, epididymis, and spermatozoa[J]. Mol Reprod Dev, 2007, 74(9): 1132-1140. DOI:10.1002/mrd.20631
[15]
GOSWAMI S, KORRODI-GREGÓRIO L, SINHA N, et al. Regulators of the protein phosphatase PP1γ2, PPP1R2, PPP1R7, and PPP1R11 are involved in epididymal sperm maturation[J]. J Cell Physiol, 2019, 234(3): 3105-3118. DOI:10.1002/jcp.27130
[16]
黄向月, 熊显荣, 张磊, 等. SPATA3基因克隆及其在牦牛和犏牛睾丸中的差异表达研究[J]. 畜牧兽医学报, 2019, 51(2): 279-287.
HUANG X Y, XIONG X R, ZHANG L, et al. Cloning of SPATA3 gene and its differential expression analysis in yak and cattle-yak testis[J]. Acta Veterinaria et Zootechnica Sinica, 2019, 51(2): 279-287. (in Chinese)
[17]
ZHANG J Y, MRUK D D, CHENG C Y. Myotubularin phosphoinositide phosphatases, protein phosphatases, and kinases: their roles in junction dynamics and spermatogenesis[J]. J Cell Physiol, 2005, 204(2): 470-483. DOI:10.1002/jcp.20303
[18]
FORGIONE N, VOGL A W, VARMUZA S. Loss of protein phosphatase 1cγ (PPP1CC) leads to impaired spermatogenesis associated with defects in chromatin condensation and acrosome development: an ultrastructural analysis[J]. Reproduction, 2010, 139(6): 1021-1029. DOI:10.1530/REP-10-0063
[19]
ZIMIN A V, DELCHER A L, FLOREA L, et al. A whole-genome assembly of the domestic cow, Bos taurus[J]. Genome Biol, 2009, 10(4): R42. DOI:10.1186/gb-2009-10-4-r42
[20]
DUDIKI T, KADUNGANATTIL S, FERRARA J K, et al. Changes in carboxy methylation and tyrosine phosphorylation of protein phosphatase PP2A are associated with epididymal sperm maturation and motility[J]. PLoS One, 2015, 10(11): e0141961. DOI:10.1371/journal.pone.0141961
[21]
VAN DEN HAM R, VAN DISSEL-EMILIANI F M F, VAN PELT A M M. Expression of the scaffolding subunit A of protein phosphatase 2A during rat testicular development[J]. Biol Reprod, 2003, 68(4): 1369-1375. DOI:10.1095/biolreprod.102.004853
[22]
SILVA J V, FREITAS M J, FARDILHA M. Phosphoprotein phosphatase 1 complexes in spermatogenesis[J]. Curr Mol Pharmacol, 2014, 7(2): 136-146.
[23]
MACLEOD G, VARMUZA S. Tandem affinity purification in transgenic mouse embryonic stem cells identifies DDOST as a novel PPP1CC2 interacting protein[J]. Biochemistry, 2012, 51(48): 9678-9688. DOI:10.1021/bi3010158
[24]
PAN X Y, CHEN X, TONG X, et al. Ppp2ca knockout in mice spermatogenesis[J]. Reproduction, 2015, 149(4): 385-391. DOI:10.1530/REP-14-0231
[25]
CAI X, YU S M, MIPAM T, et al. Comparative analysis of testis transcriptomes associated with male infertility in cattleyak[J]. Theriogenology, 2017, 88: 28-42. DOI:10.1016/j.theriogenology.2016.09.047
[26]
KORNBLUTH S, FISSORE R. Vertebrate reproduction[J]. Cold Spring Harb Perspect Biol, 2015, 7(10): a006064. DOI:10.1101/cshperspect.a006064
[27]
EITENEUER A, SEILER J, WEITH M, et al. Inhibitor-3 ensures bipolar mitotic spindle attachment by limiting association of SDS22 with kinetochore-bound protein phosphatase-1[J]. EMBO J, 2014, 33(22): 2704-2720. DOI:10.15252/embj.201489054
[28]
PENG Q H, WU J G, ZHANG Y, et al. 1A6/DRIM, a novel t-UTP, activates RNA polymerase I transcription and promotes cell proliferation[J]. PLoS One, 2010, 5(12): e14244. DOI:10.1371/journal.pone.0014244
[29]
YUE F, CHENG Y, BRESCHI A, et al. A comparative encyclopedia of DNA elements in the mouse genome[J]. Nature, 2014, 515(7527): 355-364. DOI:10.1038/nature13992
[30]
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(1): 3230. DOI:10.1038/ncomms4230
[31]
JIANG Y, XIE M, CHEN W B, et al. The sheep genome illuminates biology of the rumen and lipid metabolism[J]. Science, 2014, 344(6188): 1168-1173. DOI:10.1126/science.1252806
[32]
LIU G S, WU J J, QIAO M, et al. Transcription elements AREB6 and miR-34a affect apoptosis of PAMs by regulating the expression of SS2-related gene PPP1R11[J]. Cell Cycle, 2019, 18(9): 1033-1044. DOI:10.1080/15384101.2019.1610241
[33]
HAN Y B, SONG X X, FENG H L, et al. Mutations of t-complex testis expressed gene 5 transcripts in the testis of sterile t-haplotype mutant mouse[J]. Asian J Androl, 2008, 10(2): 219-226. DOI:10.1111/j.1745-7262.2008.00323.x
[34]
ZOU S, LIU P, YU S, et al. Hsp60 expression and localization in different tissues and testis development of male cattle (cattle-yak and yak)[J]. Folia Morphol (Warsz), 2021. DOI:10.5603/FM.a2020.0126
[35]
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
[36]
LIU P G, YU S J, CUI Y, et al. Regulation by Hsp27/P53 in testis development and sperm apoptosis of male cattle (cattle-yak and yak)[J]. J Cell Physiol, 2018, 234(1): 650-660.

(编辑   郭云雁)