加入放歸個體后的栗子坪保護(hù)區(qū)大熊貓小種群遺傳結(jié)構(gòu)研究

陳子君， 黃蜂，2， 何可, 馮菲菲,古曉東,楊旭煜,戴強(qiáng), 齊敦武,楊志松

(1.西華師范大學(xué) 西南野生動植物資源保護(hù)教育部重點實驗室， 四川　南充　637002；2.四川栗子坪國家級自然保護(hù)區(qū)管理處， 四川　石棉　625400；3.四川省野生動物資源調(diào)查保護(hù)管理站， 四川　成都　610081；4.中國科學(xué)院成都生物研究所， 四川　成都　610041；5.成都大熊貓繁育研究基地, 四川　成都　610086)

?

加入放歸個體后的栗子坪保護(hù)區(qū)大熊貓小種群遺傳結(jié)構(gòu)研究

陳子君1， 黃蜂1，2， 何可1, 馮菲菲1,古曉東3,楊旭煜3,戴強(qiáng)4, 齊敦武5,楊志松1

(1.西華師范大學(xué) 西南野生動植物資源保護(hù)教育部重點實驗室， 四川南充637002；2.四川栗子坪國家級自然保護(hù)區(qū)管理處， 四川石棉625400；3.四川省野生動物資源調(diào)查保護(hù)管理站， 四川成都610081；4.中國科學(xué)院成都生物研究所， 四川成都610041；5.成都大熊貓繁育研究基地, 四川成都610086)

小相嶺山系是現(xiàn)有大熊貓種群數(shù)量最少的山系之一。國家林業(yè)局?jǐn)M通過在該山系的栗子坪國家級自然保護(hù)區(qū)放歸大熊貓“瀘欣”、“淘淘”、“張想”來復(fù)壯大熊貓小種群。為評估加入3只放歸個體后的小種群遺傳結(jié)構(gòu)，2012～2013年在研究區(qū)域采集大熊貓糞便樣品92份，鑒定出14只個體(含放歸個體)。放歸個體加入后，測得放歸個體與本地大熊貓小種群親緣關(guān)系較遠(yuǎn)。種群年齡組成為青年組8只、成年組4只、未知年齡個體2只。平均觀測雜合度(Ho)為0.557，平均期望雜合度(He)為0.581，遺傳多樣性處于中等偏上水平。5個位點上共計增加6個等位基因，平均等位基因數(shù)(MNA)達(dá)到3.667。遺傳基因簇由1簇增加為4簇，出現(xiàn)新的基因頻率。以上結(jié)果顯示放歸個體加入后該區(qū)大熊貓小種群遺傳多樣性有所提高，放歸成效初顯。

大熊貓；小種群；遺傳多樣性；放歸；栗子坪國家級自然保護(hù)區(qū)

放歸是促進(jìn)小種群復(fù)壯，恢復(fù)破碎化棲息地的有效方法[1-2]。國內(nèi)外已有不少成功案例，國內(nèi)的如普氏野馬(Equus ferus przewalskii)、麋鹿(Elaphurus davidianus )、揚子鱷(Alligator sinensis)等。國外的如加拿大盤羊(Ovis canadensis)、美洲野牛(Bison bison)、美洲黑熊(Ursus americanus)等。然而，也有少數(shù)放歸活動因近交[3]、遺傳多樣性丟失[4]或者人類干擾、天敵捕食[5]而成為失敗的嘗試。因此，及時分析放歸成效，總結(jié)經(jīng)驗、教訓(xùn)，對放歸工作的進(jìn)一步開展有著極為重要的意義。

