Experimental swap of Anopheles gambiae's assortative mating preferences demonstrates key role of X-chromosome divergence island in incipient sympatric speciation

Unravelling the role of mitochondrial DNA in hybrid incompatibility within species of the Anopheles gambiae complex

Isolation mechanisms between mosquito species of the Anopheles gambiae complex, which includes major malaria vectors, remain poorly understood. In some cases, pre-zygotic barriers have been shown to limit gene flow between species of the complex, leading to a low level of hybridisation in nature. Post-zygotic mechanisms manifest with F1 hybrid males fully sterile and F1 hybrid females with reduced fertility. Genetic approaches combined with DNA sequencing techniques have highlighted the involvement of genomic regions in hybrid incompatibility with a predominant role of the X chromosome. In addition, differences in the phenotype of F1 hybrid males have been identified depending on the directionality of the parental cross used to generate them. All these studies have focused on the interaction of nuclear DNA elements in hybrid individuals. Given the role that mitochondrial DNA plays in genetic incompatibilities within other organisms and its unique inheritance pattern, commonly maternal, we conducted a genetic study that relied on the introgression of mitochondrial DNA between Anopheles gambiae and Anopheles arabiensis. The findings indicate that the mitochondrial switch does not appear to restore the fertility of F1 hybrid males, suggesting that mitochondrial DNA may not play a role in hybrid incompatibilities in these Anopheles species.

High Levels of Admixture in Anopheles gambiae Populations from Côte d'Ivoire Revealed by Multilocus Genotyping

Anopheles coluzzii and An. gambiae-the two most recently radiated species of the An. gambiae complex and the major Afrotropical malaria vector species-are identified by markers in the X-centromeric IGS rDNA region. Putative IGS-hybrids are rarely found in the field, except in restricted areas where genomic studies have led to the hypothesis that the observed IGS-patterns are due to cryptic taxa rather than to hybridization between the two species. We investigated the genome-wide levels of admixture in two villages in Côte d'Ivoire where high levels of IGS-hybrids have been detected, confirming unparalleled high frequencies in the coastal village. Genotyping of 24 Ancestry Informative Markers (AIMs) along the three chromosomes produced discordant results between the IGS-marker and the multilocus genotype obtained for AIMs across the whole genome (29%) as well as AIMs on chromosome-X (considered to be fundamental for species reproductive isolation) only (21%). Results highlight a complicated pattern of admixture that deserves deeper genomic analyses to understand better possible underlying causes (from extensive processes of hybridization to the existence of different cryptic taxa), and stress the need of developing advanced diagnostics for An. coluzzii, An. gambiae and putative new taxa, instrumental for assessing taxon-specific epidemiological characters.

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.

Swarming Behavior in Anopheles gambiae (sensu lato): Current Knowledge and Future Outlook

Effective management of insect disease vectors requires a detailed understanding of their ecology and behavior. In Anopheles gambiae sensu lato (s.l.) (Diptera: Culicidae) mating occurs during swarming, but knowledge of their mating behavior under natural conditions is limited. Mosquitoes mate in flight over specific landmarks, known as swarm markers, at particular locations. Swarms consist of males; the females usually approach the swarm and depart following copulation. The number of mating pairs per swarm is closely associated with swarm size. The shape and height of swarm markers vary and may depend on the environmental conditions at the swarm's location. Male-male interactions in mosquito swarms with similar levels of attractive flight activity can offer a mating advantage to some individuals. Flight tone is used by mosquitoes to recognize the other sex and choose a desirable mate. Clarifying these and other aspects of mosquito reproductive behavior can facilitate the development of population control measures that target swarming sites. This review describes what is currently known about swarming behavior in Anopheles gambiae s.l., including swarm characteristics; mating within and outside of swarms, insemination in females, and factors affecting and stimulating swarming.

Asymmetric Phenotypes of Sterile Hybrid Males From Reciprocal Crosses Between Species of the Anopheles gambiae Complex

