High-throughput genotyping of a full voltage-gated sodium channel gene via genomic DNA using target capture sequencing and analytical pipeline MoNaS to discover novel insecticide resistance mutations

Insects’ voltage-gated sodium channel (VGSC) is the primary target site of pyrethroid insecticides. Various amino acid substitutions in the VGSC protein are known to confer insecticide resistance and are selected under insecticide pressure. In the genome, the VGSC gene consists of more than 30 exons sparsely distributed across a large genomic region, which often exceeds 100 kbp. Due to this complex genomic structure of gene VGSC , it is usually challenging to genotype full coding nucleotide sequences (CDSs) of VGSC from individual genomic DNA (gDNA). In this study, we designed biotinylated oligonucleotide probes via annotated CDSs of VGSC of Asian tiger mosquito, Aedes albopictus . The probe set effectively concentrated (>80,000-fold) all targeted regions of gene VGSC from pooled barcoded Illumina libraries each constructed from individual A. albopictus gDNAs. The probe set also captured all homologous VGSC CDSs except tiny exons from the gDNA of other Culicinae mosquitos, A. aegypti and Culex pipiens complex, with comparable efficiency by virtue of the high conservation of VGSC at the nucleotide level. Furthermore, we developed an automated bioinformatic analysis pipeline to genotype VGSC after capture sequencing—MoNaS ( Mo squito Na + channel mutation S earch)—which conducts mapping of reads, variant calling, and variant annotation for nonsynonymous mutations. The proposed method and our bioinformatic tool should facilitate the discovery of novel amino acid variants conferring insecticide resistance on VGSC and population genetics studies on resistance alleles with respect to the origin, selection, and migration of both clinically and agriculturally important insect pests.