Germline Cas9 promoters with improved performance for homing gene drive.

Gene drive systems could be a viable strategy to prevent pathogen transmission or suppress vector populations by propagating drive alleles with super-Mendelian inheritance. CRISPR-based homing gene drives convert wild type alleles into drive alleles in heterozygotes with Cas9 and gRNA. It is thus desirable to identify Cas9 promoters that yield high drive conversion rates, minimize the formation rate of resistance alleles in both the germline and the early embryo, and limit somatic Cas9 expression. In Drosophila, the nanos promoter avoids leaky somatic expression, but at the cost of high embryo resistance from maternally deposited Cas9. To improve drive efficiency, we test eleven Drosophila melanogaster germline promoters. Some achieve higher drive conversion efficiency with minimal embryo resistance, but none completely avoid somatic expression. However, such somatic expression often does not carry detectable fitness costs for a rescue homing drive targeting a haplolethal gene, suggesting somatic drive conversion. Supporting a 4-gRNA suppression drive, one promoter leads to a low drive equilibrium frequency due to fitness costs from somatic expression, but the other outperforms nanos, resulting in successful suppression of the cage population. Overall, these Cas9 promoters hold advantages for homing drives in Drosophila species and may possess valuable homologs in other organisms.

Resource-explicit interactions in spatial population models.

Continuous-space population models can yield significantly different results from their panmictic counterparts when assessing evolutionary, ecological, or population-genetic processes. However, the computational burden of spatial models is typically much greater than that of panmictic models due to the overhead of determining which individuals interact with one another and how strongly they interact. Though these calculations are necessary to model local competition that regulates the population density, they can lead to prohibitively long runtimes. Here, we present a novel modeling method in which the resources available to a population are abstractly represented as an additional layer of the simulation. Instead of interacting directly with one another, individuals interact indirectly via this resource layer. We find that this method closely matches other spatial models, yet can dramatically increase the speed of the model, allowing the simulation of much larger populations. Additionally, models structured in this manner exhibit other desirable characteristics, including more realistic spatial dynamics near the edge of the simulated area, and an efficient route for modeling more complex heterogeneous landscapes.

A framework for identifying fertility gene targets for mammalian pest control.

Fertility-targeted gene drives have been proposed as an ethical genetic approach for managing wild populations of vertebrate pests for public health and conservation benefit. This manuscript introduces a framework to identify and evaluate target gene suitability based on biological gene function, gene expression and results from mouse knockout models. This framework identified 16 genes essential for male fertility and 12 genes important for female fertility that may be feasible targets for mammalian gene drives and other non-drive genetic pest control technology. Further, a comparative genomics analysis demonstrates the conservation of the identified genes across several globally significant invasive mammals. In addition to providing important considerations for identifying candidate genes, our framework and the genes identified in this study may have utility in developing additional pest control tools such as wildlife contraceptives.

Incorporating ecology into gene drive modelling.

Gene drive technology, in which fast-spreading engineered drive alleles are introduced into wild populations, represents a promising new tool in the fight against vector-borne diseases, agricultural pests and invasive species. Due to the risks involved, gene drives have so far only been tested in laboratory settings while their population-level behaviour is mainly studied using mathematical and computational models. The spread of a gene drive is a rapid evolutionary process that occurs over timescales similar to many ecological processes. This can potentially generate strong eco-evolutionary feedback that could profoundly affect the dynamics and outcome of a gene drive release. We, therefore, argue for the importance of incorporating ecological features into gene drive models. We describe the key ecological features that could affect gene drive behaviour, such as population structure, life-history, environmental variation and mode of selection. We review previous gene drive modelling efforts and identify areas where further research is needed. As gene drive technology approaches the level of field experimentation, it is crucial to evaluate gene drive dynamics, potential outcomes, and risks realistically by including ecological processes.

A framework for identifying fertility gene targets for mammalian pest control.

Fertility-targeted gene drives have been proposed as an ethical genetic approach for managing wild populations of vertebrate pests for public health and conservation benefit.This manuscript introduces a framework to identify and evaluate target gene suitability based on biological gene function, gene expression, and results from mouse knockout models.This framework identified 16 genes essential for male fertility and 12 genes important for female fertility that may be feasible targets for mammalian gene drives and other non-drive genetic pest control technology. Further, a comparative genomics analysis demonstrates the conservation of the identified genes across several globally significant invasive mammals.In addition to providing important considerations for identifying candidate genes, our framework and the genes identified in this study may have utility in developing additional pest control tools such as wildlife contraceptives.

