Adaptive introgression in an African malaria mosquito coincident with the increased usage of insecticide-treated bed nets

Complex patterns of gene flow and convergence in the evolutionary history of the spiral-horned antelopes (Tragelaphini)

The Tragelaphini, also known as spiral-horned antelope, is a phenotypically diverse mammalian tribe comprising a single genus, Tragelaphus. The evolutionary history of this tribe has attracted the attention of taxonomists and molecular geneticists for decades because its diversity is characterised by conflicts between morphological and molecular data as well as between mitochondrial, nuclear and chromosomal DNA. These inconsistencies point to a complex history of ecological diversification, coupled by either phenotypic convergence or introgression. Therefore, to unravel the phylogenetic relationships among spiral-horned antelopes, and to further investigate the role of divergence and gene flow in trait evolution, we sequenced genomes for all nine accepted species of the genus Tragelaphus, including a genome each for the highly divergent bushbuck lineages (T. s. scriptus and T. s. sylvaticus). We successfully reconstructed the Tragelaphus species tree, providing genome-level support for the early Pliocene divergence and monophyly of the nyala (T. angasii) and lesser kudu (T. imberbis), the monophyly of the two eland species (T. oryx and T. derbianus) and, importantly, the monophyly of kéwel (T. s. scriptus) and imbabala (T. s. sylvaticus) bushbuck. We found strong evidence for gene flow in at least four of eight nodes on the species tree. Among the six phenotypic traits assessed here, only habitat type mapped onto the species tree without homoplasy, showing that trait evolution was the result of complex patterns of divergence, introgression and convergent evolution.

Genetic admixture drives climate adaptation in the bank vole

Genetic admixture introduces new variants at relatively high frequencies, potentially aiding rapid responses to environmental changes. Here, we evaluate its role in adaptive variation related to climatic conditions in bank voles (Clethrionomys glareolus) in Britain, using whole-genome data. Our results reveal loci showing excess ancestry from one of the two postglacial colonist populations inconsistent with overall admixture patterns. Notably, loci associated with climate adaptation exhibit disproportionate amounts of excess ancestry, highlighting the impact of admixture between colonist populations on local adaptation. The results suggest strong and localized selection on climate-adaptive loci, as indicated by steep clines and/or shifted cline centres, during population replacement. A subset, including a haemoglobin gene, is associated with oxidative stress responses, underscoring a role of oxidative stress in local adaptation. Our study highlights the important contribution of admixture during secondary contact between populations from distinct climatic refugia enriching adaptive diversity. Understanding these dynamics is crucial for predicting future adaptive capacity to anthropogenic climate change.

The Anopheles coluzzii range extends into Kenya: detection, insecticide resistance profiles and population genetic structure in relation to conspecific populations in West and Central Africa

BACKGROUND

Anopheles coluzzii is a primary vector of malaria found in West and Central Africa, but its presence has hitherto never been documented in Kenya. A thorough understanding of vector bionomics is important as it enables the implementation of targeted and effective vector control interventions. Malaria vector surveillance efforts in the country have tended to focus on historically known primary vectors. The current study sought to determine the taxonomic status of samples collected from five different malaria epidemiological zones in Kenya as well as describe the population genetic structure and insecticide resistance profiles in relation to other An. coluzzii populations.

METHODS

Mosquitoes were sampled as larvae from Busia, Kwale, Turkana, Kirinyaga and Kiambu counties, representing the range of malaria endemicities in Kenya, in 2019 and 2021 and emergent adults analysed using Whole Genome Sequencing (WGS) data processed in accordance with the Anopheles gambiae 1000 Genomes Project phase 3. Where available, historical samples from the same sites were included for WGS. Comparisons were made with An. coluzzii cohorts from West and Central Africa.

RESULTS

This study reports the detection of An. coluzzii for the first time in Kenya. The species was detected in Turkana County across all three time points from which samples were analyzed and its presence confirmed through taxonomic analysis. Additionally, there was a lack of strong population genetic differentiation between An. coluzzii from Kenya and those from the more northerly regions of West and Central Africa, suggesting they represent a connected extension to the known species range. Mutations associated with target-site resistance to DDT and pyrethroids and metabolic resistance to DDT were found at high frequencies up to 64%. The profile and frequencies of the variants observed were similar to An. coluzzii from West and Central Africa but the ace-1 mutation linked to organophosphate and carbamate resistance present in An. coluzzii from coastal West Africa was absent in Kenya.

CONCLUSIONS

These findings emphasize the need for the incorporation of genomics in comprehensive and routine vector surveillance to inform on the range of malaria vector species, and their insecticide resistance status to inform the choice of effective vector control approaches.

An optical system to detect, surveil, and kill flying insect vectors of human and crop pathogens

Sustainable and effective means to control flying insect vectors are critically needed, especially with widespread insecticide resistance and global climate change. Understanding and controlling vectors requires accurate information about their movement and activity, which is often lacking. The Photonic Fence (PF) is an optical system that uses machine vision, infrared light, and lasers to identify, track, and interdict vectors in flight. The PF examines an insect's outline, flight speed, and other flight parameters and if these match those of a targeted vector species, then a low-power, retina-safe laser kills it. We report on proof-of-concept tests of a large, field-sized PF (30 mL × 3 mH) conducted with Aedes aegypti, a mosquito that transmits dangerous arboviruses, and Diaphorina citri, a psyllid which transmits the fatal huanglongbing disease of citrus. In tests with the laser engaged, < 1% and 3% of A. aegypti and D. citri, respectfully, were recovered versus a 38% and 19% recovery when the lacer was silenced. The PF tracked, but did not intercept the orchid bee, Euglossa dilemma. The system effectively intercepted flying vectors, but not bees, at a distance of 30 m, heralding the use of photonic energy, rather than chemicals, to control flying vectors.

Yosemite toad (Anaxyrus canorus) transcriptome reveals interplay between speciation genes and adaptive introgression

Genomes are heterogeneous during the early stages of speciation, with small 'islands' of DNA appearing to reflect strong adaptive differences, surrounded by vast seas of relative homogeneity. As species diverge, secondary contact zones between them can act as an interface and selectively filter through advantageous alleles of hybrid origin. Such introgression is another important adaptive process, one that allows beneficial mosaics of recombinant DNA ('rivers') to flow from one species into another. Although genomic islands of divergence appear to be associated with reproductive isolation, and genomic rivers form by adaptive introgression, it is unknown whether islands and rivers tend to be the same or different loci. We examined three replicate secondary contact zones for the Yosemite toad (Anaxyrus canorus) using two genomic data sets and a morphometric data set to answer the questions: (1) How predictably different are islands and rivers, both in terms of genomic location and gene function? (2) Are the adaptive genetic trait loci underlying tadpole growth and development reliably islands, rivers or neither? We found that island and river loci have significant overlap within a contact zone, suggesting that some loci are first islands, and later are predictably converted into rivers. However, gene ontology enrichment analysis showed strong overlap in gene function unique to all island loci, suggesting predictability in overall gene pathways for islands. Genome-wide association study outliers for tadpole development included LPIN3, a lipid metabolism gene potentially involved in climate change adaptation, that is island-like for all three contact zones, but also appears to be introgressing (as a river) across one zone. Taken together, our results suggest that adaptive divergence and introgression may be more complementary forces than currently appreciated.

Accurate identification of Helicoverpa armigera-Helicoverpa zea hybrids using genome admixture analysis: implications for genomic surveillance

Helicoverpa armigera, the cotton bollworm moth, is one of the world's most important crop pests, and is spreading throughout the New World from its original range in the Old World. In Brazil, invasive H. armigera has been reported to hybridize with local populations of Helicoverpa zea. The correct identification of H. armigera-H. zea hybrids is important in understanding the origin, spread and future outlook for New World regions that are affected by outbreaks, given that hybridization can potentially facilitate H. zea pesticide resistance and host plant range via introgression of H. armigera genes. Here, we present a genome admixture analysis of high quality genome sequences generated from two H. armigera-H. zea F1 hybrids generated in two different labs. Our admixture pipeline predicts 48.8% and 48.9% H. armigera for the two F1 hybrids, confirming its accuracy. Genome sequences from five H. zea and one H. armigera that were generated as part of the study show no evidence of hybridization. Interestingly, we show that four H. zea genomes generated from a previous study are predicted to possess a proportion of H. armigera genetic material. Using unsupervised clustering to identify non-hybridized H. armigera and H. zea genomes, 8511 ancestry informative markers (AIMs) were identified. Their relative frequencies are consistent with a minor H. armigera component in the four genomes, however its origin remains to be established. We show that the size and quality of genomic reference datasets are critical for accurate hybridization prediction. Consequently, we discuss potential pitfalls in genome admixture analysis of H. armigera-H. zea hybrids, and suggest measures that will improve such analyses.

