基于線粒體COI基因分析廣州市15個白紋伊蚊種群的遺傳多樣性

鄭梓豪，吳珊珊，魏 勇，鐘代斌，鄭學禮

白紋伊蚊(Aedesalbopictus)是登革病毒、寨卡病毒的主要傳播媒介[1]。廣州市的首個登革熱病例報告于1978年，之后廣州市登革熱疫情就不斷出現，在2014年出現了有記錄以來最嚴重的登革熱疫情的暴發[2]。而目前針對登革熱流行的預防方法主要是對攜帶登革病毒的白紋伊蚊通過使用化學殺蟲劑進行消殺[3-5]。細胞色素c氧化酶亞基Ⅰ(COI)基因在線粒體內，其為母系遺傳且遺傳進化速率較核DNA快，比核DNA更易受遺傳漂變的影響，利用線粒體DNA多態性分析不僅可以在一定層面上了解種群的遺傳多樣性，還可以探討形成特定遺傳結構的歷史原因和地理原因，使其成為研究分子進化遺傳常用分子標記物，并能運用于建立物種鑒定的DNA條形碼系統的研究，目前已被廣泛應用于昆蟲的系統進化和種內遺傳分化的分析[6-9]。蔡燕麗等利用COI基因及SSCP技術分析了廣州市13個區域的白紋伊蚊遺傳多態性[10]。郭玉燕等利用COI基因分析了我國部分地區的白紋伊蚊遺傳多態性[11]。張瑞玲等報告了不同地理種群白紋伊蚊線粒體基因COI的遺傳多樣性分析[12]。游麗斌等報告福建省白紋伊蚊 mtDNA-COI基因的系統發育分析[13]。我們在上述文獻報道的基礎上，開展本項研究，分析白紋伊蚊的系統進化和種內遺傳分化，是這項研究的新穎之處。

本文通過采用更大的樣本量來分析廣州市不同白紋伊蚊種群的線粒體的細胞色素c氧化酶亞基Ⅰ基因(COI)片段的遺傳結構，從白紋伊蚊的COI基因多樣性，為昆蟲的系統進化和種內遺傳分化”的理論的提升方面，提供科學參考資料。

1 方 法

1.1 白紋伊蚊的采樣 本實驗只采集白紋伊蚊樣本，采集時間為2020年9月至2020年11月，采樣點選取是基于廣州市的行政區，行政街道的劃分，白紋伊蚊的孳生習性，棲息習性見圖1。成蚊采用人誘法，幼蟲采用幼蟲吸管法，并將幼蟲帶回實驗室飼養至成蚊。形態學鑒定方法參考陸寶麟和吳永厚編著的《中國重要醫學昆蟲分類與鑒別》[14]。經鑒定為白紋伊蚊的成蚊，在-20 ℃凍死后浸泡在無水乙醇中-80 ℃保存。

圖1 廣州市白紋伊蚊種群采樣地點分布圖Fig.1 Distribution of sampling sites of Ae. albopictus population in Guangzhou

1.2 基因組DNA的提取 采用DNA提取試劑盒(Omega，Insect DNA kit)提取單只成蚊的基因組DNA，-80 ℃備用。

1.3 白紋伊蚊線粒體COI基因片段的擴增和測序 擴增COI基因片段的引物參考Folmer[15]等人使用的通用引物，正向引物(5′-GGTCAACAAATCATAAAGATATTGG-3′)，反向引物(5′-TGATTTTTTGGTCACCCTG AAGTTTA-3′)。PCR反應擴增體系為25 μL，其中Premix Taq(Takara公司)12.5 μL，正向引物(10 μmol/L)1 μL，反向引物(10 μmol/L)1 μL，模板DNA 2 μL,滅菌ddH2O 8.5 μL。PCR反應條件：94 ℃ 1 min；94 ℃ 40 s，45 ℃ 40 s，72 ℃ 1 min，5個循環；94 ℃ 40 s，53 ℃ 40 s，72 ℃ 1 min，30個循環。72 ℃ 7 min。4 ℃保存。用0.5%TBE電泳緩沖液，取PCR產物3～5 μL， 1.2%瓊脂糖凝膠(含有0.5 μg/mL溴化乙錠染料)中電泳，150 V電壓30 min后，紫外燈下觀察是否有目的條帶。出現目標條帶的樣本送北京擎科新業生物技術有限公司用反向引物測序。

