#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Double drives and private alleles for localised population genetic control


Authors: Katie Willis aff001;  Austin Burt aff001
Authors place of work: Department of Life Sciences, Imperial College London, Silwood Park, Ascot, United Kingdom aff001
Published in the journal: Double drives and private alleles for localised population genetic control. PLoS Genet 17(3): e1009333. doi:10.1371/journal.pgen.1009333
Category: Research Article
doi: https://doi.org/10.1371/journal.pgen.1009333

Summary

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.

Keywords:

Genetic loci – Guide RNA – Heterozygosity – Islands – Nucleases – Pest control – Population genetics – Population size

Introduction

Gene drive is a natural phenomenon in which some genes are able to increase in frequency and spread through populations by contriving to be inherited at a greater-than-Mendelian rate [1,2]. Strong drive can cause genes to increase rapidly in frequency even if they also harm the organisms carrying them, and there is currently much effort trying to develop synthetic gene drive constructs (or gene drives) to control disease-transmitting mosquitoes and other pest populations that have thus far been difficult or impossible to manage satisfactorily [36]. If a species is harmful and subject to control measures wherever it exists, then, in principle (i.e., in the computer), highly efficient gene drive strategies can be devised that exploit natural processes of dispersal and gene flow such that relatively small inoculative releases in a few locations can lead to substantial and widespread impacts over subsequent generations [79]. However, some species are pests only in a part of their range (e.g., invasive species), and other approaches are needed.

Two broad approaches have been proposed for restricting the impact of genetic control interventions to a target population. First, one can use a strategy requiring relatively large releases, which can be restricted to the target population, with any introductions into non-target populations (by dispersal, or by accidental or unauthorised releases) being too small to have a significant impact. Potentially suitable genetic constructs include those that do not drive (e.g., dominant lethals, autosomal X-shredders, or Y-linked editors; [1012]), those that show transient drive due to a non-driving helper construct (e.g., killer-rescue systems and split drives; [1315]) or those that drive, but only if they are above some threshold frequency (e.g., various underdominant [heterozygote inferiority] strategies, tethered drives, and split drive killer-rescue systems [1619]). Some of these approaches are more efficient than others [12,20,21], but, by necessity, all of them require a non-trivial production and release effort.

Alternatively, if there are pre-existing sequence differences between target and non-target populations, it may be possible to exploit these differences with a sequence-specific nuclease-based gene drive that would only spread in the target population, in which case the small release rates and overall efficiency of low threshold gene drive approaches may be retained [22,23]. Sudweeks et al. [23] present useful modelling of this approach, considering the case where there is a locally fixed allele of an essential gene in the target population, while non-target populations carry a functional cleavage-resistant allele at some frequency. A single-locus gene drive that uses the homing reaction (i.e., sequence-specific cleavage followed by homologous repair [4,24]) to disrupt the locally fixed allele could be released into and eliminate the target population, but have little impact, or only a transient impact, on non-target populations. However, as emphasised by the authors, if the target population has even a small frequency of the resistant allele, then that allele could be rapidly selected for and the intervention fail. Single locus homing drives targeting an essential gene in order to suppress a population necessarily generate strong selection pressure for resistant sequences if these can arise [24]; one approach to this problem is to target sites where functional resistance is unlikely to arise [25], but this is difficult to engineer if the target site is chosen such that a resistant allele exists at high frequency in non-target populations of the same species.

In this paper we explore alternative two-locus “double drive” low release rate strategies to restrict population control based on pre-existing sequence differences between target and non-target populations. All our designs are based on a division of labour between the two constructs, with one imposing a reproductive load by disrupting a gene needed for survival or reproduction, and therefore responsible for the desired impact (population suppression), and the other responsible for the population restriction. The first construct can then be designed to target a well-conserved essential site where functional resistance is unlikely to arise, and selection for resistant alleles at the differentiated locus will be relatively weaker because it is not directly responsible for the reproductive load. As a result, these designs are substantially less susceptible to pre-existing resistance in the target population at the differentiated locus than single drive designs, and can even work if the differentiated locus is selectively neutral. Double drives may also be useful for population replacement. Finally, analyses of published genome sequences from island and mainland populations of the malaria mosquito Anopheles gambiae indicates that the sort of population differentiation we model can exist in nature.

Results

Simple double drives for population suppression

The simplest double drive designs we consider consist of one construct (call it α) inserted into and disrupting a haplo-sufficient female-essential gene, such that homozygous females die without reproducing while heterozygous females and all males are unaffected, and a second construct (β) inserted into a sequence that is significantly more common in the target than the non-target population(s). Both constructs are able to drive by the homing reaction but α can drive only in the presence of β, while β may either drive autonomously or rely on the presence of α. With CRISPR-based designs, α would encode its cognate gRNA, β would encode the Cas9, and either construct could encode the gRNA for the second locus (Fig 1, Designs 1 and 2). We assume the α construct has been designed such that functional resistance is not possible (e.g. by targeting sequences that are essential at the nucleotide level or by using multiple gRNAs), though non-functional resistant alleles can arise by end-joining repair [2527]. For the β construct we initially suppose its insertion site (i.e., the differentiated locus) is selectively neutral and unlinked with the α insertion site, and that differentiation is nearly complete, with the recognition sequence present at a frequency of 99% in the target population and absent in the non-target population (i.e., it is a virtually fixed private allele).

Fig. 1. Alternative double drive designs for population suppression.
Alternative double drive designs for population suppression.
Constructs α and β can drive autonomously (a) or non-autonomously (n); one or the other may encode an X-shredder; and the A target locus can be a gene that is haplo-sufficient (HS) or haplo-insufficient (HI) for female viability or fertility. fsHS—female-specific haplo-sufficient locus; fsHI—female-specific haplo-insufficient locus; diff—differentiated sequence; X-shr—X-shredder targeting an X-linked repeat.

Under these conditions, a small (0.1%) release of males carrying Design 1 constructs into the target population leads to both constructs rapidly increasing in frequency and, as a result, an increasing fraction of female zygotes are homozygous for α and die without reproducing. The population size crashes to a minimum size of 3.58e-6 (relative to the pre-release equilibrium) after 25 generations (Fig 2A). Depending on the initial population size and the biology of the species (e.g., whether there are Allee effects [28] such that the population cannot persist at small sizes), this decline could be enough to eliminate the population. However, in our simple deterministic model population elimination is not possible. Instead, the population recovers due to selection of cleavage-resistant alleles at the differentiated locus (which either pre-existed or arose due to end-joining repair), leading to loss of the β construct, followed by loss of α, allowing the wild-type allele and population fertility to be restored. By contrast, the same releases into the non-target population have minimal effect: β cannot increase in frequency (because its target site is absent), and therefore α remains rare, and population size is little affected (Fig 2B).

