Malaria protection due to sickle haemoglobin depends on parasite genotype.

Host genetic factors can confer resistance against malaria1, raising the question of whether this has led to evolutionary adaptation of parasite populations. Here we searched for association between candidate host and parasite genetic variants in 3,346 Gambian and Kenyan children with severe malaria caused by Plasmodium falciparum. We identified a strong association between sickle haemoglobin (HbS) in the host and three regions of the parasite genome, which is not explained by population structure or other covariates, and which is replicated in additional samples. The HbS-associated alleles include nonsynonymous variants in the gene for the acyl-CoA synthetase family member2-4 PfACS8 on chromosome 2, in a second region of chromosome 2, and in a region containing structural variation on chromosome 11. The alleles are in strong linkage disequilibrium and have frequencies that covary with the frequency of HbS across populations, in particular being much more common in Africa than other parts of the world. The estimated protective effect of HbS against severe malaria, as determined by comparison of cases with population controls, varies greatly according to the parasite genotype at these three loci. These findings open up a new avenue of enquiry into the biological and epidemiological significance of the HbS-associated polymorphisms in the parasite genome and the evolutionary forces that have led to their high frequency and strong linkage disequilibrium in African P. falciparum populations.

Detection of simple and complex de novo mutations with multiple reference sequences.

The characterization of de novo mutations in regions of high sequence and structural diversity from whole-genome sequencing data remains highly challenging. Complex structural variants tend to arise in regions of high repetitiveness and low complexity, challenging both de novo assembly, in which short reads do not capture the long-range context required for resolution, and mapping approaches, in which improper alignment of reads to a reference genome that is highly diverged from that of the sample can lead to false or partial calls. Long-read technologies can potentially solve such problems but are currently unfeasible to use at scale. Here we present Corticall, a graph-based method that combines the advantages of multiple technologies and prior data sources to detect arbitrary classes of genetic variant. We construct multisample, colored de Bruijn graphs from short-read data for all samples, align long-read-derived haplotypes and multiple reference data sources to restore graph connectivity information, and call variants using graph path-finding algorithms and a model for simultaneous alignment and recombination. We validate and evaluate the approach using extensive simulations and use it to characterize the rate and spectrum of de novo mutation events in 119 progeny from four Plasmodium falciparum experimental crosses, using long-read data on the parents to inform reconstructions of the progeny and to detect several known and novel nonallelic homologous recombination events.

A forward genetic screen reveals a primary role for Plasmodium falciparum Reticulocyte Binding Protein Homologue 2a and 2b in determining alternative erythrocyte invasion pathways.

Invasion of human erythrocytes is essential for Plasmodium falciparum parasite survival and pathogenesis, and is also a complex phenotype. While some later steps in invasion appear to be invariant and essential, the earlier steps of recognition are controlled by a series of redundant, and only partially understood, receptor-ligand interactions. Reverse genetic analysis of laboratory adapted strains has identified multiple genes that when deleted can alter invasion, but how the relative contributions of each gene translate to the phenotypes of clinical isolates is far from clear. We used a forward genetic approach to identify genes responsible for variable erythrocyte invasion by phenotyping the parents and progeny of previously generated experimental genetic crosses. Linkage analysis using whole genome sequencing data revealed a single major locus was responsible for the majority of phenotypic variation in two invasion pathways. This locus contained the PfRh2a and PfRh2b genes, members of one of the major invasion ligand gene families, but not widely thought to play such a prominent role in specifying invasion phenotypes. Variation in invasion pathways was linked to significant differences in PfRh2a and PfRh2b expression between parasite lines, and their role in specifying alternative invasion was confirmed by CRISPR-Cas9-mediated genome editing. Expansion of the analysis to a large set of clinical P. falciparum isolates revealed common deletions, suggesting that variation at this locus is a major cause of invasion phenotypic variation in the endemic setting. This work has implications for blood-stage vaccine development and will help inform the design and location of future large-scale studies of invasion in clinical isolates.

Cortisol and socioeconomic status in early childhood: A multidimensional assessment.

