王湘君 陳 文 林熾賢 林杰曼

(1.海南熱帶海洋學(xué)院生命科學(xué)與生態(tài)學(xué)院，兩棲爬行動(dòng)物研究省級(jí)重點(diǎn)實(shí)驗(yàn)室，三亞，572022；2.南京師范大學(xué)生命科學(xué)學(xué)院，南京，201146；3.海南熱帶海洋學(xué)院海洋科學(xué)技術(shù)學(xué)院，三亞，572022)

白唇竹葉青蛇(Trimeresurusalbolabris)，白唇腺蝮隸屬于蛇亞目(Serpentes)蝰科(Viperiae)竹葉青蛇屬，棲息于平原、丘陵或低山區(qū)，各種水域附近的灌叢和雜草上。國內(nèi)分布于澳門、福建、廣東、廣西、貴州、海南、湖南、江西、臺(tái)灣、香港、云南；國外分布于尼泊爾、印度北部、緬甸、泰國、印度支那、馬來西亞、印度尼西亞[1-2]。

Avise首次提出把動(dòng)物學(xué)、遺傳學(xué)、系統(tǒng)學(xué)、地理學(xué)和分子生物學(xué)的理論和方法相結(jié)合，構(gòu)建“種內(nèi)系統(tǒng)發(fā)育地理學(xué)”。其方法論基礎(chǔ)是動(dòng)物線粒體DNA單倍型的系統(tǒng)發(fā)育關(guān)系及其群體動(dòng)力學(xué)[3-5]。運(yùn)用親緣地理學(xué)的理論與方法，揭示種群的歷史，例如：發(fā)生、建群、發(fā)展、隔離、分化、遷移、擴(kuò)散、瓶頸和滅絕；在種內(nèi)水平上，闡述了種群演化歷程，及其與地理學(xué)的關(guān)系。其中，12S、16S、Cytb、ND2、ND4、D-loop、COI等序列，是分子標(biāo)記上運(yùn)用比較廣的幾個(gè)基因序列[6-10]。

作者通過研究白唇竹葉青蛇12S、16S、Cytb、ND4四個(gè)線粒體基因部分序列，分析諸多單倍型序列彼此之間的親緣以及變異程度關(guān)系，推斷竹葉青蛇的擴(kuò)散方向和地理分化中心。

1 材料與方法

1.1 實(shí)驗(yàn)材料

1.1.1 PCR實(shí)驗(yàn)所用標(biāo)本

采集來自海南產(chǎn)白唇竹葉青蛇共12個(gè)個(gè)體，其中3個(gè)采自樂東縣(西南部) 編號(hào)為C01，C02，C03；2個(gè)采自東方市(西部)編號(hào)為C04，C05；4個(gè)采自昌江縣(西北北部)編號(hào)為C06，C07，C08，C09；采自三亞市(南部)、陵水縣(東南部)、瓊海市(東部)各1個(gè)，分別編號(hào)為C10，C11，C12。

1.1.2 從網(wǎng)上下載的竹葉青蛇序列

表1 GenBank下載的序列

Tab.1 The sequences downloaded from GeneBank

1.2 實(shí)驗(yàn)方法

采集標(biāo)本運(yùn)用PCR法分別擴(kuò)增12S、16S、Cytb、ND4四個(gè)基因片段，并測(cè)列。采用通用引物，引物序列和PCR擴(kuò)增條件如下表。

表2 引物設(shè)計(jì)圖表

Tab.2 Design of primer

表3 PCR反應(yīng)條件

Tab.3 Reaction requirement for PCR

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

利用12S、16S、Cytb、ND4基因序列的單倍型判別相對(duì)原始及進(jìn)化的序列來確定遷徙擴(kuò)散路線。采用DNASTAR軟件ClustalX程序排列DNA同源序列后，人工逐項(xiàng)核查，將全部序列都視為無序特征。再將12S、16S、Cytb、ND4基因序列片段的seq文件分別加載到MAGA6.0軟件后，根據(jù)其堿基突變位點(diǎn)數(shù)逐步比較排查，從而篩選出較為原始的單倍型序列，作為起源中心。

2 結(jié)果與分析

2.1 基因序列的核苷酸組成

表4 竹葉青蛇屬蛇類基因序列的核苷酸組成(bp)

Tab.4 Nucleotide composition of genus Trimeresurus sequence(bp)