Haldane’s rule of speciation states that sterility or inviability affects the heterogametic sex of inter-species hybrids. Darwin’s corollary to Haldane’s rule implies that there are asymmetric phenotypes in inter-species hybrids from reciprocal crosses. Studying the phenotypes of F1 hybrids among closely related species of malaria mosquitoes can assist researchers in identifying the genetic factors and molecular mechanisms of speciation. To characterize phenotypes of sterile hybrid males in the Anopheles gambiae complex, we performed crosses between laboratory strains of An. merus and either An. gambiae or An. coluzzii. The reproductive tracts had normal external morphology in hybrid males from crosses between female An. merus and male An. gambiae or An. coluzzii. Despite being sterile, these males could copulate with females for a normal period of time and could transfer a mating plug to induce female oviposition and monogamy. In contrast, the entire reproductive tracts in hybrid males from crosses between female An. gambiae or An. coluzzii and male An. merus were severely underdeveloped. These males had atrophic testes and reduced somatic organs of the reproductive system including male accessary glands and ejaculatory duct. In addition, hybrid males with underdeveloped reproductive tracts displayed a shorter copulation time with females and failed to induce female oviposition and monogamy due to their inability to form and transfer a plug to females during mating. The asymmetry of the phenotypes associated with hybrid male sterility suggests that different genetic factors and molecular mechanisms are responsible for reproductive isolation in reciprocal crosses among species of the An. gambiae complex.

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.

Chromosome-level genome assemblies of the malaria vectors Anopheles coluzzii and Anopheles arabiensis

BACKGROUND

Anopheles coluzzii and Anopheles arabiensis belong to the Anopheles gambiae complex and are among the major malaria vectors in sub-Saharan Africa. However, chromosome-level reference genome assemblies are still lacking for these medically important mosquito species.

FINDINGS

In this study, we produced de novo chromosome-level genome assemblies for A. coluzzii and A. arabiensis using the long-read Oxford Nanopore sequencing technology and the Hi-C scaffolding approach. We obtained 273.4 and 256.8 Mb of the total assemblies for A. coluzzii and A. arabiensis, respectively. Each assembly consists of 3 chromosome-scale scaffolds (X, 2, 3), complete mitochondrion, and unordered contigs identified as autosomal pericentromeric DNA, X pericentromeric DNA, and Y sequences. Comparison of these assemblies with the existing assemblies for these species demonstrated that we obtained improved reference-quality genomes. The new assemblies allowed us to identify genomic coordinates for the breakpoint regions of fixed and polymorphic chromosomal inversions in A. coluzzii and A. arabiensis.

CONCLUSION

The new chromosome-level assemblies will facilitate functional and population genomic studies in A. coluzzii and A. arabiensis. The presented assembly pipeline will accelerate progress toward creating high-quality genome references for other disease vectors.

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.

Unveiling mosquito cryptic species and their reproductive isolation

Mosquitoes are major vectors of many infectious pathogens or parasites. Understanding cryptic species and the speciation of disease vectors has important implications for vector management, evolution and host-pathogen and/or host-parasite interactions. Currently, mosquito cryptic species have been reported in many studies, most of which focus on the reproductive isolation of cryptic species and mainly on Anopheles gambiae sensu lato complex. Emerging species within the primary malaria vector Anopheles gambiae show different ecological preferences and significant prezygotic reproductive isolation, while Aedes mariae and Aedes zammitii show postmating reproductive isolation. However, data reporting the reproductive isolation in Culex and Aedes albopictus mosquito cryptic species is absent. The lack of systematic studies leaves many questions open, such as whether cryptic species are more common in particular habitats, latitudes or taxonomic groups; what mosquito cryptic species evolutionary processes bring about reproductive isolation in the absence of morphological differentiation? How does Wolbachia infection affect in mosquitoes' reproductive isolation? In this review, we provide a summary of recent advances in the discovery and identification of sibling or cryptic species within mosquito genera.

Assessment of the ecologically dependent post-zygotic isolation between Anopheles coluzzii and Anopheles gambiae

