Schematic illustration of transmission rules and acquisition of host immunity within the compartmental ODE model

(see Figure 1—figure Supplement 1 for a detailed representation of the compartment model). (A) Rules for new infections given the host’s past infection history and current multiplicity of infection (i.e., MOI). Upon transmission of a specific parasite strain A, if the host has had an infection of strain A in the past (hands raised), a new infection will not be added to the current MOI; instead, the infection will be considered cleared and added to the total number of cleared infections; if the host is new to strain A and does not have specific immunity to it (inferred from Eq. 1), a new infection will be added (i.e., MOI increase by 1) as long as MOI does not exceed the carrying capacity of coexisting strains. (B) Rules of symptomatic infections and treatment in the different generalized immunity (G) classes. With increasing generalized immunity (G), hosts are less likely to show clinical symptoms. Hosts in G0 have a risk of death in addition to symptomatic infections; Hosts in G1 do not die from infections but show symptoms upon new infections; Hosts in G2 carry asymptomatic infections most of the time with a slight chance of showing symptoms. Symptomatic infections result in a daily treatment rate that removes the infections caused by wild-type strains. Hosts that have cleared enough number of infections will move to the next G class. Hosts will move back to a lower G class when the generalized immunity memory is slowly lost if not boosted by constant infections.

Figure 1—figure supplement 1. Compartment model of drug resistance evolution.

The frequency of resistance under varying strain diversity and vectorial capacity.

(A) The heatmap shows a nonlinear parasite prevalence response given increasing vectorial capacity and the number of strains under no drug treatment, with warmer colors representing high prevalence and cooler colors representing low prevalence. X and Y axes correspond to increasing vectorial capacity and the number of strains in logarithmic scales. White tiles indicate the highest prevalence given a fixed number of strains. (B) The heatmaps show resistance frequencies under varying strain diversity and vectorial capacity at two levels of drug treatment rate, with warmer colors representing higher resistance frequency (in this example, ssingle = 0.1, smixed = 0.9). (C) A negative relationship between parasite prevalence and resistance frequency. The color of the points indicates combinations of resistance fitness costs in hosts with resistant strains alone (ssingle) or mixed infections of resistant and wild-type strains (smixed).

Figure 2—figure supplement 1. Prevalence given the combination of vectorial capacity and the number of strains from no treatment to high treatment rate for wild-type-only infections.

Figure 2—figure supplement 2. Infectivity of a new infection as a function of the number of strains and mean immunity.

Figure 2—figure supplement 3. Relationship between parasite prevalence and resistance frequency under full treatment.

Figure 2—figure supplement 4. Relationship between parasite prevalence and resistance frequency under partial treatment.

Global distribution of chloroquine-resistant genotype (pfcrt 76T) against P. falciparum prevalence in children between 2-10 yrs old.

Sampling between 1990 and 2000 was included to ensure genotyping was performed largely before the policy switch of the first-line antimalarial drugs to ACT. Different shapes indicate samples from different continents, while shape sizes correspond to sample sizes for genotyping (see Methods for details).

Relationship between host immunity, drug treatment, and resistance evolution.

Fraction of hosts in different G classes with increasing strain diversity and the corresponding vectorial capacity indicated by white circles in Figure 1A at equilibrium or year 50 after the invasion of resistant genotypes. Hosts in drug-treated classes are indicated by stripes. Red dotted lines show the corresponding frequency of resistance. The upper panel is generated under wild-type-only infections with increasing treatment rates. The lower panel represents resistance-only infections without treatment or resistant invasion under treatments.

Temporal trajectories of resistance invasion.

Host and parasite dynamics under resistance invasion are shown for lower (nstrains = 20) A and higher (nstrains = 113) diversity B. Because drug treatment does not affect resistant parasites, they surge quickly after introduction, thus leading to more infections. Hosts recovered from a large number of new infections move into higher G classes (from year 1-8). The high host immune protection selects against resistant parasites (year 4-10). Under low diversity, wild-type parasites quickly go to extinction A. Under high diversity, the less symptomatic G2 class provides a niche for wild-type parasites to multiply (year 4-10), where the two genotypes co-exist. Meanwhile, resistant parasites dominate in hosts that are in G0 and G1 B.

