引用本文
郑梓豪, 吴珊珊, 魏勇, 钟代斌, 郑学礼. 基于线粒体COI基因分析广州市15个白纹伊蚊种群的遗传多样性[J]. 中国人兽共患病学报, ,37(11): 985-994
ZHENG Zi-hao, WU Shan-shan, WEI Yong, ZHONG Dai-bin, ZHENG Xue-li. Analysis of the genetic diversity of 15 Aedes albopictus populations in Guangzhou based on the mitochondrial COI gene [J]. Chinese Journal of Zoonoses, ,37(11): 985-994.  
Doi: 10.3969/j.issn.1002-2694.2021.00.141
Permissions
Copyright©2021, 中国人兽共患病学报
中国人兽共患病学报
基于线粒体COI基因分析广州市15个白纹伊蚊种群的遗传多样性
郑梓豪1, 吴珊珊1, 魏勇1, 钟代斌2, 郑学礼1
1. 南方医科大学公共卫生学院病原生物学系,广州 510515
2. 加州大学尔湾分校健康科学学院公共卫生项目组, 美国尔湾 92617
通讯作者：郑学礼,Email: zhengxueli2001@126.com; ORCID:0000-0003-0646-9398
收稿日期: 2021-04-21
资助项目: 国家自然科学基金项目(No. 31630011); 广州市科技计划项目(No. 201804020084)
摘要

目的 探讨广州市白纹伊蚊不同地理种群的遗传多样性、遗传分化和系统发育关系。方法 本实验于2020年9月至2020年11月期间,采集广州市11个行政区的共计15个白纹伊蚊种群。单只蚊虫提取基因组DNA,通过PCR法扩增COI基因序列并测序,获得的序列在GenBank上经过BLAST比对。BioEdit 7.2软件观察序列峰图。MEGA X软件对齐序列,分析碱基组成,构建系统发育树(NJ法)。DNAsp 6.12软件分析位点多态性,单倍型及其多样性,进行错配分布分析。Arlequin 3.5软件进行分子变异分析、进行中性检验。DAMBE 7.2软件分析序列的系统发育信号。 PopART 1.7软件构建单倍型网络图(Median joining法)。结果 广州市15个白纹伊蚊种群共获得642条序列,其长度为603 bp。A碱基加T碱基平均含量是67.5%,符合线粒体DNA的AT偏向性。单倍型分析共检出45种单倍型,其中单倍型1为优势单倍型。中性检验表明广州市部分种群不符合中性理论,但大部分种群都经历过种群扩张,而错配分布却表明仅有少部分种群经历过种群扩张。Mantel检验表明15个白纹伊蚊种群不存在地理隔离现象。种群遗传分化显示各种群间交流较频繁,且种群间的遗传分化很低,差异更多来自于个体间。结论 广州市的白纹伊蚊种群间基因交流较频繁,遗传分化小,遗传多样性偏低。这可能使不同地理种群的白纹伊蚊具有相似的媒介能力。

关键词: 白纹伊蚊; 线粒体细胞色素c氧化酶亚基Ⅰ; 种群遗传多样性
中图分类号:R384.1 文献标志码:A 文章编号:1002-2694(2021)11-0985-10
Analysis of the genetic diversity of 15 Aedes albopictus populations in Guangzhou based on the mitochondrial COI gene
ZHENG Zi-hao1, WU Shan-shan1, WEI Yong1, ZHONG Dai-bin2, ZHENG Xue-li1
1. Department of Pathogen Biology, School of Public Health, Southern Medical University, Guangzhou 510515, China
2. Program in Public Health, College of Health Sciences, University of California, Irvine 92617, USA
Corresponding author: Zheng Xue-li, Email: zhengxueli2001@126.com
Fund:Supported by the National Natural Science Foundation of China (No. 31630011); Guangzhou Science and Technology Project (No. 201804020084)
Abstract

