Search

Abstract Gene drives are a promising means of malaria control with the potential to cause sustained reductions in transmission. In real environments, however, their impacts will depend on local ecological and epidemiological factors. We develop a data-driven model to investigate the impacts of gene drives that causes vector population suppression. We simulate gene drive releases in sixteen ~ 12,000 km2 areas of west Africa that span variation in vector ecology and malaria prevalence, and estimate reductions in vector abundance, malaria prevalence and clinical cases. Average reductions in vector abundance ranged from 71.6–98.4% across areas, while impacts on malaria depended strongly on which vector species were targeted. When other new interventions including RTS,S vaccination and pyrethroid-PBO bednets were in place, at least 60% more clinical cases were averted when gene drives were added, demonstrating the benefits of integrated interventions. Our results show that different strategies for gene drive implementation may be required across different African settings.

Abstract Background Intensive deployment of insecticide based malaria vector control tools resulted in the rapid evolution of phenotypes resistant to these chemicals. Understanding this process at the genomic level is important for the deployment of successful vector control interventions. Therefore, longitudinal sampling followed by whole genome sequencing (WGS) is necessary to understand how these evolutionary processes evolve over time. This study investigated the change in genetic structure and the evolution of the insecticide resistance variants in natural populations of Anopheles gambiae over time and space from 2012 to 2017 in Burkina Faso. Methods New genomic data have been generated from An. gambiae mosquitoes collected from three villages in the western part of Burkina Faso between 2012 and 2017. The samples were whole-genome sequenced and the data used in the An. gambiae 1000 genomes (Ag1000G) project as part of the Vector Observatory. Genomic data were analysed using the analysis pipeline previously designed by the Ag1000G project. Results The results showed similar and consistent nucleotide diversity and negative Tajima’s D between An. gambiae sensu stricto (s.s.) and Anopheles coluzzii. Principal component analysis (PCA) and the fixation index (F ST ) showed a clear genetic structure in the An. gambiae sensu lato (s.l.) species. Genome-wide F ST and H12 scans identified genomic regions under divergent selection that may have implications in the adaptation to ecological changes. Novel voltage-gated sodium channel pyrethroid resistance target-site alleles (V402L, I1527T) were identified at increasing frequencies alongside the established alleles (Vgsc-L995F, Vgsc-L995S and N1570Y) within the An. gambiae s.l. populations. Organophosphate metabolic resistance markers were also identified, at increasing frequencies, within the An. gambiae s.s. populations from 2012 to 2017, including the SNP Ace1-G280S and its associated duplication. Variants simultaneously identified in the same vector populations raise concerns about the long-term efficacy of new generation bed nets and the recently organophosphate pirimiphos-methyl indoor residual spraying in Burkina Faso. Conclusion These findings highlighted the benefit of genomic surveillance of malaria vectors for the detection of new insecticide resistance variants, the monitoring of the existing resistance variants, and also to get insights into the evolutionary processes driving insecticide resistance.

Abstract Sustainable reductions in African malaria transmission require innovative tools for mosquito control. One proposal involves the use of low-threshold gene drive in Anopheles vector species, where a ‘causal pathway’ would be initiated by (i) the release of a gene drive system in target mosquito vector species, leading to (ii) its transmission to subsequent generations, (iii) its increase in frequency and spread in target mosquito populations, (iv) its simultaneous propagation of a linked genetic trait aimed at reducing vectorial capacity for Plasmodium, and (v) reduced vectorial capacity for parasites in target mosquito populations as the gene drive system reaches fixation in target mosquito populations, causing (vi) decreased malaria incidence and prevalence. Here the scope, objectives, trial design elements, and approaches to monitoring for initial field releases of such gene dive systems are considered, informed by the successful implementation of field trials of biological control agents, as well as other vector control tools, including insecticides, Wolbachia, larvicides, and attractive-toxic sugar bait systems. Specific research questions to be addressed in initial gene drive field trials are identified, and adaptive trial design is explored as a potentially constructive and flexible approach to facilitate testing of the causal pathway. A fundamental question for decision-makers for the first field trials will be whether there should be a selective focus on earlier points of the pathway, such as genetic efficacy via measurement of the increase in frequency and spread of the gene drive system in target populations, or on wider interrogation of the entire pathway including entomological and epidemiological efficacy. How and when epidemiological efficacy will eventually be assessed will be an essential consideration before decisions on any field trial protocols are finalized and implemented, regardless of whether initial field trials focus exclusively on the measurement of genetic efficacy, or on broader aspects of the causal pathway. Statistical and modelling tools are currently under active development and will inform such decisions on initial trial design, locations, and endpoints. Collectively, the considerations here advance the realization of developer ambitions for the first field trials of low-threshold gene drive for malaria vector control within the next 5 years.

