水牛睪丸生精細胞蛋白質組學及其磷酸化修飾組學聯合分析

DOI：10.3969/j.issn.2095-1191.2025.01.023

摘要：【目的】綜合分析水牛睪丸生精細胞蛋白質組學及其磷酸化修飾組學，探究水牛精子發生過程中蛋白的動態表達變化及其調控機制，為解析水牛精子發生的分子機制提供理論依據。【方法】從性成熟期水牛睪丸組織中分離生精細胞，采用液相色譜—質譜聯用技術（LC-MS/MS）構建水牛生精細胞蛋白質組和磷酸化蛋白質組表達譜。利用生物信息學軟件對數據進行聯合分析，篩選出關鍵蛋白及其磷酸化修飾位點。構建蛋白互作網絡（PPI）和激酶—底物互作網絡。采用Western blotting檢測和免疫組織化學分析驗證篩選出的關鍵蛋白。【結果】蛋白質組和磷酸化蛋白質組中分別鑒定出1785和129個蛋白，共同鑒定到78個蛋白質。蛋白質組中的高豐度蛋白包括熱休克蛋白和微管蛋白等。GO功能注釋和KEGG信號通路富集分析結果表明，RNA剪接在蛋白質組和磷酸化蛋白質組中均被顯著富集（Plt;0.05），剪接體是磷酸化蛋白富集最多的信號通路。蛋白質組結構域富集分析結果表明，RNA識別基序結構域的富集程度最高。細胞周期依賴性蛋白激酶在蛋白質組PPI和磷酸化蛋白質組激酶—底物互作網絡中均發揮關鍵作用。蛋白質組和磷酸化蛋白質組共同鑒定到的78個蛋白中，有17個蛋白與精子發生密切相關。Western blotting檢測和免疫組織化學分析結果表明，LMNA蛋白及其磷酸化形式LMNA（pS392）均能被成功鑒定，LMNA主要分布于曲細精管基底膜附近，而LMNA（pS392）主要分布于曲細精管管腔附近。【結論】水牛睪丸生精細胞高表達熱休克蛋白可能與水牛適應濕熱氣候有關，高表達微管蛋白可能與生精細胞在曲細精管內的遷移密切相關。RNA剪接和磷酸化修飾可能通過影響蛋白的功能和分布來調控精子的發生過程。

關鍵詞：生精細胞；蛋白質組；磷酸化蛋白質組；RNA剪接

中圖分類號：S823.83文獻標志碼：A文章編號：2095-1191（2025）01-0261-14

Integrated analysis of proteomics and phosphoproteomics in testicular spermatogenic cells of buffalo

HUANG Yu-lin1，WANG Chen-yang1，LIU Xu-yang1，ZHANG Peng-fei2*，FU Qiang3，ZHANG Ming3

（1College of Basic Medicine，Guangxi University of Chinese Medicine，Nanning，Guangxi 530200，China；2College of Life Science and Technology，Guangxi University，Nanning，Guangxi 530004，China；3State Key Laboratory for Conservation and Utilization of Subtropical Agro-bioresources，Guangxi University，Nanning，Guangxi 530004，China）

Abstract：【Objective】To comprehensively analyze the proteomics and immunohistochemistry-based phosphopro-teomics of buffalo testicular spermatogenic cells，and explore the dynamic expression changes of proteins and their regula-tory mechanisms during buffalo spermatogenesis，which could provide theoretical basis for analyzing the molecular mechanism of buffalo spermatogenesis.【Method】Spermatogenic cells were isolated from testicular tissue of sexually ma-ture buffaloes，and the proteome and phosphorylation proteome expression profiles of buffalo spermatogenic cells were constructed by liquid chromatography-mass spectrometry（LC-MS/MS）.The data were jointly analyzed by bioinformatics softwares to screen out key proteins and their phosphorylation modification sites.Protein interaction networks（PPI）and kinase-substrate interaction networks were constructed.Western blotting and immunohistochemical analysis were used to verify the screened key proteins.【Result】A total of 1785 and 129 proteins were identified in the proteome and phospho-proteome respectively，and 78 shared proteins were identified.High-abundance proteins in the proteome included heatshock proteins and tubulin.GO functional annotation and KEGG signaling pathway enrichment analysis results showedthat RNA splicing was significantly enriched in both the proteome and the phosphoproteome（Plt;0.05），and the spliceo-some was the signaling pathway with the most phosphorylated proteins enriched.The results of proteomic domain enrich-ment analysis showed that the RNA recognition motif domain was the most enriched.Cell cycle-dependent protein kinases played a key role in both proteome PPI and kinase-substrate interaction network in phosphoproteome.Among the 78shared proteins identified in the proteome and phosphoproteome，17 proteins were closely related to spermatogenesis.Western blotting and immunohistochemical analysis showed that both LMNA protein and the phosphorylated form LMNA（pS392）could be successfully identified.LMNA was mainly distributed near the basement membrane of theseminiferoustubules，while LMNA（pS392）was mainly distributed near the lumen of theseminiferous tubules.【Conclusion】The high expression of heat shock proteins in buffalo testicular spermatogenic cells may be related to the adaptation of buffalo to hot and humid climate，and the high expression of tubulin may be closely related to the migration of spermatogenic cellsin seminiferous tubules.RNA splicing and phosphorylation modifications may regulate spermatogenesis by affecting pro-tein function and distribution.

Key words：spermatogenic cells；proteome；phosphoproteome；RNA splicing

Foundation items：National Natural Science Foundation of China（32160789）

0引言

