基于轉錄組測序的藏豬睪丸組織性成熟相關基因鑒定

摘要：【目的】基于轉錄組測序分析篩選出藏豬睪丸組織中的性成熟相關基因，為提高藏豬的繁殖育種性能提供 參考依據(jù)。【方法】以性成熟前（2月齡）和性成熟后（24月齡）的藏公豬為研究對象，通過人工閹割法采集藏豬睪丸組織 進行轉錄組測序，依據(jù)錯誤發(fā)現(xiàn)率（FDR）lt;0.05且logFold Change|gt;1的標準篩選差異表達基因（DEGs），綜合GO功 能注釋分析和KEGG信號通路富集分析篩選出藏公豬性成熟相關候選基因，并通過實時熒光定量PCR驗證轉錄組測序結果的準確性?！窘Y果】經轉錄組測序，從性成熟前后藏豬睪丸組織樣品中篩選出5929個DEGs，其中有3822個DEGs呈上調表達、2107個DEGs呈下調表達。篩選獲得的5929個DEGs主要注釋到生精過程、纖毛運動、鞭毛精子運動、運動纖毛、精子主塊、精子頂囊、動力蛋白輕鏈結合、動態(tài)蛋白中間鏈結合等與精子發(fā)生和精子運動等有關的GO功能條目；KEGG信號通路富集分析發(fā)現(xiàn)，DEGs主要富集于受體相互作用通路、軸突引導通路、神經活性配體一受體相互作用通路、鈣信號通路等；綜合DEGs的GO功能注釋分析和KEGG信號通路富集分析結果，共篩選出9個與藏公豬性成 熟相關的基因，分別是BAX、DMRTI、FSHR、IGFIR、IZUMO1、LYZL4、LYZL6、SYCP3和TNP2。實時熒光定量PCR檢測 結果顯示，BAX、DMRTI、FSHR、IGFIR等4個基因在性成熟后藏豬睪丸組織中的相對表達量較性成熟前呈明顯 的下降趨勢，且下降幅度均在50%以上；IZUMO1、LYZL4、LYZL6、SYCP3、TNP2等5個基因的相對表達量較性成熟前 上升了1.07～30.25倍，其中IZUMO1、LYZL4和LYZL6基因的上升倍數(shù)達極顯著水平（Plt;0.01），與轉錄組測序結果基本 一致?！窘Y論】基于轉錄組測序分析從藏豬睪丸組織中篩選獲得9個與性成熟相關的DEGs，即BAX、DMRT1、FSHR、 IGFIR、IZUMOI、LYZL4、LYZL6、SYCP3和TNP2，可作為研究藏豬繁殖育種能力的主要候選基因。

關鍵詞：藏豬；睪丸組織；性成熟；精子；轉錄組測序

文章編號：2095-1191（2024）03-0680-09

中圖分類號：S828.89

文獻標志碼：A

Identification of genes related to sexual maturation in Tibetan pig testis tissue based on transcriptome sequencing

SHI Chun-yuan， ZHAO Yan-ling

（College of Animal Science， Xizang Agricultural and Animal Husbandry University， Nyingchi， Xizang 860000， China）

Abstract： [Objective]Based on transcriptome sequencing analysis， the genes related to sexual maturation were screened out in testicular tissue of Tibetan pigs so as to provide reference for improving the breeding performance of Tibetan pigs. 【Method】Tibetan boars before sexual maturity （2 months old） and after sexual maturity （24 months old） were selected as the study objects， testicular tissues of the Tibetan boars were collected by artificial castration method for tran- scriptomic sequencing， and differentially expressed genes （DEGs） were screened according to the criteria of 1 dis- covery rate （FDR）lt;0.05 and |log2 Fold Changegt;1. GO functional annotation analysis and KEGG signaling pathway enrich- ment analysis were combined to select the candidate genes related to sexual maturity in Tibetan boars. And the accuracy of transcriptome sequencing results were verified by real-time fluorescence quantitative PCR. 【Result]According to tran- scriptome sequencing， 5929 DEGs were screened from testicular tissue samples of Tibetan pigs before and after sexual maturity， among which 3822 DEGs were up-regulated and 2107 DEGs were down-regulated. The 5929 DEGs were se-lected to include spermatogenesis， ciliary motility， flagellate sperm motility， motile cilia， sperm mass， sperm apical sac， kinetic protein light chain binding， dynamic protein intermediate chain binding and other GO functional entries related to spermatogenesis and sperm motility. KEGG signaling pathway enrichment analysis showed that the DEGs was mainly enriched in receptor interaction pathway， axon-guided pathway， neural active ligand-receptor interaction pathway， calcium signaling pathway， etc. Based on the GO functional annotation analysis of DEGs and the enrichment analysis of KEGG signaling pathway， a total of nine genes related to sexual maturation in Tibetan boars were selected， which were BAX， DMRT1， FSHR， IGF1R， IZUMO1， LYZL4， LYZL6， SYCP3 and TNP2. Real-time fluorescence quantitative PCR detection results showed that the relative expression levels of the four genes， including BAX， DMRT1， FSHR and IGF1R， were decreased after sexual maturity compared with that before sexual maturity， all decreased by more than 50%.The relative expression levels of the five genes， including IZUMO1， LYZL4， LYZL6， SYCP3 and TNP2， were increased by 1.07-30.25 times compared with that before sexual maturity， and the increasing times of IZUMO1，LYZL4 and LYZL6 genes reached extremely significant level （Plt;0.01）， which was basically consistent with the transcriptome sequencing results. 【Conclusion】Nine DEGs related to sexual maturity are detected in testicular tissue of Tibetan pigs before and after sexual maturity， namely BAX， DMRT1， FSHR， IGF1R， IZUMO1， LYZL4， LYZL6， SYCP3 and TNP2， which can be used as the main candidate genes to study the breeding ability of Tibetan pigs.

