Hidden genomic features of an invasive malaria vector, Anopheles stephensi, revealed by a chromosome-level genome assembly

Genome structural variants shape adaptive success of an invasive urban malaria vector Anopheles stephensi

Global changes are associated with the emergence of several invasive species. However, the genomic determinants of the adaptive success of an invasive species in a new environment remain poorly understood. Genomic structural variants (SVs), consisting of copy number variants, play an important role in adaptation. SVs often cause large adaptive shifts in ecologically important traits, which makes SVs compelling candidates for driving rapid adaptations to environmental changes, which is critical to invasive success. To address this problem, we investigated the role SVs play in the adaptive success of Anopheles stephensi , a primary vector of urban malaria in South Asia and an invasive malaria vector in several South Asian islands and Africa. We collected whole genome sequencing data from 115 mosquitoes from invasive island populations and four locations from mainland India, an ancestral range for the species. We identified 2,988 duplication copy number variants and 16,038 deletions in these strains, with ∼50% overlapping genes. SVs are enriched in genomic regions with signatures of selective sweeps in the mainland and invasive island populations, implying a putative adaptive role of SVs. Nearly all high-frequency SVs, including the candidate adaptive variants, in the invasive island populations are present on the mainland, suggesting a major contribution of existing variation to the success of the island populations. Among the candidate adaptive SVs, three duplications involving toxin-resistance genes evolved, likely due to the widespread application of insecticides in India since the 1950s. We also identify two SVs associated with the adaptation of An. stephensi larvae to brackish water in the island and two coastal mainland populations, where the mutations likely originated. Our results suggest that existing SVs play a vital role in the evolutionary success of An. stephensi in new environmental conditions.

Establishing a Male-Positive Genetic Sexing Strain in the Asian Malaria Vector Anopheles stephensi

Genetic biocontrol interventions targeting mosquito-borne diseases require the release of male mosquitoes exclusively, as only females consume blood and transmit human pathogens. This reduces the risk of spreading pathogens while enabling effective population control. Robust sex sorting methods to enable early larval sorting in mosquitoes need to be developed to allow for scalable sex sorting for genetic biocontrol interventions. This study applies the SEPARATOR (Sexing Element Produced by Alternative RNA-splicing of A Transgenic Observable Reporter) system, previously developed for Aedes aegypti, to the Asian malaria vector Anopheles stephensi. We hypothesized that the intron from the doublesex gene in Anopheles gambiae would function in An. stephensi due to evolutionary conservation. Our results confirm that the splicing module from An. gambiae operates effectively in An. stephensi, demonstrating evolutionary conservation in sex-specific splicing events between these species. This system enables reliable positive male selection from first instar larval to pupal stages. RT-PCR analysis demonstrates that male-specific EGFP expression is dependent on doublesex sex-specific splicing events. The SEPARATOR system's independence from sex-chromosome linkage confers resistance to meiotic recombination and chromosomal rearrangements. This approach may facilitate the mass release of males, and the cross-species portability of SEPARATOR establishes it as a valuable tool for genetic biocontrol interventions across various pest species.

A decade of invasive Anopheles stephensi sequence-based identification: toward a global standard

Anopheles stephensi is an invasive malaria vector in Africa that has been implicated in malaria outbreaks in the Horn of Africa. In 10 years, it has been detected as far east as Djibouti and as far west as Ghana. Early detections were mostly incidental, but now active surveillance in Africa has been updated to include An. stephensi. Morphological identification of An. stephensi from native vectors can be challenging, thus, sequence-based assays have been used to confirm identification during initial detections. Methods of sequence-based identification of An. stephensi have varied across initial detections to date. Here, we summarize initial detections, make suggestions that could provide a standardized approach, and discuss how sequences can inform additional genomic studies beyond species identification.

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.

Detection of Spatially Localized Sounds Is Robust to Saccades and Concurrent Eye Movement-Related Eardrum Oscillations (EMREOs)

