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

“Far-West” Africa is known to be a secondary contact zone between the two major malaria vectors Anopheles coluzzii and A. gambiae. We investigated gene-flow and potentially adaptive introgression between these species along a west-to-east transect in Guinea Bissau, the putative core of this hybrid zone. To evaluate the extent and direction of gene flow, we genotyped site 702 in Intron-1 of the para Voltage-Gated Sodium Channel gene, a species-diagnostic nucleotide position throughout most of A. coluzzii and A. gambiae sympatric range. We also analyzed polymorphism in the thioester-binding domain (TED) of the innate immunity-linked thioester-containing protein 1 (TEP1) to investigate whether elevated hybridization might facilitate the exchange of variants linked to adaptive immunity and Plasmodium refractoriness. Our results confirm asymmetric introgression of genetic material from A. coluzzii to A. gambiae and disruption of linkage between the centromeric "genomic islands" of inter-specific divergence. We report that A. gambiae from the Guinean hybrid zone possesses an introgressed TEP1 resistant allelic class, found exclusively in A. coluzzii elsewhere and apparently swept to fixation in West Africa (i.e. Mali and Burkina Faso). However, no detectable fixation of this allele was found in Guinea Bissau, which may suggest that ecological pressures driving segregation between the two species in larval habitats in this region may be different from those experienced in northern and more arid parts of the species’ range. Finally, our results also suggest a genetic subdivision between coastal and inland A. gambiae Guinean populations and provide clues on the importance of ecological factors in intra-specific differentiation processes.     

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.     

STRONG ASSORTATIVE MATING BY DIET,COLOR, SIZE,AND MORPHOLOGY BUT LIMITED PROGRESS TOWARD SYMPATRIC SPECIATION IN A CLASSIC EXAMPLE: CAMEROON CRATER LAKE CICHLIDS

Models predict that sympatric speciation depends on restrictive parameter ranges, such as sufficiently strong disruptive selection and assortative mating, but compelling examples in nature have rarely been used to test these predictions. I measured the strength of assortative mating within a species complex of Tilapia in Lake Ejagham, Cameroon, a celebrated example of incipient sympatric adaptive radiation. This species complex is in the earliest stages of speciation: morphological and ecological divergence are incomplete, species differ primarily in breeding coloration, and introgression is common. I captured 27 mated pairs in situ and measured the diet, color, size, and morphology of each individual. I found strong assortative mating by color, size, head depth, and dietary source of benthic or pelagic prey along two independent dimensions of assortment. Thus, Ejagham Tilapia showed strong assortative mating most conducive to sympatric speciation. Nonetheless, in contrast to a morphologically bimodal Sarotherodon cichlid species pair in the lake, Ejagham Tilapia show more limited progress toward speciation, likely due to insufficient strength of disruptive selection on morphology estimated in a previous study (γ = 0.16). This supports the predicted dependence of sympatric speciation on strong assortment and strong disruptive selection by examining a potentially stalled example in nature.     

昆虫同域物种形成机制及检验

从生物学、生态和遗传的角度阐述昆虫同域物种形成过程中涉及到的可能性机制。昆虫同域种的分化与作用于同域初始种群的歧化选择密切相关,歧化选择间接导致种群生态特征和遗传特征的分化,促进同域近缘种群间的生殖隔离。同域物种形成的过程中涉及到性状替换、性选择、同型交配等机制。寄主专化型多见于昆虫同域种的分化过程中,一般以植食性昆虫为主。有关昆虫同域物种形成的检验机制有多种,归纳起来主要包括同型交配的检验、遗传漂流的量化、遗传分化程度和连锁不平衡(LD)的检测、杂交后代适合度的估算等。目前发现在许多昆虫种类中存在同域物种形成的可能性,但是有关其隔离机制并没有得到充分的解释。     

Resistance to pirimiphos-methyl in West African Anopheles is spreading via duplication and introgression of the Ace1 locus