The Anopheles coluzzii range extends into Kenya: Detection, insecticide resistance profiles and population genetic structure in relation to conspecific populations in West and Central Africa

Background

Anopheles coluzzii is a primary vector of malaria found in West and Central Africa, but its presence has hitherto never been documented in Kenya. A thorough understanding of vector bionomics is important as it enables the implementation of targeted and effective vector control interventions. Malaria vector surveillance efforts in the country have tended to focus on historically known primary vectors. In the current study, we sought to determine the taxonomic status of samples collected from five different malaria epidemiological zones in Kenya as well asdescribe the population genetic structure and insecticide resistance profiles in relation to other An. coluzzi populations.

Methods

Mosquitoes were sampled as larvae from Busia, Kwale, Turkana, Kirinyaga and Kiambu counties, representing the range of malaria endemicities in Kenya, in 2019 and 2021 and emergent adults analysed using Whole Genome Sequencing data processed in accordance with the Anopheles gambiae 1000 Genomes Project phase 3. Where available, historical samples from the same sites were included for WGS.

Results

This study reports the detection of Anopheles coluzzii for the first time in Kenya. The species was detected in Turkana County across all three time points sampled and its presence confirmed through taxonomic analysis. Additionally, we found a lack of strong population genetic differentiation between An. coluzzii from Kenya and those from the more northerly regions of West and Central Africa, suggesting they represent a connected extension to the known species range. Mutations associated with target-site resistance to DDT and pyrethroids and metabolic resistance to DDT were found at high frequencies of ~60%. The profile and frequencies of the variants observed were similar to An. coluzzii from West and Central Africa but the ace-1 mutation linked to organophosphate and carbamate resistance present in An. coluzzii from coastal West Africa was absent in Kenya.

Conclusions

These findings emphasise the need for the incorporation of genomics in comprehensive and routine vector surveillance to inform on the range of malaria vector species, and their insecticide resistance status to inform the choice of effective vector control approaches.

IntroUNET: Identifying introgressed alleles via semantic segmentation

A growing body of evidence suggests that gene flow between closely related species is a widespread phenomenon. Alleles that introgress from one species into a close relative are typically neutral or deleterious, but sometimes confer a significant fitness advantage. Given the potential relevance to speciation and adaptation, numerous methods have therefore been devised to identify regions of the genome that have experienced introgression. Recently, supervised machine learning approaches have been shown to be highly effective for detecting introgression. One especially promising approach is to treat population genetic inference as an image classification problem, and feed an image representation of a population genetic alignment as input to a deep neural network that distinguishes among evolutionary models (i.e. introgression or no introgression). However, if we wish to investigate the full extent and fitness effects of introgression, merely identifying genomic regions in a population genetic alignment that harbor introgressed loci is insufficient-ideally we would be able to infer precisely which individuals have introgressed material and at which positions in the genome. Here we adapt a deep learning algorithm for semantic segmentation, the task of correctly identifying the type of object to which each individual pixel in an image belongs, to the task of identifying introgressed alleles. Our trained neural network is thus able to infer, for each individual in a two-population alignment, which of those individual's alleles were introgressed from the other population. We use simulated data to show that this approach is highly accurate, and that it can be readily extended to identify alleles that are introgressed from an unsampled "ghost" population, performing comparably to a supervised learning method tailored specifically to that task. Finally, we apply this method to data from Drosophila, showing that it is able to accurately recover introgressed haplotypes from real data. This analysis reveals that introgressed alleles are typically confined to lower frequencies within genic regions, suggestive of purifying selection, but are found at much higher frequencies in a region previously shown to be affected by adaptive introgression. Our method's success in recovering introgressed haplotypes in challenging real-world scenarios underscores the utility of deep learning approaches for making richer evolutionary inferences from genomic data.

IntroUNET: identifying introgressed alleles via semantic segmentation

A growing body of evidence suggests that gene flow between closely related species is a widespread phenomenon. Alleles that introgress from one species into a close relative are typically neutral or deleterious, but sometimes confer a significant fitness advantage. Given the potential relevance to speciation and adaptation, numerous methods have therefore been devised to identify regions of the genome that have experienced introgression. Recently, supervised machine learning approaches have been shown to be highly effective for detecting introgression. One especially promising approach is to treat population genetic inference as an image classification problem, and feed an image representation of a population genetic alignment as input to a deep neural network that distinguishes among evolutionary models (i.e. introgression or no introgression). However, if we wish to investigate the full extent and fitness effects of introgression, merely identifying genomic regions in a population genetic alignment that harbor introgressed loci is insufficient-ideally we would be able to infer precisely which individuals have introgressed material and at which positions in the genome. Here we adapt a deep learning algorithm for semantic segmentation, the task of correctly identifying the type of object to which each individual pixel in an image belongs, to the task of identifying introgressed alleles. Our trained neural network is thus able to infer, for each individual in a two-population alignment, which of those individual's alleles were introgressed from the other population. We use simulated data to show that this approach is highly accurate, and that it can be readily extended to identify alleles that are introgressed from an unsampled "ghost" population, performing comparably to a supervised learning method tailored specifically to that task. Finally, we apply this method to data from Drosophila, showing that it is able to accurately recover introgressed haplotypes from real data. This analysis reveals that introgressed alleles are typically confined to lower frequencies within genic regions, suggestive of purifying selection, but are found at much higher frequencies in a region previously shown to be affected by adaptive introgression. Our method's success in recovering introgressed haplotypes in challenging real-world scenarios underscores the utility of deep learning approaches for making richer evolutionary inferences from genomic data.

Analysis of geometric morphometrics and molecular phylogeny for Anopheles species in the Republic of Korea

Human malaria, transmitted by Anopheles mosquitoes, is the most predominant mosquito-borne disease that is responsible for hundreds of thousands of deaths worldwide each year. In the Republic of Korea (ROK), there are currently several hundred malaria cases annually, mostly near the demilitarized zone (DMZ). Eight species of Anopheles mosquitoes are currently known to be present in the ROK. Similar to other major malaria vectors in Africa and India, it is very challenging to morphologically differentiate Anopheles mosquitoes in the ROK due to their extremely similar morphology. In this study, wing geometric morphometrics (WGM) were used to differentiate the eight Anopheles species collected at six locations near the DMZ, Seoul and Pyeongtaek from April-October 2021. Phylogenetic analysis was also performed using cytochrome c oxidase subunit 1 (COI), internal transcribed spacer 2 (ITS2), and tyrosine hydroxylase (TH) genes for comparison with WGM analysis and to infer evolutionary relationships. The results of cross-validation (overall accuracy = 74.8%) demonstrated that species identification using WGM alone was not possible with a high accuracy for all eight species. While phylogenetic analyses based on the COI region could not clearly distinguish some species, the analysis based on ITS2 and TH was more useful for resolving the phylogenetic correlation of the eight species. Our results may improve Anopheles species identification strategies for effective identification and control of malaria vectors in the ROK.

The ecological importance of hybridization

Hybridization as an evolutionary process has been studied in depth over the past few decades. Research has focused on its role in shaping reproductive barriers, its adaptive value, and its genomic consequences. In contrast, our knowledge of ecological dimensions of hybridization is still in its infancy, despite hybridization being an inherently ecological interaction. Using examples from various organisms, we show that hybridization can affect and be affected by non-reproductive interactions, including predation, competition, parasitism, mutualism, commensalism, and organism-environment interactions, with significant implications for community structure and ecosystem functioning. However, since these dimensions of hybridization have mostly been revealed from studies designed to decipher other evolutionary processes, we argue that much of the eco-evolutionary importance of hybridization is yet to be discovered.

Conceptual risk assessment of mosquito population modification gene-drive systems to control malaria transmission: preliminary hazards list workshops

