动态染色质变化与基因调控网络揭示鸭卵泡发生的多组学机制

《Journal of Animal Science and Biotechnology》：Comprehensive multi-omics reveals dynamic chromatin changes and gene regulatory networks during duck folliculogenesis

　　

摘要 卵泡发育是脊椎动物繁殖的先决条件，受复杂的基因组构象和调控元件精确调控。然而，禽类卵泡发育过程中颗粒细胞（Granulosa Cells, GCs）的三维基因组与调控元件相互作用的动态变化仍不清楚。在此，研究人员整合了北京鸭（Anas platyrhynchos domestica）GCs在7个发育阶段的RNA测序（RNA-seq）、ATAC测序（ATAC-seq）、CUT&Tag和Hi-C数据，构建了卵泡发育的高分辨率三维顺式调控图谱，揭示了禽类卵泡发生的染色质动力学基础。整合分析表明，H3K27ac的动态变化而非单纯的染色质可及性，与卵泡选择和成熟阶段特异性的转录增加密切相关。研究人员鉴定了与关键卵泡基因表达显著相关的增强子和超级增强子（Super-enhancers, SEs）。在三维基因组组织方面，拓扑关联结构域（Topologically Associating Domains, TADs）保持高度稳定，作为结构支架；然而，阶段特异性的边界变化与关键调节基因的转录改变相一致。此外，研究人员推断了一个包含46个核心转录因子（Transcription Factors, TFs）的推定基因调控网络（Gene Regulatory Networks, GRNs），预测其与卵泡发育密切相关。最后，比较分析强调了这些调控元件在鸟类和哺乳动物之间的保守性和物种特异性。该研究提供了一个整合的多组学资源，为鸭卵泡发育的表观基因组景观提供了新颖的见解，所得到的数据集和调控图谱为进一步研究卵泡发生的机制和理解跨物种调控差异奠定了宝贵的基础。

论文解读：鸭卵泡发生过程中染色质动态与基因调控网络的整合解析

研究背景与意义

卵泡发育是脊椎动物生殖生物学的关键环节，其精确调控依赖于复杂的基因组三维结构与表观遗传机制。尽管哺乳动物卵泡发育的表观调控已有较多研究，但禽类作为卵生动物代表，其卵巢形态、卵泡层级结构及激素分泌模式与哺乳动物存在显著差异，且禽类卵泡选择（follicle selection）与排卵过程中的染色质动态互作机制尚不明确。现有研究多局限于单一组学层面，缺乏对三维基因组与调控元件协同作用的系统性解析。为此，本研究通过多组学整合策略，构建鸭卵泡发育的三维顺式调控图谱，旨在揭示禽类卵泡发生的表观遗传基础，相关成果发表于《Journal of Animal Science and Biotechnology》。

关键技术方法

研究人员选取高产北京鸭（产蛋率>90%）的颗粒细胞（GCs），涵盖小白色卵泡（SWFs）、大白色卵泡（LWFs）、小黄色卵泡（SYFs）、F5/F3/F1等级前卵泡及排卵后卵泡（POFs）共7个阶段，每组设置生物学重复。采用RNA测序（RNA-seq）、ATAC测序（ATAC-seq）、组蛋白修饰CUT&Tag（H3K4me3、H3K4me1、H3K27ac、H3K27me3）及Hi-C技术构建多组学数据集，结合ChromHMM注释染色质状态，ROSE算法鉴定超级增强子（SEs）与超级沉默子（SSs），并通过整合三维染色质互作与基序分析构建基因调控网络（GRN）。

研究结果

1. 鸭卵泡发育过程中GCs的转录组与表观基因组图谱

通过对7个阶段GCs的多组学分析，研究人员鉴定了8,137个差异表达基因（DEGs），其中SYF至F5（卵泡选择期）和F1至POF（排卵期）的DEG比例显著高于其他阶段。K-means聚类将DEGs分为8个簇，功能富集显示簇4、5、6在层级卵泡中高表达，主要参与类固醇合成、Wnt信号通路等关键生物学过程。表观数据显示，启动子与TSS近端转录区域在转录起始位点（TSS）处富集最高，且染色质状态在相邻阶段间存在7%-25%的动态变化。

2. 卵泡发育过程中的特异性染色质状态变化

