Antigenic strain diversity predicts different biogeographic patterns of maintenance and decline of antimalarial drug resistance.

The establishment and spread of antimalarial drug resistance vary drastically across different biogeographic regions. Though most infections occur in sub-Saharan Africa, resistant strains often emerge in low-transmission regions. Existing models on resistance evolution lack consensus on the relationship between transmission intensity and drug resistance, possibly due to overlooking the feedback between antigenic diversity, host immunity, and selection for resistance. To address this, we developed a novel compartmental model that tracks sensitive and resistant parasite strains, as well as the host dynamics of generalized and antigen-specific immunity. Our results show a negative correlation between parasite prevalence and resistance frequency, regardless of resistance cost or efficacy. Validation using chloroquine-resistant marker data supports this trend. Post discontinuation of drugs, resistance remains high in low-diversity, low-transmission regions, while it steadily decreases in high-diversity, high-transmission regions. Our study underscores the critical role of malaria strain diversity in the biogeographic patterns of resistance evolution.

Drug resistance among strains of the parasites that cause malaria is a growing problem for people relying on antimalarial drugs to protect them from the disease. This phenomenon is global yet exactly how resistance emerges, spreads and persists in a population often differs greatly between regions, which can complicate malaria control projects. For example, discontinuing the use of antimalarials can lead to the frequency of resistant strains declining in an area, such as Africa, but persisting at high levels in others, including Asia and South America. Gaining resistance often leads to parasites becoming less transmissible than other strains. When antimalarials are not used, sensitive strains usually outcompete their resistant counterparts. However, prolonged use of antimalarial drugs tends to eliminate susceptible strains, allowing the previously outcompeted resistant strains to dominate. The local dynamics of antimalarial resistance are also shaped by multiple other factors such as transmission levels (how common the disease is in the region), the type of antimalarial measures used (such as drugs and mosquito nets), or previous immunity the population may have developed to specific strains. While many computational models have been developed to capture these dynamics, they usually fail to include strain diversity ­ a parameter reflecting the number of malaria strains the immune system is exposed to. This parameter is important as parasites need to escape both host immunity and drugs in order to be successful. To address this gap, He, Chaillet, and Labbé created a computational model to investigate how strain diversity, transmission levels and other related factors influence antimalarial resistance. The model was used to explore how the frequency of resistant and susceptible strains changes over time once antimalarial drugs are rolled out and then halted. These analyses show that in areas with both low strain diversity and low transmission levels, susceptible parasites are more likely to be wiped out from the population, leading to a high frequency of resistant strains that persist after drugs are discontinued. However, in high diversity and high transmission regions, susceptible strains can remain in the population. Therefore, when drug treatments are stopped, resistance levels are more likely to drop due to these parasites outcompeting the drug-resistant ones. Overall, this work demonstrates how modelling approaches that include strain diversity can help inform public health decisions aimed at reducing antimalarial resistance. In particular, they can provide important insights into the control strategies that are best suited for a specific region, suggesting that in low transmission areas intensive drug treatment may contribute to resistance. Instead, preventative strategies such as eliminating mosquitos and preventing bites with bed nets may prove more beneficial at reducing transmission rates in such areas.

Comparison of molecular surveillance methods to assess changes in the population genetics of Plasmodium falciparum in high transmission.

A major motivation for developing molecular methods for malaria surveillance is to measure the impact of control interventions on the population genetics of Plasmodium falciparum as a potential marker of progress towards elimination. Here we assess three established methods (i) single nucleotide polymorphism (SNP) barcoding (panel of 24-biallelic loci), (ii) microsatellite genotyping (panel of 12-multiallelic loci), and (iii) varcoding (fingerprinting var gene diversity, akin to microhaplotyping) to identify changes in parasite population genetics in response to a short-term indoor residual spraying (IRS) intervention. Typical of high seasonal transmission in Africa, multiclonal infections were found in 82.3% (median 3; range 1-18) and 57.8% (median 2; range 1-12) of asymptomatic individuals pre- and post-IRS, respectively, in Bongo District, Ghana. Since directly phasing multilocus haplotypes for population genetic analysis is not possible for biallelic SNPs and microsatellites, we chose ~200 low-complexity infections biased to single and double clone infections for analysis. Each genotyping method presented a different pattern of change in diversity and population structure as a consequence of variability in usable data and the relative polymorphism of the molecular markers (i.e., SNPs < microsatellites < var). Varcoding and microsatellite genotyping showed the overall failure of the IRS intervention to significantly change the population structure from pre-IRS characteristics (i.e., many diverse genomes of low genetic similarity). The 24-SNP barcode provided limited information for analysis, largely due to the biallelic nature of SNPs leading to a high proportion of double-allele calls and a view of more isolate relatedness compared to microsatellites and varcoding. Relative performance, suitability, and cost-effectiveness of the methods relevant to sample size and local malaria elimination in high-transmission endemic areas are discussed.