小相嶺山系是現(xiàn)有大熊貓種群數(shù)量最少的山系之一，種群生存力分析表明，如果與其他種群沒有遺傳交流或者無外來個體補(bǔ)充，該小種群極有可能在80年內(nèi)會趨于滅絕[6]。為提高該區(qū)大熊貓小種群生存能力，促進(jìn)大熊貓小種群的復(fù)壯，(截止到本研究開展的2013年)國家林業(yè)局在小相嶺山系的栗子坪國家級保護(hù)區(qū)先后放歸了3只大熊貓(Ailuropoda melanoleuca)：“瀘欣”已放歸超過6年，“淘淘”已放歸超過3年，“張想”也已放歸超過2年。然而，加入放歸個體后，栗子坪自然保護(hù)區(qū)大熊貓小種群遺傳結(jié)構(gòu)是否有所改善，目前仍不清楚。本研究利用野外調(diào)查收集的大熊貓糞便，通過糞便咬節(jié)測量、糞便DNA性別鑒定以及微衛(wèi)星分析，確定了放歸地大熊貓種群的年齡結(jié)構(gòu)、性比和種群遺傳結(jié)構(gòu)，并從種群遺傳上對放歸成效進(jìn)行了評估。這一評估結(jié)果有助于確定將來大熊貓放歸和小種群復(fù)壯項目的方向和重點，對其他瀕危珍稀野生動物的小種群復(fù)壯工作也有很好的參考價值。

1　研究對象及區(qū)域

圖1　研究區(qū)域區(qū)位示意圖

栗子坪自然保護(hù)區(qū)位于四川石棉縣境內(nèi)，地理位置為E 102°10′33″～102°29′07″，N 28°51′02″～29°08′42″；南北長23 km，東西寬17.8 km，總面積47 940 hm2。區(qū)內(nèi)最高海拔為4 551 m，最低海拔1 330 m，最大相對高差達(dá)3 221 m，平均相對高差大于2 000 m。保護(hù)區(qū)是小相嶺山系大熊貓種群的核心分布區(qū)和關(guān)鍵棲息地，也是小相嶺山系大熊貓種群遺傳基因交流的關(guān)鍵走廊帶和重要集中分布區(qū)。

本項目研究區(qū)域位于栗子坪自然保護(hù)區(qū)內(nèi)公益海保護(hù)站周邊，面積約為23 251.45 hm2?！盀o欣”、“淘淘”、“張想”先后放歸于該區(qū)(詳見圖1和表1)，遺傳分析所用糞便也采集自這一區(qū)域。這一區(qū)域的大熊貓棲息地與周邊大熊貓棲息地存在較為嚴(yán)重的隔離[7-8]，第四次全國大熊貓調(diào)查也認(rèn)為該區(qū)域大熊貓種群是一個相對獨立的局域種群[9]。

表1　放歸大熊貓個體信息

2　研究方法

2．1野外采樣

野外樣品于2012年4月至2013年12月采集自保護(hù)區(qū)內(nèi)公益海保護(hù)站周邊(見圖1)。 研究區(qū)域內(nèi)調(diào)查樣線的抽樣強(qiáng)度為1 條/200 hm2,采用非損傷性遺傳取樣方法在野外收集糞便樣品[10]。 在采樣同時，對所有的樣品采集點均進(jìn)行GPS定位并記錄。

2．2大熊貓年齡組劃分

測量采集到的大熊貓糞便的咬節(jié)長度，采用“離差平方和法”對咬節(jié)長度進(jìn)行聚類，參照胡錦矗的年齡組劃分標(biāo)準(zhǔn)[11]確定個體年齡段：幼年組(小于1.5 歲)，咬節(jié)值小于30 mm，糞便直徑30 mm以下；青年組(大于等于1.5歲，并且小于5.5 歲),咬節(jié)值在30～35 mm(包含30 mm) 之間, 糞便直徑30～55 mm(包含30 mm)；成年組(大于等于5.5歲，并且小于20 歲)，咬節(jié)長度在35～40 mm(包含35 mm)之間，糞便直徑55～65 mm (包含55 mm)；老年組(大于等于20 歲)，咬節(jié)長度在40 mm(包含40 mm)以上，糞便直徑在65 mm(包含65 mm)以上。

2．3糞便微衛(wèi)星DNA提取與分析

