白腹鷺CO Ⅰ和Cyt b基因序列分析及鷺屬系統(tǒng)發(fā)育研究

陳熙， 王斌， 羅旭*

(1.西南林業(yè)大學，西南地區(qū)生物多樣性保育國家林業(yè)局重點實驗室，昆明650224； 2. 云南高黎貢山國家級自然保護區(qū)瀘水管護分局，云南瀘水673100)

白腹鷺Ardeainsignis隸屬于鷺科Ardeidae鷺屬Ardea(Gill & Donsker，2014)，是一種體型較大的涉禽，高達127 cm，整體呈灰色，分布于尼泊爾、不丹、印度東北部和緬甸北部，我國云南西部(Hoyoetal. ，1992；韓聯(lián)憲等，2015)。自1994年，白腹鷺被世界自然保護聯(lián)盟(IUCN)列為瀕危(EN)物種，后升級為極危(CR) 物種，全球數(shù)量不足500只(IUCN，2016)。在《中國脊椎動物紅色名錄》中，白腹鷺的現(xiàn)狀是缺乏數(shù)據(jù)(DD)(蔣志剛等，2016)。對于該物種，不論是野外分布、種群現(xiàn)狀，還是系統(tǒng)分類、分子生物學研究，均近于空白。

關(guān)于鷺類的系統(tǒng)分類研究歷來不多，或者分類取樣十分有限(張保衛(wèi)等，2002；常青等，2003；張國萍等，2005)。鷺屬內(nèi)的物種歸屬也一直存在爭議。比如，大白鷺A.alba和中白鷺A.intermedia長期以來被置于白鷺屬Egretta(鄭作新，1964，1994；Hermannetal. ，1972；Benetal. ，1975)，但Sibley和Monoroe(1990)根據(jù)DNA雜交研究結(jié)果，將大白鷺歸入大白鷺屬Casmerodius，中白鷺則獨立成中白鷺屬Mesophoyx。但在國際鳥類聯(lián)合會(International Ornithologists Committee，IOC)的鳥類名錄和目前國內(nèi)普遍采用的鳥類分類系統(tǒng)中，大白鷺和中白鷺都置于鷺屬(鄭光美，2017)。對于白腹鷺，由于缺乏研究，中國鳥類分類系統(tǒng)一直沿用Peters的分類方式，將其置于鷺屬(Peters，1945；鄭作新，1964；鄭光美，2017)，但尚無分子系統(tǒng)學的驗證。

本研究根據(jù)IOC的分類系統(tǒng)，對鷺屬鳥類進行系統(tǒng)發(fā)育分析。分子標記采用線粒體DNA中的細胞色素b(cytochrome b，Cyt b)基因和細胞色素C氧化酶亞基Ⅰ(cytochrome oxidase subunit Ⅰ，COⅠ)基因，它們都被廣泛運用于鳥類分類與系統(tǒng)發(fā)育關(guān)系的研究(Mindell，1997；Weibel & Moore，2002；Hebertetal. ，2004；Webb & Moore，2005)。選擇鷺屬內(nèi)7種作為內(nèi)群，選擇鷺科中與鷺屬親緣關(guān)系近且有相應基因序列信息的其他6屬7種作為外群，構(gòu)建系統(tǒng)發(fā)育樹。

1 材料和方法

1.1 材料

2014年，在云南瀘水發(fā)現(xiàn)的1只白腹鷺是該物種在中國分布的第一條確切記錄(韓聯(lián)憲等，2015)。該個體于云南野生動物園進行救助，飼養(yǎng)過程中死亡，我們獲得少量實驗樣品。其余物種的Cytb和COⅠ基因序列均來自美國國家生物技術(shù)信息中心(National Center for Biotechnology Information，NCBI)(表1)。

表1 鷺科14種鳥類基因序列信息Table 1 List and gene information of 14 Ardeidae species

1.2 總DNA提取及基因擴增

取-80 ℃保存的肌肉樣品提取DNA。具體操作參考生工生物工程(上海)股份有限公司的Ezup柱式動物細胞/組織基因組DNA抽提試劑盒中使用說明書的標準抽提步驟。提取DNA后經(jīng)瓊脂糖凝膠電泳，在凝膠成像儀下觀察提取結(jié)果，若出現(xiàn)明顯條帶且與目的條帶長度相近，則說明已成功提取出樣品的總DNA，-20 ℃儲存?zhèn)溆谩?/p>

