How driving endonuclease genes can be used to combat pests and disease vectors
© Godfray et al. 2017
Published: 11 September 2017
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.
KeywordsGene drive Gene editing Endonucleases CRISPR-Cas9 Vector control Pest control Mosquitoes
What gene drive is
In the absence of selection a particular gene on a particular chromosome has a 50:50 chance of being present in an individual gamete. Gene drive occurs when a gene is able in some way to subvert Mendelian segregation so that it is overrepresented in the gametes. There is currently great interest in the possibility of using genetic technologies to drive genes through populations of harmful insect pests or vectors [1–8]. The aim of some interventions is to reduce insect densities (population suppression) while the purpose of others is to bring about a genetic change (population replacement) so that, for example, an insect is no longer capable of transmitting a disease. Ideas about utilising gene drive in pest and vector management date back to the 1960s when Hickey and Craig [9, 10] and Curtis  explored the possibility of transforming mosquito populations using self-spreading chromosomal variants, and when the influential evolutionary biologist WD Hamilton  presented models showing that genes that spread through non-Mendelian mechanisms may cause population suppression or elimination.
Many examples of natural gene drive were discovered and investigated in the 20th century , and sporadic discussion about applied gene drive continued in the evolutionary and applied biology literatures, but serious work on operationalizing the technology awaited the new century. A strong stimulus was provided by the discovery of homing endonuclease genes (HEGs) in single-celled eukaryotes , and the suggestion by Burt  that they might be used to control the major vectors of malaria and other human diseases. More recently, the discovery of the CRISPR-Cas9 system in bacteria [16, 17] and the realisation that it may be adapted to create an artificial drive mechanism in eukaryotes analogous to HEGs (zinc-finger nucleases and TALENs might similarly be employed) [3, 18] has brought many new laboratories into the field (see also Additional file 1: Note 1). At the same time the unexpected emergence of the Zika virus in South America  has underscored the relatively restricted armoury we have to defend ourselves against vector-borne diseases and stimulated work on novel approaches.
Gene drive works through competition amongst alleles at the level of the gene rather than the individual (hence they are often called ‘selfish genes’). That it works is demonstrated by mathematical models [15, 20] as well as by observations of nature  and, increasingly, proof of principle in the laboratory [22–24]. However, the logic behind the spread can be counterintuitive and alien to biologists familiar with thinking of selection acting primarily through the differential survival and reproduction of individuals, while the technical literature is difficult to access without a modelling background. The aim of this review is to summarise the theoretical population biology relevant to the application of gene drive using endonucleases for pest and vector control. While providing an entry into the modelling literature, we aim to do this without anything but the minimum of mathematics, relying chiefly on intuitive and graphical arguments. We concentrate on what we call driving endonuclease genes (DEGs), which include HEGs, modified CRISPR-Cas9 systems and other existing and likely-to-be-discovered constructs that spread in the same general way, and we chiefly focus on population suppression and replacement (see Additional file 1: Note 2 for a further way of using gene drive). We do not discuss in any detail the molecular biology of DEGs (for reviews see [3, 4, 13]) or non-DEG drive mechanisms [1, 4]. Finally, we focus on applications of DEGs to address problems arising from insect pests and vectors, the group that has attracted most attention so far, although nearly all the concepts explored here apply equally to gene drive in other sexually reproducing organisms.
Why driving endonuclease genes spread
The normal intuition that a gene spreads because it provides some advantage to the individual that carries it (typically more offspring) fails in the case of non-Mendelian genes such as DEGs, which break the normal rules of inheritance. One now has to think of the gene in competition with its alternative ‘allele’—the homologous stretch of the chromosome that does not carry a DEG. A DEG will spread when rare if an arbitrarily chosen copy of the DEG produces more copies of itself than an arbitrarily chosen example of the alternative allele . Because the DEG is initially rare, it will nearly always be in a heterozygote and the fraction of offspring carrying the DEG will be ½(1 + e), where e is the so-called homing frequency. The DEG gets into half the offspring by simple Mendelian inheritance and into a fraction e of the rest of the offspring by the ‘homing’ action of the endonuclease. A randomly chosen example of the alternative allele will nearly always be in a homozygote (because the DEG is rare) and be transmitted through Mendelian inheritance to ½ the offspring. Because ½(1 + e) is greater than ½, a DEG that has no deleterious effect on its host will always spread, and will eventually become fixed (it excludes the other allele and achieves a frequency of 1).
Natural examples of DEGs tend to be found in self-splicing introns or inteins and probably have few effects on their host’s phenotype . Values of the homing frequency (e) can approach 1 in these systems and this provides an extraordinarily strong selective advantage for the DEG allele. It rapidly goes to fixation where, in the absence of any further selection, it would accumulate mutations and become non-functional. DEGs can only persist in nature by jumping between species, which is probably why they are only found in unicellular organisms with a relatively unprotected germline [21, 26]. The aim of population replacement strategies is to use artificially engineered DEGs to drive a useful construct through a wild population, for example a sequence coding for a peptide that interferes with disease transmission. In the ideal case, this would have no fitness consequences for the host, and it would spread rapidly to fixation.
Population suppression strategies seek to reduce the density of a pest or vector, possibly to a level when the population can no longer sustain itself (population elimination). Unlike population replacement where ideally the DEG does not affect fitness, now a DEG is designed that does cause substantial harm to the population. Consider the case where the DEG is inserted in a chromosome in the middle of a functional gene whose activity it disrupts. Focus first on a fully dominant gene where functional homozygotes and heterozygotes have the same fitness (which we can define to be 1) but non-functional homozygotes have fitness 1 – s. Will the DEG spread? Let us suppose that the DEG acts after the functional gene is expressed so that the conversion of a heterozygote to a non-functional homozygote does not affect the fitness of the individual in which it occurs. When rare, the DEG is nearly always in a heterozygote and its fitness is near the wild type because the functional gene is fully dominant. It will thus begin to spread through the population for the reasons described above.
As the DEG advances through the population more and more non-functional homozygotes will be produced and this will clearly act to slow its spread. We can ask whether the DEG reaches an equilibrium frequency which we denote q. At equilibrium, by definition, the number of copies produced by arbitrarily chosen individual DEG and wild-type alleles should be equal—if this were not true then the frequency would change. An exactly equivalent condition is that the net costs and benefits for each allele due to the segregating DEG should be the same at the equilibrium frequency q (Additional file 1: Note 3). Consider first a wild-type allele: as a homozygote it is unaffected by the DEG but with probability q it is in a heterozygote and suffers the risk of conversion e (so the net effect is –qe). Now the DEG will be in a heterozygote with probability (1 – q) and gain the benefit of conversion with probability e, while with probability q it will be in a homozygote and suffer fitness costs s (so the net effect is (1 − q)e – qs). Equating the net effects for the two alleles and solving for q we obtain q = e/s.
The fact that a gene that reduces population fitness can spread is the basis for the strategy of population suppression. In this simple case of a dominant functional gene it is straightforward to calculate the genetic load the spread of the DEG places on the population . At equilibrium, the frequency of homozygote DEGs is q 2: relative population fitness is thus reduced from 1 to 1 – q 2 s. In the case where the DEG goes to fixation (s ≤ e) and so q = 1 then population fitness is 1 – s. The load is thus the same as the fitness reduction of the DEG homozygote. Where there is a polymorphism (and q = e/s) the load is q 2 s = e 2/s. The potential power of a potent DEG to reduce population fitness is thus very great. For example, were homing to be absolute (e = 1) then DEGs targeting genes that reduced population fitness to near zero (s → 1) would spread, leading to certain elimination.
We have dwelt at some length on this simplest of cases to try to provide an intuitive understanding of why DEGs spread. It is straightforward to relax the simplifying assumptions, though typically mathematical reasoning has to replace verbal arguments, and in Additional file (and for full details see ) we explore some of these issues. Potential candidate genes to be targets for DEGs include viability or fecundity genes that may be active in both sexes or just one. Theory suggests female fecundity genes may be the most effective single genes to target (Fig. 2; Additional file 1: Note 4). DEGs can still spread if they target genes where heterozygote individuals have reduced fitness, and mild heterozygote and homozygote costs are not a barrier for population replacement strategies. High heterozygote costs can lead to complex dynamic behaviours [20, 27, 28] and such genes are poorer targets when population replacement is the goal (Fig. 2; Additional file 1: Note 5). Finally, for population suppression there is a clear advantage in choosing a DEG that homes after the developmental stage during which the target gene is expressed so that the costs imposed by the DEG do not reduce the boost it gets from converting heterozygotes to homozygotes (Additional file 1: Note 6). This advantage is much less important for population replacement strategies and several recent models motivated by this type of intervention assume homing occurs early in development [28, 29].
It is also possible to explore the speed with which a DEG spreads through a population, an important question in designing deployment strategies (Additional file 1: Note 7). In Fig. 2a at the time required for a selection of DEGs to spread is plotted, which shows that rapid increases in frequency can occur over a relatively small number of generations.
Where population suppression is being attempted, multiple DEGs targeting different essential genes can be introduced into a population at the same time. The dynamics of these DEGs is largely independent (homing acts to break up linkage associations between alleles even if they are quite close together on the chromosome). Where each DEG reduces population fitness by an amount L (the load) the overall reduction is 1 – (1 – L 1) (1 – L 2) … (1 – L n ), where the subscript indexes the load of each of the n DEGs. Thus, even if any individual load is quite small, the combination of multiple drags on population fitness quickly adds up [30, 31].
Population genetic arguments can tell you the extent to which a DEG can reduce average population fitness (the genetic load, L, expressed as a fraction between 0 and 1). But for population suppression we need to know what this means for population densities and in particular whether the load is sufficient to cause population elimination, if this is the goal. Answering this question requires an understanding of the target insect’s ecology in the field, often a more challenging task than measuring genetic parameters in the laboratory.
Begin by assuming the insect lives in a constant environment and the DEG causes a substantial reduction in density such that the species is no longer subject to factors such as competition for resources that limit its density when common. Further, assume that the insect has discrete generations and that its capacity to increase when rare is R m , which is measured by the number of female offspring produced per female parent (the argument for insects with overlapping generations is essentially the same). The species will persist if, on average, every female more than replaces itself: R m >1. In the simplest case the genetic load can be measured by the reduction in the number of female offspring produced. Thus, once the DEG is established the growth rate of the population is R m (1 − L) and for elimination to occur L > 1 – 1/R m .
A number of conclusions follow immediately from this simple calculation. First, it will be much harder to eliminate a species with a great capacity to increase when rare. Second, population reduction may work best in peripheral and marginal populations where by definition the species has a lower R m : thus, one might get population elimination in these areas but suppression with persistence in the core of the range. Third, any other control measure that reduced R m at the same time as the DEG would increase the chance of elimination; provided the control measure does not differentially target individuals carrying the new construct, they should help and not hinder control. Fourth, where the load is not large enough to lead to inevitable elimination, the population may be so small that it is subject to demographic stochastic elimination—for example by ‘good’ luck all individuals may fail to breed and so the population is eliminated.
It is surprising how little we know about the magnitude of R m , even for medically and economically important insects. For the main mosquito vector of malaria in Africa (the Anopheles gambiae complex), Deredec et al.  found only one good estimate (from the famous Garki project  in the early 1970s in Nigeria) where R m ~ 9. This suggests that a load of about 0.9 (a 90% reduction in fitness) would be required for the elimination of a population with similar ecology, though R m will certainly vary across sites.
Of course, vectors and pests are not found in constant environments where R m remains constant throughout the year. However, the above arguments still hold, to a good approximation, if R m is replaced by its geometric (rather than arithmetic) mean over time (Additional file 1: Note 8) . For some species there may be a minimum population density below which elimination always occurs (an Allee effect ; Additional file 1: Note 9), which will tend to increase the likelihood of population elimination.
A DEG that substantially reduces population density, but does not cause elimination, would still play a valuable role in reducing the damage done by a pest or vector. However, it is hard to say how the DEG will affect population size because this requires an understanding of how it interacts with all the other ecological factors determining population size and, in particular, those that are density dependent—that is, they vary in magnitude as population density changes. Density-dependent factors are of particular importance because they determine the typical abundance of a species. The relative position in the insect’s life cycle where density dependence acts and where the organism suffers from the presence of the DEG can be very important in determining the effectiveness of the genetic intervention.
The important point here is that it matters where in the lifecycle the DEG acts relative to where density-dependent mortality occurs and the life stage that causes the damage. To get maximum effect, the DEG should act after density dependence but before the damaging stage, otherwise the density-dependent mortality will compensate for some of the load imposed by the DEG. In the major mosquito vectors of human disease it is generally thought that density dependence occurs in the larval stage [35–37] and obviously it is the adults that transmit disease. A DEG acting in the pupal stage or immediately after the adult emerges would be optimal. A pest insect that causes damage in an early instar and where density dependence occurs later in the lifecycle might be best targeted by a DEG that acted at the egg stage.
The density-dependent processes that determine typical population abundances may operate by affecting mortality (as in the example above) or fecundity—perhaps females in high density populations have insufficient food to produce a full complement of eggs, or grow up small and stunted with reduced fecundity. Density-dependent effects may also increase generation time and so slow population growth rates. The precise shape of the density-dependent mortality function can also have important consequences for population dynamic behaviour  (Additional file 1: Note 10). The rather artificial example described above where exactly k individuals emerge from a life stage is termed perfectly compensating density dependence, but more typically the effects of density operate more smoothly without such a sharp threshold (under-compensation). Over-compensation is said to occur when low population densities in one generation lead to very large population densities in the next generation and vice versa. Were such dynamics to characterise mosquito populations then there would be concern that introducing a DEG that did not cause elimination might under some circumstances actually increase equilibrium populations , but what we know about these insects’ biology suggests that this is just a theoretical curiosity. Finally, studies of real populations have to take into account that both density-dependent and density-independent factors vary over time and space. Implementing population suppression successfully will thus require careful study of the population ecology of the target species.
Eliminating or reducing the density of a target species will affect the dynamics of its prey, predator and competitors, and possibly the further species with which they interact. There has been little theoretical exploration of this subject to date in a specific DEG context, though community ecology provides a rich tool box to explore this question (Additional file 1: Note 11).
Many of the possible targets of DEGs are of concern because of the human, animal or plant diseases they transmit, rather than the direct harm they themselves inflict. Using DEGs to drive a gene that interferes with disease transmission through a vector population is a major goal of population replacement strategies, while population suppression aims to target the disease by lowering vector numbers rather than competence. In exactly the same way that the insect population is eliminated if a load is imposed so that each adult female has, on average, less than one female offspring, so the disease is eliminated if an initial infection leads to less than one secondary infection when the pathogen is uncommon. The number of secondary infections when the disease is rare is called the basic epidemiological number or R 0 .
Most expressions that have been derived for the R 0 of vector-borne diseases [41, 42] are the product of (i) biting insects per host and (ii) a parameter combination that includes propensity to bite, adult insect lifespan, and transmission efficiency (see Additional file 1: Note 12 for a worked example). Thus, any vector population suppression will lead to a proportionate reduction in R 0 through the first component, while population replacement strategies affect R 0 through the second. Obviously, vector elimination will lead to disease elimination (in the absence of alternative vector species), but are there circumstances when fully competent vectors persist (R m in the presence of the DEG > 1) but their abundance is so reduced that the disease goes extinct (R 0 < 1)? In a model motivated by the mosquitoes that transmit malaria, Deredec et al.  found circumstances where this might occur (high R m but low R 0), though further work including such factors as seasonality and other interventions in the model is needed to assess its relevance to disease control.
For an insect vector to transmit disease it must live long enough to bite a host and acquire the pathogen, for the pathogen to develop to a life stage in which it can be transmitted, and then for the vector to bite a second host and transmission to occur . Given the high rates of daily mortality suffered by most small invertebrates, it is likely that only a relatively small fraction of long-lived individuals is responsible for disease transmission. A DEG (or other) intervention that reduced adult vector longevity might be more effective in reducing disease incidence than its effects on vector numbers might suggest [43–45]. A DEG that was expressed and caused mortality exclusively in old individuals could be particularly effective, not only because it would kill the dangerous fraction of the population but also because if the numbers involved were relatively small then natural selection for resistance (see next section) would be lower .
One of the consequences of reduced population size and relaxed competition is that larvae may experience reduced competition for food and emerge as larger adults. If larger adults live longer then they could be more efficient at transmitting the disease (the opposite of the process described in the last paragraph), reducing or conceivably reversing the effect of the intervention . However, the situation is more complicated as smaller females may feed more often so that they may be more efficient vectors than larger individuals. Whether this occurs depends critically on the shape of the relevant density-dependent functions. More work is needed but our (unpublished) modelling suggests that quite extreme forms of prolonged lifespan are required for the size-longevity effect to counter the advantages of fewer vectors.
Resistance and recall
Any DEG that reduces the fitness of its host will face the potential evolution of resistance. In this, DEGs are no different from any other intervention that seeks to kill or impair harmful vectors and hosts.
Perhaps the most likely type of resistance allele to arise is one that does not contain the DEG recognition site. Such an allele might be present undetected at very low frequency or may arise de novo after the DEG has been released, either due to background mutation or through the action of homing itself. Homing works because double-strand breaks are usually repaired using the homologous chromosome as a template. But repairs can also be made by end joining (EJ), which as the name suggests involves direct ligation of sundered chromosome ends . If EJ leads to a perfect repair in which the DEG recognition site is completely reconstructed then there is no lasting affect and EJ can be modelled simply as a reduction in homing frequency. But if EJ destroys the DEG recognition site then an allele resistant to the DEG will have been created. Irrespective of origin, the fate of such an allele will depend on its fitness relative to the different DEG genotypes . Other types of resistance are possible; for example, any mutation not in the DEG that resulted in the disabling or blocking of the DEG could be favoured. Resistance mutations generated by EJ, including non-homologous and micro-homology-mediated EJ , have been observed in the laboratory .
These considerations suggest a series of strategies to try to avoid resistance. The probability of functional escape mutations can be reduced by ensuring that the DEG recognition site is in a conserved region of the gene where any change to the recognition site also destroys gene function. Third-base pair redundancy of course makes this difficult, but not impossible—in mosquitoes some sequences are, for unclear reasons, conserved at the nucleotide level for tens of millions of years . A critical issue is to ensure that no alleles immune to the DEG are already segregating in the target population. Genetic surveys can detect such alleles if they are at high frequency but will be much less likely to find rare variants. Again, choosing recognition sites in conserved regions of genes is likely to be the best precautionary strategy.
A further way to avoid resistance is to combine multiple DEGs in a single intervention, or to use DEGs that cause multiple double-strand chromosome breaks [15, 29, 31]. As in combination drug therapy, the probability of resistance can be reduced by requiring multiple resistance mutations to occur simultaneously. Additional file 1: Note 14 describes an elegant analysis by Marshall et al.  suggesting that resistance to single CRISPR-like DEGs is almost inevitable (if it has no fitness cost), while resistance to multiple DEGs becomes exponentially more difficult as their numbers increase. Two aspects of DEGs facilitate this approach. First DEGs, especially those based on CRISPR-Cas9, can be reprogrammed to target multiple sites . Such a strategy would fail if a mutation arose that suppressed all CRISPR activity, but then the second advantage, the availability of very different DEG systems (or CRISPR variants) that are unlikely to be suppressed by the same mutation, comes into play. A mutation that suppressed chromosome repair using the homologous chromosome could affect all DEG approaches though likely at the cost of reducing the organism’s overall fitness. Even then, there are other ways to use DEGs that might be ‘stacked’ to reduce the likelihood of resistance.
In the context of population replacement, Unckless et al.  develop a valuable theoretical framework to explore the dynamics of resistance alleles that are already present at low frequencies in a target population or that arise de novo after deployment. They assume that individuals in which homing occurs experience the fitness costs of the homozygote DEG (that is, homing occurs before expression) but that these costs are much lower than for a population suppression strategy. For replacement strategies, minimising fitness costs and targeting conserved regions for recognition sites are both important mitigation strategies. Noble et al.  (Additional file 1: Note 15) model a neat variant of a replacement strategy originally suggested by Esvelt et al. . In their scheme, part of the payload of the construct is a copy of the gene being targeted that has been recoded at the DNA sequence level so that it is immune to cleavage but still produces a viable protein product. Any fitness cost due to recognition site disruption would thus be avoided, reducing selection for resistance. Finally, when population replacement seeks to render vector individuals incapable of transmitting a pathogen, the possibility of the latter evolving resistance to the intervention must be considered.
The strength of selection for resistance means that any artificially created resistance allele would quickly spread through the population. It would thus be possible to use natural selection to reverse the effect of any costly DEG should it be desired . Such a remediation strategy might be requested by regulators, or could be used to reinstate an insect vector population after the disease it carries has gone extinct. The reconstituted population might not be genetically identical to the original wild type, but the differences could be minimised to the small number of bases in the recognition sequence. Related to this, population replacement might be achieved via a population suppression intermediate step [50, 51]. The idea here is first to impose a load on the population and then to introduce an artificial resistance allele in combination with a beneficial construct. A reason for doing this is that the end result would be the establishment of the desired construct but with the loss of the genetic element capable of drive.
Even when a DEG used for population replacement has no effect on fitness, a DEG could in principle be engineered to reinstate an allele that had been disrupted by a previous gene drive. To do this it would be necessary to prevent the original DEG from cutting the replacement allele, for example by altering the DNA sequence in a way that did not affect the protein sequence, or in a two stage process in which a sequence immune to cutting is first driven through a population followed by the desired allele once the original DEG has been lost . Alternatively, Wu et al.  suggest that a CRISPR-Cas9 DEG system may be recalled by introducing a guide RNA that targets the Cas9 gene itself (they call this mechanism CATCHA after Cas9-triggered chain ablation [of Cas9]; Additional file 1: Note 16; Additional file 2).
These elaborations on the basic theme of DEG-drive are conceptually fascinating and potentially important. Nevertheless, the reality is that regulators will almost certainly need to be convinced that ‘basic’ gene drive is efficient and safe before they countenance more complicated constructs.
There is a second way that a driving endonuclease genes can spread . In species where there are differentiated sex chromosomes the two different chromosomes in the heterogametic sex can be thought of as competing to be represented amongst the gametes. Consider a DEG that is carried by the Y chromosome but which codes for an endonuclease that targets a sequence on the X chromosome. We assume males are the heterogametic sex (XY). If the double-strand break the endonuclease causes on the X is not repaired then by destroying competing gametes it gives the Y chromosome on which it is carried an advantage. The Y chromosome spreads, and more and more individuals will come to carry it; the increasing preponderance of Y over X gametes produced by males carrying the DEG will lead to a population-wide biased sex ratio. This may be good in itself, for example in mosquitoes where only females transmit human diseases, but depending on the extent of the bias and the details of the species’ ecology it may lead to population suppression or even elimination [30, 53]. A number of examples of ‘driving Y chromosomes’ leading to male-biased sex ratios have been observed in nature, though the mechanisms involved are not known . Proof of principle of artificial DEG Y-drive has been shown in the laboratory with mosquitoes [54, 55].
The reason a driving Y spreads can be understood by comparing the number of offspring to which an arbitrarily chosen wild-type and DEG-carrying Y chromosome can expect to be transmitted . Begin by assuming the mating competitiveness and capacity of males carrying the two types of Y chromosome are the same. The wild-type Y will be present in half the offspring that its bearer sires. Because we are assuming equal competitiveness, a male carrying the modified Y will father the same number of offspring (the effect of relaxing this assumption is explored in Additional file 1: Note 17) but the Y:X ratio of its gametes will be 1:1 − e where e is the probability of successful cleavage (and hence the fraction of X chromosomes destroyed). The modified Y can thus expect to be present in 1/(2 – e) of the offspring, which is always greater than ½. In this simple case, a driving Y chromosome invariably spreads all the way, it can be shown, to fixation.
One of the attractions of gene drive for pest and vector control is that they spread across the landscape: unlike other control methods, such as insecticides or sterile-insect releases, they do not need to be applied everywhere. Modelling is important to give insights into the speed of spread and for the ability of the DEG to spread in spatially heterogeneous regions.
A critical determinant of spatial spread is the extent to which individual insects move in their lifetime. The simplest way to conceptualise this is that an insect that emerges as an adult in one place will move during its lifespan to different places in the landscape with a probability that can be described by a normal distribution centred at its birthplace. The variance of this distribution is proportional to what is called the diffusion coefficient, by analogy to diffusion in physical processes that can be described in the same way . Implicit in this description is the assumption that there is no directional movement or type of behaviour, such as rare long-distance jumps, that would give rise to non-normal distributions. Despite its biological simplicity, it is a good place to start to understand spatial spread.
A simple driving DEG that spreads through homing is also likely to spread through a homogeneous landscape at a speed proportional to the square root of the product of average insect movement and the DEG’s net advantage, though this has not been formally demonstrated. As described in Additional file 1: Note 5, some types of DEG are only spread when their frequency exceeds a threshold and here the pattern of spatial spread is through a ‘Bartonian’ rather than a ‘Fisherian’ wave  (Additional file 1: Note 20).
The assumption of homogenous populations is clearly unrealistic and results from spatial population ecology will be helpful in analysing these more realistic situations (Additional file 1: Note 21). North et al. , for example, developed a stochastic individual-based model of DEG spread (either homing or Y drive) motivated by the biology of mosquitoes. Adults move through the landscape searching for larval oviposition sites or humans on which to feed. When these resources are relatively abundant and uniformly distributed the dynamics of the system can be approximated by non-spatial or spatially homogenous models. When resources are scarcer, an introduced HEG can eliminate a local population and then itself go extinct while wild-type refugia remain. The chance of this happening increases as adult feeding sites become scarce and as the covariance between adult feeding and breeding sites increases: both effects increase the autonomy of the local populations. After the DEG goes extinct the wild-type populations may then recolonize the landscape. The risk of wild-type refuges remaining can be substantially reduced by making multiple releases of DEG-carrying mosquitoes across the environment.
The future of gene-drive technology
Modelling of gene drive is still relatively in its infancy. Beginning with stylised population genetic and then population dynamic models, theoretical work has expanded to begin to include spatial processes and the ecology of particular species at specific locations. Models have also been developed to explore new tactical possibilities presented by advances in molecular biology, in particular the discovery of the CRISP-Cas9 system.
We believe it highly likely within the next decade that a driving element based on a DEG will be available for deployment to transform or suppress the population of an injurious arthropod, most likely a mosquito vector of an important human disease. Whether deployment actually occurs will depend on whether regulators are satisfied the technology is safe to humans and poses few environmental risks. Modelling will play an important role in submissions to regulators as the genetic and population dynamic consequences of DEG release can be explored in no other way. Deployment will also require a ‘licence to operate’ from civil society, and for citizens to engage profitably in discussions it is critical for modelling outputs to made available in the most accessible forms. This review concentrates on the population biology of gene drive but the same endonucleases that can be used in gene drive may also be used in contexts where spread does not occur or is limited, for example in genetic versions of the sterile insect technique (Additional file 1: Note 22). It is likely that these will be the first modified endonucleases permitted to be released.
The field has changed enormously in the mere five years since CRISPR-Cas9 burst onto the scene  and it would be rash to predict how gene drive technology may be affected by new molecular advances in the coming years. Conceptually, CRISPR-Cas9 DEGs function largely in the same way as the first generation of homing endonucleases and much of the theory developed for this first generation of DEGs is equally applicable. But the ease of manipulation of CRISPR systems is transformational and has brought new laboratories and new funds into the field, as well as new suggestions about how DEGs may be deployed. It has also led to greatly increased press and public attention, and the field, rightly, will be expected to be transparent and responsible . This is a uniquely exciting time to work on the genetic control of pests and vectors.
We thank Anne Deredec, Philip Eckhoff, Ben Lambert and the members of the Target Malaria Consortium for valuable discussion. The authors are supported by a grant from the Bill & Melinda Gates Foundation.
HCJG drafted the review and all authors worked on the final version. All authors read and approved the final manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Sinkins SP, Gould F. Gene drive systems for insect disease vectors. Nat Rev Genet. 2006;7(6):427–35.View ArticlePubMedGoogle Scholar
- Burt A. Heritable strategies for controlling insect vectors of disease. Philos Trans R Soc B Biol Sci. 2014;369(1656):20130432.View ArticleGoogle Scholar
- Esvelt KM, Smidler AL, Catteruccia F, Church GM. Concerning RNA-guided gene drives for the alteration of wild populations. Elife. 2014;3:e03401.Google Scholar
- Champer J, Buchman A, Akbari OS. Cheating evolution: engineering gene drives to manipulate the fate of wild populations. Nat Rev Genet. 2016;17(3):146–59.View ArticlePubMedGoogle Scholar
- Gantz VM, Bier E. The dawn of active genetics. Bioessays. 2016;38(1):50–63.View ArticlePubMedGoogle Scholar
- Alphey L. Can CRISPR-Cas9 gene drives curb malaria? Nat Biotechnol. 2016;34(2):149–50.View ArticlePubMedGoogle Scholar
- Adelman ZN, Tu ZJ. Control of mosquito-borne infectious diseases: sex and gene drive. Trends Parasitol. 2016;32(3):219–29.View ArticlePubMedPubMed CentralGoogle Scholar
- Leftwich PT, Bolton M, Chapman T. Evolutionary biology and genetic techniques for insect control. Evol Appl. 2016;9(1):212–30.View ArticlePubMedGoogle Scholar
- Craig GB, Hickey WA, van de Hey RC. An inherited male-producing factor in Aedes aegypti. Science. 1960;132:1887–9.View ArticlePubMedGoogle Scholar
- Hickey WA, Craig GB. Genetic distortion of sex ratio in a mosquito, Aedes aegypti. Genetics. 1966;53:1177–96.PubMedPubMed CentralGoogle Scholar
- Curtis CF. Possible use of translocations to fix desirable genes in insect pest populations. Nature. 1968;218:368–9.View ArticlePubMedGoogle Scholar
- Hamilton WD. Extraordinary sex ratios. Science. 1967;156:477–88.View ArticlePubMedGoogle Scholar
- Burt A, Trivers R. Genes in conflict. Cambridge: Belknap Press; 2006.View ArticleGoogle Scholar
- Stoddard BL. Homing endonuclease structure and function. Q Rev Biophys. 2005;38(1):49–95.View ArticlePubMedGoogle Scholar
- Burt A. Site-specific selfish genes as tools for the control and genetic engineering of natural populations. Proc Biol Sci. 2003;270(1518):921–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Jinek M, Chylinski K, Fonfara I, Hauer M, Doudna JA, Charpentier E. A programmable dual-RNA-guided DNA endonuclease in adaptive bacterial immunity. Science. 2012;337(6096):816–21.View ArticlePubMedGoogle Scholar
- Doudna JA, Charpentier E. The new frontier of genome engineering with CRISPR-Cas9. Science. 2014;346(6213):1258096.View ArticlePubMedGoogle Scholar
- Simoni A, Siniscalchi C, Chan YS, Huen DS, Russell S, Windbichler N, et al. Development of synthetic selfish elements based on modular nucleases in Drosophila melanogaster. Nucleic Acids Res. 2014;42(11):7461–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Fauci AS, Morens DM. Zika Virus in the Americas - yet another Arbovirus threat. New Engl J Med. 2016;374(7):601–4.View ArticlePubMedGoogle Scholar
- Deredec A, Burt A, Godfray HCJ. The population genetics of using homing endonuclease genes in vector and pest management. Genetics. 2008;179(4):2013–26.View ArticlePubMedPubMed CentralGoogle Scholar
- Goddard MR, Burt A. Recurrent invasion and extinction of a selfish gene. Proc Natl Acad Sci U S A. 1999;96(24):13880–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Windbichler N, Menichelli M, Papathanos PA, Thyme SB, Li H, Ulge UY, et al. A synthetic homing endonuclease-based gene drive system in the human malaria mosquito. Nature. 2011;473(7346):212–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Hammond A, Galizi R, Kyrou K, Simoni A, Siniscalchi C, Katsanos D, et al. A CRISPR-Cas9 gene drive system-targeting female reproduction in the malaria mosquito vector Anopheles gambiae. Nat Biotechnol. 2016;34(1):78–83.View ArticlePubMedGoogle Scholar
- Gantz VM, Jasinskiene N, Tatarenkova O, Fazekas A, Macias VM, Bier E, et al. Highly efficient Cas9-mediated gene drive for population modification of the malaria vector mosquito Anopheles stephensi. Proc Natl Acad Sci U S A. 2015;112(49):E6736–43.View ArticlePubMedPubMed CentralGoogle Scholar
- Burt A, Koufopanou V. Homing endonuclease genes: the rise and fall and rise again of a selfish element. Curr Opin Genet Dev. 2004;14(6):609–15.View ArticlePubMedGoogle Scholar
- Koufopanou V, Goddard MR, Burt A. Adaptation for horizontal transfer in a homing endonuclease. Mol Biol Evol. 2002;19(3):239–46.View ArticlePubMedGoogle Scholar
- Unckless RL, Messer PW, Connallon T, Clark AG. Modeling the manipulation of natural populations by the mutagenic chain reaction. Genetics. 2015;201(2):425–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Unckless RL, Clark AG, Messer PW. Evolution of resistance against CRISPR/Cas9 gene drive. Genetics. 2017;205(2):827–41.View ArticlePubMedGoogle Scholar
- Noble C, Olejarz J, Esvelt K, Church G, Nowak M. Evolutionary dynamics of CRISPR gene drives. Sci. Advances. 2017;3(4):e1601964.Google Scholar
- Deredec A, Godfray HCJ, Burt A. Requirements for effective malaria control with homing endonuclease genes. Proc Natl Acad Sci U S A. 2011;108(43):E874–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Marshall JM, Buchman A, Sánchez HM, Akbari OS. Overcoming evolved resistance to population-suppressing homing-based gene drives. Sci Rep-Uk. 2016;in press.Google Scholar
- Molineaux L, Gramiccia G. The Garki Project. Geneva: World Health Organisation; 1980.Google Scholar
- Tuljapurkar S. Population dynamics in variable environments. Berlin: Springer-Verlag; 1990.View ArticleGoogle Scholar
- Stephens PA, Sutherland WJ, Freckleton RP. What is the Allee effect? Oikos. 1999;87(1):185–90.View ArticleGoogle Scholar
- Gimnig JE, Ombok M, Otieno S, Kaufman MG, Vulule JM, Walker ED. Density-dependent development of Anopheles gambiae (Diptera: Culicidae) larvae in artificial habitats. J Med Entomol. 2002;39(1):162–72.View ArticlePubMedGoogle Scholar
- White MT, Griffin JT, Churcher TS, Ferguson NM, Basanez MG, Ghani AC. Modelling the impact of vector control interventions on Anopheles gambiae population dynamics. Parasites Vectors. 2011;4:153.View ArticlePubMedPubMed CentralGoogle Scholar
- Muriu SM, Coulson T, Mbogo CM, Godfray HCJ. Larval density dependence in Anopheles gambiae s.s., the major African vector of malaria. J Anim Ecol. 2013;82(1):166–74.View ArticlePubMedGoogle Scholar
- May RM. Stability and complexity in model ecosystems. Princeton: Princeton University Press; 1974.Google Scholar
- Alphey N, Bonsall MB. Interplay of population genetics and dynamics in the genetic control of mosquitoes. J R Soc Interface. 2014;11(93):20131071.View ArticlePubMedPubMed CentralGoogle Scholar
- Anderson RM, May RM. Infectious diseases of humans. Oxford: Oxford University Press; 1991.Google Scholar
- Macdonald G. The epidemiology and control of malaria. Oxford: Oxford University Press; 1957.Google Scholar
- Smith DL, McKenzie FE. Statics and dynamics of malaria infection in Anopheles mosquitoes. Malaria J. 2004;3:Art. No. 13.View ArticleGoogle Scholar
- Cook PE, McMeniman CJ, O'Neill SL. Modifying insect population age structure to control vector-borne disease. Adv Exp Med Biol. 2008;627:126–40.View ArticlePubMedGoogle Scholar
- Hancock PA, Thomas MB, Godfray HCJ. An age-structured model to evaluate the potential of novel malaria-control interventions: a case study of fungal biopesticide sprays. Proc Biol Sci. 2009;276(1654):71–80.View ArticlePubMedGoogle Scholar
- Koella JC, Lynch PA, Thomas MB, Read AF. Towards evolution-proof malaria control with insecticides. Evol Appl. 2009;2(4):469–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Lyimo EO, Koella JC. Relationship between body size of adult Anopheles gambiae s.1. and infection with the malaria parasite Plasmodium falciparum. Parasitology. 1992;104:233–7.View ArticlePubMedGoogle Scholar
- Weterings E, Chen DJ. The endless tale of non-homologous end-joining. Cell Res. 2008;18(1):114–24.View ArticlePubMedGoogle Scholar
- Yajima H, Fujisawa H, Nakajima NI, Hirakawa H, Jeggo PA, Okayasu R, et al. The complexity of DNA double strand breaks is a critical factor enhancing end-resection. DNA Repair. 2013;12(11):936–46.View ArticlePubMedGoogle Scholar
- Neafsey DE, Waterhouse RM, Abai MR, Aganezov SS, Alekseyev MA, Allen JE, et al. Highly evolvable malaria vectors: the genomes of 16 Anopheles mosquitoes. Science. 2015;347(6217):1258522.View ArticlePubMedGoogle Scholar
- Beaghton A, Hammond A, Nolan T, Crisanti A, Godfray HCJ, Burt A. Requirements for driving anti-pathogen effector genes into populations of disease vectors by homing. Genetics. 2017;205:1587–96.View ArticlePubMedPubMed CentralGoogle Scholar
- Huang YX, Magori K, Lloyd AL, Gould F. Introducing desirable transgenes into insect populations using Y-linked meiotic drive--a theoretical assessment. Evolution. 2007;61(4):717–26.View ArticlePubMedGoogle Scholar
- Wu B, Luo LQ, Gao XJJ. Cas9-triggered chain ablation of cas9 as a gene drive brake. Nat Biotechnol. 2016;34(2):137–8.View ArticlePubMedPubMed CentralGoogle Scholar
- North A, Burt A, Godfray HCJ. Modelling the spatial spread of a homing endonuclease gene in a mosquito population. J Appl Ecol. 2013;50(5):1216–25.PubMedPubMed CentralGoogle Scholar
- Galizi R, Doyle LA, Menichelli M, Bernardini F, Deredec A, Burt A, et al. A synthetic sex ratio distortion system for the control of the human malaria mosquito. Nat Commun. 2014;5:3977.View ArticlePubMedPubMed CentralGoogle Scholar
- Galizi R, Hammond A, Kyrou K, Taxiarchi C, Bernardini F, O'Loughlin SM, et al. A CRISPR-Cas9 sex-ratio distortion system for genetic control. Sci Rep. 2016;6:31139.Google Scholar
- Windbichler N, Papathanos PA, Catteruccia F, Ranson H, Burt A, Crisanti A. Homing endonuclease mediated gene targeting in Anopheles gambiae cells and embryos. Nucleic Acids Res. 2007;35(17):5922–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Ovaskainen O, de Knegt HJ, del Mar Delgado M. Quantitative ecology and evolutionary biology. Oxford: Oxford University Press; 2016.View ArticleGoogle Scholar
- Fisher RA. The wave of advance of advantageous genes. Ann Hum Genet. 1937;7:355–69.Google Scholar
- Beaghton A, Beaghton PJ, Burt A. Gene drive through a landscape: Reaction-diffusion models of population suppression and elimination by a sex ratio distorter. Theor Popul Biol. 2016;108:51–69.View ArticlePubMedGoogle Scholar
- Barton NH. The dynamics of hybrid zones. Heredity. 1979;43:341–59.View ArticleGoogle Scholar
- Eckhoff PA, Wenger EA, Godfray HCJ, Burt A. Impact of mosquito gene drive on malaria elimination in a computational model with explicit spatial and temporal dynamics. Proc Natl Acad Sci U S A. 2017;114(2):E255–64.View ArticlePubMedGoogle Scholar
- Esvelt K. Gene editing can drive science to openness. Nature. 2016;534(7606):153.View ArticlePubMedGoogle Scholar