【研究意義】精子發生是一個復雜的生物學進程，伴隨時空變化，包括精原細胞有絲分裂、精母細胞減數分裂及精子形成。生精細胞作為精子發生的執行者，涵蓋了多種類型的精原細胞、精母細胞與精子細胞（Tan et al.，2023）。精原細胞為二倍體祖細胞，具有分化為精母細胞和進行有絲分裂自我更新的雙重功能；初級精母細胞經過第一次減數分裂后形成次級精母細胞，次級精母細胞在第二次減數分裂中進一步分裂為單倍體的圓形精子細胞；圓形精子細胞經歷細胞質及細胞核的系列變化，最終形成精子并釋放至生精小管（Rabbaniet al.，2022）。水牛的生精表觀遺傳周期為8.6 d，生精過程持續約38 d，屬于大型家畜中較短的生精周期（Huang et al.，2016）。目前，針對水牛雄性生殖的研究較為匱乏，且其生精過程與其他物種差異較大。因此，深入探究其分子機制具有重要意義。【前人研究進展】精子發生歷來是研究的重點，尤其是不同物種間的精子發生差異。Gan等（2013）利用基于iTRAQ的定量蛋白質組學方法，鑒定了2008種小鼠精原細胞、精母細胞及精子中的蛋白；Chen等（2017）測序分析了男性生殖細胞中的小RNA，并構建了miRNAs圖譜；Chen等（2018）采用高精度單細胞測序方法構建了小鼠精子發生過程的轉錄組圖譜；Wang等（2018）首次在單細胞水平闡明了人類精子發生過程中基因表達的調控網絡；吳朝棟（2019）利用全基因組DNA甲基化測序和NA-SEQ技術探究了圩豬睪丸甲基化和轉錄組的變化，篩選出了相關候選基因和通路；Wang等（2019a）在分子層面探究了哺乳動物減數分裂過程中染色質的三維結構；哈斯高娃（2022）通過對不同年齡段雙峰駝睪丸進行全轉錄組測序，構建了競爭性內源性RNA調控網絡，篩選出與睪丸發育和精子發生相關的候選靶基因；伏曉玉（2023）通過對藏綿羊和胡羊睪丸進行轉錄組測序與生物信息學分析，識別了與精子發生相關的差異表達基因、轉錄本及可變剪接事件；石軍（2023）利用轉錄組測序分析了牦牛睪丸組織不同發育階段的差異表達基因，篩選出了與胰島素和肽類激素相關的信號通路；鄒玥（2023）對黔北麻羊睪丸進行轉錄組分析，發現IGF1基因在睪丸發育過程中發揮重要作用；Huang等（2023）采用單細胞多組學測序揭示了男性精母細胞減數分裂重組過程中DNA去甲基化事件；Sahlu（2023）通過對荷斯坦公牛睪丸進行全轉錄組分析，構建了競爭性內源性RNA調控網絡，并篩選出了與睪丸發育相關的候選基因；馬劍峰等（2024）則通過轉錄組測序分析了大白豬和青峪豬睪丸中差異基因的表達，探討了其生物學功能；梅艷芳等（2024）通過高通量轉錄組測序發現，12月齡布萊克特黑牛精子發生過程中差異表達基因主要富集在Wnt和PI3K-Akt信號通路；石淳元和趙彥玲（2024）通過對藏野豬的睪丸進行轉錄組測序，篩選出了與性成熟相關的候選基因。【本研究切入點】蛋白作為生命活動的關鍵執行者，其翻譯后修飾對信號轉導、蛋白定位和降解等生物學過程具有重要影響，并與多種疾病密切相關（Kitamura and Galligan，2023）。磷酸化是哺乳動物體內最常見的翻譯后修飾，由500種不同的激酶催化絲氨酸、蘇氨酸和酪氨酸發生磷酸化。水牛精子發生的蛋白質組學研究主要集中在睪丸和精子水平（Holland and Ohlendieck，2015；Huang et al.，2015；Hou et al.，2019；Li etal.，2019），鮮見對生精細胞蛋白質組學及磷酸化蛋白質組學的研究報道。【擬解決的關鍵問題】結合蛋白質組學與磷酸化蛋白質組學技術，探究水牛生精細胞中的蛋白表達及翻譯后修飾，為探究精子發生機制和水牛優良種質的培育提供理論參考。

1材料與方法

1.1試驗材料

新鮮睪丸組織取自當地屠宰場性成熟期水牛，取材后，將組織置于無菌生理鹽水中，3 h內送至實驗室。所用抗體均購自美國Abcam公司；ProLongTM Gold抗熒光淬滅封片劑和磷酸化肽段富集試劑盒購自美國ThermoFisher Scientific公司；色譜柱和純化柱購自美國Waters公司。

1.2生精細胞分離與鑒定

將取回的睪丸去除被膜和白膜，使用75%酒精消毒，剪碎，放入含0.1%膠原酶Ⅳ的DMEM/F12細胞培養基中，37℃搖床孵育15min，離心去除上清液。加入含0.25%胰蛋白酶和2μg/mL DNaseⅠ的DMEM/F12細胞培養基，在37℃下消化10 min，待曲細精管完全酶解后，加入等量完全培養基終止消化。用40μm過濾篩過濾細胞，130×g離心5 min。培養基重懸后，將細胞轉移到0.1%明膠覆蓋的100 mm細胞培養皿中，在37℃培養箱中培養4h，收集懸浮細胞進行鑒定。將生精細胞滴于載玻片中央，自然晾干后，使用4%多聚甲醛在室溫下固定10 min；PBS洗滌后，加入含0.1%Triton X-100的PBS，室溫孵育10min。加入一抗DDX4，4℃下孵育12 h；加入二抗IgG-Hamp;L（Alexa Fluor?488），37℃下避光孵育1h；用PBS洗滌載玻片3次，加入含DAPI的ProLongTM Gold抗熒光淬滅封片劑，然后將載玻片置于熒光顯微鏡下采集圖像。

