Genome Biol | 時(shí)空組學(xué)技術(shù)在豬早期卵子發(fā)生階段的應用研究

2025年1月，青島農業(yè)大學(xué)葛偉教授與沈偉教授團隊在學(xué)術(shù)期刊《Genome Biology》上發(fā)表了一篇題為《Spatiotemporal dynamics of early oogenesis in pigs》的研究論文。該研究利用了單細胞RNA測序和空間轉錄組技術(shù)，解析了豬早期卵子發(fā)生中的基因表達及空間組織特征。研究發(fā)現，豬生殖細胞在卵巢中呈“皮質(zhì)到髓質(zhì)”分布，且與人類(lèi)卵子發(fā)生過(guò)程相似，表明豬是研究人類(lèi)卵子發(fā)生的良好模型。RNA速度分析揭示了顆粒細胞譜系特征，空間共定位分析和細胞通訊分析顯示皮質(zhì)與髓質(zhì)區域生殖細胞與體細胞通訊模式不同，特別是NOTCH信號通路和ECM蛋白在調控中起關(guān)鍵作用，強調了卵巢微環(huán)境對生殖細胞命運的重要性。

文章標題：Spatiotemporal dynamics of early oogenesis in pigs

期刊名稱(chēng)：Genome Biology

影響因子：10.1

合作單位：青島農業(yè)大學(xué)

發(fā)表日期：2025.01.02

研究方法：單細胞轉錄組（10x Chromium）、空間轉錄組（10x Visium）、高分辨率空間轉錄組（BMKMANU S1000）、免疫熒光

百邁客生物為該研究中的E65胚胎卵巢提供了BMKMANU S1000空間轉錄組技術(shù)服務(wù)。

 

研究背景

在發(fā)育中的哺乳動(dòng)物卵巢中，卵泡的位置對其生物學(xué)功能具有重要影響。休眠卵泡主要位于卵巢皮質(zhì)外部，對雌性生育周期長(cháng)短具有重要影響，而生長(cháng)卵泡則位于髓質(zhì)區域，對雌性生育的起始至關(guān)重要。然而，目前對胚胎階段卵巢微環(huán)境形成機制的理解尚不深入，這主要是由于該過(guò)程發(fā)生在胚胎期，涉及多種異質(zhì)細胞類(lèi)型的協(xié)同作用。

在小鼠中，原始生殖細胞從外胚層遷移到生殖脊并在性腺中增殖，隨后在特定時(shí)間點(diǎn)沿前后軸異步啟動(dòng)減數分裂。最近的研究顯示，這一過(guò)程與經(jīng)典的維甲酸信號無(wú)關(guān)，而是通過(guò)TEX14形成細胞間橋來(lái)實(shí)現。這些發(fā)現強調了生殖細胞微環(huán)境在決定其發(fā)育命運中的關(guān)鍵作用。

目前對卵巢發(fā)育中卵母細胞的研究主要集中在基因表達上，但對卵母細胞在卵泡發(fā)生過(guò)程中的細微空間定位及主要體細胞類(lèi)型的空間動(dòng)力學(xué)特征了解有限。通過(guò)空間轉錄組學(xué)和單細胞RNA測序技術(shù)，該研究揭示了豬早期卵巢發(fā)育中的皮質(zhì)-髓質(zhì)梯度，并發(fā)現人豬之間卵母細胞空間定位的保守性，進(jìn)一步強調了卵巢微環(huán)境在決定卵母細胞命運中的重要作用。

 

材料方法

研究材料：通過(guò)人工授精培育妊娠豬（三元雜交：Landrace、Large White和Duroc），分別在受精后45、55、65和75天采集胚胎階段的豬卵巢樣本。

研究方法：scRNA-seq（E45、E55、E65、E75胚胎卵巢，共計18,956個(gè)高質(zhì)量細胞，10x Genomics Chromium）；空間轉錄組（E45、E55、E65、E75胚胎卵巢，10x Visium；E65胚胎卵巢，BMKMANU S1000）

