在群體遺傳學中,有效種群大小是指一個種群中有能力繁殖下一代的個體的數量。美國遺傳學家休厄爾·賴特在兩篇標誌性的論文(Wright 1931, 1938)中引入了有效群體大小(effective population size,又做有效種群大小)這一概念。[1][2]有效群體大小是很多群體遺傳學模型中的基本參量,通常小於絕對的種群大小(N)。參見小種群大小。
有效群體大小可能有兩種定義方式,即方差有效群體大小(variance effective size)和近交有效群體大小(inbreeding effective size)。兩種定義由F-統計衍生而來並緊密相連。
另一種定義,近交有效群體大小通過近交系數相鄰兩代之間的變化來計算,此時Ne定義爲和理想群體的具有相同近交變化的種群大小。此處文字參照Kempthorne (1957)。
對於理想群體,近交系數依以下遞歸方程式計算
。
利用回歸交配指數(1 − F)代替近交指數得到近似遞歸方程式
。
每一世代的差值是
。
通過解
![{\displaystyle {\frac {P_{t+1}}{P_{t}}}=1-{\frac {1}{2N_{e}^{(F)}}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/1be7ff707d1b1aaefd08111db87bb71f6d912aa9)
可得出近交有效群體大小爲
![{\displaystyle N_{e}^{(F)}={\frac {1}{2\left(1-{\frac {P_{t+1}}{P_{t}}}\right)}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/9339434ded19901cb2be29454e8b223630184c2b)
,儘管研究者很少直接使用這個方程式。
種群大小隨時間變化。假設有t個互相不重疊的世代,那麼有效種群大小由種群的調和平均數得出:
![{\displaystyle {1 \over N_{e}}={1 \over t}\sum _{i=1}^{t}{1 \over N_{i}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/ae7418227e9991f5e2ea2b6136ff07580d3de423)
例如,若6個世代(t = 6)的種群大小N = 10, 100, 50, 80, 20, 500,則有效種群大小爲:
|
|
|
|
|
|
|
|
注意調和平均數的值小於算術平均數的大小。後者在本例中爲126.7。
尤其需要注意的是種群瓶頸效應,即有效群體大小受最低種群數量的嚴重影響。
另外一種計算家族有效種群大小Ne的方法爲Ne = (4N)/(Vk + 2) 其中Vk是種群大小的方差。家族大小變化增大是不利的,因爲這將導致家族的近親交配。
對於雌雄異體的種群,即沒有自體受精的情況,則
![{\displaystyle N_{e}=N+{\begin{matrix}{\frac {1}{2}}\end{matrix}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/c88c9e80151bdf0429da70a038284cdd24f345d3)
或者更普遍地,
![{\displaystyle N_{e}=N+{\begin{matrix}{\frac {D}{2}}\end{matrix}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/65dcc068a71b12eeaef94cac3c9c74a22dddadfe)
其中D表示雌雄異體特性,可取0(非雌雄異體)或者1(雌雄異體)。
當N很大時,Ne近似等於N,因此這通常不重要,可以忽略,即:
![{\displaystyle N_{e}=N+{\begin{matrix}{\frac {1}{2}}\approx N\end{matrix}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/1d456f8336fa20763a81a9ca762d5c95d9cd3184)
如果種群大小保持恆定,每個個體必須平均向下一代貢獻兩個配子。一個理想種群假設這遵循泊松分佈(Poisson distribution),此時每個個體貢獻配子數k的方差也等於貢獻配子的平均數,即2:
![{\displaystyle \operatorname {var} (k)={\bar {k}}=2.}](//wikimedia.org/api/rest_v1/media/math/render/svg/ccaac8057a2d95281158d83028ee2fa829f3aa7e)
然而,在自然種群中,方差通常大於此:
![{\displaystyle \operatorname {var} (k)>2.}](//wikimedia.org/api/rest_v1/media/math/render/svg/1778e61cbd67e5fc06f15292fbdc69aa8fac9c44)
則有效種群數量爲
![{\displaystyle N_{e}^{(v)}={4N-2D \over 2+\operatorname {var} (k)}}](//wikimedia.org/api/rest_v1/media/math/render/svg/c8a0b2ea1d327e91bc456f7f2d7e3193b7b94802)
注意如果k的方差小於2,則Ne大於N。生育能力造成的遺傳偏差通常造成Ne值下降。
當生物的生存時間跨越多於一個生育季節,有效種群大小需要考慮進物種的生命表(life table)。
類似地,對於具有離散年齡結構的二倍體種群的近交有效種群大小也可以計算。這首先由Johnson(1977)提出,但本文的符號更接近Emigh和Pollak(1979)。
假定生命表的基本參量和上述單倍體情況相同,但區分雌雄性,如用N0ƒ 和N0m分別表示新生的雌性和雄性個體數量(注意小寫ƒ代表雌性而大寫F表示近交)。則近交有效種群大小爲
![{\displaystyle {\begin{aligned}{\frac {1}{N_{e}^{(F)}}}={\frac {1}{4T}}\left\{{\frac {1}{N_{0}^{f}}}+{\frac {1}{N_{0}^{m}}}+\sum _{i}\left(\ell _{i+1}^{f}\right)^{2}\left(v_{i+1}^{f}\right)^{2}\left({\frac {1}{\ell _{i+1}^{f}}}-{\frac {1}{\ell _{i}^{f}}}\right)\right.\,\,\,\,\,\,\,\,&\\\left.{}+\sum _{i}\left(\ell _{i+1}^{m}\right)^{2}\left(v_{i+1}^{m}\right)^{2}\left({\frac {1}{\ell _{i+1}^{m}}}-{\frac {1}{\ell _{i}^{m}}}\right)\right\}.&\end{aligned}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/6738d459d4b8b25e956f8618a277c672f1b751ed)
- Minimum viable population