Within the Anopheles gambiae complex, the sibling species An. coluzzii and An. gambiae are undergoing sympatric speciation. These species are characterized by rare hybrids in most of their geographical distribution. A strong assortative mating mediated by spatial swarm segregation has been shown whereas no intrinsic post-zygotic barriers have been found in laboratory conditions. To test the role of the hybridisation in reproductive isolation in natural populations transplant experiment are therefore needed to establish the significance of post-zygotic barriers. Previous studies indicated that predation is one of the major forces driving ecological divergence between An. gambiae and An. coluzzii. Here we extended these studies to their hybrids. Parental species and their F1 hybrids from reciprocal crosses were generated by the forced-mating technique as follows: female An. coluzzii × male An. coluzzii; female An. coluzzii × male An. gambiae; female An. gambiae × male An. coluzzii and female An. gambiae × Male An. gambiae. First instar larvae of each group from the crossing (here after An. coluzzii, Hybrid COL/GAM, Hybrid GAM/COL and An. gambiae, respectively) were transplanted in a field experiment with predation effect. Emergence success, development time of larvae and body size of the newly emerging adults were estimated as fitness components and then compared between parental species and F1 hybrids in absence and in presence of predators. Our findings confirm that An. coluzzii had higher fitness than An. gambiae in presence of predators versus in absence of predators. Moreover, the fitness of the F1 hybrid COL/GAM whose female parent was An. coluzzii matched that of An. coluzzii while that of the F1 reciprocal hybrid GAM/COL was similar to An. gambiae.

Structural Variation of the X Chromosome Heterochromatin in the Anopheles gambiae Complex

Heterochromatin is identified as a potential factor driving diversification of species. To understand the magnitude of heterochromatin variation within the Anopheles gambiae complex of malaria mosquitoes, we analyzed metaphase chromosomes in An. arabiensis, An. coluzzii, An. gambiae, An. merus, and An. quadriannulatus. Using fluorescence in situ hybridization (FISH) with ribosomal DNA (rDNA), a highly repetitive fraction of DNA, and heterochromatic Bacterial Artificial Chromosome (BAC) clones, we established the correspondence of pericentric heterochromatin between the metaphase and polytene X chromosomes of An. gambiae. We then developed chromosome idiograms and demonstrated that the X chromosomes exhibit qualitative differences in their pattern of heterochromatic bands and position of satellite DNA (satDNA) repeats among the sibling species with postzygotic isolation, An. arabiensis, An. merus, An. quadriannulatus, and An. coluzzii or An. gambiae. The identified differences in the size and structure of the X chromosome heterochromatin point to a possible role of repetitive DNA in speciation of mosquitoes. We found that An. coluzzii and An. gambiae, incipient species with prezygotic isolation, share variations in the relative positions of the satDNA repeats and the proximal heterochromatin band on the X chromosomes. This previously unknown genetic polymorphism in malaria mosquitoes may be caused by a differential amplification of DNA repeats or an inversion in the sex chromosome heterochromatin.

Premeiotic and meiotic failures lead to hybrid male sterility in the Anopheles gambiae complex

Hybrid male sterility (HMS) contributes to speciation by restricting gene flow between related taxa. Detailed cytological characterization of reproductive organs in hybrid males is important for identifying phenotypes that can help guide searches of speciation genes. To investigate possible cellular causes of HMS, we performed crosses between closely related species of the Anopheles gambiae complex: An. merus with An. gambiae or An. coluzzii. We demonstrate that HMS in African malaria mosquitoes involves two defects in the reciprocal crosses: a premeiotic arrest of germline stem cells in degenerate testes and a failure of the reductional meiotic division of primary spermatocytes in normal-like testes. The premeiotic arrest in degenerate testes of hybrids is accompanied by a strong suppression of meiotic and postmeiotic genes. Unlike pure species, sex chromosomes in normal-like testes of F1 hybrids are largely unpaired during meiotic prophase I and all chromosomes show various degrees of insufficient condensation. Instead of entering reductional division in meiosis I, primary spermatocytes prematurely undergo an equational mitotic division producing non-motile diploid sperm. Thus, our study identified cytogenetic errors in interspecies hybrids that arise during the early stages of postzygotic isolation.