分析方法：空間轉錄組解卷積（Cell2location）；RNA速度分析（Velocyto和scVelo）；單細胞軌跡推斷（Cellrank）；同源基因轉化（biomaRt）；細胞通訊分析（CellphoneDB、biomaRt、ktplots）；蛋白互作分析；GO富集分析；基因模塊分析（Hotspot）

驗證實(shí)驗：免疫熒光；圖像分析和定量分析；Western blotting；卵巢組織體外培養

 

研究結果

1.豬早期卵子發(fā)生的單細胞圖譜的構建

為了深入揭示豬卵巢早期卵子發(fā)生的時(shí)空發(fā)育特征，研究者將單細胞RNA測序（scRNA-seq）與10× Visium空間轉錄組（ST）技術(shù)相結合，對豬卵巢發(fā)育過(guò)程中的細胞空間屬性進(jìn)行了系統分析。4個(gè)發(fā)育時(shí)間點(diǎn)（E45、E55、E65和E75；scRNA-seq各包括2個(gè)胚胎；ST：E45、E55、E65各包括2個(gè)胚胎，E75，1個(gè)胚胎），其中包括早期卵子發(fā)生的各階段。保留了18,956個(gè)高質(zhì)量細胞，共鑒定出17個(gè)聚類(lèi)，包含9種主要細胞類(lèi)型：生殖細胞、雙潛能前顆粒細胞、上皮前顆粒細胞、性腺間質(zhì)細胞、平滑肌細胞、內皮細胞、血管周?chē)毎?、T細胞、巨噬細胞。

2.豬卵巢空間轉錄組的細胞類(lèi)型解卷積

為了揭示卵巢發(fā)育過(guò)程中不同細胞類(lèi)型的空間動(dòng)態(tài)變化，研究者進(jìn)行了細胞類(lèi)型的解卷積分析。盡管scRNA-seq技術(shù)為剖析復雜組織中的細胞異質(zhì)性提供了有力工具，但它缺乏細胞的空間位置信息。因此，研究者采用10× Visium空間轉錄組（ST）技術(shù)進(jìn)行測序，以更好地展現豬卵巢發(fā)育過(guò)程中不同細胞類(lèi)型的空間動(dòng)態(tài)變化。研究者首先對ST數據進(jìn)行了質(zhì)量評估，確保數據質(zhì)量符合分析要求。隨后，利用scRNA-seq數據作為參考，通過(guò)Cell2location算法對ST數據進(jìn)行了細胞類(lèi)型的解卷積，實(shí)現了細胞類(lèi)型的高分辨率定位。

通過(guò)解卷積分析，研究者發(fā)現在E45和E55階段，細胞類(lèi)型的空間分布相對較為“隨機”。然而，隨著(zhù)發(fā)育進(jìn)展到E65階段，研究者觀(guān)察到了一個(gè)明顯的空間分布模式：生殖細胞在E45-E55階段隨機分布在豬卵巢中，而在E65-E75階段則主要定位于皮層區域，僅有少量位于髓質(zhì)區域。對于前顆粒細胞群，我們觀(guān)察到了兩種具有不同空間定位模式的細胞類(lèi)型：BPG細胞在E55之前隨機分布，并逐漸在E75時(shí)集中于髓質(zhì)區域；而EPG細胞則在整個(gè)發(fā)育過(guò)程中始終位于卵巢表面。

為進(jìn)一步驗證分析結果，研究者進(jìn)行了全組織染色實(shí)驗，使用生殖細胞標志物DDX4和前顆粒細胞標志物KRT19對E45至E75的胎兒卵巢進(jìn)行染色，實(shí)驗結果與ST數據分析結果一致。此外，研究者還分析了這些陽(yáng)性細胞在卵巢背腹軸上的空間分布，并通過(guò)熒光強度分析確認了生殖細胞和前顆粒細胞在皮層和髓質(zhì)區域的不同分布模式。這些發(fā)現為深入理解豬卵巢發(fā)育過(guò)程中細胞類(lèi)型的空間動(dòng)態(tài)變化提供了重要線(xiàn)索。

3.空間尺度解析精細尺度豬生殖細胞發(fā)育軌跡