Abstract The source of malaria vector populations that re-establish at the beginning of the rainy season is still unclear yet knowledge of mosquito behaviour is required to effectively institute control measures. Alternative hypotheses like aestivation, local refugia, migration between neighbouring sites, and long-distance migration (LDM) are stipulated to support mosquito persistence. This work assessed the malaria vector persistence dynamics and examined various studies done on vector survival via these hypotheses; aestivation, local refugia, local or long-distance migration across sub-Saharan Africa, explored a range of methods used, ecological parameters and highlighted the knowledge trends and gaps. The results about a particular persistence mechanism that supports the re-establishment of Anopheles gambiae, Anopheles coluzzii or Anopheles arabiensis in sub-Saharan Africa were not conclusive given that each method used had its limitations. For example, the Mark-Release-Recapture (MRR) method whose challenge is a low recapture rate that affects its accuracy, and the use of time series analysis through field collections whose challenge is the uncertainty about whether not finding mosquitoes during the dry season is a weakness of the conventional sampling methods used or because of hidden shelters. This, therefore, calls for further investigations emphasizing the use of ecological experiments under controlled conditions in the laboratory or semi-field, and genetic approaches, as they are known to complement each other. This review, therefore, unveils and assesses the uncertainties that influence the different malaria vector persistence mechanisms and provides recommendations for future studies.

The sterile insect technique (SIT) has been successful in controlling some pest species but is not practicable for many others due to the large number of individuals that need to be reared and released. Previous computer modelling has demonstrated that the release of males carrying a Y-linked editor that kills or sterilises female descendants could be orders of magnitude more efficient than SIT while still remaining spatially restricted, particularly if combined with an autosomal sex distorter. In principle, further gains in efficiency could be achieved by using a self-propagating double drive design, in which each of the two components (the Y-linked editor and the sex ratio distorter) boosted the transmission of the other. To better understand the expected dynamics and impact of releasing constructs of this new design we have analysed a deterministic population genetic and population dynamic model. Our modelling demonstrates that this design can suppress a population from very low release rates, with no invasion threshold. Importantly, the design can work even if homing rates are low and sex chromosomes are silenced at meiosis, potentially expanding the range of species amenable to such control. Moreover, the predicted dynamics and impacts can be exquisitely sensitive to relatively small (e.g., 25%) changes in allele frequencies in the target population, which could be exploited for sequence-based population targeting. Analysis of published Anopheles gambiae genome sequences indicates that even for weakly differentiated populations with an FST of 0.02 there may be thousands of suitably differentiated genomic sites that could be used to restrict the spread and impact of a release. Our proposed design, which extends an already promising development pathway based on Y-linked editors, is therefore a potentially useful addition to the menu of options for genetic biocontrol.

Abstract 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 (dsxF CRISPRh ) 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 dsxF CRISPRh 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 dsxF CRISPRh 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 dsxF CRISPRh 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.

Abstract 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.

