Adaptive Potential of Hybridization among Malaria Vectors: Introgression at the Immune Locus TEP1 between Anopheles coluzzii and A. gambiae in 'Far-West' Africa

HM, Antoshechkin I, Marshall JM, Akbari OS. Eliminating malaria vectors with precision-guided sterile males

Controlling the principal African malaria vector, the mosquito Anopheles gambiae, is considered essential to curtail malaria transmission. However, existing vector control technologies rely on insecticides, which are becoming increasingly ineffective. Sterile insect technique (SIT) is a powerful suppression approach that has successfully eradicated a number of insect pests, yet the A. gambiae toolkit lacks the requisite technologies for its implementation. SIT relies on iterative mass releases of nonbiting, nondriving, sterile males which seek out and mate with monandrous wild females. Once mated, females are permanently sterilized due to mating-induced refractoriness, which results in population suppression of the subsequent generation. However, sterilization by traditional methods renders males unfit, making the creation of precise genetic sterilization methods imperative. Here, we introduce a vector control technology termed precision-guided sterile insect technique (pgSIT), in A. gambiae for inducible, programmed male sterilization and female elimination for wide-scale use in SIT campaigns. Using a binary CRISPR strategy, we cross separate engineered Cas9 and gRNA strains to disrupt male-fertility and female-essential genes, yielding >99.5% male sterility and >99.9% female lethality in hybrid progeny. We demonstrate that these genetically sterilized males have good longevity, are able to induce sustained population suppression in cage trials, and are predicted to eliminate wild A. gambiae populations using mathematical models, making them ideal candidates for release. This work provides a valuable addition to the malaria genetic biocontrol toolkit, enabling scalable SIT-like confinable, species-specific, and safe suppression in the species.

HM, Antoshechkin I, Marshall JM, Akbari OS. Eliminating Malaria Vectors with Precision Guided Sterile Males

Controlling the principal African malaria vector, the mosquito Anopheles gambiae, is considered essential to curtail malaria transmission. However existing vector control technologies rely on insecticides, which are becoming increasingly ineffective. Sterile insect technique (SIT) is a powerful suppression approach that has successfully eradicated a number of insect pests, yet the A. gambiae toolkit lacks the requisite technologies for its implementation. SIT relies on iterative mass-releases of non-biting, non-driving, sterile males which seek out and mate with monandrous wild females. Once mated, females are permanently sterilized due to mating-induced refractoriness, which results in population suppression of the subsequent generation. However, sterilization by traditional methods renders males unfit, making the creation of precise genetic sterilization methods imperative. Here we develop precision guided Sterile Insect Technique (pgSIT) in the mosquito A. gambiae for inducible, programmed male-sterilization and female-elimination for wide scale use in SIT campaigns. Using a binary CRISPR strategy, we cross separate engineered Cas9 and gRNA strains to disrupt male-fertility and female-essential genes, yielding >99.5% male-sterility and >99.9% female-lethality in hybrid progeny. We demonstrate that these genetically sterilized males have good longevity, are able to induce population suppression in cage trials, and are predicted to eliminate wild A. gambiae populations using mathematical models, making them ideal candidates for release. This work provides a valuable addition to the malaria genetic biocontrol toolkit, for the first time enabling scalable SIT-like confinable suppression in the species.

Gene drive in species complexes: defining target organisms

Engineered gene drives, which bias their own inheritance to increase in frequency in target populations, are being developed to control mosquito malaria vectors. Such mosquitoes can belong to complexes of both vector and nonvector species that can produce fertile interspecific hybrids, making vertical gene drive transfer (VGDT) to sibling species biologically plausible. While VGDT to other vectors could positively impact human health protection goals, VGDT to nonvectors might challenge biodiversity ones. Therefore, environmental risk assessment of gene drive use in species complexes invites more nuanced considerations of target organisms and nontarget organisms than for transgenes not intended to increase in frequency in target populations. Incorporating the concept of target species complexes offers more flexibility when assessing potential impacts from VGDT.

Molecular characterization and genotype distribution of thioester-containing protein 1 gene in Anopheles gambiae mosquitoes in western Kenya

BACKGROUND

Evolutionary pressures lead to the selection of efficient malaria vectors either resistant or susceptible to Plasmodium parasites. These forces may favour the introduction of species genotypes that adapt to new breeding habitats, potentially having an impact on malaria transmission. Thioester-containing protein 1 (TEP1) of Anopheles gambiae complex plays an important role in innate immune defenses against parasites. This study aims to characterize the distribution pattern of TEP1 polymorphisms among populations of An. gambiae sensu lato (s.l.) in western Kenya.

