Systematic identification of plausible pathways to potential harm via problem formulation for investigational releases of a population suppression gene drive to control the human malaria vector Anopheles gambiae in West Africa.

BACKGROUND: Population suppression gene drive has been proposed as a strategy for malaria vector control. A CRISPR-Cas9-based transgene homing at the doublesex locus (dsxFCRISPRh) has recently been shown to increase rapidly in frequency in, and suppress, caged laboratory populations of the malaria mosquito vector Anopheles gambiae. Here, problem formulation, an initial step in environmental risk assessment (ERA), was performed for simulated field releases of the dsxFCRISPRh transgene in West Africa. METHODS: Building on consultative workshops in Africa that previously identified relevant environmental and health protection goals for ERA of gene drive in malaria vector control, 8 potentially harmful effects from these simulated releases were identified. These were stratified into 46 plausible pathways describing the causal chain of events that would be required for potential harms to occur. Risk hypotheses to interrogate critical steps in each pathway, and an analysis plan involving experiments, modelling and literature review to test each of those risk hypotheses, were developed. RESULTS: Most potential harms involved increased human (n = 13) or animal (n = 13) disease transmission, emphasizing the importance to subsequent stages of ERA of data on vectorial capacity comparing transgenics to non-transgenics. Although some of the pathways (n = 14) were based on known anatomical alterations in dsxFCRISPRh homozygotes, many could also be applicable to field releases of a range of other transgenic strains of mosquito (n = 18). In addition to population suppression of target organisms being an accepted outcome for existing vector control programmes, these investigations also revealed that the efficacy of population suppression caused by the dsxFCRISPRh transgene should itself directly affect most pathways (n = 35). CONCLUSIONS: Modelling will play an essential role in subsequent stages of ERA by clarifying the dynamics of this relationship between population suppression and reduction in exposure to specific potential harms. This analysis represents a comprehensive identification of plausible pathways to potential harm using problem formulation for a specific gene drive transgene and organism, and a transparent communication tool that could inform future regulatory studies, guide subsequent stages of ERA, and stimulate further, broader engagement on the use of population suppression gene drive to control malaria vectors in West Africa.

Modelling the suppression of a malaria vector using a CRISPR-Cas9 gene drive to reduce female fertility.

BACKGROUND: Gene drives based on CRISPR-Cas9 technology are increasingly being considered as tools for reducing the capacity of mosquito populations to transmit malaria, and one of the most promising options is driving endonuclease genes that reduce the fertility of female mosquitoes. In particular, there is much interest in constructs that target the conserved mosquito doublesex (dsx) gene such that the emergence of functional drive-resistant alleles is unlikely. Proof of principle that these constructs can lead to substantial population suppression has been obtained in population cages, and they are being evaluated for use in sub-Saharan Africa. Here, we use simulation modelling to understand the factors affecting the spread of this type of gene drive over a one million-square kilometre area of West Africa containing substantial environmental and social heterogeneity. RESULTS: We found that a driving endonuclease gene targeting female fertility could lead to substantial reductions in malaria vector populations on a regional scale. The exact level of suppression is influenced by additional fitness costs of the transgene such as the somatic expression of Cas9, and its deposition in sperm or eggs leading to damage to the zygote. In the absence of these costs, or of emergent drive-resistant alleles that restore female fertility, population suppression across the study area is predicted to stabilise at ~ 95% 4 years after releases commence. Small additional fitness costs do not greatly affect levels of suppression, though if the fertility of females whose offspring transmit the construct drops by more than ~ 40%, then population suppression is much less efficient. We show the suppression potential of a drive allele with high fitness costs can be enhanced by engineering it also to express male bias in the progeny of transgenic males. Irrespective of the strength of the drive allele, the spatial model predicts somewhat less suppression than equivalent non-spatial models, in particular in highly seasonal regions where dry season stochasticity reduces drive efficiency. We explored the robustness of these results to uncertainties in mosquito ecology, in particular their method of surviving the dry season and their dispersal rates. CONCLUSIONS: The modelling presented here indicates that considerable suppression of vector populations can be achieved within a few years of using a female sterility gene drive, though the impact is likely to be heterogeneous in space and time.

Modelling the potential of genetic control of malaria mosquitoes at national scale.

