Harnessing Wolbachia cytoplasmic incompatibility alleles for confined gene drive: A modeling study

Cleave and Rescue gamete killers create conditions for gene drive in plants

Gene drive elements promote the spread of linked traits and can be used to change the composition or fate of wild populations. Cleave and Rescue (ClvR) drive elements sit at a fixed chromosomal position and include a DNA sequence-modifying enzyme such as Cas9/gRNAs that disrupts endogenous versions of an essential gene and a recoded version of the essential gene resistant to cleavage. ClvR spreads by creating conditions in which those lacking ClvR die because they lack functional versions of the essential gene. Here we demonstrate the essential features of the ClvR gene drive in the plant Arabidopsis thaliana through killing of gametes that fail to inherit a ClvR that targets the essential gene YKT61. Resistant alleles, which can slow or prevent drive, were not observed. Modelling shows plant ClvRs are robust to certain failure modes and can be used to rapidly drive population modification or suppression. Possible applications are discussed.

Performance characteristics allow for confinement of a CRISPR toxin-antidote gene drive for population suppression in a reaction-diffusion model

Gene drive alleles that can bias their own inheritance could engineer populations for control of disease vectors, invasive species and agricultural pests. There are successful examples of suppression drives and confined modification drives, but developing confined suppression drives has proven more difficult. However, CRISPR-based toxin-antidote dominant embryo (TADE) suppression drive may fill this niche. It works by targeting and disrupting a haplolethal target gene in the germline with its gRNAs while rescuing this target. It also disrupts a female fertility gene by driving insertion or additional gRNAs. Here, we used a reaction-diffusion model to assess drive performance in continuous space, where outcomes can be substantially different from those in panmictic populations. We measured drive wave speed and found that moderate fitness costs or target gene disruption in the early embryo from maternally deposited nuclease can eliminate the drive's ability to form a wave of advance. We assessed the required release size, and finally we investigated migration corridor scenarios. It is often possible for the drive to suppress one population and then persist in the corridor without invading the second population, a potentially desirable outcome. Thus, even imperfect variants of TADE suppression drive may be excellent candidates for confined population suppression.

Modelling daisy quorum drive: A short-term bridge across engineered fitness valleys

Engineered gene-drive techniques for population modification and/or suppression have the potential for tackling complex challenges, including reducing the spread of diseases and invasive species. Gene-drive systems with low threshold frequencies for invasion, such as homing-based gene drive, require initially few transgenic individuals to spread and are therefore easy to introduce. The self-propelled behavior of such drives presents a double-edged sword, however, as the low threshold can allow transgenic elements to expand beyond a target population. By contrast, systems where a high threshold frequency must be reached before alleles can spread-above a fitness valley-are less susceptible to spillover but require introduction at a high frequency. We model a proposed drive system, called "daisy quorum drive," that transitions over time from a low-threshold daisy-chain system (involving homing-based gene drive such as CRISPR-Cas9) to a high-threshold fitness-valley system (requiring a high frequency-a "quorum"-to spread). The daisy-chain construct temporarily lowers the high thresholds required for spread of the fitness-valley construct, facilitating use in a wide variety of species that are challenging to breed and release in large numbers. Because elements in the daisy chain only drive subsequent elements in the chain and not themselves and also carry deleterious alleles ("drive load"), the daisy chain is expected to exhaust itself, removing all CRISPR elements and leaving only the high-threshold fitness-valley construct, whose spread is more spatially restricted. Developing and analyzing both discrete patch and continuous space models, we explore how various attributes of daisy quorum drive affect the chance of modifying local population characteristics and the risk that transgenic elements expand beyond a target area. We also briefly explore daisy quorum drive when population suppression is the goal. We find that daisy quorum drive can provide a promising bridge between gene-drive and fitness-valley constructs, allowing spread from a low frequency in the short term and better containment in the long term, without requiring repeated introductions or persistence of CRISPR elements.

Cleave and Rescue gamete killers create conditions for gene drive in plants

Gene drive elements promote the spread of linked traits, even when their presence confers a fitness cost to carriers, and can be used to change the composition or fate of wild populations. Cleave and Rescue (ClvR) drive elements sit at a fixed chromosomal position and include a DNA sequence-modifying enzyme such as Cas9/gRNAs (the Cleaver/Toxin) that disrupts endogenous versions of an essential gene, and a recoded version of the essential gene resistant to cleavage (the Rescue/Antidote). ClvR spreads by creating conditions in which those lacking ClvR die because they lack functional versions of the essential gene. We demonstrate the essential features of ClvR gene drive in the plant Arabidopsis thaliana through killing of gametes that fail to inherit a ClvR that targets the essential gene YKT61, whose expression is required in male and female gametes for their survival. Resistant (uncleavable but functional) alleles, which can slow or prevent drive, were not observed. Modeling shows plant ClvRs are likely to be robust to certain failure modes and can be used to rapidly drive population modification or suppression. Possible applications in plant breeding, weed control, and conservation are discussed.

Transgenic expression of cif genes from Wolbachia strain wAlbB recapitulates cytoplasmic incompatibility in Aedes aegypti

The endosymbiotic bacteria Wolbachia can invade insect populations by modifying host reproduction through cytoplasmic incompatibility (CI), an effect that results in embryonic lethality when Wolbachia-carrying males mate with Wolbachia-free females. Here we describe a transgenic system for recreating CI in the major arbovirus vector Aedes aegypti using CI factor (cif) genes from wAlbB, a Wolbachia strain currently being deployed to reduce dengue transmission. CI-like sterility is induced when cifA and cifB are co-expressed in testes; this sterility is rescued by maternal cifA expression, thereby reproducing the pattern of Wolbachia-induced CI. Expression of cifB alone is associated with extensive DNA damage and disrupted spermatogenesis. The strength of rescue by maternal cifA expression is dependent on the comparative levels of cifA/cifB expression in males. These findings are consistent with CifB acting as a toxin and CifA as an antitoxin, with CifA attenuating CifB toxicity in both the male germline and in developing embryos. These findings provide important insights into the interactions between cif genes and their mechanism of activity and provide a foundation for the building of a cif gene-based drive system in Ae. aegypti.