Key words： Tibetan pigs; testicular tissue; sexual maturity; sperm; transcriptome sequencing

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

0 引言

【研究意義】藏豬作為我國古老珍稀的瘦肉型豬 種，主要分布在青藏高原半農半牧區(qū)，具有四肢堅 實、脂肪沉積能力強及肉質風味鮮美等特點，但也存 在生長速度緩慢、育肥增重難、屠宰率低和繁殖能力 不強等缺陷（楊燁城等，2022），且多數(shù)藏豬為近親繁 殖，嚴重制約其規(guī)?；B(yǎng)殖的發(fā)展速度。篩選鑒定 藏公豬性成熟相關基因，有利于對藏豬遺傳特性進 行深度了解，如胰島素樣生長因子1（IGF1）和胰島素樣生長因子1受體（IGF1R）是影響動物生長發(fā)育 的重要候選基因。因此，加強藏豬遺傳特性研究 及繁殖相關基因鑒定，對開展其種質資源保護及品 種改良具有重要意義?！厩叭搜芯窟M展】至今，轉錄 組學在豬業(yè)科學上的研究已得到高度重視和廣泛 應用（易凡利等，2018；馬圍圍等，2021；韋明邦等， 2022）。如基于轉錄組測序分析篩選出豬遺傳候選 基因（甘麥鄰等，2020；張福平等，2021）、豬脂肪沉積 關鍵基因（孫菁希等，2021；翟麗維等，2022）、豬黑色 素生成相關基因（覃海等，2022）、豬骨骼肌發(fā)育相關 基因（熊和麗等，2023）等，在豬相關病毒性疾病研究中也有應用（梁婉等，2019；王璐夢等，2021）。這些 研究通過轉錄組測序獲得豬種在某些性狀方面的特 殊基因，后續(xù)通過對這些特殊基因進行深入挖掘，可 為生產性能研究提供重要的理論依據(jù)。在藏豬上， 聶靖茹等（2022）基于RNA-Seq技術對迪慶藏豬背最長肌組織進行轉錄組測序，結果篩選出ATF3、CCL2、CCN1、EGR1、FOXO1等12個影響大型迪慶藏豬不同生長階段肌肉生長的差異表達基因（Dif-ferentially expressed genes，DEGs）;宋明坤等（2022） 利用背最長肌轉錄組數(shù)據(jù)篩選藏公豬與大約克夏公豬差異表達IncRNAs和DEGs，發(fā)現(xiàn)差異表達 lncRNAs可能通過調控HESI、LEP、TGFB2和PER2 等基因的表達而影響豬肉品質；董新星等（2023）、王 琳等（2023）通過轉錄組測序篩選并鑒定出不同生長 階段迪慶藏豬背、腹脂代謝的DEGs?！颈狙芯壳腥?點】藏豬作為我國珍貴的物種遺傳資源，在生產性 狀、繁殖性能及相關基因功能等方面已得到深入研 究（易凡利等，2018；董世雄等，2022；劉思源等， 2022；韓雨晴等，2023），但針對藏豬性成熟前后睪丸 組織轉錄組測序的研究鮮見報道。【擬解決的關鍵問 題】以性成熟前（2月齡）和性成熟后（24月齡）的藏 公豬為研究對象，基于轉錄組測序分析篩選出藏豬 睪丸組織中的性成熟相關基因，為提高藏豬的繁殖 育種性能提供參考依據(jù)。

1材料與方法

1.1試驗材料

以飼養(yǎng)于西藏林芝（海拔3000m）牧場且營養(yǎng) 良好、體況相似、繁殖機能正常的藏豬為研究對象， 隨機選取性成熟前（2月齡）和性成熟后（24月齡）的 藏公豬各3頭，按照《實驗動物管理條例》規(guī)定，通過人工閹割法采集藏豬睪丸組織，將表皮劃開取黃豆 粒大小樣本，生理鹽水沖洗后放入裝有RNA保存液的凍存管中，編號[性成熟前（TPA）：TPO1、TP02、TP03；性成熟后（TPB）：TP04、TP05、TP06]，然后置于-80℃冰箱保存?zhèn)溆?。RNATRIzol提取試劑盒（T102096）購自南京諾唯贊生物科技股份有限公司；TruSeq Stranded Total RNA with Ribo-Zero Gold（RS122-2301）購自美國Illumina公司；Agencourt AMPure XP（A63881）購自美國Beckman Coulter公司;Bioana- lyzer 2100 DNA-1000 Kit（5067-1504）購自美國Aglient公司;SuperScript Ⅱ Reverse Transcriptase（18064014） 購自美國Invitrogen公司;TransScript?All-in-One First- Strand cDNA Synthesis SuperMix for qPCR （AT341- 02）和PerfectStartTM Green qPCR SuperMix（AQ601）購自北京全式金生物技術股份有限公司。主要儀 器設備：ST16R冷凍離心機（美國ThermoFisher Sci- entific公司），NanoDrop 2000紫外分光光度計（美國 ThermoFisher Scientific公司），Centrifuge 5418R臺式 離心機（德國Eppendorf公司），Tanon 2500凝膠成像系統(tǒng)（上海蓓米爾生物科技有限公司），MyCyclerPCR儀（美國Bio-Rad公司），Qubit2.0熒光定量儀（美國 ThermoFisher Scientific公司），Bioanalyzer生化分析儀（安捷倫科技有限公司）。

