Modeling the Dynamics of a Non-Limited and a Self-Limited Gene Drive System in Structured Aedes aegypti Populations

Recently there have been significant advances in research on genetic strategies to control populations of disease-vectoring insects. Some of these strategies use the gene drive properties of selfish genetic elements to spread physically linked anti-pathogen genes into local vector populations. Because of the potential of these selfish elements to spread through populations, control approaches based on these strategies must be carefully evaluated to ensure a balance between the desirable spread of the refractoriness-conferring genetic cargo and the avoidance of potentially unwanted outcomes such as spread to non-target populations. There is also a need to develop better estimates of the economics of such releases. We present here an evaluation of two such strategies using a biologically realistic mathematical model that simulates the resident Aedes aegypti mosquito population of Iquitos, Peru. One strategy uses the selfish element Medea, a non-limited element that could permanently spread over a large geographic area; the other strategy relies on Killer-Rescue genetic constructs, and has been predicted to have limited spatial and temporal spread. We simulate various operational approaches for deploying these genetic strategies, and quantify the optimal number of released transgenic mosquitoes needed to achieve definitive spread of Medea-linked genes and/or high frequencies of Killer-Rescue-associated elements. We show that for both strategies the most efficient approach for achieving spread of anti-pathogen genes within three years is generally to release adults of both sexes in multiple releases over time. Even though females in these releases should not transmit disease, there could be public concern over such releases, making the less efficient male-only release more practical. This study provides guidelines for operational approaches to population replacement genetic strategies, as well as illustrates the use of detailed spatial models to assist in safe and efficient implementation of such novel genetic strategies.     

Designing gene drives to limit spillover to non-target populations

The prospect of utilizing CRISPR-based gene-drive technology for controlling populations has generated much excitement. However, the potential for spillovers of gene-drive alleles from the target population to non-target populations has raised concerns. Here, using mathematical models, we investigate the possibility of limiting spillovers to non-target populations by designing differential-targeting gene drives, in which the expected equilibrium gene-drive allele frequencies are high in the target population but low in the non-target population. We find that achieving differential targeting is possible with certain configurations of gene-drive parameters, but, in most cases, only under relatively low migration rates between populations. Under high migration, differential targeting is possible only in a narrow region of the parameter space. Because fixation of the gene drive in the non-target population could severely disrupt ecosystems, we outline possible ways to avoid this outcome. We apply our model to two potential applications of gene drives—field trials for malaria-vector gene drives and control of invasive species on islands. We discuss theoretical predictions of key requirements for differential targeting and their practical implications.     

Development of a targeted transgenesis strategy in highly differentiated cells: a powerful tool for functional genomic analysis

Functional genomic analysis is a challenging step in the so-called post-genomic field. Identification of potential targets using large-scale gene expression analysis requires functional validation to identify those that are physiologically relevant. Genetically modified cell models are often used for this purpose allowing up- or down-expression of selected targets in a well-defined and if possible highly differentiated cell type. However, the generation of such models remains time-consuming and expensive. In order to alleviate this step, we developed a strategy aimed at the rapid and efficient generation of genetically modified cell lines with conditional, inducible expression of various target genes. Efficient knock-in of various constructs, called targeted transgenesis, in a locus selected for its permissibility to the tet inducible system, was obtained through the stimulation of site-specific homologous recombination by the meganuclease I-SceI. Our results demonstrate that targeted transgenesis in a reference inducible locus greatly facilitated the functional analysis of the selected recombinant cells. The efficient screening strategy we have designed makes possible automation of the transfection and selection steps. Furthermore, this strategy could be applied to a variety of highly differentiated cells.     

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

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

Site-specific selfish genes as tools for the control and genetic engineering of natural populations

Site-specific selfish genes exploit host functions to copy themselves into a defined target DNA sequence, and include homing endonuclease genes, group II introns and some LINE-like transposable elements. If such genes can be engineered to target new host sequences, then they can be used to manipulate natural populations, even if the number of individuals released is a small fraction of the entire population. For example, a genetic load sufficient to eradicate a population can be imposed in fewer than 20 generations, if the target is an essential host gene, the knockout is recessive and the selfish gene has an appropriate promoter. There will be selection for resistance, but several strategies are available for reducing the likelihood of it evolving. These genes may also be used to genetically engineer natural populations, by means of population-wide gene knockouts, gene replacements and genetic transformations. By targeting sex-linked loci just prior to meiosis one may skew the population sex ratio, and by changing the promoter one may limit the spread of the gene to neighbouring populations. The proposed constructs are evolutionarily stable in the face of the mutations most likely to arise during their spread, and strategies are also available for reversing the manipulations.     

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

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

Inverse Medea as a novel gene drive system for local population replacement: a theoretical analysis