Hearing is an active process, and recent studies show that even the ear is affected by cognitive states or motor actions. One example are movements of the eardrum induced by saccadic eye movements, known as "eye movement-related eardrum oscillations" (EMREOs). While these are systematically shaped by the direction and size of saccades, the consequences of saccadic eye movements and their resulting EMREOs for hearing remain unclear. We here studied their implications for the detection of near-threshold clicks in human participants. Across three experiments, sound detection was not affected by their time of presentation relative to saccade onset, by saccade amplitude or direction. While the EMREOs were shaped by the direction and amplitude of the saccadic movement, inducing covert shifts in spatial attention did not affect the EMREO, suggesting that this signature of active sensing is restricted to overt changes in visual focus. Importantly, in our experiments, fluctuations in the EMREO amplitude were not related to detection performance, at least when monaural cues are sufficient. Hence, while eye movements may shape the transduction of acoustic information, the behavioral implications remain to be understood.SIGNIFICANCE STATEMENT Previous studies suggest that oculomotor behavior may influence how we perceive spatially localized sounds. Recent work has introduced a new perspective on this question by showing that eye movements can directly modulate the eardrum. Yet, it remains unclear whether this signature of active hearing accounts for behavioral effects. We here show that overt but not covert changes in visual attention modulate the eardrum, but these modulations do not interfere with the detection of sounds. Our results provide a starting point to obtain a deeper understanding about the interplay of oculomotor behavior and the active ear.

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.

Tracking Adaptive Pathways of Invasive Insects: Novel Insight from Genomics

Despite the huge human and economic costs of invasive insects, which are the main group of invasive species, their environmental impacts through various mechanisms remain inadequately explained in databases and much of the invasion biology literature. High-throughput sequencing technology, especially whole-genome sequencing, has been used as a powerful method to study the mechanisms through which insects achieve invasion. In this study, we reviewed whole-genome sequencing-based advances in revealing several important invasion mechanisms of invasive insects, including (1) the rapid genetic variation and evolution of invasive populations, (2) invasion history and dispersal paths, (3) rapid adaptation to different host plant ranges, (4) strong environmental adaptation, (5) the development of insecticide resistance, and (6) the synergistic damage caused by invasive insects and endosymbiotic bacteria. We also discussed prevention and control technologies based on whole-genome sequencing and their prospects.

Identification of two insecticide resistance markers in Ethiopian Anopheles stephensi mosquitoes using a multiplex amplicon sequencing assay

Since its first detection in 2012 in Djibouti, Anopheles stephensi has invaded and established in the Horn of Africa, and more recently Nigeria. The expansion of this vector poses a significant threat to malaria control and elimination efforts. Integrated vector management is the primary strategy used to interrupt disease transmission; however, growing insecticide resistance is threatening to reverse gains in global malaria control. We present a next-generation amplicon-sequencing approach, for high-throughput monitoring of insecticide resistance genes (ace1, GSTe2, vgsc and rdl), species identification and characterization of genetic diversity (its2 and cox1) in An. stephensi. Ninety-five An. stephensi mosquitoes, collected in Ethiopia, were screened, identifying 104 SNPs, including the knock-down mutation L958F (L1014F in Musca domestica), and for the first time in this vector species, the A296S substitution (A301S in Drosophila melanogaster) in the rdl locus. Two other amino acid substitutions (ace1-N177D, GSTe2-V189L) were also identified but have not been previously implicated in insecticide resistance. Genetic diversity in the mitochondrial cox1 gene revealed shared haplotypes between Ethiopian An. stephensi with samples from Pakistan, Sudan, and Djibouti. Overall, we present a reliable, cost-effective strategy using amplicon-sequencing to monitor known insecticide resistance mutations, with the potential to identify new genetic variants, to assist in the high-throughput surveillance of insecticide resistance in An. stephensi populations.

Mapping current and future thermal limits to suitability for malaria transmission by the invasive mosquito Anopheles stephensi

BACKGROUND

Anopheles stephensi is a malaria-transmitting mosquito that has recently expanded from its primary range in Asia and the Middle East, to locations in Africa. This species is a competent vector of both Plasmodium falciparum and Plasmodium vivax malaria. Perhaps most alarming, the characteristics of An. stephensi, such as container breeding and anthropophily, make it particularly adept at exploiting built environments in areas with no prior history of malaria risk.