The hypothalamic-pituitary-adrenal (HPA) axis is sensitive to early life stress, with enduring consequences for biological stress vulnerability and health (Gunnar & Talge, 2008). Low socioeconomic status (SES) is associated with dysregulation of the stress hormone cortisol in early childhood. However, a mechanistic understanding of this association is lacking. Multidimensional assessment of both SES and cortisol is needed to characterize the intricate relations between SES and cortisol function in early childhood. We assessed parent-reported family income, parent education, occupational prestige, neighborhood risk, food insecurity, and household chaos for 12-month-old infants (N = 90) and 3.5-year-old children (N = 91). Hair cortisol concentration (HCC) was obtained from parent and child, indexing chronic biological stress, and diurnal salivary cortisol was measured in the children. Controlling for parent HCC, parent education uniquely predicted infant and child HCC and, in addition, neighborhood risk uniquely predicted infant HCC. Household chaos predicted bedtime salivary cortisol concentration (SCC) for both infants and children, and infant daily cortisol output. Food insecurity was associated with flattened cortisol slope in 3.5-year-old children. Parental sensitivity did not mediate relations between SES and cortisol. Results highlight the utility of SES measures that index unpredictable and unsafe contexts, such as neighborhood risk, food insecurity, and household chaos.

Genomic Analysis of Plasmodium vivax in Southern Ethiopia Reveals Selective Pressures in Multiple Parasite Mechanisms.

The Horn of Africa harbors the largest reservoir of Plasmodium vivax in the continent. Most of sub-Saharan Africa has remained relatively vivax-free due to a high prevalence of the human Duffy-negative trait, but the emergence of strains able to invade Duffy-negative reticulocytes poses a major public health threat. We undertook the first population genomic investigation of P. vivax from the region, comparing the genomes of 24 Ethiopian isolates against data from Southeast Asia to identify important local adaptions. The prevalence of the Duffy binding protein amplification in Ethiopia was 79%, potentially reflecting adaptation to Duffy negativity. There was also evidence of selection in a region upstream of the chloroquine resistance transporter, a putative chloroquine-resistance determinant. Strong signals of selection were observed in genes involved in immune evasion and regulation of gene expression, highlighting the need for a multifaceted intervention approach to combat P. vivax in the region.

Indels, structural variation, and recombination drive genomic diversity in Plasmodium falciparum.

The malaria parasite Plasmodium falciparum has a great capacity for evolutionary adaptation to evade host immunity and develop drug resistance. Current understanding of parasite evolution is impeded by the fact that a large fraction of the genome is either highly repetitive or highly variable and thus difficult to analyze using short-read sequencing technologies. Here, we describe a resource of deep sequencing data on parents and progeny from genetic crosses, which has enabled us to perform the first genome-wide, integrated analysis of SNP, indel and complex polymorphisms, using Mendelian error rates as an indicator of genotypic accuracy. These data reveal that indels are exceptionally abundant, being more common than SNPs and thus the dominant mode of polymorphism within the core genome. We use the high density of SNP and indel markers to analyze patterns of meiotic recombination, confirming a high rate of crossover events and providing the first estimates for the rate of non-crossover events and the length of conversion tracts. We observe several instances of meiotic recombination within copy number variants associated with drug resistance, demonstrating a mechanism whereby fitness costs associated with resistance mutations could be compensated and greater phenotypic plasticity could be acquired.

Analysis of Plasmodium falciparum diversity in natural infections by deep sequencing.

Malaria elimination strategies require surveillance of the parasite population for genetic changes that demand a public health response, such as new forms of drug resistance. Here we describe methods for the large-scale analysis of genetic variation in Plasmodium falciparum by deep sequencing of parasite DNA obtained from the blood of patients with malaria, either directly or after short-term culture. Analysis of 86,158 exonic single nucleotide polymorphisms that passed genotyping quality control in 227 samples from Africa, Asia and Oceania provides genome-wide estimates of allele frequency distribution, population structure and linkage disequilibrium. By comparing the genetic diversity of individual infections with that of the local parasite population, we derive a metric of within-host diversity that is related to the level of inbreeding in the population. An open-access web application has been established for the exploration of regional differences in allele frequency and of highly differentiated loci in the P. falciparum genome.