BACKGROUND: The persistence of malaria in large parts of sub-Saharan Africa has motivated the development of novel tools to complement existing control programmes, including gene-drive technologies to modify mosquito vector populations. Here, we use a stochastic simulation model to explore the potential of using a driving-Y chromosome to suppress vector populations in a 106 km2 area of West Africa including all of Burkina Faso. RESULTS: The consequence of driving-Y introductions is predicted to vary across the landscape, causing elimination of the target species in some regions and suppression in others. We explore how this variation is determined by environmental conditions, mosquito behaviour, and the properties of the gene-drive. Seasonality is particularly important, and we find population elimination is more likely in regions with mild dry seasons whereas suppression is more likely in regions with strong seasonality. CONCLUSIONS: Despite the spatial heterogeneity, we suggest that repeated introductions of modified mosquitoes over a few years into a small fraction of human settlements may be sufficient to substantially reduce the overall number of mosquitoes across the entire geographic area.

Modelling the persistence of mosquito vectors of malaria in Burkina Faso.

BACKGROUND: Populations of the Anopheles gambiae complex are found during the rainy season throughout West Africa, even in arid areas with long dry seasons during which mosquitoes appear to be absent. Several hypotheses have been proposed to explain this apparent paradox, including aestivation, dispersal between neighbouring settlements, and long distance migration using high-altitude wind currents. METHODS: An individual-based, spatially explicit model of mosquito populations was developed for a region of West Africa centred on, and including all of, Burkina Faso. Populations associated with human settlements were linked by dispersal and the model incorporated geospatial data on the distribution of settlements, water bodies and rainfall. RESULTS: Local dispersal (at rates consistent with experimental data) was necessary to explain observed patterns of rainy season populations across all of the simulation area, but by itself failed to account for the presence of populations in the arid North (the Sahel). The presence of rare dry-season larval sites could explain these northern populations, but seems inconsistent with field surveys. Aestivation by female mosquitoes explained rainy-season populations in all but the very sparsest and driest areas of human habitation, while long-distance migration based on annual wind patterns could account for all observed populations. CONCLUSIONS: Modelling studies such as this can help assess the potential validity of different hypotheses and suggest priority areas for experimental study. In particular, the results highlight a shortage of empirical research on mosquito dispersal between neighbouring settlements, which may be critically important to the continued presence of many mosquito populations in West Africa. Further research that establishes the extent to which mosquitoes aestivate, and migrate using high altitude winds, is also much needed to understand Sahelian mosquito populations.

The dynamics of disease in a metapopulation: The role of dispersal range.

The establishment and spread of a disease within a metapopulation is influenced both by dynamics within each population and by the host and pathogen spatial processes through which they are connected. We develop a spatially explicit metapopulation model to investigate how the form of host and disease dispersal jointly influence the probability of disease establishment and invasion. We show that diseases are more likely to establish if both the host and the disease tend to disperse locally, since the former leads to the spatial aggregation of host populations in the environment while the latter facilitates the pathogen's exploitation of this spatial pattern. In contrast, local pathogen dispersal is likely to reduce the probability of subsequent disease spread because it increases the spatial segregation of infected and uninfected populations. The effects of local dispersal on disease dynamics are less pronounced when the pathogen spreads through the movement of infected hosts and more pronounced when pathogen dispersal is independent (for example through airborne viruses) though the details of host and pathogen biology can be important. These spatial effects tend to be more pronounced if the sites available for host occupation are themselves spatially aggregated.

Anopheles gambiae (s.l.) is found where few are looking: assessing mosquito diversity and density outside inhabited areas using diverse sampling methods.

BACKGROUND: One of the promising current approaches to curb malaria lies in genetic vector control, the implementation of which will require an improved understanding of the movement of genetic constructs among mosquito populations. To predict potential gene flow from one area to another, it is important to begin to understand mosquito dynamics outside of the commonly-sampled village areas, and thus how genes may move between villages. This study assessed the presence and relative abundance of mosquitoes in a 6-km corridor between two villages in western Burkina Faso. METHODS: The area surrounding the villages was mapped and the road between them was used as the basis of a transect along which to sample. Five collection points were placed along this transect. To investigate both larval and adult mosquito presence, multiple sampling approaches were used surrounding each point: searching for larval sites in an area of 500 m radius, swarm sampling, human landing catches (HLC), CDC light traps and backpack aspiration catches of potential resting sites. Sampling took place twice: in September and October 2015. RESULTS: Adult mosquitoes from six species of Anopheles and three other genera were found along the whole transect. Anopheles gambiae (s.l.) was the most abundant followed by Anopheles nili and Anopheles coustani. Larvae of Anopheles spp. were found in small pools of surface water along the whole transect, though their presence increased with human proximity. HLC and aspiration were the most efficient methods of collecting adult mosquitoes along the whole transect, indicating that there are both host-seeking and resting mosquitoes well away from core village areas. In contrast, swarms of male mosquitoes, thought to be the principle mating locations of Anopheles spp. mosquitoes in West Africa, were only found close to the core village areas. CONCLUSIONS: This preliminary study indicates that Anopheles spp. mosquitoes are both present and breeding in low human-density areas along transit axes and provides both a relative evaluation of methods for use in these areas and evidence that gene flow between Sahelian population centres is likely. More robust and structured studies are nevertheless needed to come with stronger conclusions.

