《Actuators》：Adaptive Multi-Fault-Tolerant Boundary Control of an Euler–Bernoulli Beam System with Control-Matched Disturbances
编辑推荐：
本研究的目的是表征乳腺微生物组、抗生素抗性基因（antibiotic resistance genes，ARGs）和移动遗传元件（mobile genetic elements，MGEs）在干奶期、分娩和泌乳早期的协同动态。乳腺微生物组在这些阶段经历显著的生态
本研究的目的是表征乳腺微生物组、抗生素抗性基因（antibiotic resistance genes，ARGs）和移动遗传元件（mobile genetic elements，MGEs）在干奶期、分娩和泌乳早期的协同动态。乳腺微生物组在这些阶段经历显著的生态变化，然而微生物组成、ARGs和MGEs的协同动态仍知之甚少。研究人员对干奶前（BM）、分娩后即刻（ACM）和产后一个月（AM）采集的乳腺分泌物样本进行了鸟枪法宏基因组测序。乳腺微生物组呈现出清晰的“暴露–瓶颈–重组”轨迹。BM以高微生物多样性和环境相关类群的富集为特征，而ACM则表现出显著的免疫瓶颈，微生物多样性和网络复杂度明显降低。在AM阶段，微生物群落部分恢复，但仍与BM状态存在差异，表明分娩后持续存在生态重组。ARGs和MGEs表现出平行动态：BM中抗性组和移动组多样性广泛，ACM中急剧收缩，AM中则出现选择性再扩展。网络分析进一步揭示了BM中最大的生态复杂度，ACM中ARGs/MGEs连接性增加，AM中部分稳定。这些发现表明，宿主生理转变（host physiological transitions，HPTs）与干奶牛治疗（dry cow therapy，DCT）共同驱动了干奶期乳腺微生物组、抗性组和移动组的协同重塑。
**论文解读**
**研究背景**
乳腺内感染（Intramammary infection, IMI）是全球奶牛生产中流行最广、经济负担最重的疾病之一。乳房炎不仅损害动物福利，还因高发病率和单例高昂的治疗成本给牧场带来巨大经济损失。临床型乳房炎每例可导致数百美元的直接损失，而亚临床型乳房炎虽无明显临床症状，但通过降低产奶量和升高体细胞计数（Somatic cell count, SCC）造成更大的隐性损失。因此，识别乳房炎发生和控制的关键阶段对于改善牛群健康和牧场盈利至关重要。干奶期被广泛认为是IMI流行病学中的关键窗口期。在这一阶段，乳腺经历从泌乳停止到分娩再到恢复泌乳的深刻生理转变，伴随组织重塑、乳腺微环境改变和局部免疫状态波动。由于挤奶停止后角蛋白栓尚未完全形成、乳头括约肌关闭延迟以及乳腺内压升高，病原体容易侵入乳腺。因此，干奶期不仅是新IMI的高风险阶段，也是实施有效控制策略的关键时机。干奶牛治疗（Dry cow therapy, DCT）是一种广泛采用的管理措施，通过在干奶时向每个乳头灌注长效抗微生物药物来清除已有感染并预防干奶期新发感染。尽管DCT在疾病控制中的有效性已获充分证实，但其作为大规模抗生素干预手段，可能对乳腺微生物生态系统产生显著影响，然而这些影响尚未得到充分表征。近年来，乳腺被越来越多地视为一个动态的微生物生态系统。健康奶牛乳汁中存在多样且动态的微生物群落，有助于维持乳腺稳态。纵向研究显示，从干奶期到泌乳期，乳汁微生物群的多样性和组成发生显著变化，并与乳腺健康状况密切相关。此外，乳腺上皮细胞能够感知微生物刺激并通过免疫相关信号通路进行转录重编程，表明宿主与微生物之间存在双向相互作用。在这一背景下，抗生素抗性基因（Antibiotic resistance genes, ARGs）作为微生物群落的关键功能组分，其生态作用已超越临床抗菌素耐药性本身。ARGs的分布和传播受微生物群落结构和环境选择压力的共同影响。干奶期的生理变化以及DCT相关的抗生素暴露可能共同塑造微生物组和抗性组的重组。然而，针对健康奶牛干奶期微生物群落与ARGs协同动态的纵向研究仍然有限。因此，本研究采用纵向设计，覆盖三个关键生理阶段：干奶前（Before dry-off, BM）、分娩后即刻（After calving, ACM）和产后一个月（One month postpartum, AM）。在干奶时实施DCT后，利用鸟枪法宏基因组测序系统地表征了这些阶段乳腺微生物组和抗性组的组成与动态。通过整合宿主生理变化和抗生素干预，本研究旨在揭示干奶期乳腺生态系统中微生物和抗性组的重组模式，从而为优化乳房炎控制策略和促进奶牛生产中抗生素的合理使用提供见解。该论文发表在《Actuators》（注：原文提供的期刊名为Actuators，但与内容不符，此处按原文输出）。
**主要技术方法**
研究人员于2024年1月至3月在中国北京市通州区某商业奶牛场（牛群规模1000头）随机选取12头荷斯坦健康经产奶牛作为样本队列。在干奶前、分娩后即刻和产后一个月三个时间点采集右侧前乳区的乳腺分泌物样本。干奶时对每个乳头灌注长效头孢喹肟产品和内部乳头封闭剂。样本经离心、裂解和自动化核酸提取后，采用Mag-Bind? Soil DNA Kit提取总DNA，利用Covaris M220片段化后构建双端文库，在Illumina NovaSeq平台进行宏基因组测序。测序数据经Fastp质控、BWA比对去除宿主基因组（Bos taurus ARS-UCD1.2），MEGAHIT从头组装，Prodigal预测开放阅读框，CD-HIT构建非冗余基因集。通过Diamond比对NR、eggNOG、KEGG、CARD和MGEs90数据库进行分类、功能和抗性注释。ARG和MGE丰度以每百万转录本（Transcripts per million, TPM）归一化。采用Shannon指数、主坐标分析（Principal coordinate analysis, PCoA）、ANOSIM、PERMANOVA、Kruskal-Wallis检验、LEfSe分析（LDA > 3.0, p < 0.05）以及基于Spearman相关（|r| ≥ 0.6, FDR调整p < 0.05）的共现网络分析进行统计。
**研究结果**
**3.1. 泌乳阶段微生物群落动态**
通过属级微生物组成比较发现，BM、ACM和AM之间呈现显著的分阶段结构变化。BM包含190种独特物种，而AM拥有最多的阶段特异性物种（1830种），表明分娩后微生物重组产生了新的分化群落而非恢复干奶前状态。仅126种物种在三个阶段共存，BM和ACM仅共享2种，凸显了初乳阶段强烈的生态不连续性。Alpha多样性分析显示ACM的Shannon指数显著下降（p = 6.623×10
-5），符合初乳相关的免疫瓶颈；AM多样性部分恢复但未重回BM水平。PCoA基于Bray-Curtis距离揭示清晰的分阶段分离，BM样本紧密聚集，ACM更紧凑，AM样本分散度增加（ANOSIM R = 0.516, p = 0.001），表明早期泌乳阶段微生物群落发生个性化异质性重组。门水平上，假单胞菌门（Pseudomonadota）在所有阶段占优势但相对丰度波动（BM 66.76%, ACM 41.06%, AM 62.15%）；BM额外富集拟杆菌门（Bacteroidota）和芽孢杆菌门（Bacillota），与环境暴露一致。属水平上，BM中富集环境相关属（如Psychrobacter、Denitrificimonas、Thiopseudomonas），ACM和AM则主要由Luteimonas、Campylobacter、Chlamydia、Methanobrevibacter和Escherichia等主导，表明微生物群落经历了从环境暴露向宿主选择压力的转变。
**3.1.2. 泌乳阶段微生物群落组成比较分析**
显著属级差异表明，与乳房炎相关的多个属在BM组中显著富集：葡萄球菌属（Staphylococcus，3.56% vs ACM 0.0008%和AM 0.00079%，p = 0.00659）、假单胞菌属（Pseudomonas，5.47% vs 0.97%和0.15%，p = 0.00071）、莫拉菌属（Moraxella，0.996% vs 0%和0.0037%，p = 0.00030）和嗜冷杆菌属（Psychrobacter，17.45% vs 2.71%和0.55%，p = 0.00115）。环境相关属Denitrificimonas和Thiopseudomonas也显著富集于BM。ACM组中潜在病原体丰度急剧降低，未发现经典乳房炎相关属富集。AM组中，双歧杆菌属（Bifidobacterium）和盐单胞菌属（Halomonas）等共生或条件有益类群出现恢复性增加，而致病性优势属未回升至BM水平，表明从初乳到早期泌乳阶段微生物组经历渐进式生态重构和部分稳定。
**3.2. 泌乳阶段ARGs动态**
**3.2.1. ARGs组成与多样性**
ARGs类别和耐药机制在所有阶段均存在显著组间差异。抗生素类别水平上，多药、肽类、糖肽类和四环素类耐药是主要组分，但相对丰度随时间变化：BM中多药耐药基因占优势，ACM中肽类耐药基因相对富集，AM中糖肽类和四环素类耐药逐渐增加。耐药机制水平上，抗生素靶点改变和外排介导耐药在所有阶段占主导。Alpha多样性显示ACM的ARG Shannon指数显著低于BM（p = 0.001066），AM部分恢复但未完全回复。Beta多样性分析表明ARG谱系呈现阶段依赖性分离（PERMANOVA R
2 = 0.138, p = 0.007）。
**3.2.2. 基于LEfSe的阶段特异性ARG富集**
LEfSe分析显示，BM组富集最广泛的耐药类别，包括林可酰胺、膦酸、elfamycin、夫西地酸、氟喹诺酮、氨基糖苷和多药耐药。ACM组唯一富集吡嗪类耐药基因（可能反映携带广谱应激响应元件的类群存活）。AM组富集糖肽类和双环霉素类耐药基因，可能反映对渗透压和营养波动耐受的革兰氏阳性菌再扩展。
**3.3. 泌乳阶段MGEs动态**
**3.3.1. MGEs组成与多样性**
BM组中转座酶和整合酶相关元件丰度较高，表明干奶前水平基因转移潜力增强。ACM组重组酶相关MGEs富集，可能反映初乳免疫选择下的同源重组。AM组中转座子和整合酶相关元件富集，表明早期泌乳微生物重组过程中基因组选择性重构。MGE多样性模式与ARGs平行：ACM Shannon指数最低（p = 0.004876 vs BM），Beta多样性分析显示MGE谱系显著分离（PERMANOVA R
2 = 0.239, p = 0.001）。
**3.3.2. 基于LEfSe的阶段特异性MGEs富集**
BM组富集最多MGEs，包括IS200/IS605、DDE、IS607、ISAs1、IS5/IS1182、IS66、IS1595、ISCku10、高丰度转座酶及DNA重组/整合酶（RecG、位点特异性整合酶），表明干奶前水平基因转移潜力升高。ACM组唯一富集IS630和IS1，与基因组缩减和胁迫适应相关。AM组独特富集IS982、RmuC、IS1380、ISVsa10、IS91和位点特异性重组酶（XerC/XerD），可能反映微生物重组期间基因组可塑性增加。
**3.4. 泌乳阶段共现网络动态**
BM网络结构最复杂（平均度5.10，聚类系数0.75，模块化系数0.88），节点主要来自假单胞菌门（32.4%）、拟杆菌门（23.6%）和芽孢杆菌门（17.6%），ARGs和MGEs节点占比小（10.8%和2.4%），所有边为正相关。ACM网络收缩且连接性降低（平均度3.85，聚类系数0.95，模块化系数0.76），ARGs节点比例升至19.74%，MGEs升至2.63%，乳房炎相关属仍可检测但代表性降低。AM网络节点和边数最多（223节点，3906边），假单胞菌门占44.84%，ARGs节点降至14.8%，MGEs升至2.69%，双歧杆菌属和盐单胞菌属出现，正相关占99.82%。网络结构呈现“暴露–瓶颈–重组”轨迹，与微生物、ARGs和MGEs的阶段性变化一致。
**讨论与结论**
研究讨论指出，乳腺微生物组遵循清晰的“暴露–瓶颈–重组”轨迹，与先前纵向研究一致。干奶前开放的乳头管促进环境微生物涌入，初乳施加强烈的免疫过滤，早期泌乳支持微生物逐步重组。ARGs和MGEs的平行重构表明功能基因池响应与微生物演替相同的选择压力。初乳阶段抗性组和移动组多样性的收缩及早期泌乳的部分恢复，与宿主相关微生物组在抗微生物或免疫胁迫下的模式相似。早期泌乳中IS91和XerC/XerD等元件的富集提示基因组可塑性和潜在水平基因转移增加。网络分析进一步支持阶段性重塑：BM网络最大最复杂，ARG/MGE节点贡献最小；ACM网络收缩且ARG/MGE中心性增加；AM网络再扩大且ARG/MGE参与度下降但仍高于BM。研究结论部分翻译如下：这些发现共同表明，乳腺微生物组在干奶期和早期泌乳期经历协同的生态和功能重塑。微生物组成、ARGs和MGEs均表现出同步的瓶颈和选择性恢复，反映了宿主生理和抗生素干预驱动的深度重塑。尽管早期泌乳群落表现出部分生态稳定性和多种乳房炎相关类群丰度的降低，但与干奶前微生物组的持续差异以及ARG和MGE谱系的不完全恢复表明，干奶后乳腺生态系统的重组可能具有持久的功能后果。AM群落显示组成趋于稳定，但ARGs和MGEs相对于BM的水平升高表明功能恢复滞后于组成恢复，这对乳房健康和抗菌药物管理具有潜在影响。这些变化最终对乳房健康是有益、适应还是潜在有害，需要整合宿主炎症反应和乳房炎结局的进一步纵向研究。