Abstract Background Measures of linkage disequilibrium (LD) play a key role in a wide range of applications from disease association to demographic history estimation. The true population LD cannot be measured directly and instead can only be inferred from genetic samples, which are unavoidably subject to measurement error. Previous studies of r 2 (a measure of LD), such as the bias due to finite sample size and its variance, were based on the special case that the true population-wise LD is zero. These results generally do not hold for non-zero rtrue2 $$ {r}_{true}^2 $$ values, which are more common in real genetic data. Results This work generalises the estimation of r 2 to all levels of LD, and for both phased and unphased data. First, we provide new formulae for the effect of finite sample size on the observed r 2 values. Second, we find a new empirical formula for the variance of the observed r 2, equals to 2E[r 2](1 − E[r 2])/n, where n is the diploid sample size. Third, we propose a new routine, Constrained ML, a likelihood-based method to directly estimate haplotype frequencies and r 2 from diploid genotypes under Hardy-Weinberg Equilibrium. While serving the same purpose as the pre-existing Expectation-Maximisation algorithm, the new routine can have better convergence and is simpler to use. A new likelihood-ratio test is also introduced to test for the absence of a particular haplotype. Extensive simulations are run to support these findings. Conclusion Most inferences on LD will benefit from our new findings, from point and interval estimation to hypothesis testing. Genetic analyses utilising r 2 information will become more accurate as a result.

CRISPR-based homing gene drives can be designed to disrupt essential genes whilst biasing their own inheritance, leading to suppression of mosquito populations in the laboratory. This class of gene drives relies on CRISPR-Cas9 cleavage of a target sequence and copying ('homing') therein of the gene drive element from the homologous chromosome. However, target site mutations that are resistant to cleavage yet maintain the function of the essential gene are expected to be strongly selected for. Targeting functionally constrained regions where mutations are not easily tolerated should lower the probability of resistance. Evolutionary conservation at the sequence level is often a reliable indicator of functional constraint, though the actual level of underlying constraint between one conserved sequence and another can vary widely. Here we generated a novel adult lethal gene drive (ALGD) in the malaria vector Anopheles gambiae, targeting an ultra-conserved target site in a haplosufficient essential gene (AGAP029113) required during mosquito development, which fulfils many of the criteria for the target of a population suppression gene drive. We then designed a selection regime to experimentally assess the likelihood of generation and subsequent selection of gene drive resistant mutations at its target site. We simulated, in a caged population, a scenario where the gene drive was approaching fixation, where selection for resistance is expected to be strongest. Continuous sampling of the target locus revealed that a single, restorative, in-frame nucleotide substitution was selected. Our findings show that ultra-conservation alone need not be predictive of a site that is refractory to target site resistance. Our strategy to evaluate resistance in vivo could help to validate candidate gene drive targets for their resilience to resistance and help to improve predictions of the invasion dynamics of gene drives in field populations.

Abstract 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.

Synthetic gene drive constructs could, in principle, provide the basis for highly efficient interventions to control disease vectors and other pest species. This efficiency derives in part from leveraging natural processes of dispersal and gene flow to spread the construct and its impacts from one population to another. However, sometimes (for example, with invasive species) only specific populations are in need of control, and impacts on non-target populations would be undesirable. Many gene drive designs use nucleases that recognise and cleave specific genomic sequences, and one way to restrict their spread would be to exploit sequence differences between target and non-target populations. In this paper we propose and model a series of low threshold double drive designs for population suppression, each consisting of two constructs, one imposing a reproductive load on the population and the other inserted into a differentiated locus and controlling the drive of the first. Simple deterministic, discrete-generation computer simulations are used to assess the alternative designs. We find that the simplest double drive designs are significantly more robust to pre-existing cleavage resistance at the differentiated locus than single drive designs, and that more complex designs incorporating sex ratio distortion can be more efficient still, even allowing for successful control when the differentiated locus is neutral and there is up to 50% pre-existing resistance in the target population. Similar designs can also be used for population replacement, with similar benefits. A population genomic analysis of CRISPR PAM sites in island and mainland populations of the malaria mosquito Anopheles gambiae indicates that the differentiation needed for our methods to work can exist in nature. Double drives should be considered when efficient but localised population genetic control is needed and there is some genetic differentiation between target and non-target populations.