Genomic Analysis Reveals a Common Breakpoint in Amplifications of the Plasmodium vivax Multidrug Resistance 1 Locus in Thailand.

In regions of coendemicity for Plasmodium falciparum and Plasmodium vivax where mefloquine is used to treat P. falciparum infection, drug pressure mediated by increased copy numbers of the multidrug resistance 1 gene (pvmdr1) may select for mefloquine-resistant P. vivax Surveillance is not undertaken routinely owing in part to methodological challenges in detection of gene amplification. Using genomic data on 88 P. vivax samples from western Thailand, we identified pvmdr1 amplification in 17 isolates, all exhibiting tandem copies of a 37.6-kilobase pair region with identical breakpoints. A novel breakpoint-specific polymerase chain reaction assay was designed to detect the amplification. The assay demonstrated high sensitivity, identifying amplifications in 13 additional, polyclonal infections. Application to 132 further samples identified the common breakpoint in all years tested (2003-2015), with a decline in prevalence after 2012 corresponding to local discontinuation of mefloquine regimens. Assessment of the structure of pvmdr1 amplification in other geographic regions will yield information about the population-specificity of the breakpoints and underlying amplification mechanisms.

K13-propeller polymorphisms in Plasmodium falciparum parasites from sub-Saharan Africa.

Mutations in the Plasmodium falciparum K13-propeller domain have recently been shown to be important determinants of artemisinin resistance in Southeast Asia. This study investigated the prevalence of K13-propeller polymorphisms across sub-Saharan Africa. A total of 1212 P. falciparum samples collected from 12 countries were sequenced. None of the K13-propeller mutations previously reported in Southeast Asia were found, but 22 unique mutations were detected, of which 7 were nonsynonymous. Allele frequencies ranged between 1% and 3%. Three mutations were observed in >1 country, and the A578S was present in parasites from 5 countries. This study provides the baseline prevalence of K13-propeller mutations in sub-Saharan Africa.

Genomic variation during culture adaptation of genetically complex Plasmodium falciparum clinical isolates.

Experimental studies on the biology of malaria parasites have mostly been based on laboratory-adapted lines, but there is limited understanding of how these may differ from parasites in natural infections. Loss-of-function mutants have previously been shown to emerge during culture of some Plasmodium falciparum clinical isolates in analyses focusing on single-genotype infections. The present study included a broader array of isolates, mostly representing multiple-genotype infections, which are more typical in areas where malaria is highly endemic. Genome sequence data from multiple time points over several months of culture adaptation of 28 West African isolates were analysed, including previously available sequences along with new genome sequences from additional isolates and time points. Some genetically complex isolates eventually became fixed over time to single surviving genotypes in culture, whereas others retained diversity, although proportions of genotypes varied over time. Drug resistance allele frequencies did not show overall directional changes, suggesting that resistance-associated costs are not the main causes of fitness differences among parasites in culture. Loss-of-function mutants emerged during culture in several of the multiple-genotype isolates, affecting genes (including AP2-HS, EPAC and SRPK1) for which loss-of-function mutants were previously seen to emerge in single-genotype isolates. Parasite clones were derived by limiting dilution from six of the isolates, and sequencing identified de novo variants not detected in the bulk isolate sequences. Interestingly, several of these were nonsense mutants and frameshifts disrupting the coding sequence of EPAC, the gene with the largest number of independent nonsense mutants previously identified in laboratory-adapted lines. Analysis of genomic identity by descent to explore relatedness among clones revealed co-occurring non-identical sibling parasites, illustrative of the natural genetic structure within endemic populations.

Drug resistance and vaccine target surveillance of Plasmodium falciparum using nanopore sequencing in Ghana.

Malaria results in over 600,000 deaths annually, with the highest burden of deaths in young children living in sub-Saharan Africa. Molecular surveillance can provide important information for malaria control policies, including detection of antimalarial drug resistance. However, genome sequencing capacity in malaria-endemic countries is limited. We designed and implemented an end-to-end workflow to detect Plasmodium falciparum antimalarial resistance markers and diversity in the vaccine target circumsporozoite protein (csp) using nanopore sequencing in Ghana. We analysed 196 clinical samples and showed that our method is rapid, robust, accurate and straightforward to implement. Importantly, our method could be applied to dried blood spot samples, which are readily collected in endemic settings. We report that P. falciparum parasites in Ghana are mostly susceptible to chloroquine, with persistent sulfadoxine-pyrimethamine resistance and no evidence of artemisinin resistance. Multiple single nucleotide polymorphisms were identified in csp, but their significance is uncertain. Our study demonstrates the feasibility of nanopore sequencing for malaria genomic surveillance in endemic countries.