從測(cè)定結(jié)果可以看出，白唇竹葉青蛇的堿基數(shù)A+T的含量都明顯高于堿基數(shù)G+C的含量，這與脊椎動(dòng)物線粒體DNA的L鏈堿基組成相一致[11-12]。A的含量明顯偏高，而G明顯偏低，與其他蛇類此片段堿基組成比例相當(dāng)[11，13-14]。

2.2 白唇竹葉青蛇擴(kuò)散路線

樣本產(chǎn)自海南島的8個(gè)市縣，按地理距離和地理相似性被劃分在5個(gè)地理單元內(nèi)地理分布(圖1)：昌江縣、東方市為西部單元(Western unit：Changjiang，Dongfang)，瓊海市為東部單元(Easternmost unit：Qionghai)；海口市為北部單元(Northern unit：Haikou)，樂東縣、三亞市、陵水縣為南部單元(Southern unit：Ledong，Sanya，Lingshui)，白沙縣為中部單元(Central unit：Baisha)。

基于12S、16S、Cytb、ND4的基因序列片段，采用DNASTAR軟件ClustalX程序排列DNA同源序列后，再通過人工詳細(xì)核查，將全部序列都視為無序特征。再將12S、16S、Cytb、ND4基因序列片段的seq文件分別加載到Megalin軟件后，根據(jù)其堿基突變位點(diǎn)數(shù)逐步比較排查，從而篩選出較為原始的單倍型序列，再根據(jù)所示遺傳距離圖，逐步有選擇性的局部序列比對(duì)和原始序列的比對(duì)，尋找突變位點(diǎn)、突變堿基數(shù)目，嘗試判斷突變方向，繪制出單倍型簡(jiǎn)約網(wǎng)絡(luò)圖[15]，如圖2(a)、3(a)、4(a)、5(a)所示。

2.2.1 基于12S基因，白唇竹葉青蛇擴(kuò)散路線

基于12S基因下，獲得分別產(chǎn)自海南島西部、東部、南部單元共12條白唇竹葉青蛇基因序列片段，經(jīng)詳細(xì)的剪輯比對(duì)排查發(fā)現(xiàn)393個(gè)堿基位點(diǎn)中有15個(gè)突變位點(diǎn)，共定義了9個(gè)單倍型。由圖2(a)-2(b)可知，基于12S基因下的單倍型網(wǎng)絡(luò)呈現(xiàn)出較清晰的星狀聚類圖。其中產(chǎn)自海南東方的C04和產(chǎn)自海南昌江的C05、C08的白唇竹葉青為同一單倍型，產(chǎn)自海南樂東的C01、C03及產(chǎn)自海南三亞的C10為同一單倍型，其他序列均為“孤序列”，即：?jiǎn)为?dú)一個(gè)序列對(duì)應(yīng)一個(gè)單倍型[15]。編號(hào)為C04、C06、C08的白唇竹葉青蛇相較原始且占主導(dǎo)優(yōu)勢(shì)，因此為中央單倍型，其他單倍型均以中央單倍型進(jìn)行擴(kuò)散進(jìn)化。其中編號(hào)C11、C02與中央單倍型相較遠(yuǎn)，均相差4個(gè)單倍型，因采集地點(diǎn)及個(gè)數(shù)較少，大多數(shù)的過渡類群的單倍體無法獲取，均以黑點(diǎn)表示未觀察到的中間單倍型。

根據(jù)圖2(a)，繪制出基于12S基因下，白唇竹葉青蛇擴(kuò)散路線圖2(b)。由圖可知，海南產(chǎn)白唇竹葉青蛇總體呈現(xiàn)出由西部向東部擴(kuò)散的趨勢(shì)，其起源中心分別是海南省的昌江縣、東方市這兩個(gè)西線市縣。根據(jù)單倍型可知，C07、C09可能由C04、C06、C08擴(kuò)散進(jìn)化得到，但C07、C09與C06、C08均產(chǎn)自昌江縣，地理位置相較為近，因此C07由C06、C08擴(kuò)散進(jìn)化得到的可能性大。同理，C04與C05均產(chǎn)自東方市，C05由C04擴(kuò)散進(jìn)化得到的可能性大。C05是C04擴(kuò)散進(jìn)化到C01、C03、C10的過渡類群，C07是C06、C09擴(kuò)散進(jìn)化到C01、C03、C10的過渡類群。