The field-testing and eventual adoption of genetically-engineered mosquitoes (GEMs) to control vector-borne pathogen transmission will require them meeting safety criteria specified by regulatory authorities in regions where the technology is being considered for use and other locales that might be impacted. Preliminary risk considerations by researchers and developers may be useful for planning the baseline data collection and field research used to address the anticipated safety concerns. Part of this process is to identify potential hazards (defined as the inherent ability of an entity to cause harm) and their harms, and then chart the pathways to harm and evaluate their probability as part of a risk assessment. The University of California Malaria Initiative (UCMI) participated in a series of workshops held to identify potential hazards specific to mosquito population modification strains carrying gene-drive systems coupled to anti-parasite effector genes and their use in a hypothetical island field trial. The hazards identified were placed within the broader context of previous efforts discussed in the scientific literature. Five risk areas were considered i) pathogens, infections and diseases, and the impacts of GEMs on human and animal health, ii) invasiveness and persistence of GEMs, and interactions of GEMs with target organisms, iii) interactions of GEMs with non-target organisms including horizontal gene transfer, iv) impacts of techniques used for the management of GEMs and v) evolutionary and stability considerations. A preliminary hazards list (PHL) was developed and is made available here. This PHL is useful for internal project risk evaluation and is available to regulators at prospective field sites. UCMI project scientists affirm that the subsequent processes associated with the comprehensive risk assessment for the application of this technology should be driven by the stakeholders at the proposed field site and areas that could be affected by this intervention strategy.

Mitochondrial introgression and mito-nuclear discordance obscured the closely related species boundaries in Cletus Stål from China (Heteroptera: Coreidae)

Accurate taxonomy and delimitation are of great importance for pest control strategies and management programs. Here, we focus on Cletus (Insecta: Hemiptera: Coreidae), which includes many crop pests. The species boundaries still conflict and only cytochrome c oxidase subunit I (COI) barcoding has been previously used for molecular studies. We generated new mitochondrial genome and nuclear genome-wide SNPs to explore the species boundaries of 46 Cletus samples from China using multiple species delimitation approaches. All results recovered a monophyly with high support, except for two closely related species in clade I - C. punctiger and C. graminis. Mitochondrial data demonstrated admixture in clade I, while genome-wide SNPs unambiguously identified two separate species, which were confirmed by morphological classification. Inconsistent nuclear and mitochondrial data indicated mito-nuclear discordance. Mitochondrial introgression is the most likely explanation, and more extensive sampling and more comprehensive data are needed to ascertain a pattern. Accurate species delimitation will shed light on species status; thus, an accurate taxonomy is of particular concern, as there is a pressing need to implement precise control of agricultural pests and to perform further research on diversification.

The fitness of an introgressing haplotype changes over the course of divergence and depends on its size and genomic location

The genomic era has made clear that introgression, or the movement of genetic material between species, is a common feature of evolution. Examples of both adaptive and deleterious introgression exist in a variety of systems. What is unclear is how the fitness of an introgressing haplotype changes as species diverge or as the size of the introgressing haplotype changes. In a simple model, we show that introgression may more easily occur into parts of the genome which have not diverged heavily from a common ancestor. The key insight is that alleles from a shared genetic background are likely to have positive epistatic interactions, increasing the fitness of a larger introgressing block. In regions of the genome where few existing substitutions are disrupted, this positive epistasis can be larger than incompatibilities with the recipient genome. Further, we show that early in the process of divergence, introgression of large haplotypes can be favored more than introgression of individual alleles. This model is consistent with observations of a positive relationship between recombination rate and introgression frequency across the genome; however, it generates several novel predictions. First, the model suggests that the relationship between recombination rate and introgression may not exist, or may be negative, in recently diverged species pairs. Furthermore, the model suggests that introgression that replaces existing derived variation will be more deleterious than introgression at sites carrying ancestral variants. These predictions are tested in an example of introgression in Drosophila melanogaster, with some support for both. Finally, the model provides a potential alternative explanation to asymmetry in the direction of introgression, with expectations of higher introgression from rapidly diverged populations into slowly evolving ones.

Evidence of hybridization between genetically distinct Baltic cod stocks during peak population abundance(s)

Range expansions can lead to increased contact of divergent populations, thus increasing the potential of hybridization events. Whether viable hybrids are produced will most likely depend on the level of genomic divergence and associated genomic incompatibilities between the different entities as well as environmental conditions. By taking advantage of historical Baltic cod (Gadus morhua) otolith samples combined with genotyping and whole genome sequencing, we here investigate the genetic impact of the increased spawning stock biomass of the eastern Baltic cod stock in the mid 1980s. The eastern Baltic cod is genetically highly differentiated from the adjacent western Baltic cod and locally adapted to the brackish environmental conditions in the deeper Eastern basins of the Baltic Sea unsuitable for its marine counterparts. Our genotyping results show an increased proportion of eastern Baltic cod in western Baltic areas (Mecklenburg Bay and Arkona Basin)-indicative of a range expansion westwards-during the peak population abundance in the 1980s. Additionally, we detect high frequencies of potential hybrids (including F1, F2 and backcrosses), verified by whole genome sequencing data for a subset of individuals. Analysis of mitochondrial genomes further indicates directional gene flow from eastern Baltic cod males to western Baltic cod females. Our findings unravel that increased overlap in distribution can promote hybridization between highly divergent populations and that the hybrids can be viable and survive under specific and favourable environmental conditions. However, the observed hybridization had seemingly no long-lasting impact on the continuous separation and genetic differentiation between the unique Baltic cod stocks.

Localizing Post-Admixture Adaptive Variants with Object Detection on Ancestry-Painted Chromosomes

Gene flow between previously differentiated populations during the founding of an admixed or hybrid population has the potential to introduce adaptive alleles into the new population. If the adaptive allele is common in one source population, but not the other, then as the adaptive allele rises in frequency in the admixed population, genetic ancestry from the source containing the adaptive allele will increase nearby as well. Patterns of genetic ancestry have therefore been used to identify post-admixture positive selection in humans and other animals, including examples in immunity, metabolism, and animal coloration. A common method identifies regions of the genome that have local ancestry "outliers" compared with the distribution across the rest of the genome, considering each locus independently. However, we lack theoretical models for expected distributions of ancestry under various demographic scenarios, resulting in potential false positives and false negatives. Further, ancestry patterns between distant sites are often not independent. As a result, current methods tend to infer wide genomic regions containing many genes as under selection, limiting biological interpretation. Instead, we develop a deep learning object detection method applied to images generated from local ancestry-painted genomes. This approach preserves information from the surrounding genomic context and avoids potential pitfalls of user-defined summary statistics. We find the method is robust to a variety of demographic misspecifications using simulated data. Applied to human genotype data from Cabo Verde, we localize a known adaptive locus to a single narrow region compared with multiple or long windows obtained using two other ancestry-based methods.

Introgression of a synthetic sex ratio distortion transgene into different genetic backgrounds of Anopheles coluzzii

The development of genetically modified mosquitoes (GMM) and their subsequent field release offers innovative approaches for vector control of malaria. A non-gene drive self-limiting male-bias Ag(PMB)1 strain has been developed in a 47-year-old laboratory G3 strain of Anopheles gambiae s.l. When Ag(PMB)1 males are crossed to wild-type females, expression of the endonuclease I-PpoI during spermatogenesis causes the meiotic cleavage of the X chromosome in sperm cells, leading to fertile offspring with a 95% male bias. However, World Health Organization states that the functionality of the transgene could differ when inserted in different genetic backgrounds of Anopheles coluzzii which is currently a predominant species in several West-African countries and thus a likely recipient for a potential release of self-limiting GMMs. In this study, we introgressed the transgene from the donor Ag(PMB)1 by six serial backcrosses into two recipient colonies of An. coluzzii that had been isolated in Mali and Burkina Faso. Scans of informative Single Nucleotide Polymorphism (SNP) markers and whole-genome sequencing analysis revealed a nearly complete introgression of chromosomes 3 and X, but a remarkable genomic divergence in a large region of chromosome 2 between the later backcrossed (BC6) transgenic offspring and the recipient paternal strains. These findings suggested to extend the backcrossing breeding strategy beyond BC6 generation and increasing the introgression efficiency of critical regions that have ecological and epidemiological implications through the targeted selection of specific markers. Disregarding differential introgression efficiency, we concluded that the phenotype of the sex ratio distorter is stable in the BC6 introgressed An. coluzzii strains.

Rapid adaptation in a fast-changing world: Emerging insights from insect genomics