Antigenic strain diversity predicts different biogeographic patterns of maintenance and decline of anti-malarial drug resistance.

The establishment and spread of anti-malarial drug resistance vary drastically across different biogeographic regions. Though most infections occur in Sub-Saharan Africa, resistant strains often emerge in low-transmission regions. Existing models on resistance evolution lack consensus on the relationship between transmission intensity and drug resistance, possibly due to overlooking the feedback between antigenic diversity, host immunity, and selection for resistance. To address this, we developed a novel compartmental model that tracks sensitive and resistant parasite strains, as well as the host dynamics of generalized and antigen-specific immunity. Our results show a negative correlation between parasite prevalence and resistance frequency, regardless of resistance cost or efficacy. Validation using chloroquine-resistant marker data supports this trend. Post discontinuation of drugs, resistance remains high in low-diversity, low-transmission regions, while it steadily decreases in high-diversity, high-transmission regions. Our study underscores the critical role of malaria strain diversity in the biogeographic patterns of resistance evolution.

Performance of SNP barcodes to determine genetic diversity and population structure of Plasmodium falciparum in Africa.

Panels of informative biallelic single nucleotide polymorphisms (SNPs) have been proposed to be an economical method to fast-track the population genetic analysis of Plasmodium falciparum in malaria-endemic areas. Whilst used successfully in low-transmission areas where infections are monoclonal and highly related, we present the first study to evaluate the performance of these 24- and 96-SNP molecular barcodes in African countries, characterised by moderate-to-high transmission, where multiclonal infections are prevalent. For SNP barcodes it is generally recommended that the SNPs chosen i) are biallelic, ii) have a minor allele frequency greater than 0.10, and iii) are independently segregating, to minimise bias in the analysis of genetic diversity and population structure. Further, to be standardised and used in many population genetic studies, these barcodes should maintain characteristics i) to iii) across various iv) geographies and v) time points. Using haplotypes generated from the MalariaGEN P. falciparum Community Project version six database, we investigated the ability of these two barcodes to fulfil these criteria in moderate-to-high transmission African populations in 25 sites across 10 countries. Predominantly clinical infections were analysed, with 52.3% found to be multiclonal, generating high proportions of mixed-allele calls (MACs) per isolate thereby impeding haplotype construction. Of the 24- and 96-SNPs, loci were removed if they were not biallelic and had low minor allele frequencies in all study populations, resulting in 20- and 75-SNP barcodes respectively for downstream population genetics analysis. Both SNP barcodes had low expected heterozygosity estimates in these African settings and consequently biased analyses of similarity. Both minor and major allele frequencies were temporally unstable. These SNP barcodes were also shown to identify weak genetic differentiation across large geographic distances based on Mantel Test and DAPC. These results demonstrate that these SNP barcodes are vulnerable to ascertainment bias and as such cannot be used as a standardised approach for malaria surveillance in moderate-to-high transmission areas in Africa, where the greatest genomic diversity of P. falciparum exists at local, regional and country levels.

Genomic analysis of two phlebotomine sand fly vectors of Leishmania from the New and Old World.