Evolution under a model of functionally buffered deleterious mutations can lead to positive selection in protein-coding genes.

Selective pressures on DNA sequences often result in departures from neutral evolution that can be captured by the McDonald-Kreitman (MK) test. However, the nature of such selective forces often remains unknown to experimentalists. Amino acid fixations driven by natural selection in protein-coding genes are commonly associated with a genetic arms race or changing biological purposes, leading to proteins with new functionality. Here, we evaluate the expectations of population genetic patterns under a buffering mechanism driving selective amino acids to fixation, which is motivated by an observed phenotypic rescue of otherwise deleterious nonsynonymous substitutions at bag of marbles (bam) and Sex lethal (Sxl) in Drosophila melanogaster. These two genes were shown to experience strong episodic bursts of natural selection potentially due to infections of the endosymbiotic bacteria Wolbachia observed among multiple Drosophila species. Using simulations to implement and evaluate the evolutionary dynamics of a Wolbachia buffering model, we demonstrate that selectively fixed amino acid replacements will occur, but that the proportion of adaptive amino acid fixations and the statistical power of the MK test to detect the departure from an equilibrium neutral model are both significantly lower than seen for an arms race/change-in-function model that favors proteins with diversified amino acids. We find that the observed selection pattern at bam in a natural population of D. melanogaster is more consistent with an arms race model than with the buffering model.

SLiM 4: Multispecies Eco-Evolutionary Modeling.

AbstractThe SLiM software framework for genetically explicit forward simulation has been widely used in population genetics. However, it has been largely restricted to modeling only a single species, which has limited its broader utility in evolutionary biology. Indeed, to our knowledge no general-purpose, flexible modeling framework exists that provides support for simulating multiple species while also providing other key features, such as explicit genetics and continuous space. The lack of such software has limited our ability to model higher biological levels such as communities, ecosystems, coevolutionary and eco-evolutionary processes, and biodiversity, which is crucial for many purposes, from extending our basic understanding of evolutionary ecology to informing conservation and management decisions. We here announce the release of SLiM 4, which fills this important gap by adding support for multiple species, including ecological interactions between species such as predation, parasitism, and mutualism, and illustrate its new features with examples.

Anopheles homing suppression drive candidates exhibit unexpected performance differences in simulations with spatial structure.

Recent experiments have produced several Anopheles gambiae homing gene drives that disrupt female fertility genes, thereby eventually inducing population collapse. Such drives may be highly effective tools to combat malaria. One such homing drive, based on the zpg promoter driving CRISPR/Cas9, was able to eliminate a cage population of mosquitoes. A second version, purportedly improved upon the first by incorporating an X-shredder element (which biases inheritance towards male offspring), was similarly successful. Here, we analyze experimental data from each of these gene drives to extract their characteristics and performance parameters and compare these to previous interpretations of their experimental performance. We assess each suppression drive within an individual-based simulation framework that models mosquito population dynamics in continuous space. We find that the combined homing/X-shredder drive is actually less effective at population suppression within the context of our mosquito population model. In particular, the combined drive often fails to completely suppress the population, instead resulting in an unstable equilibrium between drive and wild-type alleles. By contrast, otherwise similar drives based on the nos promoter may prove to be more promising candidates for future development than originally thought.

Deep cis-regulatory homology of the butterfly wing pattern ground plan.

Butterfly wing patterns derive from a deeply conserved developmental ground plan yet are diverse and evolve rapidly. It is poorly understood how gene regulatory architectures can accommodate both deep homology and adaptive change. To address this, we characterized the cis-regulatory evolution of the color pattern gene WntA in nymphalid butterflies. Comparative assay for transposase-accessible chromatin using sequencing (ATAC-seq) and in vivo deletions spanning 46 cis-regulatory elements across five species revealed deep homology of ground plan-determining sequences, except in monarch butterflies. Furthermore, noncoding deletions displayed both positive and negative regulatory effects that were often broad in nature. Our results provide little support for models predicting rapid enhancer turnover and suggest that deeply ancestral, multifunctional noncoding elements can underlie rapidly evolving trait systems.