Changes in frequency of resistance after the first-line drug is changed.

Each trajectory represents the combination of variables indicated by the white circles in Figure 1A. Color from cool to warm represents increasing diversity in strains. Here the usage of the drug, to which parasites have developed resistance, is reduced to 0.52, 0.52, 0.52, 0.52, 0.21, 0.21, 0.21, 0.21, 0, 0, 0, 0, 0, 0, 0, 0 each year following the change in the treatment regime. The trajectory of reduction in resistant drug usage follows the usage survey in Western Kenya from 2003 to 2018 (Hemming-Schroeder et al., 2018).

Figure 6—figure supplement 1. Percentage of reduction in resistance after one year.

Changes in frequency of resistant genotypes across different biogeographic regions.

Each circle represents one studied sample (at least 20 infected hosts) from one geographic location. Circles connected by dotted lines represent longitudinal samples from the same study. After the policy switch in first-line antimalarial drugs, frequencies of resistance decreased gradually in Africa, but maintained high in Asia, Oceania, and South America despite the policy change for more than 20 years. CQ: chloroquine; SP: sulfadoxine-pyrimethamine; MQ: mefloquine; AQ: amodiaquine; PQ: primaquine; QN-TET: quinine+tetracycline; ACT: artemisinin-based combination therapy.

Relationship between parasite prevalence and resistance frequency for the generalized-immunity-only model.

Paths are connected from low vectorial capacity to high vectorial capacity. Colors represent different combinations of single-genotype infection cost and mixed-genotype infection cost of resistant parasites.

Figure 8—figure supplement 1. Relationship between host immunity, drug treatment, and resistance evolution for the generalized-immunity-only model.

Source of drug policy data.

Epidemiological parameters used for numerically solving ODEs. All rates are measured per day; time is measured in days.

Compartment model of drug resistance evolution.

(A) The number of hosts and movements are tracked in different generalized immunity classes (G), to-gether with their drug treatment states (treated, D; untreated, U); (B) wild-type (P W) and resistant parasite (P R) population sizes are tracked in different host immunity classes; (C) Changes in total immunity (T I, total number of cleared infections) per G class are followed. See Appendix 1 for a detailed explanation of the ODE system.

Prevalence given the combination of vectorial capacity and the number of strains from no treatment to high treatment rate for wild-type-only infections.

Grey areas indicate that transmission is eliminated.

Infectivity of a new infection as a function of the number of strains and mean immunity.

(Total immunity divided by the number of hosts per G class) (see Eq. 1. ssingle : 0.1; smixed : 0.9.

Relationship between parasite prevalence and resistance frequency under full treatment (daily treatment rate d1 = 0.2).

Each subgraph represents the combination of resistance fitness costs in hosts with resistant strains alone (ssingle) and mixed-genotype infections of resistant and wild-type strains (smixed), as well as the efficacy of resistance . Color indicates vectorial capacity.

Relationship between parasite prevalence and resistance frequency under partial treatment (daily treatment rate d1 = 0.2).

Each subgraph represents the combination of resistance fitness costs in hosts with resistant strains alone (ssingle) or mixed-genotype infections of resistant and wild-type strains (smixed). Color indicates vectorial capacity, as well as the efficacy of resistance .

Percentage of reduction in resistance after one year of policy change in drug treatment as a function of vectorial capacity and the number of strains under different combinations of resistance costs (ssingle ;smixed).

Relationship between host immunity, drug treatment, and resistance evolution for the generalized-immunity-only model.

Note that in the generalized-immunity-only model, there is no strain diversity. The only parameter that determines transmission intensity is vectorial capacity. In general, prevalence (blue dotted line) increases as vectorial capacity increases despite hosts increasingly concentrating in G2 class (A). The fraction of resistant parasites decreases initially with increasing vectorial capacity, but rises again as high transmission results in a higher proportion of G2 hosts in the drug-treated class (B).