PCR反應體系為25 μL：5 μL DNA模板+12.5 μL Mix+5.5 μL ddH2O+2 μL引物(上下游引物各1 μL)，其中DNA模板和引物的密度均為10-5mol·L-1。反應在Applied Biosystems型PCR儀中進行。COⅠ基因片段擴增選用的引物對為BirdF1(TTCTCCAACCACAAAGACATTGGCAC)和BirdR1(ACGTGGGAGATAATTCCAAATCCTG)(Kerretal. ，2007)。PCR反應條件為：94 ℃ 4 min；94 ℃1 min，48 ℃ 1 min，72 ℃ 1 min，5個循環(huán)；94 ℃1 min，55 ℃ 1 min，72 ℃ 1 min，35個循環(huán)；72 ℃7 min。Cytb基因片段擴增選用的引物對為H16067(AGCCTTCAATCTTTGGCTTACAAG)和L14731(AATTGCATCCCACTTAATCGA)(Saetreetal. ，2001)。PCR反應條件為：94 ℃ 4 min；94 ℃ 1 min，48 ℃ 1 min， 72 ℃1.5 min，35個循環(huán)；72 ℃ 7 min。在凝膠成像儀下確認PCR產(chǎn)物是目標條帶后，送生工生物工程(上海)股份有限公司進行雙向測序。

1.3 DNA序列數(shù)據(jù)整理及分析

將測序結(jié)果用SEQMAN進行拼接并輔以人工校正。結(jié)合序列圖(abI格式文件)，剪去兩端信號較雜亂的序列，再查看2條鏈的重疊區(qū)，即在序列的內(nèi)部進行人工校正。一般情況下，為保證序列的準確性，重疊區(qū)域的堿基至少≥50 bp(黃原，1998)。將拼接校正過的DNA序列在NCBI上進行BLAST相似性搜索，以確保獲得的序列為目的基因片段。

運用Clustal W(Thompsonetal. ，1997)將所有物種的序列進行比對，再運用MEGA6.0對已經(jīng)比對好的序列計算堿基組成情況、保守位點、變異位點、簡約信息位點等，用Kimura 2-parameter(K2P)法計算遺傳距離。選取除鷺屬之外的其他6屬7種為外群，將2段基因合并，應用MrBayes3.1(Huelsenbeck & Ronquist，2003)構(gòu)建貝葉斯樹，運行4個馬爾可夫鏈，共運行5 000 000代，至平均標準離差(sd)=0.004 315(只有當sd<0.01時說明達到收斂，才可結(jié)束運行)。在FigTree 1.4.2中查看和修飾系統(tǒng)樹。

2 研究結(jié)果

2.1 白腹鷺基因序列特征

白腹鷺Cytb基因序列共有1 163個位點(表2)。堿基含量中G最高，C最低，A+T(51.7%)高于G+C(48.3%)。嘧啶含量(39.7%)明顯低于嘌呤含量(60.3%)。第一位點和第二位點的A+T均高于G+C，第三位點的A+T(48.0%)低于G+C(52.0%)。在密碼子的不同位點，堿基頻率有不同偏倚，第二位點上4種堿基含量較接近，偏倚程度最?。坏谌稽c上G含量增加到46.8%，C含量只有5.2%，偏倚最突出。

白腹鷺COⅠ基因序列共有722個位點(表2)。C含量最高，G含量最低。A+T(52.5%)高于G+C(47.5%)。嘧啶含量(57.0%)高于嘌呤含量(43.0%)。第一位點的A+T(41.5%)低于G+C(58.5%)，第二位點和第三位點的A+T均高于G+C。第一位點上堿基含量比較平均；第二位點上T高達41.0%，A只有14.9%；第三位點上A含量最高，G最少，堿基偏倚現(xiàn)象在第三位點最明顯。

2.2 鷺屬基因序列特征及種間遺傳距離

本研究獲得白腹鷺的Cytb基因片段，合并來自NCBI的另外5種鷺屬鳥類的同源區(qū)序列(NCBI中未查詢到黑冠白頸鷺A.cocoi的Cytb基因信息)。經(jīng)過比對排列后共有822個位點。除了中白鷺A+T的含量(49.3%)低于G+C的含量(50.7%)外，其余5種的A+T的含量均高于G+C的含量(表3)。保守位點650個，可變位點172個，簡約信息位點47個，單變異位點125個。物種間序列位點變異表現(xiàn)為轉(zhuǎn)換多于顛換，R(轉(zhuǎn)換/顛換)的平均值為3.46。

表2 白腹鷺Cyt b和CO Ⅰ基因序列的各位點堿基組成Table 2 Sequence composition of Cyt b and CO Ⅰ genes of Ardea insignis