Fitness effects of CRISPR endonucleases in Drosophila melanogaster populations.

Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/Cas9 provides a highly efficient and flexible genome editing technology with numerous potential applications ranging from gene therapy to population control. Some proposed applications involve the integration of CRISPR/Cas9 endonucleases into an organism's genome, which raises questions about potentially harmful effects to the transgenic individuals. One example for which this is particularly relevant are CRISPR-based gene drives conceived for the genetic alteration of entire populations. The performance of such drives can strongly depend on fitness costs experienced by drive carriers, yet relatively little is known about the magnitude and causes of these costs. Here, we assess the fitness effects of genomic CRISPR/Cas9 expression in Drosophila melanogaster cage populations by tracking allele frequencies of four different transgenic constructs that allow us to disentangle 'direct' fitness costs due to the integration, expression, and target-site activity of Cas9, from fitness costs due to potential off-target cleavage. Using a maximum likelihood framework, we find that a model with no direct fitness costs but moderate costs due to off-target effects fits our cage data best. Consistent with this, we do not observe fitness costs for a construct with Cas9HF1, a high-fidelity version of Cas9. We further demonstrate that using Cas9HF1 instead of standard Cas9 in a homing drive achieves similar drive conversion efficiency. These results suggest that gene drives should be designed with high-fidelity endonucleases and may have implications for other applications that involve genomic integration of CRISPR endonucleases.

Experimental demonstration of tethered gene drive systems for confined population modification or suppression.

BACKGROUND: Homing gene drives hold great promise for the genetic control of natural populations. However, current homing systems are capable of spreading uncontrollably between populations connected by even marginal levels of migration. This could represent a substantial sociopolitical barrier to the testing or deployment of such drives and may generally be undesirable when the objective is only local population control, such as suppression of an invasive species outside of its native range. Tethered drive systems, in which a locally confined gene drive provides the CRISPR nuclease needed for a homing drive, could provide a solution to this problem, offering the power of a homing drive and confinement of the supporting drive. RESULTS: Here, we demonstrate the engineering of a tethered drive system in Drosophila, using a regionally confined CRISPR Toxin-Antidote Recessive Embryo (TARE) drive to support modification and suppression homing drives. Each drive was able to bias inheritance in its favor, and the TARE drive was shown to spread only when released above a threshold frequency in experimental cage populations. After the TARE drive had established in the population, it facilitated the spread of a subsequently released split homing modification drive (to all individuals in the cage) and of a homing suppression drive (to its equilibrium frequency). CONCLUSIONS: Our results show that the tethered drive strategy is a viable and easily engineered option for providing confinement of homing drives to target populations.

A homing suppression gene drive with multiplexed gRNAs maintains high drive conversion efficiency and avoids functional resistance alleles.

Gene drives are engineered alleles that can bias inheritance in their favor, allowing them to spread throughout a population. They could potentially be used to modify or suppress pest populations, such as mosquitoes that spread diseases. CRISPR/Cas9 homing drives, which copy themselves by homology-directed repair in drive/wild-type heterozygotes, are a powerful form of gene drive, but they are vulnerable to resistance alleles that preserve the function of their target gene. Such resistance alleles can prevent successful population suppression. Here, we constructed a homing suppression drive in Drosophila melanogaster that utilized multiplexed gRNAs to inhibit the formation of functional resistance alleles in its female fertility target gene. The selected gRNA target sites were close together, preventing reduction in drive conversion efficiency. The construct reached a moderate equilibrium frequency in cage populations without apparent formation of resistance alleles. However, a moderate fitness cost prevented elimination of the cage population, showing the importance of using highly efficient drives in a suppression strategy, even if resistance can be addressed. Nevertheless, our results experimentally demonstrate the viability of the multiplexed gRNAs strategy in homing suppression gene drives.

Predicting the genomic resolution of bulk segregant analysis.

