Improved reference genome of Aedes aegypti informs arbovirus vector control

Long non-coding RNAs as promising targets for controlling disease vector mosquitoes

Hematophagous female mosquitoes are important vectors of numerous devastating human diseases, posing a major public health threat. Effective prevention and control of mosquito-borne diseases rely considerably on progress in understanding the molecular mechanisms of various life activities, and accordingly, the molecules that regulate the various life activities of mosquitoes are potential targets for implementing future vector control strategies. Many long non-coding RNAs (lncRNAs) have been identified in mosquitoes and significant progress has been made in determining their functions. Here, we present a comprehensive overview of the research advances on mosquito lncRNAs, including their molecular identification, function, and interaction with other non-coding RNAs, as well as their synergistic regulatory roles in mosquito life activities. We also highlight the potential roles of competitive endogenous RNAs in mosquito growth and development, as well as in insecticide resistance and virus-host interactions. Insights into the biological functions and mechanisms of lncRNAs in mosquito life activities, viral replication, pathogenesis, and transmission will contribute to the development of novel drugs and safe vaccines.

Challenges of Robust RNAi-Mediated Gene Silencing in Aedes Mosquitoes

In this study, we report the complexities and challenges associated with achieving robust RNA interference (RNAi)-mediated gene knockdown in the mosquitoes Aedes aegypti and Aedes albopictus, a pivotal approach for genetic analysis and vector control. Despite RNAi's potential for species-specific gene targeting, our independent efforts to establish oral delivery of RNAi for identifying genes critical for mosquito development and fitness encountered significant challenges, failing to reproduce previously reported potent RNAi effects. We independently evaluated a range of RNAi-inducing molecules (siRNAs, shRNAs, and dsRNAs) and administration methods (oral delivery, immersion, and microinjection) in three different laboratories. We also tested various mosquito strains and utilized microorganisms for RNA delivery. Our results reveal a pronounced inconsistency in RNAi efficacy, characterized by minimal effects on larval survival and gene expression levels in most instances despite strong published effects for the tested targets. One or multiple factors, including RNase activity in the gut, the cellular internalization and processing of RNA molecules, and the systemic dissemination of the RNAi signal, could be involved in this variability, all of which are barely understood in mosquitoes. The challenges identified in this study highlight the necessity for additional research into the underlying mechanisms of mosquito RNAi to develop more robust RNAi-based methodologies. Our findings emphasize the intricacies of RNAi application in mosquitoes, which present a substantial barrier to its utilization in genetic control strategies.

Odorant receptors for floral- and plant-derived volatiles in the yellow fever mosquito, Aedes aegypti (Diptera: Culicidae)

Adult mosquitoes require regular sugar meals, including nectar, to survive in natural habitats. Both males and females locate potential sugar sources using sensory proteins called odorant receptors (ORs) activated by plant volatiles to orient toward flowers or honeydew. The yellow fever mosquito, Aedes aegypti (Linnaeus, 1762), possesses a large gene family of ORs, many of which are likely to detect floral odors. In this study, we have uncovered ligand-receptor pairings for a suite of Aedes aegypti ORs using a panel of environmentally relevant, plant-derived volatile chemicals and a heterologous expression system. Our results support the hypothesis that these odors mediate sensory responses to floral odors in the mosquito's central nervous system, thereby influencing appetitive or aversive behaviors. Further, these ORs are well conserved in other mosquitoes, suggesting they function similarly in diverse species. This information can be used to assess mosquito foraging behavior and develop novel control strategies, especially those that incorporate mosquito bait-and-kill technologies.

Effects of Arboviral Infections on Transposable Element Transcript Levels in Aedes aegypti

Transposable elements are mobile repeated sequences found in all genomes. Transposable elements are controlled by RNA interference pathways in most organisms, and this control involves the PIWI-interacting RNA pathway and the small interfering RNA pathway, which is also known to be the first line of antiviral defense in invertebrates. Using Drosophila, we recently showed that viral infections result in the modulation of transposable element transcript levels through modulation of the small RNA repertoire. The Aedes aegypti mosquito is of particular interest because almost half of its genome is made of transposable elements, and it is described as a major vector of viruses (such as the dengue [DENV], Zika [ZIKV], and chikungunya [CHIKV] arboviruses). Moreover, Aedes mosquitoes are unique among insects in that the PIWI-interacting RNA pathway is also involved in the somatic antiviral response, in addition to the transposable element control and PIWI-interacting RNA pathway genes expanded in the mosquito genome. For these reasons, we studied the impacts of viral infections on transposable element transcript levels in A. aegypti samples. We retrieved public datasets corresponding to RNA-seq data obtained from viral infections by DENV, ZIKV, and CHIKV in various tissues. We found that transposable element transcripts are moderately modulated following viral infection and that the direction of the modulation varies greatly across tissues and viruses. These results highlight the need for an in-depth investigation of the tightly intertwined interactions between transposable elements and viruses.

Antibodies to Aedes aegypti D7L salivary proteins as a new serological tool to estimate human exposure to Aedes mosquitoes

Introduction

Aedes spp. are the most prolific mosquito vectors in the world. Found on every continent, they can effectively transmit various arboviruses, including the dengue virus which continues to cause outbreaks worldwide and is spreading into previously non-endemic areas. The lack of widely available dengue vaccines accentuates the importance of targeted vector control strategies to reduce the dengue burden. High-throughput tools to estimate human-mosquito contact and evaluate vector control interventions are lacking. We propose a novel serological tool that allows rapid screening of human cohorts for exposure to potentially infectious mosquitoes.

Methods

We tested 563 serum samples from a longitudinal pediatric cohort study previously conducted in Cambodia. Children enrolled in the study were dengue-naive at baseline and were followed biannually for dengue incidence for two years. We used Western blotting and enzyme-linked immunosorbent assays to identify immunogenic Aedes aegypti salivary proteins and measure total anti-Ae. aegypti IgG.

Results

We found a correlation (rs=0.86) between IgG responses against AeD7L1 and AeD7L2 recombinant proteins and those to whole salivary gland homogenate. We observed seasonal fluctuations of AeD7L1+2 IgG responses and no cross-reactivity with Culex quinquefasciatus and Anopheles dirus mosquitoes. The baseline median AeD7L1+2 IgG responses for young children were higher in those who developed asymptomatic versus symptomatic dengue.

Discussion

The IgG response against AeD7L1+2 recombinant proteins is a highly sensitive and Aedes specific marker of human exposure to Aedes bites that can facilitate standardization of future serosurveys and epidemiological studies by its ability to provide a robust estimation of human-mosquito contact in a high-throughput fashion.

Roles of a newly lethal cuticular structural protein, AaCPR100A, and its upstream interaction protein, G12-like, in Aedes aegypti

Mosquitoes form a vital group of vector insects, which can transmit various diseases and filarial worms. The cuticle is a critical structure that protects mosquitoes from adverse environmental conditions and penetration resistance. Thus, cuticle proteins can be used as potential targets for controlling the mosquito population. In the present study, we found that AaCPR100A is a structural protein in the soft cuticle, which has flexibility and elasticity allowing insects to move or fly freely, of Aedes aegypti. RNA interference (RNAi) of AaCPR100A caused high mortality in Aedes aegypti larvae and adults and significantly decreased the egg hatching rate. Transmission electron microscopy (TEM) analysis revealed that the larval microstructure had no recognizable endocuticle in AaCPR100A-deficient mosquitoes. A yeast two-hybrid assay was performed to screen proteins interacting with AaCPR100A. We verified that the G12-like protein had the strongest interaction with AaCPR100A using yeast two-hybrid and GST pull-down assays. Knockdown of G12-like transcription resulted in high mortality in Ae. aegypti larvae, but not in adults. Interestingly, RNAi of G12-like rescued the high mortality of adults caused by decreased AaCPR100A expression. Additionally, adults treated with G12-like dsRNA were found to be sensitive to low temperature, and their eggshell formation and hatching were decreased. Overall, our results demonstrated that G12-like may interacts with AaCPR100A, and both G12-like and AaCPR100A are involved in Ae. aegypti cuticle development and eggshell formation. AaCPR100A and G12-like can thus be considered newly potential targets for controlling the Ae. aegypti mosquito.

Structural basis of ligand specificity and channel activation in an insect gustatory receptor