1.4 線粒體COI基因片段分析 將序列在NCBI上的GenBank經過BLAST比對，與白紋伊蚊線粒體COI基因(KX266719.1)匹配度在99%以上的確認為白紋伊蚊(分子生物學鑒定)。用BioEdit 7.2.軟件觀察序列峰圖，用MEGA X[16]基于MUSCLE算法比對序列，使目標序列等長并截去兩端測序不準確部分，并統計序列的堿基組成。線粒體COI基因片段的位點多態性和DNA多態性分析，錯配分布分析使用DNAsp 6.12[17]軟件進行操作。分子變異分析，各種群間的固定系數(Fst)和基因流(Nm)的計算，Mantel檢驗，中性檢驗使用Arlequin 3.5軟件[18]。系統發育信號檢測使用DAMBE 7.2[19]軟件進行操作。如果序列未飽和，說明適合系統發育分析，則使用MEGA X軟件構建系統發育樹，統計方法采用鄰接法(Neighbor-Joining)[20]，Bootstrap值設置為1 000，替換模型使用Kimura 2-parameter雙參數法[21]，空位處理使用成對刪除(Pairwise deletion)，其余選項默認。構建系統發育樹時從GenBank中獲得了一株埃及伊蚊的序列(外群)和多個不同地理株的白紋伊蚊序列，序列號分別是：NC_035159.1(埃及伊蚊)、KF406541.1(巴基斯坦)、AB738203.1(日本)、MF148287.1(馬來西亞1)、MF148279.1(馬來西亞2)、JQ235749.1(云南)、KP843400.1(泰國)、GQ143719.1(澳大利亞)。單倍型網絡圖使用PopART[22]繪制，適用于分析種內數據的Median joining network[23]方法。

2 結 果

2.1 蚊蟲樣本采集結果 經過了形態學鑒定和分子生物學鑒定，并刪去部分峰圖質量不佳的樣本序列后，最終確定了642只白紋伊蚊用于本次實驗。采樣地點見圖1和表1。白云區包括南方醫科大學(SMU)和松南閣(SNG)兩個白紋伊蚊種群，海珠區包括南洲花苑(NZHY)和素社(SS)兩個白紋伊蚊種群，荔灣區包括東漖(DJ)、廣州醫科大學附屬第三醫院(GY)、鶴翔(HX)3個白紋伊蚊種群，天河區(天河山莊，TH)、越秀區(盤福社區，YX)、從化區(從化區第七中學，CH)、番禺區(小平村，PY)、花都區(東鏡村，HD)、黃埔區(新田村，HP)、南沙區(潭洲公園，NS)、增城區(南崗村，ZC)各只有1個白紋伊蚊種群。

表1 廣州市15個白紋伊蚊種群采樣點詳細信息Tab.1 Information about sampling sites among 15 Ae. albopictus populations in Guangzhou