Fig. 2. Performance of double drives for population suppression.
Performance of double drives for population suppression.
(A, B) Timecourse for Design 1 in target and non-target populations, assuming 1% and 100% pre-existing resistance at the differentiated locus, respectively. In the target population the α and β constructs increase in frequency together (blue and red solid lines), causing the number of females to decline. If the population is not eliminated, then eventually the resistant b allele replaces β, followed by the wild-type A allele replacing α, allowing the population to recover. In the non-target population both constructs remain rare and the reduction in female numbers remains small. (C, D) Timecourse for Design 3 assuming 20% and 100% pre-existing resistance in the target and non-target populations. (E, F) Timecourse for Design 5 assuming 50% and 100% pre-existing resistance in the target and non-target populations.

Because the spread of construct α in the target population depends on β, and therefore will be affected by the association between them, it might be expected that close linkage between the two constructs may increase construct spread and the extent of population suppression. Close linkage has been observed to affect the dynamics of other two-locus drive systems [15]. Furthermore, because population recovery (if it occurs) will be due to the evolution of resistance at the differentiated locus, additional improvements might be expected by using an essential gene as the differentiated locus while designing β to have minimal fitness effects (e.g., by containing a recoded, cleavage-resistant version of the essential gene [14,26,29], or by being inserted in an artificial intron [30]). End-joining repair will then tend to produce non-functional resistance alleles, increasing the load, and functional resistance at the differentiated locus will be slower to evolve, relying instead on pre-existing resistant alleles. Moreover, this effect may be stronger if the essential gene is haplo-insufficient than if it is haplo-sufficient, as found with some other gene drive designs [31]. Both these expectations about linkage and using an essential differentiated gene are met individually, and, in combination, can reduce the minimum population size achieved by many orders of magnitude (Fig 3; see also S1 Fig for the separate effect of each modification). If it is not possible to have close linkage, then the maximum level of suppression can also be increased by releasing the two constructs in different males rather than in the same males, which allows β to escape some of the fitness costs of α and get to a higher frequency than it otherwise would, though at the cost of the impact being delayed, and separate releases perform worse than combined releases when linkage is tight (S2 Fig).

Fig. 3. Timecourse for the relative number of females over time for Designs 1 and 2.
Timecourse for the relative number of females over time for Designs 1 and 2.
Solid lines are for β in a neutral locus unlinked to the α construct, and dashed lines for β as a neutral insertion in a haplo-insufficient essential gene closely linked (r = 0.01) to the α construct. In all cases there is 1% pre-existing resistance at the differentiated locus. Also shown for comparison are results for a single construct drive targeting a haplo-sufficient female-specific viability gene, assuming 1% pre-existing functional resistance and that end-joining repair produces only nonfunctional resistant alleles (S).

Design 2, which has the same components as Design 1, but arranged differently such that homing of the β construct only occurs in the presence of α, has dynamics qualitatively similar to Design 1, but quantitatively different (S3 Fig). Interestingly, if the two constructs are unlinked then the extent of suppression is less than with Design 1, but if they are closely linked then the suppression can be greater (Figs 3 and S1). For comparison we also model a single drive homing into a differentiated female-essential gene with 1% pre-existing resistance. The maximum extent of suppression is much less than with any of the double drives considered, because selection for resistance is much stronger, being directly at the fitness-determining locus (Fig 3).

Coping with higher frequencies of pre-existing resistance

Though these simple double drive designs work well with 1% pre-existing target site resistance at the differentiated locus, performance declines rapidly after that. For example, if there is 10% pre-existing resistance, then even the best of these designs (Design 2 with close linkage and the differentiated locus being haplo-insufficient) only suppresses the target population to a minimum of 2.38e-4 (Fig 4). In some situations the target population may not have a private allele with frequency over 90% and alternative approaches would need to be considered. One possibility is to increase the load imposed on the population by the α construct by adding to it an X-shredder locus that destroys the X-chromosome during spermatogenesis such that α now produces a male-biased sex ratio as well as killing homozygous females (Fig 1, Designs 3 and 4). Since population productivity in many species depends on the number of females, population size may thereby be further reduced. A single drive based on these components has previously been constructed in Anopheles gambiae by Simoni et al. [32]. Our modelling indicates that adding an X-shredder to a double drive gives a quantitative improvement in the dynamics, and even pre-existing resistance frequencies of 20% are compatible with good control, while still having minimal effect on non-target populations (Fig 2C and 2D).

Fig. 4. Minimum population size for each of the 5 designs as a function of the pre-existing frequency of resistance at the differentiated locus.
Minimum population size for each of the 5 designs as a function of the pre-existing frequency of resistance at the differentiated locus.
Solid lines are for the baseline case (r = 0.5, β in a neutral locus), while dashed lines are for the improved case (r = 0.01, β as a neutral insertion in a haplo-insufficient essential gene).

Even more robust control can be obtained by adding the X-shredder to the β construct and having the α construct drive autonomously in males and cause dominant sterility or lethality in females (e.g., target a female-specific haplo-insufficient locus; Fig 1, Design 5). As the Cas9 is encoded by α, homing of β also occurs only in males. The dynamics in this case are somewhat different from the others: the X-shredder does not function to directly increase the load, but instead it allows the α construct to spread in the population, because it will end up more often in males (where it homes), and less often in females (where it is a dead end). The male bias also protects the β construct from the female lethality produced by the α construct, and so selection against β is much weaker than in the previous designs, and resistance evolves more slowly (compare the rate of spread of the resistant b allele in Fig 2E to that in Fig 2A and 2C). As a result the design is able to perform well even with pre-existing resistance of up to 50%, but still not spread in the non-target population (Fig 2E and 2F). Moreover, if the population is not eliminated, it can nevertheless be suppressed for many generations. For example, with 50% pre-existing resistance the minimum population size reached is 2.15e-3, and the population remains below 5% of its pre-intervention size for 63 generations; with close linkage (r = 0.01), then the corresponding values are 8.79e-7 and 147 generations. A comparison of the maximum extent of suppression as a function of the pre-existing resistance frequency for the different designs is shown in Fig 4 (see also S4 Fig). Note that none of the modifications considered (linkage, use of an essential differentiated gene, or separate releases) has a qualitative effect on dynamics in the non-target population, as β is still unable to increase in frequency, and impacts on population size remain small (S5 Fig).