Effect of water source and feed regime on development and phenotypic quality in Anopheles gambiae (s.l.): prospects for improved mass-rearing techniques towards release programmes

Background

In many malaria-endemic sub-Saharan countries, insecticide resistance poses a threat to existing mosquito control measures, underscoring the need for complementary control methods such as sterile and/or genetically-modified mosquito release programmes. The sibling species Anopheles gambiae and An. coluzzii are responsible for malaria transmission in most of this region. In their natural habitat, these species generally breed in clean, soft water and it is believed that divergent preference in their larval breeding sites have played a role in their speciation process. Mosquito release programmes rely on the rearing of mosquitoes at high larval densities. Current rearing protocols often make use of deionised water regardless of the strain reared. They also depend on a delicate balance between the need for adequate feeding and the negative effect of toxic ammonia and food waste build-up on mosquito development, making managing and improving water quality in the insectary imperative.

Methods

Here, we investigated the impact of water source and feed regimes on emergence rate and phenotypic quality of mosquitoes in the insectary. First-instar larvae of An. gambiae (Kisumu strain) and An. coluzzii (Mopti and VK3 strains) were reared in three water sources with varying degrees of hardness (deionised, mineral and a mix of the two), with a daily water change. Larvae were fed daily using two standardised feeding regimes, solution and powder feed.

Results

Water source had a significant impact on mosquito size and development time for all strains. Earlier emergence of significantly larger mosquitoes was observed in mineral water with the smallest mosquitoes developing later from deionised water. Wing-length was significantly longer in mineral, mixed water and in powder feed, irrespective of sex, strains or water types. Deionised water was the least favourable for mosquito quality across all strains.

Conclusions

Mineral water and powder feed should be used in rearing protocols to improve mosquito quality where the optimal quality of mosquitoes is desired. Although results obtained were not significant for improved mosquito numbers, the phenotypic quality of mosquitoes reared was significantly improved in mineral water and mix water. Further studies are recommended on the impact mineral water has on other fitness traits such as longevity, fecundity and mating competitiveness.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3465-0) contains supplementary material, which is available to authorized users.

A High-Quality De novo Genome Assembly from a Single Mosquito Using PacBio Sequencing

A high-quality reference genome is a fundamental resource for functional genetics, comparative genomics, and population genomics, and is increasingly important for conservation biology. PacBio Single Molecule, Real-Time (SMRT) sequencing generates long reads with uniform coverage and high consensus accuracy, making it a powerful technology for de novo genome assembly. Improvements in throughput and concomitant reductions in cost have made PacBio an attractive core technology for many large genome initiatives, however, relatively high DNA input requirements (~5 µg for standard library protocol) have placed PacBio out of reach for many projects on small organisms that have lower DNA content, or on projects with limited input DNA for other reasons. Here we present a high-quality de novo genome assembly from a single Anopheles coluzzii mosquito. A modified SMRTbell library construction protocol without DNA shearing and size selection was used to generate a SMRTbell library from just 100 ng of starting genomic DNA. The sample was run on the Sequel System with chemistry 3.0 and software v6.0, generating, on average, 25 Gb of sequence per SMRT Cell with 20 h movies, followed by diploid de novo genome assembly with FALCON-Unzip. The resulting curated assembly had high contiguity (contig N50 3.5 Mb) and completeness (more than 98% of conserved genes were present and full-length). In addition, this single-insect assembly now places 667 (>90%) of formerly unplaced genes into their appropriate chromosomal contexts in the AgamP4 PEST reference. We were also able to resolve maternal and paternal haplotypes for over 1/3 of the genome. By sequencing and assembling material from a single diploid individual, only two haplotypes were present, simplifying the assembly process compared to samples from multiple pooled individuals. The method presented here can be applied to samples with starting DNA amounts as low as 100 ng per 1 Gb genome size. This new low-input approach puts PacBio-based assemblies in reach for small highly heterozygous organisms that comprise much of the diversity of life.