1.3蛋白提取與烷基化修飾

用PBS洗滌生精細胞3次，加入10 mL蛋白裂解緩沖液（8 mol/L尿素、20 mmol/LHEPES、2.5 mmol/L焦磷酸鈉、1 mmol/Lβ-甘油磷酸酯、1 mmol/L釩酸鈉和1%蛋白酶抑制劑，pH 8.0）；使用超聲細胞破碎儀對樣品進行超聲處理，檢測蛋白質濃度；每個樣品取5 mg蛋白，加入1.25 mol/L DTT，55℃溫育30 min；加入樣品體積1/10的碘乙酰胺，室溫避光孵育15 min；加入3倍樣品體積的20 mmol/L HEPES （pH 8.0）稀釋烷基化蛋白溶液；加入樣品體積1/50的1 mg/mL胰蛋白酶，37℃下消化12h。

1.4肽段純化

酶解后的肽段利用Sep-Pak C18純化柱純化。向樣品中加入20%TFA，于冰上孵育15 min；離心，將懸浮液轉移到新的15 mL離心管中；使用5 mL乙腈激活Sep-Pak C18色譜柱，并用1、3、6 mL 0.1%TFA連續清洗；加入樣品，采用1、5、6 mL 0.1%TFA清洗色譜柱；利用2 mL色譜柱清洗液（0.1%TFA，5%乙腈）清洗一次；加入2 mL洗脫液（0.1%TFA，40%乙腈）收集多肽，重復3次；收集洗脫液，于-80℃冰箱冷凍過夜。純化后的肽段5%用于蛋白質組分析，95%用于磷酸化肽段富集和磷酸化蛋白質組分析。

1.5磷酸化肽段富集

利用TiO?磷酸肽段富集試劑盒富集肽段。組裝磷酸化柱，使用20μL洗滌緩沖液清洗，3000×g離心2 min；使用20μL結合/校正緩沖液在3000×g下離心2 min；將樣品加入磷酸化柱中，1000×g離心5 min；分別加入結合/平衡緩沖液、平衡緩沖液、洗滌緩沖液和水，3000×g離心2 min。磷酸化樣品用50μL洗脫緩沖液在1000×g下洗脫5min。

1.6液相色譜分離

采用E2695高pH反相液相色譜法，利用Waters X-Bridge C18色譜柱分離多肽。流動相A為2%的乙腈水溶液，流動相B為98%的乙腈水溶液，pH 10.0，洗脫時間80 min。在0.5 mL/min的流速下，流動相B的濃度在10～55 min內從0增加到35%。收集10～55 min的餾分，最后合并為12個餾分。

1.7質譜分析

使用EASY-nLC 1000納升級液相色譜系統對肽段進行分離。多肽溶解于溶劑A（2%乙腈，0.1%甲酸）中，然后用溶劑B（98%乙腈，0.1%甲酸）洗脫，流速為250 nL/min，洗脫程序：2%～6%溶劑B 5 min，6%～15%溶劑B 27 min，15%～23%溶劑B 17 min，23%～98%溶劑B 2 min，98%～2%溶劑B 5 min和2%溶劑B 5min。對于磷酸化蛋白質組的肽段，使用溶劑B以250 nL/min的流速洗脫，歷時90 min，洗脫程序：0～13%溶劑B 4 min、13%～20%溶劑B 40 min、20%～32%溶劑B 36 min、32%～98%溶劑B 5 min和98%溶劑B 5min。質譜儀參數設置：一級質譜分辨率70000 m/z，一級質譜掃描范圍300～2000 m/z，二級質譜分辨率17500m/z，碰撞模式為高能碰撞解離（HCD），歸一化碰撞能量27。

1.8生物信息學分析

質譜數據分析采用Proteome Discoverer 2.1和UniProt蛋白數據庫（https：//www.uniprot.org/tax-onomy/9913）。檢索條件設定：一級質量誤差：±20 PPM；二級質量誤差：±0.02 Da；靜態修飾：氨基甲酰化（C）乙基；動態修飾：蛋氨酸（M）氧化；酶消化方法：胰蛋白酶；最大缺失位點：2；蛋白和肽段鑒定的假陽性率（FDR）：1%。建立表達譜后，對蛋白質組和磷酸化蛋白質組進行GO功能注釋（https：//david.ncifcrf.gov/tools.jsp）和KEGG信號通路富集（https：//www.genome.jp/kegg）分析。使用STRING 11.0（https：//string-db.org/）獲取蛋白互作網絡（PPI），并將數據導入Cytoscape 3.2.0進行可視化。采用精子發生在線1.0數據庫（https：//mcg.ustc.edu.cn/bsc/spermgenes/simple.php）富集與精子發生相關的蛋白。

1.9免疫組織化學染色

睪丸組織使用4%多聚甲醛固定24h，采用梯度酒精、二甲苯脫水并用石蠟包埋；將樣本切成石蠟切片，放入50℃烘箱中烘烤12h；切片采用梯度酒精脫蠟，并使用3%甲醇過氧化氫孵育30min；采用pH 6.0的檸檬酸鈉緩沖液進行抗原修復，PBST沖洗；1%Triton X-100-PBST溶液孵育30 min；5%山羊血清封閉1 h；加入一抗DDX4、LMNA和LMNA（pS392）于4℃下孵育12 h，PBST洗滌3次；加入二抗IgG-Hamp;L（Alexa Fluor?488）于37℃下避光孵育1 h，PBST洗滌3次；加入含DAPI的ProLongTM Gold抗熒光淬滅封片劑，于熒光顯微鏡下采集圖像。

1.10 Western blotting檢測