1.2轉錄組測序

轉錄組測序分析委托上海歐易生物醫(yī)學科技有 限公司完成。采用TRIzol法提取藏公豬睪丸組織總RNA，利用NanoDrop 2000分光光度計檢測其濃度及OD250/OD230，以1%瓊脂糖凝膠電泳檢測RNA完整性。質檢合格的總RNA去除rRNA后進行片段化處理，逆轉錄生成cDNA，末端補平poly（A）尾及加接頭，PCR擴增后上機測序，測序原始數(shù)據(jù)經生物分 析前預處理后，進行基因組比對、轉錄本拼接及基因 表達定量分析。

1.3DEGs篩選與分析

通過HtSeq-count獲取每個樣本中比對到蛋白編碼基因上的Reads，再根據(jù)FPKM公式計算基因表達量（FPKM值）。FPKM（A）=（比對到基因A的 Fragments/比對到所有基因的Fragmentsx基因A長 度×10°）。然后計算差異倍數(shù)，并采用NB（負二項分布檢驗）進行差異顯著性檢驗，最終根據(jù)差異倍 數(shù)及差異顯著性檢驗結果篩選[DEGs錯誤發(fā)現(xiàn)率 （FDR）lt;0.05且|log2Fold Change|gt;1]，并對篩選獲得 的DEGs進行火山圖（Volcano Plot）分析、GO功能注釋分析和KEGG信號通路富集分析，根據(jù)分析結果篩選出藏公豬性成熟相關候選基因。

1.4實時熒光定量PCR驗證

使用 TransScript?All-in-One First-Strand cDNA Synthesis SuperMix for qPCR試劑盒將質檢合格的 總RNA反轉錄合成cDNA。反轉錄體系10.0μL：總 RNA 0.5 μg， 5×TransScript? All-in-One First-Strand cDNA Synthesis SuperMix for qPCR 2.0 uL， gDNA Remover 0.5 uL，Nuclease-free H20補足至10.0 uL。反轉錄程序：42℃15min，85℃5s。實時熒光定量PCR擴增引物（表1）由上海歐易生物醫(yī)學科技有限 公司設計，并委托北京擎科新業(yè)生物技術有限公司 合成。使用PerfectStartTM Green qPCR SuperMix 試 劑盒在LightCycler?480Ⅱ型熒光定量PCR儀上完 成實時熒光定量PCR驗證。反應體系10.0μL：2×PerfectStartTMGreen qPCR SuperMix 5.0 uL，上、下游引物（10umol/L）各0.2uL，cDNA模板1.0uL， Nuclease-free H2O3.6μL。擴增程序：94℃預變性30s；94℃5s，60℃30s，進行45個循環(huán)。以HMBS為內參基因，采用2-^Ac法計算目的基因相對表達量，采 用SPSS 19.0進行單因素方差分析（One-way ANOVA）。 2結果與分析

2.1轉錄組測序分析結果

對6份藏豬睪丸組織樣品的轉錄組測序共獲 得94.26 Gb的有效堿基（Clean bases），其中14.61～16.64 Gb的Clean bases占比較高，Q30分布范圍在93.28%～93.76％，GC含量為48.86%～50.03%（表2）。將有效序列（Clean reads）與參考基因組進行比對， 結果發(fā)現(xiàn)6個樣本的轉錄組測序數(shù)據(jù)比對率均 在90.00%以上，Clean reads的準確性和可靠性為90.07%～91.16%。此外，每個樣本原始堿基（Raw bases）過濾后獲得的Clean bases不少于14.50 Gb，且Q30gt;93.00%、GC含量gt;48.50%，說明轉錄組測序數(shù)據(jù)質量 較高，可用于后續(xù)研究。

2.2DEGs篩選結果

以性成熟前的藏豬睪丸組織樣品（TP01、TP02、TP03）為對照組，性成熟后的藏豬睪丸組織樣品（TP04、TP05、TP06）為試驗組，依據(jù)FDRlt;0.05且|log2Fold Change|gt;1的標準，共篩選出5929個DEGs，其中有3822個DEGs呈上調表達、2107個DEGs呈下調表達。性成熟前后藏豬睪丸組織DEGs火山圖如 圖1所示。

2.3DEGs的GO功能注釋分析結果

DEGs的GO功能注釋分析結果（圖2）顯示：在 生物學過程（Biological process）中，DEGs主要注 釋到生精過程（Spermatogenesis）、纖毛運動（Cilium movement）、鞭毛精子運動（Flagellated sperm moti- lity）、多細胞生物發(fā)育（Multicellular organism deve- lopment）、纖毛組裝（Cilium assembly）、細胞分化 （Cell differentiation）及精細胞發(fā)育（Spermatid deve- lopment）等GO功能條目；在細胞組分（Cellular com- ponent）中，DEGs主要注釋到基因絲（Axoneme）、運 動纖毛（Motile cilium）、纖毛（Cilium）、精子主塊（Sperm principal piecee）、精子頂囊（Acrosomal ve- sicle）、微管（Microtubule）和睫狀基體（Ciliary basal body）等GO功能條目；在分子功能（Molecular funca- tion）中，DEGs主要注釋到動力蛋白輕鏈結合 （Dynein light chain binding）、動態(tài)蛋白中間鏈結合（Dynein intermediate chain binding）、金屬內肽酶活 性（Metalloendopeptidase activity）、微管結合（Microtu- bule binding）、鈣離子結合（Calcium ion binding）等 GO功能條目。

