由玉巖周軍英杜洋李瑩曹長雷

(1.圈養(yǎng)野生動物技術(shù)北京市重點實驗室，北京動物園，北京，100044；2.中國動物園協(xié)會，北京，100037；3.長江師范學(xué)院生命科學(xué)與技術(shù)學(xué)院，重慶，408100)

北京動物園圈養(yǎng)平原斑馬種群統(tǒng)計學(xué)與遺傳學(xué)分析

由玉巖1*周軍英杜 洋李 瑩曹長雷3

(1.圈養(yǎng)野生動物技術(shù)北京市重點實驗室，北京動物園，北京，100044；2.中國動物園協(xié)會，北京，100037；3.長江師范學(xué)院生命科學(xué)與技術(shù)學(xué)院，重慶，408100)

稿件運行過程

斑馬； 動物園； 圈養(yǎng)種群； 繁育； 遺傳管理

圈養(yǎng)野生動物的目的在于使物種能夠長期生存，期望能夠在有限的空間和動物種群數(shù)量條件下，實現(xiàn)動物種群良性的可持續(xù)發(fā)展[1]，我們有必要把物種作為動態(tài)的進化實體進行保育并維持物種的遺傳健康[2]。早在20世紀70年代末，美國動物園最先發(fā)現(xiàn)，圈養(yǎng)野生動物的近親繁殖會影響動物種群的可持續(xù)發(fā)展[3]，因此，圈養(yǎng)野生動物種群的遺傳學(xué)管理便應(yīng)運而生。圈養(yǎng)種群的遺傳學(xué)管理是避免近親繁殖、保持種群較高基因多樣性水平非常必要的一種管理手段[3]，是動物園保護物種健康發(fā)展的需要[4]。通過對圈養(yǎng)野生動物種群的遺傳學(xué)管理，可以為不同物種建立真實的譜系，以確定特定個體對當(dāng)前和未來種群的遺傳重要性，并監(jiān)測遺傳多樣性的時空變化情況[5]。我國近些年有關(guān)圈養(yǎng)野生動物種群的遺傳學(xué)管理的報道多集中于珍稀瀕危鳥類[6-8]和靈長類[9-10]，而針對動物園種群較大，且繁育良好的食草動物往往予以忽視。然而，經(jīng)過近幾十年的繁衍，部分物種在繁育過程中逐步出現(xiàn)了種群衰退現(xiàn)象，如死胎、流產(chǎn)、幼崽生存率低等問題凸顯，同時，部分種群的管理也出現(xiàn)了問題，其中，平原斑馬(Equusburchelli)就是一個最為顯著的例子。

平原斑馬亦稱草原斑馬或普通斑馬。被列入《世界自然保護聯(lián)盟》(IUCN)2015年瀕危物種紅色名錄ver 3.1，野外種群現(xiàn)狀穩(wěn)定，該物種的生存尚未受到威脅。根據(jù)2002年IUCN的記錄，野外平原斑馬數(shù)量是66萬只，但目前數(shù)量仍在緩慢下降，原因主要來自人類活動的影響，如野生動物棲息地的占用、畜牧業(yè)的發(fā)展以及被家畜占用的水源地等原因。平原斑馬廣泛分布于非洲東部、中部和南部，棲息于水草豐盛的熱帶草原。斑馬是高度社群性的動物，其家庭結(jié)構(gòu)多為每個家庭有1只雄性斑馬，最多6只雌性斑馬及它們的子女組成小群生活[11]。而一些尚未結(jié)婚、交配的雄性斑馬則會自己單獨生活，或是跟其他雄性一起生活，直至它們有能力去挑戰(zhàn)有家室的雄性斑馬[11]。在繁殖季節(jié)，斑馬表現(xiàn)得十分緊張而活躍，雄獸之間互相會毫不客氣地進行激烈的爭斗，打斗的方式為互相碰擊頸部、用嘴咬、用前蹄踢等[12]。野外生境中，每只獲勝的雄獸每年要交配數(shù)只雌獸。雌獸的妊娠期約為11～13個月，每隔3 a生產(chǎn)1次[12]，每胎產(chǎn)1仔，幼仔出生后不久，即可站立和走路，哺乳期約為6個月，3.5～4歲時性成熟。壽命約為20～30 a[12]。