METHODS

Anopheles gambiae adult and larvae were collected using pyrethrum spray catches (PSC) and plastic dippers respectively from Homa Bay, Kakamega, Bungoma, and Kisumu counties between 2017 and 2020. Collected adults and larvae reared to the adult stage were morphologically identified and then identified to sibling species by PCR. TEP1 alleles were determined in 627 anopheles mosquitoes using restriction fragment length polymorphisms-polymerase chain reaction (RFLP-PCR) and to validate the TEP1 genotyping results, a representative sample of the alleles was sequenced.

RESULTS

Two TEP1 alleles (TEP1*S1 and TEP1*R2) and three corresponding genotypes (*S1/S1, *R2/S1, and *R2/R2) were identified. TEP1*S1 and TEP1*R2 with their corresponding genotypes, homozygous *S1/S1 and heterozygous *R2/S1 were widely distributed across all sites with allele frequencies of approximately 80% and 20%, respectively both in Anopheles gambiae and Anopheles arabiensis. There was no significant difference detected among the populations and between the two mosquito species in TEP1 allele frequency and genotype frequency. The overall low levels in population structure (FST = 0.019) across all sites corresponded to an effective migration index (Nm = 12.571) and low Nei's genetic distance values (< 0.500) among the subpopulation. The comparative fixation index values revealed minimal genetic differentiation between species and high levels of gene flow among populations.

CONCLUSION

Genotyping TEP1 has identified two common TEP1 alleles (TEP1*S1 and TEP1*R2) and three corresponding genotypes (*S1/S1, *R2/S1, and *R2/R2) in An. gambiae s.l. The TEP1 allele genetic diversity and population structure are low in western Kenya.

The origin of island populations of the African malaria mosquito, Anopheles coluzzii

Anopheles coluzzii is a major malaria vector throughout its distribution in west-central Africa. Here we present a whole-genome study of 142 specimens from nine countries in continental Africa and three islands in the Gulf of Guinea. This sample set covers a large part of this species' geographic range. Our population genomic analyses included a description of the structure of mainland populations, island populations, and connectivity between them. Three genetic clusters are identified among mainland populations and genetic distances (F_ST) fits an isolation-by-distance model. Genomic analyses are applied to estimate the demographic history and ancestry for each island. Taken together with the unique biogeography and history of human occupation for each island, they present a coherent explanation underlying levels of genetic isolation between mainland and island populations. We discuss the relationship of our findings to the suitability of São Tomé and Príncipe islands as candidate sites for potential field trials of genetic-based malaria control strategies.

Evaluating the probability of CRISPR-based gene drive contaminating another species

The probability D that a given clustered regularly interspaced short palindromic repeats (CRISPR)-based gene drive element contaminates another, nontarget species can be estimated by the following Drive Risk Assessment Quantitative Estimate (DRAQUE) Equation: D = h y b + t r a n s f × e x p r e s s × c u t × f l a n k × i m m u n e × n o n e x t i n c t with hyb = probability of hybridization between the target species and a nontarget species; transf = probability of horizontal transfer of a piece of DNA containing the gene drive cassette from the target species to a nontarget species (with no hybridization); express = probability that the Cas9 and guide RNA genes are expressed; cut = probability that the CRISPR-guide RNA recognizes and cuts at a DNA site in the new host; flank = probability that the gene drive cassette inserts at the cut site; immune = probability that the immune system does not reject Cas9-expressing cells; nonextinct = probability of invasion of the drive within the population. We discuss and estimate each of the seven parameters of the equation, with particular emphasis on possible transfers within insects, and between rodents and humans. We conclude from current data that the probability of a gene drive cassette to contaminate another species is not insignificant. We propose strategies to reduce this risk and call for more work on estimating all the parameters of the formula.

Agricultural insect hybridization and implications for pest management