2.4DEGs的KEGG信號通路富集分析結果

對DEGs進行KEGG信號通路分析，并利用超幾何分布檢測評估各通路中DEGs的富集情況，得到前20條KEGG信號通路富集分析氣泡圖。如圖3所示，DEGs主要富集于受體相互作用通路（ssc04512：ECM-receptor interaction）、軸突引導通路（ssc04360： Axon guidance）、神經活性配體一受體相互作用通路 （ssc04080： Neuroactive ligand-receptor interaction）、鈣信號通路（ssc04020：Calcium signa-ling pathway）、 心律失常性右心室心肌通路（ssc05412：Arrhythmo- genic right ventricular cardiomyopathy）、肥厚性心肌 病通路（ssc05410：Hypertrophic cardiomyopathy）、化學致癌-DNA加合物通路（ssc05204：Chemical carcino-genesis-DNA adducts）、礦物吸收通路（ssc04978：Mine-ral absorption）、PPAR信號通路（ssc03320：PPAR sig-naling pathway）及消化吸收通路（ssc04974：Protein digestion and absorption）等。

2.5藏公豬性成熟相關基因篩選結果

綜合DEGs的GO功能注釋分析和KEGG信號 通路富集分析結果，共篩選出9個與藏公豬性成熟 相關的基因，分別是凋亡蛋白基因（BAX）、雙性和mab-3相關轉錄因子1（DMRT1）、卵泡刺激素受體基因（FSHR）、胰島素樣生長因子1受體基因（IGFIR）、精卵融合必需蛋白1基因（IZUMO1）、人源類溶菌酶蛋白4基因（LYZL4）、人源類溶菌酶蛋白6基因（LYZL6）、突觸復合蛋白3基因（SYCP3）和過渡蛋白2基因（TNP2）。其中，BAX基因在線粒體凋亡過程 中發(fā)揮作用；DMRTI基因通過控制睪丸發(fā)育和生殖細胞增殖，在雄性性別決定和分化中起關鍵作用； FSHR基因通過cAMP的產生激活下游PI3K-Akt和ERK1/ERK2信號通路；IGF1R是胰島素樣生長因子 1受體a鏈，介導胰島素樣生長因子1作用的受體酪 氨酸激酶；IZUMO1基因是受精的必需因子；LYZL4基因可能參與受精，而LYZL6基因可能與受精過程 中精卵質膜的黏附與融合有關；SYCP3是雄性生 殖細胞減數(shù)分裂過程中著絲粒配對的必需蛋白； TNP2基因在精細胞伸長與濃縮過程中發(fā)揮重要 作用。

2.6實時熒光定量PCR驗證結果

以性成熟前的藏豬睪丸組織樣品（TP01、TP02、TP03）為對照組，性成熟后的藏豬睪丸組織樣品（TP04、TP05、TP06）為試驗組，通過實時熒光定量PCR檢測篩選出的9個藏公豬性成熟相關基因，設性成熟前藏豬睪丸組織的相對表達量為1，以2-Ac法 計算目的基因在性成熟后藏豬睪丸組織的相對表達 量，結果（圖4）顯示，BAX、DMRT1、FSHR、IGFIR等 4個基因在性成熟后藏豬睪丸組織中的相對表達量 較性成熟前呈明顯的下降趨勢，且下降幅度均在50%以上；IZUMO1、LYZL4、LYZL6、SYCP3、TNP2等 5個基因的相對表達量較性成熟前上升了1.07～30.25倍，其中IZUMO1、LYZL4和LYZL6基因的上升倍數(shù)達極顯著水平（Plt;0.01）。同時，計算9個基因的表達差異倍數(shù)（表3），發(fā)現(xiàn)實時熒光定量PCR檢測結果與轉錄組測序結果基本一致，即BAX、DMRT1、 FSHR、IGF1R等4個基因呈下調表達，IZUMO1、LYZL4、LYZL6、SYCP3、TNP2等5個基因呈上調表達，進一步 證實轉錄組測序結果的準確性。

3討論

本研究對性成熟前后的藏公豬睪丸組織進行轉 錄組測序，共篩選獲得5929個DEGs，其中有3822個DEGs呈上調表達、2107個DEGs呈下調表達。經GO功能注釋分析和KEGG信號通路富集分析，發(fā)現(xiàn) 5929個DEGs主要注釋到與精子發(fā)生和精子運動等 有關的GO功能條目，且主要富集在與精子形成相關的經典信號通路上。綜合DEGs的GO功能注釋分析和KEGG信號通路富集分析結果，最終篩選出9個與藏公豬性成熟相關的基因，分別是BAX、DMRT1、FSHR、IGF1R、IZUMO1、 LYZL4、LYZL6、 SYCP3和TNP2。這9個候選基因在前人的相關研究中已有報道。LYZL4和LYZL6基因參與細胞壁大 分子分解代謝過程及細胞外空間通路，李忠邦和閻 萍（2019）研究發(fā)現(xiàn)LYZL4、LYZL6基因在牦牛睪丸/附睪中高表達，LYZL4基因在牦牛性成熟前存在于 間質細胞中，牦牛性成熟后則存在于精原細胞和支 持細胞中，有助于精子產生及其成熟過程；LYZL6基 因表達于睪丸和附睪中，對受精過程起重要作用。 本研究發(fā)現(xiàn)，在卵巢類固醇生成通路中FSHR和IGFIR基因在性成熟后個體睪丸組織中高表達。FSHR基因作用于精子的產生過程，且對性腺細胞有營養(yǎng)支持作用，可通過FSHR調節(jié)卵巢功能；IGFIR 基因是影響動物生長發(fā)育的主要候選基因，在免疫 調節(jié)及肌肉和骨骼發(fā)育生長過程中發(fā)揮重要作用 （Ding et al.，2022）。蘇少鋒等（2019）研究證實，IGFIR基因主要存在于睪丸間質細胞的細胞質中。龍亞麗等（2022）研究表明，F(xiàn)SHR基因主要表達于支 持細胞、間質細胞及各種生精細胞中，與曲細精管和 精子形成有密切關系。

