Genetic Variation in Insect Vectors: Death of Typology?

Anticipating evolutionary responses of mosquito mating systems to population suppression with mass-reared males

Strategies that rely on the mass release of males to suppress mosquito populations will exert selective pressure on natural mating systems. Here, we investigate how mass releases might affect the mating behaviors of wild target populations. We highlight gaps in our understanding of both variation in these aspects of mosquito behavior and the evolutionary forces that maintain variation within and between populations. We provide a mathematical framework for integrating mosquito mating ecology into models of population suppression. Given that these strategies are being increasingly deployed, anticipating and managing evolutionary responses of target population behavior should be a priority for research.

Insights into the mitochondrial cytochrome oxidase I (mt-COI) gene and wing morphometrics of Anopheles baimaii (Diptera: Culicidae) in malaria-endemic islands of Thailand

Anopheles baimaii (Diptera: Culicidae) significantly contributes to the transmission of parasites causing malaria in Southeast Asia and South Asia. This study examined the morphological (wing shape) and molecular (mitochondrial gene) variations of An. baimaii in four of Thailand's border islands, and also investigated the presence of Plasmodium parasites in these mosquitoes. No Plasmodium infections were detected in the samples. Significant differences in wing shape were observed in most island populations (p < 0.05). A single-linkage tree, constructed using Mahalanobis distances, clustered the populations into two groups based on geographical locations. Genetic variation in An. baimaii was also analyzed through cytochrome c oxidase subunit I (COI) gene sequences. This analysis identified 22 segregating sites and a low nucleotide diversity of 0.004. Furthermore, 18 distinct haplotypes were identified, indicating a high haplotype diversity of 0.825. Neutrality tests for the overall population revealed a significantly negative Fu's Fs value (-5.029), indicating a population expansion. In contrast, Tajima's D yielded a negative value (-1.028) that did not reach statistical significance. The mismatch distribution analysis exhibited a bimodal pattern, and the raggedness index was 0.068, showing no significant discrepancy (p = 0.485) between observed and expected distributions. Pairwise genetic differentiation assessments demonstrated significant differences between all populations (p < 0.05). These findings provide valuable insights into the COI gene and wing morphometric variations in An. baimaii across Thailand's islands, offering critical information for understanding the adaptations of this malaria vector and guiding future comprehensive research.

Genetic diversity and phylogeography of Phlebotomus argentipes (Diptera: Psychodidae, Phlebotominae), using COI and ND4 mitochondrial gene sequences

BACKGROUND

Phlebotomus argentipes complex is the primary vector for cutaneous leishmaniasis, a burgeoning health concern in contemporary Sri Lanka, where effective vector control is important for proper disease management. Understanding the genetic diversity of the P. argentipes population in Sri Lanka is vital before implementing a successful vector control program. Various studies have indicated that genetic divergence, caused by genetic drift or selection, can significantly influence the vector capacity of arthropod species. To devise innovative control strategies for P. argentipes, exploring genetic diversity and phylogeography can offer valuable insights into vector competence, key genetic trait transfer, and impact on disease epidemiology. The primary objective is to analyze the genetic diversity and phylogeography of the P. argentipes complex in Sri Lanka, based on two mitochondrial genomic regions in modern representatives of P. argentipes populations.

METHODOLOGY