北京動物園飼養(yǎng)斑馬的歷史可以追溯到1971年，1971～2013年，北京動物園圈養(yǎng)及繁育斑馬數(shù)量累計達到132只，在過去42 a中，北京動物園繁育的斑馬后代已經(jīng)遍布全國各大動物園。由于早期對斑馬的飼養(yǎng)和繁殖并沒有考慮到遺傳管理問題，致使早期斑馬的飼養(yǎng)繁育只注重數(shù)量增長，并未考慮種群的可持續(xù)發(fā)展及種群繁育力衰退問題。近些年斑馬的流產(chǎn)率非常高，本研究通過查詢北京動物園歷史檔案，建立斑馬譜系，運用譜系分析軟件進行種群遺傳學(xué)分析，并最終制定科學(xué)的種群管理策略，為斑馬的科學(xué)繁育提供依據(jù)，對圈養(yǎng)斑馬種群的健康發(fā)展具有重要的意義。

1 研究方法

1.1譜系調(diào)查

本研究對北京動物園自1971年以來飼養(yǎng)繁育的斑馬情況進行了調(diào)查，對歷史圈養(yǎng)斑馬的引進、輸出、出生、死亡、繁育后代等譜系信息進行了查詢與總結(jié)。

1.2數(shù)據(jù)分析

譜系數(shù)據(jù)錄入及分析使用ISIS(International Species Information System)開發(fā)的SPARKS v1.6軟件；種群統(tǒng)計學(xué)和遺傳學(xué)分析使用PMx v1.0軟件。通過PMx v1.0軟件分析并明確奠基者數(shù)量、潛在奠基者數(shù)量及奠基者編號，對奠基者后代數(shù)量進行統(tǒng)計。通過PMx v1.0軟件對種群數(shù)量變化情況及年齡結(jié)構(gòu)進行分析，對不同時期斑馬的死亡率進行統(tǒng)計，比較不同月份斑馬的出生及死亡情況，并計算種群整體基因多樣性保存比率、潛在的基因多樣性保存比率、奠基者基因等量值、平均親緣關(guān)系值、世代長度、有效種群數(shù)量，以及有效種群數(shù)量與種群數(shù)量的比值(Ne/N)等。

有效種群大小(Ne)是指一個種群中能將其基因連續(xù)傳遞到下一代的個體平均數(shù)，實際上相當(dāng)于理想狀態(tài)下(種群大小的穩(wěn)定，世代之間沒有重疊，所有個人以相同的概率給下一代貢獻配子)種群的大小。Ne=4NmNf/(Nm+Nf)，其中，Nm表示雄性個體數(shù)，Nf表示雌性個體數(shù)，Nm和Nf應(yīng)是指具有正常繁殖、生育能力和不致引起近緣變異和突變的數(shù)量。有效種群數(shù)量與種群數(shù)量的比值(Ne/N)，是判斷種群遺傳管理有效性的重要參數(shù)，Ne/N越接近1則表明有更多的個體適應(yīng)繁育后代，反之則有更少的個體適宜繁育后代。

通過使用PMx v1.0軟件計算MSI配對系數(shù)(MSI配對表是根據(jù)配對動物的平均親緣關(guān)系值、雄雌動物平均親緣關(guān)系值的差距、產(chǎn)生后代的近親繁殖系數(shù)、配對動物祖先未知因素的總量進行計算的一組數(shù)值，通過該表的繪制，可以明確斑馬的最佳配對，數(shù)值以1～7表示，1～4為適宜配對，>4為不適宜配對)。

2 結(jié)果

2.1圈養(yǎng)種群數(shù)量概況

1971～2013年，北京動物園斑馬圈養(yǎng)及繁育累計132只。剔除死胎、流產(chǎn)、死亡、來源和去向不明等個體，最終用于種群譜系分析的個體共計56只，其中雄性38只，雌性18只。

對1971～2013年北京動物園圈養(yǎng)斑馬種群繁育數(shù)量變化情況進行了統(tǒng)計，并將剔除前和剔除后的數(shù)據(jù)進行了比較(圖1)。剔除數(shù)據(jù)后，種群數(shù)量總體呈穩(wěn)定趨勢，平均增長率λ為0.98。

2.2種群年齡結(jié)構(gòu)與性比