A genetic linkage map for the salmon louse (Lepeophtheirus salmonis): evidence for high male:female and inter-familial recombination rate differences

A salmon louse (Lepeophtheirus salmonis salmonis) genetic linkage map was constructed to serve as a genomic resource for future investigations into the biology of this important marine parasitic copepod species, and to provide insights into the inheritance patterns of genetic markers in this species. SNP genotyping of 8 families confirmed the presence of 15 linkage groups based upon the assignment of 93,773 markers. Progeny sample size weight adjusted map sizes in males (with the exception of SL12 and SL15) ranged in size from 96.50 cM (SL11) to 134.61 cM (SL06), and total combined map steps or bins ranged from 143 (SL09) to 203 (SL13). The SL12 male map was the smallest linkage group with a weight-averaged size of 3.05 cM with 6 recombination bins. Male:female specific recombination rate differences are 10.49:1 and represent one of the largest reported sex-specific differences for any animal species. Recombination ratio differences (M:F) ranged from 1.0 (SL12) to 29:1 (SL15). The number of markers exhibiting normal Mendelian segregation within the sex linkage group SL15 was extremely low (N = 80) in comparison to other linkage groups genotyped [range: 1459 (SL12)-10206 markers (SL05)]. Re-evaluation of Mendelian inheritance patterns of markers unassigned to any mapping parent according to hemizygous segregation patterns (models presented) identified matches for many of these markers to hemizygous patterns. The greatest proportion of these markers assigned to SL15 (N increased to 574). Inclusion of the hemizygous markers revised SL15 sex-specific recombination rate differences to 28:1. Recombination hot- and coldspots were identified across all linkage groups with all linkage groups possessing multiple peaks. Nine of 13 linkage groups evaluated possessed adjacent domains with hot-coldspot transitional zones. The most common pattern was for one end of the linkage to show elevated recombination in addition to internal regions. For SL01 and SL06, however, a terminal region with high recombination was not evident while a central domain possessing extremely high-recombination levels was present. High levels of recombination were weakly coupled to higher levels of SNP variation within domains, but this association was very strong for the central domains of SL01 and SL06. From the pooled paternal half-sib lots (several virgin females placed with 1 male), only 1 or two surviving family lots were obtained. Surviving families possessed parents where both the male and female possessed either inherently low or high recombination rates. This study provides insight into the organization of the sea louse genome, and describes large differences in recombination rate that exist among individuals of the same sex, and between the sexes. These differences in recombination rate may be coupled to the capabilities of this species to adapt to environmental and pharmaceutical treatments, given that family survivorship appears to be enhanced when parents have similar recombination levels.

Chemosensory gene expression in olfactory organs of the anthropophilic Anopheles coluzzii and zoophilic Anopheles quadriannulatus

BACKGROUND

Anopheles (An.) coluzzii, one of Africa's primary malaria vectors, is highly anthropophilic. This human host preference contributes greatly to its ability to transmit malaria. In contrast, the closely related An. quadriannulatus prefers to feed on bovids and is not thought to contribute to malaria transmission. The diverged preference for host odor profiles between these sibling species is likely reflected in chemosensory gene expression levels in the olfactory organs. Therefore, we compared the transcriptomes of the antennae and maxillary palps between An. coluzzii and An. quadriannulatus, focusing on the major chemosensory gene families.

RESULTS

While chemosensory gene expression is strongly correlated between the two species, various chemosensory genes show significantly enhanced expression in one of the species. In the antennae of An. coluzzii the expression of six olfactory receptors (Ors) and seven ionotropic receptors (Irs) is considerably enhanced, whereas 11 Ors and 3 Irs are upregulated in An. quadriannulatus. In the maxillary palps, leaving aside Irs with very low level of expression, one Ir is strongly enhanced in each species. In addition, we find divergence in odorant binding protein (Obp) gene expression, with several highly expressed Obps being enhanced in the antennae and palps of An. coluzzii. Finally, the expression of several gustatory receptors (Grs) in the palps appears to be species-specific, including a homolog of a sugar-sensing Drosophila Gr.