Gustatory receptors (GRs) are critical for insect chemosensation and are potential targets for controlling pests and disease vectors, making their structural investigation a vital step toward such applications. We present structures of Bombyx mori Gr9 (BmGr9), a fructose-gated cation channel, in agonist-free and fructose-bound states. BmGr9 forms a tetramer similar to distantly related insect odorant receptors (ORs). Upon fructose binding, BmGr9's channel gate opens through helix S7b movements. In contrast to ORs, BmGr9's ligand-binding pocket, shaped by a kinked helix S4 and a shorter extracellular S3-S4 loop, is larger and solvent accessible in both agonist-free and fructose-bound states. Also, unlike ORs, fructose binding by BmGr9 involves helix S5 and a pocket lined with aromatic and polar residues. Structure-based sequence alignments reveal distinct patterns of ligand-binding pocket residue conservation in GR subfamilies associated with different ligand classes. These data provide insight into the molecular basis of GR ligand specificity and function.

A high heterozygosity genome assembly of Aedes albopictus enables the discovery of the association of PGANT3 with blood-feeding behavior

The Asian tiger mosquito, Aedes albopictus, is a global invasive species, notorious for its role in transmitting dangerous human arboviruses such as dengue and Chikungunya. Although hematophagous behavior is repulsive, it is an effective strategy for mosquitoes like Aedes albopictus to transmit viruses, posing a significant risk to human health. However, the fragmented nature of the Ae. albopictus genome assembly has been a significant challenge, hindering in-depth biological and genetic studies of this mosquito. In this research, we have harnessed a variety of technologies and implemented a novel strategy to create a significantly improved genome assembly for Ae. albopictus, designated as AealbF3. This assembly boasts a completeness rate of up to 98.1%, and the duplication rate has been minimized to 1.2%. Furthermore, the fragmented contigs or scaffolds of AealbF3 have been organized into three distinct chromosomes, an arrangement corroborated through syntenic plot analysis, which compared the genetic structure of Ae. albopictus with that of Ae. aegypti. Additionally, the study has revealed a phylogenetic relationship suggesting that the PGANT3 gene is implicated in the hematophagous behavior of Ae. albopictus. This involvement was preliminarily substantiated through RNA interference (RNAi) techniques and behavioral experiment. In summary, the AealbF3 genome assembly will facilitate new biological insights and intervention strategies for combating this formidable vector of disease. The innovative assembly process employed in this study could also serve as a valuable template for the assembly of genomes in other insects characterized by high levels of heterozygosity.

Quantifying Mosquito Host Preference

The most dangerous mosquito species for human health are those that blood feed preferentially and frequently on humans (anthropophilic mosquitoes). These include prolific disease vectors such as the African malaria mosquito Anopheles gambiae and yellow fever mosquito Aedes aegypti The chemosensory basis for anthropophilic behavior exhibited by these disease vectors, as well as the factors that drive interindividual differences in human attractiveness to mosquitoes, remain largely uncharacterized. Here, we concisely review established methods to quantify mosquito interspecific and intraspecific host preference in the laboratory, as well as semi-field and field environments. Experimental variables for investigator consideration during assays of mosquito host preference across these settings are highlighted.

The updated genome of the Hungarian population of Aedes koreicus

Vector-borne diseases pose a potential risk to human and animal welfare, and understanding their spread requires genomic resources. The mosquito Aedes koreicus is an emerging vector that has been introduced into Europe more than 15 years ago but only a low quality, fragmented genome was available. In this study, we carried out additional sequencing and assembled and characterized the genome of the species to provide a background for understanding its evolution and biology. The updated genome was 1.1 Gbp long and consisted of 6099 contigs with an N50 value of 329,610 bp and a BUSCO score of 84%. We identified 22,580 genes that could be functionally annotated and paid particular attention to the identification of potential insecticide resistance genes. The assessment of the orthology of the genes indicates a high turnover at the terminal branches of the species tree of mosquitoes with complete genomes, which could contribute to the adaptation and evolutionary success of the species. These results could form the basis for numerous downstream analyzes to develop targets for the control of mosquito populations.

Demographic history and the efficacy of selection in the globally invasive mosquito Aedes aegypti

Aedes aegypti is the main vector species of yellow fever, dengue, zika and chikungunya. The species is originally from Africa but has experienced a spectacular expansion in its geographic range to a large swath of the world, the demographic effects of which have remained largely understudied. In this report, we examine whole-genome sequences from 6 countries in Africa, North America, and South America to investigate the demographic history of the spread of Ae. aegypti into the Americas its impact on genomic diversity. In the Americas, we observe patterns of strong population structure consistent with relatively low (but probably non-zero) levels of gene flow but occasional long-range dispersal and/or recolonization events. We also find evidence that the colonization of the Americas has resulted in introduction bottlenecks. However, while each sampling location shows evidence of a past population contraction and subsequent recovery, our results suggest that the bottlenecks in America have led to a reduction in genetic diversity of only ~35% relative to African populations, and the American samples have retained high levels of genetic diversity (expected heterozygosity of ~0.02 at synonymous sites) and have experienced only a minor reduction in the efficacy of selection. These results evoke the image of an invasive species that has expanded its range with remarkable genetic resilience in the face of strong eradication pressure.

A genotyping array for the globally invasive vector mosquito, Aedes albopictus

BACKGROUND

Although whole-genome sequencing (WGS) is the preferred genotyping method for most genomic analyses, limitations are often experienced when studying genomes characterized by a high percentage of repetitive elements, high linkage, and recombination deserts. The Asian tiger mosquito (Aedes albopictus), for example, has a genome comprising up to 72% repetitive elements, and therefore we set out to develop a single-nucleotide polymorphism (SNP) chip to be more cost-effective. Aedes albopictus is an invasive species originating from Southeast Asia that has recently spread around the world and is a vector for many human diseases. Developing an accessible genotyping platform is essential in advancing biological control methods and understanding the population dynamics of this pest species, with significant implications for public health.

METHODS

We designed a SNP chip for Ae. albopictus (Aealbo chip) based on approximately 2.7 million SNPs identified using WGS data from 819 worldwide samples. We validated the chip using laboratory single-pair crosses, comparing technical replicates, and comparing genotypes of samples genotyped by WGS and the SNP chip. We then used the chip for a population genomic analysis of 237 samples from 28 sites in the native range to evaluate its usefulness in describing patterns of genomic variation and tracing the origins of invasions.

RESULTS

Probes on the Aealbo chip targeted 175,396 SNPs in coding and non-coding regions across all three chromosomes, with a density of 102 SNPs per 1 Mb window, and at least one SNP in each of the 17,461 protein-coding genes. Overall, 70% of the probes captured the genetic variation. Segregation analysis found that 98% of the SNPs followed expectations of single-copy Mendelian genes. Comparisons with WGS indicated that sites with genotype disagreements were mostly heterozygotes at loci with WGS read depth < 20, while there was near complete agreement with WGS read depths > 20, indicating that the chip more accurately detects heterozygotes than low-coverage WGS. Sample sizes did not affect the accuracy of the SNP chip genotype calls. Ancestry analyses identified four to five genetic clusters in the native range with various levels of admixture.

CONCLUSIONS

The Aealbo chip is highly accurate, is concordant with genotypes from WGS with high sequence coverage, and may be more accurate than low-coverage WGS.

Finding divergent sequences of homomorphic sex chromosomes via diploidized nanopore-based assembly from a single male

Although homomorphic sex chromosomes can have non-recombining regions with elevated sequence divergence between its complements, such divergence signals can be difficult to detect bioinformatically. If found in genomes of e.g. insect pests, these sequences could be targeted by the engineered genetic sexing and control systems. Here, we report an approach that can leverage long-read nanopore sequencing of a single XY male to identify divergent regions of homomorphic sex chromosomes. Long-read data are used for de novo genome assembly that is diploidized in a way that maximizes sex-specific differences between its haploid complements. We show that the correct assembly phasing is supported by the mapping of nanopore reads from the male's haploid Y-bearing sperm cells. The approach revealed a highly divergent region (HDR) near the centromere of the homomorphic sex chromosome of Aedes aegypti, the most important arboviral vector, for which there is a great interest in creating new genetic control tools. HDR is located ~5Mb downstream of the known male-determining locus on chromosome 1 and is significantly enriched for ovary-biased genes. While recombination in HDR ceased relatively recently (~1.4 MYA), HDR gametologs have divergent exons and introns of protein coding genes, and most lncRNA genes became X-specific. Megabases of previously invisible sex-linked sequences provide new putative targets for engineering the genetic systems to control this deadly mosquito. Broadly, our approach expands the toolbox for studying cryptic structure of sex chromosomes.

Insight into Hyalomma anatolicum biology by comparative genomics analyses