利用ChromHMM整合四种组蛋白修饰与染色质可及性，定义了15种染色质状态，包括启动子（TssA等）、增强子（EnhA等）、转录区域（TxFlnk等）及沉默区域（Repr等）。研究发现，染色质可及性与H3K4me3在卵泡选择后显著升高，H3K27ac与H3K27me3则在F1阶段（排卵前）急剧上升，而H3K4me1整体变化较小，揭示了不同组蛋白修饰在鸭卵泡发育中的特异性动态模式。

3. H3K27ac变化与GCs基因表达动态的关联

在卵泡选择与成熟阶段，上调基因伴随H3K27ac与H3K4me1水平升高，下调基因则呈相反趋势，而染色质可及性与H3K4me3在上下调基因中均表现为升高。SYF至F5阶段鉴定到3,313个差异H3K27ac峰，其中875个基因（FC=1.56）同时显示染色质可及性增加，包括ESR1、BMPR1B、WNT4等已知卵泡选择相关基因，富集于激素调节与脂质代谢通路。F1阶段H3K27ac上调基因则富集于胶原重塑与类固醇合成通路，表明H3K27ac协同染色质开放驱动阶段特异性基因表达。

4. 增强子与超级增强子（SEs）与基因表达的关联

通过ROSE算法鉴定到平均964个SEs与1,040个SSs，SE靶基因表达水平显著高于普通增强子（TEs）与沉默子（TSs/SSs）。阶段特异性SEs的靶基因聚类显示，Cluster 3与Wnt、PPAR信号通路相关，SYF特异性SEs关联AKT1、INHBB等基因，F1特异性SEs则关联INHBA、IGF1R等排卵相关基因，表明SEs通过调控关键基因参与卵泡选择与成熟。

5. 三维结构变化在转录调控中的调节作用

Hi-C分析显示，44%的基因组区室（A/B compartments）保持稳定，49%-52%位于活跃A区室，但存在25%的A→B转换与22%的B→A转换。拓扑关联结构域（TADs）在发育过程中高度保守（97.9%稳定），平均大小约296 kb。尽管TAD边界总体稳定，但SYF与F5、F5与F1阶段分别存在90个与60个特异性边界，其邻近基因富集于MAPK、PI3K等卵泡发育通路。代表性基因ZEB2在F1/F5阶段获得TAD边界并伴随表达升高，提示三维基因组通过局部精细调控而非全局重排影响转录。

6. 整合多组学揭示多层顺式调控景观

通过整合转录因子（TF）基序富集与增强子-基因关联，构建包含8,132个增强子、232个TFs及3,376个基因的GRN，筛选出46个核心TFs。其中卵泡选择相关TFs（如FOXP1、FOXO1）在SYF阶段高表达，排卵相关TFs（如PPARα、NR5A2）在F1阶段富集。结合GWAS数据，发现27个与产蛋性状相关的位点与GCs增强子重叠，例如rs7:13532237可能通过调控CYP17A1影响激素合成。进一步揭示三种顺式调控机制：近端增强子（平均距离启动子82,248 bp）、远端染色质环介导的增强子（平均距离359,411 bp）及TAD边界重组。

7. 鸟类与哺乳动物卵泡发育的比较

跨物种分析显示，30.44%的鸭增强子在鸡中序列保守，其中31.54%具功能保守性，靶基因表达模式相似。7个TFs（GATA3、GATA6、FOXK1、FOXK2、FOXO3、TRPS1、ZFX）在鸭、鸡、人、小鼠、猪中保守，其结合基序密度在不同物种间存在差异，提示调控元件的进化速率与物种特异性适应相关。

讨论与结论

本研究首次系统构建了鸭卵泡发育的三维顺式调控图谱，阐明H3K27ac动态与染色质可及性协同驱动阶段特异性转录的机制，鉴定了46个核心TFs及关键增强子元件，并揭示TAD结构稳定性与局部边界重组的调控平衡。跨物种比较发现禽类与哺乳动物在卵泡发育调控中存在保守通路与物种特异性适应，为理解脊椎动物生殖进化的表观遗传基础提供了新视角。研究结果为禽类繁殖性能改良及卵泡发育异常机制解析提供了重要的多组学资源。

结论

综上所述，本研究首次提供了鸭卵泡发育的三维顺式调控图谱，鉴定了禽类卵泡发生过程中的关键增强子、转录因子及相互作用调控网络，并揭示了鸟类与哺乳动物在卵泡发育调控中的保守转录因子。这些发现为理解颗粒细胞支持卵泡发育的分子机制提供了宝贵的资源。