CONCLUSIONS

A considerable number of Ors and Irs are differentially expressed between these two closely related species with diverging host preference. These chemosensory genes could play a role in the human host preference of the malaria vector An. coluzzii. Additionally, divergence in Obp expression between the two species suggests a possible role of these odor carrier proteins in determining host preference. Finally, divergence in chemosensory expression in the palps may point towards a possible role for the maxillary palps in host differentiation.

Pre-copula acoustic behaviour of males in the malarial mosquitoes Anopheles coluzzii and Anopheles gambiae s.s. does not contribute to reproductive isolation

We reveal that males of two members of the Anopheles gambiae s.l. species complex, Anopheles coluzzii and Anopheles gambiae s.s. (hereafter A. gambiae), which are both malaria vectors, perform a stereotypical acoustic behaviour in response to pure tones at frequencies that encompass the frequency range of the female's flight-tones. This behaviour resembles that described for Culex quinquefasciatus and consists of phonotactic flight initiated by a steep increase in wing-beat frequency (WBF) followed by rapid frequency modulation (RFM) of WBF when in close proximity to the sound source. RFM was elicited without acoustic feedback or the presence of a live female, but it appears to be a stereotypic behaviour in the immediate lead up to copula formation. RFM is an independent and different behavioural process from harmonic convergence interactions used by male-female pairs for mate recognition at earlier stages of mating. Acoustic threshold for RFM was used to plot behavioural audiograms from free-flying A coluzzii and A gambiae males. These audiograms were almost identical (minima ∼400 Hz) and encompassed the WBF ranges of A coluzzii (378-601 Hz) and A gambiae (373-590 Hz) females, indicating that males of the two species share similar frequency tuning and range. Furthermore, no differences were found between the two species in their WBFs, RFM behaviour or harmonic convergence ratios. These results indicate that assortative mating between A coluzzii and A gambiae is unlikely to be based on male-specific acoustic behaviours during RFM. The significance of these findings in relation to possible mechanisms for assortative mating is discussed.

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.

A combination of sexual and ecological divergence contributes to rearrangement spread during initial stages of speciation

Chromosomal rearrangements between sympatric species often contain multiple loci contributing to assortative mating, local adaptation and hybrid sterility. When and how these associations arise during the process of speciation remains a subject of debate. Here, we address the relative roles of local adaptation and assortative mating on the dynamics of rearrangement evolution by studying how a rearrangement covaries with sexual and ecological trait divergence within a species. Previously, a chromosomal rearrangement that suppresses recombination on the Z (sex) chromosome was identified in European corn borer moths (Ostrinia nubilalis). We further characterize this recombination suppressor and explore its association with variation in sex pheromone communication and seasonal ecological adaptation in pairs of populations that are divergent in one or both of these characteristics. Direct estimates of recombination suppression in pedigree mapping families indicated that more than 39% of the Z chromosome (encompassing up to ~10 megabases and ~300 genes) resides within a nonrecombining unit, including pheromone olfactory receptor genes and a major quantitative trait locus that contributes to ecotype differences (Pdd). Combining direct and indirect estimates of recombination suppression, we found that the rearrangement was occasionally present between sexually isolated strains (E vs. Z) and between divergent ecotypes (univoltine vs. bivoltine). However, it was only consistently present when populations differed in both sexual and ecological traits. Our results suggest that independent of the forces that drove the initial establishment of the rearrangement, a combination of sexual and ecological divergence is required for rearrangement spread during speciation.

Chromosomal inversions and ecotypic differentiation in Anopheles gambiae: the perspective from whole-genome sequencing