Hyalomma anatolicum is an obligatory blood-sucking ectoparasite and contributes to the transmission of Crimean-Congo haemorrhagic fever (CCHF) virus, Theileria spp. and Babesia spp. Progress in exploring the adaptive strategy of this ectoparasite and developing tools to fight it has been hindered by the lack of a complete genome. Herein, we assembled the genome using diverse sources of data from multiple sequencing platforms and annotated the 1.96 Gb genome of Hy. anatolicum. Comparative genome analyses and the predicted protein encoding genes reveal unique facets of this genome, including gene family expansion associated with blood feeding and digestion, multi-gene families involved in detoxification, a great number of neuropeptides and corresponding receptors regulating tick growth, development, and reproduction, and glutathione S-transferase genes playing roles in insecticide resistance and detoxification of multiple xenobiotic factors. This high quality reference genome provides fundamental data for obtaining insights into a variety of aspects of tick biology and developing novel strategies to fight notorious tick vectors of human and animal pathogens.

Extensive variation and strain-specificity in dengue virus susceptibility among African Aedes aegypti populations

African populations of the mosquito Aedes aegypti are usually considered less susceptible to infection by human-pathogenic flaviviruses than globally invasive populations found outside Africa. Although this contrast has been well documented for Zika virus (ZIKV), it is unclear to what extent it is true for dengue virus (DENV), the most prevalent flavivirus of humans. Addressing this question is complicated by substantial genetic diversity among DENV strains, most notably in the form of four genetic types (DENV1 to DENV4), that can lead to genetically specific interactions with mosquito populations. Here, we carried out a survey of DENV susceptibility using a panel of seven field-derived Ae. aegypti colonies from across the African range of the species and a colony from Guadeloupe, French West Indies as non-African reference. We found considerable variation in the ability of African Ae. aegypti populations to acquire and replicate a panel of six DENV strains spanning the four DENV types. Although African Ae. aegypti populations were generally less susceptible than the reference non-African population from Guadeloupe, in several instances some African populations were equally or more susceptible than the Guadeloupe population. Moreover, the relative level of susceptibility between African mosquito populations depended on the DENV strain, indicating genetically specific interactions. We conclude that unlike ZIKV susceptibility, there is no clear-cut dichotomy in DENV susceptibility between African and non-African Ae. aegypti. DENV susceptibility of African Ae. aegypti populations is highly heterogeneous and largely governed by the specific pairing of mosquito population and DENV strain.

An expanded neurogenetic toolkit to decode olfaction in the African malaria mosquito Anopheles gambiae

Anopheles gambiae uses its sense of smell to hunt humans. We report a two-step method yielding cell-type-specific driver lines for enhanced neuroanatomical and functional studies of its olfactory system. We first integrated a driver-responder-marker (DRM) system cassette consisting of a linked T2A-QF2 driver, QUAS-GFP responder, and a gut-specific transgenesis marker into four chemoreceptor genes (Ir25a, Ir76b, Gr22, and orco) using CRISPR-Cas9-mediated homology-directed repair. The DRM system facilitated rapid selection of in-frame integrations via screening for GFP+ olfactory sensory neurons (OSNs) in G1 larval progeny, even at genomic loci such as orco where we found the transgenesis marker was not visible. Next, we converted these DRM integrations into T2A-QF2 driver-marker lines by Cre-loxP excision of the GFP responder, making them suitable for binary use in transcuticular calcium imaging. These cell-type-specific driver lines tiling key OSN subsets will support systematic efforts to decode olfaction in this prolific malaria vector.

Niemann-Pick Type C2 Proteins in Aedes aegypti: Molecular Modelling and Prediction of Their Structure-Function Relationships

Aedes aegypti is a major vector that transmits arboviruses through the saliva injected into the host. Salivary proteins help in uninterrupted blood intake and enhance the transmission of pathogens. We studied Niemann-Pick Type C2 (NPC2) proteins, a superfamily of saliva proteins that play an important role in arbovirus infections. In vertebrates, a single conserved gene encodes for the NPC2 protein that functions in cholesterol trafficking. Arthropods, in contrast, have several genes that encode divergent NPC2 proteins. We compared the sequences of 20 A. aegypti NPC2 proteins to the cholesterol-binding residues of human and bovine, and fatty-acid-binding residues of ant NPC2 protein. We identified four mosquito NPC2 proteins as potential sterol-binding proteins. Two of these proteins (AAEL006854 and/or AAEL020314) may play a key role in ecdysteroid biosynthesis and moulting. We also identified one mosquito NPC2 protein as a potential fatty-acid-binding protein. Through molecular modelling, we predicted the structures of the potential sterol- and fatty-acid-binding proteins and compared them to the reference proteins.

Targeting sex determination to suppress mosquito populations

Each year, hundreds of millions of people are infected with arboviruses such as dengue, yellow fever, chikungunya, and Zika, which are all primarily spread by the notorious mosquito Aedes aegypti. Traditional control measures have proven insufficient, necessitating innovations. In response, here we generate a next-generation CRISPR-based precision-guided sterile insect technique (pgSIT) for Ae. aegypti that disrupts genes essential for sex determination and fertility, producing predominantly sterile males that can be deployed at any life stage. Using mathematical models and empirical testing, we demonstrate that released pgSIT males can effectively compete with, suppress, and eliminate caged mosquito populations. This versatile species-specific platform has the potential for field deployment to effectively control wild populations of disease vectors.

Neofunctionalization driven by positive selection led to the retention of the loqs2 gene encoding an Aedes specific dsRNA binding protein

BACKGROUND

Mosquito borne viruses, such as dengue, Zika, yellow fever and Chikungunya, cause millions of infections every year. These viruses are mostly transmitted by two urban-adapted mosquito species, Aedes aegypti and Aedes albopictus. Although mechanistic understanding remains largely unknown, Aedes mosquitoes may have unique adaptations that lower the impact of viral infection. Recently, we reported the identification of an Aedes specific double-stranded RNA binding protein (dsRBP), named Loqs2, that is involved in the control of infection by dengue and Zika viruses in mosquitoes. Preliminary analyses suggested that the loqs2 gene is a paralog of loquacious (loqs) and r2d2, two co-factors of the RNA interference (RNAi) pathway, a major antiviral mechanism in insects.

RESULTS

Here we analyzed the origin and evolution of loqs2. Our data suggest that loqs2 originated from two independent duplications of the first double-stranded RNA binding domain of loqs that occurred before the origin of the Aedes Stegomyia subgenus, around 31 million years ago. We show that the loqs2 gene is evolving under relaxed purifying selection at a faster pace than loqs, with evidence of neofunctionalization driven by positive selection. Accordingly, we observed that Loqs2 is localized mainly in the nucleus, different from R2D2 and both isoforms of Loqs that are cytoplasmic. In contrast to r2d2 and loqs, loqs2 expression is stage- and tissue-specific, restricted mostly to reproductive tissues in adult Ae. aegypti and Ae. albopictus. Transgenic mosquitoes engineered to express loqs2 ubiquitously undergo developmental arrest at larval stages that correlates with massive dysregulation of gene expression without major effects on microRNAs or other endogenous small RNAs, classically associated with RNA interference.

CONCLUSIONS

Our results uncover the peculiar origin and neofunctionalization of loqs2 driven by positive selection. This study shows an example of unique adaptations in Aedes mosquitoes that could ultimately help explain their effectiveness as virus vectors.

The chromosome-scale genome assembly for the West Nile vector Culex quinquefasciatus uncovers patterns of genome evolution in mosquitoes

BACKGROUND

Understanding genome organization and evolution is important for species involved in transmission of human diseases, such as mosquitoes. Anophelinae and Culicinae subfamilies of mosquitoes show striking differences in genome sizes, sex chromosome arrangements, behavior, and ability to transmit pathogens. However, the genomic basis of these differences is not fully understood.

METHODS

In this study, we used a combination of advanced genome technologies such as Oxford Nanopore Technology sequencing, Hi-C scaffolding, Bionano, and cytogenetic mapping to develop an improved chromosome-scale genome assembly for the West Nile vector Culex quinquefasciatus.

RESULTS

We then used this assembly to annotate odorant receptors, odorant binding proteins, and transposable elements. A genomic region containing male-specific sequences on chromosome 1 and a polymorphic inversion on chromosome 3 were identified in the Cx. quinquefasciatus genome. In addition, the genome of Cx. quinquefasciatus was compared with the genomes of other mosquitoes such as malaria vectors An. coluzzi and An. albimanus, and the vector of arboviruses Ae. aegypti. Our work confirms significant expansion of the two chemosensory gene families in Cx. quinquefasciatus, as well as a significant increase and relocation of the transposable elements in both Cx. quinquefasciatus and Ae. aegypti relative to the Anophelines. Phylogenetic analysis clarifies the divergence time between the mosquito species. Our study provides new insights into chromosomal evolution in mosquitoes and finds that the X chromosome of Anophelinae and the sex-determining chromosome 1 of Culicinae have a significantly higher rate of evolution than autosomes.