Phlebotomine sand flies are of global significance as important vectors of human disease, transmitting bacterial, viral, and protozoan pathogens, including the kinetoplastid parasites of the genus Leishmania, the causative agents of devastating diseases collectively termed leishmaniasis. More than 40 pathogenic Leishmania species are transmitted to humans by approximately 35 sand fly species in 98 countries with hundreds of millions of people at risk around the world. No approved efficacious vaccine exists for leishmaniasis and available therapeutic drugs are either toxic and/or expensive, or the parasites are becoming resistant to the more recently developed drugs. Therefore, sand fly and/or reservoir control are currently the most effective strategies to break transmission. To better understand the biology of sand flies, including the mechanisms involved in their vectorial capacity, insecticide resistance, and population structures we sequenced the genomes of two geographically widespread and important sand fly vector species: Phlebotomus papatasi, a vector of Leishmania parasites that cause cutaneous leishmaniasis, (distributed in Europe, the Middle East and North Africa) and Lutzomyia longipalpis, a vector of Leishmania parasites that cause visceral leishmaniasis (distributed across Central and South America). We categorized and curated genes involved in processes important to their roles as disease vectors, including chemosensation, blood feeding, circadian rhythm, immunity, and detoxification, as well as mobile genetic elements. We also defined gene orthology and observed micro-synteny among the genomes. Finally, we present the genetic diversity and population structure of these species in their respective geographical areas. These genomes will be a foundation on which to base future efforts to prevent vector-borne transmission of Leishmania parasites.

Neutral vs. non-neutral genetic footprints of Plasmodium falciparum multiclonal infections.

At a time when effective tools for monitoring malaria control and eradication efforts are crucial, the increasing availability of molecular data motivates their application to epidemiology. The multiplicity of infection (MOI), defined as the number of genetically distinct parasite strains co-infecting a host, is one key epidemiological parameter for evaluating malaria interventions. Estimating MOI remains a challenge for high-transmission settings where individuals typically carry multiple co-occurring infections. Several quantitative approaches have been developed to estimate MOI, including two cost-effective ones relying on molecular data: i) THE REAL McCOIL method is based on putatively neutral single nucleotide polymorphism loci, and ii) the varcoding method is a fingerprinting approach that relies on the diversity and limited repertoire overlap of the var multigene family encoding the major Plasmodium falciparum blood-stage antigen PfEMP1 and is therefore under selection. In this study, we assess the robustness of the MOI estimates generated with these two approaches by simulating P. falciparum malaria dynamics under three transmission conditions using an extension of a previously developed stochastic agent-based model. We demonstrate that these approaches are complementary and best considered across distinct transmission intensities. While varcoding can underestimate MOI, it allows robust estimation, especially under high transmission where repertoire overlap is extremely limited from frequency-dependent selection. In contrast, THE REAL McCOIL often considerably overestimates MOI, but still provides reasonable estimates for low and moderate transmission. Regardless of transmission intensity, results for THE REAL McCOIL indicate that an inaccurate tail at high MOI values is generated, and that at high transmission, an apparently reasonable estimated MOI distribution can arise from some degree of compensation between overestimation and underestimation. As many countries pursue malaria elimination targets, defining the most suitable approach to estimate MOI based on sample size and local transmission intensity is highly recommended for monitoring the impact of intervention programs.

Europe as a bridgehead in the worldwide invasion history of grapevine downy mildew, Plasmopara viticola.

Europe is the historical cradle of viticulture, but grapevines (Vitis vinifera) have been increasingly threatened by pathogens of American origin. The invasive oomycete Plasmopara viticola causes downy mildew, one of the most devastating grapevine diseases worldwide. Despite major economic consequences, its invasion history remains poorly understood. We analyzed a comprehensive dataset of ∼2,000 samples, collected from the most important wine-producing countries, using nuclear and mitochondrial gene sequences and microsatellite markers. Population genetic analyses revealed very low genetic diversity in invasive downy mildew populations worldwide and little evidence of admixture. All the invasive populations originated from only one of the five native North American lineages, the one parasitizing wild summer grape (V. aestivalis). An approximate Bayesian computation-random forest approach allowed inferring the worldwide invasion scenario of P. viticola. After an initial introduction into Europe, invasive European populations served as a secondary source of introduction into vineyards worldwide, including China, South Africa, and twice independently, Australia. Only the invasion of Argentina probably represents a tertiary introduction, from Australia. Our findings provide a striking example of a global pathogen invasion resulting from secondary dispersal of a successful invasive population. Our study will also help designing quarantine regulations and efficient breeding for resistance against grapevine downy mildew.

