Consider first a memoryless population model described by the usual branching process with a given mean reproduction matrix on a finite space of types. Motivated by the consequences of atavism in Evolutionary Biology, we are interested in a modification of the dynamics where individuals keep full memory of their forebears and procreation involves the reactivation of a gene picked at random on the ancestral lineage. By comparing the spectral radii of the two mean reproduction matrices (with and without memory), we observe that, on average, the model with memory always grows at least as fast as the model without memory. The proof relies on analyzing a biased Markov chain on the space of memories, and the existence of a unique ergodic law is demonstrated through asymptotic coupling.
論文ID : 2501.00306タイトル : On a population model with memory著者 : Jean Bertoin(チューリッヒ大学数学研究所、スイス)分類 : math.PR(確率論)発表日 : 2024年12月31日(arXiv プレプリント)論文リンク : https://arxiv.org/abs/2501.00306 本論文は、有限型空間上で与えられた平均繁殖行列を持つ標準分枝過程により記述される記憶のない種族モデルを考察する。進化生物学における隔世遺伝現象に着想を得て、著者は修正された動力学モデルを研究する。このモデルでは、個体が祖先の完全な記憶を保持し、繁殖過程は祖先系統からランダムに選択された遺伝子の再活性化を伴う。2つの平均繁殖行列(記憶ありと記憶なし)のスペクトル半径を比較することにより、平均的に記憶ありモデルは常に記憶なしモデルと同程度かそれ以上の速度で成長することが観察される。証明は記憶空間上の偏りを持つマルコフ連鎖の分析に依存し、漸近的結合により唯一遍歴性の存在が示される。
論文が研究する中心的な問題は、種族モデルに遺伝記憶機構を導入することが種族の成長に有利かどうかである。具体的には、個体が祖先の遺伝情報を「記憶」でき、繁殖時にこれらの情報をランダムに再活性化できる場合、種族の平均成長率はどのように変化するのか。
研究の動機は進化生物学における隔世遺伝 (atavism)現象に由来する:
遺伝記憶 :DNA内に過去の事象の記憶が保存され、表現型上消失した形質はDNA内に完全には消失していない遺伝子休眠 :ある遺伝子は多くの世代にわたって休眠状態を保つ可能性がある形質の再出現 :遠い祖先の形質が再び現れることがあり、通常は突然変異または遺伝子制御の欠陥により休眠遺伝子が活性化される著者は、自然選択と環境変化がない場合でも、祖先の遺伝物質をむやみに再活性化する可能性だけが、記憶なしモデルと比較して常に有益であることを指摘する。これは隔世遺伝の種族への利点が適応性の向上だけでなく、前世代の記憶を用いた繁殖から生じる本質的な利益にも由来することを示唆している。
理論的結果 :記憶を持つ分枝過程のスペクトル半径が常に記憶なしモデルのスペクトル半径以上であることを証明した(定理1.1:r ˉ ≥ r \bar{r} \geq r r ˉ ≥ r )数学的枠組み :記憶を持つ種族モデルを分析するための完全な数学的枠組みを確立し、作用素論とマルコフ連鎖法を含む技術的革新 :偏りを持つマルコフ連鎖と漸近的結合技術により唯一遍歴性を証明した生物学的洞察 :遺伝記憶機構の本質的な利点を数学的に証明し、進化生物学に理論的支援を提供した状態空間 :有限型集合 S S S 繁殖核 :π ( s , ⋅ ) \pi(s, \cdot) π ( s , ⋅ ) 、与えられた親の型における子孫の型の分布を指定平均繁殖行列 :m ∈ R + S × S m \in \mathbb{R}_+^{S \times S} m ∈ R + S × S 、ここで
m ( s , t ) = ∑ n ∈ N S n t π ( s , n ) m(s,t) = \sum_{n \in \mathbb{N}^S} n_t \pi(s,n) m ( s , t ) = ∑ n ∈ N S n t π ( s , n ) スペクトル半径 :r r r (Perron-Frobenius固有値)記憶空間 :S N = { s = ( s 0 , s 1 , s 2 , … ) } S^{\mathbb{N}} = \{s = (s_0, s_1, s_2, \ldots)\} S N = { s = ( s 0 , s 1 , s 2 , … )} 、ここでs j s_j s j はj j j 世代前の祖先の型を表す記憶活性化 :確率測度τ = ( τ ( j ) : j ≥ 0 ) \tau = (\tau(j) : j \geq 0) τ = ( τ ( j ) : j ≥ 0 ) は第j j j 世代の祖先記憶を活性化する確率を記述繁殖機構 :個体はまずτ \tau τ に従って祖先の世代T T T を選択し、その後π ( s T , ⋅ ) \pi(s_T, \cdot) π ( s T , ⋅ ) に従って子孫を産生作用素定義 :
m ˉ f ( s ) = ∑ j = 0 ∞ τ ( j ) ∑ t ∈ S m ( s j , t ) f ( t s ) \bar{m}f(s) = \sum_{j=0}^{\infty} \tau(j) \sum_{t \in S} m(s_j, t)f(ts) m ˉ f ( s ) = ∑ j = 0 ∞ τ ( j ) ∑ t ∈ S m ( s j , t ) f ( t s ) 記号空間S ˉ = S × S N \bar{S} = S \times S^{\mathbb{N}} S ˉ = S × S N を導入し、ここで( s , s ) (s, \mathbf{s}) ( s , s ) は個体の活性化された型と完全な記憶を記録する。転移核を構成する:
Q f ( s , s ) = ∑ j = 0 ∞ τ ( j ) ∑ t ∈ S m ˉ ( s j , t ) f ( s j , t s ) Qf(s, \mathbf{s}) = \sum_{j=0}^{\infty} \tau(j) \sum_{t \in S} \bar{m}(s_j, t)f(s_j, t\mathbf{s}) Q f ( s , s ) = ∑ j = 0 ∞ τ ( j ) ∑ t ∈ S m ˉ ( s j , t ) f ( s j , t s )
ここでm ˉ ( s , t ) = m ( s , t ) h ( t ) r h ( s ) \bar{m}(s,t) = \frac{m(s,t)h(t)}{rh(s)} m ˉ ( s , t ) = r h ( s ) m ( s , t ) h ( t ) は正規化された転移行列である。
作用素の冪とマルコフ連鎖の関係を確立する:
m ˉ k f ( s ) = r k E s [ f ( X k ) × ∏ j = 1 k h ( Y j ) h ( X j ) ] \bar{m}^k f(s) = r^k \mathbb{E}_s\left[f(X_k) \times \prod_{j=1}^k \frac{h(Y_j)}{h(X_j)}\right] m ˉ k f ( s ) = r k E s [ f ( X k ) × ∏ j = 1 k h ( X j ) h ( Y j ) ]
状態空間が無限であり連鎖が完全な履歴記憶を保つため、従来の強Feller性質は失効する。漸近的結合を採用する:
初期化段階 :記憶が共通の接頭辞を持たない場合、2つの連鎖は独立に進化統合段階 :記憶が共通の接頭辞を持つ場合、共通部分を使用して同じ型を活性化陳述 :r ˉ ≥ r \bar{r} \geq r r ˉ ≥ r 、すなわち記憶ありモデルのスペクトル半径は記憶なしモデルのスペクトル半径以上である。
証明の概要 :
many-to-one公式を通じて問題を積関数∏ j = 1 k h ( Y j ) h ( X j ) \prod_{j=1}^k \frac{h(Y_j)}{h(X_j)} ∏ j = 1 k h ( X j ) h ( Y j ) の分析に変換 Birkhoff遍歴定理を利用し、以下を証明する必要がある:
lim inf k → ∞ 1 k ∑ j = 1 k ( log h ( Y j ) − log h ( X j ) ) ≥ 0 \liminf_{k \to \infty} \frac{1}{k} \sum_{j=1}^k (\log h(Y_j) - \log h(X_j)) \geq 0 lim inf k → ∞ k 1 ∑ j = 1 k ( log h ( Y j ) − log h ( X j )) ≥ 0 重要な点は偏りを持つ連鎖の唯一遍歴性と不変測度の性質を証明することである 偏りを持つマルコフ連鎖Q Q Q は唯一遍歴的である、すなわち唯一の不変確率測度σ \sigma σ が存在する。
不変測度σ \sigma σ の周辺分布は以下を満たす:
σ ( { ( s , s ) : s = t } ) = σ ( { ( s , s ) : s j = t } ) = ϱ ( t ) h ( t ) \sigma(\{(s, \mathbf{s}) : s = t\}) = \sigma(\{(s, \mathbf{s}) : s_j = t\}) = \varrho(t)h(t) σ ({( s , s ) : s = t }) = σ ({( s , s ) : s j = t }) = ϱ ( t ) h ( t )
無限状態空間 :記憶空間S N S^{\mathbb{N}} S N は無限である強Feller性質の失効 :完全な記憶保持のため、従来の分析ツールは適用不可既約性の欠如 :古典的マルコフ連鎖理論の標準的仮定が満たされない漸近的結合 :2つの軌道が最終的に無限に接近するが完全には一致しない結合を構成コンパクト化技術 :疑似距離d ( s , s ′ ) = a ∣ s ∧ s ′ ∣ d(s,s') = a_{|s \wedge s'|} d ( s , s ′ ) = a ∣ s ∧ s ′ ∣ を通じて状態空間をコンパクト化多層的分析 :作用素論、遍歴理論、結合技術を組み合わせS = { a , b } S = \{a,b\} S = { a , b } 、m = ( 1 1 1 2 ) m = \begin{pmatrix} 1 & 1 \\ 1 & 2 \end{pmatrix} m = ( 1 1 1 2 ) 、τ ( 0 ) = u \tau(0) = u τ ( 0 ) = u 、τ ( 1 ) = 1 − u \tau(1) = 1-u τ ( 1 ) = 1 − u を考える。
記憶なしモデルのスペクトル半径:r = 3 + 5 2 ≈ 2.618 r = \frac{3+\sqrt{5}}{2} \approx 2.618 r = 2 3 + 5 ≈ 2.618 記憶ありモデルのスペクトル半径:r ˉ ( u ) \bar{r}(u) r ˉ ( u ) はu u u の関数であり、すべてのu ∈ ( 0 , 1 ) u \in (0,1) u ∈ ( 0 , 1 ) に対してr ˉ ( u ) ≥ r \bar{r}(u) \geq r r ˉ ( u ) ≥ r が成立 論文はr ˉ ( u ) \bar{r}(u) r ˉ ( u ) がu u u とともに変化する図を提供し、不等式が通常は厳密であることを示す。
平均繁殖行列m m m が平衡的である場合(すなわち∑ t ∈ S m ( s , t ) \sum_{t \in S} m(s,t) ∑ t ∈ S m ( s , t ) がs s s に依存しない)、r ˉ = r \bar{r} = r r ˉ = r であり、この場合記憶は追加的な利点を提供しない。
論文は本研究を以下の分野と関連付ける:
強化Galton-Watson過程 :Bertoinとmaleinの最近の研究3,4 との類似性ペナルティ付きマルコフ過程 :Champagnatら7 の研究と関連多型分枝過程 :古典理論1,12 の拡張に基づく遺伝記憶は常に有益である :選択圧がない場合でも、記憶機構は種族の成長率を向上させることができる数学的機構 :利点は多産的祖先の遺伝物質が現在の種族内に大きな部分で存在し、したがってより再活性化される可能性が高いことから生じる理論的完全性 :記憶を持つ分枝過程を分析するための完全な数学的枠組みを提供したモデルの単純化 :環境変化と自然選択を考慮していない記憶活性化の公平性 :活性化確率が記憶内容に依存しないと仮定技術的限界 :古典的結果(完全なPerron-Frobenius定理など)は無界支持の場合に確立が困難記憶内容に依存する活性化機構を考慮 環境変化と選択圧を導入 有限記憶長の最適化問題を研究 他の生物学的状況への応用を探索 理論的厳密性 :証明技術は精巧であり、特に漸近的結合の応用が優れている生物学的関連性 :隔世遺伝現象に数学理論的支援を提供技術的革新 :無限記憶の場合における従来のマルコフ連鎖理論の困難を克服結果の明確性 :主要定理は簡潔明瞭で、直感的な生物学的解釈を持つ実用的応用の限界 :モデルは過度に理想化されており、実際の生物系との乖離が大きい計算の複雑性 :実際の計算では、無限記憶がもたらす複雑性は禁止的である可能性があるパラメータ感度 :記憶活性化分布τ \tau τ の選択に関する指導が不足している理論的貢献 :確率論における分枝過程理論に新しい次元を追加学際的価値 :数学と進化生物学を結合方法論的意義 :漸近的結合技術は無限記憶を持つ他の確率過程に適用可能進化生物学における理論的モデリング 履歴依存性を持つ確率系の分析 遺伝的アルゴリズムと進化計算の理論的基礎研究 論文は分枝過程理論、マルコフ連鎖理論、結合技術、進化生物学など複数の分野にわたる14篇の関連文献を引用し、研究に堅実な理論的基礎を提供している。