TNP2基因在精子的形成過程中發(fā)揮重要作用。藺蕙等（2021）研究表明，在性成熟牦牛睪丸組織中 的TNP2基因相對表達量明顯高于其他發(fā)育階段，且 在圓形精子細胞和長形精子細胞中的含量較高。Amjad等（2021）研究發(fā)現(xiàn)，TNP2基因新型單核苷酸多態(tài)性與男性不育相關，隨著睪丸組織中TNP2基因 表達譜的增加，無精子癥患者恢復生精的概率提高。BAX、SYCP3和TNP2基因在脫氧核糖核酸結合及精 子DNA縮合通路中高表達（Jebur et al.，2023）。BAX基因作為誘導凋亡因子，主要作用于線粒體膜，通過促進細胞色素釋放，導致生精細胞凋亡；當BAX基因表達下降時，導致生精細胞大量存在而產生非BAX凋亡過程（Xi et al.，2017）。DMRT1基因與雄性性別 決定相關，參與雄性哺乳動物的性腺發(fā)育過程，當 DMRTI基因缺失時通過激活維甲酸信號通路，導致 干細胞向卵巢細胞分化，同時能調節(jié)精原細胞有絲 分裂及減數(shù)分裂（Zhang et al.，2016）。Li等（2018）研究表明，DMRTI基因在性成熟前后綿羊睪丸組織 中均有表達，且在性成熟后綿羊睪丸組織中的表達 水平極顯著高于性成熟前。

精子與卵子融合是通過特定的膜蛋白配子相互 作用而完成，包括IZUMO1和受精介導蛋白（JUNO）等特定的膜蛋白。精子IZUMO基因依賴配子融合，會影響小鼠的雄性生育能力（Saito et al.，2019）。有研究證明，隨著年齡的增長，IZUMO1基因表達量逐漸上調，且在性成熟前后差異顯著（Kalwar et al.，2019b）。SYCP3是一種DNA結合蛋白，位于突觸復 合物的側向元件上，在同源染色體配對中發(fā)揮重要 作用，是精子發(fā)生減數(shù)分裂的必需蛋白（Wang et al.， 2011）。精子細胞更新與分化的動態(tài)過程發(fā)生在生 精小管內，分別受有絲分裂和減數(shù)分裂的調節(jié)。在 減數(shù)分裂過程中，基因重組是由突觸復合物的組裝 引起，而突觸復合物涉及組成蛋白的參與，如SYCP3等。Kalwar等（2019a）研究表明，SYCP3基因在成熟牦牛睪丸組織中高表達，進一步證實SYCP3是精子 發(fā)生、誘導及同源重組的必需蛋白。上述研究結論 為今后深入開展藏豬繁殖育種能力相關研究提供了 重要的基礎數(shù)據(jù)。

4結論

基于轉錄組測序從藏豬睪丸組織中篩選獲得9個與性成熟相關的DEGs，即BAX、DMRTI、FSHR、 IGF1R、IZUMO1、LYZL4、LYZL6、SYCP3和TNP2，可 為作為研究藏豬繁殖育種能力的主要候選基因。

參考文獻：

董世雄，王浩田，姚雨民，徐士軍，王向東，張健，商鵬.2022.

藏豬與大約克夏豬Mxl基因CDS區(qū)多態(tài)性及生物信息 學分析[J].黑龍江畜牧獸醫(yī)，（4）：35-41.[DongSX，Wang H T， Yao Y M，Xu S J，Wang X D，Zhang J，Shang P. 2022. Polymorphism and bioinformatics analysis of Mxl gene CDS region in Tibetan pigs and Yorkshire pigs[J]. Heilong- jiang Animal Science and Veterinary Medicine， （4）：35-41.] doi：10.13881/j.cnki.hljxmsy.2021.05. 0037.

董新星，王琳，劉金橋，嚴達偉，張博，鄧俊，張浩，馬黎.2023.

大型迪慶藏豬不同生長階段腹脂代謝功能基因篩選及 其調控網絡分析[J].云南農業(yè)大學學報（自然科學），38

（1） ： 61-70. [Dong X X， Wang L， Liu J Q， Yan D W， Zhang B， Deng J，Zhang H， Ma L. 2023. Analysis of functional genes screening and regulatory networks of abdominal fat metabolism of large Diqing Tibetan pigs at different growth stages [J]. Journal of Yunnan Agricultural University（Natural Science），38（1）：61-70.] doi： 10.12101/j.issn.

1004-390X（n）.202204047.

甘麥鄰，沈林園，楊東麗，王定國，張順華，朱礪.2020.轉錄組

