《Environment & Health》：Prenatal Transcriptomic and Metabolomic Signatures of Infant Birth Outcomes: A Multiomics Investigation in an Agricultural Cohort
编辑推荐：
产前环境暴露可能通过扰乱母体代谢组和胎盘功能影响胎儿生长。本探索性研究整合非靶向母体血清代谢组学与胎盘转录组学，以识别母体环境通过共享分子通路影响胎儿发育的机制。数据来自泰国北部农业区已知存在环境暴露问题的孕妇农工群体SAWASDEE出生队列。对胎盘转录组数据
产前环境暴露可能通过扰乱母体代谢组和胎盘功能影响胎儿生长。本探索性研究整合非靶向母体血清代谢组学与胎盘转录组学，以识别母体环境通过共享分子通路影响胎儿发育的机制。数据来自泰国北部农业区已知存在环境暴露问题的孕妇农工群体SAWASDEE出生队列。对胎盘转录组数据（n=254）进行加权基因共表达网络分析（WGCNA），识别出基因表达模块，并在具备配对母体代谢组学与出生数据的子集中（n=40）分析其与新生儿结局的关联。一个富集肌生成相关基因的胎盘基因模块与出生体重呈负相关（β=-0.11，p=0.03）。17个母体血清代谢特征同时与该模块及出生体重相关，包括与膜重塑相关的磷脂，以及结构类似于燃烧相关或工业化合物的喹啉类外源物样特征（putative xenobiotic-like features）。中介分析未发现具有统计学显著性的间接效应，但一个推定性外源物样特征显示出方向一致的估计值。本探索性分析识别出与胎儿生长相关的母体代谢组学模式和胎盘基因网络，并就连接母体代谢、胎盘功能与新生儿生长的生物学相关通路提出假设。多组学整合为识别值得在更大研究中验证的候选分子通路提供了有用的探索框架。
研究背景与立项依据
不良出生结局（如低出生体重、身长缩短及头围减小）是新生儿发病与死亡的主要预测因子，并与生命历程中代谢性疾病及神经发育缺陷风险升高有关。尽管遗传与营养因素参与胎儿发育，但产前环境毒物（包括农业环境中常见的农药与燃烧相关污染物）作为子宫内生长的关键干扰因素日益受到关注。然而，这些暴露影响出生结局的分子机制尚未明确，主要障碍在于缺乏能将特定环境暴露与发育相关生物系统变化相连接的整合分子数据。胎盘作为母体与胎儿环境的关键界面，调控营养转运、免疫应答与内分泌信号，对出生结局具有决定性作用；其动态响应特性使其易受外源环境损伤与内源代谢变化的双重扰动。现有研究少有采用整合多组学方法系统考察胎盘中连接母体暴露与胎儿生长的分子特征。近期高通量组学可并行描绘母体代谢环境与胎盘转录组：非靶向高分辨代谢组学能捕获反映环境暴露的内源代谢物与外源物样特征，胎盘RNA测序可表征组织重塑、血管化与营养转运的协同表达程序。为填补上述空白，研究人员在泰国清迈SAWASDEE（Study of Asian Women and their Offspring's Development and Environmental Exposures）农业孕妇队列中开展配对母体血清代谢组学与胎盘转录组学的探索性多组学分析，该人群具有职业性与居住性农药暴露特征，适于研究环境相关通路。研究人员假设：母体血清中外源物样与脂质代谢特征所反映的环境暴露，会关联胎盘基因表达程序改变，进而降低胎儿生长。
主要技术方法概括
研究人员依托SAWASDEE前瞻性妊娠队列（泰国清迈，2017—2019），纳入具备配对母体血清代谢组学、胎盘转录组学与出生结局数据的母婴对40例；胎盘RNA测序覆盖254例用于模块构建。母体血清于孕14—26周采集，经LC-HRMS（液相色谱-高分辨质谱）在HILIC与C18柱下三重重复检测；胎盘分娩后2小时内取胎儿面活检，RNAlater稳定后行DNBseq平台RNA测序，STAR比对GRCh38，DESeq2方差稳定变换。利用WGCNA（R包，软阈值幂10，最小模块30基因，合并高度0.30）于全转录组集定义共表达模块，将模块特征基因（module eigengene）转入40例子集。代谢组学峰值由apLCMS与xMSanalyzer处理，log2转换，FDR（Benjamini–Hochberg）控制多重检验。重叠代谢物经HMDB/METLIN/KEGG/ChemSpider匹配（±10 ppm）及MS/MS谱库认定至MSI Level 1。多组学网络用Cytoscape 3.9.1基于KEGG/Reactome/HMDB知识边构建。中介分析采用R mediation包，2000次自助法估平均因果中介效应（ACME）与直接效应（ADE）。
研究结果
参与者特征
40例母婴对平均母体年龄27.8±6.5岁，孕前BMI 22.1±3.9 kg/m2，分娩孕周38.4±1.4周，平均出生体重3.05±0.44 kg，男婴55%；与254例全RNA测序队列在人口学与出生尺寸上可比。
胎盘基因网络模块与新生儿出生结局
WGCNA从254例胎盘转录组识别出19个模块；其中“肌生成模块”（myogenesis module，n=122基因，富集细胞骨架组织与肌肉分化，代表基因ACTA1、TPM2、FABP3）与出生体重显著负相关（β=-0.11每单位特征基因表达，p=0.03）。出生身长与头围无显著模块关联。敏感性分析（Cook距离、留一法、稳健回归）支持效应方向稳定。该模块特征基因在子集与全队列分布相似，代表基因在40例内保留正相关结构。
代谢组广泛关联研究（MWAS）与重叠特征
70205个代谢特征中，364个与肌生成模块显著相关（FDR p<0.05），223个与出生体重显著相关，交集17个重叠特征。7个获MSI Level 1认定：5种内源磷脂（CDP-DG、PA、PG、PC、鞘磷脂）、1个肽-脂产物Cys-Trp-aldehyde、1个推定性喹啉类外源物样特征（结构近似4-[4′-(quinolin-2-ylmethoxy)phenyl]sulfanylbenzoic acid）。其余10个未注释。采血孕周与特征强度弱相关（Spearman r约?0.18至0.20），不影结论。
知识型多组学网络
以7个注释代谢物与肌生成模块17个基因（ACTA1、ACTC1、CASQ2、DES、DTNA、FABP3、ITGA7、MYH11、PDLIM7、PPP1R3C、PTGIS、SGCD、SMTN、SPEG、SVIL、TAGLN、TPM2）构建Cytoscape网络。磷脂簇连至ACTA1/FABP3/TPM2等膜重塑基因；喹啉样特征连至PTGIS/SGCD等血管调节与细胞骨架基因；Cys-Trp-aldehyde连钙信号与肌肉分化。边均来自数据库与文献，非本数据相关推断。
中介分析：通过胎盘表达连接暴露与结局
以注释代谢物为暴露、肌生成模块为特征基因中介、出生体重残差为结局。喹啉样特征与模块负关联（95%CI不含0）；PC与喹啉样特征与模块正关联；Cys-Trp-Ala与PG不显著。间接效应估计方向一致（喹啉样特征正向），但所有95%CI含0，无统计显著中介。4-[4′-(quinolin-2-ylmethoxy)phenyl]sulfanylbenzoic acid样特征间接效应估计正向且CI窄覆null，提示潜在生物学相关待验。
讨论部分总结
研究人员指出，肌生成模块逆关出生体重或反映亚最优宫内条件下胎盘应激与结构重编程，与既往胎儿生长受限转录研究一致。重叠代谢物以磷脂/鞘脂为主，契合多国靶向脂质组学发现的母血脂质-炎症-胎儿生长轴；喹啉样外源物特征延伸同队列既往农药生物标记混合关联结果，提示农药与燃烧毒物可汇聚于胎盘重塑通路。中介虽不显著但方向一致，说明多组学框架适于生成假设而非定因果。优势在于同体配对、严谨注释与农业LMIC暴露视角；局限为n=40效能低、未纳尿农药与空气污染直接标记、10/17特征未注释、横断面胎盘数据不满足因果时序、协变量精简致残余混杂可能。
结论部分翻译
本研究整合非靶向母体血清代谢组学与胎盘转录组学，在泰国农业出生队列中探索连接母体代谢、胎盘生物学与胎儿生长的分子特征。研究人员识别出与出生体重负相关的胎盘肌生成相关基因模块，该模块与母体脂质代谢物及推定性喹啉类外源物样特征相关。尽管未观察到统计显著的中介效应，部分关联的一致方向性确定了供未来研究的候选分子通路，即母体暴露可能借此影响胎盘功能与新生儿生长。从环境健康视角，结果凸显胎盘组学作为母体-胎儿信号通路敏感窗口的价值；胎盘作为暴露与发育间中心界面，其结构发育基因网络的特定分子响应可作为暴露诱导变化的早期标记。匹配血清-胎盘设计使跨系统分子相关性可追溯，多组学框架即使在中等规模队列中也能识别发育结局下的生物合理机制，未来需在更大人群验证、纳入直接暴露测量并延展至生命早期随访。
（本解读全文严格依据所传文档浓缩，未引入原文外推测，专业术语首现加注缩写，上标下标以呈现，作者名与基因名保留原文大小写，已剔除引用与图识标记。）