CONCLUSION

The improved Cx. quinquefasciatus genome assembly uncovered new details of mosquito genome evolution and has the potential to speed up the development of novel vector control strategies.

Odorant receptors for floral- and plant-derived volatiles in the yellow fever mosquito, Aedes aegypti (Diptera: Culicidae)

Adult mosquitoes require regular sugar meals, usually floral nectar, to survive and flourish in natural habitats. Both males and females locate potential sugar sources using sensory proteins called odorant receptors activated by plant volatiles that facilitate orientation toward flowers or honeydew. The Yellow Fever mosquito, Aedes aegypti (Linnaeus, 1762), possesses a large repertoire of odorant receptors, many of which are likely to support floral odor detection and nectar-seeking. In this study, we have employed a heterologous expression system and the two-electrode voltage clamping technique to identify environmentally relevant chemical compounds that activate specific odorant receptors. Importantly, we have uncovered ligand-receptor pairings for a suite of Aedes aegypti odorant receptors likely to mediate appetitive or aversive behavioral responses, thus shaping a critical aspect of the life history of a medically important mosquito. Moreover, the high degree of conservation of these receptors in other disease-transmitting species suggests common mechanisms of floral odor detection. This knowledge can be used to further investigate mosquito foraging behavior to either enhance existing, or develop novel, control strategies, especially those that incorporate mosquito bait-and-kill or attractive toxic sugar bait technologies.

Influence of dengue virus serotypes on the abundance of Aedes aegypti insect-specific viruses (ISVs)

A comprehensive understanding of the virome in mosquito vectors is crucial for assessing the potential transmission of viral agents, designing effective vector control strategies, and advancing our knowledge of insect-specific viruses (ISVs). In this study, we utilized Oxford Nanopore Technologies metagenomics to characterize the virome of Aedes aegypti mosquitoes collected in various regions of Colombia, a country hyperendemic for dengue virus (DENV). Analyses were conducted on groups of insects with previous natural DENV infection (DENV-1 and DENV-2 serotypes), as well as mosquito samples that tested negative for virus infection (DENV-negative). Our findings indicate that the Ae. aegypti virome exhibits a similar viral composition at the ISV family and species levels in both DENV-positive and DENV-negative samples across all study sites. However, differences were observed in the relative abundance of viral families such as Phenuiviridae, Partitiviridae, Flaviviridae, Rhabdoviridae, Picornaviridae, Bromoviridae, and Virgaviridae, depending on the serotype of DENV-1 and DENV-2. In addition, ISVs are frequently found in the core virome of Ae. aegypti, such as Phasi Charoen-like phasivirus (PCLV), which was the most prevalent and showed variable abundance in relation to the presence of specific DENV serotypes. Phylogenetic analyses of the L, M, and S segments of the PCLV genome are associated with sequences from different regions of the world but show close clustering with sequences from Brazil and Guadeloupe, indicating a shared evolutionary relationship. The profiling of the Ae. aegypti virome in Colombia presented here improves our understanding of viral diversity within mosquito vectors and provides information that opens the way to possible connections between ISVs and arboviruses. Future studies aimed at deepening our understanding of the mechanisms underlying the interactions between ISVs and DENV serotypes in Ae. aegypti could provide valuable information for the design of effective vector-borne viral disease control and prevention strategies.IMPORTANCEIn this study, we employed a metagenomic approach to characterize the virome of Aedes aegypti mosquitoes, with and without natural DENV infection, in several regions of Colombia. Our findings indicate that the mosquito virome is predominantly composed of insect-specific viruses (ISVs) and that infection with different DENV serotypes (DENV-1 and DENV-2) could lead to alterations in the relative abundance of viral families and species constituting the core virome in Aedes spp. The study also sheds light on the identification of the genome and evolutionary relationships of the Phasi Charoen-like phasivirus in Ae. aegypti in Colombia, a widespread ISV in areas with high DENV incidence.

Mitochondrial DNA variation in Aedes aegypti (Diptera: Culicidae) mosquitoes from Jeddah, Saudi Arabia

Wolbachia (Hertig 1936) (Rickettsiales: Ehrlichiaceae) has emerged as a valuable biocontrol tool in the fight against dengue by suppressing the transmission of the virus through mosquitoes. Monitoring the dynamics of Wolbachia is crucial for evaluating the effectiveness of release programs. Mitochondrial (mtDNA) markers serve as important tools for molecular tracking of infected mitochondrial backgrounds over time but require an understanding of the variation in release sites. In this study, we investigated the mitochondrial lineages of Aedes aegypti (Linnaeus 1762) in Jeddah, Saudi Arabia, which is a prospective release site for the "wAlbBQ" Wolbachia-infected strain of this mosquito species. We employed a combination of comprehensive mitogenomic analysis (including all protein-coding genes) and mtDNA marker analysis (cox1 and nad5) using data collected from Jeddah. We combined our mitogenome and mtDNA marker data with those from previous studies to place mitochondrial variation in Saudi Arabia into a broader global context. Our findings revealed the presence of 4 subclades that can be broadly categorized into 2 major mitochondrial lineages. Ae. aegypti mosquitoes from Jeddah belonged to both major lineages. Whilst mitogenomic data offered a higher resolution for distinguishing Jeddah mosquitoes from the wAlbBQ strain, the combination of cox1 and nad5 mtDNA markers alone proved to be sufficient. This study provides the first important characterization of Ae. aegypti mitochondrial lineages in Saudi Arabia and offers essential baseline information for planning future molecular monitoring efforts during the release of Wolbachia-infected mosquitoes.

On the Origin and Evolution of the Mosquito Male-determining Factor Nix

The mosquito family Culicidae is divided into 2 subfamilies named the Culicinae and Anophelinae. Nix, the dominant male-determining factor, has only been found in the culicines Aedes aegypti and Aedes albopictus, 2 important arboviral vectors that belong to the subgenus Stegomyia. Here we performed sex-specific whole-genome sequencing and RNAseq of divergent mosquito species and explored additional male-inclusive datasets to investigate the distribution of Nix. Except for the Culex genus, Nix homologs were found in all species surveyed from the Culicinae subfamily, including 12 additional species from 3 highly divergent tribes comprising 4 genera, suggesting Nix originated at least 133 to 165 million years ago (MYA). Heterologous expression of 1 of 3 divergent Nix open reading frames (ORFs) in Ae. aegypti resulted in partial masculinization of genetic females as evidenced by morphology and doublesex splicing. Phylogenetic analysis suggests Nix is related to femaleless (fle), a recently described intermediate sex-determining factor found exclusively in anopheline mosquitoes. Nix from all species has a conserved structure, including 3 RNA-recognition motifs (RRMs), as does fle. However, Nix has evolved at a much faster rate than fle. The RRM3 of both Nix and fle are distantly related to the single RRM of a widely distributed and conserved splicing factor transformer-2 (tra2). The RRM3-based phylogenetic analysis suggests this domain in Nix and fle may have evolved from tra2 or a tra2-related gene in a common ancestor of mosquitoes. Our results provide insights into the evolution of sex determination in mosquitoes and will inform broad applications of mosquito-control strategies based on manipulating sex ratios toward nonbiting males.

Evolution of the insect Hox gene cluster: Comparative analysis across 243 species

The Hox gene cluster is an iconic example of evolutionary conservation between divergent animal lineages, providing evidence for ancient similarities in the genetic control of embryonic development. However, there are differences between taxa in gene order, gene number and genomic organisation implying conservation is not absolute. There are also examples of radical functional change of Hox genes; for example, the ftz, zen and bcd genes in insects play roles in segmentation, extraembryonic membrane formation and body polarity, rather than specification of anteroposterior position. There have been detailed descriptions of Hox genes and Hox gene clusters in several insect species, including important model systems, but a large-scale overview has been lacking. Here we extend these studies using the publicly-available complete genome sequences of 243 insect species from 13 orders. We show that the insect Hox cluster is characterised by large intergenic distances, consistently extreme in Odonata, Orthoptera, Hemiptera and Trichoptera, and always larger between the 'posterior' Hox genes. We find duplications of ftz and zen in many species and multiple independent cluster breaks, although certain modules of neighbouring genes are rarely broken apart suggesting some organisational constraints. As more high-quality genomes are obtained, a challenge will be to relate structural genomic changes to phenotypic change across insect phylogeny.