Many researchers have questioned the ability of biota to adapt to rapid anthropogenic environmental shifts. Here, we synthesize emerging genomic evidence for rapid insect evolution in response to human pressure. These new data reveal diverse genomic mechanisms (single locus, polygenic, structural shifts; introgression) underpinning rapid adaptive responses to a variety of anthropogenic selective pressures. While the effects of some human impacts (e.g. pollution; pesticides) have been previously documented, here we highlight startling new evidence for rapid evolutionary responses to additional anthropogenic processes such as deforestation. These recent findings indicate that diverse insect assemblages can indeed respond dynamically to major anthropogenic evolutionary challenges. Our synthesis also emphasizes the critical roles of genomic architecture, standing variation and gene flow in maintaining future adaptive potential. Broadly, it is clear that genomic approaches are essential for predicting, monitoring and responding to ongoing anthropogenic biodiversity shifts in a fast-changing world.

Hybridization between closely related songbirds is related to human habitat disturbance

Human habitat disturbances can promote hybridization between closely related, but typically reproductively isolated, species. We explored whether human habitat disturbances are related to hybridization between two closely related songbirds, black-capped and mountain chickadees, using both genomic and citizen science data sets. First, we genotyped 409 individuals from across both species' ranges using reduced-representation genome sequencing and compared measures of genetic admixture to a composite measure of human landscape disturbance. Then, using eBird observations, we compared human landscape disturbance values for sites where phenotypically diagnosed hybrids were observed to locations where either parental species was observed to determine whether hybrid chickadees are reported in more disturbed areas. We found that hybridization between black-capped and mountain chickadees positively correlates with human habitat disturbances. From genomic data, we found that (1) hybrid index (HI) significantly increased with habitat disturbance, (2) more hybrids were sampled in disturbed habitats, (3) mean HIs were higher in disturbed habitats versus wild habitats, and (4) hybrids were detected in habitats with significantly higher disturbance values than parentals. Using eBird data, we found that both hybrid and black-capped chickadees were significantly more disturbance-associated than mountain chickadees. Surprisingly, we found that nearly every black-capped chickadee we sampled contained some proportion of hybrid ancestry, while we detected very few mountain chickadee backcrosses. Our results highlight that hybridization between black-capped and mountain chickadees is widespread, but initial hybridization is rare (few F1s were detected). We conclude that human habitat disturbances can erode pre-zygotic reproductive barriers between chickadees and that post-zygotic isolation is incomplete. Understanding what becomes of recently hybridizing species following large-scale habitat disturbances is a new, but pressing, consideration for successfully preserving genetic biodiversity in a rapidly changing world.

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.

Revisiting the Baltimore-Bullock's Oriole hybrid zone reveals changing plumage colour in Bullock's Orioles

Hybrid zones are dynamic areas where populations of two or more interbreeding species may change through an influx of novel genetic material resulting from hybridization or selection on standing genetic variation. Documenting changes in populations through time, however, is challenging because repeated samples are often missing or because long-term storage can affect trait morphologies, especially colour traits that may fade through time. We document a change in carotenoid-based orange breast feathers of Bullock's Orioles (Icterus bullockii) from the Great Plains hybrid zone, USA. Contemporary Bullock's Orioles are more orange than historic individuals from the same location sampled approximately 60 years ago. Spectrophotometry revealed that contemporary Bullock's Orioles resemble orange colour profiles of Baltimore Orioles (I. galbula), the species with which they hybridize. Fading or changes in diet hypotheses do not appear to explain the shift in colour we report for Bullock's Orioles. We propose that these changes in colour are facilitated through introgression with Baltimore Orioles, and favoured by females that choose brighter, more orange males. Our study highlights the long memory of natural history collections and how they offer new insights to the dynamic roll of hybrid zones in trait evolution between interacting species.

A need for standardized reporting of introgression: Insights from studies across eukaryotes

With the rise of affordable next-generation sequencing technology, introgression-or the exchange of genetic materials between taxa-has become widely perceived to be a ubiquitous phenomenon in nature. Although this claim is supported by several keystone studies, no thorough assessment of the frequency of introgression across eukaryotes in nature has been performed to date. In this manuscript, we aim to address this knowledge gap by examining patterns of introgression across eukaryotes. We collated a single statistic, Patterson's D, which can be used as a test for introgression across 123 studies to further assess how taxonomic group, divergence time, and sequencing technology influence reports of introgression. Overall, introgression has mostly been measured in plants and vertebrates, with less attention given to the rest of the Eukaryotes. We find that the most frequently used metrics to detect introgression are difficult to compare across studies and even more so across biological systems due to differences in study effort, reporting standards, and methodology. Nonetheless, our analyses reveal several intriguing patterns, including the observation that differences in sequencing technologies may bias values of Patterson's D and that introgression may differ throughout the course of the speciation process. Together, these results suggest the need for a unified approach to quantifying introgression in natural communities and highlight important areas of future research that can be better assessed once this unified approach is met.

Perfect association between spatial swarm segregation and the X-chromosome speciation island in hybridizing Anopheles coluzzii and Anopheles gambiae populations

The sibling species An. coluzzii and An. gambiae s.s. are major malaria vectors thought to be undergoing sympatric speciation with gene flow. In the absence of intrinsic post-zygotic isolation between the two taxa, speciation is thought possible through the association of assortative mating and genomic regions protected from gene flow by recombination suppression. Such genomic islands of speciation have been described in pericentromeric regions of the X, 2L and 3L chromosomes. Spatial swarm segregation plays a major role in assortative mating between sympatric populations of the two species and, given their importance for speciation, genes responsible for such pre-mating reproductive barriers are expected to be protected within divergence islands. In this study 2063 male and 266 female An. coluzzii and An. gambiae s.s. individuals from natural swarms in Burkina Faso, West Africa were sampled. These were genotyped at 16 speciation island SNPs, and characterized as non-hybrid individuals, F₁ hybrids or recombinant F_1+n backcrossed individuals. Their genotypes at each speciation island were associated with their participation in An. coluzzii and An. gambiae-like swarms. Despite extensive introgression between the two species, the X-island genotype of non-hybrid individuals (37.6%), F₁ hybrids (0.1%) and F_1+n recombinants (62.3%) of either sex perfectly associated to each swarm type. Associations between swarm type and the 3L and 2L speciation islands were weakened or broken down by introgression. The functional demonstration of a close association between spatial segregation behaviour and the X speciation island lends further support to sympatric speciation models facilitated by pericentric recombination suppression in this important species complex.

Genome-wide Association Study Reveals New Loci Associated With Pyrethroid Resistance in Aedes aegypti

Genome-wide association studies (GWAS) use genetic polymorphism across the genomes of individuals with distinct characteristics to identify genotype-phenotype associations. In mosquitoes, complex traits such as vector competence and insecticide resistance could benefit from GWAS. We used the Aedes aegypti 50k SNP chip to genotype populations with different levels of pyrethroid resistance from Northern Brazil. Pyrethroids are widely used worldwide to control mosquitoes and agricultural pests, and their intensive use led to the selection of resistance phenotypes in many insects including mosquitoes. For Ae. aegypti, resistance phenotypes are mainly associated with several mutations in the voltage-gated sodium channel, known as knockdown resistance (kdr). We phenotyped those populations with the WHO insecticide bioassay using deltamethrin impregnated papers, genotyped the kdr alleles using qPCR, and determined allele frequencies across the genome using the SNP chip. We identified single-nucleotide polymorphisms (SNPs) directly associated with resistance and one epistatic SNP pair. We also observed that the novel SNPs correlated with the known kdr genotypes, although on different chromosomes or not in close physical proximity to the voltage gated sodium channel gene. In addition, pairwise comparison of resistance and susceptible mosquitoes from each population revealed differentiated genomic regions not associated with pyrethroid resistance. These new bi-allelic markers can be used to genotype other populations along with kdr alleles to understand their worldwide distribution. The functional roles of the genes near the newly discovered SNPs require new studies to determine if they act synergistically with kdr alleles or reduce the fitness cost of maintaining resistant alleles.

Recurring adaptive introgression of a supergene variant that determines social organization

Introgression has been proposed as an essential source of adaptive genetic variation. However, a key barrier to adaptive introgression is that recombination can break down combinations of alleles that underpin many traits. This barrier might be overcome in supergene regions, where suppressed recombination leads to joint inheritance across many loci. Here, we study the evolution of a large supergene region that determines a major social and ecological trait in Solenopsis fire ants: whether colonies have one queen or multiple queens. Using coalescent-based phylogenies built from the genomes of 365 haploid fire ant males, we show that the supergene variant responsible for multiple-queen colonies evolved in one species and repeatedly spread to other species through introgressive hybridization. This finding highlights how supergene architecture can enable a complex adaptive phenotype to recurrently permeate species boundaries.