Molecular signatures of sexual communication in the phlebotomine sand flies.

Phlebotomine sand flies employ an elaborate system of pheromone communication wherein males produce pheromones that attract other males to leks (thus acting as an aggregation pheromone) and females to the lekking males (sex pheromone). In addition, the type of pheromone produced varies among populations. Despite the numerous studies on sand fly chemical communication, little is known of their chemosensory genome. Chemoreceptors interact with chemicals in an organism's environment to elicit essential behaviors such as the identification of suitable mates and food sources. Thus, they play important roles during adaptation and speciation. Major chemoreceptor gene families, odorant receptors (ORs), gustatory receptors (GRs) and ionotropic receptors (IRs) together detect and discriminate the chemical landscape. Here, we annotated the chemoreceptor repertoire in the genomes of Lutzomyia longipalpis and Phlebotomus papatasi, major phlebotomine vectors in the New World and Old World, respectively. Comparison with other sequenced Diptera revealed a large and unique expansion where over 80% of the ~140 ORs belong to a single, taxonomically restricted clade. We next conducted a comprehensive analysis of the chemoreceptors in 63 L. longipalpis individuals from four different locations in Brazil representing allopatric and sympatric populations and three sex-aggregation pheromone types (chemotypes). Population structure based on single nucleotide polymorphisms (SNPs) and gene copy number in the chemoreceptors corresponded with their putative chemotypes, and corroborate previous studies that identified multiple populations. Our work provides genomic insights into the underlying behavioral evolution of sexual communication in the L. longipalpis species complex in Brazil, and highlights the importance of accounting for the ongoing speciation in central and South American Lutzomyia that could have important implications for vectorial capacity.

Radiation with reticulation marks the origin of a major malaria vector.

Advances in genomics have led to an appreciation that introgression is common, but its evolutionary consequences are poorly understood. In recent species radiations the sharing of genetic variation across porous species boundaries can facilitate adaptation to new environments and generate novel phenotypes, which may contribute to further diversification. Most Anopheles mosquito species that are of major importance as human malaria vectors have evolved within recent and rapid radiations of largely nonvector species. Here, we focus on one of the most medically important yet understudied anopheline radiations, the Afrotropical Anopheles funestus complex (AFC), to investigate the role of introgression in its diversification and the possible link between introgression and vector potential. The AFC comprises at least seven morphologically similar species, yet only An. funestus sensu stricto is a highly efficient malaria vector with a pan-African distribution. Based on de novo genome assemblies and additional whole-genome resequencing, we use phylogenomic and population genomic analyses to establish species relationships. We show that extensive interspecific gene flow involving multiple species pairs has shaped the evolutionary history of the AFC since its diversification. The most recent introgression event involved a massive and asymmetrical movement of genes from a distantly related AFC lineage into An. funestus, an event that predated and plausibly facilitated its subsequent dramatic geographic range expansion across most of tropical Africa. We propose that introgression may be a common mechanism facilitating adaptation to new environments and enhancing vectorial capacity in Anopheles mosquitoes.

Human Migration and the Spread of the Nematode Parasite Wuchereria bancrofti.

The human disease lymphatic filariasis causes the debilitating effects of elephantiasis and hydrocele. Lymphatic filariasis currently affects the lives of 90 million people in 52 countries. There are three nematodes that cause lymphatic filariasis, Brugia malayi, Brugia timori, and Wuchereria bancrofti, but 90% of all cases of lymphatic filariasis are caused solely by W. bancrofti (Wb). Here we use population genomics to reconstruct the probable route and timing of migration of Wb strains that currently infect Africa, Haiti, and Papua New Guinea (PNG). We used selective whole genome amplification to sequence 42 whole genomes of single Wb worms from populations in Haiti, Mali, Kenya, and PNG. Our results are consistent with a hypothesis of an Island Southeast Asia or East Asian origin of Wb. Our demographic models support divergence times that correlate with the migration of human populations. We hypothesize that PNG was infected at two separate times, first by the Melanesians and later by the migrating Austronesians. The migrating Austronesians also likely introduced Wb to Madagascar where later migrations spread it to continental Africa. From Africa, Wb spread to the New World during the transatlantic slave trade. Genome scans identified 17 genes that were highly differentiated among Wb populations. Among these are genes associated with human immune suppression, insecticide sensitivity, and proposed drug targets. Identifying the distribution of genetic diversity in Wb populations and selection forces acting on the genome will build a foundation to test future hypotheses and help predict response to current eradication efforts.

