Fitness costs in the presence and absence of insecticide use explains abundance of two common Aedes aegypti kdr resistance alleles found in the Americas.

Aedes aegypti is the vector of viruses such as chikungunya, dengue, yellow fever and Zika that have a critical impact on human health. Control of adult mosquitoes is widely done using pyrethroids, but resistance has reduced the effectiveness of this class of insecticides. Resistance to pyrethroids in mosquitoes is commonly due to mutations in the voltage-gated sodium channel (Vgsc) gene (these mutations are known as knockdown resistance, kdr). In the Americas and the Caribbean, the most common kdr alleles are 410L+1016I+1534C and 1534C. In this study, we conducted a population cage experiment to evaluate changes in the allele and genotype frequencies of the 410L+1016I+1534C allele by crossing two congenic strains; one carrying the 410L+1016I+1534C and another with the 1534C allele. Changes in allele frequencies were measured over 10 generations in the absence of insecticide exposure. We also applied one cycle of selection with deltamethrin at F9 to evaluate the changes in allele and genotype frequencies. Our findings indicate that fitness costs were higher with the 410L+1016I+1534C allele, relative to the 1534C allele, in the absence of deltamethrin exposure, but that the 410L+1016I+1534C allele provides a stronger advantage when exposed to deltamethrin relative to the 1534C allele. Changes in genotype frequencies were not in Hardy-Weinberg equilibrium and could not be explained by drift. Our results suggest the diametrically opposed fitness costs in the presence and absence of insecticides is a reason for the variations in frequencies between the 410L+1016I+1534C and 1534C alleles in field populations.

Frequencies and distribution of kdr and Ace alleles that cause insecticide resistance in house flies in the United States.