Sympatric Populations of the Anopheles gambiae Complex in Southwest Burkina Faso Evolve Multiple Diverse Resistance Mechanisms in Response to Intense Selection Pressure with Pyrethroids

Simple Summary

Targeting mosquitoes with insecticides is one of the most effective methods to prevent malaria transmission. Although numbers of malaria cases have declined substantially this century, this pattern is not universal and Burkina Faso has one of the highest burdens of malaria; it is also a hotspot for the evolution of insecticide resistance in malaria vectors. We have established laboratory colonies from multiple species within the An. gambiae complex, the most efficient group of malaria vectors in the world, from larval collections in southwest Burkina Faso. Using bioassays with different insecticides widely used to control public health pests, we provide a profile of insecticide resistance in each of these colonies and, using molecular tools, reveal the genetic changes underpinning this resistance. We show that, whilst many resistance mechanisms are shared between species, there are some important differences which may affect resistance to current and future insecticide classes. The complexity, and diversity of resistance mechanisms highlights the importance of screening any potential new insecticide intended for use in malaria control against a wide range of populations. These stable laboratory colonies provide a valuable resource for insecticide discovery, and for further studies on the evolution and dispersal of insecticide resistance within and between species.

Abstract

Pyrethroid resistance in the Anopheles vectors of malaria is driving an urgent search for new insecticides that can be used in proven vector control tools such as insecticide treated nets (ITNs). Screening for potential new insecticides requires access to stable colonies of the predominant vector species that contain the major pyrethroid resistance mechanisms circulating in wild populations. Southwest Burkina Faso is an apparent hotspot for the emergence of pyrethroid resistance in species of the Anopheles gambiae complex. We established stable colonies from larval collections across this region and characterised the resistance phenotype and underpinning genetic mechanisms. Three additional colonies were successfully established (1 An. coluzzii, 1 An. gambiae and 1 An. arabiensis) to add to the 2 An. coluzzii colonies already established from this region; all 5 strains are highly resistant to pyrethroids. Synergism assays found that piperonyl butoxide (PBO) exposure was unable to fully restore susceptibility although exposure to a commercial ITN containing PBO resulted in 100% mortality. All colonies contained resistant alleles of the voltage gated sodium channel but with differing proportions of alternative resistant haplotypes. RNAseq data confirmed the role of P450s, with CYP6P3 and CYP6Z2 elevated in all 5 strains, and identified many other resistance mechanisms, some found across strains, others unique to a particular species. These strains represent an important resource for insecticide discovery and provide further insights into the complex genetic changes driving pyrethroid resistance.

Modelling spatiotemporal trends in the frequency of genetic mutations conferring insecticide target-site resistance in African mosquito malaria vector species

BACKGROUND

Resistance in malaria vectors to pyrethroids, the most widely used class of insecticides for malaria vector control, threatens the continued efficacy of vector control tools. Target-site resistance is an important genetic resistance mechanism caused by mutations in the voltage-gated sodium channel (Vgsc) gene that encodes the pyrethroid target-site. Understanding the geographic distribution of target-site resistance, and temporal trends across different vector species, can inform strategic deployment of vector control tools.

RESULTS

We develop a Bayesian statistical spatiotemporal model to interpret species-specific trends in the frequency of the most common resistance mutations, Vgsc-995S and Vgsc-995F, in three major malaria vector species Anopheles gambiae, An. coluzzii, and An. arabiensis over the period 2005-2017. The models are informed by 2418 observations of the frequency of each mutation in field sampled mosquitoes collected from 27 countries spanning western and eastern regions of Africa. For nine selected countries, we develop annual predictive maps which reveal geographically structured patterns of spread of each mutation at regional and continental scales. The results show associations, as well as stark differences, in spread dynamics of the two mutations across the three vector species. The coverage of ITNs was an influential predictor of Vgsc allele frequencies, with modelled relationships between ITN coverage and allele frequencies varying across species and geographic regions. We found that our mapped Vgsc allele frequencies are a significant partial predictor of phenotypic resistance to the pyrethroid deltamethrin in An. gambiae complex populations.

CONCLUSIONS

Our predictive maps show how spatiotemporal trends in insecticide target-site resistance mechanisms in African An. gambiae vary across individual vector species and geographic regions. Molecular surveillance of resistance mechanisms will help to predict resistance phenotypes and track their spread.

Predictors of genomic differentiation within a hybrid taxon

Hybridization is increasingly recognized as an important evolutionary force. Novel genetic methods now enable us to address how the genomes of parental species are combined in hybrid lineages. However, we still do not know the relative importance of admixed proportions, genome architecture and local selection in shaping hybrid genomes. Here, we take advantage of the genetically divergent island populations of Italian sparrow on Crete, Corsica and Sicily to investigate the predictors of genomic variation within a hybrid taxon. We test if differentiation is affected by recombination rate, selection, or variation in ancestry proportions. We find that the relationship between recombination rate and differentiation is less pronounced within hybrid lineages than between the parent species, as expected if purging of minor parent ancestry in low recombination regions reduces the variation available for differentiation. In addition, we find that differentiation between islands is correlated with differences in signatures of selection in two out of three comparisons. Signatures of selection within islands are correlated across all islands, suggesting that shared selection may mould genomic differentiation. The best predictor of strong differentiation within islands is the degree of differentiation from house sparrow, and hence loci with Spanish sparrow ancestry may vary more freely. Jointly, this suggests that constraints and selection interact in shaping the genomic landscape of differentiation in this hybrid species.

Phylogenetic analysis of the Neotropical Albitarsis Complex based on mitogenome data

BACKGROUND

Some of the most important malaria vectors in South America belong to the Albitarsis Complex (Culicidae; Anophelinae; Anopheles). Understanding the origin, nature, and geographical distribution of species diversity in this important complex has important implications for vector incrimination, control, and management, and for modelling future responses to climate change, deforestation, and human population expansion. This study attempts to further explore species diversity and evolutionary history in the Albitarsis Complex by undertaking a characterization and phylogenetic analysis of the mitogenome of all 10 putative taxa in the Albitarsis Complex.

METHODS

Mitogenome assembly and annotation allowed for feature comparison among Albitarsis Complex and Anopheles species. Selection analysis was conducted across all 13 protein-coding genes. Maximum likelihood and Bayesian inference methods were used to construct gene and species trees, respectively. Bayesian methods were also used to jointly estimate species delimitation and species trees.

RESULTS

Gene composition and order were conserved across species within the complex. Unique signatures of positive selection were detected in two species-Anopheles janconnae and An. albitarsis G-which may have played a role in the recent and rapid diversification of the complex. The COI gene phylogeny does not fully recover the mitogenome phylogeny, and a multispecies coalescent-based phylogeny shows that considerable uncertainty exists through much of the mitogenome species tree. The origin of divergence in the complex dates to the Pliocene/Pleistocene boundary, and divergence within the distinct northern South American clade is estimated at approximately 1 million years ago. Neither the phylogenetic trees nor the delimitation approach rejected the 10-species hypothesis, although the analyses could not exclude the possibility that four putative species with scant a priori support (An. albitarsis G, An. albitarsis H, An. albitarsis I, and An. albitarsis J), represent population-level, rather than species-level, splits.

CONCLUSION

The lack of resolution in much of the species tree and the limitations of the delimitation analysis warrant future studies on the complex using genome-wide data and the inclusion of additional specimens, particularly from two putative species, An. albitarsis I and An. albitarsis J.

Prevalence and Adaptive Impact of Introgression

Alleles that introgress between species can influence the evolutionary and ecological fate of species exposed to novel environments. Hybrid offspring of different species are often unfit, and yet it has long been argued that introgression can be a potent force in evolution, especially in plants. Over the last two decades, genomic data have increasingly provided evidence that introgression is a critically important source of genetic variation and that this additional variation can be useful in adaptive evolution of both animals and plants. Here, we review factors that influence the probability that foreign genetic variants provide long-term benefits (so-called adaptive introgression) and discuss their potential benefits. We find that introgression plays an important role in adaptive evolution, particularly when a species is far from its fitness optimum, such as when they expand their range or are subject to changing environments.

The genetic architecture of target-site resistance to pyrethroid insecticides in the African malaria vectors Anopheles gambiae and Anopheles coluzzii