圖1 海南島白唇竹葉青蛇采集樣點(diǎn)Fig.1 Map of the sampling localities of Hainan Trimeresurus albolabris 注：黃色三角形編號(hào)1～10分別對(duì)應(yīng)B01～B10號(hào)樣本蛇；紅色三角形1～12分別對(duì)應(yīng)C01～C12號(hào)樣本蛇 Note：Yellow triangle No.1-10 corresponds to sample snake B01-B10，red triangle No.1-12 corresponds to sample snake C01-C12，respectively

圖2(a) 基于12S基因，白唇竹葉青蛇統(tǒng)計(jì)簡(jiǎn)約網(wǎng)絡(luò)圖Fig.2 (a) The statistics simple network diagram of Trimeresurusalbolabris in Hainan based on 12S gene 注：未觀察到的中間單倍型用黑點(diǎn)表示，每一步變異用或者較長或者較短的連接單倍型的線段表示，包括觀察到的和未觀察到的單倍型，圓圈的大小大致表示共享該單倍型的個(gè)體的數(shù)量[15] Note：Unobserved intermediate haplotypes are denoted by black dots，and each variation is denoted by lines of either longer or shorter connected haplotypes，including observed and unobserved haplotypes.The size of the circle roughly denotes the number of individuals sharing the haplotype[15]

圖2(b) 基于12S基因，白唇竹葉青蛇擴(kuò)散路線圖Fig.2 (b) The spread route map of Trimeresurus albolabris in Hainan based on the 12S gene 注：紅色小三角形1～12分別對(duì)應(yīng)C01～C12號(hào)樣本蛇 Note：Red small triangle No.1-12 corresponds to sample snake C01-C12

圖3(a) 基于16S基因，白唇竹葉青蛇統(tǒng)計(jì)簡(jiǎn)約網(wǎng)絡(luò)圖Fig.3 (a) The statistics simple network diagram of Trimeresurus albolabris in Hainan based on 16S gene 注：未觀察到的中間單倍型用黑點(diǎn)表示，每一步變異用或者較長或者較短的連接單倍型的線段表示，包括觀察到的和未觀察到的單倍型，圓圈的大小大致表示共享該單倍型的個(gè)體的數(shù)量[15] Note：Unobserved intermediate haplotypes are denoted by black dots，and each variation is denoted by lines of either longer or shorter connected haplotypes，including observed and unobserved haplotypes.The size of the circle roughly denotes the number of individuals sharing the haplotype[15]

圖3(b) 基于16S基因，白唇竹葉青蛇擴(kuò)散路線圖Fig.3 (b) The spread route map of Trimeresurus albolabris in Hainan based on the 16S gene 注：紅色小三角形1～12分別對(duì)應(yīng)C01～C12號(hào)樣本蛇 Note：Red small triangle No.1-12 corresponds to sample snake C01-C12

2.2.2 基于16S基因，白唇竹葉青蛇擴(kuò)散路線

基于16S基因，獲得分別產(chǎn)自海南島西部、東部、南部單元共9條白唇竹葉青蛇基因序列片段及產(chǎn)自中國香港的1條白唇竹葉青蛇基因序列片段。經(jīng)詳細(xì)的比對(duì)排查發(fā)現(xiàn)358個(gè)堿基位點(diǎn)中有16個(gè)突變位點(diǎn)，共定義了7個(gè)單倍型。由圖3(a)可知，基于16S基因下的單倍型網(wǎng)絡(luò)與基于12S基因下的單倍型網(wǎng)絡(luò)相似，均呈現(xiàn)出較清晰的星狀聚類圖。其中產(chǎn)自海南樂東的C01、東方的C04、海南昌江的C08的白唇竹葉青為同一單倍型，且相較其他白唇竹葉青為原始，占主導(dǎo)優(yōu)勢(shì)，因此為中央單倍型，其他單倍型均以中央單倍型進(jìn)行擴(kuò)散進(jìn)化。產(chǎn)自海南瓊海的C12與產(chǎn)自香港的HK為同一單倍型，其與中央單倍型相距遠(yuǎn)，相隔7個(gè)單倍類型。其他序列均為“孤序列”(singletons)，其中C11與中央單倍型差9個(gè)單倍型。因采集地點(diǎn)及個(gè)數(shù)較少，所有的過渡類群的單倍體無法獲取，均以黑點(diǎn)表示未觀察到的中間單倍型。