METHODS

In this paper, global maps of thermal transmission suitability and people at risk (PAR) for malaria transmission by An. stephensi were created, under current and future climate. Temperature-dependent transmission suitability thresholds derived from recently published species-specific thermal curves were used to threshold gridded, monthly mean temperatures under current and future climatic conditions. These temperature driven transmission models were coupled with gridded population data for 2020 and 2050, under climate-matched scenarios for future outcomes, to compare with baseline predictions for 2020 populations.

RESULTS

Using the Global Burden of Disease regions approach revealed that heterogenous regional increases and decreases in risk did not mask the overall pattern of massive increases of PAR for malaria transmission suitability with An. stephensi presence. General patterns of poleward expansion for thermal suitability were seen for both P. falciparum and P. vivax transmission potential.

CONCLUSIONS

Understanding the potential suitability for An. stephensi transmission in a changing climate provides a key tool for planning, given an ongoing invasion and expansion of the vector. Anticipating the potential impact of onward expansion to transmission suitable areas, and the size of population at risk under future climate scenarios, and where they occur, can serve as a large-scale call for attention, planning, and monitoring.

Diversity, duplication, and genomic organization of homeobox genes in Lepidoptera

Homeobox genes encode transcription factors with essential roles in patterning and cell fate in developing animal embryos. Many homeobox genes, including Hox and NK genes, are arranged in gene clusters, a feature likely related to transcriptional control. Sparse taxon sampling and fragmentary genome assemblies mean that little is known about the dynamics of homeobox gene evolution across Lepidoptera or about how changes in homeobox gene number and organization relate to diversity in this large order of insects. Here we analyze an extensive data set of high-quality genomes to characterize the number and organization of all homeobox genes in 123 species of Lepidoptera from 23 taxonomic families. We find most Lepidoptera have around 100 homeobox loci, including an unusual Hox gene cluster in which the lab gene is repositioned and the ro gene is next to pb A topologically associating domain spans much of the gene cluster, suggesting deep regulatory conservation of the Hox cluster arrangement in this insect order. Most Lepidoptera have four Shx genes, divergent zen-derived loci, but these loci underwent dramatic duplication in several lineages, with some moths having over 165 homeobox loci in the Hox gene cluster; this expansion is associated with local LINE element density. In contrast, the NK gene cluster content is more stable, although there are differences in organization compared with other insects, as well as major rearrangements within butterflies. Our analysis represents the first description of homeobox gene content across the order Lepidoptera, exemplifying the potential of newly generated genome assemblies for understanding genome and gene family evolution.

Visualization of the Linear and Spatial Organization of Chromosomes in Mosquitoes

Mosquitoes are vectors of dangerous human diseases such as malaria, dengue, Zika, West Nile fever, and lymphatic filariasis. Visualization of the linear and spatial organization of mosquito chromosomes is important for understanding genome structure and function. Utilization of chromosomal inversions as markers for population genetics studies yields insights into mosquito adaptation and evolution. Cytogenetic approaches assist with the development of chromosome-scale genome assemblies that are useful tools for studying mosquito biology and for designing novel vector control strategies. Fluorescence in situ hybridization is a powerful technique for localizing specific DNA sequences within the linear chromosome structure and within the spatial organization of the cell nucleus. Here, we introduce protocols used in our laboratories for chromosome visualization and their application in mosquitoes.

Tandem duplication of a genomic region encoding glutathione S-transferase epsilon-2 and -4 genes in DDT-resistant Anopheles stephensi strain from India