目前，對斑馬現(xiàn)存活個體種群年齡結(jié)構(gòu)的分析顯示(圖2)，種群中育齡個體數(shù)量較多，0～3歲幼年個體數(shù)量不足，種群數(shù)量呈緩慢穩(wěn)定增長趨勢。歷史圈養(yǎng)斑馬種群總體性別比例為雄性：雌性=1.17∶1(性別不詳個體未列入該分析，下同)。現(xiàn)存活斑馬種群性別比例為雄性：雌性=2.11∶1，適宜繁殖年齡個體性別比例為1.76∶1，性別比例嚴重偏離斑馬一雄多雌的種群現(xiàn)狀。

圖1 1971～2013年北京動物園圈養(yǎng)斑馬種群數(shù)量變化情況Fig.1 Changes of captive population of zebra in Beijing Zoo from 1971 to 2013

圖2 北京動物園斑馬圈養(yǎng)種群年齡和性別結(jié)構(gòu)Fig.2 Age and sex structure of captive zebra population in Beijing Zoo

2.3圈養(yǎng)種群的繁殖、死亡情況

譜系數(shù)據(jù)統(tǒng)計結(jié)果顯示，斑馬的繁殖年齡通常多在5～21歲間，據(jù)記載斑馬最早繁殖年齡為2歲，最晚繁殖年齡為22歲。斑馬的生理壽命約為30歲，而北京動物園圈養(yǎng)斑馬最高存活年齡為32歲。雌性繁育周期平均1.11 a/胎。1歲前是幼齡斑馬的死亡高峰期。在對1971～2013年北京動物園圈養(yǎng)斑馬月出生、死亡數(shù)量情況的調(diào)查結(jié)果表明(圖3)，42 a間，斑馬在1月出生數(shù)量最多，累計達到35只個體；而6～7月死亡數(shù)量最多，累計達到15只，3月斑馬死亡數(shù)量較6～7月為少，累計死亡5只個體。

通過不同時期斑馬死亡率的統(tǒng)計結(jié)果(表1)顯示，幼年斑馬在30日齡死亡率較高，為11%，雌性、雄性30日齡死亡率相同，均為11%；剛出生斑馬的死亡率(含流產(chǎn))為18%，雄性死亡率為13%，雌性死亡率為24%；未成年個體的死亡率為4%，死亡個體全部為雄性；繁殖期和繁殖后期死亡個體數(shù)量明顯下降。

圖3 1971～2013年北京動物園圈養(yǎng)繁育斑馬月出生、死亡數(shù)量Fig.3 The born and death number of zebra from 1971 to 2013 in Beijing Zoo

表1不同時期斑馬死亡率統(tǒng)計

Tab.1 The zebra mortality rate in different period

由于2000年前，斑馬的死亡病例還不是十分健全，僅對2000年之后的當(dāng)年新生個體的死亡原因進行了統(tǒng)計(表2)。結(jié)果表明，新生個體死亡率為40%，且近一半的死亡原因為死胎和流產(chǎn)，以2004年死亡率最高為66.67%。

表2 北京動物園2000～2012年0～2歲斑馬死亡原因統(tǒng)計

Tab.2 The death reason statistics of 0 to 2 years old the zebra in Beijing Zoo

2.4 圈養(yǎng)種群遺傳學(xué)分析

北京動物園斑馬圈養(yǎng)種群主要來源于不同時期引進的7只奠基者(表3)，其中5只個體已經(jīng)死亡，現(xiàn)今種群中還有2只存活的奠基者已經(jīng)繁育了后代，存活個體處于育齡年齡的雄性個體23只，雌性13只，種群只保留了野生個體83%的基因多樣性。該種群的奠基者基因等量值是9，而奠基者的后代數(shù)明顯不均衡(圖4)，基因庫87%的基因僅來自于3只奠基者，其中T132號(♂，已故)個體的后代數(shù)最多，T62次之。目前，存活奠基者為T62(♂)，T73(輸出，♂)，T74(輸出，♂)。由于T62后代數(shù)過多，已不適宜再繼續(xù)產(chǎn)生后代，而T73和T74已經(jīng)輸出，不能滿足園區(qū)內(nèi)個體的繁育需求。其中，僅T132和T133為長期的配偶，其余奠基者均與非奠基者進行了配對。種群中，潛在奠基者的數(shù)量僅有2只，占整個種群數(shù)量的2.13%，如果能夠使?jié)撛诘旎叩幕蛉谌敕N群，則種群的基因多樣性水平可提升至94%。種群的平均親緣關(guān)系值(MK)是0.17，近似于半同胞間的近親繁殖。

