《Biomolecules》：Comparative ChIP-Seq Analysis of ERα Binding Induced by Estradiol, Zearalenone, and Bisphenol A in Human Ovarian Cells
编辑推荐：
专业术语翻译
外源性雌激素（xenoestrogens）是一类通过结合雌激素受体（estrogen receptors，ERs）模拟内源性雌激素效应的合成或天然化合物。全球范围内人类对外源性雌激素的暴露水平日益增加，引发了关于其内分泌干扰效应的重大关切。本研究
专业术语翻译
外源性雌激素（xenoestrogens）是一类通过结合雌激素受体（estrogen receptors，ERs）模拟内源性雌激素效应的合成或天然化合物。全球范围内人类对外源性雌激素的暴露水平日益增加，引发了关于其内分泌干扰效应的重大关切。本研究利用染色质免疫沉淀测序（ChIP-seq）整合转录组分析，研究了人卵巢细胞在暴露于雌二醇（estradiol，E2）、玉米赤霉烯酮（zearalenone，ZEA）和双酚A（bisphenol A，BPA）后由ERα介导的基因组及转录水平响应。研究结果表明，ZEA诱导的ERα依赖性响应与E2高度相似，而BPA则引发显著减弱的响应。E2或ZEA处理细胞中的ERα信号强度高于BPA处理细胞，表明BPA暴露导致ERα占据率降低。研究人员鉴定了641个在任何处理条件下均被ERα占据的基因组区域，定义了一组核心共享结合位点。与mRNA-seq数据集整合分析显示，162个ERα结合位点位于109个上调基因转录起始位点（transcription start sites，TSS）±50 kb范围内。其中，28个基因在所有处理条件下均被诱导表达，构成人卵巢细胞中的核心雌激素反应基因集。对六个代表性基因（GREB1、RASGRP1、CISH、RBBP8、KCNF1和LONRF2）的RT-qPCR验证证实，E2和ZEA暴露比BPA暴露诱导更强且更持久的转录激活；这种诱导可被ERα拮抗剂MPP抑制。研究结果表明，BPA激活ERα染色质结合及下游转录程序的能力较低，主要影响人卵巢细胞中高亲和力的雌激素反应位点。
论文解读文章
一、研究背景与科学问题
外源性雌激素（xenoestrogens）是环境中广泛存在的一类内分泌干扰化合物，涵盖天然来源的植物雌激素和真菌毒素，以及人工合成的工业化学品、药物和个人护理产品成分。这类化合物通过结合雌激素受体（estrogen receptors，ERs）模拟内源性雌激素的生物学效应，进而干扰内分泌系统的正常功能。双酚A（bisphenol A，BPA）是合成外源性雌激素的典型代表，广泛用于聚碳酸酯塑料、环氧树脂、牙科密封剂、医疗器械以及食品和饮料容器的内涂层制造。BPA易于从塑料制品中浸出并污染储存的食品和液体，人体可通过胃肠道吸收、吸入和皮肤接触等多种途径暴露于BPA，且已在外周血、尿液和羊水等多种人体生物样本中被普遍检出。玉米赤霉烯酮（zearalenone，ZEA）是由镰刀菌属（Fusarium species）真菌产生的霉菌毒素，兼具真菌雌激素（mycoestrogen）的属性。ZEA化学稳定性高，能够耐受食品加工过程，导致面包、麦芽、面粉等谷物制品频繁受到污染，同时还可通过肉、奶等动物源性食品间接进入人体。
尽管BPA和ZEA的人体暴露水平因塑料污染加剧和气候变化促进真菌生长而持续上升，目前科学文献对这类化合物分子作用机制的认识仍然有限，特别是在基因组事件层面。现有研究多基于动物模型，基于人类组织或细胞的研究证据相对匮乏。在此背景下，研究人员利用人卵巢细胞模型，系统比较了生理性雌激素E2与两种代表性外源性雌激素ZEA和BPA在ERα介导的基因组结合和转录调控层面的异同，旨在阐明外源性雌激素发挥雌激素效应的基因组机制。该研究基于研究人员前期工作——已证明ZEA和BPA处理ERα阳性卵巢癌细胞系可显著促进细胞增殖和迁移，并通过转录组分析得到进一步验证——现进一步从全基因组水平揭示其分子机制。
二、研究设计与核心技术方法
研究人员选用PEO1人卵巢上皮细胞系作为实验模型，该细胞系购自Merck公司（德国达姆施塔特；ECACC，英国索尔兹伯里），其ERα表达及雌激素敏感性已在前期研究中得到确证。细胞在含10%胎牛血清（fetal bovine serum，FBS）、1% L-谷氨酰胺、100 μg/mL链霉素和100 U/mL青霉素的RPMI1640培养基中、37°C、90%湿度和5% CO?条件下培养。为降低培养基中酚红和雌激素的干扰效应，在药物处理前将培养基更换为含5%葡聚糖-活性炭处理胎牛血清（dextran-coated charcoal-treated FBS，DCC-FBS）的无酚红RPMI1640培养基；在ChIP-seq分析中该饥饿处理时间延长至72 h以增强ERα对雌激素处理的敏感性。随后分别以E2（10 nM）、ZEA（10 nM）和BPA（100 nM）处理细胞，溶剂处理组作为对照。
研究采用的核心技术方法是染色质免疫沉淀测序（chromatin immunoprecipitation followed by sequencing，ChIP-seq），使用抗ERα抗体（Sigma，目录号06-935）进行免疫沉淀，DNA纯化后构建测序文库，在Illumina NextSeq 500平台完成75 bp单端测序，所有实验条件均设两个生物学重复。测序数据通过MACS2进行峰值（peak）识别，参数设置为q值截断值0.001并启用亚峰解析（--call-summits），使用ENCODE黑名单去除人工制品。每样本获得1200万至1600万条测序读数，超过96%的读数可比对至人类参考基因组GRCh38.p13。ChIP-seq数据与研究人员前期生成的RNA-seq转录组数据（相同实验条件下的结果）进行整合分析。基因表达验证采用逆转录定量PCR（RT-qPCR），以GAPDH为内参基因，采用2?ΔΔCt法计算相对表达量。此外，研究还运用了HOMER进行从头基序发现（de novo motif discovery），利用Network Analyst平台基于STRING相互作用组和Reactome数据库进行蛋白质-蛋白质相互作用（protein-protein interaction，PPI）网络分析和功能富集分析。
三、主要研究结果
3.1 利用ChIP-Seq鉴定E2、ZEA和BPA激活细胞中的ERα结合区域
研究人员对PEO1细胞分别给予溶剂、E2、ZEA和BPA处理1小时后进行ChIP-seq分析。全基因组峰值鉴定显示，E2和ZEA处理样本中检测到的ERα结合峰数量显著多于BPA处理样本。各配体处理组的重复间共有峰值（consensus peak set）分别为：对照组63个、E2处理组2214个、ZEA处理组3200个和BPA处理组656个。从头基序富集分析表明，经典雌激素反应元件（estrogen response element，ERE）基序GGTCAnnnTGACC在所有处理条件下（除对照组外）均被识别为排名第一的富集基序，证实了ERα特异性DNA结合的序列基础。
3.2 E2和ZEA比BPA引发更强健的ERα介导反应
比较分析显示，三种配体处理存在大量重叠的ERα结合区域：641个ERα结合区域在所有三种处理条件下均被检测到，构成核心共享结合位点集。值得注意的是，98%的BPA响应峰在E2和/或ZEA处理中同样存在，仅有7个峰为BPA处理所独有。此外还鉴定出E2与ZEA共有的1450个峰、E2独有的123个峰和ZEA独有的1105个峰。信号强度分析揭示了两个关键发现：其一，在三种处理共有的结合位点上，E2或ZEA处理细胞的ERα信号强度显著高于BPA处理细胞，表明BPA在共享位点的ERα结合能力较弱；其二，同一处理条件下，核心共享峰集的ERα信号强度普遍高于各处理特异峰集，提示这些共享调控区域的ERα结合更为强健和稳定。对照组中也可检测到低强度的ERα峰，这与已有报道中ERα可在配体刺激前预先招募至特定基因组位点的发现一致。
3.3 ChIP-Seq数据与转录组谱的相关性分析
研究人员将ChIP-seq数据与前期在相同实验条件下获得的RNA-seq数据（鉴定了365个受任一配体显著上调的基因，log2FC≥1）进行整合分析，在上调基因转录起始位点（transcription start site，TSS）±50 kb窗口内鉴定出162个ERα结合区域，涉及109个上调基因，其中17个ERα结合区域同时与两个基因相关联，表明ERα参与E2、ZEA和/或BPA对这109个基因的诱导表达。信号强度分析显示，这109个基因的表达水平在E2处理样本中最高、ZEA处理样本与之相当、BPA处理后较低，ERα结合信号和基因诱导倍数（fold-change，FC）均呈现一致趋势。进一步分析表明，由三种配体共同诱导的基因在E2或ZEA处理下的倍数变化显著高于仅由E2和ZEA共同诱导或仅由单一配体诱导的基因。通过RNA-seq与ChIP-seq数据整合，最终鉴定出28个同时受E2、ZEA和BPA诱导且具有ERα结合支持的基因，构成人卵巢细胞的核心雌激素反应基因集，其中AGT、DHRS2、GREB1、KCNF1、LONRF2、PDGFRL、PIPOX、RASGRP1、RBBP8和SLC7A2的TSS±50 kb区域内存在一个以上ERα峰。
3.4 E2和ZEA比BPA诱导更强的转录反应且随时间保持稳定并可被ERα拮抗剂抑制
研究人员从28个核心基因中选取六个代表性基因（GREB1、RASGRP1、CISH、RBBP8、KCNF1和LONRF2）进行RT-qPCR验证，在处理后2、6和24小时检测基因表达动态。结果显示，三种化合物均能快速诱导这六个基因的表达，但E2和ZEA诱导的基因表达水平显著高于BPA，且这种差异在后续时间点保持不变，表明BPA较弱的转录诱导能力具有时间持续性。作为阴性对照的ACTRT3基因（已知不参与ERα介导的雌激素信号传导）在E2、ZEA或BPA处理后均未表现出诱导表达，验证了实验体系的特异性。进一步使用选择性ERα拮抗剂MPP（methyl-piperidinopyrazole）进行预处理，发现1 μM MPP能显著抑制E2、ZEA和BPA对这六个基因的诱导表达，证实这些转录反应确实由ERα介导。
3.5 雌激素反应基因在细胞周期调控中高度互联
研究人员利用Network Analyst平台基于STRING相互作用组数据库对109个基因进行PPI网络分析，构建的网络包含395个节点、456条边和23个种子基因，呈现广泛的相互作用。网络分析鉴定出NOL6、FOS和NOD2为中枢（hub）基因，它们具有最高的节点度（node degree）和介数中心性（betweenness centrality），提示这些基因在调控网络中具有重要地位。基于Reactome数据库的功能富集分析显示，这109个基因显著富集于细胞周期调控相关通路，提示其可能参与细胞增殖和细胞周期控制机制。对28个核心基因的基因集富集分析则发现，氨基酸转运相关通路是最显著富集的生物学过程，雌激素反应性基因表达程序也位列富集通路前列。
四、讨论与结论
本研究表明，ZEA诱导的ERα基因组结合效率和转录反应与生理性E2相当，与前期研究中ZEA诱导细胞增殖和迁移能力与E2相当的结果一致，进一步支持ZEA作为强效雌激素干扰物的定位。鉴于ZEA对谷物食品的污染广泛且水平较高，并已在人体血清、卵泡液、尿液、子宫内膜癌组织及结直肠癌患者血液样本中被检出，其健康风险值得高度关注。相比之下，BPA诱导的ERα介导反应在量级和ERα基因组占据程度上均显著减弱，且其较弱的转录诱导能力随时间保持稳定，这与BPA作为弱雌激素的既有认识一致。BPA在ERα结合信号最强的位点仍能诱导可检测的转录反应，提示BPA可能优先增强高雌激素反应性基因的表达。BPA雌激素效力较低的可能机制包括：BPA对ERα和ERβ的结合亲和力显著低于E2；BPA与ER配体结合域的相互作用方式不同导致受体激活效率下降；BPA-ERα复合物的转录共激活因子招募和核小体定位发生改变从而限制下游转录激活。
研究人员整合ChIP-seq和RNA-seq数据鉴定了109个可能受ERα调控的雌激素反应基因，其中28个为高雌激素反应性核心基因。部分基因（如GREB1、RBBP8、SLC7A2等）在BPA处理乳腺癌细胞系中的诱导已被报道，GREB1、NPY1R、PLAT等基因的活性与乳腺癌ER阳性状态和/或内分泌治疗耐药相关。功能富集分析揭示细胞周期调控、氨基酸转运和核糖体生物发生通路的显著富集，与雌激素和外源性雌激素促进细胞增殖的已知效应一致，提示氨基酸转运和核糖体生物发生的上调可能为增强的蛋白质合成能力提供支持，从而共同构成雌激素相关肿瘤发生发展的潜在机制。
研究结论：ZEA诱导ERα结合基因组区域的效率与生理性配体E2相当；BPA激活ERα依赖性转录程序的能力较弱，仅最高度雌激素反应性的位点对该外源性雌激素表现出可检测的转录反应。这是首次在人类卵巢细胞中表征BPA基因组效应、在任何实验模型系统中表征ZEA基因组效应的研究。研究同时鉴定出28个能被ZEA和BPA持续诱导的高雌激素反应性基因。需要指出的是，ERα介导的雌激素信号传导具有高度的细胞类型依赖性，本研究仅限于单一卵巢细胞系，研究结果的普适性仍需在其他卵巢细胞类型中进一步验证；系统研究ERα辅因子相互作用、ERα占据动态和染色质可及性变化将提供重要的机制见解，是未来研究的重要方向。