A total of 159 P. argentipes specimens were collected from five endemic areas of cutaneous leishmaniasis and identified morphologically. Two mitochondrial regions (Cytochrome c oxidase subunit I (COI) and NADH dehydrogenase subunit 4 (ND4) were amplified using the total DNA and subsequently sequenced. Partial sequences of those mitochondrial genes were utilized to analyze genetic diversity indices and to explore phylogenetic and phylogeographic relationships.

PRINCIPAL FINDINGS

Among five sampling locations, the highest genetic diversity for COI and ND4 was observed in Hambantota (Hd-0.749, π-0.00417) and Medirigiriya (Hd-0.977, π-0.01055), respectively. Phylogeographic analyses conducted using COI sequences and GenBank retrieved sequences demonstrated a significant divergence of P. argentipes haplotypes found in Sri Lanka. Results revealed that they have evolved from the Indian ancestral haplotype due to historical- geographical connections of the Indian subcontinent with Sri Lanka.

CONCLUSIONS

Utilizing high-mutation-rate mitochondrial genes, such as ND4, can enhance the accuracy of genetic variability analysis in P. argentipes populations in Sri Lanka. The phylogeographical analysis of COI gene markers in this study provides insights into the historical geographical relationship between India and P. argentipes in Sri Lanka. Both COI and ND4 genes exhibited consistent genetic homogeneity in P. argentipes in Sri Lanka, suggesting minimal impact on gene flow. This homogeneity also implies the potential for horizontal gene transfer across populations, facilitating the transmission of genes associated with traits like insecticide resistance. This dynamic undermines disease control efforts reliant on vector control strategies.

Shedding light on the controversial taxonomic status of Culicoides jamaicensis and Culicoides paolae (Diptera: Ceratopogonidae): an overseas trip among continents

Culicoides biting midges (Diptera: Ceratopogonidae) are small bloodsucking flies that act as vectors for various pathogens of medical and veterinary importance. This study aimed to examine, using a comprehensive approach, the controversial taxonomic status of 2 Culicoides species that are currently distributed in the Neotropical (Culicoides jamaicensis Edwards) and Palearctic (Culicoides paolae Boorman) areas and possess unique and distinctive features. Previous investigations based on morphological analysis have suggested that these 2 species may be synonyms. Our work updated the current geographical distribution of both species and analyzed new specimens from different geographic origins, together with publicly available sequences. We used 2 universal genetic markers (COI and 28S) to test this hypothesis. Our study reveals evidence that C. paolae and C. jamaicensis belong to the same species due to the following statements: (i) similar morphological features; (ii) low interspecific genetic variation; (iii) association with a single genetic cluster; (iv) inclusion within the subgenus Drymodesmyia, which has only been recorded in the New World; and (v) occurrence in habitats with moderate temperatures. We recommend that European and African specimens of C. paolae be considered from now on as C. jamaicensis. Our comprehensive approach shed new light on the taxonomic status of these 2 Culicoides species and has implications for future studies on their biology and ecology.

Genetic Diversity and Population Genetic Structure of Aedes albopictus in the Yangtze River Basin, China

Aedes albopictus is an indigenous primary vector of dengue and Zika viruses in China. Understanding the population spatial genetic structure, migration, and gene flow of vector species is critical to effectively preventing and controlling vector-borne diseases. The genetic variation and population structure of Ae. albopictus populations collected from 22 cities along the Yangtze River Basin were investigated with nine microsatellite loci and the mitochondrial CoxI gene. The polymorphic information content (PIC) values ranged from 0.534 to 0.871. The observed number of alleles (Na) values ranged from 5.455 to 11.455, and the effective number of alleles (Ne) values ranged from 3.106 to 4.041. The Shannon Index (I) ranged from 1.209 to 1.639. The observed heterozygosity (Ho) values ranged from 0.487 to 0.545. The FIS value ranged from 0.047 to 0.212. All Ae. albopictus populations were adequately allocated to three clades with significant genetic differences. Haplotype 2 is the most primitive molecular type and forms 26 other haplotypes after one or more site mutations. The rapid expansion of high-speed rail, aircraft routes and highways along the Yangtze River Basin have accelerated the dispersal and communication of mosquitoes, which appears to have contributed to inhibited population differentiation and promoted genetic diversity among Ae. albopictus populations.

Modifying mosquitoes to suppress disease transmission: Is the long wait over?

For more than 50 years it has been a dream of medical entomologists and public health workers to control diseases like malaria and dengue fever by modifying, through genetics and other methods, the arthropods that transmit them to humans. A brief synopsis of the history of these efforts as applied to mosquitoes is presented; none proved to be effective in reducing disease prevalence. Only in the last few years have novel approaches been developed or proposed that indicate the long wait may be over. Three recent developments are particularly promising: CRISPR-Cas9 driven genetic modification, shifting naturally occurring allele frequencies, and microbe-based modifications. The last is the furthest along in implementation. Dengue fever incidence has been reduced between 40% and 96% in 4 different regions of the world where Wolbachia-infected Aedes aegypti have been established in the field. It is not yet clear how sustainable such control programs will prove to be, but there is good reason for optimism. In light of this, the time is ripe for reinvigorated research on vectors, especially genetics. Vector-borne diseases primarily affect under-developed countries and thus have not received the attention they deserve from wealthier countries with well-developed and funded biomedical research establishments.

RP, G. S, K. AAP, K. SM, A. SP, P. R, V. S, Dasgupta S, Krishnan J, Ishtiaq F. Island biogeography and human practices drive ecological connectivity in mosquito species richness in the Lakshadweep Archipelago

Mosquitoes are globally distributed and adapted to a broad range of environmental conditions. As obligatory hosts of many infectious pathogens, mosquito abundance and distribution are primarily determined by the presence and quality of larval habitats. To understand the dynamics and productivity of larval habitats in changing island environments, we conducted a four-month mosquito survey across ten inhabited islands in the Lakshadweep archipelago. Using fine-resolution larval habitat mapping, we recorded 7890 mosquitoes representing 13 species and 7 genera. Of these, four species comprised 95% of the total collections—Aedes albopictus (Stegomyia) was the dominant species followed by Armigeres subalbatus, Culex quinquefasciatus and Malaya genurostris. We found larval species richness was positively associated with the island area and mosquito larval richness (Chao1 estimator) was higher in artificial habitats than in natural habitats. Furthermore, mosquito species composition did not deteriorate with distance between islands. Mosquito abundance by species was associated with microclimatic variables—pH and temperature. We detected co-existence of multiple species at a micro-habitat level with no evidence of interactions like competition or predation. Our study analyzed and identified the most productive larval habitats –discarded plastic container and plastic drums contributing to high larval indices predicting dengue epidemic across the Lakshadweep islands. Our data highlight the need to devise vector control strategies by removal of human-induced plastic pollution (household waste) which is a critical driver of disease risk.

Predicting the Spread of Vector-Borne Diseases in a Warming World

Predicting how climate warming affects vector borne diseases is a key research priority. The prevailing approach uses the basic reproductive number (R0) to predict warming effects. However, R0 is derived under assumptions of stationary thermal environments; using it to predict disease spread in non-stationary environments could lead to erroneous predictions. Here, we develop a trait-based mathematical model that can predict disease spread and prevalence for any vector borne disease under any type of non-stationary environment. We parameterize the model with trait response data for the Malaria vector and pathogen to test the latest IPCC predictions on warmer-than-average winters and hotter-than-average summers. We report three key findings. First, the R0 formulation commonly used to investigate warming effects on disease spread violates the assumptions underlying its derivation as the dominant eigenvalue of a linearized host-vector model. As a result, it overestimates disease spread in cooler environments and underestimates it in warmer environments, proving its predictions to be unreliable even in a constant thermal environment. Second, hotter-than-average summers both narrow the thermal limits for disease prevalence, and reduce prevalence within those limits, to a much greater degree than warmer-than-average winters, highlighting the importance of hot extremes in driving disease burden. Third, while warming reduces infected vector populations through the compounding effects of adult mortality, and infected host populations through the interactive effects of mortality and transmission, uninfected vector populations prove surprisingly robust to warming. This suggests that ecological predictions of warming-induced reductions in disease burden should be tempered by the evolutionary possibility of vector adaptation to both cooler and warmer climates.

Intraspecific Transcriptome Variation and Sex-Biased Expression in Anopheles arabiensis

The magnitude and functional patterns of intraspecific transcriptional variation in the anophelines, including those of sex-biased genes underlying sex-specific traits relevant for malaria transmission, remain understudied. As a result, how changes in expression levels drive adaptation in these species is poorly understood. We sequenced the female, male, and larval transcriptomes of three populations of Anopheles arabiensis from Burkina Faso. One-third of the genes were differentially expressed between populations, often involving insecticide resistance-related genes in a sample type-specific manner, and with the females showing the largest number of differentially expressed genes. At the genomic level, the X chromosome appears depleted of differentially expressed genes compared with the autosomes, chromosomes harboring inversions do not exhibit evidence for enrichment of such genes, and genes that are top contributors to functional enrichment patterns of population differentiation tend to be clustered in the genome. Further, the magnitude of variation for the sex expression ratio across populations did not substantially differ between male- and female-biased genes, except for some populations in which male-limited expressed genes showed more variation than their female counterparts. In fact, female-biased genes exhibited a larger level of interpopulation variation than male-biased genes, both when assayed in males and females. Beyond uncovering the extensive adaptive potential of transcriptional variation in An. Arabiensis, our findings suggest that the evolutionary rate of changes in expression levels on the X chromosome exceeds that on the autosomes, while pointing to female-biased genes as the most variable component of the An. Arabiensis transcriptome.

Sexual selection theory meets disease vector control: Testing harmonic convergence as a "good genes" signal in Aedes aegypti mosquitoes

Background

The mosquito Aedes aegypti is a medically important, globally distributed vector of the viruses that cause dengue, yellow fever, chikungunya, and Zika. Although reproduction and mate choice are key components of vector population dynamics and control, our understanding of the mechanisms of sexual selection in mosquitoes remains poor. In “good genes” models of sexual selection, females use male cues as an indicator of both mate and offspring genetic quality. Recent studies in Ae. aegypti provide evidence that male wingbeats may signal aspects of offspring quality and performance during mate selection in a process known as harmonic convergence. However, the extent to which harmonic convergence may signal overall inherent quality of mates and their offspring remains unknown.

Methodology/Principal findings

To examine this, we measured the relationship between acoustic signaling and a broad panel of parent and offspring fitness traits in two generations of field-derived Ae. aegypti originating from dengue-endemic field sites in Thailand. Our data show that in this population of mosquitoes, harmonic convergence does not signal male fertility, female fecundity, or male flight performance traits, which despite displaying robust variability in both parents and their offspring were only weakly heritable.

Conclusions/Significance

Together, our findings suggest that vector reproductive control programs should treat harmonic convergence as an indicator of some, but not all aspects of inherent quality, and that sexual selection likely affects Ae. aegypti in a trait-, population-, and environment-dependent manner.

Mosquitoes transmit numerous pathogens that disproportionately impact developing countries. The mosquito Aedes aegypti, studied here, transmits viruses that cause neglected tropical diseases such as dengue, yellow fever, chikungunya, and Zika. Disease prevention programs rely heavily upon mosquito vector control. To successfully interrupt disease transmission, several control methods depend upon the ability of laboratory-modified male mosquitoes to successfully mate with wild females to suppress or replace natural populations. However, our understanding of what determines mating success in mosquitoes is far from complete. Our study addresses the question of whether female Ae. aegypti mosquitoes use male acoustic signals to select higher quality mates and improve their offspring’s fitness. We find that acoustic signals do not serve as universal indicators of fitness. Further, the fitness metrics we measured were only weakly heritable, suggesting that females that mate with high quality males do not necessarily produce fitter offspring. Our study provides a nuanced understanding of mate choice, mating acoustic signals, and parent and offspring reproductive fitness in a key disease-transmitting mosquito species. These discoveries improve our grasp of sexual selection in mosquitoes and can be leveraged by the vector control community to improve vitally important disease prevention programs.

The Sterile Insect Technique: Success and Perspectives in the Neotropics

The sterile insect technique (SIT), an environmentally friendly means of control, is currently used against plant, animal, and human pests under the area-wide integrated pest management. It consists in the mass production, sterilization, and release of insects in an affected area where sterile males mate with wild females leading to no reproduction. Here, we review SIT in the Neotropics and focus on particular recent successful cases of eradication of the Mediterranean fruit fly, Ceratitis capitata (Wiedemann), as well as effective programs used against the Mexican fruit fly Anastrepha ludens (Loew), the New World screwworm fly Cochliomyia hominivorax (Coquerel)), and the Cactus moth Cactoblastis cactorum (Berg). We examine when SIT does not work and innovations that have made SIT more efficient and also highlight complimentary techniques that can be used in conjunction. We address potential candidate species that could be controlled through SIT, for example Philornis downsi Dodge & Aitken. Finally, we consider the impact of climate change in the context of the use of the SIT against these pests. Given the recent dramatic decline in insect biodiversity, investing in environmentally friendly means of pest control should be a priority. We conclude that SIT should be promoted in the region, and leadership and political will is needed for continued success of SIT in the Neotropics.