The glutathione S-transferases (GST) genes are a multigene family of enzymes involved in the metabolism of endogenous and xenobiotic compounds by catalysing the conjugation of the reduced form of glutathione to the substrate. The epsilon class of GST (GSTe), unique to arthropods, is known to be involved in the detoxification process of several classes of insecticides, and GSTe2 in particular is known to have DDT dehydrochlorinase activity. This communication reports a tandem duplication of a genomic region encoding GSTe2 and GSTe4 genes in a laboratory-colonized DDT-resistant Anopheles stephensi. We identified duplication breakpoints and the organization of gene duplication through Sanger sequencing performed on long-PCR products. Manual annotation of sequences revealed a tandemly-arrayed duplication of a 3.62 kb segment of GST epsilon gene clusters comprised of five genes: a partial GSTe1, GSTe2, GSTe2-pseudogene, GSTe4 and partial GSTe5, interconnected by a conserved 2.42 kb DNA insert segment major part of which is homologous to a genomic region located on a different chromosome. The tandemly duplicated array contained a total of two GSTe2 and three GSTe4 functional paralog genes. Read-depth coverage and split-read analysis of Illumina-based whole-genome sequence reads confirmed the presence of duplication in the corresponding region of the genome. The increased gene dose in mosquitoes as a result of the GSTe gene-duplication may be an adaptive process to increase levels of detoxifying enzymes to counter insecticide pressure.

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.

Selective targeting of biting females to control mosquito-borne infectious diseases

Mosquitoes are vectors for a number of infectious diseases. Only females feed on blood to provision for their embryos and, in doing so, transmit pathogens to the associated vertebrate hosts. Therefore, sex is an important phenotype in the context of genetic control programs, both for sex separation in the rearing facilities to avoid releasing biting females and for ways to distort the sex ratio towards nonbiting males. We review recent progress in the fundamental knowledge of sex determination and sex chromosomes in mosquitoes and discuss new methods to achieve sex separation and sex ratio distortion to help control mosquito-borne infectious diseases. We conclude by suggesting a few critical areas for future research.

Transcriptomic and small RNA response to Mayaro virus infection in Anopheles stephensi mosquitoes

Mayaro virus (MAYV) is an arboviral pathogen in the genus Alphavirus that is circulating in South America with potential to spread to naïve regions. MAYV is also one of the few viruses with the ability to be transmitted by mosquitoes in the genus Anopheles, as well as the typical arboviral transmitting mosquitoes in the genus Aedes. Few studies have investigated the infection response of Anopheles mosquitoes. In this study we detail the transcriptomic and small RNA responses of An. stephensi to infection with MAYV via infectious bloodmeal at 2, 7, and 14 days post infection (dpi). 487 unique transcripts were significantly regulated, 78 putative novel miRNAs were identified, and an siRNA response is observed targeting the MAYV genome. Gene ontology analysis of transcripts regulated at each timepoint shows a number of proteases regulated at 2 and 7 dpi, potentially representative of Toll or melanization pathway activation, and repression of pathways related to autophagy and apoptosis at 14 dpi. These findings provide a basic understanding of the infection response of An. stephensi to MAYV and help to identify host factors which might be useful to target to inhibit viral replication in Anopheles mosquitoes.

Mayaro virus (MAYV) is a mosquito-borne Alphavirus responsible for outbreaks in South America and the Caribbean. In this study we infected Anopheles stephensi with MAYV and sequenced mRNA and small RNA to understand how MAYV infection impacts gene transcription and the expression of small RNAs in the mosquito vector. Genes involved with innate immunity and signaling pathways related to cell death are regulated in response to MAYV infection of An. stephensi, we also discovered novel miRNAs and describe the expression patterns of miRNAs, siRNAs, and piRNAs following bloodmeal ingestion. These results suggest that MAYV does induce a molecular response to infection in its mosquito vector species.

Anopheles mosquitoes reveal new principles of 3D genome organization in insects

Chromosomes are hierarchically folded within cell nuclei into territories, domains and subdomains, but the functional importance and evolutionary dynamics of these hierarchies are poorly defined. Here, we comprehensively profile genome organizations of five Anopheles mosquito species and show how different levels of chromatin architecture influence each other. Patterns observed on Hi-C maps are associated with known cytological structures, epigenetic profiles, and gene expression levels. Evolutionary analysis reveals conservation of chromatin architecture within synteny blocks for tens of millions of years and enrichment of synteny breakpoints in regions with increased genomic insulation. However, in-depth analysis shows a confounding effect of gene density on both insulation and distribution of synteny breakpoints, suggesting limited causal relationship between breakpoints and regions with increased genomic insulation. At the level of individual loci, we identify specific, extremely long-ranged looping interactions, conserved for ~100 million years. We demonstrate that the mechanisms underlying these looping contacts differ from previously described Polycomb-dependent interactions and clustering of active chromatin.