A tale of caution: How endogenous viral elements affect virus discovery in transcriptomic data

Large-scale metagenomic and -transcriptomic studies have revolutionized our understanding of viral diversity and abundance. In contrast, endogenous viral elements (EVEs), remnants of viral sequences integrated into host genomes, have received limited attention in the context of virus discovery, especially in RNA-Seq data. EVEs resemble their original viruses, a challenge that makes distinguishing between active infections and integrated remnants difficult, affecting virus classification and biases downstream analyses. Here, we systematically assess the effects of EVEs on a prototypical virus discovery pipeline, evaluate their impact on data integrity and classification accuracy, and provide some recommendations for better practices. We examined EVEs and exogenous viral sequences linked to Orthomyxoviridae, a diverse family of negative-sense segmented RNA viruses, in 13 genomic and 538 transcriptomic datasets of Culicinae mosquitoes. Our analysis revealed a substantial number of viral sequences in transcriptomic datasets. However, a significant portion appeared not to be exogenous viruses but transcripts derived from EVEs. Distinguishing between transcribed EVEs and exogenous virus sequences was especially difficult in samples with low viral abundance. For example, three transcribed EVEs showed full-length segments, devoid of frameshift and nonsense mutations, exhibiting sufficient mean read depths that qualify them as exogenous virus hits. Mapping reads on a host genome containing EVEs before assembly somewhat alleviated the EVE burden, but it led to a drastic reduction of viral hits and reduced quality of assemblies, especially in regions of the viral genome relatively similar to EVEs. Our study highlights that our knowledge of the genetic diversity of viruses can be altered by the underestimated presence of EVEs in transcriptomic datasets, leading to false positives and altered or missing sequence information. Thus, recognizing and addressing the influence of EVEs in virus discovery pipelines will be key in enhancing our ability to capture the full spectrum of viral diversity.

Antibodies to Aedes aegypti D7L salivary proteins as a new serological tool to estimate human exposure to Aedes mosquitoes

Introduction

Aedes spp. are the most prolific mosquito vectors in the world. Found on every continent, they can effectively transmit various arboviruses, including the dengue virus which continues to cause outbreaks worldwide and is spreading into previously non-endemic areas. The lack of widely available dengue vaccines accentuates the importance of targeted vector control strategies to reduce the dengue burden. High-throughput sensitive tools to estimate human-mosquito contact and evaluate vector control interventions are lacking. We propose a novel serological tool that allows rapid screening of large human cohorts for exposure to potentially infectious mosquitoes and effective targeting of vector control.

Methods

We tested 563 serum samples from a longitudinal pediatric cohort study previously conducted in Cambodia. Children enrolled in the study were dengue-naïve at baseline and were followed biannually for dengue incidence for two years. We used Western blotting and enzyme-linked immunosorbent assays to identify the most immunogenic Aedes aegypti salivary proteins and measure total anti- Ae. Aegypti IgG.

Results

We found a strong correlation (r s =0.86) between the combined IgG responses against AeD7L1 and AeD7L2 recombinant proteins and those to whole salivary gland homogenate. We observed seasonal fluctuations of AeD7L1+2 IgG responses, corresponding to Aedes spp. abundance in the region, and no cross-reactivity with Culex quinquefasciatus and Anopheles dirus mosquitoes. The baseline median AeD7L1+2 IgG responses for young children were higher in those who developed asymptomatic dengue versus those who developed symptomatic dengue.

Conclusion

The IgG response against AeD7L1+2 recombinant proteins is a highly sensitive and Aedes specific marker of human exposure to Aedes bites that can facilitate standardization of future serosurveys and epidemiological studies by its ability to provide a robust estimation of human-mosquito contact in a high-throughput fashion.

Intrinsic factors driving mosquito vector competence and viral evolution: a review

Mosquitoes are responsible for the transmission of numerous viruses of global health significance. The term "vector competence" describes the intrinsic ability of an arthropod vector to transmit an infectious agent. Prior to transmission, the mosquito itself presents a complex and hostile environment through which a virus must transit to ensure propagation and transmission to the next host. Viruses imbibed in an infectious blood meal must pass in and out of the mosquito midgut, traffic through the body cavity or hemocoel, invade the salivary glands, and be expelled with the saliva when the vector takes a subsequent blood meal. Viruses encounter physical, cellular, microbial, and immunological barriers, which are influenced by the genetic background of the mosquito vector as well as environmental conditions. Collectively, these factors place significant selective pressure on the virus that impact its evolution and transmission. Here, we provide an overview of the current state of the field in understanding the mosquito-specific factors that underpin vector competence and how each of these mechanisms may influence virus evolution.

Structure of an insect gustatory receptor

Gustatory Receptors (GRs) are critical for insect chemosensation and are potential targets for controlling pests and disease vectors. However, GR structures have not been experimentally determined. We present structures of Bombyx mori Gr9 (BmGr9), a fructose-gated cation channel, in agonist-free and fructose-bound states. BmGr9 forms a tetramer similar to distantly related insect Olfactory Receptors (ORs). Upon fructose binding, BmGr9's ion channel gate opens through helix S7b movements. In contrast to ORs, BmGR9's ligand-binding pocket, shaped by a kinked helix S4 and a shorter extracellular S3-S4 loop, is larger and solvent accessible in both agonist-free and fructose-bound states. Also unlike ORs, fructose binding by BmGr9 involves helix S5 and a binding pocket lined with aromatic and polar residues. Structure-based sequence alignments reveal distinct patterns of ligand-binding pocket residue conservation in GR subfamilies associated with distinct ligand classes. These data provide insight into the molecular basis of GR ligand specificity and function.

Extensive variation and strain-specificity in dengue virus susceptibility among African Aedes aegypti populations

African populations of the mosquito Aedes aegypti are usually considered less susceptible to infection by human-pathogenic flaviviruses than globally invasive populations found outside Africa. Although this contrast has been well documented for Zika virus (ZIKV), it is unclear to what extent it is true for dengue virus (DENV), the most prevalent flavivirus of humans. Addressing this question is complicated by substantial genetic diversity among DENV strains, most notably in the form of four genetic types (DENV1 to DENV4), that can lead to genetically specific interactions with mosquito populations. Here, we carried out a continent-wide survey of DENV susceptibility using a panel of field-derived Ae. aegypti colonies from across the African range of the species and a colony from Guadeloupe, French West Indies as non-African reference. We found considerable variation in the ability of African Ae. aegypti populations to acquire and replicate a panel of six DENV strains spanning the four DENV types. Although African Ae. aegypti populations were generally less susceptible than the reference non-African population from Guadeloupe, in several instances some African populations were equally or more susceptible than the Guadeloupe population. Moreover, the relative level of susceptibility between African mosquito populations depended on the DENV strain, indicating genetically specific interactions. We conclude that unlike ZIKV susceptibility, there is no clear-cut dichotomy in DENV susceptibility between African and non-African Ae. aegypti. DENV susceptibility of African Ae. aegypti populations is highly heterogeneous and largely governed by the specific pairing of mosquito population and DENV strain.

Sequencing and Analysis of the Mitochondrial Genome of Aedes aegypti (Diptera: Culicidae) from the Brazilian Amazon Region

Aedes aegypti is a mosquito native to the African continent, which is now widespread in the tropical and subtropical regions of the world. In many regions, it represents a major challenge to public health, given its role in the cycle of transmission of important arboviruses, such as Dengue, Zika, and Chikungunya. Considering the epidemiological importance of Ae. aegypti, the present study sequenced the partial mitochondrial genome of a sample collected in the municipality of Balsas, in the Brazilian state of Maranhão, followed by High Throughput Sequencing and phylogenetic analyses. The mitochondrial sequence obtained here was 15,863 bp long, and contained 37 functional subunits (thirteen PCGs, twenty-two tRNAs and two rRNAs) in addition to a partial final portion rich in A+T. The data obtained here contribute to the enrichment of our knowledge of the taxonomy and evolutionary biology of this prominent disease vector. These findings represent an important advancement in the understanding of the characteristics of the populations of northeastern Brazil and provide valuable insights into the taxonomy and evolutionary biology of this prominent disease vector.

In silico identification, characterization, and expression analysis of RNA recognition motif (RRM) containing RNA-binding proteins in Aedes aegypti