本研究獲得白腹鷺的COⅠ基因片段，合并來自NCBI的另外6種鷺屬鳥類的同源區(qū)序列，經(jīng)過比對排列后共有599個位點。A+T的含量均高于G+C的含量(表3)。保守位點485個，可變位點111個，簡約信息位點60個，單變異位點51個。物種間序列位點變異表現(xiàn)為轉(zhuǎn)換多于顛換，R(轉(zhuǎn)換/顛換)的平均值為10.04。

基于COⅠ基因計算鷺屬各物種間的遺傳距離，白腹鷺與草鷺A.purpurea的遺傳距離最小(0.062)，與大白鷺的遺傳距離最大(0.093)(表4)。

表3 鷺屬7種的Cyt b和CO Ⅰ基因序列組成Table 3 Base composition of Cyt b andCO Ⅰ genes in the 7 Ardea species

表4 鷺屬7種的種間遺傳距離Table 4 Genetic distance between the 7 Ardea species

2.3 鷺屬系統(tǒng)進化樹

基于COⅠ和Cytb基因片段構(gòu)建貝葉斯系統(tǒng)發(fā)育樹顯示(圖1)，鷺屬鳥類構(gòu)成單系(支持率96%)。鷺屬7種鳥類中，最先分化的是大白鷺，其次是白腹鷺、草鷺。蒼鷺A.cinerea、大藍鷺A.herodias、中白鷺、黑冠白頸鷺聚為一支(支持率97%)。白腹鷺位于鷺屬分支內(nèi)，單獨聚為一支，無姐妹類群。

3 討論

3.1 白腹鷺的線粒體基因序列特征

線粒體基因序列的堿基組成研究較多，對74種鳥類線粒體基因組序列的分析表明，密碼子第二位的G、C含量波動較小，第三位的波動幅度大，這與第三位的C、T含量變化有關(guān)(高英凱等，2009)。白腹鷺Cytb和COⅠ基因序列的嘧啶和嘌呤含量有明顯差異，同樣是由第三位點上G、C含量變化引起，這與第三位點的變異通常不改變氨基酸編碼、受選擇壓力小有關(guān)(鐘東等，2002)。

在鳥類線粒體基因序列中，堿基組成存在一定的偏向性是普遍存在的現(xiàn)象(Arnaiz-Villenaetal. ，2001；Webb & Moore，2005)。偏倚情況可能因物種不同而不同，對24種雀類的研究表明，Cytb基因序列密碼子第一位點出現(xiàn)的頻率相似，第二位點G的頻率明顯低于其他3種堿基，第三位點偏倚G、T(Arnaiz-Villenaetal. ，2001)。在本研究中，Cytb基因序列的第一、二位點均無明顯偏倚，但第三位點明顯偏倚G，在該位點上的堿基偏倚與上述雀類研究相似。白腹鷺COⅠ基因序列的堿基偏倚和Cytb基因不同，第一位點沒有明顯的偏倚，第二位點偏倚T，第三位點偏倚A。

3.2 白腹鷺的系統(tǒng)分類地位

根據(jù)本文結(jié)果，鷺屬7種聚為一支，該節(jié)點有較高的支持率(96%)，支持IOC分類系統(tǒng)中把白腹鷺、大白鷺、中白鷺均放入鷺屬中的處理(Gill &Donsker，2014)，這與目前國內(nèi)最新的鳥類分類系統(tǒng)相同(鄭光美，2017)。在鷺屬中，白腹鷺是較早分化的物種，遺傳距離分析表明和白腹鷺親緣關(guān)系最近的是草鷺。然而，Rasmussen和Anderton(2005)提及大嘴鷺Ardeasumatrana外部形態(tài)特征和白腹鷺最為接近，但大嘴鷺的物種地位尚不明確，抑或是白腹鷺的某些地理種群產(chǎn)生了形態(tài)分化，因此，白腹鷺的系統(tǒng)分類地位還需進一步研究確認。

圖1 基于Cytb+COⅠ基因序列構(gòu)建的鷺科14種鳥類的貝葉斯樹
Fig.1 Bayes tree of 14 species from family Ardeidae based onCytbandCOⅠ genes

：

常青， 張保衛(wèi)， 金宏. 2003. 從12s RNA基因序列推測鷺類13種鳥類的系統(tǒng)發(fā)生關(guān)系[J]. 動物學報， 49(2)： 205-210.

高英凱， 苗正旺， 蘇小茜， 等. 2009. 74種鳥類線粒體基因組堿基組成及特征分析[J]. 云南農(nóng)業(yè)大學學報， 24(1)： 51-58.