學技術及其在豬肉質遺傳研究中的應用[J].豬業(yè)科學，

37（2）：94-98. [Gan M L，Shen L Y， Yang D L，Wang D G， Zhang S H，Zhu L. 2020. Application of transcriptomics to improve pork quality traits [J]. Swine Industry Science，37

（2）：94-98.] doi：10.3969/j.issn.1673-5358.2020.02.026.

韓雨晴，薄素雪，嚴飛飛，段夢琪，強巴央宗，商鵬.2023.不同 品種豬PVALB基因多態(tài)性及其在肌肉組織的表達差異 分析[J].河南農業(yè)科學，52（5）：150-155.[Han YQ，Bo S X， Yan F F， Duan M Q， Chamba Yangzom， Shang P. 2023. Analysis of porcine PVALB gene polymorphism and diffe- rential expression in muscle tissues of different pig breeds [J]. Journal of Henan Agricultural Sciences， 52 （5）： 150- 155.] doi：10.15933/j.cnki.1004-3268.2023.05.017.

李忠邦，閻萍.2019.不同月齡牦牛睪丸中LYZL4、LYZL6的基

因表達及蛋白質定位[J].河南農業(yè)科學，48（7）：116-121.

[Li Z B， Yan P. 2019. Gene expression and protein localization of LYZL4 and LYZL6 in testes of different months

old yak[J]. Journal of Henan Agricultural Sciences，48（7）：

116-121.] doi： 10.15933/j.cnki.1004-3268.2019.07.017.

梁婉，周丹娜，彭忠，楊克禮，段正贏，郭銳，劉威，吳斌，陳文 杰，孟麗，田永祥.2019.轉錄組測序技術在豬繁殖與呼 吸綜合征病毒研究中的應用[J].中國畜牧獸醫(yī)，46（1）： 37-44. [Liang W， Zhou D N， Peng Z， Yang K L， Duan Z Y， Guo R， Liu W，Wu B， Chen W J， Meng L， Tian Y X. 2019. Application of transcriptome sequencing technology in porcine reproductive and respiratory syndrome virus research[J]. China Animal Husbandry amp; Veterinary Medi- cine，46（1）：37-44.] doi：10.16431/j.cnki.1671-7236.2019.

01.005.

藺蕙，潘陽陽，韓小紅，蔡得琪，王萌，盧國杰，柯良備，王瑞 剛，王學義，余四九.2021.牦牛TNP2基因分子特征及其

在不同發(fā)育階段睪丸中的表達[J].畜牧與獸醫(yī)，53（11）： 6-13. [Lin H，Pan Y Y，Han X H，Cai D Q，Wang M，Lu G J，Ke L B， Wang R G， Wang X Y， Yu S J. 2021. Molecular characteristics of TNP2 gene in yaks and its quantitation in testis at different developmental stages [J]. Animal Hus- bandry amp; Veterinary Medicine，53（11）：6-13.]

劉思源，何辰慶，陳軍軍，葉幼榮，張健，商鵬，強巴央宗.

2022.藏豬和大約克夏豬CXCL10基因多態(tài)性及表達差

異分析[J].中國畜牧雜志，58（11）：100-104.[LiuSY， He C Q， Chen J J， Ye Y R， Zhang J， Shang P， Chamba Yangzom. 2022. Analysis of CXCL10 gene polymorphism and expression differences between Tibetan and Yorkshire

pigs[J]. Chinese Journal of Animal Science， 58 （11） ： 100-

104.] doi：10.19556/j.0258-7033.20211112-12.

龍亞麗，王軍乾，潘陽陽，王萌，趙凌，馬文斌，張虔，王靖雷， 余四九，王立斌.2022.牦牛睪丸發(fā)育過程中INHB、

FSHR和LHR的表達差異性分析[J].動物醫(yī)學進展，43

（12） ： 43-49. [Long L Y， Wang J Q， Pan Y Y， Wang M， Zhao L， Ma W B，Zhang Q，Wang J L， Yu S J，Wang L B.

2022. Differential analysis of expressions of INH?B， FSHR and LHR during development of testes in yaks [J].

Progress in Veterinary Medicine， 43 （12） ： 43-49.] doi： 10.

3969/j.issn.1007-5038.2022.12.008.

馬圍圍，宋尚橋，曾素先，孫翠翠，李鑫，王睿敏，嚴瑾，黎宗 強.2021.豬肌肉的轉錄組學研究進展[J].基因組學與

應用生物學，40（1）：159-164.[Ma W W，Song SQ，Zeng S X，Sun C C，Li X， Wang R M， Yan J，Li Z Q. 2021. Prog- ress in research on porcine muscle mRNA transcriptome [J]. Genomics and Applied Biology，40（1）： 159-164.] doi：

10.13417/j.gab.040.000159.

聶靖茹，張博，馬黎，張浩，李國美，嚴達偉，董新星.2022.大 型迪慶藏豬不同生長階段肌肉生長差異基因及其調控 通路分析[J].中國農業(yè)大學學報，27（6）：132-144.[Nie J R，Zhang B，Ma L，Zhang H，Li G M，Yan D W，Dong X

X. 2022. Screening of differential genes for muscle growth and regulatory pathways analysis in large Diqing Tibetan pig at different growth stages[J]. Journal of China Agricultural University，27（6）：132-144.] doi：10.11841/j.issn.1007-4333.2022.06.12.