AbstractDNA sequences that are exactly conserved over long evolutionary time scales have been observed in a variety of taxa. Such sequences are likely under strong functional constraint and they have been useful in the field of comparative genomics for identifying genome regions with regulatory function. A potential new application for these ultra-conserved elements (UCEs) has emerged in the development of gene drives to control mosquito populations. Many gene drives work by recognizing and inserting at a specific target sequence in the genome, often imposing a reproductive load as a consequence. They can therefore select for target sequence variants that provide resistance to the drive. Focusing on highly conserved, highly constrained sequences lowers the probability that variant, gene drive-resistant alleles can be tolerated. Here, we search for conserved sequences of 18 bp and over in an alignment of 21 AnophelesAnophelesAnophelesAnophelesAnopheles

Homing-based gene drives use a germline source of nuclease to copy themselves at specific target sites in a genome and bias their inheritance. Such gene drives can be designed to spread and deliberately suppress populations of malaria mosquitoes by impairing female fertility. However, strong unintended fitness costs of the drive and a propensity to generate resistant mutations can limit a gene drive's potential to spread. Alternative germline regulatory sequences in the drive element confer improved fecundity of carrier individuals and reduced propensity for target site resistance. This is explained by reduced rates of end-joining repair of DNA breaks from parentally deposited nuclease in the embryo, which can produce heritable mutations that reduce gene drive penetrance. We tracked the generation and selection of resistant mutations over the course of a gene drive invasion of a population. Improved gene drives show faster invasion dynamics, increased suppressive effect and later onset of target site resistance. Our results show that regulation of nuclease expression is as important as the choice of target site when developing a robust homing-based gene drive for population suppression.

Multiple national and international trends and drivers are radically changing what biological security means for the United Kingdom (UK). New technologies present novel opportunities and challenges, and globalisation has created new pathways and increased the speed, volume and routes by which organisms can spread. The UK Biological Security Strategy (2018) acknowledges the importance of research on biological security in the UK. Given the breadth of potential research, a targeted agenda identifying the questions most critical to effective and coordinated progress in different disciplines of biological security is required. We used expert elicitation to generate 80 policy-relevant research questions considered by participants to have the greatest impact on UK biological security. Drawing on a collaboratively-developed set of 450 questions, proposed by 41 experts from academia, industry and the UK government (consulting 168 additional experts) we subdivided the final 80 questions into six categories: bioengineering; communication and behaviour; disease threats (including pandemics); governance and policy; invasive alien species; and securing biological materials and securing against misuse. Initially, the questions were ranked through a voting process and then reduced and refined to 80 during a one-day workshop with 35 participants from a variety of disciplines. Consistently emerging themes included: the nature of current and potential biological security threats, the efficacy of existing management actions, and the most appropriate future options. The resulting questions offer a research agenda for biological security in the UK that can assist the targeting of research resources and inform the implementation of the UK Biological Security Strategy. These questions include research that could aid with the mitigation of Covid-19, and preparation for the next pandemic. We hope that our structured and rigorous approach to creating a biological security research agenda will be replicated in other countries and regions. The world, not just the UK, is in need of a thoughtful approach to directing biological security research to tackle the emerging issues.