為了深入理解生殖細胞在空間環(huán)境中的發(fā)育過(guò)程，研究者提取了生殖細胞進(jìn)行重新聚類(lèi)。通過(guò)對生殖細胞群的精細分析，研究者發(fā)現了五個(gè)對應的細胞簇，分別是有絲分裂生殖細胞（FGC_mitotic）、表達MECOM和ATM的前減數生殖細胞（Oogonia_STRA8）、表達SYCP1和DMC1的減數分裂生殖細胞（Oogonia_meiotic）、表達FIGLA和NANOS1的早期卵母細胞（Pre_oocyte）以及表達ZP4和ZP3的卵母細胞（Oocyte）。

為了深入了解早期豬卵母細胞生成過(guò)程中模塊基因的功能富集情況，研究人員采用Hotspot算法，利用生殖細胞的scRNA-seq數據識別了相關(guān)基因模塊，共鑒定出11個(gè)功能各異的基因模塊。通過(guò)將單個(gè)模塊投影到UMAP圖中，確定了具有簇特異性的信息性基因模塊。剖析了早期豬卵母細胞生成過(guò)程中的基因表達譜后，研究人員進(jìn)一步在空間環(huán)境中研究了這一過(guò)程。通過(guò)對ST芯片上的生殖細胞進(jìn)行解卷積，揭示了早期卵母細胞生成的皮質(zhì)-髓質(zhì)梯度。為了驗證上述分析，研究人員使用生殖細胞標記DDX4和減數分裂生殖細胞細線(xiàn)期標記γH2AX進(jìn)行了全組織染色實(shí)驗，以定位在不同階段啟動(dòng)減數程序的生殖細胞。評估了C-M軸的熒光強度后發(fā)現在E45和E55卵巢中，γH2AX信號沿C-M軸的水平相似，而在E65至E75卵巢中，γH2AX信號在內皮層達到峰值，在髓質(zhì)區域則減少。

為了進(jìn)一步驗證分析，研究人員還在同一E65卵巢組織上使用了更高分辨率（BMKMANU S1000；Spot點(diǎn)大小約為10微米）的平臺進(jìn)行了ST測序。結果與E65時(shí)的Visium ST解卷積一致，使用S1000平臺在E65時(shí)可視化ZP3陽(yáng)性卵母細胞的空間位置也證實(shí)了其在內皮層的表達。

綜上所述，研究人員成功地在單細胞分辨率下再現了豬生殖細胞的發(fā)育過(guò)程，并闡明了豬卵巢組織中生殖細胞的空間分布模式。

4.跨物種分析揭示了豬與人類(lèi)之間生殖細胞基因表達程序及空間位置模式的保守性

為了深入理解哺乳動(dòng)物卵母細胞生成的空間調控機制，研究人員進(jìn)行了跨物種分析。研究團隊首先獲取了妊娠后19周人類(lèi)卵巢的空間轉錄組（ST）數據，并將其與豬E65時(shí)期的ST數據進(jìn)行了比較，以分析卵母細胞生成各階段生殖細胞的空間位置模式。與人類(lèi)情況一致的是，生殖細胞標記DDX4與5mC和5hmC的共染色顯示，在E55和E75豬卵巢中，5mC和5hmC信號在DDX4陽(yáng)性細胞中不可檢測，且主要在周?chē)韵袤w細胞中表達。這些結果進(jìn)一步揭示，除了保守的基因表達程序外，人類(lèi)和豬在早期卵母細胞生成過(guò)程中的表觀(guān)遺傳重編程模式也是保守的。

為了研究豬和人類(lèi)早期卵母細胞生成過(guò)程中皮質(zhì)-髓質(zhì)梯度是否保守，研究人員基于ST數據探索了卵母細胞生成各階段生殖細胞的空間位置模式。對于早期卵母細胞階段的生殖細胞，這些細胞在卵巢內皮質(zhì)區域有清晰定位，而在人類(lèi)中，這些細胞主要位于內皮質(zhì)和外髓質(zhì)區域，且在內皮質(zhì)區域豐度更高。到了卵母細胞階段，豬和人類(lèi)中這些生殖細胞的位置模式相似，主要位于外髓質(zhì)區域。通過(guò)評估生殖細胞到卵巢表面的相對距離（以卵巢寬度標準化），觀(guān)察到豬和人類(lèi)早期卵母細胞生成過(guò)程均顯示出生殖細胞的皮質(zhì)-髓質(zhì)梯度定位。這些數據不僅揭示了豬和人類(lèi)早期卵母細胞生成過(guò)程中保守的基因表達程序和表觀(guān)遺傳重編程模式，還強調了兩種物種在整個(gè)早期卵母細胞生成階段空間動(dòng)態(tài)的保守性。