覃海，陳偉，畢歡，袁巍，周萍.2022.基于轉錄組測序對劍白 香豬黑色素生成相關基因的篩選與分析[J].中國獸醫(yī)學

報，42（11）：2261-2271.[Qin H，Chen W，Bi H，Yuan W， Zhou P. 2022. Screening and analysis of genes related to melanin production in Jianbaixiang pig based on transcrip- tome sequencing[J]. Chinese Journal of Veterinary Science，

42（11）：2261-2271.] doi：10.16303/j.cnki.1005-4545.2022.

11.21.

宋明坤，薛明明，張力戈，顏鐸，夏寧，商鵬，李新建，韓雪蕾， 喬瑞敏，李秀領，李明，王克君.2022.影響藏豬與大約克 夏豬肌肉品質的lncRNA的篩選及功能分析[J].畜牧獸 醫(yī)學報，53（1）：53-65.[Song MK，Xue MM，Zhang LG， Yan D，Xia N，Shang P，Li X J，Han X L，Qiao R M， Li X L， Li M， Wang K J. 2022. Screening and functional analy- sis of IncRNA affecting muscle quality in Tibetan and Yorkshire pigs[J]. Acta Veterinaria et Zootechnica Sinica，

53 （1） ： 53-65. ] doi ： 10.11843/j.issn.0366-6964.2022.01.

006.

蘇少鋒，趙一萍，李蓓，白東義，宋海燕，趙若陽，芒來.2019.

不同發(fā)育時期蒙古馬睪丸組織中GH和IGF-I及其受體

的表達[J].中國農業(yè)大學學報，24（1）：52-60.[SuSF，

Zhao Y P，Li B，Bai D Y， Song H Y，Zhao R Y，Mang L. 2019. Expression profiles of GH， IGF-I and their corresponding receptors in testicular tissue of Mongolian horses at different developmental stages[J]. Journal of China Agricultural University， 24 （1） ： 52-60. ] doi ： 10.11841/j.issn.

1007-4333.2019.01.08.

孫菁希，逄世龍，孫敬春，褚瑰燕.2021.基于轉錄組測序鑒定

豬肌內脂肪沉積相關生物學通路及候選基因[J].中國畜

牧獸醫(yī)，48（12）：4508-4519.[SunJX，Pang SL，SunJ C， Chu G Y. 2021. Identification of biological pathways and candidate genes related to porcine intramuscular fat deposition using transcriptome sequencing [J]. China Ani- mal Husbandry amp; Veterinary Medicine，48（12）：4508-4519.]

doi：10.16431/j.cnki.1671-7236.2021.12.019.

王琳，馬黎，張博，鄧俊，張浩，歐陽曉芳，嚴達偉，董新星.

2023.大型迪慶藏豬不同生長階段背脂與腹脂脂質代謝 差異基因及調控網絡分析[J].畜牧獸醫(yī)學報，54（2）： 520-533. [Wang L， Ma L， Zhang B， Deng J， Zhang H， Ouyang X F， Yan D W， Dong X X. 2023. Key genes and regulatory network analysis of lipid metabolism diffe- rences between back fat and abdominal fat of large Diqing Tibetan pigs at different growth stages[J]. Acta Veterinaria et Zootechnica Sinica， 54 （2） ： 520-533. ] doi ： 10.11843/j.

issn.0366-6964.2023.02.010.

王璐夢，孫秀秀，朱陽華，王金坡，張萬坡，胡薛英，谷長勤，程 國富.2021.人工感染豬繁殖與呼吸綜合征病毒豬腦組 織的轉錄組學分析[J].中國獸醫(yī)科學，51（3）：368-373.

[Wang L M，Sun X X，Zhu Y H，Wang J P，Zhang W P， Hu X Y， Gu C Q， Cheng G F. 2021. Transcriptome analysis of porcine brain tissue artificially infected with porcine reproductive and respiratory syndrome virus[J]. Chinese Veterinary Science，51 （3） ： 368-373.] doi： 10.16656/j.issn.

1673-4696.2021.0016.

韋明邦，段夢琪，羅布白珍，肖青青，徐士軍，陳軍軍，葉幼榮， 張健，商鵬.2022.藏豬與大約克豬EDN1基因多態(tài)性與

表達差異分析[J].中國畜牧雜志，58（3）：79-83.[Wei M

B， Duan M Q， Luobubaizhen，Xiao Q Q，Xu S J， Chen J J， Ye Y R， Zhang J， Shang P. 2022. Analysis of EDNl gene polymorphism and expression difference between Tibetan pig and Gram pig[J]. Chinese Journal of Animal Science，

58（3）：79-83.] doi：10.19556/j.0258-7033.20210611-02.

熊和麗，張彥，張斌，相德才，常雅潔，劉韶娜，趙智勇.2023.

基于全基因組重測序分析藏豬骨骼肌發(fā)育相關基因[J].

中國農業(yè)大學學報，28（3）：83-93.[Xiong HL，Zhang Y， Zhang B，Xiang D C，Chang Y J， Liu S N，Zhao Z Y. 2023. Analysis of genes related to skeletal muscle development in Tibetan pigs based on whole genome sequencing[J]. Journal of China Agricultural University， 28 （3） ： 83-93.] doi： 10.11841/j.issn.1007-4333.2023.03.08.

楊燁城，黃秋艷，辛海云，李寶紅，杜宗亮，孟繁明，朱向星，李 劍豪，王塑天，唐冬生.2022.藏豬和杜藏豬生長、肉質性 狀及脂肪沉積相關基因的表達差異分析[J].南方農業(yè)學