Population genomics reveal distinct and diverging populations of An. minimus in Cambodia.

Anopheles minimus is an important malaria vector throughout its wide geographic range across Southeast Asia. Genome sequencing could provide important insights into the unique malaria transmission dynamics in this region, where many vector species feed and rest outdoors. We describe results from a study using Illumina deep whole-genome sequencing of 302 wild-caught An. minimus collected from three Cambodian provinces over several years (2010, 2014, 2016) and seasons to examine the level of population structure and genetic diversity within this species. These specimens cluster into four distinct populations of An. minimus s.s., with two populations overlapping geographically. We describe the underlying genetic diversity and divergence of these populations and investigated the genetic variation in genes known to be involved in insecticide resistance. We found strong signals of selection within these An. minimus populations, most of which were present in the two Northeastern Cambodian populations and differ from those previously described in African malaria vectors. Cambodia is the focus of the emergence and spread of drug-resistant malaria parasites, so understanding the underlying genetic diversity and resilience of the vectors of these parasites is key to implementing effective malaria control and elimination strategies. These data are publicly available as part of the MalariaGEN Vector Observatory, an open access resource of genome sequence data.

An open dataset of Plasmodium vivax genome variation in 1,895 worldwide samples.

This report describes the MalariaGEN Pv4 dataset, a new release of curated genome variation data on 1,895 samples of Plasmodium vivax collected at 88 worldwide locations between 2001 and 2017. It includes 1,370 new samples contributed by MalariaGEN and VivaxGEN partner studies in addition to previously published samples from these and other sources. We provide genotype calls at over 4.5 million variable positions including over 3 million single nucleotide polymorphisms (SNPs), as well as short indels and tandem duplications. This enlarged dataset highlights major compartments of parasite population structure, with clear differentiation between Africa, Latin America, Oceania, Western Asia and different parts of Southeast Asia. Each sample has been classified for drug resistance to sulfadoxine, pyrimethamine and mefloquine based on known markers at the dhfr, dhps and mdr1 loci. The prevalence of all of these resistance markers was much higher in Southeast Asia and Oceania than elsewhere. This open resource of analysis-ready genome variation data from the MalariaGEN and VivaxGEN networks is driven by our collective goal to advance research into the complex biology of P. vivax and to accelerate genomic surveillance for malaria control and elimination.

An open dataset of Plasmodium falciparum genome variation in 7,000 worldwide samples.

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed.  Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination.

Intrinsic multiplication rate variation and plasticity of human blood stage malaria parasites.

Pathogen multiplication rate is theoretically an important determinant of virulence, although often poorly understood and difficult to measure accurately. We show intrinsic asexual blood stage multiplication rate variation of the major human malaria parasite Plasmodium falciparum to be associated with blood-stage infection intensity in patients. A panel of clinical isolates from a highly endemic West African population was analysed repeatedly during five months of continuous laboratory culture, showing a range of exponential multiplication rates at all timepoints tested, mean rates increasing over time. All isolates had different genome sequences, many containing within-isolate diversity that decreased over time in culture, but increases in multiplication rates were not primarily attributable to genomic selection. New mutants, including premature stop codons emerging in a few isolates, did not attain sufficiently high frequencies to substantially affect overall multiplication rates. Significantly, multiplication rate variation among the isolates at each of the assayed culture timepoints robustly correlated with parasite levels seen in patients at clinical presentation, indicating innate parasite control of multiplication rate that contributes to virulence.

Mapping imported malaria in Bangladesh using parasite genetic and human mobility data.