圖4(a) 基于Cyt b基因，白唇竹葉青蛇統(tǒng)計(jì)簡(jiǎn)約網(wǎng)絡(luò)圖Fig.4(a)The statistics simple network diagram of Trimeresurusalbolabris in Hainan based on Cyt b gene 注：未觀察到的中間單倍型用黑點(diǎn)表示，每一步變異用或者較長或者較短的連接單倍型的線段表示，包括觀察到的和未觀察到的單倍型，圓圈的大小大致表示共享該單倍型的個(gè)體的數(shù)量[15]. Note：Unobserved intermediate haplotypes are denoted by black dots,and each variation is denoted by lines of either longer or shorter connected haplotypes,including observed and unobserved haplotypes.The size of the circle roughly denotes the number of individuals sharing the haplotype[15]

圖4(b) 基于Cyt b基因，白唇竹葉青蛇擴(kuò)散路線圖Fig.4(b)The spread route map of Trimeresurus albolabris in Hainan based on the Cyt b gene 注：黃色三角形編號(hào)1～10分別對(duì)應(yīng)B01～B10號(hào)樣本蛇，紅色三角形1～12分別對(duì)應(yīng)C01～C12號(hào)樣本蛇。 Note：Yellow triangle No. 1-10 corresponds to sample snakes B01-B10 and red triangle No.1-12 corresponds to sample snakes C01-C12,respectively

圖5(a) 基于ND4基因，白唇竹葉青蛇統(tǒng)計(jì)簡(jiǎn)約網(wǎng)絡(luò)圖Fig.5(a)The statistics simple network diagram of Trimeresurusalbolabris in Hainan based on ND4 gene 注：未觀察到的中間單倍型用黑點(diǎn)表示，每一步變異用或者較長或者較短的連接單倍型的線段表示，包括觀察到的和未觀察到的單倍型，圓圈的大小大致表示共享該單倍型的個(gè)體的數(shù)量[13] Note：Unobserved intermediate haplotypes are denoted by black dots,and each variation is denoted by lines of either longer or shorter connected haplotypes,including observed and unobserved haplotypes.The size of the circle roughly denotes the number of individuals sharing the haplotype[13]

圖5(b) 基于ND4基因，白唇竹葉青蛇擴(kuò)散路線圖Fig.5(b)The spread route map of Trimeresurus albolabrisin Hainan based on the ND4 gene 注：黃色三角形編號(hào)1～10分別對(duì)應(yīng)B01～B10號(hào)樣本蛇，紅色三角形1～12分別對(duì)應(yīng)C01～C12號(hào)樣本蛇 Note：Yellow triangle No.1-10 corresponds to sample snakes B01-1B10, red triangle No.1-12 corresponds to sample snakes C01-C12,respectively

根據(jù)圖3(a)，繪制出基于16S基因下，白唇竹葉青蛇擴(kuò)散路線圖3(b)。由圖可知，海南產(chǎn)白唇竹葉青蛇總體呈現(xiàn)出由西部向東部擴(kuò)散的趨勢(shì)，其起源中心分別是海南省的昌江縣、東方市、樂東縣這3個(gè)西線市縣，這與基于12S基因下，所得的起源中心為昌江縣、東方市相似。根據(jù)單倍型可知，C09可能由C01、C04、C08擴(kuò)散進(jìn)化得到，但C09與C08均產(chǎn)自昌江縣，地理位置相較為近，因此C09由C08擴(kuò)散進(jìn)化得到的可能性大。同理，C02與C01均產(chǎn)自樂東縣，地理位置近，C02由C01擴(kuò)散進(jìn)化得到的可能性大。

2.2.3 基于Cytb基因，白唇竹葉青蛇擴(kuò)散路線

基于Cytb基因下，獲得分別產(chǎn)自海南島東部、南部、北部單元共13條白唇竹葉青蛇基因序列片段及產(chǎn)自中國香港的1條白唇竹葉青蛇基因序列片段。經(jīng)詳細(xì)的比對(duì)排查發(fā)現(xiàn)308個(gè)堿基位點(diǎn)中有10個(gè)突變位點(diǎn)，共定義了5個(gè)單倍型。由圖4(a)可知，基于Cytb基因下的單倍型網(wǎng)絡(luò)與基于12S、16S基因下的單倍型網(wǎng)絡(luò)相似，均呈現(xiàn)出較清晰的星狀聚類圖。其中產(chǎn)自海南樂東縣的C01、白沙縣的B10、?？谑械腂01、B02、B03及產(chǎn)自陵水縣的B04、B05、B06、B07、B08、B09為同一單倍型，且相較其他白唇竹葉青為原始，占主導(dǎo)優(yōu)勢(shì)，因此為中央單倍型，其他單倍型均以中央單倍型進(jìn)行擴(kuò)散進(jìn)化。產(chǎn)自海口市的B03與產(chǎn)自香港的HK為同一單倍型，其與中央單倍型相隔1個(gè)單倍類型。其他序列均為“孤序列”0，其中C10與中央單倍型相距最遠(yuǎn)差6個(gè)單倍型。因采集地點(diǎn)及個(gè)數(shù)較少，所有的過渡類群的單倍體無法獲取，均以黑點(diǎn)表示未觀察到的中間單倍型。