2.2 線粒體COI基因堿基組成分析 廣州市15個種群一共642條核苷酸序列，長度603 bp。A、T、C、G四種堿基組成統計結果顯示(表2)：第1位和第2位堿基A加T的比例比G加C的比例稍高，第3位堿基的A和T的含量遠比G和C的高，A加T的總含量比例更是高達89.6%。4種堿基的平均含量分別為A(28.7%)、T(38.8%)、C(16.9%)、G(15.6%)，A堿基和T堿基的平均含量是67.5%，G堿基和C堿基的平均含量是32.5%，表明存在AT偏向性，符合線粒體DNA的特性。根據AT偏斜值公式：(A-T)/(A+T)和GC偏斜值公式：(G-C)/(G+C)，得出廣州市白紋伊蚊種群的AT偏斜值為：-0.150，GC偏斜值為-0.04。對所有序列基于K2P計算種內距離，發現其范圍在0～0.67%，平均為：0.12%。不同采樣點之間的白紋伊蚊種群的種內距離范圍在0.05%(HX、PY)～0.16%(NZHY、ZC)。

表2 廣州市15個白紋伊蚊種群線粒體COI基因核苷酸組成Tab.2 Nucleotide composition of mitochondrial COI gene in 15 Ae. albopictus populations in Guangzhou

2.3 線粒體COI基因序列特征 位點多態性分析結果顯示：603個位點中保守位點有564個，變異位點有39個(共有40種突變類型)，單態位點有15個，簡約信息位點有24個，無插入和缺失，37處出現轉換，3處出現顛換，轉換/顛換比為7.9，轉換遠比顛換多。

DNA多樣性指數(表3)分別為：廣州市整體的白紋伊蚊種群的單倍型多樣性(Haplotype diversity，Hd)是0.499，核苷酸多樣性(Nucleotide diversity，Pi)是0.001 08，平均核苷酸差異數(Average number of nucleotide difference，k)是0.653。各白紋伊蚊種群單倍型分析如下：單倍型多樣性在0.236(HX)～0.725(DJ)，核苷酸多樣性在0.000 40(HX)～0.001 54(DJ)。平均核苷酸差異數在0.243(HX)～0.931(DJ)。DNA多態性分析結果(表4)顯示：642條序列共檢出單倍型有45個。其中特有單倍型有33個，共享單倍型有12個，共享單倍型占所有單倍型的26.7%。種群中包含10個個體以上的單倍型有5種，分別是Hap-1(450)、Hap-3(20)、Hap-4(57)、Hap-10(16)、Hap-15(13)，其中Hap-1為所有種群共享,Hap-3、Hap-4為部分種群共享。

表3 廣州市15個白紋伊蚊種群線粒體COI基因遺傳多樣性Tab.3 Genetic diversity of mitochondrial COI gene of 15 Ae. albopictus populations in Guangzhou

表4 廣州市15個白紋伊蚊種群線粒體COI基因單倍型分布及頻率Tab.4 Distribution and frequency of mitochondrial COI gene haplotypes in 15 Ae. albopictus populations in Guangzhou

2.4 15個地理種群的遺傳結構 分子變異分析(Analysis of molecular variance, AMOVA)結果表明：15個種群間COⅠ基因的遺傳差異占4.09%，而種群內部的差異占95.91%。固定系數(Fixation Index，Fst)的值為0.040 9，P<0.01，認為廣州市不同種群之間的白紋伊蚊的遺傳差異主要來自種群個體內部，而各個種群之間的白紋伊蚊種群則沒有表現出大的差異性。兩兩種群間的固定系數Fst(表5)在-0.02～0.13之間，SS和HP種群間的Fst值最小，PY/HD和SNG種群間的Fst值最大。除去部分Nm值為負值外，所有地點兩兩配對的Nm值都大于1，而且還有部分地點的兩兩配對的Nm值大于10，說明廣州市11個行政區的的白紋伊蚊都有不同程度的基因交流。

表5 廣州市白紋伊蚊種群之間的固定系數(Fst)統計值(左下部分)和基因流(Nm)統計值(右上部分)Tab.5 Statistic of fixation(Fst, bottom left) and Gene flow(Nm, top right) of Ae. Albopictus population in Guangzhou

Mantel檢驗的回歸系數為0.000 127，相關系數為0.072 0，P值為0.354。不能認為地理距離和遺傳距離(Fst)這兩者存在相關性，即廣州市15個白紋伊蚊種群不存在地理隔離現象。