The genome trilogy of Anopheles stephensi, an urban malaria vector, reveals structure of a locus associated with adaptation to environmental heterogeneity

Anopheles stephensi is the most menacing malaria vector to watch for in newly urbanising parts of the world. Its fitness is reported to be a direct consequence of the vector adapting to laying eggs in over-head water tanks with street-side water puddles polluted by oil and sewage. Large frequent inversions in the genome of malaria vectors are implicated in adaptation. We report the genome assembly of a strain of An. stephensi of the type-form, collected from a construction site from Chennai (IndCh) in 2016. The genome reported here with a L50 of 4, completes the trilogy of high-resolution genomes of strains with respect to a 16.5 Mbp 2Rb genotype in An. stephensi known to be associated with adaptation to environmental heterogeneity. Unlike the reported genomes of two other strains, STE2 (2R+^b/2Rb) and UCI (2Rb/2Rb), IndCh is found to be homozygous for the standard form (2R+^b/2R+^b). Comparative genome analysis revealed base-level details of the breakpoints and allowed extraction of 22,650 segregating SNPs for typing this inversion in populations. Whole genome sequencing of 82 individual mosquitoes from diverse geographical locations reveal that one third of both wild and laboratory populations maintain the heterozygous genotype of 2Rb. The large number of SNPs can be tailored to 1740 exonic SNPs enabling genotyping directly from transcriptome sequencing. The genome trilogy approach accelerated the study of fine structure and typing of an important inversion in An. stephensi, putting the genome resources for this understudied species on par with the extensively studied malaria vector, Anopheles gambiae. We argue that the IndCh genome is relevant for field translation work compared to those reported earlier by showing that individuals from diverse geographical locations cluster with IndCh, pointing to significant convergence resulting from travel and commerce between cities, perhaps, contributing to the survival of the fittest strain.

Spontaneous mutation rate estimates for the principal malaria vectors Anopheles coluzzii and Anopheles stephensi

Using high-depth whole genome sequencing of F0 mating pairs and multiple individual F1 offspring, we estimated the nuclear mutation rate per generation in the malaria vectors Anopheles coluzzii and Anopheles stephensi by detecting de novo genetic mutations. A purpose-built computer program was employed to filter actual mutations from a deep background of superficially similar artifacts resulting from read misalignment. Performance of filtering parameters was determined using software-simulated mutations, and the resulting estimate of false negative rate was used to correct final mutation rate estimates. Spontaneous mutation rates by base substitution were estimated at 1.00 × 10⁻⁹ (95% confidence interval, 2.06 × 10⁻¹⁰—2.91 × 10⁻⁹) and 1.36 × 10⁻⁹ (95% confidence interval, 4.42 × 10⁻¹⁰—3.18 × 10⁻⁹) per site per generation in A. coluzzii and A. stephensi respectively. Although similar studies have been performed on other insect species including dipterans, this is the first study to empirically measure mutation rates in the important genus Anopheles, and thus provides an estimate of µ that will be of utility for comparative evolutionary genomics, as well as for population genetic analysis of malaria vector mosquito species.

Extending and Running the Mosquito Small RNA Genomics Resource Pipeline

The Mosquito Small RNA Genomics (MSRG) resource is a repository of analyses on the small RNA transcriptomes of mosquito cell cultures and somatic and gonadal tissues. This resource allows for comparing the regulation dynamics of small RNAs generated from transposons and viruses across mosquito species. This chapter covers the procedures to set up the MSRG resource pipeline as a new installation by detailing the necessary collection of genome reference and annotation files and lists of microRNAs (miRNAs) hairpin sequences, transposon repeats consensus sequences, and virus genome sequences. Proper execution of the MSRG resource pipeline yields outputs amenable to biologists to further analyze with desktop and spreadsheet software to gain insights into the balance between arthropod endogenous small RNA populations and the proportions of virus-derived small RNAs that include Piwi-interacting RNAs (piRNAs) and endogenous small interfering RNAs (siRNAs).