No persistent circadian oscillator at genome resolution: pseudo-coherence in gut microbiome dynamics

技術要約：ゲノム解像度における持続的な概日振動子の不在

問題提起
腸内細菌叢の日周リズムに関する支配的な解釈は、結合位相振動子という概念的枠組みに依存している。この見解は、微生物コミュニティが共通の時計に位相ロックされた固有または宿主同調型の振動子を含んでおり、それが観察されるリズム性、反同期クラスター、およびスペクトルピークを生成していると仮定している。しかし、この解釈は、システムがホップ分岐（Hopf bifurcation）の近傍で動作しているか、あるいは固有の振動モードを有しているという強い物理的仮定に基づいている。

本論文は、ゲノム解像度（メタゲノム構築ゲノム、すなわちMAGレベル）の腸内細菌叢ダイナミクスの時間-周波数構造が、振動子の存在しない安定した確率論的レジームに由来する可能性があるのではないか、という問いを立てることで、この仮定に異議を唱える。具体的には、「非正規過渡増幅（non-normal transient amplification）」——すなわち、相互作用行列の幾何学的特性であり、非直交な固有ベクトルによって摂動が減衰する前に一時的に増幅される現象——が、基礎となる振動子なしに、集団的な振る舞いを擬似的なコヒーレント・パターンへと組織化できるのかを調査する。

方法論
著者らは、2匹のマウス（マウスAおよびマウスB）から得られた、1時間間隔の2週間の時系列データ[45]を再解析した。解析パイプラインは、以下の4つの「擬似コヒーレンス（pseudo-coherence）」のシグネチャーをテストするように設計された。ここで、擬似コヒーレンスとは、幾何学的増幅が確率的なゆらぎを低次元の反応部分空間へと再形成するレジームを指す。

1. 局所ヤコビ行列の推定: 4つの連続する1時間ごとの観測値を用い、最小二乗更新法によって局所ヤコビ行列（$\hat{A}_k$）を推定した。高次元行列全体を分析するのではなく、交換子 $B_k = [\hat{A}_k, \hat{A}_k^\top]$ を対角化することで、ランク2の部分空間を抽出した。この部分空間は2つのモードを定義する：すなわち、確率的な強制力を吸収する「非正規モード（non-normal mode）」と、ゆらぎが一時的に増幅される「反応モード（reaction mode）」である。
2. 非正規性の診断: システムが反応的レジームにあるか（$K/K_c > 1$）を判断するために、幾何学的閾値（$K_c$）に対する非正規性指数（$K$）を算出した。また、これらのモードの「サポート（support）」（$s_r, s_n$）を追跡し、増幅がどの程度広範にMAG（Metagenome-Assembled Genome）間に分布しているかを測定した。
3. 位相およびクラスター分析: 非パラメトリック位相決定（NPPD）法を用いて各MAGに位相を割り当て、クラスター化した。極めて重要な点として、位相情報に依存しない、交換子から導出された反応モード成分（$r_i$）の符号のみに基づく「位相非依存型」のクラスター回復手法を開発した。
4. スペクトルおよび共分散テスト:
   • 時間-周波数解析: モルレット・ウェーブレット・スケーログラムを生成し、持続的なリッジ（固定された時間-周波数構造）の有無を探索した。
   • サロゲート・テスト: 限界スペクトルのみから期待される値よりも観測された低周波成分のパワーが大きいかどうかをテストするため、振幅調整フーリエ変換（AAFT）サロゲートを生成し、偽発見を制御した。
   • 時間反転対称性: ラグ付き共分散（$C(\tau)$）を分析し、非対称性（$I(\tau)$）を調べた。これは不可逆性と循環的な確率流のシグネチャーである。
   • エントロピー生成: $\tau \to 0$ におけるラグ付き不均衡の傾きから、局所エントロピー生成のプロキシ（$\Sigma_{local}$）を導出した。

主な結果
解析の結果、両個体において一貫した4つの知見が得られ、これらは擬似コヒーレンス仮説を強力に支持し、持続的振動子仮説を否定するものであった。