中性檢驗(表6)顯示所有種群的兩種檢驗的統計值都為負值。Tajima’sD檢驗的統計值的P值

表6 廣州市白紋伊蚊種群中性檢驗Tab.6 Neutral test of Ae. albopictus population in Guangzhou

小于0.05的種群有7個，分別是：SMU、SS、TH、YX、HP、NS、ZC，表明廣州市部分白紋伊蚊種群不符合中性理論；Fu’sFS檢驗的統計值的P值小于0.05的種群有12個，分別是：SMU、NZHY、SS、DJ、GY、HX、TH、YX、PY、HP、NS、ZC，表明了廣州市的大部分種群(12/15)都在歷史上發生過擴張。

如果錯配分布(Mismatch distribution)的實際觀察結果與期望值的某個假定模型(恒定或收縮)相吻合，則說明群體狀態與期望的假定模型一致，反之，則說明背離期望值的假定模型，群體處于擴張狀態(曲線為單峰)。在圖2中可見DJ、ZC種群的分布為單峰分布，與期望的分布不一致，認為這些種群經歷過擴張，SMU、NZHY、GY種群的分布在圖形中有小凸起，結合中性檢驗的結果來看，可能是曾經經歷過種群擴張，其余種群的分布與期望分布基本一致，認為這些種群的發展較平穩。

圖2 廣州市白紋伊蚊線粒體COI基因的錯配分布Fig.2 Mismatch distribution of mitochondrial COI gene of Ae. albopictus in Guangzhou

2.5 廣州市白紋伊蚊單倍型的系統發育關系 DAMBE的序列飽和度分析會給出Iss和Iss.C兩個統計值，當Iss的統計值比Iss.C的統計值小，且具有統計學意義的時候，可以認為序列沒有達到飽和。結果為：Iss的值為0.001和0.002，Iss.C的值在0.716～0.804之間，P值小于0.01，表明序列沒有達到飽和，可以用來建立系統發育樹。各單倍型的聚類關系結果如圖3所示，埃及伊蚊與所有的單倍型分開，單倍型34、單倍型3獨自成為一支。單倍型16和單倍型17聚成為一個小支，接著其余單倍型則分成3個大支，其他地理種群的白紋伊蚊則分布在這3個大分支。巴基斯坦種群和單倍型21聚成一個小支，日本種群和單倍型15、單倍型5、單倍型22聚成一個小分支，云南種群和單倍型23聚成一個小支。

圖3 廣州市15個白紋伊蚊種群線粒體COI基因各單倍型間的聚類關系(NJ法)Fig.3 Clustering relationship among mitochondrial COI gene haplotypes of 15 Ae. albopictus populations in Guangzhou (NJ method)

在單倍型網絡(圖4)可見單倍型1(被全部15個種群所共享)為優勢單倍型，基本所有的單倍型都以單倍型1為中心，說明單倍型1可能是祖先序列，其他單倍型則可能從單倍型1演化而來。而單倍型4，單倍型15也成為一個小中心。各單倍型之間只相差一個突變步驟，而單倍型1和單倍型37之間相差兩個突變步驟。

3 討 論

本次實驗采樣點包括廣州市11個行政區，共15個白紋伊蚊種群。廣州市15個白紋伊蚊種群的單倍型多樣性在0.256～0.759之間，核苷酸多樣性在0.000 50～0.001 59之間，平均核苷酸差異數在0.289～1.003 6之間，其中核苷酸多樣性和廣州市2011年的白紋伊蚊樣品(0.000 6)相符合，略高于2017年福建省的白紋伊蚊樣品(0.001 26)，低于深圳市的白紋伊蚊樣品(0.016 7)[13]。

Mantel檢驗的相關系數為0.033 4，但是沒有統計學意義。原因可能有：①基因流(Nm)的計算結果表明廣州市的白紋伊蚊種群存在較頻繁的基因交流；②分子變異分析表明種群間的遺傳距離很小。所以導致了遺傳距離和地理距離之間沒有相關性。

