A theory of resistance to multiplexed gene drive demonstrates the significant role of weakly deleterious natural genetic variation.

Evolution of resistance is a major barrier to successful deployment of gene-drive systems to suppress natural populations, which could greatly reduce the burden of many vector-borne diseases. Multiplexed guide RNAs (gRNAs) that require resistance mutations in all target cut sites are a promising antiresistance strategy since, in principle, resistance would only arise in unrealistically large populations. Using stochastic simulations that accurately model evolution at very large population sizes, we explore the probability of resistance due to three important mechanisms: 1) nonhomologous end-joining mutations, 2) single-nucleotide mutants arising de novo, or 3) single-nucleotide polymorphisms preexisting as standing variation. Our results explore the relative importance of these mechanisms and highlight a complexity of the mutation-selection-drift balance between haplotypes with complete resistance and those with an incomplete number of resistant alleles. We find that this leads to a phenomenon where weakly deleterious naturally occurring variants greatly amplify the probability of multisite resistance compared to de novo mutation. This key result provides design criterion for antiresistance multiplexed systems, which, in general, will need a larger number of gRNAs compared to de novo expectations. This theory may have wider application to the evolution of resistance or evolutionary rescue when multiple changes are required before selection can act.

Gene drive designs for efficient and localisable population suppression using Y-linked editors.

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.

Considerations for first field trials of low-threshold gene drive for malaria vector control.

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.

Double drives and private alleles for localised population genetic control.

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.

Regulating the expression of gene drives is key to increasing their invasive potential and the mitigation of resistance.

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.

Resistance to a CRISPR-based gene drive at an evolutionarily conserved site is revealed by mimicking genotype fixation.

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.

Potential persistence mechanisms of the major Anopheles gambiae species complex malaria vectors in sub-Saharan Africa: a narrative review.

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.

Gene drives and population persistence vs elimination: The impact of spatial structure and inbreeding at low density.

Synthetic gene drive constructs are being developed to control disease vectors, invasive species, and other pest species. In a well-mixed random mating population a sufficiently strong gene drive is expected to eliminate a target population, but it is not clear whether the same is true when spatial processes play a role. In species with an appropriate biology it is possible that drive-induced reductions in density might lead to increased inbreeding, reducing the efficacy of drive, eventually leading to suppression rather than elimination, regardless of how strong the drive is. To investigate this question we analyse a series of explicitly solvable stochastic models considering a range of scenarios for the relative timing of mating, reproduction, and dispersal and analyse the impact of two different types of gene drive, a Driving Y chromosome and a homing construct targeting an essential gene. We find in all cases a sufficiently strong Driving Y will go to fixation and the population will be eliminated, except in the one life history scenario (reproduction and mating in patches followed by dispersal) where low density leads to increased inbreeding, in which case the population persists indefinitely, tending to either a stable equilibrium or a limit cycle. These dynamics arise because Driving Y males have reduced mating success, particularly at low densities, due to having fewer sisters to mate with. Increased inbreeding at low densities can also prevent a homing construct from eliminating a population. For both types of drive, if there is strong inbreeding depression, then the population cannot be rescued by inbreeding and it is eliminated. These results highlight the potentially critical role that low-density-induced inbreeding and inbreeding depression (and, by extension, other sources of Allee effects) can have on the eventual impact of a gene drive on a target population.

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.

Robust Estimation of Recent Effective Population Size from Number of Independent Origins in Soft Sweeps.

Estimating recent effective population size is of great importance in characterizing and predicting the evolution of natural populations. Methods based on nucleotide diversity may underestimate current day effective population sizes due to historical bottlenecks, whereas methods that reconstruct demographic history typically only detect long-term variations. However, soft selective sweeps, which leave a fingerprint of mutational history by recurrent mutations on independent haplotype backgrounds, holds promise of an estimate more representative of recent population history. Here, we present a simple and robust method of estimation based only on knowledge of the number of independent recurrent origins and the current frequency of the beneficial allele in a population sample, independent of the strength of selection and age of the mutation. Using a forward-time theoretical framework, we show the mean number of origins is a function of θ=2Nµ and current allele frequency, through a simple equation, and the distribution is approximately Poisson. This estimate is robust to whether mutants preexisted before selection arose and is equally accurate for diploid populations with incomplete dominance. For fast (e.g., seasonal) demographic changes compared with time scale for fixation of the mutant allele, and for moderate peak-to-trough ratios, we show our constant population size estimate can be used to bound the maximum and minimum population size. Applied to the Vgsc gene of Anopheles gambiae, we estimate an effective population size of roughly 6×107, and including seasonal demographic oscillations, a minimum effective population size >3×107, and a maximum <6×109, suggesting a mean ∼109.