RNA-binding proteins (RBPs) are the proteins that bind RNAs and regulate their functioning. RBPs in mosquitoes are gaining attention due to their ability to bind flaviviruses and regulate their replication and transmission. Despite their relevance, RBPs in mosquitoes are not explored much. In this study, we screened the whole genome of Aedes aegypti, the primary vector of several pathogenic viruses, and identified the proteins containing RNA recognition motif (RRM), the most abundant protein domain in eukaryotes. Using several in silico strategies, a total of 135 RRM-containing RBPs were identified in Ae. aegypti. The proteins were characterized based on their available annotations and the sequence similarity with Drosophila melanogaster. Ae. aegypti RRM-containing RBPs included serine/arginine-rich (SR) proteins, polyadenylate-binding proteins (PABP), heteronuclear ribonucleoproteins (hnRNP), small nuclear ribonucleoproteins (snRNP), splicing factors, eukaryotic initiation factors, transformers, and nucleolysins. Phylogenetic analysis revealed that the proteins and the domain organization are conserved among Ae. aegypti, Bombyx mori, and Drosophila melanogaster. However, the gene length and the intron-exon organization varied across the insect species. Expression analysis of the genes encoding RBPs using publicly available RNA sequencing data for different developmental time points of the mosquito life cycle starting from the ovary and eggs up to the adults revealed stage-specific expression with several genes preferentially expressed in early embryonic stages and blood-fed female ovaries. This is the first database for the Ae. aegypti RBPs that can serve as the reference base for future investigations. Stage-specific genes can be further explored to determine their role in mosquito growth and development with a focus on developing novel mosquito control strategies.

Advances in genome sequencing reveal changes in gene content that contribute to arthropod macroevolution

Current sequencing technology allows for the relatively affordable generation of highly contiguous genomes. Technological advances have made it possible for researchers to investigate the consequences of diverse sorts of genomic variants, such as gene gain and loss. With the extraordinary number of high-quality genomes now available, we take stock of how these genomic variants impact phenotypic evolution. We take care to point out that the identification of genomic variants of interest is only the first step in understanding their impact. Painstaking lab or fieldwork is still required to establish causal relationships between genomic variants and phenotypic evolution. We focus mostly on arthropod research, as this phylum has an impressive degree of phenotypic diversity and is also the subject of much evolutionary genetics research. This article is intended to both highlight recent advances in the field and also to be a primer for learning about evolutionary genetics and genomics.

De novo genome assembly of the invasive mosquito species Aedes japonicus and Aedes koreicus

BACKGROUND

Recently, two invasive Aedes mosquito species, Ae. japonicus and Ae. koreicus, are circulating in several European countries posing potential health risks to humans and animals. Vector control is the main option to prevent mosquito-borne diseases, and an accurate genome sequence of these mosquitoes is essential to better understand their biology and to develop effective control strategies.

METHODS

A de novo genome assembly of Ae. japonicus (Ajap1) and Ae. koreicus (Akor1) has been produced based on a hybrid approach that combines Oxford Nanopore long-read and Illumina short-read data. Their quality was ascertained using various metrics. Masking of repetitive elements, gene prediction and functional annotation was performed.

RESULTS

Sequence analysis revealed a very high presence of repetitive DNA and, among others, thermal adaptation genes and insecticide-resistance genes. Through the RNA-seq analysis of larvae and adults of Ae. koreicus and Ae. japonicus exposed to different temperatures, we also identified genes showing a differential temperature-dependent activation.

CONCLUSIONS

The assembly of Akor1 and Ajap1 genomes constitutes the first updated collective knowledge of the genomes of both mosquito species, providing the possibility of understanding key mechanisms of their biology such as the ability to adapt to harsh climates and to develop insecticide-resistance mechanisms.

Automated analysis of feeding behaviors of females of the mosquito Aedes aegypti using a modified flyPAD system

Mosquitoes present a global health challenge due to their ability to transmit human and animal pathogens upon biting and blood feeding. The investigation of tastants detected by mosquitoes and their associated feeding behaviors is needed to answer physiological and ecological questions that could lead to novel control methods. A high-throughput system originally developed for research in fruit flies feeding behavior, the flyPAD, was adapted and tested for behaviors associated with the interaction or consumption of liquid diets offered to females of the mosquito Aedes aegypti Liverpool strain. Females were given water, sucrose solution and sheep blood in choice and non-choice assays. The volume ingested was evaluated with fluorescein. The placement of the system on a heated surface allowed blood consumption, and without females puncturing a membrane. The flyPAD system recorded nine feeding behavioral variables, of which the number of sips and number of activity bouts correlated with meal volume ingested for both sucrose solution and blood. The adaptation to mosquitoes of the flyPAD system differentiated feeding behavior variables between two feeding deterrents, capsaicin, and caffeine. The flyPAD has potential to quickly assess diverse tastants in both sucrose and blood and may contribute to characterizing more precisely their mode of action.

Aedes aegypti microbiome composition covaries with the density of Wolbachia infection

BACKGROUND

Wolbachia is a widespread bacterial endosymbiont that can inhibit vector competency when stably transinfected into the mosquito, Aedes aegypti, a primary vector of the dengue virus (DENV) and other arboviruses. Although a complete mechanistic understanding of pathogen blocking is lacking, it is likely to involve host immunity induction and resource competition between Wolbachia and DENV, both of which may be impacted by microbiome composition. The potential impact of Wolbachia transinfection on host fitness is also of importance given the widespread release of mosquitos infected with the Drosophila melanogaster strain of Wolbachia (wMel) in wild populations. Here, population-level genomic data from Ae. aegypti was surveyed to establish the relationship between the density of wMel infection and the composition of the host microbiome.

RESULTS

Analysis of genomic data from 172 Ae. aegypti females across six populations resulted in an expanded and quantitatively refined, species-level characterization of the bacterial, archaeal, and fungal microbiome. This included 844 species of bacteria across 23 phyla, of which 54 species were found to be ubiquitous microbiome members across these populations. The density of wMel infection was highly variable between individuals and negatively correlated with microbiome diversity. Network analyses revealed wMel as a hub comprised solely of negative interactions with other bacterial species. This contrasted with the large and highly interconnected network of other microbiome species that may represent members of the midgut microbiome community in this population.

CONCLUSION

Our bioinformatic survey provided a species-level characterization of Ae. aegypti microbiome composition and variation. wMel load varied substantially across populations and individuals and, importantly, wMel was a major hub of a negative interactions across the microbiome. These interactions may be an inherent consequence of heightened pathogen blocking in densely infected individuals or, alternatively, may result from antagonistic Wolbachia-incompatible bacteria that could impede the efficacy of wMel as a biological control agent in future applications. The relationship between wMel infection variation and the microbiome warrants further investigation in the context of developing wMel as a multivalent control agent against other arboviruses. Video Abstract.

Decoupled evolution of the Sex Peptide gene family and Sex Peptide Receptor in Drosophilidae

Across internally fertilising species, males transfer ejaculate proteins that trigger wide-ranging changes in female behaviour and physiology. Much theory has been developed to explore the drivers of ejaculate protein evolution. The accelerating availability of high-quality genomes now allows us to test how these proteins are evolving at fine taxonomic scales. Here, we use genomes from 264 species to chart the evolutionary history of Sex Peptide (SP), a potent regulator of female post-mating responses in Drosophila melanogaster. We infer that SP first evolved in the Drosophilinae subfamily and has followed markedly different evolutionary trajectories in different lineages. Outside of the Sophophora-Lordiphosa, SP exists largely as a single-copy gene with independent losses in several lineages. Within the Sophophora-Lordiphosa, the SP gene family has repeatedly and independently expanded. Up to seven copies, collectively displaying extensive sequence variation, are present in some species. Despite these changes, SP expression remains restricted to the male reproductive tract. Alongside, we document considerable interspecific variation in the presence and morphology of seminal microcarriers that, despite the critical role SP plays in microcarrier assembly in D. melanogaster, appear to be independent of changes in the presence/absence or sequence of SP. We end by providing evidence that SP's evolution is decoupled from that of its receptor, SPR, in which we detect no evidence of correlated diversifying selection. Collectively, our work describes the divergent evolutionary trajectories that a novel gene has taken following its origin and finds a surprisingly weak coevolutionary signal between a supposedly sexually antagonistic protein and its receptor.

Targeting Sex Determination to Suppress Mosquito Populations