This study aimed to investigate the genetic diversity, genetic differentiation, and phylogenetic relationship of the mitochondrial cytochrome c oxidase subunit I (COI) gene among 15 Aedes albopictus populations in Guangzhou. A total of 15 Ae. albopictus populations were collected from 11 districts of Guangzhou from September 2020 to November 2020. Genomic DNA was extracted from a single mosquito. The COI gene sequences were amplified by PCR methodology and sequenced afterwards. Sequences were blasted against GenBank. Chromatograms of sequences were observed with BioEdit 7.2 software. Sequence alignment, nucleotide composition analysis, and phylogenetic tree construction were conducted with MEGA X. Analysis of polymorphic sites, haplotype morphism and mismatch distribution was performed with DNAsp 6.12. A figure of the mismatch distribution was drawn with R 4.0.3 using the “ggplot2” and “cowplot” packages. The results of the molecular variance test and neutral hypothesis test were obtained with Arlequin 3.5. Substitution saturation was checked with DAMBE 7.2 and the haplotype network was drawn with PopART 1.7. In total, 642 sequences were obtained from the 15 Ae. albopictus populations in Guangzhou, and the length of the COI gene was 603 bp, with an average A+T content of 67.5%, which conforms to the mitochondrial AT bias. Forty-five haplotypes were detected and haplotype 1 was the dominant haplotype. Haplotype diversity (Hd) was 0.499 and nucleotide diversity (Π) was 0.00108. Neutrality tests showed that half of the Ae. albopictus populations (7/15) did not fit the neutrality theory and most (12/15) had experienced population expansion. However, the mismatch distribution showed that only a few of the Ae. albopictus populations had undergone population growth. The Mantel test showed that 15 Ae. albopictus populations were not geographically isolated. Population differentiation revealed a common occurrence of gene flow and greater variation within populations than between populations. In conclusion, gene flow has commonly occurred in Ae. albopictus populations in Guangzhou. Genetic differentiation and genetic diversities among populations were not significant. Therefore, Ae. albopictus populations in Guangzhou might have similar vector competence. The control of Ae. albopictus necessarily needs effort applied to all aspects of its biology.

Key words: Aedes albopictus; mitochondrial cytochrome C oxidase subunit Ⅰ; population genetic diversities