表3 北京動物園圈養(yǎng)斑馬種群遺傳學(xué)參數(shù)

Tab.3 Genetic parameters of captive zebra population in Beijing Zoo

圖4 北京動物園歷史圈養(yǎng)斑馬種群奠基者后代數(shù)量Fig.4 Founder contribution and descendants of zebra of captive population in Beijing Zoo

2004年繁殖的數(shù)量最多為9只，2010年和2012年分別繁育后代7只，2001年和2005年繁育后代6只，其余年份繁育后代數(shù)在1～4只不等。后代年數(shù)量變化呈現(xiàn)增減式波動現(xiàn)象。根據(jù)配對適宜度指數(shù)(MSI)評估現(xiàn)存圈養(yǎng)斑馬種群情況，結(jié)果表明大部分都是不利于遺傳健康的配對繁育的后代。根據(jù)配對適宜度指數(shù)(MSI)表的分析結(jié)果(表4)，適合配對的遺傳優(yōu)勢個體已經(jīng)很少。

表4 圈養(yǎng)斑馬MSI配對表

Tab.4 Breeding pairs with various MSI scores in captive zebra

注：本MSI配對表僅列出了適宜配對的個體；“-”表示數(shù)值已經(jīng)>7，不適宜再配對

Note：The MSI pairing table only suitable pairing of individuals；“-”indicates values>7 which unsuitably to pair

3 討論

3.1 斑馬種群存在問題

1971～2013年，北京動物園斑馬圈養(yǎng)及繁育累計132只，圈養(yǎng)斑馬種群總數(shù)量呈持續(xù)穩(wěn)定趨勢，但如果不做有效的種群管理，種群數(shù)量將面臨下降。至今，經(jīng)由北京動物園繁育并輸出的后代已經(jīng)遍布全國，隨著斑馬種群的發(fā)展，問題也逐步凸顯。

鑒于北京動物園現(xiàn)存斑馬和歷史斑馬性別比來看，過多雄性后代的出現(xiàn)，加劇了斑馬種群中雄性個體對等級地位的競爭，是導(dǎo)致接近繁殖年齡雄性死亡率上升的直接原因。在動物園有限的空間內(nèi)，過多的雄性競爭將會導(dǎo)致僅有少數(shù)雄性個體能夠繁育后代，而那些潛在遺傳優(yōu)勢個體無法將其優(yōu)勢基因輸入到該種群中，出現(xiàn)了單一雄性后代數(shù)量過多，形成了單一性別個體的不足或過多而導(dǎo)致后代遺傳多樣性比例的不均衡發(fā)展的現(xiàn)象出現(xiàn)[13]。此外，數(shù)代的近親繁殖更加劇了遺傳多樣性的降低和稀有等位基因的丟失。因此，整個斑馬種群在更多雄性降生和近親繁殖的情況下逐步的衰退。

遺傳管理有效性指數(shù)是有效數(shù)量與種群數(shù)量的比率(Ne/N)，該指數(shù)是對未來種群遺傳狀況預(yù)測的重要指標(biāo)[13]。目前，斑馬種群中Ne/N的比率接近0.06，這一結(jié)果表明現(xiàn)存斑馬的有效種群已經(jīng)很少，再加之大量早期奠基者的死亡和潛在奠基者未能有效參與到種群的繁育之中，給現(xiàn)有的斑馬種群發(fā)展帶來了更大的問題。影響Ne/N最主要的因素是種群大小的波動，其次為家系大小的變化和性比失衡[13]。鑒于斑馬屬一雄多雌的家族體系，后代中更多雄性的出現(xiàn)會在一定程度上影響了有效種群大小[14]。目前，北京動物園圈養(yǎng)斑馬種群Ne/N的值很低，主要有兩方面的原因。首先，種群的數(shù)量始終維持穩(wěn)定；其次，種群中存在某些個體后代過多的情況。對北京動物園斑馬種群而言，過少的有效種群以及單一建群者后代過多，直接導(dǎo)致了斑馬種群遺傳多樣性的降低；更多近交繁育的出現(xiàn)，更加劇了斑馬種群的近交衰退現(xiàn)象，使得斑馬后代因死胎和流產(chǎn)的死亡的個體占了死亡總數(shù)的近一半比例。