Evolutionary stability and impact of fitness costs

We now explore the consequences of relaxing two assumptions that have been implicit thus far in our modelling. First, we have assumed that our various constructs remain intact after release. In fact, mutations that destroy the function of one component or another will be expected to arise as the constructs spread through a population, particularly as homing may be associated with a higher mutation rate than normal DNA replication [3335]. For components that contribute directly to their construct’s spread, one would expect that loss-of-function mutations would remain rare in the population and have little effect, whereas for other components (e.g., the X-shredder, especially in Designs 3 and 4), such mutations may be actively selected for. To investigate we allowed homing-associated loss-of-function mutations to occur in each component of each construct. Mutation rates of 10e-3 have a small but significant impact on the performance of the three designs with an X-shredder, due to the accumulation of mutant constructs missing that component, while mutation rates of 10e-4 have negligible impact for all designs (Figs 5A, S6 and S7).

Fig. 5. The impact of evolutionary stability and added fitness costs on the performance of each of the 5 designs for population suppression.
The impact of evolutionary stability and added fitness costs on the performance of each of the 5 designs for population suppression.
(A) The effect of loss-of-function mutations on the minimum population size (number of adult females) reached. Solid lines are for the baseline case of no mutations, and dashed lines are with each component of each construct having a mutation rate of 10e-3 per homing event. Note that the effect is only visible for designs with an X-shredder, and if the mutation rate was 10e-4, the results for all designs would be virtually indistinguishable from the solid lines. (B, C) Contour plots showing combinations of fitness costs and pre-existing frequency of resistance giving a minimum population size of 10e-4 for different double drive designs. For (B) the costs are reductions in female fitness due to somatic expression of the nuclease targeting the A locus, and for (C) the costs are reductions in male fitness due to the X-shredder. Vertical lines indicate the cost is irrelevant, either because heterozygous females in any case have fitness 0 (Design 5 in (B)), or because the designs do not include an X-shredder (Designs 1 and 2 in (C)).

Second, we have assumed thus far that the genetic constructs have little unintended impact on survival or reproduction. Experiments with An. gambiae have revealed at least two unintended fitness costs can occur, a reduced fitness of homing heterozygous females due to somatic expression of the nuclease [25,32], and reduced fitness of males expressing an X-shredder, possibly due to paternal deposition of the nuclease and/or reduced sperm production [36]. The first of these costs is not relevant to Design 5 (because heterozygous females die anyway), and the second is not relevant to Designs 1 and 2 (because they do not use an X-shredder), but in other contexts, as expected, these costs reduce performance, requiring a lower frequency of pre-existing resistance in order to achieve a particular level of suppression (Fig 5B and 5C).

Population replacement

Gene drive can be used not only for population suppression but also to introduce a new desirable ‘cargo’ gene into a target population for population replacement or modification–for example, a gene reducing a mosquito’s ability to transmit a pathogen [37,38]. In double drive designs for population replacement the α construct would carry the cargo and homing by α would require β, while that by β could be either autonomous or depend on α (analogous to Designs 1 and 2 for population suppression; Fig 6A). Both α and β could be inserted into neutral sites, or into essential genes in such a way as to minimise fitness effects. We have modelled these approaches assuming, for purposes of illustration, the cargo imposes a dominant 20% fitness cost on females, and find that, again, such double drives can spread rapidly through target populations even when there is significant pre-existing resistance, and would not spread in non-target populations fixed for the resistant allele (Fig 6B and 6C). Unless there is virtually no pre-existing resistance at the differentiated locus, double drives can keep the frequency of the cargo gene above 95% much longer than a single drive construct targeting a differentiated locus, either neutral or essential (Fig 6D). In the single locus case selection rapidly increases the frequency of a pre-existing functional resistant allele because there is both significant variation and significant fitness differences (arising from the cost of the cargo gene) at the same locus. By contrast, in the double drive case there is one locus at which there are fitness difference (due to presence/absence of the cargo) but very little variation (initially none, and arising only after release due to rare end-joining and loss-of-cargo events), and another (differentiated) locus at which there is pre-existing variation but much smaller fitness differences (arising only due to the statistical correlation between alleles at the two loci). Finally, as with double drives for population suppression, the protection provided by a double drive for replacement can be extended even further if the two loci are tightly linked and α is inserted in an essential gene, with end-joining repair producing nonfunctional alleles (Fig 6D). In the latter case the only source of functional cleavage-resistant cargo-less alleles are the constructs that lose their cargo during homing, and, since this is assumed to occur at a much lower rate than end-joining repair, more generations are needed for such alleles to become common and the duration of protection is extended.

Fig. 6. Double drives for population replacement.
Double drives for population replacement.
(A) Alternative double drive designs. (B) Timecourse of allele frequencies for Design 2r in a target population assuming 50% pre-existing resistance; “cargoless α construct” refers to α constructs that have lost their cargo gene. The dynamics for Design 1r are qualitatively similar. (C) Allele frequencies for Designs 1r and 2r in a non-target population with 100% pre-existing resistance, assuming insertion of both α and β into neutral unlinked loci (solid lines); or α as a neutral insertion into a haplo-insufficient essential gene closely linked to β (r = 0.01) (dashed lines). (D) Duration of at least 95% of adult females carrying a cargo gene (whether the rest of the construct is functional or not) as a function of the pre-existing frequency of resistance at the B locus for double drives 1r and 2r, where solid and dashed lines are as in (C). Also shown for comparison are results for a single drive (S) carrying the cargo at a neutral A locus (solid line) or a haplo-insufficient essential gene assuming end-joining repair produces only nonfunctional resistant alleles (dashed lines). Note that results for insertion of β into a haplo-insufficient essential gene would be virtually indistinguishable from the solid lines (C, D). All plots assume 20% fitness cost of the cargo on females and a homing-associated loss-of-function mutation rate of 10e-3.

PAM site analysis in An. gambiae