Mosquito Sexual Selection and Reproductive Control Programs

The field of mosquito mating biology has experienced a considerable expansion in the past decade. Recent work has generated many key insights about specific aspects of mating behavior and physiology. Here, we synthesize these findings and classify swarming mosquito systems as polygynous. Male mating success is highly variable in swarms and evidence suggests that it is likely determined by both scramble competition between males and female choice. Incorporating this new understanding will improve both implementation and long-term stability of reproductive control tools.

Dispersal patterns and population genetic structure of Aedes albopictus (Diptera: Culicidae) in three different climatic regions of China

Background

Aedes albopictus is an indigenous primary vector for dengue and Zika viruses in China. Compared with its insecticide resistance, biology and vector competence, little is known about its genetic variation, which corresponds to environmental variations. Thus, the present study examines how Ae. albopictus varies among different climatic regions in China and deciphers its potential dispersal patterns.

Methods

The genetic variation and population structure of 17 Ae. albopictus populations collected from three climatic regions of China were investigated with 11 microsatellite loci and the mitochondrial coxI gene.

Results

Of 44 isolated microsatellite markers, 11 pairs were chosen for genotyping analysis and had an average PIC value of 0.713, representing high polymorphism. The number of alleles was high in each population, with the n_e value increasing from the temperate region (3.876) to the tropical region (4.144). Twenty-five coxI haplotypes were detected, and the highest diversity was observed in the tropical region. The mean Ho value (ca. 0.557) of all the regions was significantly lower than the mean He value (ca. 0.684), with nearly all populations significantly departing from HWE and displaying significant population expansion (p value < 0.05). Two genetically isolated groups and three haplotype clades were evaluated via STRUCTURE and haplotype phylogenetic analyses, and the tropical populations were significantly isolated from those in the other regions. Most genetic variation in Ae. albopictus was detected within populations and individuals at 31.40 and 63.04%, respectively, via the AMOVA test, and a relatively significant positive correlation was observed among only the temperate populations via IBD analysis (R² = 0.6614, p = 0.048). Recent dispersions were observed among different Ae. albopictus populations, and four major migration trends with high gene flow (Nm > 0.4) were reconstructed between the tropical region and the other two regions. Environmental factors, especially temperature and rainfall, may be the leading causes of genetic diversity in different climatic regions.