北京動物園圈養(yǎng)斑馬種群在衰退的過程中，新生及幼年個體的死亡率也隨之升高。北京動物園斑馬死亡的主要原因是死胎、流產(chǎn)、外傷及疾病等原因。斑馬死亡多集中在早期，表現(xiàn)為出生時及30日齡前死亡率較高，早期流產(chǎn)或死胎是主要的死亡原因；30日齡前死亡個體均表現(xiàn)為體弱，抵抗力低下，極少有因斑馬奶水不足死亡的個體，此期間斑馬幼崽自身免疫力可能是導(dǎo)致其死亡的主要原因。近親的斑馬，如父與女，母與子，兄弟與姐妹間的繁育致使后代逐代親緣系數(shù)增高，這可能也是致使免疫缺陷凸顯的誘因之一。通常情況下，家庭中的雄性斑馬會驅(qū)逐單身雄性斑馬，以保證雌性普通斑馬有著良好的小環(huán)境孕育后代而生活不會受到騷擾[12]。本研究中繁殖前期(5歲前)死亡斑馬均為雄性，多為打斗和外傷致死；另一方面也是由于場地限制，沒有及時將幼年斑馬進行分群。當(dāng)斑馬形成穩(wěn)定種群后，死亡率便直線下降，在幼年斑馬接近繁殖年齡時，應(yīng)盡早分群，以避免更多的爭斗。

3.2 斑馬種群管理建議

鑒于國內(nèi)早期可追溯來源的斑馬幾乎均源于北京動物園繁育的后代，已不適宜在國內(nèi)動物園尋找更適宜的配對個體，引進新的種源迫在眉睫。由于飼養(yǎng)場地和引進新個體難度的限制，根據(jù)上述研究結(jié)果，我們引進2只無親緣關(guān)系的雄性斑馬，計算新增個體后的MSI配對表，合理進行配對選擇，以促進種群血緣調(diào)配。將原有種群進行重新組和，建立了2個斑馬種群，通過新個體的引進，至少可以在3～5 a內(nèi)維持斑馬種群遺傳多樣性在90%以上的可持續(xù)發(fā)展目標(biāo)。增加奠基者的繁殖機會，減少近親繁育，避免近親繁育，延長單一雌性斑馬世代時間，育齡雄性個體分群飼養(yǎng)，努力每5～10 a引進新的個體或進行合作繁殖。

致謝：感謝北京動物園張金國、盧巖、張耀華、王澤重對圈養(yǎng)斑馬種群管理工作的大力支持，感謝飼養(yǎng)隊各位同仁對相關(guān)數(shù)據(jù)和資料的提供。

[1] 周軍英.中國動物園白鶴圈養(yǎng)種群統(tǒng)計學(xué)與遺傳學(xué)分析[J].四川動物，2013，32(4)：492-497.

[2] 黃宏文，康明.保育遺傳學(xué)導(dǎo)論[M].北京：科學(xué)出版社，2005：154-269.

[3] Kleiman A，Thompson L.Wild mammals in captivity principles and techniques[M].Chicago：The University of Chicago Press，1996：277-281.

[4] 于澤英.動物園圈養(yǎng)種群的遺傳學(xué)管理[J].野生動物，2004，25(4)：41-42.

[5] Haig S M，Ballou J D.Pedigree analyses in wild populations[M]// Beissinger S R，McCullough D R.Population viability analysis.Chicago：University of Chicago Press，2002：388-405.

[6] 吳菡.上海動物園圈養(yǎng)東方白鸛種群動態(tài)研究[J].現(xiàn)代農(nóng)業(yè)科技，2013(12)：221-224.

[7] 周軍英.中國動物園冠鶴圈養(yǎng)種群現(xiàn)狀分析[J].野生動物，2013，34(2)：115-119.

[8] 周軍英，盧雁平，張金國.中國動物園圈養(yǎng)鶴類現(xiàn)狀調(diào)查及分析[J].野生動物，2012，33(4)：225-230.