To explore whether the type of population differentiation assumed in our modelling can exist in nature, we analysed published genome sequence data on An. gambiae mosquitoes from the Ag1000G project [39]. The Ag1000G dataset includes sequences from 16 mainland African populations and from populations on Mayotte and Bioko, two islands 500km off the east and 30km off the west coast of Africa, respectively. Note that in presenting this analysis we are not advocating the use of double drives on these islands, and merely wish to investigate whether the requisite differentiation can be found on island populations. For our analysis we focussed on potential PAM sequences (NGG or CCN), on the logic that a construct would be unlikely to mutate to recognise a new PAM, whereas this could occur for a protospacer. The entire dataset includes 57 million polymorphic sites, which we screened for PAM sites present in the island population and at a frequency <10%, <5%, or absent from all other populations. In Mayotte, for PAM sequences that were completely private to the island (i.e., not found in any other population), only 1 of them had no pre-existing resistance (i.e., was found in all 48 sequences from the island), whereas 25 had pre-existing resistance less than 20%, and 353 had pre-existing resistance less than 50%. PAM sequences with small but nonzero frequencies on the mainland were even more abundant (Fig 7). Bioko island is not as differentiated as Mayotte from the mainland populations, and the sample size is smaller (18 sequences), but still there are some potential candidate sites.

Fig. 7. Frequency of PAM sites in island populations of An.
Frequency of PAM sites in island populations of An.
gambiae. Numbers of PAM sites (NGG or CCN) with varying frequencies of resistance in samples of An. gambiae from two oceanic islands (blue map circles): (A) Bioko island (n = 18 sequences, from 9 individuals), and (B) Mayotte island (n = 48 sequences, from 24 individuals), where the PAM site frequency in each non-target population (black map circles) is <10%, <5% or 0% (i.e., target site resistance is >90%, >95% or 100%). GG or CC dinucleotides which varied by at least one base were considered to be resistant.

Discussion

Given that some of the most promising gene drive approaches for population control use (CRISPR-based) sequence-specific nucleases, an obvious way to limit their spread and impact is to exploit sequence differences between target and non-target populations. In this paper we have proposed using a double drive design, here defined as one that uses two constructs, inserted at different locations in the genome, both of which can increase in frequency, at least initially, and which interact such that the transmission of at least one of them depends on the other. Previously published examples that fit this definition include those for 2-locus under-dominance [16,19,40,41], and Medusa [42], tethered [17], integral [30], and transcomplementing [43] gene drives. As with single-construct gene drives, these various proposed designs differ in purpose (suppression vs. modification), release rate needed to initiate spread (low vs high threshold), and the molecular basis for the superMendelian inheritance (homing, toxin-antidote interactions, or a combination of the two), and the suggested rationales for adopting these designs over single drives include allowing more localised population control and a more modular product development pipeline. The requirement that both constructs can increase in frequency over time excludes split drives [14,4446] and killer-rescue systems [13,47], in which only one of the two components increases in frequency. In our proposed designs there is a division of labour between the two constructs, with one responsible for the desired impact (suppression or replacement) and the other for the population restriction, such that together they act as a double drive in the target population and as a split drive in non-target populations.

Note that if there are multiple populations of the same species requiring control, each with a different private allele, the same α construct could be used in each case, with only a change in the insertion site of the β construct and the corresponding gRNA. This flexibility may be particularly useful when the α construct requires significant optimisation [30]. Moreover, the same strategies may also be useful for controlled suppression of a target population even when there is no concern about non-target populations (and therefore no need to target a private allele): by appropriate choice of construct components and insertion sites a form of “planned obsolescence” could be achieved, with a wider and more predictable range of suppression profiles (e.g., depth and duration of suppression) possible than with conventional single drives [see also [48]].

For designs that involve interacting insertions at two or more loci, their population genetic dynamics and impact will usually depend on the statistical correlation between the constructs, and therefore also on the degree of linkage between them [12]. As previously demonstrated for split cleave-and-rescue designs by Oberhofer et al. [15], the degree of linkage between constructs can therefore be used as a tunable parameter to control the dynamics. They found that the expected impacts of a release were stronger and longer-lasting with closely linked constructs than with distantly linked ones, and we found much the same with our double drive designs, though there are some differences in detail between the systems. In particular, with split cleave-and-rescue designs, which do not rely on homing, if there is complete linkage (i.e., no meiotic recombination) then the system behaves the same as a single locus construct, whereas that is not the case for our homing-based double drives, where constructs can be separated if one homes and the other does not, even if there is never any meiotic crossing over between them. Thus, in our model, setting r = 0 (while allowing separate homing events) does not reproduce the single drive dynamics. If the insertions are physically very close to each other, then there may be some direct mechanistic interaction between them (e.g., binding of one nuclease complex preventing the other from binding, or resection of DNA during the repair process leading to co-homing of the two insertions), or simultaneous cleavage may lead to a large deletion. We have assumed that our constructs are far enough apart as not to interact in this way (e.g., r = 0.01 corresponds on average to between 600kb and 1Mb in An. gambiae [49]).

We have considered a range of double drive designs of increasing robustness, as judged by their ability to cope with an increasing frequency of pre-existing resistance at the differentiated locus. The simplest designs do not have any component beyond those needed for any CRISPR-based construct, and so should be widely applicable [43]. More powerful constructs can be made by adding an X-shredding sex ratio distorter to the load-inducing construct; these have been most effectively demonstrated in An. gambiae mosquitoes [11,50], but may also work more broadly [51]. Note that the optimal timing of homing and X-shredding during gametogenesis may be different, requiring different and compatible control sequences, which will need to be taken into account in construct design [50,52]. The ability to control gRNA expression in a tissue-specific manner would be helpful in this regard. In other species there are other ways to distort the sex ratio [5355], and it would be interesting to model whether these alternatives would be expected to have the same impact as an X-shredder in the context of a double drive. Potentially an even simpler way to increase the load may be to include additional gRNAs in the α construct that cleave and knock out the function of other female fertility genes elsewhere in the genome [31,56]. The effect of such gRNAs would depend on the heterozygous and homozygous fitness effects of the mutations caused and, again, on the degree of linkage with the target site, and further modelling would be needed to investigate whether such a strategy is worthwhile. The most powerful design we considered targets a female-specific haplo-insufficient gene, or otherwise causes dominant female sterility or lethality. Such genes are not common, but there are some possible candidates [5760], and our modelling motivates the search for others. Finally, performance (in terms of being able to cope with ever higher frequencies of pre-existing resistance) could presumably also be improved by using a third construct, to construct a triple drive, though modelling would be required to explore the implications of the many different configurations this extension would allow.