使用QIAGEN公司的QIAamp DNA Stool kit試劑盒進(jìn)行大熊貓糞便DNA的提取，性別鑒定PCR、微衛(wèi)星PCR以及結(jié)果判定參考胡義波[12]、詹祥江[10]的方法，我們挑選了12對微衛(wèi)星引物(Ame -μ10, Ame -μ11, Ame -μ22, Ame -μ13, Ame -μ15, Ame -μ24, Ame -μ26, Ame -μ27以及AY161179, AY161195, AY161213, AY161217)來對大熊貓糞便DNA進(jìn)行微衛(wèi)星擴(kuò)增[12]。

因大熊貓糞便中DNA含量較低，為避免因引物選取、PCR擴(kuò)增、數(shù)據(jù)讀取等引起的誤差，我們利用MICRO-CHECKER軟件[13]檢測基因分型數(shù)據(jù)中是否存在無效等位基因(Null alleles)、大片段等位基因丟失(Large allele dropout)及鬼影峰誤判(Stutter scoring error)等分型錯誤。利用GIMLET 1．3．3軟件[14]計算了12個座位的聯(lián)合P(ID)值，以確定選擇的微衛(wèi)星引物數(shù)量是否能將影子效應(yīng)(Shadow effect)[15]減少到一個不顯著的水平。

2．4個體鑒定與種群數(shù)量估計

利用MStools插件[16]，尋找基因分型數(shù)據(jù)中相匹配的基因型。判斷為同一個個體的標(biāo)準(zhǔn)是：所有座位上的基因型都相同，或者只有一個座位上的一個等位基因存在差別[17]。由于我們采樣分析是針對放歸成效的研究，糞便采集也集中在放歸大熊貓主要活動的區(qū)域(依據(jù)放歸大熊貓佩戴的GPS項圈確定)，可以認(rèn)為出現(xiàn)在此區(qū)域外或是很少出現(xiàn)在該區(qū)域的個體被捕捉概率很低，因此我們選擇CAPWIRE軟件里TIRM(Two Innate Rates Model)模型[10]進(jìn)行種群數(shù)量估計。

2．5種群遺傳多樣性及親緣關(guān)系分析

利用ARLEQUIN 3．5[18]和GENEPOP 3．4[19]計算期望雜合度(He, expected heterozygosity)、表觀雜合度(Ho, observed heterozygosity)、等位基因數(shù)(A)，測算總體和各個位點是否符合哈溫平衡(Hardy-Weinberg)。利用FSTAT 2．9．3．2[20]來檢驗位點是否存在連鎖不平衡(Linkage disequilibrium)。利用STRUCTURE方法[21]檢測研究區(qū)域所有個體中存在的遺傳基因簇。利用Kingroup v2．08[22]來測算個體間的親緣關(guān)系遠(yuǎn)近，其中個體間的親緣關(guān)系值(pairwise relatedness, r)參考Wang[23]的研究。

3　結(jié)果

3．1種群數(shù)量、性比及年齡結(jié)構(gòu)

本研究從92 份樣品(含糞便樣品及放歸個體血樣)中得到52個可靠基因型。TIRM模型分析表明，研究區(qū)域大熊貓種群數(shù)量(含3 只放歸個體)為25 只(95%置信范圍：14～28)。種群數(shù)量最低估計值(lowCI)為14 只(與實驗室分析鑒定結(jié)果一致)，因此研究區(qū)域中大熊貓種群最保守的數(shù)量應(yīng)該是14 只。14只個體分布在公益海保護(hù)站區(qū)域，海拔集中在2 291～3 510 m(見圖2)，分別是8 只雌性、6只雄性，雌雄比例為4∶3。其中“瀘欣”、“淘淘”和“張想”為放歸個體，年齡有準(zhǔn)確記錄，無需通過糞便咬節(jié)確定年齡，另有兩只個體未能獲得有效咬節(jié)數(shù)據(jù)，無法判斷年齡，其余個體的年齡段均根據(jù)咬節(jié)法測得。結(jié)果表明，本地種群的年齡組成為：青年組6 只，成年組3 只，未知年齡個體2只(見表2)。從年齡結(jié)構(gòu)上來看，青年組個體數(shù)明顯多于成年組個體數(shù)，年齡組成為增長型。 加入放歸個體后的大熊貓種群年齡組成為：青年組8 只、成年組4只，未知年齡個體2只(見表2)，與放歸前比較，年齡組結(jié)構(gòu)的增長型特征更為明顯。