Each year, hundreds of millions of people are infected with arboviruses such as dengue, yellow fever, chikungunya, and Zika, which are all primarily spread by the notorious mosquito Aedes aegypti. Traditional control measures have proven insufficient, necessitating innovations. In response, here we generate a next generation CRISPR-based precision-guided sterile insect technique (pgSIT) for Aedes aegypti that disrupts genes essential for sex determination and fertility, producing predominantly sterile males that can be deployed at any life stage. Using mathematical models and empirical testing, we demonstrate that released pgSIT males can effectively compete with, suppress, and eliminate caged mosquito populations. This versatile species-specific platform has the potential for field deployment to effectively control wild populations of disease vectors.

Efficient sex separation by exploiting differential alternative splicing of a dominant marker in Aedes aegypti

Only female mosquitoes consume blood giving them the opportunity to transmit deadly human pathogens. Therefore, it is critical to remove females before conducting releases for genetic biocontrol interventions. Here we describe a robust sex-sorting approach termed SEPARATOR (Sexing Element Produced by Alternative RNA-splicing of A Transgenic Observable Reporter) that exploits sex-specific alternative splicing of an innocuous reporter to ensure exclusive dominant male-specific expression. Using SEPARATOR, we demonstrate reliable sex selection from early larval and pupal stages in Aedes aegypti, and use a Complex Object Parametric Analyzer and Sorter (COPAS) to demonstrate scalable high-throughput sex-selection of first instar larvae. Additionally, we use this approach to sequence the transcriptomes of early larval males and females and find several genes that are sex-specifically expressed. SEPARATOR can simplify mass production of males for release programs and is designed to be cross-species portable and should be instrumental for genetic biocontrol interventions.

Dishevelled Has Anti-Viral Activity in Rift Valley Fever Virus Infected Aedes aegypti

Mosquitoes in the genera Aedes and Culex are vectors of Rift Valley fever virus (RVFV), which emerges in periodic epidemics in Africa and Saudi Arabia. Factors that influence the transmission dynamics of RVFV are not well characterized. To address this, we interrogated mosquito host-signaling responses through analysis of differentially expressed genes (DEGs) in two mosquito species with marked differences in RVFV vector competence: Aedes aegypti (Aae, low competence) and Culex tarsalis (Cxt, high competence). Mosquito-host transcripts related to three different signaling pathways were investigated. Selected genes from the Wingless (Wg, WNT-beta-catenin) pathway, which is a conserved regulator of cell proliferation and differentiation, were assessed. One of these, dishevelled (DSH), differentially regulates progression/inhibition of the WNT and JNK (c-Jun N-terminal Kinase) pathways. A negative regulator of the JNK-signaling pathway, puckered, was also assessed. Lastly, Janus kinase/signal transducers and activators of transcription (JAK-STAT) are important for innate immunity; in this context, we tested domeless levels. Here, individual Aae and Cxt were exposed to RVFV MP-12 via oral bloodmeals and held for 14 days. Robust decreases in DEGs in both Aae and Cxt were observed. In particular, Aae DSH expression, but not Cxt DSH, was correlated to the presence/absence of viral RNA at 14 days post-challenge (dpc). Moreover, there was an inverse relationship between the viral copy number and aaeDSH expression. DSH silencing resulted in increased viral copy numbers compared to controls at 3 dpc, consistent with a role for aaeDSH in antiviral immunity. Analysis of cis-regulatory regions for the genes of interest revealed clues to upstream regulation of these pathways.

The evolution and mutational robustness of chromatin accessibility in Drosophila

BACKGROUND

The evolution of genomic regulatory regions plays a critical role in shaping the diversity of life. While this process is primarily sequence-dependent, the enormous complexity of biological systems complicates the understanding of the factors underlying regulation and its evolution. Here, we apply deep neural networks as a tool to investigate the sequence determinants underlying chromatin accessibility in different species and tissues of Drosophila.

RESULTS

We train hybrid convolution-attention neural networks to accurately predict ATAC-seq peaks using only local DNA sequences as input. We show that our models generalize well across substantially evolutionarily diverged species of insects, implying that the sequence determinants of accessibility are highly conserved. Using our model to examine species-specific gains in accessibility, we find evidence suggesting that these regions may be ancestrally poised for evolution. Using in silico mutagenesis, we show that accessibility can be accurately predicted from short subsequences in each example. However, in silico knock-out of these sequences does not qualitatively impair classification, implying that accessibility is mutationally robust. Subsequently, we show that accessibility is predicted to be robust to large-scale random mutation even in the absence of selection. Conversely, simulations under strong selection demonstrate that accessibility can be extremely malleable despite its robustness. Finally, we identify motifs predictive of accessibility, recovering both novel and previously known motifs.

CONCLUSIONS

These results demonstrate the conservation of the sequence determinants of accessibility and the general robustness of chromatin accessibility, as well as the power of deep neural networks to explore fundamental questions in regulatory genomics and evolution.

Current Status of Omics Studies Elucidating the Features of Reproductive Biology in Blood-Feeding Insects

Female insects belonging to the genera Anopheles, Aedes, Glossina, and Rhodnius account for the majority of global vector-borne disease mortality. In response to mating, these female insects undergo several molecular, physiological, and behavioral changes. Studying the dynamic post-mating molecular responses in these insects that transmit human diseases can lead to the identification of potential targets for the development of novel vector control methods. With the continued advancements in bioinformatics tools, we now have the capability to delve into various physiological processes in these insects. Here, we discuss the availability of multiple datasets describing the reproductive physiology of the common blood-feeding insects at the molecular level. Additionally, we compare the male-derived triggers transferred during mating to females, examining both shared and species-specific factors. These triggers initiate post-mating genetic responses in female vectors, affecting not only their reproductive success but also disease transmission.

Time-series analysis of transcriptomic changes due to permethrin exposure reveals that Aedes aegypti undergoes detoxification metabolism over 24 h

Insecticide resistance is a multifaceted response and an issue across taxa. Aedes aegypti, the mosquito that vectors Zika, dengue, chikungunya, and yellow fever, demonstrates high levels of pyrethroid resistance across the globe, presenting a challenge to public health officials. To examine the transcriptomic shifts across time after exposure to permethrin, a 3'Tag-Seq analysis was employed on samples 6, 10, and 24 h after exposure along with controls. Differential expression analysis revealed significant shifts in detoxifying enzymes and various energy-producing metabolic processes. These findings indicate significant alterations in gene expression associated with key energy mobilization pathways within the system. These changes encompass a coordinated response involving lipolysis, beta-oxidation, and the citric acid cycle, required for the production of energetic molecules such as ATP, NADH, NADPH, and FADH. These findings highlight a complex interplay of metabolic processes that may have broader implications for understanding insect physiology and response to environmental stimuli. Among the upregulated detoxifying enzymes are cytochrome P450s, glutathione s-transferases and peroxidases, and ATP-binding cassette transporters. Additionally, eight heat shock genes or genes with heat shock domains exhibit the highest fold change across time. Twenty-four hours after exposure, samples indicate a global downregulation of these processes, though principal component analysis suggests lasting signatures of the response. Understanding the recovery response to insecticide exposure provides information on possible new genetic and synergist targets to explore.

Eggs of the mosquito Aedes aegypti survive desiccation by rewiring their polyamine and lipid metabolism

Upon water loss, some organisms pause their life cycles and escape death. While widespread in microbes, this is less common in animals. Aedes mosquitoes are vectors for viral diseases. Aedes eggs can survive dry environments, but molecular and cellular principles enabling egg survival through desiccation remain unknown. In this report, we find that Aedes aegypti eggs, in contrast to Anopheles stephensi, survive desiccation by acquiring desiccation tolerance at a late developmental stage. We uncover unique proteome and metabolic state changes in Aedes embryos during desiccation that reflect reduced central carbon metabolism, rewiring towards polyamine production, and enhanced lipid utilisation for energy and polyamine synthesis. Using inhibitors targeting these processes in blood-fed mosquitoes that lay eggs, we infer a two-step process of desiccation tolerance in Aedes eggs. The metabolic rewiring towards lipid breakdown and dependent polyamine accumulation confers resistance to desiccation. Furthermore, rapid lipid breakdown is required to fuel energetic requirements upon water reentry to enable larval hatching and survival upon rehydration. This study is fundamental to understanding Aedes embryo survival and in controlling the spread of these mosquitoes.

An 11-point time course midgut transcriptome across 72 h after bloodfeeding provides detailed temporal resolution of transcript expression in the arbovirus vector, Aedes aegypti