The proposed strategy requires that there be a differentiated locus between target and non-target populations. It need not be an essential gene, and could even be selectively neutral. Our focus has been on using so-called private alleles–sequences that are present (but not necessarily fixed) in the target population, and absent (or of negligible frequency) in non-target populations. Our analysis of PAM sites in An. gambiae indicates that appropriately differentiated sites may exist in island populations of this species, though our analysis must be considered preliminary: the dataset does not include mainland sites in closest proximity to the island populations, where differentiation may be lower, and we have not considered potential polymorphism in the protospacer sequence (which, if present, may require the use of multiple gRNAs). We have focussed on nucleotide variation at PAM sites on the assumption that a construct is unlikely to mutate to recognise a new PAM; structural variation in the protospacer region may also be an appropriate basis for geographically restricting double drive spread. We have also not attempted to determine whether the observed differentiation is due solely to mutation and drift, or if selection may be involved as well.

Note that the single drives modelled by Sudweeks et al. [23] require the opposite type of differentiation: sequences that are fixed in the target population, even if not private (i.e., even if found at appreciable frequencies in the non-target population [61]). In this latter scenario the challenge is not so much to have an impact on the target population as to not have an impact on the non-target population. What constitutes “acceptable non-impact” may differ widely from one use case to another and must be assessed on a case-by-case basis: in some circumstances spread of the construct and a transient decline in population size followed by recovery may be acceptable, whereas in others any significant spread of the construct may be unacceptable, regardless of impact on population size. Designs with non-autonomous homing of the β construct (Designs 2, 4, and 5) should be less likely to increase in frequency in the non-target population, and may therefore be preferable. We have focused in this paper on differentiated loci on autosomes, but note that for Design 5 the X-shredder is required for the spread of the α construct and, in principle, one could achieve population-restricted spread if the shredder targeted a population-specific sequence on the X chromosome (rather than inserting it into a population-specific autosomal sequence). In many species the X chromosome shows greater population differentiation than autosomes [62], so this alternative may be useful. Finally, if there are no private alleles in the target population, it may be worthwhile considering a two-step approach of first introducing a private allele into a population and then using that allele to control the population [22]. The ability of double drives to exploit private alleles that are selectively neutral and that have a frequency of only 50% (suppression) or 20% (modification) potentially makes this approach more feasible than would otherwise be the case.

In this paper we have used a simple high-level modelling framework in which the generations are discrete, the population is well mixed, and dynamics are deterministic. This framework is appropriate for strategic models aiming to identify candidate approaches that are worthy of further investigation. For any specific use case the appropriate tactical models would need to be developed that incorporate more biological detail, including spatial and stochastic effects. In spatially distributed populations with local mating the statistical association between alleles at the two loci may differ from that in our well mixed model, and the quantitative dynamics thereby affected. Issues of evolutionary stability and the breakdown of constructs can also be more important in such models, as previously demonstrated for single drive homing constructs for population replacement [63]. Such extensions will be particularly important when the goal is to eliminate the target population, which is not possible in our deterministic models. Instead, we have reported the minimum relative population size achieved, which is expected to be related to the size of a population that could be eliminated, but determining the precise connection will require bespoke modelling tailored to a specific situation. Further extensions would be needed to allow for on-going movement between target and non-target populations–if there is on-going immigration into the target population, and this cannot be stopped, then it may not be possible to eliminate the target population with a single release of a double drive. Nevertheless, such a release may be sufficient to suppress the population to such an extent that it can be controlled by other means, including recurrent releases of the same constructs. If one is able to achieve an initial release rate of 1% into a target population, and that suppresses the population by a factor of 1000, then the same releases going forwards will constitute a 10-fold inundation, and self-limiting genetic approaches may be sufficient.

Methods