The molecular mechanisms and genetic architecture that facilitate adaptive radiation of lineages remain elusive. Polymorphic chromosomal inversions, due to their recombination-reducing effect, are proposed instruments of ecotypic differentiation. Here, we study an ecologically diversifying lineage of Anopheles gambiae, known as the Bamako chromosomal form based on its unique complement of three chromosomal inversions, to explore the impact of these inversions on ecotypic differentiation. We used pooled and individual genome sequencing of Bamako, typical (non-Bamako) An. gambiae and the sister species Anopheles coluzzii to investigate evolutionary relationships and genomewide patterns of nucleotide diversity and differentiation among lineages. Despite extensive shared polymorphism and limited differentiation from the other taxa, Bamako clusters apart from the other taxa, and forms a maximally supported clade in neighbour-joining trees based on whole-genome data (including inversions) or solely on collinear regions. Nevertheless, F_ST outlier analysis reveals that the majority of differentiated regions between Bamako and typical An. gambiae are located inside chromosomal inversions, consistent with their role in the ecological isolation of Bamako. Exceptionally differentiated genomic regions were enriched for genes implicated in nervous system development and signalling. Candidate genes associated with a selective sweep unique to Bamako contain substitutions not observed in sympatric samples of the other taxa, and several insecticide resistance gene alleles shared between Bamako and other taxa segregate at sharply different frequencies in these samples. Bamako represents a useful window into the initial stages of ecological and genomic differentiation from sympatric populations in this important group of malaria vectors.

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.

Genome-Wide Divergence in the West-African Malaria Vector Anopheles melas

Anopheles melas is a member of the recently diverged An. gambiae species complex, a model for speciation studies, and is a locally important malaria vector along the West-African coast where it breeds in brackish water. A recent population genetic study of An. melas revealed species-level genetic differentiation between three population clusters. An. melas West extends from The Gambia to the village of Tiko, Cameroon. The other mainland cluster, An. melas South, extends from the southern Cameroonian village of Ipono to Angola. Bioko Island, Equatorial Guinea An. melas populations are genetically isolated from mainland populations. To examine how genetic differentiation between these An. melas forms is distributed across their genomes, we conducted a genome-wide analysis of genetic differentiation and selection using whole genome sequencing data of pooled individuals (Pool-seq) from a representative population of each cluster. The An. melas forms exhibit high levels of genetic differentiation throughout their genomes, including the presence of numerous fixed differences between clusters. Although the level of divergence between the clusters is on a par with that of other species within the An. gambiae complex, patterns of genome-wide divergence and diversity do not provide evidence for the presence of pre- and/or postmating isolating mechanisms in the form of speciation islands. These results are consistent with an allopatric divergence process with little or no introgression.

An Evolution-Based Screen for Genetic Differentiation between Anopheles Sister Taxa Enriches for Detection of Functional Immune Factors

Nucleotide variation patterns across species are shaped by the processes of natural selection, including exposure to environmental pathogens. We examined patterns of genetic variation in two sister species, Anopheles gambiae and Anopheles coluzzii, both efficient natural vectors of human malaria in West Africa. We used the differentiation signature displayed by a known coordinate selective sweep of immune genes APL1 and TEP1 in A. coluzzii to design a population genetic screen trained on the sweep, classified a panel of 26 potential immune genes for concordance with the signature, and functionally tested their immune phenotypes. The screen results were strongly predictive for genes with protective immune phenotypes: genes meeting the screen criteria were significantly more likely to display a functional phenotype against malaria infection than genes not meeting the criteria (p = 0.0005). Thus, an evolution-based screen can efficiently prioritize candidate genes for labor-intensive downstream functional testing, and safely allow the elimination of genes not meeting the screen criteria. The suite of immune genes with characteristics similar to the APL1-TEP1 selective sweep appears to be more widespread in the A. coluzzii genome than previously recognized. The immune gene differentiation may be a consequence of adaptation of A. coluzzii to new pathogens encountered in its niche expansion during the separation from A. gambiae, although the role, if any of natural selection by Plasmodium is unknown. Application of the screen allowed identification of new functional immune factors, and assignment of new functions to known factors. We describe biochemical binding interactions between immune proteins that underlie functional activity for malaria infection, which highlights the interplay between pathogen specificity and the structure of immune complexes. We also find that most malaria-protective immune factors display phenotypes for either human or rodent malaria, with broad specificity a rarity.