Making waves: Comparative analysis of gene drive spread characteristics in a continuous space model

With their ability to rapidly increase in frequency, gene drives can be used to modify or suppress target populations after an initial release of drive individuals. Recent advances have revealed many possibilities for different types of drives, and several of these have been realized in experiments. These drives have advantages and disadvantages related to their ease of construction, confinement and capacity to be used for modification or suppression. Though characteristics of these drives have been explored in modelling studies, assessment in continuous space environments has been limited, often focusing on outcomes rather than fundamental properties. Here, we conduct a comparative analysis of many different gene drive types that have the capacity to form a wave of advance in continuous space using individual-based simulations in continuous space. We evaluate the drive wave speed as a function of drive performance and ecological parameters, which reveals substantial differences between drive performance in panmictic versus spatial environments. In particular, we find that suppression drive waves are uniquely vulnerable to fitness costs and undesired CRISPR cleavage activity in embryos by maternal deposition. Some drives, however, retain robust performance even with widely varying efficiency parameters. To gain a better understanding of drive waves, we compare their panmictic performance and find that the rate of wild-type allele removal is correlated with drive wave speed, though this is also affected by other factors. Overall, our results provide a useful resource for understanding the performance of drives in spatially continuous environments, which may be most representative of potential drive deployment in many relevant scenarios.

wMel replacement of dengue-competent mosquitoes is robust to near-term change

Rising temperatures are impacting the range and prevalence of mosquito-borne diseases. A promising biocontrol technology replaces wild mosquitoes with those carrying the virus-blocking Wolbachia bacterium. Because the most widely used strain, wMel, is adversely affected by heat stress, we examined how global warming may influence wMel-based replacement. We simulated interventions in two locations with successful field trials using Coupled Model Intercomparison Project Phase 5 climate projections and historical temperature records, integrating empirical data on wMel's thermal sensitivity into a model of Aedes aegypti population dynamics to evaluate introgression and persistence over one year. We show that in Cairns, Australia, climatic futures necessitate operational adaptations for heatwaves exceeding two weeks. In Nha Trang, Vietnam, projected heatwaves of three weeks and longer eliminate wMel under the most stringent assumptions of that symbiont's thermal limits. We conclude that this technology is generally robust to near-term (2030s) climate change. Accelerated warming may challenge this in the 2050s and beyond.

Modeling emergence of Wolbachia toxin-antidote protein functions with an evolutionary algorithm

Evolutionary algorithms (EAs) simulate Darwinian evolution and adeptly mimic natural evolution. Most EA applications in biology encode high levels of abstraction in top-down population ecology models. In contrast, our research merges protein alignment algorithms from bioinformatics into codon based EAs that simulate molecular protein string evolution from the bottom up. We apply our EA to reconcile a problem in the field of Wolbachia induced cytoplasmic incompatibility (CI). Wolbachia is a microbial endosymbiont that lives inside insect cells. CI is conditional insect sterility that operates as a toxin antidote (TA) system. Although, CI exhibits complex phenotypes not fully explained under a single discrete model. We instantiate in-silico genes that control CI, CI factors (cifs), as strings within the EA chromosome. We monitor the evolution of their enzymatic activity, binding, and cellular localization by applying selective pressure on their primary amino acid strings. Our model helps rationalize why two distinct mechanisms of CI induction might coexist in nature. We find that nuclear localization signals (NLS) and Type IV secretion system signals (T4SS) are of low complexity and evolve fast, whereas binding interactions have intermediate complexity, and enzymatic activity is the most complex. Our model predicts that as ancestral TA systems evolve into eukaryotic CI systems, the placement of NLS or T4SS signals can stochastically vary, imparting effects that might impact CI induction mechanics. Our model highlights how preconditions and sequence length can bias evolution of cifs toward one mechanism or another.

Assessment of distant-site rescue elements for CRISPR toxin-antidote gene drives

Gene drive is a genetic engineering technology that can enable super-mendelian inheritance of specific alleles, allowing them to spread through a population. New gene drive types have increased flexibility, offering options for confined modification or suppression of target populations. Among the most promising are CRISPR toxin-antidote gene drives, which disrupt essential wild-type genes by targeting them with Cas9/gRNA. This results in their removal, increasing the frequency of the drive. All these drives rely on having an effective rescue element, which consists of a recoded version of the target gene. This rescue element can be at the same site as the target gene, maximizing the chance of efficient rescue, or at a distant site, which allows useful options such as easily disrupting another essential gene or increasing confinement. Previously, we developed a homing rescue drive targeting a haplolethal gene and a toxin-antidote drive targeting a haplosufficient gene. These successful drives had functional rescue elements but suboptimal drive efficiency. Here, we attempted to construct toxin-antidote drives targeting these genes with a distant-site configuration from three loci in Drosophila melanogaster. We found that additional gRNAs increased cut rates to nearly 100%. However, all distant-site rescue elements failed for both target genes. Furthermore, one rescue element with a minimally recoded sequence was used as a template for homology-directed repair for the target gene on a different chromosomal arm, resulting in the formation of functional resistance alleles. Together, these results can inform the design of future CRISPR-based toxin-antidote gene drives.