One strategy to control mosquito-borne diseases, such as malaria and dengue fever, on a regional scale is to use gene drive systems to spread disease-refractory genes into wild mosquito populations. The development of a synthetic Medea element that has been shown to drive population replacement in laboratory Drosophila populations has provided encouragement for this strategy but has also been greeted with caution over the concern that transgenes may spread into countries without their consent. Here, we propose a novel gene drive system, inverse Medea, which is strong enough to bring about local population replacement but is unable to establish itself beyond an isolated release site. The system consists of 2 genetic components--a zygotic toxin and maternal antidote--which render heterozygous offspring of wild-type mothers unviable. Through population genetic analysis, we show that inverse Medea will only spread when it represents a majority of the alleles in a population. The element is best located on an autosome and will spread to fixation provided any associated fitness costs are dominant and to very high frequency otherwise. We suggest molecular tools that could be used to build the inverse Medea system and discuss its utility for a confined release of transgenic mosquitoes.     

Genomic differentiation and patterns of gene flow between two long‐tailed tit species (Aegithalos)

Patterns of heterogeneous genomic differentiation have been well documented between closely related species, with some highly differentiated genomic regions (“genomic differentiation islands”) spread throughout the genome. Differential levels of gene flow are proposed to account for this pattern, as genomic differentiation islands are suggested to be resistant to gene flow. Recent studies have also suggested that genomic differentiation islands could be explained by linked selection acting on genomic regions with low recombination rates. Here, we investigate genomic differentiation and gene‐flow patterns for autosomes using RAD‐seq data between two closely related species of long‐tailed tits (Aegithalos bonvaloti and A. fuliginosus) in both allopatric and contact zone populations. The results confirm recent or ongoing gene flow between these two species. However, there is little evidence that the genomic regions that were found to be highly differentiated between the contact zone populations are resistant to gene flow, suggesting that differential levels of gene flow is not the cause of the heterogeneous genomic differentiation. Linked selection may be the cause of genomic differentiation islands between the allopatric populations with no or very limited gene flow, but this could not account for the heterogeneous genomic differentiation between the contact zone populations, which show evidence of recent or ongoing gene flow.     

Gene drive systems in mosquitoes: rules of the road

Population replacement strategies for controlling transmission of mosquito-borne diseases call for the introgression of antipathogen effector genes into vector populations. It is anticipated that these genes, if present at high enough frequencies, will impede transmission of the target pathogens and result in reduced human morbidity and mortality. Recent laboratory successes in the development of virus- and protozoan-resistant mosquito strains make urgent research of gene drive systems capable of moving effector genes into wild populations. A systematic approach to developing safe and effective gene drive systems that includes defining the requirements of the system, identifying naturally occurring or synthetic genetic mechanisms for gene spread upon which drive systems can be based and the successful adaptation of a mechanism to a drive system, should mitigate concerns about using genetically engineered mosquitoes for disease control.     

What Magnitude Are Observed Non-Target Impacts from Weed Biocontrol?

A systematic review focused by plant on non-target impacts from agents deliberately introduced for the biological control of weeds found significant non-target impacts to be rare. The magnitude of direct impact of 43 biocontrol agents on 140 non-target plants was retrospectively categorized using a risk management framework for ecological impacts of invasive species (minimal, minor, moderate, major, massive). The vast majority of agents introduced for classical biological control of weeds (>99% of 512 agents released) have had no known significant adverse effects on non-target plants thus far; major effects suppressing non-target plant populations could be expected to be detectable. Most direct non-target impacts on plants (91.6%) were categorized as minimal or minor in magnitude with no known adverse long-term impact on non-target plant populations, but a few cacti and thistles are affected at moderate (n = 3), major (n = 7) to massive (n = 1) scale. The largest direct impacts are from two agents (Cactoblastis cactorum on native cacti and Rhinocyllus conicus on native thistles), but these introductions would not be permitted today as more balanced attitudes exist to plant biodiversity, driven by both society and the scientific community. Our analysis shows (as far as is known), weed biological control agents have a biosafety track record of >99% of cases avoiding significant non-target impacts on plant populations. Some impacts could have been overlooked, but this seems unlikely to change the basic distribution of very limited adverse effects. Fewer non-target impacts can be expected in future because of improved science and incorporation of wider values. Failure to use biological control represents a significant opportunity cost from the certainty of ongoing adverse impacts from invasive weeds. It is recommended that a simple five-step scale be used to better communicate the risk of consequences from both action (classical biological control) and no action (ongoing impacts from invasive weeds).     

Testing transposable elements as genetic drive mechanisms using Drosophila P element constructs as a model system