Vector population control using insecticides is a key element of current strategies to prevent malaria transmission in Africa. The introduction of effective insecticides, such as the organophosphate pirimiphos-methyl, is essential to overcome the recurrent emergence of resistance driven by the highly diverse Anopheles genomes. Here, we use a population genomic approach to investigate the basis of pirimiphos-methyl resistance in the major malaria vectors Anopheles gambiae and A. coluzzii. A combination of copy number variation and a single non-synonymous substitution in the acetylcholinesterase gene, Ace1, provides the key resistance diagnostic in an A. coluzzii population from Côte d’Ivoire that we used for sequence-based association mapping, with replication in other West African populations. The Ace1 substitution and duplications occur on a unique resistance haplotype that evolved in A. gambiae and introgressed into A. coluzzii, and is now common in West Africa primarily due to selection imposed by other organophosphate or carbamate insecticides. Our findings highlight the predictive value of this complex resistance haplotype for phenotypic resistance and clarify its evolutionary history, providing tools to for molecular surveillance of the current and future effectiveness of pirimiphos-methyl based interventions.     

Conditions for sympatric speciation: A diploid model incorporating habitat fidelity and non-habitat assortative mating

Summary Three types of genes have been proposed to promote sympatric speciation: habitat preference genes, assortative mating genes and habitat-based fitness genes. Previous computer models have analysed these genes separately or in pairs. In this paper we describe a multilocus model in which genes of all three types are considered simultaneously. Our computer simulations show that speciation occurs in complete sympatry under a broad range of conditions. The process includes an initial diversification phase during which a slight amount of divergence occurs, a quasi-equilibrium phase of stasis during which little or no detectable divergence occurs and a completion phase during which divergence is dramatic and gene flow between diverging habitat morphs is rapidly eliminated. Habitat preference genes and habitat-specific fitness genes become associated when assortative mating occurs due to habitat preference, but interbreeding between individuals adapted to different habitats occurs unless habitat preference is almost error free. However, nonhabitat assortative mating, when coupled with habitat preference can eliminate this interbreeding. Even when several loci contribute to the probability of expression of non-habitat assortative mating and the contributions of individual loci are small, gene flow between diverging portions of the population can terminate within less than 1000 generations.     

Incipient speciation driven by hypertrophied lips in Midas cichlid fishes?

Sympatric speciation has been debated in evolutionary biology for decades. Although it has gained in acceptance recently, still only a handful of empirical examples are seen as valid (e.g. crater lake cichlids). In this study, we disentangle the role of hypertrophied lips in the repeated adaptive radiations of Nicaraguan crater lake cichlid fish. We assessed the role of disruptive selection and assortative mating during the early stages of divergence and found a functional trade‐off in feeding behaviour between thick‐ and thin‐lipped ecotypes, suggesting that this trait is a target of disruptive selection. Thick‐lipped fish perform better on nonevasive prey at the cost of a poorer performance on evasive prey. Using enclosures in the wild, we found that thick‐lipped fish perform significantly better in rocky than in sandy habitats. We found almost no mixed pairs during two breeding seasons and hence significant assortative mating. Genetic differentiation between ecotypes seems to be related to the time since colonization, being subtle in L. Masaya (1600 generations ago) and absent in the younger L. Apoyeque (<600 generations ago). Genome‐wide differentiation between ecotypes was higher in the old source lakes than in the young crater lakes. Our results suggest that hypertrophied lips might be promoting incipient sympatric speciation through divergent selection (ecological divergence in feeding performance) and nonrandom mating (assortative mating) in the young Nicaraguan crater lakes. Nonetheless, further manipulative experiments are needed in order to confirm the role of hypertrophied lips as the main cue for assortative mating.     

Reproductive Isolation among Sympatric Molecular Forms of An. gambiae from Inland Areas of South-Eastern Senegal