表2　本研究大熊貓種群年齡組組成

3．2本地小種群遺傳結(jié)構(gòu)

MICRO-CHECKER檢驗表明，本研究基因分型數(shù)據(jù)結(jié)果未受無效等位基因、大片段等位基因丟失及鬼影帶誤判等分型錯誤的顯著干擾。GIMLET分析顯示，12 個微衛(wèi)星座位的聯(lián)合區(qū)分率很高，兩只個體具有相同基因型的無偏概率P(ID)為3．286×10-9，即使為雙胞胎，判斷錯誤的概率P(sib) 也僅為0．04778%，小于0．01[15]。哈溫平衡檢驗結(jié)果顯示研究區(qū)域本地大熊貓小種群未偏離哈溫平衡(P>0．05)。連鎖不平衡分析結(jié)果顯示本研究采用的12 個位點，對當(dāng)?shù)卮笮茇埿》N群而言不存在連鎖遺傳。12個微衛(wèi)星位點分析結(jié)果顯示研究區(qū)域本地大熊貓小種群的平均觀測雜合度(Ho)為0．635，平均期望雜合度(He)為0．608。平均等位基因數(shù)(MNA)為3．455(見表3)。STRUCTURE 分析結(jié)果顯示K=1時(圖3A)，LnP(D)具有最大的平均值，圖3B表明研究區(qū)域本地大熊貓小種群整體上是一個遺傳基因簇。

表3　研究區(qū)域本地大熊貓小種群(代號A)與加入放歸個體后的種群(代號B)遺傳多樣性信息

圖3　LnP(D)值及遺傳基因簇圖圖A:本地小種群LnP(D)值；圖C:放歸個體加入后的種群LnP(D)值；圖B圖D分別為K=1、K=4的遺傳基因簇圖

3．3加入放歸個體后小種群遺傳結(jié)構(gòu)

3．3．1放歸個體和本地大熊貓親緣關(guān)系

親緣關(guān)系分析結(jié)果顯示(見表4)，“瀘欣”、“淘淘”、“張想”和整個本地小種群的親緣關(guān)系都較遠(yuǎn)，與其它13 只個體的親緣系數(shù)(r)平均值平均為-0．07、-0．21、-0．23。

雖然“瀘欣”和lzp6、lzp7、lzp10、lzp11、“張想”的親緣系數(shù)(r)>0，但是這五只個體均為雌性，不存在因親緣關(guān)系近而帶來的近交風(fēng)險?！疤蕴浴焙蚻zp2、lzp7、“張想”的親緣系數(shù)(r)分別為0．05、0．10、0．10，但是r<0．125，“淘淘”和這三只個體的親緣關(guān)系水平低于雙重表親(Double first cousins)，不存在明顯的因親緣關(guān)系較近而帶來的近交風(fēng)險。“張想”和lzp6、“淘淘”、“瀘欣”的親緣系數(shù)(r)分別為0．03、0．10、0．11，r<0．125，“張想”和這三只個體的親緣關(guān)系水平也低于雙重表親(Double first cousins)，不存在明顯的因親緣關(guān)系近而帶來的近交風(fēng)險。因lzp7是雌性大熊貓，雖然它和“張想”親緣系數(shù)(r)為0．30，屬于親子/全同胞(parent-offspring/full sibs)的水平，但也不存在因親緣關(guān)系近的近交風(fēng)險。

表4　研究區(qū)域大熊貓個體親緣關(guān)系

