Abstract
Insecticide resistance in the bed bug Cimex lectularius is poorly understood due to the lack of genome sequences for resistant strains. In Japan, we identified a resistant strain of C. lectularius that exhibits a higher pyrethroid resistance ratio compared to many previously discovered strains. We sequenced the genomes of the pyrethroid-resistant and susceptible strains using long-read sequencing, resulting in the construction of highly contiguous genomes (N50 of the resistant strain: 2.1 Mb and N50 of the susceptible strain: 1.5 Mb). Gene prediction was performed by BRAKER3, and the functional annotation was performed by the Fanflow4insects workflow. Next, we compared their amino acid sequences to identify gene mutations, identifying 729 mutated transcripts that were specific to the resistant strain. Among them, those defined previously as resistance genes were included. Additionally, enrichment analysis implicated DNA damage response, cell cycle regulation, insulin metabolism, and lysosomes in the development of pyrethroid resistance. Genome editing of these genes can provide insights into the evolution and mechanisms of insecticide resistance. This study expanded the target genes to monitor allele distribution and frequency changes, which will likely contribute to the assessment of resistance levels. These findings highlight the potential of genome-wide approaches to understand insecticide resistance in bed bugs.