Resistance to pyrethroid insecticides is a major concern for malaria vector control. Pyrethroids target the voltage-gated sodium channel (VGSC), an essential component of the mosquito nervous system. Substitutions in the amino acid sequence can induce a resistance phenotype. We use whole-genome sequence data from phase 2 of the Anopheles gambiae 1000 Genomes Project (Ag1000G) to provide a comprehensive account of genetic variation in the Vgsc gene across 13 African countries. In addition to known resistance alleles, we describe 20 other non-synonymous nucleotide substitutions at appreciable population frequency and map these variants onto a protein model to investigate the likelihood of pyrethroid resistance phenotypes. Thirteen of these novel alleles were found to occur almost exclusively on haplotypes carrying the known L995F kdr (knock-down resistance) allele and may enhance or compensate for the L995F resistance genotype. A novel mutation I1527T, adjacent to a predicted pyrethroid-binding site, was found in tight linkage with V402L substitutions, similar to allele combinations associated with resistance in other insect species. We also analysed genetic backgrounds carrying resistance alleles, to determine which alleles have experienced recent positive selection, and describe ten distinct haplotype groups carrying known kdr alleles. Five of these groups are observed in more than one country, in one case separated by over 3000 km, providing new information about the potential for the geographical spread of resistance. Our results demonstrate that the molecular basis of target-site pyrethroid resistance in malaria vectors is more complex than previously appreciated, and provide a foundation for the development of new genetic tools for insecticide resistance management.

An introgression breakthrough left by an anthropogenic contact between two ascidians

Human-driven translocations of species have diverse evolutionary consequences such as promoting hybridization between previously geographically isolated taxa. This is well illustrated by the solitary tunicate, Ciona robusta, native to the North East Pacific and introduced in the North East Atlantic. It is now co-occurring with its congener Ciona intestinalis in the English Channel, and C. roulei in the Mediterranean Sea. Despite their long allopatric divergence, first and second generation crosses showed a high hybridization success between the introduced and native taxa in the laboratory. However, previous genetic studies failed to provide evidence of recent hybridization between C. robusta and C. intestinalis in the wild. Using SNPs obtained from ddRAD-sequencing of 397 individuals from 26 populations, we further explored the genome-wide population structure of the native Ciona taxa. We first confirmed results documented in previous studies, notably (i) a chaotic genetic structure at regional scale, and (ii) a high genetic similarity between C. roulei and C. intestinalis, which is calling for further taxonomic investigation. More importantly, and unexpectedly, we also observed a genomic hotspot of long introgressed C. robusta tracts into C. intestinalis genomes at several locations of their contact zone. Both the genomic architecture of introgression, restricted to a 1.5 Mb region of chromosome 5, and its absence in allopatric populations suggest introgression is recent and occurred after the introduction of the non-native species. Overall, our study shows that anthropogenic hybridization can be effective in promoting introgression breakthroughs between species at a late stage of the speciation continuum.

The genomic consequences of hybridization

In the past decade, advances in genome sequencing have allowed researchers to uncover the history of hybridization in diverse groups of species, including our own. Although the field has made impressive progress in documenting the extent of natural hybridization, both historical and recent, there are still many unanswered questions about its genetic and evolutionary consequences. Recent work has suggested that the outcomes of hybridization in the genome may be in part predictable, but many open questions about the nature of selection on hybrids and the biological variables that shape such selection have hampered progress in this area. We synthesize what is known about the mechanisms that drive changes in ancestry in the genome after hybridization, highlight major unresolved questions, and discuss their implications for the predictability of genome evolution after hybridization.

Deep ancestral introgression shapes evolutionary history of dragonflies and damselflies

Introgression is an important biological process affecting at least 10% of the extant species in the animal kingdom. Introgression significantly impacts inference of phylogenetic species relationships where a strictly binary tree model cannot adequately explain reticulate net-like species relationships. Here we use phylogenomic approaches to understand patterns of introgression along the evolutionary history of a unique, non-model insect system: dragonflies and damselflies (Odonata). We demonstrate that introgression is a pervasive evolutionary force across various taxonomic levels within Odonata. In particular, we show that the morphologically "intermediate" species of Anisozygoptera (one of the three primary suborders within Odonata besides Zygoptera and Anisoptera), which retain phenotypic characteristics of the other two suborders, experienced high levels of introgression likely coming from zygopteran genomes. Additionally, we find evidence for multiple cases of deep inter-superfamilial ancestral introgression.

Genomic introgression from a distant congener in the Levant fritillary butterfly, Melitaea acentria

Introgressive hybridization is more common in nature than previously thought, and its role and creative power in evolution is hotly discussed but not completely understood. Introgression occurs more frequently in sympatry between recently diverged taxa, or when the speciation process has not yet been completed. However, there are relatively few documented cases of hybridization that erodes reproductive barriers between distantly related species. Here, we use whole genome and mitochondrial data to examine how introgression from a distant congener affects pattern of genetic differentiation in the Levant fritillary butterfly Melitaea acentria. We show that this local taxon has evolved as a peripatric geographic isolate of the widespread Melitaea persea, and that there has been significant unidirectional gene flow from the sympatric, nonclosely related Melitaea didyma to M. acentria. We found direct evidence of ongoing sporadic hybridization between M. didyma and M. acentria, which are separated by at least 5 million years of independent evolution. Elevated differentiation and lower level of introgression on the sex Z chromosome compared to autosomes suggest that the Z chromosome has accumulated loci acting as intrinsic postzygotic barriers. Our results show that introgression from M. didyma has been an additional source of nucleotide diversity in the M. acentria population, providing material for drift and selection.

Adaptive introgression of the beta-globin cluster in two Andean waterfowl

Introgression of beneficial alleles has emerged as an important avenue for genetic adaptation in both plant and animal populations. In vertebrates, adaptation to hypoxic high-altitude environments involves the coordination of multiple molecular and cellular mechanisms, including selection on the hypoxia-inducible factor (HIF) pathway and the blood-O2 transport protein hemoglobin (Hb). In two Andean duck species, a striking DNA sequence similarity reflecting identity by descent is present across the ~20 kb β-globin cluster including both embryonic (HBE) and adult (HBB) paralogs, though it was yet untested whether this is due to independent parallel evolution or adaptive introgression. In this study, we find that identical amino acid substitutions in the β-globin cluster that increase Hb-O2 affinity have likely resulted from historical interbreeding between high-altitude populations of two different distantly-related species. We examined the direction of introgression and discovered that the species with a deeper mtDNA divergence that colonized high altitude earlier in history (Anas flavirostris) transferred adaptive genetic variation to the species with a shallower divergence (A. georgica) that likely colonized high altitude more recently possibly following a range shift into a novel environment. As a consequence, the species that received these β-globin variants through hybridization might have adapted to hypoxic conditions in the high-altitude environment more quickly through acquiring beneficial alleles from the standing, hybrid-origin variation, leading to faster evolution.

Exchanging fluids The sociocultural implications of microbial, cultural, and ethnic admixture in Latin America

Knowledge of evolutionary influences on patterns of human mating, social interactions, and differential health is increasing, yet these insights have rarely been applied to historical analyses of human population dynamics. The genetic and evolutionary forces behind biases in interethnic mating and in the health of individuals of different ethnic groups in Latin America and the Caribbean since the European colonization of America are still largely ignored. We discuss how historical and contemporary sociocultural interactions and practices are strongly influenced by population-level evolutionary forces. Specifically, we discuss the historical implications of functional (de facto) polygyny, sex-biased admixture, and assortative mating in Latin America. We propose that these three evolutionary mechanisms influenced mating patterns, shaping the genetic and cultural landscape across Latin America and the Caribbean. Further, we discuss how genetic differences between the original populations that migrated at different times into Latin America contributed to their accommodation to and survival in the different local ecologies and interethnic interactions. Relevant medical and social implications follow from the genetic and cultural changes reviewed.

Gene Flow in Volant Vertebrates: Species Biology, Ecology and Climate Change

Gene flow, the exchange of genetic material between populations is an important biological process, which shapes and maintains biodiversity. The successful movement of individuals between populations depends on multiple factors determined by species biology and the environment. One of the most important factors regulating gene flow is the ability to move, and flight allows individuals to easily move across geographical barriers. Volant vertebrates are found on some of the remotest islands and contribute significantly to the biodiversity and ecosystem. The availability of next-generation sequencing data for non-model animals has substantially improved our understanding of gene flow and its consequences, allowing us to look at fine-scale patterns. However, most of our understanding regarding gene flow comes from the temperate regions and the Neotropics. The lack of studies from species-rich Asia is striking. In this review, we outline the importance of gene flow and the factors affecting gene flow, especially for volant vertebrates. We especially discuss research studies from tropical biomes of South and Southeast Asia, highlight the lacuna in literature and provide an outline for future studies in this species-rich region.