根據(jù)圖4(a)，繪制出基于Cytb基因下，白唇竹葉青蛇擴(kuò)散路線圖4(b)。由圖可知，海南產(chǎn)白唇竹葉青蛇總體呈現(xiàn)出由西部向東部擴(kuò)散的趨勢(shì)，其起源中心分別是海南省的?？谑?、白沙縣、樂東縣、陵水縣這4個(gè)市縣，起源中心相對(duì)分散，這可能與沒測(cè)序出產(chǎn)自海南島西部單元的蛇序列有關(guān)。根據(jù)單倍型親緣地理位置可知，B03由B01、B02擴(kuò)散進(jìn)化得到的可能性比較大，C11由B04、B05、B06、B07、B08、B09擴(kuò)散進(jìn)化得到的可能性大。

2.2.4 基于ND4基因，白唇竹葉青蛇擴(kuò)散路線

基于ND4基因下，獲得分別產(chǎn)自海南島西部、東部、南部、北部單元共21條白唇竹葉青蛇基因序列片段及產(chǎn)自中國香港的1條白唇竹葉青蛇基因序列片段。經(jīng)詳細(xì)的比對(duì)排查發(fā)現(xiàn)666個(gè)堿基位點(diǎn)中有13個(gè)突變位點(diǎn)，共定義了9個(gè)單倍型。由圖5(a)可知，基于ND4基因下的單倍型網(wǎng)絡(luò)與基于12S、16S、Cytb基因下的單倍型網(wǎng)絡(luò)相似，均呈現(xiàn)出較清晰的星狀聚類圖。其中產(chǎn)自海南東方市的C05、昌江縣的C08、C09為同一單倍型，且相較其他白唇竹葉青為原始，占主導(dǎo)優(yōu)勢(shì)，為中央單倍型，其他單倍型均以中央單倍型進(jìn)行擴(kuò)散進(jìn)化。B01、B02、B04、B05、B07、B08、B09、B10、C12、HK為同一單倍型，其與中央單倍型相差4個(gè)單倍類型且是中央單倍型擴(kuò)散進(jìn)化為B03的過渡類群?jiǎn)伪缎汀03、C10、C11為同一單倍型其與中央單倍型相隔2個(gè)單倍類型且是中央單倍型擴(kuò)散進(jìn)化為C02的過渡類群的單倍型。C06是中央單倍型擴(kuò)散進(jìn)化為C01的過渡類群的單倍型。因采集地點(diǎn)及個(gè)數(shù)較少，一部分的過渡類群的單倍體無法獲取，均以黑點(diǎn)表示未觀察到的中間單倍型。

根據(jù)圖5(a)，繪制出基于ND4基因下，白唇竹葉青蛇擴(kuò)散路線圖5(b)。由圖可知，海南產(chǎn)白唇竹葉青蛇總體呈現(xiàn)出由西部向東部擴(kuò)散的趨勢(shì)，其起源中心分別是海南省的昌江縣、東方市這兩個(gè)西線市縣。根據(jù)單倍型及地理親緣關(guān)系可知，C06、C07由C08、C09擴(kuò)散進(jìn)化得到的可能性較大，C04由C05擴(kuò)散進(jìn)化得到的可能性大，C02由C03擴(kuò)散進(jìn)化得到的可能性大，B03由B01、B02擴(kuò)散進(jìn)化得到的可能性大。

3 討論與結(jié)論

根據(jù)系統(tǒng)發(fā)育的分枝格局來推測(cè)種群的擴(kuò)散線路和方向，有很多成功的示例[16]，特別適合于對(duì)島嶼種群和擴(kuò)散能力較弱的兩棲爬行動(dòng)物種群的研究[17]。

