[論文レビュー] A General Formula for the Generation Time

François Bienvenu, Lloyd Demetrius+1 more|arXiv (Cornell University)|Jul 25, 2013

Gene Regulatory Network AnalysisBiochemistry, Genetics and Molecular Biology参考文献 20被引用数 20

ひとこと要約

本稿では、有限マルコフ連鎖フレームワークを用いて、構造的集団における生成期間の一般式を提案し、重み付き有向グラフにおける生殖遷移への最初の帰還時間の平均として定義する。主な結果は、任意の原始的行列が表すライフコースに適用可能な簡単な閉形式の式である：$ T = 1 / ackslashsum_{[j o i]ackslashin R} e_{ij} $、ここで $ e_{ij} $ は生殖弧の弾力性を表す。

ABSTRACT

We show that the generation time -- a notion usually described in a biological context -- can be defined in a general way as a return time in a conveniently constructed finite Markov chain. The simple formula we obtain agrees with previous results derived for structured populations projected in discrete time, and allows to define the generation time of any process described by a weighted directed graph whose matrix is primitive.

研究の動機と目的

• 年齢区分モデルにとどまらず、任意の構造的集団に一般化された生成期間の概念を定義すること。
• 生成期間を確率過程における生殖遷移への最初の帰還時間の平均として定義すること。
• 任意の原始的行列モデルに適用可能な、簡単で普遍的な生成期間の公式を導出すること。
• 従来の人口動態モデルにおける生成期間の複雑な公式を統一・簡略化すること。
• 生成期間と人口成長率に対する生殖パラメータへの弾力性との関係を確立すること。

提案手法

• ライフコースの重み付き有向グラフから有限マルコフ連鎖を構築し、ノードを段階、弧を遷移として定義する。
• 生殖弧 $ R $ を新たに個体を生成する遷移（例：有性生物における精巣融合）として特定する。
• 弧上のマルコフ連鎖の遷移行列 $ \tilde{\mathbf{p}} $ を、生殖（$ R $）と非生殖（$ S $）の集合に分割する。
• 弧レベルのマルコフ連鎖の定常分布 $ \tilde{\boldsymbol{\pi}} $ を計算し、それをスケーリングして弧上の正規化分布 $ \boldsymbol{\varpi} $ を得る。
• 再生に至る最初の通過時間分布を導出し、生成期間を帰還時間の期待値として定義する。
• 行列の弾力性理論を用いて、主要固有値 $ \lambda $ が生殖および非生殖遷移の変化に対してどのように感度を示すかを表現する。

実験結果

リサーチクエスチョン

• RQ1 年齢区分モデルに限らず、任意の構造的集団モデルに対して生成期間を一般化して定義可能か？
• RQ2 複雑なライフコースにおける世代間の平均期間の背後にある数学的構造は何か？
• RQ3 人口成長率 $ \lambda $ に対する生殖遷移への弾力性と生成期間の関係は何か？
• RQ4 集団モデルのグラフ構造に基づいて、普遍的かつ単純な生成期間の公式を導出可能か？
• RQ5 弧上の定常分布が生殖イベントへの平均帰還時間の決定に果たす役割は何か？

主な発見

• 生成期間は簡単な式 $ T = 1 / \backslashsum_{[j\to i]\backslashin R} e_{ij} $ で与えられ、ここで $ e_{ij} $ は生殖遷移の弾力性を表す。
• この式は、年齢区分モデル（例：レスリー行列）に対する従来の結果を一般化し、識別可能な生殖弧を持つ任意の原始的行列モデルに適用可能である。
• 導出過程から、生成期間がライフコースグラフの弧上のマルコフ連鎖における生殖遷移への最初の帰還時間の平均に一致することが示された。
• すべての生殖弧に共通乗数 $ c $ をかけるときの人口成長率 $ \lambda $ に対する弾力性は $ e_\lambda(c) = 1/T $ であり、生成期間が直接的に人口動態感受性に関連している。
• 非生殖遷移については、弾力性の和が $ e_\lambda(d) = 1 - 1/T $ となる。これは、長寿命種では若年期生存率の感度が低く、成人期生存率の感度が高いことを示している。
• 本手法は、サイズ構造的モデルやメタポピュレーションモデルを含む多様な生物学的システムにおいて、生成期間を統一的に計算するフレームワークを提供する。