Detecting the population dynamics of an autosomal sex ratio distorter transgene in malaria vector mosquitoes.

The development of genetically modified (GM) mosquitoes and their subsequent field release offers innovative and cost-effective approaches to reduce mosquito-borne diseases, such as malaria. A sex-distorting autosomal transgene has been developed recently in G3 mosquitoes, a laboratory strain of the malaria vector Anopheles gambiae s.l. The transgene expresses an endonuclease called I-PpoI during spermatogenesis, which selectively cleaves the X chromosome to result in ~95% male progeny. Following the World Health Organization guidance framework for the testing of GM mosquitoes, we assessed the dynamics of this transgene in large cages using a joint experimental modelling approach.We performed a 4-month experiment in large, indoor cages to study the population genetics of the transgene. The cages were set up to mimic a simple tropical environment with a diurnal light-cycle, constant temperature and constant humidity. We allowed the generations to overlap to engender a stable age structure in the populations. We constructed a model to mimic the experiments, and used the experimental data to infer the key model parameters.We identified two fitness costs associated with the transgene. First, transgenic adult males have reduced fertility and, second, their female progeny have reduced pupal survival rates. Our results demonstrate that the transgene is likely to disappear in <3 years under our confined conditions. Model predictions suggest this will be true over a wide range of background population sizes and transgene introduction rates. Synthesis and applications. Our study is in line with the World Health Organization guidance recommendations in regard to the development and testing of GM mosquitoes. Since the transgenic sex ratio distorter strain (Ag(PMB)1) has been considered for genetic vector control of malaria, we recorded the dynamics of this transgene in indoor-large cage populations and modelled its post-release persistence under different scenarios. We provide a demonstration of the self-limiting nature of the transgene, and identified new fitness costs that will further reduce the longevity of the transgene after its release. Finally, our study has showcased an alternative and effective statistical method for characterizing the phenotypic expression of a transgene in an insect pest population.

Large-cage assessment of a transgenic sex-ratio distortion strain on populations of an African malaria vector.

BACKGROUND: Novel transgenic mosquito control methods require progressively more realistic evaluation. The goal of this study was to determine the effect of a transgene that causes a male-bias sex ratio on Anopheles gambiae target populations in large insectary cages. METHODS: Life history characteristics of Anopheles gambiae wild type and Ag(PMB)1 (aka gfp124L-2) transgenic mosquitoes, whose progeny are 95% male, were measured in order to parameterize predictive population models. Ag(PMB)1 males were then introduced at two ratios into large insectary cages containing target wild type populations with stable age distributions and densities. The predicted proportion of females and those observed in the large cages were compared. A related model was then used to predict effects of male releases on wild mosquitoes in a west African village. RESULTS: The frequency of transgenic mosquitoes in target populations reached an average of 0.44 ± 0.02 and 0.56 ± 0.02 after 6 weeks in the 1:1 and in the 3:1 release ratio treatments (transgenic male:wild male) respectively. Transgenic males caused sex-ratio distortion of 73% and 80% males in the 1:1 and 3:1 treatments, respectively. The number of eggs laid in the transgenic treatments declined as the experiment progressed, with a steeper decline in the 3:1 than in the 1:1 releases. The results of the experiment are partially consistent with predictions of the model; effect size and variability did not conform to the model in two out of three trials, effect size was over-estimated by the model and variability was greater than anticipated, possibly because of sampling effects in restocking. The model estimating the effects of hypothetical releases on the mosquito population of a West African village demonstrated that releases could significantly reduce the number of females in the wild population. The interval of releases is not expected to have a strong effect. CONCLUSIONS: The biological data produced to parameterize the model, the model itself, and the results of the experiments are components of a system to evaluate and predict the performance of transgenic mosquitoes. Together these suggest that the Ag(PMB)1 strain has the potential to be useful for reversible population suppression while this novel field develops.