Abstract Background The use of gene drive systems to manipulate populations of malaria vectors is currently being investigated as a method of malaria control. One potential system uses driving endonuclease genes (DEGs) to spread genes that impose a genetic load. Previously, models have shown that the introduction of DEG-bearing mosquitoes could suppress or even extinguish vector populations in spatially-heterogeneous environments which were constant over time. In this study, a stochastic spatially-explicit model of mosquito ecology is combined with a rainfall model which enables the generation of a variety of daily precipitation patterns. The model is then used to investigate how releases of a DEG that cause a bias in population sex ratios towards males are affected by seasonal or random rainfall patterns. The parameters of the rainfall model are then fitted using data from Bamako, Mali, and Mbita, Kenya, to evaluate release strategies in similar climatic conditions. Results In landscapes with abundant resources and large mosquito populations the spread of a DEG is reliable, irrespective of variability in rainfall. This study thus focuses mainly on landscapes with low density mosquito populations where the spread of a DEG may be sensitive to variation in rainfall. It is found that an introduced DEG will spread into its target population more reliably in wet conditions, yet an established DEG will have more impact in dry conditions. In strongly seasonal environments, it is thus preferable to release DEGs at the onset of a wet season to maximize their spread before the following dry season. If the variability in rainfall has a substantial random component, there is a net increase in the probability that a DEG release will lead to population extinction, due to the increased impact of a DEG which manages to establish in these conditions. For Bamako, where annual rainfall patterns are characterized by a long dry season, it is optimal to release a DEG at the start of the wet season, where the population is growing fastest. By contrast release timing is of lower importance for the less seasonal Mbita. Conclusion This analysis suggests that DEG based methods of malaria vector control can be effective in a wide range of climates. In environments with substantial temporal variation in rainfall, careful timing of releases which accounts for the temporal variation in population density can substantially improve the probability of mosquito suppression or extinction.

Abstract Driving endonuclease genes (DEGs) spread through a population by a non-Mendelian mechanism. In a heterozygote, the protein encoded by a DEG causes a double-strand break in the homologous chromosome opposite to where its gene is inserted and when the break is repaired using the homologue as a template the DEG heterozygote is converted to a homozygote. Some DEGs occur naturally while several classes of endonucleases can be engineered to spread in this way, with CRISPR-Cas9 based systems being particularly flexible. There is great interest in using driving endonuclease genes to impose a genetic load on insects that vector diseases or are economic pests to reduce their population density, or to introduce a beneficial gene such as one that might interrupt disease transmission. This paper reviews both the population genetics and population dynamics of DEGs. It summarises the theory that guides the design of DEG constructs intended to perform different functions. It also reviews the studies that have explored the likelihood of resistance to DEG phenotypes arising, and how this risk may be reduced. The review is intended for a general audience and mathematical details are kept to a minimum.

Abstract Background The introduction of new malaria control interventions has often led to the evolution of resistance, both of the parasite to new drugs and of the mosquito vector to new insecticides, compromising the efficacy of the interventions. Recent progress in molecular and population biology raises the possibility of new genetic-based interventions, and the potential for resistance to evolve against these should be considered. Here, population modelling is used to determine the main factors affecting the likelihood that resistance will evolve against a synthetic, nuclease-based driving Y chromosome that produces a male-biased sex ratio. Methods A combination of deterministic differential equation models and stochastic analyses involving branching processes and Gillespie simulations is utilized to assess the probability that resistance evolves against a driving Y that otherwise is strong enough to eliminate the target population. The model considers resistance due to changes at the target site such that they are no longer cleaved by the nuclease, and due to trans-acting autosomal suppressor alleles. Results The probability that resistance evolves increases with the mutation rate and the intrinsic rate of increase of the population, and decreases with the strength of drive and any pleiotropic fitness costs of the resistant allele. In seasonally varying environments, the time of release can also affect the probability of resistance evolving. Trans-acting suppressor alleles are more likely to suffer stochastic loss at low frequencies than target site resistant alleles. Conclusions As with any other intervention, there is a risk that resistance will evolve to new genetic approaches to vector control, and steps should be taken to minimize this probability. Two design features that should help in this regard are to reduce the rate at which resistant mutations arise, and to target sequences such that if they do arise, they impose a significant fitness cost on the mosquito.