蛋白樣品經5%SDS-PAGE（80 V，30 min）和12%SDS-PAGE（120 V，90 min）電泳分析后，在100 V 300 mA條件下轉膜1 h；室溫下使用5%脫脂奶粉溶液封閉1 h，以阻止非特異性結合；1×TBST洗滌3次，加入一抗LMNA和LMNA（pS392），于4℃下孵育12h；1×TBST洗滌3次，加入二抗，室溫下孵育1h；加入ECL發光顯影液，在凝膠顯影儀下曝光觀察。

2結果與分析

2.1 DDX4定位與生精細胞鑒定結果

DDX4是一種生殖細胞特異性蛋白，僅在睪丸生精細胞中表達（Kimetal.，2015）。水牛睪丸免疫組織化學染色結果（圖1-A）顯示，DDX4分布于曲細精管，定位于細胞質，與生精細胞在曲細精管中的位置相符。對分離的生精細胞進行免疫熒光染色鑒定，結果（圖1-B）顯示，DDX4分布于細胞質，經統計分析，生精細胞的陽性率為95%。

2.2生精細胞蛋白質組理化性質預測結果

在生精細胞蛋白質組中，共鑒定出11809個特異性肽段，對應1785個蛋白。蛋白理化性質預測結果顯示，蛋白理論等電點（pI）主要分布在5～10（圖2-A），蛋白分子量主要分布在20～80 Da（圖2-B），蛋白序列覆蓋率主要分布在40%以內（圖2-C）。所有鑒定出的蛋白至少有一個特異性肽段（圖2-D）。

2.3高豐度蛋白功能分析結果

肽匹配圖譜（PSM）可作為蛋白的半定量指標，PSM值越大表示蛋白豐度越高。在生精細胞蛋白質組中，發現了10種高豐度蛋白（表1），其中HSPA2、HSP90AA1、HSP90B1和HSPA5為熱休克蛋白，具有分子伴侶蛋白活性，能增強細胞抗應激能力，這可能與水牛適應濕熱氣候有關。微管蛋白家族中TUBB4B、LOC100295712、TUBB5和TUBA1A組成了微管，微管作為細胞骨架的主要組成部分，在維持細胞內物質的形狀、運動和運輸方面發揮著重要作用，可能與生精細胞在曲細精管內的遷移密切相關。血紅蛋白亞基β（HBB）和延伸因子1-α1（EEF1A1）負責調節抗病毒先天免疫反應（Yang et al.，2019），并在mRNA翻譯過程中促進多肽鏈的延伸。

2.4生精細胞蛋白質組GO功能注釋分析結果

GO功能注釋分析按細胞組分（Cellular compo-nent，CC）、生物學過程（Biological process，BP）和分子功能（Molecular function，MF）三大類別進行。生精細胞蛋白質組中有1740個蛋白具有功能信息。在生物學過程中，顯著富集的GO功能條目包括RNA剪接（RNA splicing）、翻譯（Translation）、細胞定位（Cellular localization）、氧化還原過程（Oxidation-reduction process）、蛋白折疊（Protein folding）、細胞酰胺代謝過程（Cellular amide metabolic process）等（Plt;0.05，下同）（圖3-A）。在細胞組分中，細胞器（Organelle）的富集程度最高（圖3-B），可能是由于細胞活性增加導致生精細胞不斷分化所致。在分子功能中，顯著富集的GO功能條目包括分子結合和水解酶活性（Hydrolase activity）等，分子結合中RNA結合（RNA binding）的富集程度最高（圖3-C）。此外，在生物學過程中，還發現了精卵識別、有性生殖、精子發生、雄配子生成、精母細胞分化、生殖細胞發育、配子生成、受精、生殖和精子運動等許多與生殖相關的過程（圖4），這些過程涉及大量功能蛋白和調節因子。

2.5生精細胞蛋白質組結構域富集分析結果

使用DAVID V6.8軟件對水牛生精細胞蛋白質組進行結構域分析，結果（圖5）僅展示前10個被顯著富集的結構域，其中，RNA識別基序結構域（RNA recognition motif domain）富集程度最高，其能識別多種RNA和蛋白，參與多種生物功能，通常與轉錄后基因調控相關，表明精子發生過程中調控機制的復雜性。

2.6生精細胞蛋白質組PPI分析結果

將蛋白質組數據導入STRING數據庫，再將得到的PPI數據導入Cytoscape 3.2.0進行可視化網絡分析。結果顯示，PPI由1094個節點和13763條邊組成，其中多數節點存在2個以上的互作關系（圖6）。選取節點間中心度最大且大于平均節點度的前10個節點蛋白為關鍵蛋白（表2），其中，熱休克蛋白A8（HSPA8）、細胞周期蛋白依賴性激酶1（CDK1）和DNA定向RNA聚合酶II亞基C（POLR2C）是PPI分析中最關鍵的3個蛋白。

2.7生精細胞磷酸化蛋白質組GO功能注釋分析結果

在磷酸化蛋白質組中，共鑒定出129種磷酸化蛋白、201個磷酸化肽段和417個磷酸化位點。其中，絲氨酸（Ser）、蘇氨酸（Thr）和酪氨酸（Tyr）位點的數量分別為328、79和10個。在生精細胞磷酸化蛋白質組中，有123個蛋白質具有注釋信息。在生物學過程中，RNA剪接是最顯著富集的GO功能條目，其他富集的GO功能條目還包括基因表達（Gene expression）、細胞周期（Cell cycle）、有絲分裂運動（Mitotic cytokinesis）、氮化合物過程（Nitrogen com-pound metabolic process）、繁殖（Reproduction）和蛋白質靶向（Protein targeting）等（圖7-A），顯示了磷酸化蛋白功能的多樣性。在細胞組分中，細胞核（Nucleus）為最顯著富集的GO功能條目，表明鑒定出的磷酸化蛋白主要位于細胞核中（圖7-B）。在分子功能中，顯著富集到的GO功能條目有分子結合和酶活性（圖7-C）相關條目。