5.基于空間轉錄組學(xué)（ST）技術(shù)解析豬卵巢發(fā)育過(guò)程中兩種顆粒細胞譜系的時(shí)空特性

在發(fā)育中的卵巢中，顆粒細胞與生殖細胞相互作用，形成了卵巢的基本功能單位，即卵泡。因此，研究人員接下來(lái)提取了顆粒細胞譜系，并使用UMAP進(jìn)行了細胞聚類(lèi)分析，共鑒定出10個(gè)細胞簇。分析了潛在時(shí)間基因表達動(dòng)態(tài)和顆粒細胞前體發(fā)育軌跡中wave II型顆粒細胞（PreGC_II）的命運決定。與PreGC_I相比，簇6和7中的基因在PreGC_II命運決定中起關(guān)鍵作用，如WNT6、GREM1和SFRP4，這些結果共同強調了WNT信號分子在豬顆粒細胞前體命運決定中的保守作用。此外，ALDH1A2的上調表明，醛脫氫酶家族成員1A2可能在wave II顆粒細胞發(fā)育中發(fā)揮重要作用。

使用空間轉錄組學(xué)（ST）技術(shù)對wave I型顆粒細胞前體（PreGC_I）和wave II型顆粒細胞前體（PreGC_II）的標志性基因進(jìn)行了空間定位。結果發(fā)現，LHX9（PreGC_I的標志）和GREM1（PreGC_II的標志）的信號從E45到E75在卵巢中呈現特定的空間分布模式。在E65和E75，GREM1的信號明顯環(huán)繞著(zhù)DDX4（生殖細胞的標志）的信號，表明卵巢卵泡組裝的早期階段，wave II型顆粒細胞（來(lái)源于皮質(zhì)區域的顆粒細胞前體）在卵巢卵泡組裝過(guò)程中發(fā)揮了關(guān)鍵作用。

接下來(lái)，研究人員分析了顆粒細胞譜系的空間特性，進(jìn)行了細胞類(lèi)型反卷積，并觀(guān)察到卵巢表面上皮（OSE）細胞群主要位于豬卵巢的表面，并隨著(zhù)發(fā)育的進(jìn)行細胞數量增加。在E45時(shí)，I型顆粒細胞前體（PreGC_I）細胞廣泛分布在髓質(zhì)中，部分位于皮質(zhì)中；在E65時(shí)，觀(guān)察到I型顆粒細胞前體（PreGC_I）和II型顆粒細胞前體（PreGC_II）分別呈現出特定的皮質(zhì)和髓質(zhì)分布模式。通過(guò)免疫熒光分析進(jìn)一步驗證了這些發(fā)現，上述結果共同揭示了豬卵巢發(fā)育過(guò)程中顆粒細胞譜系的空間時(shí)間特性。

6.豬卵巢發(fā)育期間主要體細胞類(lèi)型的空間特性

除了生殖細胞和顆粒細胞外，研究人員還進(jìn)一步以更高分辨率分析了GI（間質(zhì)細胞）和Sm（平滑肌細胞）這兩種體細胞在豬卵巢中的空間分布模式，因為這兩種細胞類(lèi)型在豬卵巢中顯示出特征性分布，表明它們可能在塑造卵巢微環(huán)境方面發(fā)揮作用，這對早期卵子發(fā)生至關(guān)重要。

將GI和Sm的標記基因整合到空間轉錄組學(xué)中，觀(guān)察到這些體細胞的空間位置隨著(zhù)從E45到E75的發(fā)育進(jìn)展而發(fā)生變化。通過(guò)RNA速度分析確定的分化譜系細胞和增殖狀態(tài)細胞在E45和E55時(shí)的卵巢中隨機分布，而在E65時(shí)，觀(guān)察到這些細胞在卵巢皮質(zhì)中呈現出明顯的空間聚集模式，特別是在生殖細胞附近。這些發(fā)現共同揭示了豬卵巢發(fā)育早期體細胞空間位置模式的動(dòng)態(tài)變化，表明它們在形成對早期卵子發(fā)生至關(guān)重要的卵巢微環(huán)境中發(fā)揮潛在作用。