[9] 廖輝，袁梨，殷毓中，等.中國動物園大猩猩圈養(yǎng)種群現(xiàn)狀分析[J].野生動物學(xué)報，2015，36(4)：378-381.

[10] 于澤英.川金絲猴圈養(yǎng)種群現(xiàn)狀分析[J].動物學(xué)雜志，2004，39(4)：45-49.

[11] 遲俊，徐暉，于明斌.半散放條件下普通斑馬的種群特征[J].獸類學(xué)報，2014，34(4)：419-422.

[12] 孫忻.斑馬簡介[J].生物學(xué)通報，2014，49(1)：8-11.

[13] Kleiman D G，Thompson K V，Baer C K.Wild mammals in captivity：principles and techniques for zoo management[M].Chicago：University of Chicago Press，2010：160-246.

[14] Frankham R.Effective population size/adult population size ratios in wildlife：a review[J].Genetical Research，1995，66(2)：95-107.

Plains zebra； Zoo； Captive population； Breeding； Genetic management

根據(jù)1971～2013年北京動物園圈養(yǎng)平原斑馬的數(shù)據(jù)，應(yīng)用種群管理軟件SPARKS v1.6和PMx v1.0對北京動物園斑馬圈養(yǎng)種群進行數(shù)據(jù)整理、統(tǒng)計學(xué)和遺傳學(xué)分析。結(jié)果顯示該種群數(shù)量呈現(xiàn)穩(wěn)定緩慢增長的狀態(tài)。目前，北京動物園斑馬種群面臨的問題主要為性別比例偏離了斑馬固有的一雄多雌的種群結(jié)構(gòu)，雄性個體較多；基因多樣性丟失較快。奠基者的后代數(shù)明顯不均衡，導(dǎo)致現(xiàn)今基因多樣性下降至83%，并存在大量半同胞間的近親繁殖，且現(xiàn)有可配對個體嚴重不足。1歲前是幼齡斑馬的死亡高峰期，流產(chǎn)、死胎為死亡主要原因，其次為外傷致死。建議在北京動物園引進新個體，建立2個健康的斑馬圈養(yǎng)種群，做好配對選擇及過多育齡雄性的分群，以減少競爭引發(fā)的死亡，增加潛在奠基者的繁殖機會，從而盡快擴大種群數(shù)量，并減少后代過多奠基者個體的繁殖機會。

Demographic and Genetics Analysis of a Captive Population ofCommon Zebra(Equus burchelli)at Beijing Zoo

You Yuyan1*Zhou Junying2Du Yang1Li Ying1Cao Changlei3

(1.Beijing Key Laboratory of Captive Wildlife Technologies，Beijing Zoo，Beijing，100044，China；2.Chinese Association of Zoological Gardens，Beijing，100037，China；3.Yangze Normal University，Chongqing，408100，China)

We conducted a data collation and demographic statistical and genetic analysis based on data collected on captive common zebra(Equusburchellii)at Beijing Zoo from 1971 to 2013.We used the population management software SPARKS v1.6 and PMx v1.0.The zebra population increased continuously.The main problem for management of the zebra population is that the sex ratio deviated yielding a surplus of males，and rapid loss of genetic diversity due to polygamy.The numbers of offspring from founders is unbalanced，which has led to a decline in genetic diversity to 83%.There are many inbred descendants among half-sib relatives and a seriously shortage of individuals suitable for breeding.The mortality peak of young zebras is <1 year of age.The main causes of mortality are abortion，fetal death，follow by traumatic death.We recommend introduction of new individuals to establish two healthy captive populations of zebra.Careful selection of individuals for mating and separation of excessive numbers of mature males could reduce death caused by competition while increasing the breeding opportunities for potential founders.These measures would enable expansion of the zebra population as soon as possible while reducing the breeding opportunities for founder whose descendants are already numerous.

國家自然科學(xué)基金項目( No． 30770361) 資助

由玉巖，女，35 歲，博士，高級工程師; 主要從事野生動物保護工作。

* 通訊作者: 由玉巖，E - mail: youyy351@163. com

2016-09-14

Q959.843

A

修回日期：2016-10-21

發(fā)表日期：2017-05-10

2310-1490(2017)02-222-06