Anopheles gambiae and Anopheles coluzzii are the primary mosquito vectors of human malaria in West Africa. Both of these closely related species efficiently transmit the disease, although they display ecological differences. Previous work showed that A. coluzzii displays distinct genetic patterns in genes important for mosquito immunity. Here, we use this genetic pattern as a filter to examine a panel of potential immune genes, and show that the genetic pattern is strongly predictive for genes that play a functional role in immunity when tested with malaria parasites.

Reticulate Speciation and Barriers to Introgression in the Anopheles gambiae Species Complex

Speciation as a process remains a central focus of evolutionary biology, but our understanding of the genomic architecture and prevalence of speciation in the face of gene flow remains incomplete. The Anopheles gambiae species complex of malaria mosquitoes is a radiation of ecologically diverse taxa. This complex is well-suited for testing for evidence of a speciation continuum and genomic barriers to introgression because its members exhibit partially overlapping geographic distributions as well as varying levels of divergence and reproductive isolation. We sequenced 20 genomes from wild A. gambiae s.s., Anopheles coluzzii, Anopheles arabiensis, and compared these with 12 genomes from the “GOUNDRY” subgroup of A. gambiae s.l. Amidst a backdrop of strong reproductive isolation, we find strong evidence for a speciation continuum with introgression of autosomal chromosomal regions among species and subgroups. The X chromosome, however, is strongly differentiated among all taxa, pointing to a disproportionately large effect of X chromosome genes in driving speciation among anophelines. Strikingly, we find that autosomal introgression has occurred from contemporary hybridization between A. gambiae and A. arabiensis despite strong divergence (∼5× higher than autosomal divergence) and isolation on the X chromosome. In addition to the X, we find strong evidence that lowly recombining autosomal regions, especially pericentromeric regions, serve as barriers to introgression secondarily to the X. We show that speciation with gene flow results in genomic mosaicism of divergence and introgression. Such a reticulate gene pool connecting vector taxa across the speciation continuum has important implications for malaria control efforts.

Does extreme asymmetric dominance promote hybridization between Anopheles coluzzii and Anopheles gambiae s.s. in seasonal malaria mosquito communities of West Africa?

Background

Anopheles gambiae s.s. and An. coluzzii are two of the most important malaria vector species in sub-Saharan Africa. These recently-diverged sibling species do not exhibit intrinsic post-zygotic barriers to reproduction and are thought to be separated by strong assortative mating combined with selection against hybrids. At present, little is known about the ecological conditions that determine hybridization and introgression between these cryptic taxa.

Methods

Swarm segregation and assortative mating were studied in Western Burkina Faso in the villages of Vallée du Kou (VK7) and Soumousso which differed in terms of which sibling species was much rarer than the other, and in Bana where both occurred in similar proportions. Swarms and pairs in copula were collected and genotyped, the proportion of intra and interspecific mating determined, and interspecific sperm transfer checked genetically. Females were collected through larval and adult indoor collections and genotyped or sexed-and-genotyped via a novel multiplex PCR.

Results

A total of 3,687 males and 220 females were collected and genotyped from 109 swarms. Only 3 swarms were composed of males from both species, and these were from the village of VK7 where An. gambiae s.s. was comparatively rare. Mixed-species pairs captured in copula were only detected in that area and made for 3.62 % and 100 % of mating pairs involving An. coluzzii and An. gambiae s.s. individuals, respectively. The high An. gambiae s.s. cross-mating rate was mirrored by high rates of hybridizations estimated from female larvae and adults indoor collections. This contrasted with Soumousso, where despite being much less common than An. gambiae s.s., An. coluzzii males did not form mixed swarms, females were not found in interspecific swarms or copula and hybridization rates were low in both sibling species.

Conclusions

These data suggest that ecological conditions leading to rare An. gambiae s.s. in populations dominated by An. coluzzii may promote a breakdown of spatial swarm segregation and assortative mating between the two species. The lower overall hybridization rates observed at the larval and adult indoor stages compared to cross-mating rates support the idea that post-mating selection processes acting against hybrids may occur mostly prior to and/or at the early larval stages.