2.8生精細胞磷酸化蛋白質組KEGG信號通路富集分析結果

在鑒定到的129個磷酸化蛋白中，有62個被成功映射到不同的信號通路上。剪接體（Spliceosome）是磷酸化蛋白富集最多的信號通路，與GO功能注釋分析結果相符，表明這些蛋白在RNA剪接過程中發揮了重要作用。其他富集的信號通路包括代謝途徑（Metabolic pathway）、RNA轉運（RNA transport）、次生代謝產物的生物合成（Biosynthesis of secondarymetabolite）、細胞周期（Cell cycle）和細胞衰老（Cel-lular senescence）等（圖8）。

2.9生精細胞磷酸化蛋白質組激酶—底物互作網絡分析結果

使用iGPS構建生精細胞激酶及其底物的磷酸化位點網絡，根據已確定的磷酸化位點，預測激酶與底物間的互作關系。結果（圖9）顯示，共有451個激酶—底物互作關系，其中，與359個磷酸化位點對應的36個底物受205個預測激酶的調控。細胞周期蛋白依賴性激酶7（CDK7）調控的底物數量最多，其活性最高。

2.10 GO功能注釋聯合分析結果

從蛋白質組和磷酸化蛋白質組的GO功能注釋結果中，選取前20個顯著富集的生物學過程進行聯合分析。結果（圖10）顯示，2種組學的共同信號通路均與RNA剪接相關。此外，蛋白質組中單獨富集的生物學過程包括蛋白折疊、氧化還原過程、細胞酰胺代謝過程等；磷酸化蛋白質組中單獨富集的生物學過程為細胞器組織（Organelle organization）、染色體組織（Chromosome organization）、細胞衰老的正調控（Positive regulation of cell aging）等。在磷酸化蛋白質組中，參與細胞衰老的蛋白包括高遷移率組A1（HMGA1）、Sirtuin 1（SIRT1）和LaminA/C（LMNA）。

2.11 PPI聯合分析結果

蛋白質組和磷酸化蛋白質組共同鑒定到78個蛋白，分別占蛋白質組的4.4%，占磷酸化蛋白質組的60.5%。將數據導入STRING數據庫，再將得到的PPI數據導入Cytoscape 3.2.0進行可視化網絡分析。結果顯示（圖11），該網絡由56個節點和115條邊組成。在該網絡中發現了細胞周期蛋白依賴性激酶1（CDK1）、真核翻譯起始因子4 Gamma 1（EIF4G1）、小核核糖核蛋白U5亞基200（SNRNP200）、腫瘤蛋白p53結合蛋白1（TP53BP1）和通用轉錄因子IIF亞基1（TFIIF）5個關鍵蛋白（表3）。

2.12精子發生相關蛋白富集結果

利用精子發生在線1.0數據庫檢索蛋白質組和磷酸化蛋白質組共同鑒定到的78個蛋白，結果（表4）顯示，有17個蛋白在數據庫中被檢索到。

2.13數據驗證

為驗證蛋白質組學和磷酸化蛋白質組學數據的可靠性，選擇LMNA蛋白進行分析。Western blotting檢測結果（圖12-A）顯示，LMNA蛋白及其磷酸化形式LMNA（pS392）均被成功鑒定。免疫組織化學染色結果（圖12-B）顯示，LMNA主要分布于曲細精管基底膜附近，而LMNA（pS392）主要分布于曲細精管管腔附近。同一蛋白在不同亞細胞位置的分布差異，間接表明磷酸化修飾可能影響蛋白的功能和分布。

3討論

精子發生是一個復雜且嚴格調控的過程，涉及數千種蛋白的協調表達，以確保雄性配子的持續產生，而生精細胞是精子發生過程的主要參與者。本研究對生精細胞中的1785個蛋白和129個磷酸化蛋白進行了鑒定和生物信息學分析。在生精細胞蛋白質組中，發現了10種高豐度蛋白，有4個蛋白屬于熱休克蛋白家族（HSPA2、HSP90AA1、HSP90B1和HSPA5），熱休克蛋白作為免疫系統的調節因子，在免疫抑制和免疫刺激中發揮重要作用（Hristova，2012），提示精子發生可能通過大量熱休克蛋白對各種刺激做出反應，以維持精子發生內部環境的平衡。此外，4種高豐度蛋白屬于微管蛋白家族（TUBB4B、LOC100295712、TUBB5和TUBA1A），微管蛋白是微管的組成部分，微管是細胞分裂和遷移等細胞過程的基礎。微管動力學參與維持細胞形態和功能，并在生殖細胞有絲分裂和減數分裂中發揮重要作用（O’Donnell and O’Bryan，2014），這與精子發生需要持續的細胞分化相契合。