Bulk segregant analysis is a technique for identifying the genetic loci that underlie phenotypic trait differences. The basic approach is to compare two pools of individuals from the opposing tails of the phenotypic distribution, sampled from an interbred population. Each pool is sequenced and scanned for alleles that show divergent frequencies between the pools, indicating potential association with the observed trait differences. Bulk segregant analysis has already been successfully applied to the mapping of various quantitative trait loci in organisms ranging from yeast to maize. However, these studies have typically suffered from rather low mapping resolution, and we still lack a detailed understanding of how this resolution is affected by experimental parameters. Here, we use coalescence theory to calculate the expected genomic resolution of bulk segregant analysis for a simple monogenic trait. We first show that in an idealized interbreeding population of infinite size, the expected length of the mapped region is inversely proportional to the recombination rate, the number of generations of interbreeding, and the number of genomes sampled, as intuitively expected. In a finite population, coalescence events in the genealogy of the sample reduce the number of potentially informative recombination events during interbreeding, thereby increasing the length of the mapped region. This is incorporated into our model by an effective population size parameter that specifies the pairwise coalescence rate of the interbreeding population. The mapping resolution predicted by our calculations closely matches numerical simulations and is surprisingly robust to moderate levels of contamination of the segregant pools with alternative alleles. Furthermore, we show that the approach can easily be extended to modifications of the crossing scheme. Our framework will allow researchers to predict the expected power of their mapping experiments, and to evaluate how their experimental design could be tuned to optimize mapping resolution.

Modeling CRISPR gene drives for suppression of invasive rodents using a supervised machine learning framework.

Invasive rodent populations pose a threat to biodiversity across the globe. When confronted with these invaders, native species that evolved independently are often defenseless. CRISPR gene drive systems could provide a solution to this problem by spreading transgenes among invaders that induce population collapse, and could be deployed even where traditional control methods are impractical or prohibitively expensive. Here, we develop a high-fidelity model of an island population of invasive rodents that includes three types of suppression gene drive systems. The individual-based model is spatially explicit, allows for overlapping generations and a fluctuating population size, and includes variables for drive fitness, efficiency, resistance allele formation rate, as well as a variety of ecological parameters. The computational burden of evaluating a model with such a high number of parameters presents a substantial barrier to a comprehensive understanding of its outcome space. We therefore accompany our population model with a meta-model that utilizes supervised machine learning to approximate the outcome space of the underlying model with a high degree of accuracy. This enables us to conduct an exhaustive inquiry of the population model, including variance-based sensitivity analyses using tens of millions of evaluations. Our results suggest that sufficiently capable gene drive systems have the potential to eliminate island populations of rodents under a wide range of demographic assumptions, though only if resistance can be kept to a minimal level. This study highlights the power of supervised machine learning to identify the key parameters and processes that determine the population dynamics of a complex evolutionary system.

Design and analysis of CRISPR-based underdominance toxin-antidote gene drives.

CRISPR gene drive systems offer a mechanism for transmitting a desirable transgene throughout a population for purposes ranging from vector-borne disease control to invasive species suppression. In this simulation study, we assess the performance of several CRISPR-based underdominance gene drive constructs employing toxin-antidote (TA) principles. These drives disrupt the wild-type version of an essential gene using a CRISPR nuclease (the toxin) while simultaneously carrying a recoded version of the gene (the antidote). Drives of this nature allow for releases that could be potentially confined to a desired geographic location. This is because such drives have a nonzero-invasion threshold frequency required for the drive to spread through the population. We model drives which target essential genes that are either haplosufficient or haplolethal, using nuclease promoters with expression restricted to the germline, promoters that additionally result in cleavage activity in the early embryo from maternal deposition, and promoters that have ubiquitous somatic expression. We also study several possible drive architectures, considering both "same-site" and "distant-site" systems, as well as several reciprocally targeting drives. Together, these drive variants provide a wide range of invasion threshold frequencies and options for both population modification and suppression. Our results suggest that CRISPR TA underdominance drive systems could allow for the design of flexible and potentially confinable gene drive strategies.

Detection of hard and soft selective sweeps from Drosophila melanogaster population genomic data.

Whether hard sweeps or soft sweeps dominate adaptation has been a matter of much debate. Recently, we developed haplotype homozygosity statistics that (i) can detect both hard and soft sweeps with similar power and (ii) can classify the detected sweeps as hard or soft. The application of our method to population genomic data from a natural population of Drosophila melanogaster (DGRP) allowed us to rediscover three known cases of adaptation at the loci Ace, Cyp6g1, and CHKov1 known to be driven by soft sweeps, and detected additional candidate loci for recent and strong sweeps. Surprisingly, all of the top 50 candidates showed patterns much more consistent with soft rather than hard sweeps. Recently, Harris et al. 2018 criticized this work, suggesting that all the candidate loci detected by our haplotype statistics, including the positive controls, are unlikely to be sweeps at all and that instead these haplotype patterns can be more easily explained by complex neutral demographic models. They also claim that these neutral non-sweeps are likely to be hard instead of soft sweeps. Here, we reanalyze the DGRP data using a range of complex admixture demographic models and reconfirm our original published results suggesting that the majority of recent and strong sweeps in D. melanogaster are first likely to be true sweeps, and second, that they do appear to be soft. Furthermore, we discuss ways to take this work forward given that most demographic models employed in such analyses are necessarily too simple to capture the full demographic complexity, while more realistic models are unlikely to be inferred correctly because they require a large number of free parameters.