Conclusions

Continuous dispersion contributes to the genetic communication of Ae. albopictus populations across different climatic regions, and environmental factors, especially temperature and rainfall, may be the leading causes of genetic variation.

Graphical abstract

Genomic differentiation and intercontinental population structure of mosquito vectors Culex pipiens pipiens and Culex pipiens molestus

Understanding the population structure and mechanisms of taxa diversification is important for organisms responsible for the transmission of human diseases. Two vectors of West Nile virus, Culex pipiens pipiens and Cx. p. molestus, exhibit epidemiologically important behavioral and physiological differences, but the whole-genome divergence between them was unexplored. The goal of this study is to better understand the level of genomic differentiation and population structures of Cx. p. pipiens and Cx. p. molestus from different continents. We sequenced and compared the whole genomes of 40 individual mosquitoes from two locations in Eurasia and two in North America. Principal Component, ADMIXTURE, and neighbor joining analyses of the nuclear genomes identified two major intercontinental, monophyletic clusters of Cx. p. pipiens and Cx. p. molestus. The level of genomic differentiation between the subspecies was uniform along chromosomes. The ADMIXTURE analysis determined signatures of admixture in Cx. p. pipens populations but not in Cx. p. molestus populations. Comparison of mitochondrial genomes among the specimens showed a paraphyletic origin of the major haplogroups between the subspecies but a monophyletic structure between the continents. Thus, our study identified that Cx. p. molestus and Cx. p. pipiens represent different evolutionary units with monophyletic origin that have undergone incipient ecological speciation.

