Dynamics of a combined medea-underdominant population transformation system
BMC Evolutionary Biology volume 14, Article number: 98 (2014)
Transgenic constructs intended to be stably established at high frequencies in wild populations have been demonstrated to “drive” from low frequencies in experimental insect populations. Linking such population transformation constructs to genes which render them unable to transmit pathogens could eventually be used to stop the spread of vector-borne diseases like malaria and dengue.
Generally, population transformation constructs with only a single transgenic drive mechanism have been envisioned. Using a theoretical modelling approach we describe the predicted properties of a construct combining autosomal Medea and underdominant population transformation systems. We show that when combined they can exhibit synergistic properties which in broad circumstances surpass those of the single systems.
With combined systems, intentional population transformation and its reversal can be achieved readily. Combined constructs also enhance the capacity to geographically restrict transgenic constructs to targeted populations. It is anticipated that these properties are likely to be of particular value in attracting regulatory approval and public acceptance of this novel technology.
Curbing the spread of vector borne diseases such as malaria or dengue is possible by eliminating the transmission capabilities of the insect vectors. One of the many approaches to achieve this is population transformation of vector species. In the most commonly discussed application of population transformation the aim is to introduce transgenes into insect populations which render them refractory to spreading diseases. Usually the technique seeks to use evolutionary principles to establish such transgenes at high frequency in populations through the release of genetically transformed stocks (also called population replacement, ). Synthetic disease refractory genes have already been developed for human malaria, dengue fever and avian malaria [2–5]. However, to stably transform insect populations with transgenes that are not selectively advantageous it will be necessary to link refractory transgenes to systems that drive them to high frequency in a population [1, 6–8]. Three transgenic population transformation systems have been shown to be effective in laboratory populations of insects. One is a homing endonuclease based system (HEG), which works by converting heterozygotes to homozygotes . The remaining two systems work by reducing the average fitness of heterozygotes and are: Medea [10, 11] and a bi-allelic form of underdominance . Here we explore theoretically a mono-allelic form of underdominance the implementation of which has to date not been published.
While most studies examine the theoretical properties of transgenic constructs embodying single drive mechanisms [7–10, 13], the observation that “most of them have specific characteristics that make them less than ideal” led Huang et al. 2007  to explore combinations. They demonstrated that certain combinations resulted in enhanced properties relative to single systems while others had the opposite effect. Here we take an analogous approach for autosomal Medea and mono-allelic underdominance constructs (not examined in Huang et al. 2007 ). We provide a rigorous and flexible analytical framework to explore salient properties across the entire parameter space. Intuitively, the inclusion within a single transgenic construct of more than one drive mechanism provides a degree of resilience to either mutations in the transgenic construct or to drive-resistance alleles which may exist in target population. While the value of this desirable functional redundancy is not analytically explored here, it does however provides an additional motivation for analyzing the properties of combined systems. Similar to Huang et al. the motivation for the analysis presented here comes from the realization that intuitive predictions about combined systems can be misleading and that identifying the parameter space where synergistic enhancements occur can motivate technical developments, including the development of mono-allelic form of underdominance.
We briefly summarise the previously known properties of Medea and mono-allelic underdominant systems separately. Then we look in turn at each of the properties of interest and determine if the combined model performs better than each of the techniques independently. The discussion focuses on the impact of a combined system and provides an assessment of its strengths and weaknesses.
Natural Maternal effect dominant embryonic arrest (Medea) alleles were first discovered in Tribolium flour beetles  and have also been reported in the mouse [16, 17]. They derive their ability to invade populations by maternally induced lethality of wildtype offspring not inheriting a Medea allele (Figure 1) . Thus the wildtype homozygous offspring of the heterozygous mother die with a certain probability d. Despite the mechanism(s) by which natural Medea elements exert their maternal effect remaining unknown, Chen et al.  were able to generate a synthetic system (M e d e amyd88) which mimics their evolutionary properties.
To date, the only published Medea construct (M e d e amyd88) has been inserted on an autosomal chromosome in D. melanogaster[10, 11]. Autosomal Medea insertions unlike sex-chromosome insertions  exhibit a high-frequency stable equilibrium when the transgenic construct is associated with any fitness cost (see Figure 2a). As described previously [7, 8, 10, 13, 18] this stable equilibrium results in the persistence of wildtype alleles in populations transformed with autosomal Medea constructs (Figure 2a, though if a linked refractory gene is dominant this is likely to prove unproblematic from the perspective of target disease control).
When the heterozygote is less fit than both the possible homozygotes then we have a case of underdominance. However there are only a few examples where alleles at a given locus have been robustly inferred to exhibit underdominance . In a random mating, Hardy-Weinberg population, rarer alleles have larger sojourn times in the heterozygote state, consequently where an underdominant construct is rare it will mostly be in this unfit genotype. Due to the inherently unstable nature of underdominance, if the construct exceeds a threshold value through releases of sufficient homozygotes it is predicted to proceed to fixation within the population (Figure 2b). Intentional underdominant population transformation is inherently reversible where it is realistically possible to release sufficient wildtype individuals to traverse the unstable equilibrium in the lower frequency direction. However, underdominant constructs can be viewed as unappealing when transforming large populations due to the high release numbers required to initiate population transformation (Figure 2b) [20, 21]. The mono-allelic underdominance modeled here describes the situation where there is a transgenic allele at a single autosomal locus (the site of the transgenic construct integration). We have only examined situations where an insert is underdominant in both sexes. A recent publication  describing the development of a single locus bi-allelic form of underdominance where there are two functionally distinct transgenic alleles is not applicable to the mono-allelic underdominance analysis described here.
Medea and underdominance in a single transgenic construct
Here we explore the properties of combining both Medea and underdominance in a single transgenic construct on an autosome. As single locus transgenic underdominance effective in both sexes cannot by definition be configured on sex chromosomes we have modeled both systems on autosomes to permit the most direct comparison between single and combined systems. By combining systems, some properties will be discounted, remain the same or synergistically enhanced. We find a broad parameter space where the applied properties of single systems can be argued to be synergistically enhanced. The principle criteria being: (i) lower transformation threshold, (ii) faster population transformation and (iii) enhanced spatial stability of the transformed population.
Methods and results
Genotype fitnesses and expected dynamics
The recursion dynamics are analysed for genotype frequencies as maternal-effect killing violates the Hardy-Weinberg principle. With Medea the action of selection on wildtype homozygotes depends not only on their current state but also on the maternal genotype. Here we have the three genotypes, wildtype homozygous, transgenic homozygous and the heterozygous represented by ++, MM, and M+ respectively. We set the fitness of the wildtype homozygote, ++, to 1. The relative fitnesses of the MM homozygote and the M+ heterozygote are given by ν and ω. The parameter d measures the degree of lethality of homozygous wildtype offspring from Medea carrying mothers, from no Medea effect (d=0) to complete lethality (d=1). Using Table 1 we calculate the expected frequencies of all three genotypes in the next generation as,
where x, y, and z are the frequencies of MM, M+, and ++ respectively in the current generation and x′, y′, and z′ are the expected frequencies in the next generation (in  differences in fitness were ascribed to differences in maternal fecundity rather than zygotic genotypes as is done here). The total contribution from all genotypes in the population (i.e., the average fitness) is given by . It is the sum of the right hand sides of the set of Eqs. (1) . Another way to view the recursion equations is , where f x is the average fitness of the MM genotype . Equating the fitnesses of the three genotypes helps us to solve for the fixed points of this dynamical system (see Appendix). For d=1 there can be an unstable internal equilibrium (Appendix Eq. (A.3)). From the point of view of reversibility it is ideal to have this equilibrium as close as possible to one-half (see Figure 3). This is possible when the fitness values of the heterozygote and the Medea homozygote sum up to unity (see Appendix Eq. (A.5)) as can be seen in Figure 4. The fitnesses of the systems described here are assumed to relate only to the drive mechanism (i.e. without linked refractory genes). In Figure 2 we illustrate how selecting a Medea construct with appropriate parameters results in the combined system having an internal equilibrium closer to the ideal one-half. The release thresholds are determined by the unstable fixed points of the system. As illustrated in Figure 4, combined systems have the potential to be engineered towards an optimal unstable equilibrium value of 0.5 (Eqs. (A.4) and (A.5)). A release threshold substantially smaller than 0.5 would make the construct unappealing from the point of view of reversibility (Figure 3). However, if the size of the target population is large and the capacity to reverse population transformation is not important, then a ‘Medea only’ construct would be the most efficient approach. Solely from the perspective of initiating population transformation, ‘underdominance only’ is disadvantageous as it requires multigenerational release of a large numbers of individuals (albeit smaller than the numbers required for sterile male release programs ).
It is generally appreciated that once releases commence, population transformation should occur as rapidly as possible and proceed to complete fixation. This minimizes the possibility of selection for insects resistant to the transformed construct. Furthermore the pathogen itself could evolve mechanisms to evade the effects of the linked refractory genes. A rapid and complete fixation of the transgenic construct and elimination of the pathogen minimizes both possibilities (neither of which are explicitly modeled here). Clearly, releasing as many individuals as is feasible is an effective way to speed population transformation [21, 25]. We show that the time taken to achieve population transformation can also be reduced by combining two systems, where even very weak Medea (d≤0.2) has a large impact on the speed of transformation. (see Figure 5 and 0.65 starting frequency). The acceleration can also occur during reversal of population transformation. Knowledge of this effect will permit the design of efficient release strategies for both the initiation of population transformation and for its rapid reversal.
Population structure dynamics
We consider a simple two-deme model of population structure, where two populations of large and equal size are coupled by a symmetrical fraction of migrants m between the populations in each generation. Considering asymmetries in population sizes, migration needs to be dealt with separately, as in . Also migration dynamics with an explicitly set spatial system has been recently assessed  (albeit not for a combined system). In population i the expected genotype frequency of genotype k after migration is , where gk,i is the frequency of the kth genotype in population i and gk,j is the kth genotype frequency in population j. These adjusted genotype frequencies can then be substituted into Eqs. (1).
We initialize the two populations where the Medea allele is almost fixed in one and almost lost in the other. The recursions were performed for different migration rates, slowly incremented in units of 10−3. The migration rate where the difference in allele frequencies between the two populations fell below 1% (thus assuming populations to have reached an equilibrium), was recorded as the critical migration rate. At lower than critical migration rates the combined systems will not spread far from a successfully transformed zone, and will be resistant to loss by immigration. We evaluated the critical migration rate which allow the transformation of a local population stably (Figure 6). For a varying heterozygote fitness ω ranging from 0.01 to 0.95 we consistently see that having a Medea construct provides more geographical stability as compared to a system without Medea even as we move from a system with directional selection against Medea to an underdominant system. Figure 6 shows that a combined system has a higher geographic stability in terms of limiting the unintentional transformation of adjacent wildtype populations for a wide range of values of transgenic homozygote fitness ν (‘Medea only’ exhibits limited geographic control if fitness costs of being transgenic are high ). Interestingly, in combined systems geographic stability does not increase monotonically with respects to ν. This can result in maximal geographic stability for combined systems at intermediate values of ν (Figure 6). The levels of sustained migration, which maintain geographic stability, can be surprisingly high and of an order expected between highly interconnected demes rather than between isolated populations . In addition to the obvious regulatory benefits of robust geographic stability, this property can be exploited to limit the number of transgenic individuals in which unintended mutational events can occur or to lower the probability that pathogens evolve resistance to linked refractory genes.
Discussion and conclusion
In the theoretical analysis of combined population transformation systems Huang et al. considered the combination of a transgenic two-locus form of underdominance (termed engineered underdominance  with two other natural phenomena (Wolbachia and sex-linked meiotic drive). Both Wolbachia and sex-linked meiotic drive were demonstrated to have the potential to significantly impact the feasibility and dynamics of population transformation in both positive and negative ways. It was clearly shown that intuitive expectations of combined systems could be misleading and that mathematical modeling was essential in identifying potentially useful combinations and parameter values (most notably those relating to genotypic fitness). An excellent example is the Huang et al. theoretical analysis of the two-locus form of engineered underdominance which has been only recently realised .
Here we have followed an analogous approach to explore the properties of combining two currently developed transgenic drive systems within a single autosomal construct. The underdominant and Medea systems are assumed to be physically interspersed in a manner that maximizes the probability that they remain linked (e.g. in a configuration analogous to that shown in Figure 2). The described modeling framework has allowed us to identify a broad parameter space where combined systems can in some circumstances outperform single systems in terms of (i) optimizing release thresholds (Figures 3, 4, 7) (ii) increasing the speed of population transformation and (Figure 5) (iii) enhancing the geographic stability of population transformation (Figure 6). In addition, the reliance on two distinct mechanisms for population transformation could reduce the probability that resistance to the transgenic construct arises in the target insects. If however, long term selective pressures within successfully transformed target populations would result in the loss of the underdominance mechanism, this essentially leaves a ‘Medea only’ construct at high frequency. This ‘Medea only’ construct would be impractical to remove (unless it was associated with a high fitness cost) and could spread to adjacent populations. Conversely, loss of the Medea mechanism from a combined construct has a considerably smaller impact on reversibility and stability (Figures 2 and 6). Recognizing that the loss of Medea is preferable to loss of underdominance, it would be prudent to engineer underdominance which is more mutationally stable than Medea (duplicating the underdominant mechanism would be one simple strategy). It is also noteworthy that many of the synergistic enhancements ascribed to combined systems are to a significant extent shared by Medea constructs inserted on sex-chromosomes . Consequently, depending on the empirical properties of autosomal versus sex-chromosome inserts the relative merits of both approaches would warrant evaluation within the specific objectives of a given program.
It has been assumed throughout that fitness costs are directly associated with the drive mechanism or mechanisms in a transgenic construct, however it is also likely that additional costs will also be associated with anti-Plasmodial or anti-viral genes included as part of a working construct. The analytical framework described here will permit the prediction of the properties of combined systems loaded with such disease refectory genes. The fitness cost of refractory genes has in some, but not all, circumstances been estimated to be quite high . Consequently the illustrative parameters values used in Figure 2 may represent plausible values for ‘loaded’ constructs (though the framework presented here allows exploration of the entire range of parameters). The immediate practical use of this method could help protect D. melanogaster from an unintended species wide Medea transformation if combined with underdominance for testing in the lab. The most likely application of population transformation is in species of the genera Anopheles and Aedes which act as devastating disease vectors . Within these genera there are significant differences in dispersal capacities estimated at various locations, in some instances individuals migrate hundreds of meters over their lifetime . Consequently, the capacity to restrict transgenic constructs to particular populations is likely to be considered of high value. Various configurations of underdominance have been proposed as representing the most likely system to maintain geographic stability [12, 29]. Geographic stability is generally achieved by maximizing the fitness of transgenic homozygotes fitness ν. However where this is not possible due to cost arising from the underdominant drive mechanism or of refractory genes, our analysis indicates that maximal geographic stability can be achieved by combining systems for intermediate values of ν (Figure 6). Exploitation of this phenomena, in addition to the value of functional redundancy in drive mechanisms, could provide a valuable practical incentive to explore combined drive systems experimentally.
Average genotype fitnesses and calculating the equilibria
The frequencies of the genotypes in the next generation are given by x′, y′ and z′. In equilibrium we have x′=x, y′=y and z′=z. However the expressions for the next generation frequencies are rational functions given by, , and where the fitnesses of the genotypes are given by,
Now in equilibrium the frequencies of the genotypes do not change over generations but it is a consequence of their average fitnesses being the same. Hence we can deduce the equilibria of the system just be equating the average fitnesses. This is just another way of writing x′=x, y′=y and z′=z, which reduces to . Considering the average fitness of the genotypes in a pairwise fashion, two genotypes are neither increasing or decreasing relative to each other if their average fitnesses are equal, e.g., f x =f y . This argument is obvious when we view the system in continuous time. While the recursion equations predict the dynamics of the system in the next time step, one at a time, we can explore the complete dynamics by analysing the analogous differential equations given by,
where the time derivative of a variable is given by and so forth for y and z. From the form of these differential equations the equilibrial solutions are evident, either when the frequencies are zero (vertices of the simplexes in Figure 2) or when the bracketed terms are zero. Since the genotype frequencies sum up to 1, we can solve for just two frequencies. The solutions obtained though are complicated expressions with a possibility of imaginary roots.
Assuming complete Medea lethality (d=1), the equilibrium of the system is given by,
When it exists () then it is always unstable. Of particular interest is the case where the Medea allele frequency (given by ) at the unstable equilibrium is 0.5,
as this is ideal from the point of view of reversibility (see Figure 3). In order to cross an unstable equilibrium threshold, to ultimately transform a population, releases have to be made of a minimum size of p/(1−p) relative to the wild population size. To cross this boundary and then recross it (i.e. if we wish to reverse a completely transformed population) requires two releases with a minimum combined size of 1/p(1−p)−2. This function approaches positive infinity at p=0 and p=1 and has a minimum at p=1/2 with the release ratio being twice that of the wild population (Figure 3). Thus, an unstable threshold of p=1/2 is ideal from the perspective of population transformation and reversibility and is still much lower than release sizes used in successful applications of the sterile insect technique. Substituting the equilibrium values in Eqs. (A.3) into Eq. (A.4) gives
at d=1 (Figure 7). For t<1, there may exist two internal equilibria, the lower allele frequency one is unstable and the higher frequency one is stable, given by
However, in case of underdominance if Eq. (A.5) holds then only the unstable internal equilibrium exists at p=1/2. This is graphically illustrated in Figure 7.
Sinkins SP, Gould F: Gene drive systems for insect disease vectors. Nat Rev Gen. 2006, 7: 427-435. 10.1038/nrg1870.
Ito J, Ghosh A, Moreira LA, Wimmer EA, Jacobs-Lorena M: Transgenic anopheline mosquitoes impaired in transmission of a malaria parasite. Nature. 2002, 417: 452-455. 10.1038/417452a.
Franz AWE, Sanchez-Vargas I, Adelman ZN, Blair CD, Beaty BJ, James AA, Olson KE: Engineering RNA interference-based resistance to dengue virus type-2 in genetically-modified aedes aegypti. Proc Natl Acad Sci USA. 2006, 103: 4198-4203. 10.1073/pnas.0600479103.
Jasinskiene N, Coleman J, Ashikyan A, Salampessy M, Marinotti O: Genetic control of malaria parasite transmission: threshold levels for infection in an avian model system. Am J Trop Med Hyg. 2007, 76: 1072-1078.
Corby-Harris V, Drexler A, Watkins de Jong, Antonova Y, Pakpour N, Ziegler R, Ramberg F, Lewis E, Brown JM, Luckhart S, Riehle MA: Activation of Akt signaling reduces the prevalence and intensity of malaria parasite infection and lifespan in Anopheles stephensi mosquitoes. PLoS Pathog. 2010, 6: 1001003-10.1371/journal.ppat.1001003.
Burt A, Trivers R: Genes in Conflict. 2006, Cambridge: Harvard University Press
Hay BA, Chen C-H, Ward CM, Huang H, Su JT, Guo M: Engineering the genomes of wild insect popualtions: Challenges, and opportunities provided by synthetic medea selfish genetic elements. J Insect Physiol. 2010, 56: 1402-1413. 10.1016/j.jinsphys.2010.05.022.
Marshall JM: The effect of gene drive on containment of transgenic mosquitoes. J Theor Biol. 2009, 258: 250-265. 10.1016/j.jtbi.2009.01.031.
Windbichler N, Menichelli M, Papathanos PAPA, Thyme SBSB, Li H, Ulge UYUY, Hovde BTBT, Baker D, Monnat RJRJ, Burt A, Crisanti A: A synthetic homing endonuclease-based gene drive system in the human malaria mosquito. Nature. 2011, 473 (7): 212-215.
Chen CH, Huang H, Ward CM, Su JT, Schaeffer LV, Guo M, Hay BA: A synthetic maternal-effect selfish genetic element drives population replacement in drosophila. Science. 2007, 316: 597-600. 10.1126/science. 1138595.
Akbari OS, Chen C-H, Marshall JM, Huang H, Antoshechkin I, Hay BA: Novel synthetic medea selfish genetic elements drive population replacement in drosophila; a theoretical exploration of medea-dependent population suppression. ACS Synthetic Biol. 2012, doi:10.1021/sb300079h,
Akbari OS, Matzen KD, Marshall JM, Huang H, Ward CM, Hay BA: A synthetic gene drive system for local, reversible modification and suppression of insect populations. Curr Biol. 2013, 23 (8): 671-677. 10.1016/j.cub.2013.02.059.
Ward CM, Su JT, Huang Y, Lloyd AL, Gould F, Hay BA: Medea selfish genetic elements as tools for altering traits of wild populations: a theoretical analysis. Evolution. 2010, 65 (4): 1149-1162.
Huang Y, Magoria K, Lloyd AL, Gould F: Introducing transgenes into insect populations using combined gene-drive strategies: modeling and analysis. Insect Biochem Mol Biol. 2007, 37: 1054-1063. 10.1016/j.ibmb.2007.06.002.
Beeman RW, Friesen KS, Denell RE: Maternal-effect selfish genes in flour beetles. Science. 1992, 256: 89-92. 10.1126/science.1566060.
Peters LL, Barker JE: Novel inheritance of the murine severe combined anemia and thrombocytopenia (scat) phenotype. Cell. 1993, 74: 135-142. 10.1016/0092-8674(93)90301-6.
Weichenhan D, Traut W, Kunze B, Winking H: Distortion of mendelian recovery ratio for a mouse hsr is caused by maternal and zygotic effects. Genet Res. 1996, 68: 125-129. 10.1017/S0016672300034017.
Wade MJ, Beeman RW: The population dynamics of maternal-effect selfish genes. Genetics. 1994, 138: 1309-1314.
Haldane JBS: Selection against heterozygosis in man. Ann Hum Genet. 1941, 154 (1): 333-340.
Rasgon JL: Multi-locus assortment (mla) for transgene dispersal and elimination in mosquito populations. PloS one. 2009, 4 (6): 5833-10.1371/journal.pone.0005833.
Altrock PM, Traulsen A, Reeves RG, Reed FA: Using underdominance to bi-stably transform local populations. J Theor Biol. 2010, 267: 62-75. 10.1016/j.jtbi.2010.08.004.
Hofbauer J, Schuster P, Sigmund K: Game dynamics in mendelian populations. Biol Cybern. 1982, 43: 51-57. 10.1007/BF00337287.
Cressman R: Evolutionary Dynamics and Extensive Form Games. 2003, Cambridge: MIT Press
Dyck VA, Hendrichs J, Robinson AS: Sterile Insect Technique: Principles and Practice in Area-wide Integrated Pest Management. 2005, Dordrecht: Springer
Curtis CF: Possible use of translocations to fix desirable genes in insect pest populations. Nature. 1968, 218: 368-369. 10.1038/218368a0.
Altrock PM, Traulsen A, Reed FA: Stability properties of underdominance in finite subdivided populations. PLoS Comput Biol. 2011, 7: 1002260-10.1371/journal.pcbi.1002260.
Huang Y, Lloyd AL, Legros M, Gould F: Gene-drive into insect populations with age and spatial structure: a theoretical assessment. Evol Appl. 2011, 4 (3): 415-428. 10.1111/j.1752-4571.2010.00153.x.
Wang J: Application of the one-migrant-per-generation rule to conservation and management. Conserv Biol. 2004, 18 (2): 332-343. 10.1111/j.1523-1739.2004.00440.x.
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. 10.1006/jtbi.2001.2357.
Isaacs AT, Jasinskiene N, Tretiakov M, Thiery I, Zettor A, Bourgouin C, James AA: Transgenic anopheles stephensi coexpressing single-chain antibodies resist plasmodium falciparum development. Proc Natl Acad Sci USA. 2012, 109 (28): 1299-1230.
Gillies MT: Studies on the dispersion and survival of Anopheles gambiae Giles in East Africa, by means of marking and release experiments. Bull Entomol Res. 1961, 52 (01): 99-10.1017/S0007485300055309.
We thank P. M. Altrock, J. F. Baines, J. Denton, and A. Traulsen, for discussions and V. L. Reed for comments on the manuscript. CSG is supported by the Emmy-Noether program of the Deutsche Forschungsgemeinschaft. RGR is supported by grant RE-3062/2-1 of the Deutsche Forschungsgemeinschaft. FAR is supported by the College of Natural Sciences, University of Hawai‘i at Mānoa, by the Max Planck Society, and by a grant from the Hawai‘i Community Foundation.
The authors declare that they have no competing interests.
CSG, RGR and FAR conceived the project. All authors developed the model. CSG and FAR performed simulations. All authors analysed the results and wrote the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Gokhale, C.S., Reeves, R.G. & Reed, F.A. Dynamics of a combined medea-underdominant population transformation system. BMC Evol Biol 14, 98 (2014). https://doi.org/10.1186/1471-2148-14-98