House flies (Musca domestica L) are nuisances and vectors of pathogens between and among humans and livestock. Population suppression has been accomplished for decades with pyrethroids and acetylcholinesterase (AChE) inhibitors, but recurrent selection has led to increased frequency of alleles conferring resistance to those two classes of active ingredients (Geden et al., 2021). A common mechanism of resistance to both classes involves an altered target site (mutations in Voltage gated sodium channel (Vgsc) for pyrethroids or in Ace for AChE inhibitors). As part of ongoing efforts to understand the origin, spread and evolution of insecticide resistance alleles in house fly populations, we sampled flies in 11 different US states, sequenced, and then estimated frequencies of the Vgsc and Ace alleles. There was substantial variation in frequencies of the four common knockdown resistance alleles (kdr (L1014F), kdr-his (L1014H), super-kdr (M918T + L10414F) and 1B (T929I + L1014F) across the sampled states. The kdr allele was found in all 11 states and was the most common allele in four of them. The super-kdr allele was detected in only six collections, with the highest frequencies found in the north, northeast and central United States. The kdr-his allele was the most common allele in PA, NC, TN and TX. In addition, a novel super-kdr-like mutation in mutually exclusive exon 17a was found. The overall frequencies of the different Ace alleles, which we name based on the amino acid present at the mutation sites (V260L, A316S, G342A/V and F407Y), varied considerably between states. Five Ace alleles were identified: VAGF, VAVY, VAGY, VAAY and VSAY. Generally, the VSAY allele was the most common in the populations sampled. The susceptible allele (VAGF) was found in all populations, ranging in frequency from 3% (KS) to 41% (GA). Comparisons of these resistance allele frequencies with those previously found suggests a dynamic interaction between the different alleles, in terms of levels of resistance they confer and likely fitness costs they impose in the absence of insecticides.

A Globally Distributed Insecticide Resistance Allele Confers a Fitness Cost in the Absence of Insecticide in Aedes aegypti (Diptera: Culicidae), the Yellow Fever Mosquito.

The cosmopolitan mosquito Aedes aegypti is a vector of harmful arboviruses. Pyrethroid insecticides are used to reduce adult populations and prevent the spread of disease. Pyrethroids target the insect voltage-gated sodium channel (VGSC). Collectively, mutations in Vgsc that confer resistance are referred to as knock-down resistance or kdr. There are numerous kdr mutations found in A. aegypti Vgsc, and there is co-occurrence of some mutations. Full-length cDNA sequences have identified nine known kdr (e.g., 1534C) alleles. The 1534C allele is among the most common kdr alleles, but allele frequencies between populations vary considerably. We used the 1534C:RK strain, which has the 1534C (kdr) allele in the genetic background of the insecticide susceptible Rockefeller (ROCK) strain, and conducted population cage experiments to assess the potential intrinsic fitness cost of the 1534C allele relative to the susceptible allele (F1534) in the ROCK strain. Individuals were genotyped across generations using allele specific PCR. A fitness cost of the 1534C allele was detected across seven generations of mosquitos reared in the absence of insecticide selection pressure. The decrease in allele frequency was not due to drift. Comparison of our results to previous studies suggests that the magnitude of the fitness cost of kdr alleles in the absence of insecticide is disconnected from the level of resistance they confer, and that the fitness costs of different kdr alleles can be variable.

A chromosome-level assembly of the widely used Rockefeller strain of Aedes aegypti, the yellow fever mosquito.

Aedes aegypti is the vector of important human diseases, and genomic resources are crucial in facilitating the study of A. aegypti and its ecosystem interactions. Several laboratory-acclimated strains of this mosquito have been established, but the most used strain in toxicology studies is "Rockefeller," which was originally collected and established in Cuba 130 years ago. A full-length genome assembly of another reference strain, "Liverpool," was published in 2018 and is the reference genome for the species (AaegL5). However, genetic studies with the Rockefeller strain are complicated by the availability of only the Liverpool strain as the reference genome. Differences between Liverpool and Rockefeller have been known for decades, particularly in the expression of genes relevant to mosquito behavior and vector control (e.g. olfactory). These differences indicate that AaegL5 is likely not fully representative of the Rockefeller genome, presenting potential impediments to research. Here, we present a chromosomal-level assembly and annotation of the Rockefeller genome and a comparative characterization vs the Liverpool genome. Our results set the stage for a pan-genomic approach to understanding evolution and diversity within this important disease vector.

Out from under the wing: reconceptualizing the insect wing gene regulatory network as a versatile, general module for body-wall lobes in arthropods.

Body plan evolution often occurs through the differentiation of serially homologous body parts, particularly in the evolution of arthropod body plans. Recently, homeotic transformations resulting from experimental manipulation of gene expression, along with comparative data on the expression and function of genes in the wing regulatory network, have provided a new perspective on an old question in insect evolution: how did the insect wing evolve? We investigated the metamorphic roles of a suite of 10 wing- and body-wall-related genes in a hemimetabolous insect, Oncopeltus fasciatus. Our results indicate that genes involved in wing development in O. fasciatus play similar roles in the development of adult body-wall flattened cuticular evaginations. We found extensive functional similarity between the development of wings and other bilayered evaginations of the body wall. Overall, our results support the existence of a versatile development module for building bilayered cuticular epithelial structures that pre-dates the evolutionary origin of wings. We explore the consequences of reconceptualizing the canonical wing-patterning network as a bilayered body-wall patterning network, including consequences for long-standing debates about wing homology, the origin of wings and the origin of novel bilayered body-wall structures. We conclude by presenting three testable predictions that result from this reconceptualization.

Something old, something new, something borrowed, something red: the origin of ecologically relevant novelties in Hemiptera.

Comparative transcriptomics, applied in an evolutionary context, has transformed the possibilities for studying phenotypic evolution in non-model taxa. We review recent discoveries about the development of novel, ecologically relevant phenotypes in hemipteran insects. These discoveries highlight the diverse genomic substrates of novelty: 'something old', when novelty results from changes in the regulation of existing genes or gene duplication; 'something new', wherein lineage-restricted genes contribute to the evolution of new phenotypes; and 'something borrowed', showcasing contributions of horizontal gene transfer to the evolution of novelty, including carotenoid synthesis (resulting in 'something red'). These findings show the power and flexibility of comparative transcriptomic approaches for expanding beyond the 'toolkit' model for the evolution of development. We conclude by raising questions about the relationship between new genes and new traits and outlining a research framework for answering them in Hemiptera.

Co-option of wing-patterning genes underlies the evolution of the treehopper helmet.

Understanding the origin of novelty is a key question in evolutionary developmental biology. In arthropods, the body wall has served as a repeated source of morphological novelty. In treehoppers, an ancestrally flat part of the dorsal body wall (the pronotum) was transformed into a three-dimensional structure (the helmet), which was subsequently moulded by natural selection into diverse shapes. Here, we test three hypotheses for the developmental origin of the helmet by comparing body-region transcriptomes in a treehopper and a leafhopper that retains more ancestral morphology. In leafhoppers, pronotal gene expression is most similar to that of its serial homologue, the mesonotum. By contrast, in treehoppers, helmet gene expression is most similar to that of wings, supporting the wing-patterning network co-option hypothesis for the origin of the helmet. These results suggest that serial homologues may diverge evolutionarily through replacement of, rather than tinkering with, their ancestrally shared patterning network.

EnTAP: Bringing faster and smarter functional annotation to non-model eukaryotic transcriptomes.

EnTAP (Eukaryotic Non-Model Transcriptome Annotation Pipeline) was designed to improve the accuracy, speed, and flexibility of functional gene annotation for de novo assembled transcriptomes in non-model eukaryotes. This software package addresses the fragmentation and related assembly issues that result in inflated transcript estimates and poor annotation rates of protein-coding transcripts. Following filters applied through assessment of true expression and frame selection, open-source tools are leveraged to functionally annotate the reduced set of translated proteins. Downstream features include fast similarity search across five repositories, protein domain assignment, orthologous gene family assessment, and Gene Ontology (GO) term assignment. The final annotation integrates across multiple databases and selects an optimal assignment from a combination of weighted metrics describing similarity search score, taxonomic relationship, and informativeness. Researchers have the option to include additional filters to identify and remove contaminants, identify associated pathways, and prepare the transcripts for enrichment analysis. This fully featured pipeline is easy to install, configure, and runs significantly faster than comparable annotation packages. EnTAP is optimized to generate extensive functional information for the gene space of organisms with limited or poorly characterized genomic resources.