Genetic footprint of population fragmentation and contemporary collapse in a freshwater cetacean.

Understanding demographic trends and patterns of gene flow in an endangered species is crucial for devising conservation strategies. Here, we examined the extent of population structure and recent evolution of the critically endangered Yangtze finless porpoise (Neophocaena asiaeorientalis asiaeorientalis). By analysing genetic variation at the mitochondrial and nuclear microsatellite loci for 148 individuals, we identified three populations along the Yangtze River, each one connected to a group of admixed ancestry. Each population displayed extremely low genetic diversity, consistent with extremely small effective size (≤92 individuals). Habitat degradation and distribution gaps correlated with highly asymmetric gene-flow that was inefficient in maintaining connectivity between populations. Genetic inferences of historical demography revealed that the populations in the Yangtze descended from a small number of founders colonizing the river from the sea during the last Ice Age. The colonization was followed by a rapid population split during the last millennium predating the Chinese Modern Economy Development. However, genetic diversity showed a clear footprint of population contraction over the last 50 years leaving only ~2% of the pre-collapsed size, consistent with the population collapses reported from field studies. This genetic perspective provides background information for devising mitigation strategies to prevent this species from extinction.

Genetic analysis reveals efficient sexual spore dispersal at a fine spatial scale in Armillaria ostoyae, the causal agent of root-rot disease in conifers.

Armillaria ostoyae (sometimes named Armillaria solidipes) is a fungal species causing root diseases in numerous coniferous forests of the northern hemisphere. The importance of sexual spores for the establishment of new disease centres remains unclear, particularly in the large maritime pine plantations of southwestern France. An analysis of the genetic diversity of a local fungal population distributed over 500 ha in this French forest showed genetic recombination between genotypes to be frequent, consistent with regular sexual reproduction within the population. The estimated spatial genetic structure displayed a significant pattern of isolation by distance, consistent with the dispersal of sexual spores mostly at the spatial scale studied. Using these genetic data, we inferred an effective density of reproductive individuals of 0.1-0.3 individuals/ha, and a second moment of parent-progeny dispersal distance of 130-800 m, compatible with the main models of fungal spore dispersal. These results contrast with those obtained for studies of A. ostoyae over larger spatial scales, suggesting that inferences about mean spore dispersal may be best performed at fine spatial scales (i.e. a few kilometres) for most fungal species.

Genetic signature of a range expansion and leap-frog event after the recent invasion of Europe by the grapevine downy mildew pathogen Plasmopara viticola.

Biologic invasions can have important ecological, economic and social consequences, particularly when they involve the introduction and spread of plant invasive pathogens, as they can threaten natural ecosystems and jeopardize the production of human food. Examples include the grapevine downy mildew, caused by the oomycete Plasmopara viticola, an invasive species native to North America, introduced into Europe in the 1870s. We investigated the introduction and spread of this invasive pathogen, by analysing its genetic structure and diversity in a large sample from European vineyards. Populations of P. viticola across Europe displayed little genetic diversity, consistent with the occurrence of a bottleneck at the time of introduction. Bayesian coalescent analyses revealed a clear population expansion signal in the genetic data. We detected a weak, but significant, continental-wide population structure, with two geographically and genetically distinct clusters in Western and Eastern European vineyards. Approximate Bayesian computation, analyses of clines of genetic diversity and of isolation-by-distance patterns provided evidence for a wave of colonization moving in an easterly direction across Europe. This is consistent with historical reports, first mentioning the introduction of the disease in Bordeaux vineyards (France) and sub-sequently documenting its rapid spread across Europe. This initial introduction in the west was probably followed by a 'leap-frog' event into Eastern Europe, leading to the formation of the two genetic clusters we detected. This study shows that recent population genetics methods within the Bayesian and coalescence frameworks are extremely powerful for increasing our understanding of pathogen population dynamics and invasion histories.