白纹伊蚊(Aedes albopictus)是登革病毒、寨卡病毒的主要传播媒介[1]。广州市的首个登革热病例报告于1978年, 之后广州市登革热疫情就不断出现, 在2014年出现了有记录以来最严重的登革热疫情的暴发[2]。而目前针对登革热流行的预防方法主要是对携带登革病毒的白纹伊蚊通过使用化学杀虫剂进行消杀[3, 4, 5]。细胞色素c氧化酶亚基Ⅰ (COI)基因在线粒体内, 其为母系遗传且遗传进化速率较核DNA快, 比核DNA更易受遗传漂变的影响, 利用线粒体DNA多态性分析不仅可以在一定层面上了解种群的遗传多样性, 还可以探讨形成特定遗传结构的历史原因和地理原因, 使其成为研究分子进化遗传常用分子标记物, 并能运用于建立物种鉴定的DNA条形码系统的研究, 目前已被广泛应用于昆虫的系统进化和种内遗传分化的分析[6, 7, 8, 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’ s D检验的统计值的P值小于0.05的种群有7个, 分别是:SMU、SS、TH、YX、HP、NS、ZC, 表明广州市部分白纹伊蚊种群不符合中性理论; Fu’ s FS检验的统计值的P值小于0.05的种群有12个, 分别是:SMU、NZHY、SS、DJ、GY、HX、TH、YX、PY、HP、NS、ZC, 表明了广州市的大部分种群(12/15)都在历史上发生过扩张。

表6 广州市白纹伊蚊种群中性检验 Tab.6 Neutral test of Ae. albopictus population in Guangzhou

如果错配分布(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之间相差两个突变步骤。

图4 基于Median joining network方法构建的单倍型网络图
注:单倍型间连线上的每条短线代表一个突变步骤。Fig.4 Haplotype network diagram based on the medium joining network method

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’ s FS对近期种群的扩张更为敏感[24]。当Tajima’ D检验为负值且有统计学意义时, 可认为不符合中性理论, 否则认为符合中性理论, 当Fu’ s FS检验为负值且有统计学意义时, 认为种群在近期经历过扩张, 否则认为种群的发展平稳[25]。

错配分布(Mismatch distribution)是关于样品中所有DNA序列之间成对差异数(Pairwise difference)的分布, 而来自高度隔离的亚种群的合并样本中的平均序列差异会增加。在长期处于种群平衡状态的种群中, 分布变得参差不齐, 而不断增长的种群则产生了平滑且具有峰值的错配分布, 高峰的位置反映了种群增长的时间[26]。错配分布中并没有太多的种群出现种群扩张的迹象。出现种群扩张迹象的部分种群的单峰也不明显。张瑞玲等[12]有研究显示:较高的单倍型多样性和较低的核苷酸多样性表明白纹伊蚊种群可能经历过“ 瓶颈效应” , 之后伴随着种群快速扩张的历史事件, 但由于扩张历史不长, 核苷酸变异积累不充分。广州市白纹伊蚊的单倍型多样性远大于核苷酸多样性, 这说明广州市的15个不同地点的白纹伊蚊种群可能经历过“ 瓶颈效应” 后发生了种群快速扩张, 但扩张历史不长, 核苷酸变异积累不充分。

种群分化系数(Fst), 又叫固定系数(Fixation index), 表示种群之间的遗传分化程度, 范围在0~1之间, 值越大表明种群的分化程度越高。一般认为Fst< 0.05表明种群出现轻度分化, 0.05< Fst< 0.15表明种群出现中度分化, 0.15< Fst< 0.25表明种群出现高度分化[27]。分子变异分析表明广州市各个行政区的白纹伊蚊种群遗传差异很小, 遗传差异主要是来自于个体之间。两两白纹伊蚊种群的种群分化系数表明松南阁种群、越秀种群、南沙种群和番禺种群和其余大多数种群存在中度分化, 其余则是轻度分化, 尚未发现广州市的白纹伊蚊种群间出现高度分化。基因流(Nm)表示基因在群体中的流动, Nm越大说明群体之间交流越频繁, Nm小于1时表明种群间基因交流极少, 存在很大分化, 当Nm的值大于1时可以抑制群体间的分化[28]。 而且各个白纹伊蚊种群间的基因交流较为频繁(Nm值都大于1), 所以各种群之间的分化程度只有低到中度分化。同时, Mantel检验也表明广州市的白纹伊种群不存在地理隔离现象, 进一步说明它们的各地的白纹伊蚊存在较为频繁的基因交流, 遗传差异的主要来自于个体间而不是群体间。

系统发育树中各单倍型的聚类也没有表现出明显的地域特征, 同样提示各个行政区的白纹伊蚊种群遗传差异小。白纹伊蚊虽然飞行距离有限, 但是其蚊卵能够耐受恶劣环境, 通过便利的交通工具间接地促进了白纹伊蚊的扩散[29]。同时也有研究表明白纹伊蚊是侵袭性很强的物种, 甚至在白纹伊蚊能取代埃及伊蚊而成为优势种群[30, 31]。同时由于广州市重要的经济地位, 使得广州市的交通运输和经济贸易运输十分发达, 同样也可能会促进白纹伊蚊的扩散。

白纹伊蚊可携带并传播多种虫媒病毒, 其作为一种侵袭性很强的物种, 对人类的健康造成了不可忽视的影响。本研究通过对广州市15个白纹伊蚊种群的COI基因进行分析, 发现了广州市的白纹伊蚊种群处于群体扩张阶段, 频繁的基因交流使得广州市不同行政区之间的白纹伊蚊种群间分化水平不高, 没有发现任意两两白纹伊蚊种群出现高度分化的情况, 大部分种群都属于低度分化状态, 各个种群间的遗传多样性不高。这些使得不同地理种群的白纹伊蚊具有相似的媒介能力[12]。

利益冲突: 无

引用本文格式:郑梓豪, 吴珊珊, 魏勇, 等.基于线粒体COI基因分析广州市15个白纹伊蚊种群的遗传多样性[J].中国人兽共患病学报, 2021, 37(11):985-994. DOI:10.3969/j.issn.1002-2694.2021.00.141

编辑:张智芳

参考文献
[1] Moyes CL, Vontas J, Martins AJ, et al. Contemporary status of insecticide resistance in the major Aedes vectors of arboviruses infecting humans[J]. PLoS Negl Trop Dis, 2017, 11(7): e5625. DOI: 10.1371/journal.pntd.0005625 [本文引用:1]
[2] 梁琳琳, 李晓宁, 罗雷, 等. 广州市2016—2017年白纹伊蚊密度及登革热传播风险分析[J]. 中华卫生杀虫药械, 2019, 25(2): 148-152. [本文引用:1]
[3] 彭丽兰, 解锐历, 言慧, 等. 2017年广州市白纹伊蚊对5种杀虫剂的抗药性调查[J]. 寄生虫与医学昆虫学报, 2018, 25(4): 206-211. DOI: 10.3969/j.issn.1005-0507.2018.04.003 [本文引用:1]
[4] 蔡松武, 段金花, 卢文成, 等. 广州市城区白纹伊蚊抗药性现状与预防对策探讨[J]. 华南预防医学, 2004, 30(4): 40-41. DOI: 10.3969/j.issn.1005-0507.2018.04.003 [本文引用:1]
[5] 李美青, 刘永权, 黄燕惠, 等. 广州市城乡不同季节白纹伊蚊密度、带病毒情况及抗药性监测[J]. 中国媒介生物学及控制杂志, 2013, 24(2): 108-111. [本文引用:1]
[6] 胡坤敏, 陈韶红, 艾琳, 等. 豫皖闽浙4省溪蟹并殖吸虫囊蚴核糖体ITS2和线粒体COI基因序列分析[J]. 中国寄生虫学与寄生虫病杂志, 2020, 38(1): 87-94. DOI: 10.12140/j.issn.1000-7423.2020.01.013 [本文引用:1]
[7] 赵健, 胡劭骥, 和秋菊, 等. 基于线粒体COⅠ和ND1基因分析三尾凤蝶遗传多样性[J]. 福建农林大学学报(自然科学版), 2017, 46(4): 387-391. [本文引用:1]
[8] 周慧, 李迪强, 张于光, 等. 藏羚羊mtDNA D-loop区遗传多样性研究[J]. 遗传, 2006, 28: 299-305. DOI: 10.3321/j.issn:0253-9772.2006.03.009 [本文引用:1]
[9] 王刚. 基于DNA条形码的我国主要蚊虫分子分类系统的建立[D]. 北京: 中国人民解放军军事医学科学院, 2011. [本文引用:1]
[10] 蔡燕丽, 罗雷, 郑学礼. 广州市白纹伊蚊细胞色素C氧化酶亚基Ⅰ基因遗传多态性分析[J]. 中国人兽共患病学报, 2012, 28(8): 781-785. DOI: 10.3969/j.issn.1002-2694.2012.08.004 [本文引用:1]
[11] 郭玉燕, 罗雷, 宋璋瑶, 等. 我国部分地区28种蚊虫的线粒体COⅠ基因序列分析[J]. 中国寄生虫学与寄生虫病杂志, 2017, 35(3): 280-287. [本文引用:1]
[12] 张瑞玲, 姚广琴, 潘晓倩, 等. 不同地理种群白纹伊蚊线粒体基因COI的遗传多样性分析[J]. 中国人兽共患病学报, 2017, 33(4): 316-320. DOI: 10.3969/j.issn.1002-2694.2017.04.005 [本文引用:3]
[13] 游丽斌, 李玉伟, 王金章, 等. 福建省白纹伊蚊mtDNA-COI基因的系统发育分析[J]. 中国病原生物学杂志, 2019, 14(5): 521-524. DOI: 10.13350/j.cjpb.190506 [本文引用:2]
[14] 陆宝麟, 吴厚永. 中国重要医学昆虫分类与鉴别[M]. 郑州: 河南科学技术出版社, 2003: 3299-3302. [本文引用:1]
[15] Folmer O, Black M, Hoeh W, et al. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates[J]. Mol Mar Biol Biotechnol, 1994, 3(5): 294-299. [本文引用:1]
[16] Kumar S, Stecher G, Li M, et al. MEGA X: molecular evolutionary genetics analysis across computing platforms[J]. Mol Biol Evol, 2018, 35(6): 1547-1549. DOI: 10.1093/molbev/msy096 [本文引用:1]
[17] Rozas J, Ferrer-Mata A, Sánchez-Delbarrio J C, et al. DnaSP 6: DNA sequence polymorphism analysis of large data sets[J]. Mol Biol Evol, 2017, 34(12): 3299-3302. DOI: 10.1093/molbev/msx248 [本文引用:1]
[18] Excoffier L, Lischer HE. Arlequin suite ver 3. 5: a new series of programs to perform population genetics analyses under Linux and Windows[J]. Mol Ecol Resour, 2010, 10(3): 564-567. DOI: 10.1111/j.1755-0998.2010.02847.x [本文引用:1]
[19] Xia X. DAMBE7: new and improved tools for data analysis in molecular biology and evolution[J]. Mol Biol Evol, 2018, 35(6): 1550-1552. DOI: 10.1093/molbev/msy073 [本文引用:1]
[20] Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees[J]. Mol Biol Evol, 1987, 4(4): 406-425. DOI: 10.1093/oxfordjournals.molbev.a040454 [本文引用:1]
[21] Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences[J]. J Mol Evol, 1980, 16(2): 111-120. DOI: 10.1007/BF01731581 [本文引用:1]
[22] Leigh J, Bryant D. popart: full-feature software for haplotype network construction[J]. Methods Ecol Evol, 2015, 6: 1110-1116. DOI: 10.1111/2041-210X.12410 [本文引用:1]
[23] Band elt HJ, Forster P, Röhl A. Median-joining networks for inferring intraspecific phylogenies[J]. Mol Biol Evol, 1999, 16(1): 37-48. DOI: 10.1093/oxfordjournals.molbev.a026036 [本文引用:1]
[24] 苏思博, 黄松, 董丙君. 东北三省棕黑锦蛇(Elaphe schrenckii)与赤峰锦蛇(Elaphe anomala)基于Cyt b基因的序列多态性及遗传多样性[J]. 长春师范大学学报, 2019, 38(04): 76-86. DOI: 10.3969/j.issn.1008-178X.2019.02.016 [本文引用:1]
[25] 常雪莲, 钟代斌, 李小聪, 等. 基于mtDNA-COI基因序列分析我国中华按蚊种群的遗传结构[J]. 南方医科大学学报, 2015, 35(02): 234-238. DOI: 10.3969/j.issn.1673-4254.2015.02.15 [本文引用:1]
[26] Harpending HC. Signature of ancient population growth in a low-resolution mitochondrial DNA mismatch distribution[J]. Hum Biol, 1994, 66(4): 591-600. [本文引用:1]
[27] 郭颂, 凌锋, 王金娜, 等. 浙江省不同地理株白纹伊蚊mtDNA-COⅠ基因多态性研究[J]. 中国人兽共患病学报, 2016, 32(02): 133-136. DOI: 10.3969/j.issn.1002-2694.2016.02.007 [本文引用:1]
[28] 郑宇婷, 姜进勇, 周红宁, 等. 利用ssr标记浅析云南登革热重点地区埃及伊蚊种群遗传特征[J]. 中国人兽共患病学报, 2021, 37(02): 176-182. DOI: 10.3969/j.issn.1002-2694.2021.00.014 [本文引用:1]
[29] Bonizzoni M, Gasperi G, Chen X, et al. The invasive mosquito species Aedes albopictus: current knowledge and future perspectives[J]. Trends Parasitol, 2013, 29(9): 460-468. DOI: 10.1016/j.pt.2013.07.003 [本文引用:1]
[30] 杜建伟, 潘先海. 中国登革热流行概况与流行特征[J]. 中华流行病学杂志, 2010(12): 1429-1433. DOI: 10.3760/cma.j.issn.0254-6450.2010.12.025 [本文引用:1]
[31] 宁丹, 孙九峰, 彭志强, 等. 广东省登革热流行概况与流行特征[J]. 华南预防医学, 2017, 43(04): 368-372. DOI: 10.13217/j.scjpm.2017.0368 [本文引用:1]