The use of transposable elements (TEs) as genetic drive mechanisms was explored using Drosophila melanogaster as a model system. Alternative strategies, employing autonomous and nonautonomous P element constructs were compared for their efficiency in driving the ry⁺ allele into populations homozygous for a ry^- allele at the genomic rosy locus. Transformed flies were introduced at 1%, 5%, and 10% starting frequencies to establish a series of populations that were monitored over the course of 40 generations, using both phenotypic and molecular assays. The transposon-borne ry⁺ marker allele spread rapidly in almost all populations when introduced at 5% and 10% seed frequencies, but 1% introductions frequently failed to become established. A similar initial rapid increase in frequency of the ry⁺ transposon occurred in several control populations lacking a source of transposase. Constructs carrying ry⁺ markers also increased to moderate frequencies in the absence of selection on the marker. The results of Southern and in situ hybridization studies indicated a strong inverse relationship between the degree of conservation of construct integrity and transposition frequency. These finding have relevance to possible future applications of transposons as genetic drive mechanisms. This revised version was published online in July 2006 with corrections to the Cover Date.     

nanos-Driven expression of piggyBac transposase induces mobilization of a synthetic autonomous transposon in the malaria vector mosquito,Anopheles stephensi

Transposons are a class of selfish DNA elements that can mobilize within a genome. If mobilization is accompanied by an increase in copy number (replicative transposition), the transposon may sweep through a population until it is fixed in all of its interbreeding members. This introgression has been proposed as the basis for drive systems to move genes with desirable phenotypes into target species. One such application would be to use them to move a gene conferring resistance to malaria parasites throughout a population of vector mosquitos. We assessed the feasibility of using the piggyBac transposon as a gene-drive mechanism to distribute anti-malarial transgenes in populations of the malaria vector, Anopheles stephensi. We designed synthetic gene constructs that express the piggyBac transposase in the female germline using the control DNA of the An. stephensi nanos orthologous gene linked to marker genes to monitor inheritance. Two remobilization events were observed with a frequency of one every 23 generations, a rate far below what would be useful to drive anti-pathogen transgenes into wild mosquito populations. We discuss the possibility of optimizing this system and the impetus to do so.     

Use of Wolbachia to drive nuclear transgenes through insect populations

Wolbachia is an inherited intracellular bacterium found in many insects of medical and economic importance. The ability of many strains to spread through populations using cytoplasmic incompatibility, involving sperm modification and rescue, provides a powerful mechanism for driving beneficial transgenes through insect populations, if such transgenes could be inserted into and expressed by Wolbachia. However, manipulating Wolbachia in this way has not yet been achieved. Here, we demonstrate theoretically an alternative mechanism whereby nuclear rather than cytoplasmic transgenes could be driven through populations, by linkage to a nuclear gene able to rescue modified sperm. The spread of a 'nuclear rescue construct' occurs as long as the Wolbachia show imperfect maternal transmission under natural conditions and/or imperfect rescue of modified sperm. The mechanism is most efficient when the target population is already infected with Wolbachia at high frequency, whether naturally or by the sequential release of Wolbachia-infected individuals and subsequently the nuclear rescue construct. The results provide a potentially powerful addition to the few insect transgene drive mechanisms that are available.     

Modelling Allee effects in a transgenic mosquito population during range expansion

ABSTRACT

Mosquitoes are vectors for many diseases that cause significant mortality and morbidity. As mosquito populations expand their range, they may undergo mate-finding Allee effects such that their ability to successfully reproduce becomes difficult at low population density. With new technology, creating target specific gene modification may be a viable method for mosquito population control. We develop a mathematical model to investigate the effects of releasing transgenic mosquitoes into newly established, low-density mosquito populations. Our model consists of two life stages (aquatic and adults), which are divided into three genetically distinct groups: heterogeneous and homogeneous transgenic that cause female infertility and a homogeneous wild type. We perform analytical and numerical analyses on the equilibria to determine the level of saturation needed to eliminate mosquitoes in a given area. This model demonstrates the potential for a gene drive system to reduce the spread of invading mosquito populations.     

Introducing desirable transgenes into insect populations using Y-linked meiotic drive - a theoretical assessment

The use of genetic drive mechanisms to replace native mosquito genotypes with individuals bearing antipathogen transgenes is a potential strategy for repressing insect transmission of human diseases such as malaria and dengue. Antipathogen transgenes have been developed and tested, but efficient gene drive mechanisms are lacking. Here we theoretically assess the feasibility of introducing antipathogen genes into wild Aedes aegypti populations by using a naturally occurring meiotic drive system. We consider the release of males having both a Y-linked meiotic drive gene and an X-linked drive-insensitive response allele to which an antipathogen gene is linked. We use mathematical models and computer simulations to determine how the post-introduction dynamics of the antipathogen gene are affected by specific genetic characteristics of the system. The results show that when the natural population is uniformly sensitive to the meiotic drive gene, the antipathogen gene may be driven close to fixation if the fitness costs of the drive gene, the insensitive response allele, and the antipathogen gene are low. However, when the natural population has a small proportion of an X-linked insensitive response allele or an autosomal gene that strongly reduces the effect of the drive gene, the antipathogen gene does not spread if it has an associated fitness cost. Our modeling results provide a theoretical foundation for further experimental tests.     

