中华按蚊代谢抗性相关基因多态性分析

目的测定中华按蚊现场样本的宏基因组,分析代谢解毒酶相关基因的多态性。方法在山东济宁曹县和北湖区用人帐诱捕采集按蚊,基于形态和分子特征鉴定确认蚊种,进行宏基因组Illumina测序。从NCBI数据库下载属于CYP450s、GSTs和Coe基因家族的中华按蚊8个代谢解毒酶基因序列。应用DNASTAR将宏基因组测序结果与代谢解毒酶相关基因序列进行对比,在SeqMan Pro软件上分别进行读取和SNP分析,根据开放阅读框确定SNP的位置,并判别其为同义突变SNP或非同义突变SNP。结果分子鉴定的山东济宁曹县和北湖区中华按蚊样本分别为31只和12只,分为2组进行宏基因组测序。Illumina测序每组样本产生的reads数近200 Mb,平均长度101 bp,GC含量平均为41%(曹县)和42%(北湖区),mapping到8个代谢解毒酶基因序列上的reads数为391~2810条,每个位点的平均覆盖度的范围为23.72~144.93。本研究中华按蚊8个代谢抗性基因共有135个SNP位点,其中109个导致同义突变,占80.74%,26个非同义突变SNP位点(19.26%);每Kb的SNPs位点数为8.73~29.94;非同义突变SNPs在密码子的位置主要是第1位,占50.00%(13/26),第2位和第3位分别占26.92%(7/26)和23.08%(6/26)。结论Illumina测序获得中华按蚊宏基因组序列中代谢抗性相关基因具有相当程度的多态性,可为进一步研究中华按蚊的抗性分子机制提供重要数据。