RNA剪接是蛋白質組和磷酸化蛋白質組中富集最多的生物學過程。Zhang等（2019）研究發現，雙酚A（BPA）對RNA剪接的阻斷會導致后代雄性小鼠睪丸發育不良；Song等（2020）研究發現，睪丸發育受可變剪接的調控水平高于其他組織；Wu等（2022）研究表明，RNA剪接的調控是男性生育的必要條件，BUD31介導的精子干細胞內的可變剪接可以調節哺乳動物雄性生殖細胞的自我更新和分化。提示RNA剪接在生精細胞中的存在是維持精子發生正常進行的必要條件。此外，在高豐度蛋白中發現了大量熱休克蛋白家族成員，而在生精細胞蛋白質組PPI分析中同樣發現了關鍵蛋白HSPA8。Wang等（2019b）研究發現，急性熱應激會上調雞睪丸中HSPA8蛋白的表達，提示HSPA8可在外界刺激下維持正常的精子發生。生精細胞磷酸化蛋白質組KEGG信號通路富集分析結果顯示，剪接體是主要的調控途徑。Wu等（2016）研究表明，剪接體缺陷會影響人類精原細胞的分化，并導致非梗阻性無精子癥；Alikhani等（2017）研究發現，生精細胞匱乏的人類睪丸中參與剪接體的蛋白質表達顯著下調。Kim等（2001）研究發現，CDK7存在于減數分裂過程中，并在精母細胞中高度表達。本研究發現，CDK7調控的底物數量最多，活性最高，且CDK1處于PPI網絡中心，表明細胞周期依賴激酶在水牛生精細胞的分化過程中高度活躍。Liu等（2017）研究表明，SIRT1基因敲除小鼠對精子發生過程中生精細胞的早期有絲分裂和減數分裂無明顯影響，但在頂體發生過程存在缺陷，導致圓形精子不育的表型；Wang等（2018）對人類睪丸細胞的研究發現，HMGA1在精原細胞基底層附近活躍，在其他生精細胞中逐漸下調。本研究發現，細胞衰老的正調控僅在磷酸化蛋白質組中被富集，涉及的蛋白有HMGA1、SIRT1和LMNA，與Liu等（2017）、Wang等（2018）的研究結果相符。Alsheimer等（2004）研究發現，LMNA基因敲除的小鼠在精子發生過程中出現障礙，具體表現為精母細胞中性染色體配對存在嚴重缺陷，導致青春期小鼠精母細胞大量凋亡；Shen等（2014）通過siRNA干擾技術降低了LMNA的表達，能增加小鼠精子頂體分裂，并顯著增加精子頭部和尾部的畸形率；Milbradt等（2016）研究表明，LMNA的第22和392位絲氨酸位點的磷酸化是有絲分裂過程中核纖層短暫解體的關鍵因素。在本研究中，LMNA及其磷酸化形式LMNA（pS392）在細胞中均有表達，表明磷酸化可能會改變蛋白在精子發生中的作用。

4結論

水牛睪丸生精細胞高表達熱休克蛋白可能與水牛適應濕熱氣候有關，高表達微管蛋白可能與生精細胞在曲細精管內的遷移密切相關。RNA剪接和磷酸化修飾可能通過影響蛋白的功能和分布來調控精子的發生過程。

參考文獻（References）：

伏曉玉.2023.藏羊和湖羊睪丸的組織學研究及與精子發生相關基因的鑒定［D］.蘭州：甘肅農業大學［Fu X Y.2023.Histological study of testes and identification of genes associa-ted with spermatogenesis in Tibetan sheep and Hu sheep［D］.Lanzhou：Gansu Agricultural Univer-sity.］doi：10.27025/d.cnki.ggsnu.2023.000228.

哈斯高娃.2022.基于全轉錄組測序篩選雙峰駝睪丸發育及精子發生相關基因［D］.呼和浩特：內蒙古農業大學.［Haskova.2022.Whole transcriptome sequencing forscreening of genes related to testicular development andspermatogene-sis in Bactrian camels［D］.Hohhot：Inner Mongolia Agricultural University.］doi：10.27229/d.cnki.gnmnu.2022.001368.

馬劍峰，陳思宇，師愉倩，沈林園，朱礪，甘麥鄰.2024.轉錄組測序解析瘦肉型和脂肪型公豬睪丸的基因表達模式差異［J］.中國畜牧雜志，60（8）：191-198.［Ma J F，Chen S Y，Shi Y Q，Shen L Y，Zhu L，Gan M L.2024.Transcrip-tome sequencing resolves differences in gene expression patterns in the testes of lean and fat boars［J］.Chinese Jour-nal of Animal Science，60（8）：191-198.］doi：10.19556/j.0258-7033.20240430-05.

梅艷芳，劉瑞莉，韓明軒，張文琦，陳振鵬，王秀源，苗秀萍，肖超柱，張任政，董雅娟，柏學進.2024.基于轉錄組測序技術挖掘影響布萊凱特黑牛睪丸大小的功能基因［J］.華北農學報，39（1）：191-200.［Mei Y F，Liu R L，Han M X，Zhang W Q，Chen Z P，Wang X Y，Miao X P，Xiao C Z，Zhang R Z，Dong Y J，Bai X J.2024.Analysis of reproduc-tion function gene of testicle in different size Black cattle based on transcriptome sequencing［J］.Acta Agriculturae Boreali-Sinica，39（1）：191-200.］doi：10.7668/hbnxb.2019 4430.

石淳元，趙彥玲.2024.基于轉錄組測序的藏豬睪丸組織性成熟相關基因鑒定［J］.南方農業學報，55（3）：680-688.［Shi C Y，Zhao Y L.2024.Identification of genes relatedto sexual maturation in Tibetan pig testis tissue based ontranscriptome sequencing［J］.Journal of Southern Agricul-ture，55（3）：680-688.］doi：10.3969/j.issn.2095-1191.2024.03.009.

石軍.2023.基于轉錄組分析不同年齡牦牛睪丸差異基因及GRB10和GRB14基因功能探究［D］.蘭州：甘肅農業大學.［Shi J.2023.Transcriptomics-based analysis of the diffe-rential expressed genes and functional exploration of GRB10 and GRB14 genes in yaks testises with different ages［D］.Lanzhou：Gansu Agricultural University.］doi：10.27025/d.cnki.ggsnu.2023.000177.