Biological invasions, the expansion of agricultural frontiers, and climate change favor encounters of divergent lineages of animals and plants, increasing the likelihood of hybridization. However, hybridization of insect species and its consequences for agroecosystems have not received sufficient attention. Gene exchange between distinct and distant genetic pools can improve the survival and reproduction of insect pests, and threaten beneficial insects in disturbed agricultural environments. Hybridization may be the underlying explanation for the recurrent pest outbreaks and control failures in putative hybrid zones, as suspected for bollworm, corn borer, whiteflies, and stink bugs. Reliable predictions of the types of changes that can be expected in pest insect genomes and fitness, and of their impacts on the fate of species and populations remain elusive. Typical steps in pest management, such as insect identification, pest monitoring, and control are likely affected by gene flow and adaptive introgression mediated by hybridization, and we do not have ways to respond to or mitigate the problem. To address the adverse effects of farming intensification and global trade, we must ensure that current integrated pest management programs incorporate up-to-date monitoring and diagnostic tools. The rapid identification of hybrids, quantification of levels of introgression, and in-depth knowledge of what genes have been transferred may help to explain and predict insect population outbreaks and control failures in the future. © 2019 Society of Chemical Industry.

The genomic impact of historical hybridization with massive mitochondrial DNA introgression

BACKGROUND

The extent to which selection determines interspecific patterns of genetic exchange enlightens the role of adaptation in evolution and speciation. Often reported extensive interspecific introgression could be selection-driven, but also result from demographic processes, especially in cases of invasive species replacements, which can promote introgression at their invasion front. Because invasion and selective sweeps similarly mold variation, population genetics evidence for selection can only be gathered in an explicit demographic framework. The Iberian hare, Lepus granatensis, displays in its northern range extensive mitochondrial DNA introgression from L. timidus, an arctic/boreal species that it replaced locally after the last glacial maximum. We use whole-genome sequencing to infer geographic and genomic patterns of nuclear introgression and fit a neutral model of species replacement with hybridization, allowing us to evaluate how selection influenced introgression genome-wide, including for mtDNA.

RESULTS

Although the average nuclear and mtDNA introgression patterns contrast strongly, they fit a single demographic model of post-glacial invasive replacement of timidus by granatensis. Outliers of elevated introgression include several genes related to immunity, spermatogenesis, and mitochondrial metabolism. Introgression is reduced on the X chromosome and in low recombining regions.

CONCLUSIONS

General nuclear and mtDNA patterns of introgression can be explained by purely demographic processes. Hybrid incompatibilities and interplay between selection and recombination locally modulate levels of nuclear introgression. Selection promoted introgression of some genes involved in conflicts, either interspecific (parasites) or possibly cytonuclear. In the latter case, nuclear introgression could mitigate the potential negative effects of alien mtDNA on mitochondrial metabolism and male-specific traits.

Dissecting functional components of reproductive isolation among closely related sympatric species of the Anopheles gambiae complex

Explaining how and why reproductive isolation evolves and determining which forms of reproductive isolation have the largest impact on the process of population divergence are major goals in the study of speciation. By studying recent adaptive radiations in incompletely isolated taxa, it is possible to identify barriers involved at early divergence before other confounding barriers emerge after speciation is complete. Sibling species of the Anopheles gambiae complex offer opportunities to provide insights into speciation mechanisms. Here, we studied patterns of reproductive isolation among three taxa, Anopheles coluzzii, An. gambiae s.s. and Anopheles arabiensis, to compare its strength at different spatial scales, to dissect the relative contribution of pre‐ versus postmating isolation, and to infer the involvement of ecological divergence on hybridization. Because F1 hybrids are viable, fertile and not uncommon, understanding the dynamics of hybridization in this trio of major malaria vectors has important implications for how adaptations arise and spread across the group, and in planning studies of the safety and efficacy of gene drive as a means of malaria control. We first performed a systematic review and meta‐analysis of published surveys reporting on hybrid prevalence, showing strong reproductive isolation at a continental scale despite geographically restricted exceptions. Second, we exploited our own extensive field data sets collected at a regional scale in two contrasting environmental settings, to assess: (i) levels of premating isolation; (ii) spatio/temporal and frequency‐dependent dynamics of hybridization, (iii) relationship between reproductive isolation and ecological divergence and (iv) hybrid viability penalty. Results are in accordance with ecological speciation theory predicting a positive association between the strength of reproductive isolation and degree of ecological divergence, and indicate that postmating isolation does contribute to reproductive isolation among these species. Specifically, only postmating isolation was positively associated with ecological divergence, whereas premating isolation was correlated with phylogenetic distance.

The implementation of long-lasting insecticidal bed nets has differential effects on the genetic structure of the African malaria vectors in the Anopheles gambiae complex in Dielmo, Senegal

BACKGROUND