As the major vector for dengue, Zika, yellow fever, and chikungunya viruses, the mosquito Aedes aegypti is one of the most important insects in public health. These viruses are transmitted by bloodfeeding, which is also necessary for the reproduction of the mosquito. Thus, the midgut plays an essential role in mosquito physiology as the center for bloodmeal digestion and as an organ that serves as the first line of defense against viruses. Despite its importance, transcriptomic dynamics with fine temporal resolution across the entire digestion cycle have not yet been reported. To fill this gap, we conducted a transcriptomic analysis of A. aegypti female midguts across a 72-h bloodmeal digestion cycle for 11 time points, with a particular focus on the first 24 h. PCA analysis confirmed that 72 h is indeed a complete digestion cycle. Cluster and GO enrichment analysis showed the orchestrated modulation of thousands of genes to accomplish the midgut's role as the center for digestion, as well as nutrient transport with a clear progression with sequential emphasis on transcription, translation, energy production, nutrient metabolism, transport, and finally, autophagy by 24-36 h. We further determined that many serine proteases are robustly expressed as if to prepare for unexpected physiological challenges. This study provides a powerful resource for the analysis of genomic features that coordinate the rapid and complex transcriptional program induced by mosquito bloodfeeding.

Multifaceted contributions of Dicer2 to arbovirus transmission by Aedes aegypti

Arthropod-borne viruses (arboviruses) transmitted by Aedes aegypti mosquitoes are an increasing threat to global health. The small interfering RNA (siRNA) pathway is considered the main antiviral immune pathway of insects, but its effective impact on arbovirus transmission is surprisingly poorly understood. Here, we use CRISPR-Cas9-mediated gene editing in vivo to mutate Dicer2, a gene encoding the RNA sensor and key component of the siRNA pathway. The loss of Dicer2 enhances early viral replication and systemic viral dissemination of four medically significant arboviruses (chikungunya, Mayaro, dengue, and Zika viruses) representing two viral families. However, Dicer2 mutants and wild-type mosquitoes display overall similar levels of vector competence. In addition, Dicer2 mutants undergo significant virus-induced mortality during infection with chikungunya virus. Together, our results define a multifaceted role for Dicer2 in the transmission of arboviruses by Ae. aegypti mosquitoes and pave the way for further mechanistic investigations.

Comparative Transcriptomic Analysis of Insecticide-Resistant Aedes aegypti from Puerto Rico Reveals Insecticide-Specific Patterns of Gene Expression

Aedes aegypti transmits major arboviruses of public health importance, including dengue, chikungunya, Zika, and yellow fever. The use of insecticides represents the cornerstone of vector control; however, insecticide resistance in Ae. aegypti has become widespread. Understanding the molecular basis of insecticide resistance in this species is crucial to design effective resistance management strategies. Here, we applied Illumina RNA-Seq to study the gene expression patterns associated with resistance to three widely used insecticides (malathion, alphacypermethrin, and lambda-cyhalothrin) in Ae. aegypti populations from two sites (Manatí and Isabela) in Puerto Rico (PR). Cytochrome P450s were the most overexpressed detoxification genes across all resistant phenotypes. Some detoxification genes (CYP6Z7, CYP28A5, CYP9J2, CYP6Z6, CYP6BB2, CYP6M9, and two CYP9F2 orthologs) were commonly overexpressed in mosquitoes that survived exposure to all three insecticides (independent of geographical origin) while others including CYP6BY1 (malathion), GSTD1 (alpha-cypermethrin), CYP4H29 and GSTE6 (lambda-cyhalothrin) were uniquely overexpressed in mosquitoes that survived exposure to specific insecticides. The gene ontology (GO) terms associated with monooxygenase, iron binding, and passive transmembrane transporter activities were significantly enriched in four out of six resistant vs. susceptible comparisons while serine protease activity was elevated in all insecticide-resistant groups relative to the susceptible strain. Interestingly, cuticular-related protein genes (chinase and chitin) were predominantly downregulated, which was also confirmed in the functional enrichment analysis. This RNA-Seq analysis presents a detailed picture of the candidate detoxification genes and other pathways that are potentially associated with pyrethroid and organophosphate resistance in Ae. aegypti populations from PR. These results could inform development of novel molecular tools for detection of resistance-associated gene expression in this important arbovirus vector and guide the design and implementation of resistance management strategies.

The evolution of borneol repellency in culicine mosquitoes

Insects have developed remarkable adaptations to effectively interact with plant secondary metabolites and utilize them as cues to identify suitable hosts. Consequently, humans have used aromatic plants for centuries to repel mosquitoes. The repellent effects of plant volatile compounds are mediated through olfactory structures present in the antennae, and maxillary palps of mosquitoes. Mosquito maxillary palps contain capitate-peg sensilla, which house three olfactory sensory neurons, of which two are mainly tuned to either carbon dioxide or octenol - two animal host odorants. However, the third neuron, which expresses the OR49 receptor, has remained without a known ecologically-relevant odorant since its initial discovery. In this study, we used odorant mixtures and terpenoid-rich Cannabis essential oils to investigate the activation of OR49. Our results demonstrate that two monoterpenoids, borneol and camphor, selectively activate OR49, and OR9-expressing neurons, as well as the MD3 glomerulus in the antennal lobe. We confirm that borneol repels female mosquitoes, and knocking out the gene encoding the OR49 receptor suppresses the response of the corresponding olfactory sensory neuron. Importantly, this molecular mechanism of action is conserved across culicine mosquito species, underscoring its significance in their olfactory systems.

In Silico Characterisation of the Aedes aegypti Gustatory Receptors

Aedes aegypti, also known as the dengue mosquito or the yellow fewer mosquito, is the vector of dengue, chikungunya, Zika, Mayaro and yellow fever viruses. The A. aegypti genome contains an array of gustatory receptor (GR) proteins that are related to the recognition of taste. In this study, we performed in silico molecular characterization of all 72 A. aegypti GRs reported in the latest version of A. aegypti genome AaegL5. Phylogenetic analysis classified the receptors into three major clads. Multiple GRs were found to encode multiple transcripts. Physicochemical attributes such as the aliphatic index, hydropathicity index and isoelectric point indicated that A. aegypti gustatory receptors are highly stable and are tailored to perform under a variety of cellular environments. Analysis for subcellular localization indicated that all the GRs are located either in the extracellular matrix or the plasma membrane. Results also indicated that the GRs are distributed mainly on chromosomes 2 and 3, which house 22 and 49 GRs, respectively, whereas chromosome 1 houses only one GR. NCBI-CDD analysis showed the presence of a highly conserved 7tm_7 chemosensory receptor protein superfamily that includes gustatory and odorant receptors from insect species Anopheles gambiae and Drosophila melanogaster. Further, three significantly enriched ungapped motifs in the protein sequence of all 72 A. aegypti gustatory receptors were found. High-quality 3D models for the tertiary structures were predicted with significantly higher confidence, along with ligand-binding residues. Prediction of S-nitrosylation sites indicated the presence of target cysteines in all the GRs with close proximity to the ligand-bindings sites within the 3D structure of the receptors. In addition, two highly conserved motifs inside the GR proteins were discovered that house a tyrosine (Y) and a cysteine (C) residue which may serve as targets for NO-mediated tyrosine nitration and S-nitrosylation, respectively. This study will help devise strategies for functional genomic studies of these important receptor molecules in A. aegypti and other mosquito species through in vitro and in vivo studies.

Midgut transcriptomic responses to dengue and chikungunya viruses in the vectors Aedes albopictus and Aedes malayensis

Dengue (DENV) and chikungunya (CHIKV) viruses are among the most preponderant arboviruses. Although primarily transmitted through the bite of Aedes aegypti mosquitoes, Aedes albopictus and Aedes malayensis are competent vectors and have an impact on arbovirus epidemiology. Here, to fill the gap in our understanding of the molecular interactions between secondary vectors and arboviruses, we used transcriptomics to profile the whole-genome responses of A. albopictus to CHIKV and of A. malayensis to CHIKV and DENV at 1 and 4 days post-infection (dpi) in midguts. In A. albopictus, 1793 and 339 genes were significantly regulated by CHIKV at 1 and 4 dpi, respectively. In A. malayensis, 943 and 222 genes upon CHIKV infection, and 74 and 69 genes upon DENV infection were significantly regulated at 1 and 4 dpi, respectively. We reported 81 genes that were consistently differentially regulated in all the CHIKV-infected conditions, identifying a CHIKV-induced signature. We identified expressed immune genes in both mosquito species, using a de novo assembled midgut transcriptome for A. malayensis, and described the immune architectures. We found the JNK pathway activated in all conditions, generalizing its antiviral function to Aedines. Our comprehensive study provides insight into arbovirus transmission by multiple Aedes vectors.