Estimating linkage disequilibrium from genotypes under Hardy-Weinberg equilibrium.

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 r2 (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 [Formula: see text] values, which are more common in real genetic data. RESULTS: This work generalises the estimation of r2 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 r2 values. Second, we find a new empirical formula for the variance of the observed r2, equals to 2E[r2](1 - E[r2])/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 r2 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 r2 information will become more accurate as a result.

Impact of mosquito gene drive on malaria elimination in a computational model with explicit spatial and temporal dynamics.

The renewed effort to eliminate malaria and permanently remove its tremendous burden highlights questions of what combination of tools would be sufficient in various settings and what new tools need to be developed. Gene drive mosquitoes constitute a promising set of tools, with multiple different possible approaches including population replacement with introduced genes limiting malaria transmission, driving-Y chromosomes to collapse a mosquito population, and gene drive disrupting a fertility gene and thereby achieving population suppression or collapse. Each of these approaches has had recent success and advances under laboratory conditions, raising the urgency for understanding how each could be deployed in the real world and the potential impacts of each. New analyses are needed as existing models of gene drive primarily focus on nonseasonal or nonspatial dynamics. We use a mechanistic, spatially explicit, stochastic, individual-based mathematical model to simulate each gene drive approach in a variety of sub-Saharan African settings. Each approach exhibits a broad region of gene construct parameter space with successful elimination of malaria transmission due to the targeted vector species. The introduction of realistic seasonality in vector population dynamics facilitates gene drive success compared with nonseasonal analyses. Spatial simulations illustrate constraints on release timing, frequency, and spatial density in the most challenging settings for construct success. Within its parameter space for success, each gene drive approach provides a tool for malaria elimination unlike anything presently available. Provided potential barriers to success are surmounted, each achieves high efficacy at reducing transmission potential and lower delivery requirements in logistically challenged settings.

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.

Gene drive for population genetic control: non-functional resistance and parental effects.

Gene drive is a natural process of biased inheritance that, in principle, could be used to control pest and vector populations. As with any form of pest control, attention should be paid to the possibility of resistance evolving. For nuclease-based gene drive aimed at suppressing a population, resistance could arise by changes in the target sequence that maintain function, and various strategies have been proposed to reduce the likelihood that such alleles arise. Even if these strategies are successful, it is almost inevitable that alleles will arise at the target site that are resistant to the drive but do not restore function, and the impact of such sequences on the dynamics of control has been little studied. We use population genetic modelling of a strategy targeting a female fertility gene to demonstrate that such alleles may be expected to accumulate, and thereby reduce the reproductive load on the population, if nuclease expression per se causes substantial heterozygote fitness effects or if parental (especially paternal) deposition of nuclease either reduces offspring fitness or affects the genotype of their germline. All these phenomena have been observed in synthetic drive constructs. It will, therefore, be important to allow for non-functional resistance alleles in predicting the dynamics of constructs in cage populations and the impacts of any field release.

Self-limiting population genetic control with sex-linked genome editors.

In male heterogametic species the Y chromosome is transmitted solely from fathers to sons, and is selected for based only on its impacts on male fitness. This fact can be exploited to develop efficient pest control strategies that use Y-linked editors to disrupt the fitness of female descendants. With simple population genetic and dynamic models we show that Y-linked editors can be substantially more efficient than other self-limiting strategies and, while not as efficient as gene drive approaches, are expected to have less impact on non-target populations with which there is some gene flow. Efficiency can be further augmented by simultaneously releasing an autosomal X-shredder construct, in either the same or different males. Y-linked editors may be an attractive option to consider when efficient control of a species is desired in some locales but not others.

The use of driving endonuclease genes to suppress mosquito vectors of malaria in temporally variable environments.

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.

How driving endonuclease genes can be used to combat pests and disease vectors.

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.

Vector control with driving Y chromosomes: modelling the evolution of resistance.

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.

A synthetic homing endonuclease-based gene drive system in the human malaria mosquito.

Genetic methods of manipulating or eradicating disease vector populations have long been discussed as an attractive alternative to existing control measures because of their potential advantages in terms of effectiveness and species specificity. The development of genetically engineered malaria-resistant mosquitoes has shown, as a proof of principle, the possibility of targeting the mosquito's ability to serve as a disease vector. The translation of these achievements into control measures requires an effective technology to spread a genetic modification from laboratory mosquitoes to field populations. We have suggested previously that homing endonuclease genes (HEGs), a class of simple selfish genetic elements, could be exploited for this purpose. Here we demonstrate that a synthetic genetic element, consisting of mosquito regulatory regions and the homing endonuclease gene I-SceI, can substantially increase its transmission to the progeny in transgenic mosquitoes of the human malaria vector Anopheles gambiae. We show that the I-SceI element is able to invade receptive mosquito cage populations rapidly, validating mathematical models for the transmission dynamics of HEGs. Molecular analyses confirm that expression of I-SceI in the male germline induces high rates of site-specific chromosomal cleavage and gene conversion, which results in the gain of the I-SceI gene, and underlies the observed genetic drive. These findings demonstrate a new mechanism by which genetic control measures can be implemented. Our results also show in principle how sequence-specific genetic drive elements like HEGs could be used to take the step from the genetic engineering of individuals to the genetic engineering of populations.