3．3．2加入放歸個體后的遺傳結(jié)構(gòu)

根據(jù)哈溫平衡檢驗分析結(jié)果，加入放歸個體后的大熊貓小種群仍未偏離哈溫平衡(p>0．05)。連鎖不平衡分析結(jié)果顯示對于加入放歸個體后的大熊貓小種群，也不存在連鎖遺傳。12個微衛(wèi)星位點分析結(jié)果顯示：加入放歸個體后種群平均觀測雜合度(Ho)為0．557，平均期望雜合度(He)為0．581(見表3)。相比Zhang等[24]的研究，加入放歸個體后的遺傳多樣性仍處于中等偏上的水平。平均等位基因數(shù)(MNA)增加到3．667，表現(xiàn)在Ame-μ11(Locus1)、Ame-μ22(Locus3)、Ame-μ24(Locus7)、Ame-μ13(Locus10)和Ame-μ27(Locus12)5個位點上，共計增加6 個等位基因。 STRUCTURE分析結(jié)果表明，K=4時(見圖3C)，LnP(D)具有最大的平均值，遺傳基因簇由一簇變?yōu)榱怂拇?見圖3D)，這表明放歸個體加入后，出現(xiàn)了新的基因頻率。

4　討 論

三只放歸個體加入當(dāng)?shù)胤N群后，研究區(qū)域大熊貓小種群的等位基因數(shù)有所增加，同時也出現(xiàn)了新的基因頻率，這表明該區(qū)大熊貓小種群的遺傳多樣性確實呈增長趨勢。遺傳多樣性的高低在一定程度上反映了大熊貓種群的進(jìn)化潛力[25]，因此放歸個體的加入有助于研究區(qū)域大熊貓小種群的種群遺傳多樣性的提高，將對該種群的復(fù)壯帶來積極影響。

本研究分析結(jié)果顯示放歸個體與研究區(qū)域本地大熊貓小種群親緣關(guān)系較遠(yuǎn)。楊波等[26]的研究表明：親緣關(guān)系越遠(yuǎn)，釋放個體與本地種群組成的混合種群內(nèi)近交率會越低，從而減少近交衰退的風(fēng)險，更有利于改善野生種群遺傳水平。因此這三只放歸個體的加入，能有效提高研究區(qū)域本地大熊貓小種群的遺傳多樣性水平。

迄今，“瀘欣”、“淘淘”、“張想”已放歸分別超過6年、3年、2年，它們在大熊貓放歸地——栗子坪國家級自然保護(hù)區(qū)生存狀況良好。年齡組結(jié)構(gòu)分析顯示，加入放歸個體后，研究區(qū)域大熊貓小種群年齡結(jié)構(gòu)得到進(jìn)一步改善。此外，該區(qū)域的紅外相機(jī)監(jiān)測結(jié)果顯示，“瀘欣”已成功帶崽，表明放歸個體已經(jīng)初步融入當(dāng)?shù)胤N群，放歸工作取得初步成效。若其余兩只放歸個體也能成功繁殖，將進(jìn)一步促進(jìn)該地區(qū)大熊貓小種群的復(fù)壯。

目前，小相嶺山系大熊貓種群數(shù)量仍然很少，成功放歸三只大熊貓遠(yuǎn)不足以復(fù)壯當(dāng)?shù)匦》N群。針對這一問題，該山系還需進(jìn)一步持續(xù)開展大熊貓放歸工作。而在放歸個體的選擇上，需要比較放歸個體和當(dāng)?shù)匦》N群的遺傳背景，優(yōu)先選擇能有效提高當(dāng)?shù)匦》N群遺傳多樣性的個體。

[1]ZHU L F, ZHAN X J, MENG T, et al． Landscape features influence gene flow as measured by cost-distance and genetic analyses: a case study for giant pandas in the Daxiangling and Xiaoxiangling mountains [J]． BMC genetics, 2010(1):72．