報，53（4）：926-934.[Yang Y C，HuangQY，Xin H Y，Li B H， Du Z L， Meng F M， Zhu X X， Li J H， Wang S T， Tang D S. 2022. Analysis of differences in growth， meat quality traits and expression of genes related to fat deposi- tion between Tibetan and Duzang pigs[J]. Journal of Sou-

thern Agriculture ， 53 （4） ： 926-934. ] doi ： 10.3969/j.issn.

2095-1191.2022.04.005.

易凡利，王嘉福，張福平，阮亦麒，唐靚婷，毛寧，冉雪琴.

2018.5個繁殖候選基因在香豬卵巢和子宮中的差異表 達[J].中國畜牧雜志，54（5）：30-33.[Yi F L，WangJF， Zhang F P， Ruan Y Q， Tang L T， Mao N，Ran X Q. 2018. Differential expression of five reproductive candidate genes in the ovary and uterus of Xiang pig[J]. Chinese Journal of Animal Science， 54 （5） ： 30-33.] doi： 10.19556/

j.0258-7033.2018-05-030.

翟麗維，趙延輝，李文軍，邢凱，王楚端.2022.系統(tǒng)分析多組 織轉錄組鑒定影響豬脂肪沉積的關鍵基因[J].畜牧獸醫(yī)

學報，53（6）：1702-1711.[Zhai LW，Zhao YH，Li WJ， Xing K，Wang C D. 2022. System analysis of multi tissue transcriptome to identify key genes affecting porcine fat deposition [J]. Acta Veterinaria et Zootechnica Sinica， 53

（6）：1702-1711.] doi：10.11843/j.issn.0366-6964.2022.06.

005.

張福平，唐靚婷，王嘉福，冉雪琴，李升，黃世會.2021.基于轉

錄組測序篩選影響從江香豬產仔數(shù)的候選基因[J].南方

農業(yè)學報，52（4）：847-856.[Zhang F P，TangLT，WangJ F， Ran X Q， Li S， Huang S H. 2021. Selection of candidate genes affecting litter size of Congjiang Xiang pig by tran- scriptome sequencing[J]. Journal of Southern Agriculture，

52 （4）： 847-856.] doi： 10.3969/j.issn.2095-1191.2021.04.

001.

Amjad S， Mushtaq S， Rehman R， Zahid N， Munir A， Siddiqui P Q R. 2021. Spermatozoa retrieval in azoospermia and expression profile of JMJD1A， TNP2， and PRM2 in a subset of the Karachi population [J]. Andrology， 9 （6） ： 19341942. doi：10.1111/andr.13076.

Ding N，Tian D H，Li X，Zhang Z C，Tian F，Liu S J，Han B Y， Liu D H， Zhao K. 2022.Genetic polymorphisms of IGF1 and IGFIR genes and their effects on growth traits in Hulun Buir sheep[J]. Genes， 13 （4） ： 666. doi： 10.3390/ GENES13040666.

Jebur M I I， Dastmalchi N， Banamolaei P， Safaralizadeh R.

2023. Polymorphisms and expression levels of TNP2，

SYCP3， and AZFa genes in patients with azoospermia [J].

Clinical and Experimental Reproductive Medicine， 50（4） ：

253-261. doi： 10.5653/CERM.2023.06219.

Kalwar Q，Chu M，Ahmad A A，Ma X M，Zhang R Z，Ma F L， Xie J P，Ding X Z，Wu X Y，Bao P G， Yan P. 2019a. Molecular cloning and characterization of SYCP3 and TSEG2 genes in the testicles of sexually mature and immature yak [J]. Genes， 10（11）：867. doi：10.3390/genes10110867.

Kalwar Q，Ding X Z， Ahmad A A， Chu M， Wu X Y， Bao P J， Yan P. 20196. Expression analysis of IZUMO1 gene during testicular development of Datong yak （Bos grunniens）[J].

Animals，9（6）：292. doi： 10.3390/ani9060292.

Li T T， Lu Z K， Luo R R， Gao J F， Zhao X X， Ma Y J. 2018. Expression and cellular localization of double sex and mab-3 related transcription factor 1 in testes of postnatal SmallTail Han sheep at different developmental stages[J]. Gene，

642：467-473. doi：10.1016/j.gene.2017.11.053.

Saito T， Wada I， Inoue N. 2019. Sperm IZUMO1-dependent gamete fusion influences male fertility in mice[J]. International Journal of Molecular Sciences， 20 （19） ： 4809. doi：

10.3390/ijms20194809.

Wang S， Pan Z， Zhang Q， Xie Z， Liu H， Li Q. 2011. Differential mRNA expression and promoter methylation status of SYCP3 gene in testes of yaks and cattle-yaks[J]. Reproduction in Domestic Animals， 47 （3） ： 455-462. doi： 10.1111/

j.1439-0531.2011.01902.×.

Xi H，F(xiàn)an X，Zhang Z， Liang Y， Li Q，He J. 2017. Bax and Bcl-2 are involved in the apoptosis induced by local testicular heating in the boar testis[J]. Reproduction in Domestic Animals，52（3）：359-365. doi： 10.1111/rda.12904.

Zhang T，Oatley J，Bardwell V J，Zarkower D. 2016. DMRT1 is required for mouse spermatogonial stem cell maintenance and replenishment [J]. PLoS Genetics， 12 （9） ： e1006293.

doi： 10.1371/journal.pgen.1006293.

（責任編輯

蘭宗寶）