For countries aiming for malaria elimination, travel of infected individuals between endemic areas undermines local interventions. Quantifying parasite importation has therefore become a priority for national control programs. We analyzed epidemiological surveillance data, travel surveys, parasite genetic data, and anonymized mobile phone data to measure the spatial spread of malaria parasites in southeast Bangladesh. We developed a genetic mixing index to estimate the likelihood of samples being local or imported from parasite genetic data and inferred the direction and intensity of parasite flow between locations using an epidemiological model integrating the travel survey and mobile phone calling data. Our approach indicates that, contrary to dogma, frequent mixing occurs in low transmission regions in the southwest, and elimination will require interventions in addition to reducing imported infections from forested regions. Unlike risk maps generated from clinical case counts alone, therefore, our approach distinguishes areas of frequent importation as well as high transmission.

Evolution and expansion of multidrug-resistant malaria in southeast Asia: a genomic epidemiology study.

BACKGROUND: A multidrug-resistant co-lineage of Plasmodium falciparum malaria, named KEL1/PLA1, spread across Cambodia in 2008-13, causing high rates of treatment failure with the frontline combination therapy dihydroartemisinin-piperaquine. Here, we report on the evolution and spread of KEL1/PLA1 in subsequent years. METHODS: For this genomic epidemiology study, we analysed whole genome sequencing data from P falciparum clinical samples collected from patients with malaria between 2007 and 2018 from Cambodia, Laos, northeastern Thailand, and Vietnam, through the MalariaGEN P falciparum Community Project. Previously unpublished samples were provided by two large-scale multisite projects: the Tracking Artemisinin Resistance Collaboration II (TRAC2) and the Genetic Reconnaissance in the Greater Mekong Subregion (GenRe-Mekong) project. By investigating genome-wide relatedness between parasites, we inferred patterns of shared ancestry in the KEL1/PLA1 population. FINDINGS: We analysed 1673 whole genome sequences that passed quality filters, and determined KEL1/PLA1 status in 1615. Before 2009, KEL1/PLA1 was only found in western Cambodia; by 2016-17 its prevalence had risen to higher than 50% in all of the surveyed countries except for Laos. In northeastern Thailand and Vietnam, KEL1/PLA1 exceeded 80% of the most recent P falciparum parasites. KEL1/PLA1 parasites maintained high genetic relatedness and low diversity, reflecting a recent common origin. Several subgroups of highly related parasites have recently emerged within this co-lineage, with diverse geographical distributions. The three largest of these subgroups (n=84, n=79, and n=47) mostly emerged since 2016 and were all present in Cambodia, Laos, and Vietnam. These expanding subgroups carried new mutations in the crt gene, which arose on a specific genetic background comprising multiple genomic regions. Four newly emerging crt mutations were rare in the early period and became more prevalent by 2016-17 (Thr93Ser, rising to 19·8%; His97Tyr to 11·2%; Phe145Ile to 5·5%; and Ile218Phe to 11·1%). INTERPRETATION: After emerging and circulating for several years within Cambodia, the P falciparum KEL1/PLA1 co-lineage diversified into multiple subgroups and acquired new genetic features, including novel crt mutations. These subgroups have rapidly spread into neighbouring countries, suggesting enhanced fitness. These findings highlight the urgent need for elimination of this increasingly drug-resistant parasite co-lineage, and the importance of genetic surveillance in accelerating malaria elimination efforts. FUNDING: Wellcome Trust, Bill & Melinda Gates Foundation, UK Medical Research Council, and UK Department for International Development.

Origins of the current outbreak of multidrug-resistant malaria in southeast Asia: a retrospective genetic study.