Population genetic structure and predominance of cyclical parthenogenesis in the bird cherry-oat aphid Rhopalosiphum padi in England

Genetic diversity is the determinant for pest species' success and vector competence. Understanding the ecological and evolutionary processes that determine the genetic diversity is fundamental to help identify the spatial scale at which pest populations are best managed. In the present study, we present the first comprehensive analysis of the genetic diversity and evolution of Rhopalosiphum padi, a major pest of cereals and a main vector of the barley yellow dwarf virus (BYDV), in England. We have used a genotyping-by-sequencing approach to study whether (a) there is any underlying population genetic structure at a national and regional scale in this pest that can disperse long distances; (b) the populations evolve as a response to environmental change and selective pressures; and (c) the populations comprise anholocyclic lineages. Individual R. padi were collected using the Rothamsted Insect Survey's suction-trap network at several sites across England between 2004 and 2016 as part of the RIS long-term nationwide surveillance. Results identified two genetic clusters in England that mostly corresponded to a North-South division, although gene flow is ongoing between the two subpopulations. These genetic clusters do not correspond to different life cycle types, and cyclical parthenogenesis is predominant in England. Results also show that there is dispersal with gene flow across England, although there is a reduction between the northern and southern sites with the south-western population being the most genetically differentiated. There is no evidence for isolation by distance and other factors such as primary host distribution, uncommon in the south and absent in the south-west, could influence the dispersal patterns. Finally, results also show no evidence for the evolution of the R. padi population, and it is demographically stable despite the ongoing environmental change. These results are discussed in view of their relevance to pest management and the transmission of BYDV.