韓聯(lián)憲， 何臣相， 王斌， 等. 2015. 云南發(fā)現(xiàn)白腹鷺[J]. 四川動物， 34(2)： 281.

黃原. 2012. 分子系統(tǒng)發(fā)生學[M]. 北京： 科學出版社.

蔣志剛， 江建平， 王躍招， 等. 2016. 中國脊椎動物紅色名錄[J]. 生物多樣性， 24(5)： 500-551.

張保衛(wèi)， 常青， 魏輔文. 2002. 鷺科鳥類分類及系統(tǒng)學研究進展[J]. 動物學雜志， 37(3)： 84-88.

張國萍， 王蔚， 朱世杰， 等. 2005. 鸛形目12種鳥類核c-mos基因和線粒體12s rRNA基因序列分析及其系統(tǒng)發(fā)生研究[J]. 四川動物， 24(4)： 500-506.

鄭光美. 2017. 中國鳥類分類與分布名錄(第三版)[M]. 北京： 科學出版社.

鄭作新. 1964. 中國鳥類系統(tǒng)檢索[M]. 北京： 科學出版社： 9-12.

鄭作新. 1994. 中國鳥類種和亞種分類名錄大全[M]. 北京： 科學出版社： 9-11.

鐘東， 趙貴軍， 張振書， 等. 2002. 基因組內(nèi)堿基分布整體均衡與局部不均衡的研究進展[J]. 遺傳， 24(3)： 351-355.

Arnaiz-Villena A， Guillén J， Ruiz-del-Valle V，etal. 2001. Phylogeography of crossbills， bullfinches， grosbeaks， and rosefinches[J]. Cellular and Molecular Life Sciences， 58(8)： 1159-1166.

Ben FK， Edward CD， Martin WW. 1975. A field guide to the birds of south-east Asia[M]. Boston： Houghton Mifflin Press： 41-47.

Gill F， Donsker D. 2014. IOC world bird list (v 4.3)[EB/OL]. [2017-11-07]. DOI： 10.14344/IOC.ML.4.3.

Hebert PDN， Stoeckle MY， Zemlak TS，etal. 2004. Identification of birds through DNA barcodes[J]. PLoS Biology， 2(10)： 1657-1663.

Hermann H， Richard F， John P. 1972. The birds of Britain and Europe with north African and the middle East[M]. London： Collins Press： 34-39.

Hoyo JD， Elliott A， Sargatal J. 1992. Handbook of the birds of the world vol. 5[M]. Spain： Lynx Edicions Press： 407.

Huelsenbeck JP， Ronquist F. 2003. MrBayes： Bayesian inference of phylogenetic trees[J]. Bioinformatics， 17： 754-755.

IUCN. 2016. The IUCN red list of threatened species[EB/OL]. (2016-10-01)[2017-11-07]. http：//www.iucnredlist.org.

Kerr KCR， Stoeckle MY， Dove CJ，etal. 2007. Comprehensive DNA barcode coverage of north American birds[J]. Molecular Ecology Notes， 7(4)： 535-543.

Mindell DP. 1997. Avian molecular evolution and systematics[M]. San Diego： Academic Press： 84-86.

Peters JL. 1945. Checklist of birds of the world vol. 5[M]. Cambridge， MA： Harvard University Press.

Rasmussen PC， Anderton JC. 2005. Birds of south Asia： the ripley guide[M]. Washington， Barcelona： Smithsonian Institution and Lynx Edicions.

Saetre GP， Borge T， Lindell J，etal. 2001. Speciation， introgressive hybridization and nonlinear rate of molecular evolution in flycatchers[J]. Molecular Ecology， 10(3)： 737-749.

Sibley CG， Monoroe BLJ. 1990. Distribution and taxonomy of birds of the world[M]. New Haven & London： Yale University Press： 302-310.

Thompson JD， Gibson TJ， Plewniak F，etal. 1997. The Clustal X windows interface： flexible strategies for multiple sequence alignment aided by quality analysis tools[J]. Nucleic Acids Resreach， 24： 4876-4882.

Webb DM， Moore WS. 2005. A phylogenetic analysis of woodpeckers and their allies using 12S， Cyt b and COⅠ nucleotide sequences (Class Aves； Order Piciformes)[J]. Molecular Phylogenetics Evolution， 36： 233-248.

Weibel AC， Moore WS. 2002. Molecular phylogeny of a cosmopolitan group of woodpeckers (GenusPicoides) based onCOⅠ andCytbmitochondrial gene sequences[J]. Molecular Phylogenetics Evolution， 22： 65-75.