The Anopheles gambiae species complex includes at least seven morphologically indistinguishable species, one of which, Anopheles gambiae sensu stricto, is the primary mosquito vector responsible for the transmission of malaria across sub-Saharan Africa. Sympatric ecological diversification of An. gambiae s.s. is in progress within this complex, leading to the emergence of at least two incipient species (the M and S molecular forms now recognized as good species and named An. coluzzii and An. gambiae respectively) that show heterogeneous levels of divergence in most parts of Africa. However, this process seems to have broken down in coastal areas of West Africa at the extreme edge of the distribution. We undertook a longitudinal study to describe An. gambiae s.s. populations collected from two inland transects with different ecological characteristics in south-eastern Senegal. Analysis of samples collected from 20 sites across these two transects showed the M and S molecular forms coexisted at almost all sampled sites. Overall, similar hybridization rates (2.16% and 1.86%) were recorded in the two transects; sites with relatively high frequencies of M/S hybrids (up to 7%) were clustered toward the north-western part of both transects, often near urban settings. Estimated inbreeding indices for this putative speciation event varied spatially (range: 0.52–1), with hybridization rates being generally lower than expected under panmictic conditions. Such observations suggest substantial reproductive isolation between the M and S molecular forms, and further support the ongoing process of speciation in these inland areas. According to a recent reclassification of the An. gambiae complex, the M and S molecular forms from this zone correspond to An. coluzzii and An. gambiae, respectively. There is considerable evidence that these molecular forms differ in their behavioural and ecological characteristics. Detailed study of these characteristics will allow the development and implementation of better insect control strategies for combating malaria.     

Waiting for sympatric speciation

While it now appears likely that sympatric speciation is possible, its generality remains contentious. If it really is rare, then most natural populations must not fit the assumptions of sympatric speciation theory. A better understanding of these assumptions may help identify when sympatric speciation is or is not likely. This paper investigates two such assumptions: that genetic variation for stringent assortative mating is not limiting and that females are not penalized for mating assortatively. Simulations demonstrate that the speed of sympatric speciation is very sensitive to the population's capacity for stringent assortative mating and is potentially extremely slow. The rapid divergence often thought to be a hallmark of sympatric speciation may only occur in a restricted area of parameter space.     

Prominent intraspecific genetic divergence within Anopheles gambiae sibling species triggered by habitat discontinuities across a riverine landscape

The Anopheles gambiae complex of mosquitoes includes malaria vectors at different stages of speciation, whose study enables a better understanding of how adaptation to divergent environmental conditions leads to evolution of reproductive isolation. We investigated the population genetic structure of closely related sympatric taxa that have recently been proposed as separate species (An. coluzzii and An. gambiae), sampled from diverse habitats along the Gambia river in West Africa. We characterized putatively neutral microsatellite loci as well as chromosomal inversion polymorphisms known to be associated with ecological adaptation. The results revealed strong ecologically associated population subdivisions within both species. Microsatellite loci on chromosome‐3L revealed clear differentiation between coastal and inland populations, which in An. coluzzii is reinforced by a unusual inversion polymorphism pattern, supporting the hypothesis of genetic divergence driven by adaptation to the coastal habitat. A strong reduction of gene flow was observed between An. gambiae populations west and east of an extensively rice‐cultivated region apparently colonized exclusively by An. coluzzii. Notably, this ‘intraspecific’ differentiation is higher than that observed between the two species and involves also the centromeric region of chromosome‐X which has previously been considered a marker of speciation within this complex, possibly suggesting that the two populations may be at an advanced stage of differentiation triggered by human‐made habitat fragmentation. These results confirm ongoing ecological speciation within these most important Afro‐tropical malaria vectors and raise new questions on the possible effect of this process in malaria transmission.     

Genetic Structure of a Local Population of the Anopheles gambiae Complex in Burkina Faso

Members of the Anopheles gambiae species complex are primary vectors of human malaria in Africa. Population heterogeneities for ecological and behavioral attributes expand and stabilize malaria transmission over space and time, and populations may change in response to vector control, urbanization and other factors. There is a need for approaches to comprehensively describe the structure and characteristics of a sympatric local mosquito population, because incomplete knowledge of vector population composition may hinder control efforts. To this end, we used a genome-wide custom SNP typing array to analyze a population collection from a single geographic region in West Africa. The combination of sample depth (n = 456) and marker density (n = 1536) unambiguously resolved population subgroups, which were also compared for their relative susceptibility to natural genotypes of Plasmodium falciparum malaria. The population subgroups display fluctuating patterns of differentiation or sharing across the genome. Analysis of linkage disequilibrium identified 19 new candidate genes for association with underlying population divergence between sister taxa, A. coluzzii (M-form) and A. gambiae (S-form).     

Positive assortative mating between recently described sympatric morphs of Icelandic sticklebacks