Mosquitoes belonging to the Anopheles gambiae complex are the main vectors of malaria in sub-Saharan Africa. Among these, An. gambiae, Anopheles coluzzii and Anopheles arabiensis are the most efficient vectors and are largely distributed in sympatric locations. However, these species present ecological and behavioural differences that impact their vectorial capacity and complicate vector-control efforts, mainly based on long-lasting insecticidal bed nets (LLINs) and indoor residual spraying (IRS). In this study, the genetic structure of these three species in a Senegalese village (Dielmo) was investigated using microsatellite data in samples collected in 2006 before implementation of LLINs, in 2008, when they were introduced, and in 2010, 2 years after the use of LLINs.

RESULTS

In this study 611 individuals were included, namely 136 An. coluzzii, 101 An. gambiae, 6 An. coluzzii/An. gambiae hybrids and 368 An. arabiensis. According to the species, the effect of the implementation of LLINs in Dielmo is differentiated. Populations of the sister species An. coluzzii and An. gambiae regularly experienced bottleneck events, but without significant inbreeding. The Fst values suggested in 2006 a breakdown of assortative mating resulting in hybrids, but the introduction of LLINs was followed by a decrease in the number of hybrids. This suggests a decrease in mating success of hybrids, ecological maladaptation, or a lesser probability of mating between species due to a decrease in An. coluzzii population size. By contrast, the introduction of LLINs has favoured the sibling species An. arabiensis. In this study, some spatial and temporal structuration between An. arabiensis populations were detected, especially in 2008, and the higher genetic diversity observed could result from a diversifying selection.

CONCLUSIONS

This work demonstrates the complexity of the malaria context and shows the need to study the genetic structure of Anopheles populations to evaluate the effectiveness of vector-control tools and successful management of malaria vector control.

Massive introgression drives species radiation at the range limit of Anopheles gambiae

Impacts of introgressive hybridisation may range from genomic erosion and species collapse to rapid adaptation and speciation but opportunities to study these dynamics are rare. We investigated the extent, causes and consequences of a hybrid zone between Anopheles coluzzii and Anopheles gambiae in Guinea-Bissau, where high hybridisation rates appear to be stable at least since the 1990s. Anopheles gambiae was genetically partitioned into inland and coastal subpopulations, separated by a central region dominated by A. coluzzii. Surprisingly, whole genome sequencing revealed that the coastal region harbours a hybrid form characterised by an A. gambiae-like sex chromosome and massive introgression of A. coluzzii autosomal alleles. Local selection on chromosomal inversions may play a role in this process, suggesting potential for spatiotemporal stability of the coastal hybrid form and providing resilience against introgression of medically-important loci and traits, found to be more prevalent in inland A. gambiae.

Immersive virtual reality in computational chemistry: Applications to the analysis of QM and MM data

The role of Virtual Reality (VR) tools in molecular sciences is analyzed in this contribution through the presentation of the Caffeine software to the quantum chemistry community. Caffeine, developed at Scuola Normale Superiore, is specifically tailored for molecular representation and data visualization with VR systems, such as VR theaters and helmets. Usefulness and advantages that can be gained by exploiting VR are here reported, considering few examples specifically selected to illustrate different level of theory and molecular representation.

The last bastion? X chromosome genotyping of Anopheles gambiae species pair males from a hybrid zone reveals complex recombination within the major candidate 'genomic island of speciation'

Speciation with gene flow may be aided by reduced recombination helping to build linkage between genes involved in the early stages of reproductive isolation. Reduced recombination on chromosome X has been implicated in speciation within the Anopheles gambiae complex, species of which represent the major Afrotropical malaria vectors. The most recently diverged, morphologically indistinguishable, species pair, A. gambiae and Anopheles coluzzii, ubiquitously displays a 'genomic island of divergence' spanning over 4 Mb from chromosome X centromere, which represents a particularly promising candidate region for reproductive isolation genes, in addition to containing the diagnostic markers used to distinguish the species. Very low recombination makes the island intractable for experimental recombination studies, but an extreme hybrid zone in Guinea Bissau offers the opportunity for natural investigation of X-island recombination. SNP analysis of chromosome X hemizygous males revealed: (i) strong divergence in the X-island despite a lack of autosomal divergence; (ii) individuals with multiple-recombinant genotypes, including likely double crossovers and localized gene conversion; (iii) recombination-driven discontinuity both within and between the molecular species markers, suggesting that the utility of the diagnostics is undermined under high hybridization. The largely, but incompletely protected nature of the X centromeric genomic island is consistent with a primary candidate area for accumulation of adaptive variants driving speciation with gene flow, while permitting some selective shuffling and removal of genetic variation.