7.空間共定位分析揭示了皮質(zhì)區和髓質(zhì)區存在截然不同的微環(huán)境，強調了其在調控生殖細胞命運中的關(guān)鍵作用

在闡明了發(fā)育中的豬卵巢中主要細胞類(lèi)型的時(shí)空特征后，研究人員接下來(lái)研究了空間微環(huán)境的變化如何影響發(fā)育中卵巢的細胞命運決定。首先對細胞類(lèi)型豐度進(jìn)行了非負矩陣分解（NMF），以研究發(fā)育中卵巢中細胞的空間共定位。NMF分解的應用成功地表征了位于皮質(zhì)和髓質(zhì)的生殖細胞niche，并揭示了與不同空間來(lái)源的生殖細胞共定位的體細胞。為了全面理解細胞-細胞通信模式的功能性，研究人員使用CellphoneDB進(jìn)行GO富集分析，識別了不同區域生殖細胞與體細胞之間的配體-受體（L-R）對，發(fā)現皮質(zhì)和髓質(zhì)區域之間細胞-細胞通信模式的功能富集結果不同。

此外，研究人員觀(guān)察到皮質(zhì)區域的L-R對顯著(zhù)富集了NOTCH信號通路，PreGC_ii0和PreGC_ii8表達NOTCH信號介導因子NOTCH2，而其配體DLK1和JAG1則在向減數分裂進(jìn)程中的生殖細胞中表達，結果表明NOTCH2在豬早期生殖細胞命運決定中發(fā)揮著(zhù)關(guān)鍵作用。在髓質(zhì)區域細胞-細胞通信分析揭示了在卵細胞階段的生殖細胞表達了參與原始卵泡組裝的配體，包括BMP4和BDNF，而它們相應的受體，包括BMPR1A、BMPR1B和SORT1，則由PreGC i4和PreGC i7表達。研究人員發(fā)現髓質(zhì)區域的體細胞表達了一系列細胞外基質(zhì)（ECM）蛋白，而皮質(zhì)區域則不是這樣，表明ECM蛋白在調節豬卵母細胞生成中發(fā)揮著(zhù)關(guān)鍵作用。

為了進(jìn)一步驗證細胞通訊分析，研究人員接下來(lái)探討了在使用從E55卵巢皮質(zhì)和髓質(zhì)中分離的卵巢組織進(jìn)行體外條件下，補充N(xiāo)OTCH信號抑制劑DAPT和基質(zhì)金屬蛋白酶抑制劑多西環(huán)素是否會(huì )影響生殖細胞命運。從E55胎兒卵巢中解剖出皮質(zhì)和髓質(zhì)組織，將這些組織進(jìn)行體外培養，并在用DAPT處理10天后觀(guān)察到皮質(zhì)卵巢組織中減數分裂標志物STRA8的表達水平顯著(zhù)降低，而髓質(zhì)卵巢組織則沒(méi)有這種情況。相反，在用基質(zhì)金屬蛋白酶抑制劑多西環(huán)素處理10天后，觀(guān)察到卵泡形成標志物L(fēng)HX8的表達水平降低，而在皮質(zhì)組織中未觀(guān)察到可比的影響。綜上所述，這些數據進(jìn)一步支持了在空間背景下的細胞-細胞通訊分析，并強調了卵巢微環(huán)境在調節生殖細胞命運中的重要作用。

 

研究總結

綜上所述，該研究基于單細胞和空間轉錄組測序，深入探討了豬早期卵子發(fā)生過(guò)程中的時(shí)空基因表達譜和卵巢微環(huán)境的空間組織，以及這些特征在人類(lèi)中的保守性。這些發(fā)現不僅豐富了我們對豬繁殖性狀形成復雜機制的認識，還為生殖醫學(xué)領(lǐng)域的研究提供了新的視角和方法。