[2]ZHU L F, ZHANG S N, GU X D, et al． Significant genetic boundaries and spatial dynamics of giant pandas occupying fragmented habitat across southwest China [J]． Molecular Ecology, 2011, 20(6):1122-1132．

[3]BIEBACH I, KELLER L F． Inbreeding in reintroduced populations: the effects of early reintroduction history and contemporary processes [J]． Conservation Genetics, 2010(2):527-538．

[4]ALCAIDE M．, NEGRO J J．, SERRANO D, et al． Captive breeding and reintroduction of the lesser kestrel Falco naumanni: a genetic analysis using microsatellites [J]． Conservation Genetics, 2010(11): 331-338．

[5]GREY-ROSS R, DOWNS C T, KIRKMAN K．． Reintroduction Failure of Captive-Bred Oribi (Ourebia ourebi) [J]． South African Journal of Wildlife Research, 2009, 39(1): 34-38．

[6]朱磊, 吳攀文, 張洪, 等． 小相嶺山系大熊貓種群生存力分析 [J]． 西華師范大學(xué)學(xué)報( 自然科學(xué)版), 2008, 29(2): 112-116．

[7]ZHU L F, ZHAN X J, WU H, et al． Conservation implications of drastic reductions in the smallest and most isolated populations of giant pandas [J]． Conservation Biology, 2010, 24(5): 1299-1306．

[8]青菁, 胥池, 楊彪, 等． 小相嶺山系大熊貓廊道規(guī)劃[J]． 生態(tài)學(xué)報, 2016, 36(4): 1-9．

[9]四川省林業(yè)廳． 四川的大熊貓[M]． 成都：四川科學(xué)技術(shù)出版社，2015.

[10]詹祥江．利用非損傷性遺傳取樣研究大熊貓的種群數(shù)量和擴(kuò)散模式 [D]． 中國科學(xué)院動物研究院, 2006．

[11]胡錦矗． 從野外大熊貓的糞便估計年齡及其種群年齡結(jié)構(gòu)的研究 [J]． 獸類學(xué)報, 1987, 7 (2): 81-84．

[12]胡義波． 涼山山系大熊貓的種群歷史、景觀與空間遺傳格局研究 [D]． 北京：中國科學(xué)院研究生院, 2008．

[13]VAN OOSTERHOUT C, HUTCHINSON W F, WILLS D P M, et al． Micro-checker: software for identifying and correcting genotyping errors in microsatellite data [J]． Molecular Ecology Notes, 2004(3): 535-538．

[14]VALIERE N． GIMLET: a computer program for analysing genetic individual identification data [J]． Molecular Ecology Notes, 2002, 2: 377-379．

[15]MILLS L S, CITTA J J, LAIR K P, et al． Estimating animal abundance using noninvasive DNA sampling: Promise and pitfalls [J]． Ecological Applications, 2000, 10(1):283-294．

[16]PARK S D E． Trypanotolerance in west African cattle and the population genetic effects of selection [D]． Dublin：University of Dublin, 2001．

[17]BELLEMAIN E, SWENSON J E, TALLMON D, et al． Estimating population size of elusive animals with DNA from hunter-collected feces: Four methods for brown bears [J]． Conservation Biology, 2005, 19(1): 150-161．

[18]EXCOFFIER L, LAVAL G, SCHNEIDER S． Arlequin (version 3．0): an integrated software package for population genetics data analysis [J]． Evolutionary Bioinformatics Online, 2005, 1: 47-50．

[19]RAYMOND M, ROUSSET F． Genepop (Version1．2): Population Genetics Software for Exact Tests and Ecumenicism [J]． Journal of Heredity, 1995,86(3):248-249．

[20]GOUDET J． FSTAT, a program to estimate and test gene diversities and fixationindices (version 2．9．3)． 2001． Updated from Goudet J． FSTAT (Version 1．2): A computer program to calculate F-statistics [J]． Journal of Heredity, 1995, 86: 485-486．