1. サポートに駆動される断続的な位相整列: マクロな位相コヒーレンス（クラスター秩序パラメータ）は、断続的であることが判明した。決定的なことに、これらのコヒーレンスのピークのタイミングは、スペクトルの不安定性への近接性ではなく、反応モードの空間的サポート（$s_r$）と一致していた。サポートと秩序パラメータ間のウェーブレットコヒーレンスは、すべての周波数帯域で高かった（$0.55–0.70$）であり、Granger因果関係はサポートがコヒーレンスを駆動していることを示唆した。
2. 持続的なリッジの不在: 時間-周波数表現（スケーログラム）には、サーカディアン帯域を含むどの周波数においても、持続的なリッジは見られなかった。時間平均スペクトルでは低周波領域（1/24から1/16 $h^{-1}$の間でピーク）に弱い増強が見られたものの、これは固定された特徴ではなく、ドリフトするパワーの集中であった。サロゲート・テストにより、この低周波の過剰は統計的に有意ではあるものの、振動子のような定常的な構造を欠いていることが確認された。
3. 時間反転対称性の破れ: ラグ付き共分散は著しく非対称であり、グローバルな不均衡ピーク（$I(\tau)$）は中間ラグ（数十時間）で発生した。これは、循環的な確率流を伴う非平衡定常状態を示している。原点におけるこの不均衡の傾きは、両個体において一貫して正の準定常的なエントロピー生成の推定値を与えた。
4. 栄養学的ギルドの生物学的回復: 反応モードの符号に基づく位相非依存型クラスターリングは、元の位相ベースのクラスターと高い一致度（>93%）で、2つの明確な機能的ギルドを正常に回収した。
   • クラスター1: Bacteroidota（例：Muribaculaceae, Bacteroidaceae）が支配的であり、主要な多糖類分解者として特定された。
   • クラスター2: Bacillota A（例：Lachnospiraceae, Butyricicoccaceae）が支配的であり、二次的な短鎖脂肪酸発酵者として特定された。
     この回復は、位相情報を一切用いずに行われたため、従来観察されてきた「夜／昼」の分割は、実際にはクロスフィーディング（共生的な栄養摂取）の上流／下流の栄養学的アーキテクチャを反映しているものであることが示唆された。

意義および主張
本論文は、ゲノム解像度の腸内細菌叢ダイナミクスは、結合振動子系ではなく、安定した、強く非正規な確率論的レジームとして記述されるのが最適であると主張している。

• リズム性の再解釈: 観察される「リズム」は、固有または宿主同調型の振動子の証拠ではない。むしろ、それらは幾何学的増幅が確率的なノイズを低次元で断続的な同期様パターンへと再形成する「擬似コヒーレント」なエピソードである。
• メカニズム: システムは、直接的なクロスフィーディングのカスケード（一次分解者 $\to$ 二次発酵者）を通じて動作している。この相互作用ネットワークの非正規な構造がゆらぎを増幅し、ホップ分岐を必要とすることなく、過渡的な集団的逸脱と時間反転対称性の破れを生み出している。
• 帰無仮説: 著者らは、ゲノム解像度の微生物叢ダイナミクスに対する適切な帰無仮説は、位相ロックされた振動子ではなく、安定した非正規な確率論的システムであると提案している。
• 検証可能な予測: 本論文は、時計遺伝子のノックアウト個体（例：Bmal1, Per ノックアウト）を用いた具体的なテストを提案している。非正規モデルによれば、宿主の時計にロックされた24時間成分は弱まる可能性があるものの、過渡的増幅構造、リード・ラグの非対称性、および2つのクラスターによるギルド構成は維持されると予測される。逆に、振動子モデルであれば、リズム構造そのものが完全に崩壊すると予測される。

結論として、著者らは、この文脈における「生命」を、非正規増幅によるフラックス（流束）のレンズを通して捉えることができると述べている。そこでは、腸内細菌叢は、振動子による位相ロックではなく、方向性を持ったクロスフィーディングによって維持される非平衡定常状態の経験的な実現例となっている。