由圖可知，基于12S、16S、Cytb、ND4基因序列片段的海南產(chǎn)白唇竹葉青蛇單倍型網(wǎng)絡(luò)圖均顯示了較清晰的星狀聚類圖。暗示了局部種群可能的地理起源和較距離的遷移事件，也證明了本研究的地理采樣有一個(gè)較好的覆蓋度，但因獲取的個(gè)體數(shù)目較少以致一部分的單倍型過渡類群無法獲取，均以黑點(diǎn)表示未觀察到的中間單倍型?；诓煌幕蛐蛄械玫降钠鹪粗行拇嬖诓町悾傮w擴(kuò)散進(jìn)化水平均呈現(xiàn)出由西部單元向東部單元擴(kuò)散的趨勢(shì)。但仍無法推斷出海南產(chǎn)白唇竹葉青蛇樣本與中國大陸、東南亞白唇竹葉青蛇的擴(kuò)散方向和較為準(zhǔn)確的起源地點(diǎn)。雖然海南島的氣候地理?xiàng)l件非常適合竹葉青蛇的生活習(xí)性，但是海南島與大陸存在著無可爭(zhēng)議的地理隔離。

DNA不僅是主要的遺傳物質(zhì)，同時(shí)也是生物進(jìn)化史的重要記錄者[18]。與核基因組相比，線粒體基因沿著一條直接而又簡(jiǎn)單的路線進(jìn)化，可由后代追溯到單一祖先；因此使得該分子適合用于親緣地理學(xué)研究和種群[15，19]。ND2和D-loop區(qū)域又因?yàn)檫M(jìn)化速率快最早成為該領(lǐng)域研究熱點(diǎn)[15]。Cytb基因在線粒體中，進(jìn)化速率中等[13]，為常被用作分子標(biāo)記的典型基因。但本實(shí)驗(yàn)所獲的Cytb基因單倍型較少，待獲取更多樣本后，進(jìn)一步分析。ND4基因進(jìn)化速率相對(duì)較快[6-7，10]，且PCR實(shí)驗(yàn)條件要求很低，因此經(jīng)常用于種下分子研究。本實(shí)驗(yàn)獲取的單倍型中，個(gè)體數(shù)最多的單倍型并未處于變化的中心區(qū)域，這與采樣點(diǎn)少、樣本隨機(jī)性較強(qiáng)有關(guān)。

就本實(shí)驗(yàn)中的12S、16S、Cytb、ND4基因序列片段而言，12S基因和16S基因的效果相對(duì)較好。其處于星圖的中心單倍型個(gè)體相對(duì)較多，拓?fù)浣Y(jié)構(gòu)清晰，且兩個(gè)基因的星圖一致性較高。由于本實(shí)驗(yàn)樣本量少且討論的地域面積不大，進(jìn)化速率相對(duì)較慢的12S和16S基因[8-9]，呈現(xiàn)出更清晰的結(jié)果。

竹葉青常生活于丘陵及山區(qū)，棲于溪邊灌叢、竹林、壁巖上或石縫間，喜陰雨天氣活動(dòng)[1-2]。而海南島地勢(shì)為四周低中間高，呈現(xiàn)出彎隆形態(tài)，其是由山地、丘陵、盆地、平原構(gòu)成階梯式環(huán)形層狀地貌[20]。再根據(jù)海南島歷年的降雨量，發(fā)現(xiàn)海南東部地區(qū)的降雨量遠(yuǎn)遠(yuǎn)多于海南島西部地區(qū)[21-22]。因此，推斷海南產(chǎn)白唇竹葉青蛇由西部地區(qū)向東部地區(qū)擴(kuò)散，尋求更適宜的棲息環(huán)境?；蛟S正是因海南島適宜的氣候及地理?xiàng)l件使得白唇竹葉青蛇能夠在此大量繁殖，迅速擴(kuò)大種群數(shù)量，豐富了白唇竹葉青蛇的基因庫，但因盆地、平原等居民密集區(qū)對(duì)白唇竹葉青蛇的大量捕食也可能是造成白唇竹葉青蛇數(shù)目銳減，而山地、丘陵等居民稀疏區(qū)多發(fā)現(xiàn)竹葉青蛇的原因之一。

因此作者得出初步結(jié)論，從親緣地理學(xué)上分析，海南產(chǎn)白唇竹葉青蛇整體呈現(xiàn)出西部地區(qū)往東部地區(qū)遷徙擴(kuò)散。