Recently, models of sympatric speciation have suggested that assortative mating can develop between sympatric morphs due to divergence in an ecologically important character. For example, in sympatric pairs of threespine stickleback (Gasterosteus aculeatus L.) size-assortative mating seems to be instrumental in reproductive isolation. Here, we examine courtship behaviour and assortative mating of newly described sympatric stickleback morphs in Lake Thingvallavatn, Iceland. We find that the two morphs show strong positive assortative mating. However, the mechanism involved in mate choice does not seem to be as straightforward as in other similar systems of sympatric stickleback morphs and may involve variation in nest type.     

Evolution of assortative mating following selective introgression of pigmentation genes between two Drosophila species

Adaptive introgression is ubiquitous in animals, but experimental support for its role in driving speciation remains scarce. In the absence of conscious selection, admixed laboratory strains of Drosophila asymmetrically and progressively lose alleles from one parental species and reproductive isolation against the predominant parent ceases after 10 generations. Here, we selectively introgressed during 1 year light pigmentation genes of D. santomea into the genome of its dark sibling D. yakuba, and vice versa. We found that the pace of phenotypic change differed between the species and the sexes and identified through genome sequencing common as well as distinct introgressed loci in each species. Mating assays showed that assortative mating between introgressed flies and both parental species persisted even after 4 years (~60 generations) from the end of the selection. Those results indicate that selective introgression of as low as 0.5% of the genome can beget morphologically distinct and reproductively isolated strains, two prerequisites for the delimitation of new species. Our findings hence represent a significant step toward understanding the genome‐wide dynamics of speciation‐through‐introgression.     

Multi-species outcomes in a common model of sympatric speciation

While models of sympatric speciation are motivated in part by multi-species adaptive radiations such as the Cameroon crater lake cichlids, existing models have focused on bifurcation into a single pair of daughter species. This paper shows that a familiar model of sympatric speciation, driven by intraspecific competition and assortative mating based on ecological characters values, can yield multiple daughter species if individual niche widths are sufficiently restricted. Surprisingly, the multi-species outcome is not produced by successive bifurcation events, but by simultaneous divergence resulting in a hard polytomy. This result is sensitive to a number of assumptions, whose violation may prevent speciation. In some cases when speciation fails, the population instead ends in a state that closely resembles incipient species pairs, with an ecological polymorphism and partial reproductive isolation. However, this polymorphism is stable and does not lead to complete reproductive isolation, suggesting that empirical cases of incipient species pairs may not always end in speciation.     

Demography and selection analysis of the incipient adaptive radiation of a Hawaiian woody species

Ecological divergence in a species provides a valuable opportunity to study the early stages of speciation. We focused on Metrosideros polymorpha, a unique example of the incipient radiation of woody species, to examine how an ecological divergence continues in the face of gene flow. We analyzed the whole genomes of 70 plants collected throughout the island of Hawaii, which is the youngest island with the highest altitude in the archipelago and encompasses a wide range of environments. The continuous M. polymorpha forest stands on the island of Hawaii were differentiated into three genetic clusters, each of which grows in a distinctive environment and includes substantial genetic and phenotypic diversity. The three genetic clusters showed signatures of selection in genomic regions encompassing genes relevant to environmental adaptations, including genes associated with light utilization, oxidative stress, and leaf senescence, which are likely associated with the ecological differentiation of the species. Our demographic modeling suggested that the glaberrima cluster in wet environments maintained a relatively large population size and two clusters split: polymorpha in the subalpine zone and incana in dry and hot conditions. This ecological divergence possibly began before the species colonized the island of Hawaii. Interestingly, the three clusters recovered genetic connectivity coincidentally with a recent population bottleneck, in line with the weak reproductive isolation observed in the species. This study highlights that the degree of genetic differentiation between ecologically-diverged populations can vary depending on the strength of natural selection in the very early phases of speciation.     

Complex histories of repeated gene flow in Cameroon crater lake cichlids cast doubt on one of the clearest examples of sympatric speciation