吳朝棟.2019.圩豬初情期前后睪丸組織WGBS和RNA-seq組學分析［D］.合肥：安徽農業大學.［Wu C D.2019.Studies on the genome-wide DNA methylation changes and transcriptomic profiling of testis in Weizhu boars［D］.Hefei：Anhui Agricultural University.］doi：10.26919/d.cnki.gannu.2019.000793.

鄒玥.2023.基于轉錄組學的黔北麻羊睪丸發育相關基因篩選及表達規律研究［D］.貴陽：貴州大學.［Zou Y.2023.Transcriptomics based screening and expression of genes related to testicular development in Qianbei Ma goat［D］.Guiyang：Guizhou University.］doi：10.27047/d.cnki.ggudu.2023.000725.

Alikhani M，Mirzaei M，Sabbaghian M，Parsamatin P，Karamzadeh R，Adib S，Sodeifi N，Gilani M A S，Zabet-Moghaddam M，Parker L，Wu Y Q，Gupta V，Haynes PA，Gourabi H，Baharvand H，Salekdeh G H，Salekdeh G H.2017.Quantitative proteomic analysis of human testis reveals system-wide molecular and cellular pathways asso-ciated with non-obstructive azoospermia［J］.Journal of Proteomics，162：141-54.doi：10.1016/j.jprot.2017.02.007.

Alsheimer M，Liebe B，Sewell L，Stewart CL，Scherthan H，Benavente R.2004.Disruption of spermatogenesis in micelacking A-type lamins［J］.Journal of Cell Science，117（Pt 7）：1173-1178.doi：10.1242/jcs.00975.

Chen X X，Che D X，Zhang P F，Li X L，Yuan Q Q，Liu T T，Guo J Y，Feng T Y，Wu L G，Liao M Z，He Z P，Zeng W X.2017.Profiling of miRNAs in porcine germ cells during spermatogenesis［J］.Reproduction，154（6）：789-798.doi：10.1530/REP-17-0441.

Chen Y，Zheng Y X，Gao Y，Lin Z，Yang S M，Wang T T，Wang Q，Xie N N，Hua R，Liu M X，Sha J H，Griswold M D，Li J S，Tang F C，Tong M H.2018.Single-cell RNA-seq unco-vers dynamic processes and critical regulators in mouse spermatogenesis［J］.Cell Research，28（9）：879-896.doi：10.1038/s41422-018-0074-y.

Gan H Y，Cai T X，Lin Y W，Wu Y J，Wang X X，Yang F Q，Han C S.2013.Integrative proteomic and transcriptomic analyses reveal multiple post-transcriptional regulatory mechanisms of mouse spermatogenesis［J］.Molecularamp;Cellular Proteomics：MCP，12（5）：1144-57.doi：10.1074/mcp.M112.020123.

Holland A，Ohlendieck K.2015.Comparative profiling of the sperm proteome［J］.Proteomics，15（4）：632-648.doi：10.1002/pmic.201400032.

Hou Z，Fu Q，Huang Y L，Zhang P F，Chen F M，Li M X，Xu Z Z，Yao S，Chen D R，Zhang M.2019.Comparative pro-teomic identification buffalo spermatozoa during in vitro capacitation［J］.Theriogenology，126：303-09.doi：10.1016/j.theriogenology.2018.12.025.

Hristova I.2012.Role of heat shock proteins（Hsp）in human and mammalian fertilization and pregnancy.Part II［J］.Akusherstvo I Ginekologiia，51（6）：37-40.

Huang Y L，Fu Q，Pan H，Chen F M，Zhao X L，Wang H J，Zhang P F，Huang F L，Lu Y Q，Zhang M.2016.Spermatogenesis-associated proteins at different develop-mental stages of buffalo testicular seminiferous tubules identified by comparative proteomic analysis［J］.Proteo-mics，16（14）：2005-2018.doi：10.1002/pmic.201500547.

Huang Y L，Fu Q，Yang L，Guan J L，Pan H，Chen F M，Lu K L，Zhang M.2015.Differences between high-and low-motility buffalo sperm identified by comparative proteomics［J］.Reproduction in Domestic Animals，50（3）：443-451.doi：10.1111/rda.12511.

Huang Y P，Li L，An G，Yang X Y，Cui M M，Song X L，Lin J，Zhang X L，Yao Z K，Wan C，Zhou C，Zhao J X，Song K，Ren S F，Xia X Y，Fu X，Lan Y，Hu X S，Wang W，Wang M，Zheng Y，Miao K，Bai X C，Hutchins A P，Chang G，

Gao S，Zhao X Y.2023.Single-cell multi-omics sequen-cing of human spermatogenesis reveals a DNA demethyla-tion event associated with male meiotic recombination［J］.Nature Cell Biology，25（10）：1520-1534.doi：10.1038/s41556-023-01232-7.

Kim J M，McGaughy J T，Bogle R K，Ravnik S E.2001.Meio-tic expression of the cyclin H/Cdk7 complex in male germ cells of the mouse［J］.Biology of Reproduction，64（5）：1400-1408.doi：10.1095/biolreprod64.5.1400.

Kim J Y，Jung H J，Yoon M J.2015.VASA（DDX4）is a puta-tive marker for spermatogonia，spermatocytes and round spermatids in stallions［J］.Reproduction in Domestic Ani-mals，50（6）：1032-1038.doi：10.1111/rda.12632.

Kitamura N，Galligan J J.2023.A global view of the human post-translational modification landscape［J］.Biochemical Journal，480（16）：1241-1265.doi：10.1042/BCJ20220251.