Chromosome and Genome Divergence between the Cryptic Eurasian Malaria Vector-Species Anopheles messeae and Anopheles daciae

Chromosomal inversions are important drivers of genome evolution. The Eurasian malaria vector Anophelesmesseae has five polymorphic inversions. A cryptic species, An. daciae, has been discriminated from An. messeae based on five fixed nucleotide substitutions in the internal transcribed spacer 2 (ITS2) of ribosomal DNA. However, the inversion polymorphism in An. daciae and the genome divergence between these species remain unexplored. In this study, we sequenced the ITS2 region and analyzed the inversion frequencies of 289 Anopheles larvae specimens collected from three locations in the Moscow region. Five individual genomes for each of the two species were sequenced. We determined that An. messeae and An. daciae differ from each other by the frequency of polymorphic inversions. Inversion X1 was fixed in An. messeae but polymorphic in An. daciae populations. The genome sequence comparison demonstrated genome-wide divergence between the species, especially pronounced on the inversion-rich X chromosome (mean Fst = 0.331). The frequency of polymorphic autosomal inversions was higher in An. messeae than in An. daciae. We conclude that the X chromosome inversions play an important role in the genomic differentiation between the species. Our study determined that An. messeae and An. daciae are closely related species with incomplete reproductive isolation.

Transgenic Aedes aegypti Mosquitoes Transfer Genes into a Natural Population

In an attempt to control the mosquito-borne diseases yellow fever, dengue, chikungunya, and Zika fevers, a strain of transgenically modified Aedes aegypti mosquitoes containing a dominant lethal gene has been developed by a commercial company, Oxitec Ltd. If lethality is complete, releasing this strain should only reduce population size and not affect the genetics of the target populations. Approximately 450 thousand males of this strain were released each week for 27 months in Jacobina, Bahia, Brazil. We genotyped the release strain and the target Jacobina population before releases began for >21,000 single nucleotide polymorphisms (SNPs). Genetic sampling from the target population six, 12, and 27-30 months after releases commenced provides clear evidence that portions of the transgenic strain genome have been incorporated into the target population. Evidently, rare viable hybrid offspring between the release strain and the Jacobina population are sufficiently robust to be able to reproduce in nature. The release strain was developed using a strain originally from Cuba, then outcrossed to a Mexican population. Thus, Jacobina Ae. aegypti are now a mix of three populations. It is unclear how this may affect disease transmission or affect other efforts to control these dangerous vectors. These results highlight the importance of having in place a genetic monitoring program during such releases to detect un-anticipated outcomes.