Unravelling population structure heterogeneity within the genome of the malaria vector Anopheles gambiae

Background

Whole genome re-sequencing provides powerful data for population genomic studies, allowing robust inferences of population structure, gene flow and evolutionary history. For the major malaria vector in Africa, Anopheles gambiae, other genetic aspects such as selection and adaptation are also important. In the present study, we explore population genetic variation from genome-wide sequencing of 765 An. gambiae and An. coluzzii specimens collected from across Africa. We used t-SNE, a recently popularized dimensionality reduction method, to create a 2D-map of An. gambiae and An. coluzzii genes that reflect their population structure similarities.

Results

The map allows intuitive navigation among genes distributed throughout the so-called “mainland” and numerous surrounding “island-like” gene clusters. These gene clusters of various sizes correspond predominantly to low recombination genomic regions such as inversions and centromeres, and also to recent selective sweeps. Because this mosquito species complex has been studied extensively, we were able to support our interpretations with previously published findings. Several novel observations and hypotheses are also made, including selective sweeps and a multi-locus selection event in Guinea-Bissau, a known intense hybridization zone between An. gambiae and An. coluzzii.

Conclusions

Our results present a rich dataset that could be utilized in functional investigations aiming to shed light onto An. gambiae s.l genome evolution and eventual speciation. In addition, the methodology presented here can be used to further characterize other species not so well studied as An. gambiae, shortening the time required to progress from field sampling to the identification of genes and genomic regions under unique evolutionary processes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07722-y.

Improving mosquito control strategies with population genomics

Mosquito control strategies increasingly apply knowledge from population genomics research. This review highlights recent applications to three research domains: mosquito invasions, insecticide resistance evolution, and rear and release programs. Current research trends follow developments in reference assemblies, either as improvements to existing assemblies (particularly Aedes) or assemblies for new taxa (particularly Anopheles). With improved assemblies, studies of invasive and rear and release target populations are better able to incorporate adaptive as well as demographic hypotheses. New reference assemblies are aiding comparisons of insecticide resistance across sister taxa while helping resolve taxon boundaries amidst frequent introgression. Anopheles gene drive deployments and improved Aedes genome assemblies should lead to a convergence in research aims for Anopheles and Aedes in the coming years.

Detecting adaptive introgression in human evolution using convolutional neural networks

Studies in a variety of species have shown evidence for positively selected variants introduced into a population via introgression from another, distantly related population—a process known as adaptive introgression. However, there are few explicit frameworks for jointly modelling introgression and positive selection, in order to detect these variants using genomic sequence data. Here, we develop an approach based on convolutional neural networks (CNNs). CNNs do not require the specification of an analytical model of allele frequency dynamics and have outperformed alternative methods for classification and parameter estimation tasks in various areas of population genetics. Thus, they are potentially well suited to the identification of adaptive introgression. Using simulations, we trained CNNs on genotype matrices derived from genomes sampled from the donor population, the recipient population and a related non-introgressed population, in order to distinguish regions of the genome evolving under adaptive introgression from those evolving neutrally or experiencing selective sweeps. Our CNN architecture exhibits 95% accuracy on simulated data, even when the genomes are unphased, and accuracy decreases only moderately in the presence of heterosis. As a proof of concept, we applied our trained CNNs to human genomic datasets—both phased and unphased—to detect candidates for adaptive introgression that shaped our evolutionary history.

Inferring Adaptive Introgression Using Hidden Markov Models

Adaptive introgression-the flow of adaptive genetic variation between species or populations-has attracted significant interest in recent years and it has been implicated in a number of cases of adaptation, from pesticide resistance and immunity, to local adaptation. Despite this, methods for identification of adaptive introgression from population genomic data are lacking. Here, we present Ancestry_HMM-S, a hidden Markov model-based method for identifying genes undergoing adaptive introgression and quantifying the strength of selection acting on them. Through extensive validation, we show that this method performs well on moderately sized data sets for realistic population and selection parameters. We apply Ancestry_HMM-S to a data set of an admixed Drosophila melanogaster population from South Africa and we identify 17 loci which show signatures of adaptive introgression, four of which have previously been shown to confer resistance to insecticides. Ancestry_HMM-S provides a powerful method for inferring adaptive introgression in data sets that are typically collected when studying admixed populations. This method will enable powerful insights into the genetic consequences of admixture across diverse populations. Ancestry_HMM-S can be downloaded from https://github.com/jesvedberg/Ancestry_HMM-S/.

Novel genotyping approaches to easily detect genomic admixture between the major Afrotropical malaria vector species, Anopheles coluzzii and An. gambiae

The two most efficient and most recently radiated Afrotropical vectors of human malaria - Anopheles coluzzii and An. gambiae - are identified by single-locus diagnostic PCR assays based on species-specific markers in a 4 Mb region on chromosome-X centromere. Inherently, these diagnostic assays cannot detect interspecific autosomal admixture shown to be extensive at the westernmost and easternmost extremes of the species range. The main aim of this study was to develop novel, easy-to-implement tools for genotyping An. coluzzii and An. gambiae-specific ancestral informative markers (AIMs) identified from the Anopheles gambiae 1000 genomes (Ag1000G) project. First, we took advantage of this large set of data in order to develop a multilocus approach to genotype 26 AIMs on all chromosome arms valid across the species range. Second, we tested the multilocus assay on samples from Guinea Bissau, The Gambia and Senegal, three countries spanning the westernmost hybridization zone, where conventional species diagnostic is problematic due to the putative presence of a novel "hybrid form". The multilocus assay was able to capture patterns of admixture reflecting those revealed by the whole set of AIMs and provided new original data on interspecific admixture in the region. Third, we developed an easy-to-use, cost-effective PCR approach for genotyping two AIMs on chromosome-3 among those included in the multilocus approach, opening the possibility for advanced identification of species and of admixed specimens during routine large scale entomological surveys, particularly, but not exclusively, at the extremes of the range, where WGS data highlighted unexpected autosomal admixture.

Impact of seasonality and malaria control interventions on Anopheles density and species composition from three areas of Uganda with differing malaria endemicity

BACKGROUND

Long-lasting insecticidal nets (LLINs) and indoor residual spraying (IRS) are the malaria control interventions primarily responsible for reductions in transmission intensity across sub-Saharan Africa. These interventions, however, may have differential impact on Anopheles species composition and density. This study examined the changing pattern of Anopheles species in three areas of Uganda with markedly different transmission intensities and different levels of vector control.

METHODS

From October 2011 to June 2016 mosquitoes were collected monthly using CDC light traps from 100 randomly selected households in three areas: Walukuba (low transmission), Kihihi (moderate transmission) and Nagongera (high transmission). LLINs were distributed in November 2013 in Walukuba and Nagongera and in June 2014 in Kihihi. IRS was implemented only in Nagongera, with three rounds of bendiocarb delivered between December 2014 and June 2015. Mosquito species were identified morphologically and by PCR (Polymerase Chain Reaction).

RESULTS

In Walukuba, LLIN distribution was associated with a decline in Anopheles funestus vector density (0.07 vs 0.02 mosquitoes per house per night, density ratio [DR] 0.34, 95% CI: 0.18-0.65, p = 0.001), but not Anopheles gambiae sensu stricto (s.s.) nor Anopheles arabiensis. In Kihihi, over 98% of mosquitoes were An. gambiae s.s. and LLIN distribution was associated with a decline in An. gambiae s.s. vector density (4.00 vs 2.46, DR 0.68, 95% CI: 0.49-0.94, p = 0.02). In Nagongera, the combination of LLINs and multiple rounds of IRS was associated with almost complete elimination of An. gambiae s.s. (28.0 vs 0.17, DR 0.004, 95% CI: 0.002-0.009, p < 0.001), and An. funestus sensu lato (s.l.) (3.90 vs 0.006, DR 0.001, 95% CI: 0.0005-0.004, p < 0.001), with a less pronounced decline in An. arabiensis (9.18 vs 2.00, DR 0.15 95% CI: 0.07-0.33, p < 0.001).

CONCLUSIONS

LLIN distribution was associated with reductions in An. funestus s.l. in the lowest transmission site and An. gambiae s.s. in the moderate transmission site. In the highest transmission site, a combination of LLINs and multiple rounds of IRS was associated with the near collapse of An. gambiae s.s. and An. funestus s.l. Following IRS, An. arabiensis, a behaviourally resilient vector, became the predominant species, which may have implications for malaria vector control activities. Development of interventions targeted at outdoor biting remains a priority.