Management of orthopteran pests: a conservation perspective

Although the vast majority of orthopterans are not pests, some species have the potential to cause serious damage to human interests. Management of pest populations frequently conflicts with conservation of orthopteran species and processes, particularly when the pest species or its ecological processes are susceptible to extinction or when the pest population is coincident with non-target orthopterans. With respect to chemical control, the greatest hazards are the broad-spectrum, highly lethal properties of most agents, which can be mitigated with formulation and application methods. Biological control risks permanent, large-scale changes to orthopteran species and processes which can be minimized with bioinsecticidal and other short-lived or selective formulations and reliable host-range testing. Cultural control may have large-scale, broad-spectrum impacts to non-target orthopterans, but these hazards can be diminished by appropriate testing and monitoring. Mechanical control methods may be impractically labour intensive, but they are highly target specific and therefore warrant further consideration. Social control measures such as education, insurance and compensation programmes appear to have little direct potential for harm to orthopteran conservation, but the complex socioeconomic and, ultimately, environmental consequences of such programmes have not been assessed. The melding of orthopteran pest management and conservation requires that we perceive these insects and their ecological processes to be vital elements of sustainable agroecosystems. Our management of orthopterans (both non-target and pest populations) must focus on keeping good stewards on the land.     

The influence of antisense gene location on target gene suppression in the fission yeast Schizosaccharomyces pombe

A fission yeast model was employed to investigate the influence of antisense gene location on the efficacy of antisense RNA-mediated target gene suppression. Fission yeast transformants were generated that contained the target lacZ gene at a fixed position and a single copy antisense lacZ gene integrated into various genomic locations, including the same locus as the target gene. No significant difference in lacZ suppression was observed when the antisense gene was integrated in close proximity to the target gene locus compared with other genomic locations, indicating that target and antisense gene colocalization is not a critical factor for efficient antisense RNA-mediated gene expression in vivo. Instead, increased lacZ downregulation correlated with an increase in antisense dose, with the steady-state levels of antisense RNA being dependent on genomic position effects and transgene copy number.     

Population Genetic Structure of the Grasshopper Eyprepocnemis plorans in the South and East of the Iberian Peninsula

The grasshopper Eyprepocnemis plorans subsp. plorans harbors a very widespread polymorphism for supernumerary (B) chromosomes which appear to have arisen recently. These chromosomes behave as genomic parasites because they are harmful for the individuals carrying them and show meiotic drive in the initial stages of population invasion. The rapid increase in B chromosome frequency at intrapopulation level is thus granted by meiotic drive, but its spread among populations most likely depends on interpopulation gene flow. We analyze here the population genetic structure in 10 natural populations from two regions (in the south and east) of the Iberian Peninsula. The southern populations were coastal whereas the eastern ones were inland populations located at 260–655 m altitude. The analysis of 97 ISSR markers revealed significant genetic differentiation among populations (average G_ST =  0.129), and the Structure software and AMOVA indicated a significant genetic differentiation between southern and eastern populations. There was also significant isolation by distance (IBD) between populations. Remarkably, these results were roughly similar to those found when only the markers showing low or no dropout were included, suggesting that allelic dropout had negligible effects on population genetic analysis. We conclude that high gene flow helped this parasitic B chromosome to spread through most of the geographical range of the subspecies E. plorans plorans.     

A proposal to use plant demographic data to assess potential weed biological control agents impacts on non-target plant populations

Weed biocontrol programs aim to reduce the spread and population growth rate of the target plant while stabilizing or increasing populations of those native species considered under threat by invasive plants. This goal is not unique to weed biocontrol but applies to all other invasive plant management techniques, though such information is rarely collected. Without this information, success of management interventions can be ambiguous, and regulatory agencies, the public, policy makers, funders and land managers cannot be held accountable for chosen treatments. A fundamental reform, including use of demographic studies and long-term assessments, are essential to guide weed biocontrol programs. We propose to add use of plant demography (an assessment of how environmental factors and ecological interactions, for example competition, disease or herbivory, may affect plant populations by altering survival, growth, development and reproductive rates of plant individuals) during host specificity risk assessments of potential biological control agents. Demographic models can refine assessments of potential impacts for those plant species that experience some feeding or larval development during host specificity testing. Our proposed approach to focus on impact on plant demography instead of attack on plant individuals is useful in appropriately gauging threats potential weed biocontrol agents may pose to non-target species after field release.