Li Y，Cheng Y W，Zhu T Y，Zhang H，Li W，Guo Y S，Qi Y L，Chen X，Zhang J，Sha J H，Zhou Z M，Zhu H，Guo X J.2019.The protein phosphorylation landscape of mouse spermatids during spermiogenesis［J］.Proteomics，19（11）：e1900055.doi：10.1002/pmic.201900055.

Liu C，Song Z H，Wang LN，Yu H Y，Liu W X，Shang Y L，Xu Z L，Zhao H C，Gao F Y，Wen J M，Zhao L N，Gui Y T，Jiao J W，Gao F，Li W.2017.Sirt1 regulates acrosome bio-genesis by modulating autophagic flux during spermioge-nesis in mice［J］.Development，144（3）：441-451.doi：10.1242/dev.147074.

Milbradt J，Hutterer C，Bahsi H，Wagner S，Sonntag E，Horn A H C，Kaufer B B，Mori Y，Sticht H，Fossen T，Marschall M.2016.The prolyl isomerase Pin1 promotes the herpesvirus-induced phosphorylation-dependent disassem-bly of the nuclear lamina required for nucleocytoplasmic egress［J］.PLoS Pathog，12（8）：e1005825.doi：10.1371/journal.ppat.1005825.

O’Donnell，L，O’Bryan M K.2014.Microtubules and spermato-genesis［J］.Seminars in Cellamp;Developmental Biology，30：45-54.doi：10.1016/j.semcdb.2014.01.003.

Rabbani M，Zheng X N，Manske G L，Vargo A，Shami A N，Li J Z，Hammoud S S.2022.Decoding the spermatogenesis program：New insights from transcriptomic analyses［J］.Annual Review of Genetics，56：339-368.doi：10.1146/annurev-genet-080320-040045.

Sahlu W B.2023.Whole transcriptome sequen-cing for screen-ing of genes related to testicular development and sper-matogenesis in Holstein bull testes［D］.Beijing：Chinese Academy of Agricultural Sciences.doi：10.27630/d.cnki.gznky.2023.000899.

Shen J，Chen W，Shao B B，Qi Y J，Xia Z R，Wang F Q，Wang L，Guo X J，Huang X Y，Sha J H.2014.Lamin A/C pro-teins in the spermatid acroplaxome are essential in mouse spermiogenesis［J］.Reproduction，148（5）：479-487.doi：10.1530/REP-14-0012.

Song H，Wang L，Chen D，Li F.2020.The function of pre-mrna alternative splicing in mammal spermatogenesis［J］.Inter-national Journal of Biological Sciences，16（1）：38-48.doi：10.7150/ijbs.34422.2.

Tan X S，Zheng C H，Zhuang Y H，Jin P P，Wang F C.2023.The m6A reader PRRC2A is essential for meiosis I comple-tion during spermatogenesis［J］.Nature Communications，14（1）：1636.doi：10.1038/s41467-023-37252-y.

Wang M，Liu X X，Chang G，Chen Y D，An G，Yan L Y，Gao S，Xu Y W，Cui Y L，Dong J，Chen Y H，Fan XY，HuY Q，Song K，Zhu X H，Gao Y，Yao Z K，Bian S H，Hou Y，Lu J H，Wang R，Fan Y，Lian Y，Tang W H，Wang Y P，Liu J Q，Zhao LM，Wang LY，Liu Z T，Yuan R P，Shi Y J，Hu B Q，Ren X L，Tang F C，Zhao X Y，Qiao J.2018.Single-cellrna sequencing analysis reveals sequential cell fate transi-tion during human spermatogenesis［J］.Cell Stem Cell，23（4）：599-614.doi：10.1016/j.stem.2018.08.007.

Wang S H，Cheng C Y，Chen C J，Chan H L，Chen H H，Tang P C，Chen C F，Lee Y P，Huang S Y.2019a.Acute heat stress changes protein expression in the testes of a broiler-typestrain of taiwan country chickens［J］.Animal Biotech-nology，30（2）：129-145.doi：10.1080/10495398.2018.1446 972.

Wang Y，Wang H B，Zhang Y，Du Z H，Si W，Fan S X，Qin D D，Wang M，Duan Y C，Li L F，Jiao Y Y，Li Y Y，Wang Q

J，Shi Q H，Wu X，Xie W.2019b.Reprogramming of meio-tic chromatin architecture during spermatogenesis［J］.Mo-lecular Cell，73（3）：547-561.doi：10.1016/j.molcel.2018.11.019.

Wu D，Khan FA，Huo L J，Sun F，Huang C J.2022.Alternative splicing and MicroRNA：Epigenetic mystique in male reproduction［J］.RNA Biology，19（1）：162-175.doi：10.1080/15476286.2021.2024033.

Wu H，Sun L W，Wen Y，Liu Y J，Yu J，Mao F Y，Wang Y，Tong C，Guo X J，Hu Z B，Sha J H，Liu M X，Xia L X.2016.Major spliceosome defects cause male infertility and are associated with nonobstructive azoospermia in humans［J］.Proceedings of the National Academy of Sciences ofthe United States of America，113（15）：4134-4139.doi：10.1073/pnas.1513682113.

Yang Q，Bai S Y，Li L F，Li S，Zhang Y X，Munir M，Qiu H J.2019.Human hemoglobin subunit beta functions as a pleio-tropic regulator of RIG-I/MDA5-mediated antiviral innate immune responses［J］.Journal of Virology，93（16）：00718-00719.doi：10.1128/JVI.00718-19.

Zhang S L，Bao J L，Gong X C，Shi W Y，Zhong X H.2019.Hazards of bisphenol A—Blocks RNA splicing leading to abnormal testicular development in offspring male mice［J］.Chemosphere，230：432-39.doi：10.1016/j.chemo-sphere.2019.05.044.

（責任編輯：蘭宗寶，張博）