BACKGROUND: Antimalarial resistance is rapidly spreading across parts of southeast Asia where dihydroartemisinin-piperaquine is used as first-line treatment for Plasmodium falciparum malaria. The first published reports about resistance to antimalarial drugs came from western Cambodia in 2013. Here, we analyse genetic changes in the P falciparum population of western Cambodia in the 6 years before those reports. METHODS: We analysed genome sequence data on 1492 P falciparum samples from 11 locations across southeast Asia, including 464 samples collected in western Cambodia between 2007 and 2013. Different epidemiological origins of resistance were identified by haplotypic analysis of the kelch13 artemisinin resistance locus and the plasmepsin 2-3 piperaquine resistance locus. FINDINGS: We identified more than 30 independent origins of artemisinin resistance, of which the KEL1 lineage accounted for 140 (91%) of 154 parasites resistant to dihydroartemisinin-piperaquine. In 2008, KEL1 combined with PLA1, the major lineage associated with piperaquine resistance. By 2013, the KEL1/PLA1 co-lineage had reached a frequency of 63% (24/38) in western Cambodia and had spread to northern Cambodia. INTERPRETATION: The KEL1/PLA1 co-lineage emerged in the same year that dihydroartemisinin-piperaquine became the first-line antimalarial drug in western Cambodia and spread rapidly thereafter, displacing other artemisinin-resistant parasite lineages. These findings have important implications for management of the global health risk associated with the current outbreak of multidrug-resistant malaria in southeast Asia. FUNDING: Wellcome Trust, Bill & Melinda Gates Foundation, Medical Research Council, UK Department for International Development, and the Intramural Research Program of the National Institute of Allergy and Infectious Diseases.

Genetic markers associated with dihydroartemisinin-piperaquine failure in Plasmodium falciparum malaria in Cambodia: a genotype-phenotype association study.

BACKGROUND: As the prevalence of artemisinin-resistant Plasmodium falciparum malaria increases in the Greater Mekong subregion, emerging resistance to partner drugs in artemisinin combination therapies seriously threatens global efforts to treat and eliminate this disease. Molecular markers that predict failure of artemisinin combination therapy are urgently needed to monitor the spread of partner drug resistance, and to recommend alternative treatments in southeast Asia and beyond. METHODS: We did a genome-wide association study of 297 P falciparum isolates from Cambodia to investigate the relationship of 11 630 exonic single-nucleotide polymorphisms (SNPs) and 43 copy number variations (CNVs) with in-vitro piperaquine 50% inhibitory concentrations (IC50s), and tested whether these genetic variants are markers of treatment failure with dihydroartemisinin-piperaquine. We then did a survival analysis of 133 patients to determine whether candidate molecular markers predicted parasite recrudescence following dihydroartemisinin-piperaquine treatment. FINDINGS: Piperaquine IC50s increased significantly from 2011 to 2013 in three Cambodian provinces (2011 vs 2013 median IC50s: 20·0 nmol/L [IQR 13·7-29·0] vs 39·2 nmol/L [32·8-48·1] for Ratanakiri, 19·3 nmol/L [15·1-26·2] vs 66·2 nmol/L [49·9-83·0] for Preah Vihear, and 19·6 nmol/L [11·9-33·9] vs 81·1 nmol/L [61·3-113·1] for Pursat; all p≤10-3; Kruskal-Wallis test). Genome-wide analysis of SNPs identified a chromosome 13 region that associates with raised piperaquine IC50s. A non-synonymous SNP (encoding a Glu415Gly substitution) in this region, within a gene encoding an exonuclease, associates with parasite recrudescence following dihydroartemisinin-piperaquine treatment. Genome-wide analysis of CNVs revealed that a single copy of the mdr1 gene on chromosome 5 and a novel amplification of the plasmepsin 2 and plasmepsin 3 genes on chromosome 14 also associate with raised piperaquine IC50s. After adjusting for covariates, both exo-E415G and plasmepsin 2-3 markers significantly associate (p=3·0 × 10-8 and p=1·7 × 10-7, respectively) with decreased treatment efficacy (survival rates 0·38 [95% CI 0·25-0·51] and 0·41 [0·28-0·53], respectively). INTERPRETATION: The exo-E415G SNP and plasmepsin 2-3 amplification are markers of piperaquine resistance and dihydroartemisinin-piperaquine failures in Cambodia, and can help monitor the spread of these phenotypes into other countries of the Greater Mekong subregion, and elucidate the mechanism of piperaquine resistance. Since plasmepsins are involved in the parasite's haemoglobin-to-haemozoin conversion pathway, targeted by related antimalarials, plasmepsin 2-3 amplification probably mediates piperaquine resistance. FUNDING: Intramural Research Program of the US National Institute of Allergy and Infectious Diseases, National Institutes of Health, Wellcome Trust, Bill & Melinda Gates Foundation, Medical Research Council, and UK Department for International Development.

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.