當兩種中性檢驗的統計值接近0時，說明種群處于穩定階段，而統計值為負值并且有統計學意義的時候，認為種群在歷史上有過擴張跡象，Tajima’D更傾向于揭示古老種群發生擴張的歷史，而Fu’sFS對近期種群的擴張更為敏感[24]。當Tajima’D檢驗為負值且有統計學意義時，可認為不符合中性理論，否則認為符合中性理論，當Fu’sFS檢驗為負值且有統計學意義時，認為種群在近期經歷過擴張，否則認為種群的發展平穩[25]。

錯配分布(Mismatch distribution)是關于樣品中所有DNA序列之間成對差異數(Pairwise difference)的分布，而來自高度隔離的亞種群的合并樣本中的平均序列差異會增加。在長期處于種群平衡狀態的種群中，分布變得參差不齊，而不斷增長的種群則產生了平滑且具有峰值的錯配分布,高峰的位置反映了種群增長的時間[26]。錯配分布中并沒有太多的種群出現種群擴張的跡象。出現種群擴張跡象的部分種群的單峰也不明顯。張瑞玲等[12]有研究顯示：較高的單倍型多樣性和較低的核苷酸多樣性表明白紋伊蚊種群可能經歷過“瓶頸效應”，之后伴隨著種群快速擴張的歷史事件，但由于擴張歷史不長，核苷酸變異積累不充分。廣州市白紋伊蚊的單倍型多樣性遠大于核苷酸多樣性，這說明廣州市的15個不同地點的白紋伊蚊種群可能經歷過“瓶頸效應”后發生了種群快速擴張，但擴張歷史不長，核苷酸變異積累不充分。

種群分化系數(Fst)，又叫固定系數(Fixation index)，表示種群之間的遺傳分化程度，范圍在0～1之間，值越大表明種群的分化程度越高。一般認為Fst<0.05表明種群出現輕度分化，0.05

系統發育樹中各單倍型的聚類也沒有表現出明顯的地域特征，同樣提示各個行政區的白紋伊蚊種群遺傳差異小。白紋伊蚊雖然飛行距離有限，但是其蚊卵能夠耐受惡劣環境，通過便利的交通工具間

注：單倍型間連線上的每條短線代表一個突變步驟。圖4 基于Median joining network方法構建的單倍型網絡圖Fig.4 Haplotype network diagram based on the medium joining network method

接地促進了白紋伊蚊的擴散[29]。同時也有研究表明白紋伊蚊是侵襲性很強的物種，甚至在白紋伊蚊能取代埃及伊蚊而成為優勢種群[30-31]。同時由于廣州市重要的經濟地位，使得廣州市的交通運輸和經濟貿易運輸十分發達，同樣也可能會促進白紋伊蚊的擴散。

白紋伊蚊可攜帶并傳播多種蟲媒病毒，其作為一種侵襲性很強的物種，對人類的健康造成了不可忽視的影響。本研究通過對廣州市15個白紋伊蚊種群的COI基因進行分析，發現了廣州市的白紋伊蚊種群處于群體擴張階段，頻繁的基因交流使得廣州市不同行政區之間的白紋伊蚊種群間分化水平不高，沒有發現任意兩兩白紋伊蚊種群出現高度分化的情況，大部分種群都屬于低度分化狀態，各個種群間的遺傳多樣性不高。這些使得不同地理種群的白紋伊蚊具有相似的媒介能力[12]。

利益沖突：無

引用本文格式：鄭梓豪，吳珊珊，魏勇，等.基于線粒體COI基因分析廣州市15個白紋伊蚊種群的遺傳多樣性[J].中國人獸共患病學報，2021，37(11)：985-994. DOI:10.3969/j.issn.1002-2694.2021.00.141