Suppression gene drive in continuous space can result in unstable persistence of both drive and wild-type alleles.

Rapid evolutionary processes can produce drastically different outcomes when studied in panmictic population models vs. spatial models. One such process is gene drive, which describes the spread of "selfish" genetic elements through a population. Engineered gene drives are being considered for the suppression of disease vectors or invasive species. While laboratory experiments and modelling in panmictic populations have shown that such drives can rapidly eliminate a population, it remains unclear if these results translate to natural environments where individuals inhabit a continuous landscape. Using spatially explicit simulations, we show that the release of a suppression drive can result in what we term "chasing" dynamics, in which wild-type individuals recolonize areas where the drive has locally eliminated the population. Despite the drive subsequently reconquering these areas, complete population suppression often fails to occur or is substantially delayed. This increases the likelihood that the drive is lost or that resistance evolves. We analyse how chasing dynamics are influenced by the type of drive, its efficiency, fitness costs, and ecological factors such as the maximal growth rate of the population and levels of dispersal and inbreeding. We find that chasing is more common for lower efficiency drives when dispersal is low and that some drive mechanisms are substantially more prone to chasing behaviour than others. Our results demonstrate that the population dynamics of suppression gene drives are determined by a complex interplay of genetic and ecological factors, highlighting the need for realistic spatial modelling to predict the outcome of drive releases in natural populations.

A CRISPR homing gene drive targeting a haplolethal gene removes resistance alleles and successfully spreads through a cage population.

Engineered gene drives are being explored as a new strategy in the fight against vector-borne diseases due to their potential for rapidly spreading genetic modifications through a population. However, CRISPR-based homing gene drives proposed for this purpose have faced a major obstacle in the formation of resistance alleles that prevent Cas9 cleavage. Here, we present a homing drive in Drosophila melanogaster that reduces the prevalence of resistance alleles below detectable levels by targeting a haplolethal gene with two guide RNAs (gRNAs) while also providing a rescue allele. Resistance alleles that form by end-joining repair typically disrupt the haplolethal target gene and are thus removed from the population because individuals that carry them are nonviable. We demonstrate that our drive is highly efficient, with 91% of the progeny of drive heterozygotes inheriting the drive allele and with no functional resistance alleles observed in the remainder. In a large cage experiment, the drive allele successfully spread to all individuals within a few generations. These results show that a haplolethal homing drive can provide an effective tool for targeted genetic modification of entire populations.

The Effects of Quantitative Trait Architecture on Detection Power in Short-Term Artificial Selection Experiments.

Evolve and resequence (E&R) experiments, in which artificial selection is imposed on organisms in a controlled environment, are becoming an increasingly accessible tool for studying the genetic basis of adaptation. Previous work has assessed how different experimental design parameters affect the power to detect the quantitative trait loci (QTL) that underlie adaptive responses in such experiments, but so far there has been little exploration of how this power varies with the genetic architecture of the evolving traits. In this study, we use forward simulation to build a more realistic model of an E&R experiment in which a quantitative polygenic trait experiences a short, but strong, episode of truncation selection. We study the expected power for QTL detection in such an experiment and how this power is influenced by different aspects of trait architecture, including the number of QTL affecting the trait, their starting frequencies, effect sizes, clustering along a chromosome, dominance, and epistasis patterns. We show that all of these parameters can affect allele frequency dynamics at the QTL and linked loci in complex and often unintuitive ways, and thus influence our power to detect them. One consequence of this is that existing detection methods based on models of independent selective sweeps at individual QTL often have lower detection power than a simple measurement of allele frequency differences before and after selection. Our findings highlight the importance of taking trait architecture into account when designing and interpreting studies of molecular adaptation with temporal data. We provide a customizable modeling framework that will enable researchers to easily simulate E&R experiments with different trait architectures and parameters tuned to their specific study system, allowing for assessment of expected detection power and optimization of experimental design.