[21]PRITCHARD J K, STEPHENS M, DONNELLY P． Inference of population structure using multilocus genotype data [J]． Genetics, 2000, 155(2): 945-959．

[22]KONOVALOV D A, MANNING C, HENSHAW M T． Kingrouf: a program for pedigree relationship reconstruction and kin group assignments using genetic markers [J]． Molecular Ecology Notes, 2004(4): 779-782．

[23]WANG J． COANCESTRY: a program for simulating, estimating and analysing relatedness and inbreeding coefficients [J]． Molecular Ecology Resources, 2011, 11(1): 141-145．

[24]ZHANG B W, LI M, ZHANG Z J, et al． Genetic Viability and Population History of the Giant Panda, Putting an end to the "Evolutionary dead end"? [J]． Molecular Biology and Evolution, 2007, 24(18): 1801-1810．

[25]張澤鈞, 張陜寧, 魏輔文, 等． 移地與圈養(yǎng)大熊貓野外放歸的探討 [J]． 獸類學(xué)報, 2006, 26 (3): 292-299．

[26]楊波, 楊承忠, 涂飛云, 等． 大熊貓野化放歸中的遺傳學(xué)分析 [J]． 四川動物, 2013, 32(1):149-155．

[Abstract]The population of giant panda in Xiaoxiangling is one of the smallest . In order to rejuvenate this small panda population, the State Forestry Administration has tried to release three giant pandas(Ailuropoda melanoleuca)： “Luxin”, “Taotao” and “Zhangxiang” in the Liziping National Nature Reserve. To assess the genetic structure of this small population after reintroduction, 92 stools of giant pandas were collected in research area from 2012 to 2013, and 14 individuals(including three reintroduced individuals)were identified. After reintroduced the three individuals, measurement proved that the kinship between the reintroduced individuals and local small population was far. The population age composition was: eight individuals in youth group，four individuals in adult group and two individuals in unknown group.The mean observed heterozygosity (Ho) was 0.557. The mean expected heterozygosity (He) was 0.581. Genetic diversity was in the medium level.Five sites increased six alleles.And the mean allele number (MNA) achieved 3.667. What’s more ,the genetic cluster increased to four meaning that new gene frequency had appeared. All these results suggest that the genetic diversity of the small giantpanda population in this region is improved. The reintroduction is beginning to bear results.

[Key words]giant panda(Ailuropoda melanoleuca)； small population； genetic diversity； reintroduction； Liziping National Nature Reserve

[責(zé)任編輯劉景平]

The Genetic Structure of Small Giant Panda Population in the Liziping Nature Reserve after Reintroduction

CHEN Zi-jun1, HUANG Feng1,2, HE Ke1, FENG Fei-fei1, GU Xiao-dong3, YANG Xu-yu3,DAI Qiang4, QI Dun-wu5, YANG Zhi-song1

(1.Key Laboratory of Southwest China Wildlife Resource Conservation, China West Normal University,Nanchong, Sichuan 637002, China;2.Sichuan Liziping Nature Reserve Management Office, Shimian, Sichuan 625400, China;3.Sichuan Station of Wildlife Survey and Management, Chengdu, Sichuan 610081, China;4.Chengdu Institute of Biology, Chinese Academy of Sciences, Chengdu, Sichuan 610041, China;5.Chengdu Research Base of Giant Panda Breeding, Chengdu, Sichuan 610086, China)

Q953

A

1672-9021(2016)02-0001-07

陳子君(1989-)，女，四川綿陽人，西華師范大學(xué)生命科學(xué)學(xué)院碩士研究生，主要研究方向：大熊貓保護(hù)生物學(xué)；楊志松(1975-)，男，四川南充人，博士，西華師范大學(xué)生命科學(xué)學(xué)院副教授，碩士生導(dǎo)師，主要研究方向：大熊貓保護(hù)生物學(xué)、分子系統(tǒng)地理學(xué)、 鳥類生態(tài)學(xué)。

香港海洋公園保育基金資助項目(GP07_1213)。

2016-03-01