Control of malaria-transmitting mosquitoes using gene drives

Gene drives are selfish genetic elements that can be re-designed to invade a population and they hold tremendous potential for the control of mosquitoes that transmit disease. Much progress has been made recently in demonstrating proof of principle for gene drives able to suppress populations of malarial mosquitoes, or to make them refractory to the Plasmodium parasites they transmit. This has been achieved using CRISPR-based gene drives. In this article, I will discuss the relative merits of this type of gene drive, as well as barriers to its technical development and to its deployment in the field as malaria control. This article is part of the theme issue 'Novel control strategies for mosquito-borne diseases'.

Spatial population genomics of a recent mosquito invasion

Population genomic approaches can characterize dispersal across a single generation through to many generations in the past, bridging the gap between individual movement and intergenerational gene flow. These approaches are particularly useful when investigating dispersal in recently altered systems, where they provide a way of inferring long-distance dispersal between newly established populations and their interactions with existing populations. Human-mediated biological invasions represent such altered systems which can be investigated with appropriate study designs and analyses. Here we apply temporally restricted sampling and a range of population genomic approaches to investigate dispersal in a 2004 invasion of Aedes albopictus (the Asian tiger mosquito) in the Torres Strait Islands (TSI) of Australia. We sampled mosquitoes from 13 TSI villages simultaneously and genotyped 373 mosquitoes at genome-wide single nucleotide polymorphisms (SNPs): 331 from the TSI, 36 from Papua New Guinea (PNG) and four incursive mosquitoes detected in uninvaded regions. Within villages, spatial genetic structure varied substantially but overall displayed isolation by distance and a neighbourhood size of 232-577. Close kin dyads revealed recent movement between islands 31-203 km apart, and deep learning inferences showed incursive Ae. albopictus had travelled to uninvaded regions from both adjacent and nonadjacent islands. Private alleles and a co-ancestry matrix indicated direct gene flow from PNG into nearby islands. Outlier analyses also detected four linked alleles introgressed from PNG, with the alleles surrounding 12 resistance-associated cytochrome P450 genes. By treating dispersal as both an intergenerational process and a set of discrete events, we describe a highly interconnected invasive system.

Laboratory and microcosm experiments reveal contrasted adaptive responses to ammonia and water mineralisation in aquatic stages of the sibling species Anopheles gambiae (sensu stricto) and Anopheles coluzzii

Background

The sibling species of the malaria mosquito, Anopheles gambiae (sensu stricto) and Anopheles coluzzii co-exist in many parts of West Africa and are thought to have recently diverged through a process of ecological speciation with gene flow. Divergent larval ecological adaptations, resulting in Genotype-by-Environment (G × E) interactions, have been proposed as important drivers of speciation in these species. In West Africa, An. coluzzii tends to be associated with permanent man-made larval habitats such as irrigated rice fields, which are typically more eutrophic and mineral and ammonia-rich than the temporary rain pools exploited by An. gambiae (s.s.)

Methods

To highlight G × E interactions at the larval stage and their possible role in ecological speciation of these species, we first investigated the effect of exposure to ammonium hydroxide and water mineralisation on larval developmental success. Mosquito larvae were exposed to two water sources and increasing ammonia concentrations in small containers until adult emergence. In a second experiment, larval developmental success was compared across two contrasted microcosms to highlight G × E interactions under conditions such as those found in the natural environment.

Results

The first experiment revealed significant G × E interactions in developmental success and phenotypic quality for both species in response to increasing ammonia concentrations and water mineralisation. The An. coluzzii strain outperformed the An. gambiae (s.s.) strain under limited conditions that were closer to more eutrophic habitats. The second experiment revealed divergent crisscrossing reaction norms in the developmental success of the sibling species in the two contrasted larval environments. As expected, An. coluzzii had higher emergence rates in the rice paddy environment with emerging adults of superior phenotypic quality compared to An. gambiae (s.s.), and vice versa, in the rain puddle environment.

Conclusions

Evidence for such G × E interactions lends support to the hypothesis that divergent larval adaptations to the environmental conditions found in man-made habitats such as rice fields in An. coluzzii may have been an important driver of its ecological speciation.

Rapid parallel adaptation despite gene flow in silent crickets

Gene flow is predicted to impede parallel adaptation via de novo mutation, because it can introduce pre-existing adaptive alleles from population to population. We test this using Hawaiian crickets (Teleogryllus oceanicus) in which 'flatwing' males that lack sound-producing wing structures recently arose and spread under selection from an acoustically-orienting parasitoid. Morphometric and genetic comparisons identify distinct flatwing phenotypes in populations on three islands, localized to different loci. Nevertheless, we detect strong, recent and ongoing gene flow among the populations. Using genome scans and gene expression analysis we find that parallel evolution of flatwing on different islands is associated with shared genomic hotspots of adaptation that contain the gene doublesex, but the form of selection differs among islands and corresponds to known flatwing demographics in the wild. We thus show how parallel adaptation can occur on contemporary timescales despite gene flow, indicating that it could be less constrained than previously appreciated.

Evidence for Divergent Selection on Immune Genes between the African Malaria Vectors, Anopheles coluzzii and A. gambiae

Simple Summary

A comparison of the genomes of the African malaria vectors, Anopheles gambiae and A. coluzzii, revealed that immune genes are highly diverged. Although these two species frequently co-occur within a single site, they occur in distinct larval habitats. Our results taken in the context of known differences in the larval habitats occupied by these taxa support the hypothesis that observed genetic divergence may be driven by immune response to microbial agents specific to these habitats. Strict within species mating may have subsequently evolved in part to maintain immunocompetence which might be compromised by dysregulation of immune pathways in hybrids. We conclude that the evolution of immune gene divergence among this important group of species may serve as a useful model to explore ecological speciation in general.

Abstract

During their life cycles, microbes infecting mosquitoes encounter components of the mosquito anti-microbial innate immune defenses. Many of these immune responses also mediate susceptibility to malaria parasite infection. In West Africa, the primary malaria vectors are Anopheles coluzzii and A. gambiae sensu stricto, which is subdivided into the Bamako and Savanna sub-taxa. Here, we performed whole genome comparisons of the three taxa as well as genotyping of 333 putatively functional SNPs located in 58 immune signaling genes. Genome data support significantly higher differentiation in immune genes compared with a randomly selected set of non-immune genes among the three taxa (permutation test p < 0.001). Among the 58 genes studied, the majority had one or more segregating mutations (72.9%) that were significantly diverged among the three taxa. Genes detected to be under selection include MAP2K4 and Raf. Despite the genome-wide distribution of immune genes, a high level of linkage disequilibrium (r² > 0.8) was detected in over 27% of SNP pairs. We discuss the potential role of immune gene divergence as adaptations to the different larval habitats associated with A. gambiae taxa and as a potential force driving ecological speciation in this group of mosquitoes.

A population genomic unveiling of a new cryptic mosquito taxon within the malaria-transmitting Anopheles gambiae complex

The Anopheles gambiae complex consists of multiple morphologically indistinguishable mosquito species including the most important vectors of the malaria parasite Plasmodium falciparum in sub-Saharan Africa. Nine cryptic species have been described so far within the complex. The ecological, immunological and reproductive differences among these species will critically impact population responses to disease control strategies and environmental changes. Here, we examine whole-genome sequencing data from a longitudinal study of putative A. coluzzii in western Burkina Faso. Surprisingly, many specimens are genetically divergent from A. coluzzii and all other Anopheles species and represent a new taxon, here designated Anopheles TENGRELA (AT). Population genetic analysis suggests that the cryptic GOUNDRY subgroup, previously collected as larvae in central Burkina Faso, represents an admixed population descended from both A. coluzzii and AT. AT harbours low nucleotide diversity except for the 2La inversion polymorphism which is maintained by overdominance. It shows numerous fixed differences with A. coluzzii concentrated in several regions reflecting selective sweeps, but the two taxa are identical at standard diagnostic loci used for taxon identification, and thus, AT may often go unnoticed. We present an amplicon-based genotyping assay for identifying AT which could be usefully applied to numerous existing samples. Misidentified cryptic taxa could seriously confound ongoing studies of Anopheles ecology and evolution in western Africa, including phenotypic and genotypic surveys of insecticide resistance. Reproductive barriers between cryptic species may also complicate novel vector control efforts, for example gene drives, and hinder predictions about evolutionary dynamics of Anopheles and Plasmodium.