The basic deterministic modelling structure follows that of Burt & Deredec [12]. In brief, populations have discrete generations, mating is random, there are two life stages (juveniles and adults), and juvenile survival is density dependent according to the Beverton-Holt model with an intrinsic rate of increase (Rm) equal to 6 [56]. Genetic parameter values (rates of DNA cleavage, rates of alternative repair pathways, and the sex ratio produced by X-shredding) are as estimated from An. gambiae (S1 Table) [11,25,52]. Constructs may be inserted into a haplo-sufficient or haplo-insufficient female-essential gene (in which case gene function is disrupted), a selectively neutral sequence (in which case the insertion is also selectively neutral), or a haplo-sufficient or haplo-insufficient gene required for male and female viability (in which case the insertion is again selectively neutral, because it contains a re-coded version of the target gene [14,26,29], or is inserted in an artificial intron [30]). For constructs inserted into an essential gene we assume end-joining repair produces non-functional cleavage-resistant alleles [25,64], while for constructs inserted into selectively neutral sites the products of end-joining repair are also neutral. In all models we assume individuals with an intact CRISPR system suffer a 1% fitness cost for every different gRNA they carry as a cost of off-target cleavage, and for population replacement we assume the cargo gene imposes a 20% fitness cost on females. Both these costs are assumed to be dominant (i.e., not dosage-dependent). For simplicity, we assume all fitness costs affect survival after density dependent juvenile mortality and before censusing (e.g., as if pupae die). All results are for populations censused at the adult stage. Releases are of heterozygous adult males at 0.1% of the pre-release number of males, and if the two constructs are linked then they are assumed to be in cis; for constructs released in separate males we assume release rates of 0.05% of each. Additional details and a list of parameters and their baseline values is given in S1 Text, S1 and S2 Tables. Code for implementation of the simulations is available on GitHub (https://github.com/KatieWillis/DoubleDriveSimulator).

For the PAM site analysis we screened the Ag1000G phase II SNP data for PAM sites (GG or CC dinucleotides) showing variation between samples at one or both nucleotides. PAM site frequencies were calculated per sampling location and filtered for those present in the island population and at <10%, 5%, or absent from all other populations, excluding those containing >5% missing data in at least one sampled population. Further details are given in the S1 Text. The map in Fig 7 was produced using the cartopy python package [65].

Supporting information

S1 Text [pdf]
Supplemental methods.

S1 Table [pdf]
Model parameters and baseline values.

S2 Table [pdf]
Host gene disruption fitness costs.

S1 Fig [a]
Timecourse for the relative number of females over time for Designs 1 and 2.

S2 Fig [a]
Effect of releasing constructs in the same or different males.

S3 Fig [tif]
Example time courses for double drives for population suppression.

S4 Fig [a]
The impact of alternative designs and their variants on population suppression as a function of the pre-existing frequency of resistance.

S5 Fig [tif]
Timecourse of allele frequencies and population suppression (1-relative number of females) for Designs 1–5 in non-target populations where the resistant allele is present at 100%.

S6 Fig [tif]
Timecourses for Designs 3, 4 and 5 assuming pre-existing frequency of resistance of 1%, where homing-associated loss-of-function mutations occur for each component of each construct with probability 10e-3.

S7 Fig [tif]
Loss-of-function mutation rates of 10e-4 have minimal impact on the extent of population suppression.


Zdroje

1. Burt A, Trivers R. Genes in conflict: the biology of selfish genetic elements. Cambridge: Belknap Press of Harvard University Press; 2006.

2. Burt A, Crisanti A. Gene drive: evolved and synthetic. ACS Chem Biol. 2018;13(2):343–6. doi: 10.1021/acschembio.7b01031 WOS:000426012800008. 29400944

3. Burt A, Coulibaly M, Crisanti A, Diabate A, Kayondo JK. Gene drive to reduce malaria transmission in sub-Saharan Africa. J Resp Innov. 2018;5:S66–S80. doi: 10.1080/23299460.2017.1419410 WOS:000434459000005.

4. Raban RR, Marshall JM, Akbari OS. Progress towards engineering gene drives for population control. J Exp Biol. 2020;223. ARTN jeb208181 doi: 10.1242/jeb.208181 WOS:000541774100005. 32034041

5. Teem JL, Alphey L, Descamps S, Edgington MP, Edwards O, Gemmell N, et al. Genetic biocontrol for invasive species. Front Bioeng Biotech. 2020;8. ARTN 452 doi: 10.3389/fbioe.2020.00452 WOS:000540500500001. 32523938

6. Hay BA, Oberhofer G, Guo M. Engineering the composition and fate of wild populations with gene drive. Ann Rev Entomol. 2021;66:407–34. Epub 2020/10/10. doi: 10.1146/annurev-ento-020117-043154 33035437

7. 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 Pop Biol. 2016;108:51–69. doi: 10.1016/j.tpb.2015.11.005 WOS:000372560000005. 26704073

8. North AR, Burt A, Godfray HCJ. Modelling the potential of genetic control of malaria mosquitoes at national scale. BMC Biol. 2019;17. ARTN 26 doi: 10.1186/s12915-019-0645-5 WOS:000463655000001. 30922310

9. North AR, Burt A, Godfray HCJ. Modelling the suppression of a malaria vector using a CRISPR-Cas9 gene drive to reduce female fertility. BMC Biol. 2020;18(1). ARTN 98 doi: 10.1186/s12915-020-00834-z WOS:000561855300001. 32782000

10. Thomas DD, Donnelly CA, Wood RJ, Alphey LS. Insect population control using a dominant, repressible, lethal genetic system. Science. 2000;287(5462):2474–6. doi: 10.1126/science.287.5462.2474 WOS:000086202200041. 10741964

11. 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. doi: 10.1038/ncomms4977 WOS:000338834900001. 24915045

12. Burt A, Deredec A. Self-limiting population genetic control with sex-linked genome editors. Proc Roy Soc Lond B. 2018;285(1883). ARTN 20180776 doi: 10.1098/rspb.2018.0776 WOS:000439907900011. 30051868

13. Gould F, Huang YX, Legros M, Lloyd AL. A killer-rescue system for self-limiting gene drive of anti-pathogen constructs. Proc Roy Soc Lond B. 2008;275(1653):2823–9. doi: 10.1098/rspb.2008.0846 WOS:000260611200006. 18765342

14. Terradas G, Buchman AB, Bennett JB, Shriner I, Marshall JM, Akbari OS, et al. Inherently confinable split-drive systems in Drosophila. Biorxiv. 2020.

15. Oberhofer G, Ivy T, Hay BA. Split versions of Cleave and Rescue selfish genetic elements for measured self limiting gene drive. PLOS Genet. 2021;17(2):e1009385. Epub 2021/02/19. doi: 10.1371/journal.pgen.1009385 33600432.

16. Davis S, Bax N, Grewe P. Engineered underdominance allows efficient and economical introgression of traits into pest populations. J Theor Biol. 2001;212(1):83–98. doi: 10.1006/jtbi.2001.2357 WOS:000170918500007. 11527447

17. Dhole S, Lloyd AL, Gould F. Tethered homing gene drives: a new design for spatially restricted population replacement and suppression. Evol Appl. 2019;12(8):1688–702. doi: 10.1111/eva.12827 WOS:000484471100014. 31462923

18. Edgington MP, Harvey-Samuel T, Alphey L. Split drive killer-rescue provides a novel threshold-dependent gene drive. Sci Reports. 2020;10(1). ARTN 20520 doi: 10.1038/s41598-020-77544-7 WOS:000596296000007. 33239631

19. Champer J, Champer SE, Kim IK, Clark AG, Messer PW. Design and analysis of CRISPR-based underdominance toxin-antidote gene drives. Evol Appl. 2020. doi: 10.1111/eva.13180 WOS:000600504300001.

20. Dhole S, Vella MR, Lloyd AL, Gould F. Invasion and migration of spatially self-limiting gene drives: A comparative analysis. Evol Appl. 2018;11(5):794–808. doi: 10.1111/eva.12583 WOS:000433572400018. 29875820

21. Sanchez CHM, Bennett JB, Wu SL, Rasic G, Akbari OS, Marshall JM. Modeling confinement and reversibility of threshold-dependent gene drive systems in spatially-explicit Aedes aegypti populations. BMC Biology. 2020;18(1). ARTN 50 doi: 10.1186/s12915-020-0759-9 WOS:000535887600001. 32398005

22. Esvelt KM, Smidler AL, Catteruccia F, Church GM. Concerning RNA-guided gene drives for the alteration of wild populations. Elife. 2014;3. doi: 10.7554/eLife.03401 WOS:000209690800001. 25035423

23. Sudweeks J, Hollingsworth B, Blondel DV, Campbell KJ, Dhole S, Eisemann JD, et al. Locally fixed alleles: a method to localize gene drive to island populations. Sci Reports. 2019;9. ARTN 15821 doi: 10.1038/s41598-019-51994-0 WOS:000493716000022. 31676762

24. Burt A. Site-specific selfish genes as tools for the control and genetic engineering of natural populations. Proc Roy Soc Lond B. 2003;270(1518):921–8. doi: 10.1098/rspb.2002.2319 WOS:000182969000005. 12803906

25. Kyrou K, Hammond AM, Galizi R, Kranjc N, Burt A, Beaghton AK, et al. A CRISPR-Cas9 gene drive targeting doublesex causes complete population suppression in caged Anopheles gambiae mosquitoes. Nat Biotechnol. 2018;36(11):1062–+. doi: 10.1038/nbt.4245 WOS:000450374000019. 30247490

26. Oberhofer G, Ivy T, Hay BA. Cleave and Rescue, a novel selfish genetic element and general strategy for gene drive. Proc Natl Acad Sci USA. 2019;116(13):6250–9. doi: 10.1073/pnas.1816928116 WOS:000462382800067. 30760597

27. Champer J, Yang E, Lee E, Liu JX, Clark AG, Messer PW. A CRISPR homing gene drive targeting a haplolethal gene removes resistance alleles and successfully spreads through a cage population. Proc Natl Acad Sci USA. 2020;117(39):24377–83. doi: 10.1073/pnas.2004373117 WOS:000576672700020. 32929034

28. Courchamp F, Berek L, Gascoigne J. Allee effects in ecology and conservation. Oxford: Oxford University Press; 2008.

29. Champer J, Lee E, Yang E, Liu C, Clark AG, Messer PW. A toxin-antidote CRISPR gene drive system for regional population modification. Nat Commun. 2020;11(1). ARTN 1082 doi: 10.1038/s41467-020-14960-3 WOS:000518590600001. 32109227

30. Nash A, Urdaneta GM, Beaghton AK, Hoermann A, Papathanos PA, Christophides GK, et al. Integral gene drives for population replacement. Biol Open. 2019;8(1). ARTN bio037762 doi: 10.1242/bio.037762 WOS:000457406100012. 30498016

31. Champer J, Kim IK, Champer SE, Clark AG, Messer PW. Performance analysis of novel toxin-antidote CRISPR gene drive systems. BMC Biol. 2020;18(1):27. Epub 2020/03/14. doi: 10.1186/s12915-020-0761-2 32164660; PubMed Central PMCID: PMC7068947.

32. Simoni A, Hammond AM, Beaghton AK, Galizi R, Taxiarchi C, Kyrou K, et al. A male-biased sex-distorter gene drive for the human malaria vector Anopheles gambiae. Nat Biotechnol. 2020;38(9):1097–. doi: 10.1038/s41587-020-0658-1 WOS:000556636100002. 32764730

33. Hicks WM, Kim M, Haber JE. Increased mutagenesis and unique mutation signature associated with mitotic gene conversion. Science. 2010;329(5987):82–5. doi: 10.1126/science.1191125 WOS:000279402700037. 20595613

34. Simoni A, Siniscalchi C, Chan Y-S, Huen DS, Russell S, Windbichler N, et al. Development of synthetic selfish elements based on modular nucleases in Drosophila melanogaster. Nucl Acids Res. 2014;42(11):7461–72. doi: 10.1093/nar/gku387 WOS:000338769400062. 24803674

35. Rodgers K, McVey M. Error-prone repair of DNA double-strand breaks. J Cell Physiol. 2016;231(1):15–24. doi: 10.1002/jcp.25053 WOS:000362217800004. 26033759

36. Pollegioni P, North AR, Persampieri T, Bucci A, Minuz RL, Groneberg DA, et al. Detecting the population dynamics of an autosomal sex ratio distorter transgene in malaria vector mosquitoes. J Appl Ecol. 2020;57(10):2086–96. doi: 10.1111/1365-2664.13702 WOS:000550495400001. 33149368

37. Marshall JM, Akbari OS. Gene drive strategies for population replacement. In: Adelman ZN, editor. Genetic control of malaria and dengue. London: Academic Press; 2016. p. 169–200.

38. Adelman ZN, Basu S, Myles KM. Engineering pathogen resistance in mosquitoes. In: Adelman ZN, editor. Genetic control of malaria and dengue. London: Academic Press; 2016. p. 277–304.

39. Anopheles gambiae Genomes Consortium. Genome variation and population structure among 1142 mosquitoes of the African malaria vector species Anopheles gambiae and Anopheles coluzzii. Genome Res. 2020;30(10):1533–46. Epub 2020/09/30. doi: 10.1101/gr.262790.120 32989001; PubMed Central PMCID: PMC7605271.

40. Magori K, Gould F. Genetically engineered underdominance for manipulation of pest populations: a deterministic model. Genetics. 2006;172(4):2613–20. doi: 10.1534/genetics.105.051789 WOS:000237225800048. 16415364

41. Khamis D, El Mouden C, Kura K, Bonsall MB. Ecological effects on underdominance threshold drives for vector control. J Theor Biol. 2018;456:1–15. doi: 10.1016/j.jtbi.2018.07.024 WOS:000445321400001. 30040965

42. Marshall JM, Hay BA. Medusa: a novel gene drive system for confined suppression of insect populations. PLOS One. 2014;9(7). ARTN e102694 doi: 10.1371/journal.pone.0102694 WOS:000339614100049. 25054803

43. Del Amo VL, Bishop AL, Sanchez CHM, Bennett JB, Feng XC, Marshall JM, et al. A transcomplementing gene drive provides a flexible platform for laboratory investigation and potential field deployment. Nat Commun. 2020;11(1). ARTN 352 doi: 10.1038/s41467-019-13977-7 WOS:000514637800004. 31953404

44. DiCarlo JE, Chavez A, Dietz SL, Esvelt KM, Church GM. Safeguarding CRISPR-Cas9 gene drives in yeast. Nat Biotechnol. 2015;33(12):1250–+. doi: 10.1038/nbt.3412 WOS:000366387700016. 26571100

45. Champer J, Chung J, Lee YL, Liu C, Yang E, Wen ZX, et al. Molecular safeguarding of CRISPR gene drive experiments. Elife. 2019;8. ARTN e41439 doi: 10.7554/eLife.41439 WOS:000457468900001. 30666960

46. Li M, Yang T, Kandul NP, Bui M, Gamez S, Raban R, et al. Development of a confinable gene drive system in the human disease vector Aedes aegypti. Elife. 2020;9. ARTN e51701 doi: 10.7554/eLife.51701 WOS:000508594700001. 31960794

47. Webster SH, Vella MR, Scott MJ. Development and testing of a novel killer-rescue self-limiting gene drive system in Drosophila melanogaster. Proc Roy Soc Lond B. 2020;287(1925). ARTN 20192994 doi: 10.1098/rspb.2019.2994 WOS:000526900400001. 32292114

48. Zapletal J, Najmitabrizi N, Erraguntla M, Lawley MA, Myles KM, Adelman ZN. Making gene drive biodegradable. Phil Trans Roy Soc Lond B. 2021;376(1818). ARTN 20190804 doi: 10.1098/rstb.2019.0804 WOS:000603646400003. 33357058

49. Pombi M, Stump AD, Della Torre A, Besansky NJ. Variation in recombination rate across the X chromosome of Anopheles gambiae. Am J Trop Med Hyg. 2006;75(5):901–3. doi: 10.4269/ajtmh.2006.75.901 WOS:000242189100022. 17123984

50. 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 Reports. 2016;6. doi: 10.1038/srep31139 WOS:000380652300001. 27484623

51. Fasulo B, Meccariello A, Morgan M, Borufka C, Papathanos PA, Windbichler N. A fly model establishes distinct mechanisms for synthetic CRISPR/Cas9 sex distorters. PLOS Genet. 2020;16(3). ARTN e1008647 doi: 10.1371/journal.pgen.1008647 WOS:000524758200018. 32168334

52. Hammond A, Karlsson X, Morianou I, Kyrou K, Beaghton A, Gribble M, et al. Regulating the expression of gene drives is key to increasing their invasive potential and the mitigation of resistance. PLOS Genet. 2021;17(1):e1009321. doi: 10.1371/journal.pgen.1009321 33513149

53. Pane A, Salvemini M, Bovi PD, Polito C, Saccone G. The transformer gene in Ceratitis capitata provides a genetic basis for selecting and remembering the sexual fate. Development. 2002;129(15):3715–25. WOS:000177570800017. 12117820

54. Meccariello A, Salvemini M, Primo P, Hall B, Koskiniot P, Dalikova M, et al. Maleness-on-the-Y (MoY) orchestrates male sex determination in major agricultural fruit fly pests. Science. 2019;365(6460):1457–+. doi: 10.1126/science.aax1318 WOS:000488838600049. 31467189

55. Aryan A, Anderson MAE, Biedler JK, Qi YM, Overcash JM, Naumenko AN, et al. Nix alone is sufficient to convert female Aedes aegypti into fertile males and myo-sex is needed for male flight. Proc Natl Acad Sci USA. 2020;117(30):17702–9. doi: 10.1073/pnas.2001132117 WOS:000555848400016. 32661163

56. Deredec A, Godfray HCJ, Burt A. Requirements for effective malaria control with homing endonuclease genes. Proc Natl Acad Sci USA. 2011;108(43):E874–E80. doi: 10.1073/pnas.1110717108 WOS:000296378100004. 21976487

57. Hatanaka K, Okada M. Retarded nuclear migration in Drosophila embryos with aberrant F-actin reorganization caused by maternal mutations and by cytochalasin treatment. Development. 1991;111(4):909–20. WOS:A1991FK57600007. 1879360

58. O’Leary S, Adelman ZN. Disrupting female flight in the vector Aedes aegypti. Biorxiv. 2019.

59. Basrur NS, De Obaldia ME, Morita T, Herre M, von Heynitz RK, Tsitohay YN, et al. fruitless mutant male mosquitoes gain attraction to human odor. Elife. 2020;9. ARTN e63982 doi: 10.7554/eLife.63982 WOS:000610890200001. 33284111

60. Krzywinska E, Ferretti L, Li J, Li JC, Chen CH, Krzywinski J. femaleless controls sex determination and dosage compensation pathways in females of Anopheles mosquitoes. Curr Biol. 2021. Epub 2021/01/09. doi: 10.1016/j.cub.2020.12.014 33417880.

61. Oh KP, Shiels AB, Shiels L, Blondel DV, Campbell KJ, Saah JR, et al. Population genomics of invasive rodents on islands: genetic consequences of colonization and prospects for localized synthetic gene drive. Evol Appl. 2021. https://doi.org/10.1111/eva.13210.

62. Presgraves DC. Evaluating genomic signatures of "the large X-effect" during complex speciation. Mol Ecol. 2018;27(19):3822–30. doi: 10.1111/mec.14777 WOS:000446838400008. 29940087

63. Beaghton A, Hammond A, Nolan T, Crisanti A, Godfray HCJ, Burt A. Requirements for driving antipathogen effector genes into populations of disease vectors by homing. Genetics. 2017;205(4):1587–96. doi: 10.1534/genetics.116.197632 WOS:000401126600017. 28159753

64. Beaghton AK, Hammond A, Nolan T, Crisanti A, Burt A. Gene drive for population genetic control: non-functional resistance and parental effects. Proc Roy Soc Lond B. 2019;286(1914). ARTN 20191586 doi: 10.1098/rspb.2019.1586 WOS:000504858100006. 31662083

65. Met Office. Cartopy: a cartographic python library with a matplotlib interface. 2010–2015. http://scitools.org.uk/cartopy.


Článek vyšel v časopise

PLOS Genetics


2021 Číslo 3
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

Důležitost adherence při depresivním onemocnění
nový kurz
Autoři: MUDr. Eliška Bartečková, Ph.D.

Koncepce osteologické péče pro gynekology a praktické lékaře
Autoři: MUDr. František Šenk

Sekvenční léčba schizofrenie
Autoři: MUDr. Jana Hořínková, Ph.D.

Hypertenze a hypercholesterolémie – synergický efekt léčby
Autoři: prof. MUDr. Hana Rosolová, DrSc.

Multidisciplinární zkušenosti u pacientů s diabetem
Autoři: Prof. MUDr. Martin Haluzík, DrSc., prof. MUDr. Vojtěch Melenovský, CSc., prof. MUDr. Vladimír Tesař, DrSc.

Všechny kurzy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#