One of the most celebrated examples of sympatric speciation in nature are monophyletic radiations of cichlid fishes endemic to Cameroon crater lakes. However, phylogenetic inference of monophyly may not detect complex colonization histories involving some allopatric isolation, such as double invasions obscured by genome‐wide gene flow. Population genomic approaches are better suited to test hypotheses of sympatric speciation in these cases. Here, we use comprehensive sampling from all four sympatric crater lake cichlid radiations in Cameroon and outgroups across Africa combined with next‐generation sequencing to genotype tens of thousands of SNPs. We find considerable evidence of gene flow between all four radiations and neighboring riverine populations after initial colonization. In a few cases, some sympatric species are more closely related to outgroups than others, consistent with secondary gene flow facilitating their speciation. Our results do not rule out sympatric speciation in Cameroon cichlids, but rather reveal a complex history of speciation with gene flow, including allopatric and sympatric phases, resulting in both reproductively isolated species and incipient species complexes. The best remaining non‐cichlid examples of sympatric speciation all involve assortative mating within microhabitats. We speculate that this feature may be necessary to complete the process of sympatric speciation in nature.     

Evidence for divergent selection between the molecular forms of <Emphasis Type="Italic">Anopheles gambiae</Emphasis>: role of predation

Background  

The molecular forms of Anopheles gambiae are undergoing speciation. They are characterized by a strong assortative mating and they display partial habitat segregation. The M form is mostly found in flooded/irrigated areas whereas the S form dominates in the surrounding areas, but the ecological factors that shape this habitat segregation are not known. Resource competition has been demonstrated between species undergoing divergent selection, but resource competition is not the only factor that can lead to divergence.     

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.     

Chromosome Inversions,Genomic Differentiation and Speciation in the African Malaria Mosquito Anopheles gambiae

The African malaria vector, Anopheles gambiae, is characterized by multiple polymorphic chromosomal inversions and has become widely studied as a system for exploring models of speciation. Near complete reproductive isolation between different inversion types, known as chromosomal forms, has led to the suggestion that A. gambiae is in early stages of speciation, with divergence evolving in the face of considerable gene flow. We compared the standard chromosomal arrangement (Savanna form) with genomes homozygous for j, b, c, and u inversions (Bamako form) in order to identify regions of genomic divergence with respect to inversion polymorphism. We found levels of divergence between the two sub-taxa within some of these inversions (2Rj and 2Rb), but at a level lower than expected and confined near the inversion breakpoints, consistent with a gene flux model. Unexpectedly, we found that the majority of diverged regions were located on the X chromosome, which contained half of all significantly diverged regions, with much of this divergence located within exons. This is surprising given that the Bamako and Savanna chromosomal forms are both within the S molecular form that is defined by a locus near centromere of X chromosome. Two X-linked genes (a heat shock protein and P450 encoding genes) involved in reproductive isolation between the M and S molecular forms of A. gambiae were also significantly diverged between the two chromosomal forms. These results suggest that genes mediating reproductive isolation are likely located on the X chromosome, as is thought to be the case for the M and S molecular forms. We conclude that genes located on the sex chromosome may be the major force driving speciation between these chromosomal forms of A. gambiae.     

Sexual signal evolution and patterns of assortative mating across an intraspecific contact zone

Contact zones provide important insights into the evolutionary processes that underlie lineage divergence and speciation. Here, we use a contact zone to ascertain speciation potential in the red-eyed treefrog (Agalychnis callidryas), a brightly coloured and polymorphic frog that exhibits unusually high levels of intraspecific variation. Populations of A. callidryas differ in a number of traits, several of which are known sexual signals that mediate premating reproductive isolation in allopatric populations. Along the Caribbean coast of Costa Rica, a ~100 km contact zone, situated between two phenotypically and genetically divergent parent populations, contains multiple colour pattern phenotypes and late-generation hybrids. This contact zone provides the opportunity to examine processes that are important in the earliest stages of lineage divergence. We performed analyses of colour pattern variation in five contact zone sites and six parental sites and found complex, continuous colour variation along the contact zone. We found discordance between the geographic distribution of colour pattern and previously described genomic population structure. We then used a parental site and contact zone site to measure assortative mating and directional selection from naturally-occurring amplectant mating pairs. We found assortative mating in a parental population, but no assortative mating in the contact zone. Furthermore, we uncovered evidence of directional preference towards the adjacent parental phenotype in the contact zone population, but no directional preference in the parent population. Combined, these data provide insights into potential dynamics at the contact zone borders and indicate that incipient speciation between parent populations will be slowed.