Skip to main content

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

Abstract

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.

Background

Powerful homing gene drives can spread from low starting frequencies throughout an entire population [1,2,3,4,5]. However, this capability also renders such drives highly invasive since small levels of migration could facilitate their spread into any connected populations [6, 7]. Though potentially desirable in some applications such as global modification or elimination of disease vectors, this could also substantially increase the social and political difficulties associated with deploying a gene drive in the field due to fears of uncontrollable spread. A gene drive that suppresses invasive species or agricultural pests would also likely raise concern if it would affect these species in their native range.

Gene drive technology has improved markedly over the past few years, and several different CRISPR-based gene drives have now been demonstrated in yeast [8,9,10,11], flies [12,13,14,15,16,17,18,19,20,21,22,23,24,25,26], mosquitoes [27,28,29,30,31,32,33,34,35,36], and mice [37]. Most of these are homing drives, which have successfully achieved the modification and suppression of laboratory populations [26, 30, 36]. Several promising applications have been proposed for this technology, such as the genetic modification of Aedes and Anopheles mosquito populations by introducing genes that could reduce transmission of malaria, dengue, and other diseases [1,2,3,4,5]. Furthermore, gene drives could be used to directly suppress populations of disease vectors or invasive species [1,2,3,4,5]. Overcoming the challenge of confining a gene drive to a target population could represent an important step in bringing such approaches closer to deployment.

Several types of potentially confinable gene drives have been developed [1, 38, 39]. One class of especially promising candidates are CRISPR toxin-antidote (TA) gene drives. These drives will only spread in the population when introduced above a threshold frequency, which is determined by the parameters of the drive and is usually above zero if the drive has any imperfections, with some forms having nonzero introduction frequencies even in idealized form [40, 41]. Below the threshold, the drive frequency will tend to decline. If migrants to connected populations cannot propel the drive above this frequency, TA drives will not spread in these populations and can thereby remain confined [40, 42]. While the “migration” threshold (assuming new migrants come in each generation) will always be lower than the “introduction” threshold, the existence of one implies the existence of the other [40,41,42]. Recent studies have already used TA drives to successfully modify cage populations [42, 43]. Despite spreading more slowly than effective homing drives and also, by design, requiring larger release sizes, such confined drives could ease concerns associated with less controllable homing drives.

The development of a confined suppression drive poses a greater challenge than confined modification drives due to the need for the drives to retain enough power to spread through the population and still provide the high genetic load required for effective population suppression. A spatially confined gene drive capable of suppression, or perhaps similarly spreading a cargo gene with a high fitness cost, would therefore be of considerable interest. Three novel systems have been suggested to this end: Toxin-Antidote Dominant Embryo (TADE) suppression drives [40, 41], daisy homing drives [44], and tethered homing drive systems [45]. TADE suppression drives disrupt essential fertility genes and additionally target a haplolethal gene with Cas9, while the drive allele provides a rescue copy of the haplolethal gene [40, 41]. These drives are capable of confined population suppression, but they require a haplolethal target gene, which can make them difficult to engineer. Many transformed individuals can be lost after embryo microinjection, since many cells would possess disrupted copies of the haplolethal gene [26]. Daisy-chain homing drives consist of a series of homing elements, with each element driving the next in the chain [44]. As the system spreads through a population, each element is lost in turn until the final element is no longer able to increase in frequency [44]. Daisy-chain drives require multiple suitable target genes to avoid resistance, and because many driving elements are required, construction is complicated [46]. In addition, with higher migration rates and lower cargo fitness costs, daisy-chain drives could be more difficult to confine [46] and unable to keep intermediate levels of suppression for long if they fail to rapidly eliminate the population.

Tethered homing drive systems provide a potential alternative to solve some of these issues. A tethered homing drive system (Fig. 1) involves the modification of a population with a confined drive, along with the spread of a homing drive that is only capable of drive activity in the presence of the confined drive [45]. Thus, a tethered system potentially combines benefits of confinement with the power of a homing drive (Additional file 1: Table S1). Modeling suggests that such a system could be capable of confined suppression of a target population [45]. A tethered system for modification should also be able to spread cargo genes with much higher fitness costs than other confined drives.

Fig. 1
figure 1

Standard and tethered homing drive comparison. Populations are shown at three different times after drive release. a A standard homing drive has no threshold frequency, and it can spread to a wild-type population after a small release. b A tethered homing drive cannot increase in frequency within a wild-type population since it lacks Cas9 for drive conversion. c A TARE drive can spread after a moderate release, followed by a tethered homing drive after a small release. d A sufficiently small introduction of TARE drive to a population will not allow the TARE drive to spread, which will in turn prevent spread of the tethered homing drive

Here, we show that a tethered gene drive system can modify cage populations of Drosophila melanogaster. We use a regionally confinable Toxin-Antidote Recessive Embryo (TARE) drive targeting an essential but haplosufficient gene. This drive already contains the necessary Cas9 gene for a tethered homing drive, making it particularly suitable for use in a tethered system. We then test two split homing drives in the TARE-modified populations: The first is a modification drive targeting a haplolethal gene with two gRNAs. The second is a suppression drive targeting a haplosufficient but essential female fertility gene with four gRNAs. The TARE drive spread successfully when released well above its introduction threshold frequency and then was able to provide a Cas9 source for the homing drives, which increased in frequency according to their expected behavior.

Results

Drive construct design

The tethered drive concept involves the release of a confined drive and an incomplete homing drive, lacking a component provided by the confined system and thus essentially confining them to the same target area. Our tethered gene drive system involves a TARE drive carrying Cas9 together with two homing drives, one to demonstrate a population modification system and one to demonstrate a population suppression system. The homing drives, both constructed and tested in previous studies [26, 47], lack Cas9 and are incapable of drive conversion in wild-type individuals. Since the tethered system is modular, any engineered split homing drives will be compatible with the TARE drive, and any other confined system containing Cas9 can provide support for the homing drives.

The TARE drive is located in the first exon of the hairy locus (h) on chromosome 3L (Fig. 2A, B) and is similar to the TARE drive used to modify a cage population in a previous study [42]. It includes a recoded h sequence, Cas9 expressed by the nanos promoter, EGFP as a marker expressed by a 3xP3 promoter and SV40 3’ UTR, and two gRNAs expressed by the U6:3 promoter. Downstream from the TARE drive, the third exon of h is targeted and disrupted by the drive’s gRNA. This prevents copying of the whole drive by homology-directed repair. The drive works by disrupting wild-type alleles, creating recessive lethal alleles that are then removed from the population (Additional file 1: Figs. S1A, S2)

Fig. 2
figure 2

Tethered drive constructs. a The TARE drive includes a recoded h sequence to provide rescue for h, Cas9 with the nanos promoter/5′UTR and 3′UTR/terminator, EGFP driven by the 3xP3 promoter with a SV40 UTR, and two gRNAs expressed by the U6:3 promoter with tRNAs to separate mature gRNAs. b The TARE drive is inserted into the first exon of the wild-type h allele, located on chromosome 3L, and the drive targets the third coding exon with its gRNAs. c The haplolethal homing drive similarly includes a recoded RpL35A sequence, DsRed, and two gRNAs. d The haplolethal homing drive targets the second exon of the wild-type RpL35A allele, located on chromosome 3R. e The homing suppression drive similarly includes DsRed and four gRNAs. f The homing suppression drive targets the second exon of the wild-type yellow-g allele, located on chromosome 3L

Our haplolethal homing drive was constructed previously [26]. In short, it is a population modification homing drive that targets the haplolethal gene RpL35A with two gRNAs and provides a recoded RpL35A as rescue (Fig. 2C, D). Since the target gene is haplolethal, individuals with disrupted alleles will not be viable. Our homing suppression drive is described in detail in another study [47]. It targets a haplosufficient gene with four gRNAs (Fig. 2E, F). Both homing drives are split drives and do not contain Cas9. They are only capable of drive activity in TARE-modified individuals, which express Cas9 in the germline (Additional file 1: Fig. S1B-C).

These homing drives have been engineered to reduce the incidence and impact of resistance alleles. Because the drives use multiplexed gRNAs, the formation of functional resistance alleles will be extremely rare [25]. Nonfunctional resistance alleles will lead to nonviable individuals in the haplolethal drive and will not outcompete drive alleles in the suppression drive, although the formation of large quantities of nonfunctional resistance alleles can slow the rate of spread of either drive.

TARE drive evaluation

We crossed flies heterozygous for the TARE drive to w1118 flies to determine drive inheritance. Observation of the EGFP phenotype in the eyes of offspring from this cross was used to identify the presence of drive alleles (Fig. 3, Additional file 2: Data Set S1). Progeny of drive females showed 72% inheritance, a rate significantly different than Mendelian inheritance (p < 0.0001, Fisher’s exact test). Progeny of drive males showed 52% inheritance, which is not significantly different than Mendelian inheritance (p = 0.1511, Fisher’s exact test). These results were expected. Germline Cas9 activity occurs in both male and female drive carriers, creating disrupted h alleles, but this alone does not affect inheritance, accounting for the Mendelian rate in males. However, after fertilization, further Cas9 activity occurs in embryos from female drive carriers due to maternal deposition of Cas9 and gRNAs, even among individuals that did not receive a drive allele [40]. Any individuals with two disrupted alleles are nonviable, increasing the relative frequency of drive alleles in the progeny of drive-carrying females. The disrupted h alleles introduced by germline Cas9 activity in drive-carrying males remain useful when considering future generations, since they can eventually meet other disrupted alleles and then be removed from the population, further increasing the drive frequency (Additional file 1: Fig. S1A, S2).

Fig. 3
figure 3

TARE drive inheritance. Individuals heterozygous for the TARE drive were crossed with w1118 individuals, and EGFP expression in progeny indicated the presence of the drive. The size of the dots is proportional to the number of adult progeny from a single drive individual. The rate estimates and standard error of the mean (SEM) are indicated. The increased inheritance in the progeny of females is due to elimination of progeny carrying no drive alleles due to maternal Cas9/gRNA deposition causing disruption in wild-type alleles inherited from the male parent. Inheritance is Mendelian in males due to lack of maternal deposition, even though h alleles are still disrupted in the male germline

To determine germline and embryo cleavage rates, we crossed male and female TARE drive heterozygotes (Additional file 2: Data Set S1). The drive inheritance rate for this cross can be used with the rate for the cross between TARE females and w1118 males to estimate the rates of embryo activity and germline activity. Using this approach, we calculated an embryo cleavage rate of 63.2% and a germline cleavage rate of 88.8% for the TARE drive, assuming similar germline cut rates in males and females and embryo cutting only in the progeny of females (see Additional file 2: Data Set S1 for calculation). The germline cut rate is somewhat lower than observed in other studies of gene drive in D. melanogaster [13, 18, 20, 21, 23,24,25,26, 42, 43] but is still sufficient to support effective population modification based on previous modeling [40].

To further characterize and quantify the mechanism of our TARE drive, we crossed drive/wild-type heterozygotes with w1118 individuals. Flies were allowed to lay eggs for 3 days, and then transferred to new vials once per day for the next 2 days. In the second and third vials, eggs were counted, and pupae were counted in addition to phenotyping eclosed adults. The progeny of female drive heterozygotes had 80.7% egg-to-pupae survival, and the progeny of male drive heterozygotes had 88.2% egg-to-pupae survival (Fig. 4, Additional file 2: Data Set S1). The difference between these egg-to-pupae survival rates is significant (p = 0.0013, Fisher’s exact test), though slightly lower than expected (drive inheritance among progeny of females with egg counts was also slightly lower). In crosses between wild-type and drive-carrying individuals, the TARE drive uses embryo Cas9 activity to create nonviable genotypes in the progeny of females. The difference between egg-to-pupa survival rates indicates embryo Cas9 activity in the progeny of drive females, as expected.

Fig. 4
figure 4

Egg-to-pupae viability. Individuals heterozygous for the TARE drive were crossed with w1118 individuals, and eggs and pupae were counted as well as eclosed adults. The size of the dots is proportional to the number of eggs from a single female. The rate estimates and standard error of the mean (SEM) are indicated. As expected, survival was reduced in egg females that inherited two disrupted h alleles

Tethered homing drive evaluation

We crossed flies doubly heterozygous for the TARE drive and one of the homing drives with w1118 individuals to determine drive inheritance for each homing drive. Observation of the EGFP phenotype in the eyes of offspring identified the presence of the TARE drive, and the DsRed phenotype identified either the homing suppression drive (Fig. 5, Additional file 2: Data Set S2) or the haplolethal homing drive (Fig. 6, Additional file 2: Data Set S3). In both cases, inheritance rates were lower than the rates observed in the original studies using these drive lines [26, 47]. This is likely because the nanos-Cas9 source (part of the TARE drive) led to lower levels of Cas9 expression than the previously tested source. In addition, while the TARE drive and haplolethal homing drive are located at two distant loci, they are still on the same chromosome, and D. melanogaster chromosomes do not undergo crossovers in males during meiosis. This could result in suppression of the TARE drive by nonfunctional resistance allele formation from the haplolethal homing drive in this experimental cross, but the effect should be low since males with this drive have a low rate of germline resistance allele formation [26]. Egg-to-pupae survival was also determined for crosses between TARE heterozygotes and haplolethal homing drive heterozygotes (Additional file 1: Fig. S3). It was somewhat lower for progeny of TARE/homing drive females, as expected due to additional embryo nonviability resulting from the disruption of the haplolethal target gene, but it remained high for the progeny of TARE/homing drive males due to low resistance allele formation.

Fig. 5
figure 5

Homing suppression drive inheritance. Individuals heterozygous for both the TARE drive and the homing suppression drive were crossed with w1118 individuals. EGFP indicated the presence of the TARE drive, and DsRed indicated the presence of the homing drive. The size of the dots is proportional to the number of adult progeny from a single drive individual. Rate estimates and SEM are indicated. The TARE drive is inherited as an increased rate only in females, but the homing drives can copy themselves, resulting in an increased inheritance rate in both sexes

Fig. 6
figure 6

Haplolethal homing drive inheritance. Individuals heterozygous for both the TARE drive and the haplolethal homing drive were crossed with w1118 individuals. EGFP indicated the presence of the TARE drive, and DsRed indicated the presence of the homing drive. The size of the dots is proportional to the number of adult progeny from a single drive individual. Rate estimates and SEM are indicated. The TARE drive is inherited as an increased rate only in females, but the homing drives can copy themselves, resulting in an increased inheritance rate in both sexes

Tethered system cage study

To assess the performance of the tethered homing drive system in large cage populations, we allowed individuals homozygous for the TARE drive to lay eggs in bottles for one day. Similarly, w1118 individuals were allowed to lay eggs in another set of bottles for 1 day. Flies were removed, and the bottles were placed in four population cages. Emerging adults were considered to be generation zero, and the TARE drive was allowed to spread in each population, with all adults phenotyped for EGFP (Fig. 7, Additional file 2: Data Set S4). In all cages, the population fluctuated, averaging 3491 individuals in each cage (Additional file 1: Fig. S4).

Fig. 7
figure 7

Tethered homing system cage study. Individuals homozygous for the TARE drive were introduced into a population that was wild type at the drive site. The cage populations were followed for several discrete generations, each lasting 12 days, including 1 day of egg-laying. After the TARE drive reached fixation in two cages, individuals that were homozygous for the TARE drive and heterozygous for one of the split homing drives were released. All individuals from each generation were phenotyped for DsRed and EGFP, with positive drive carriers having either one or two drive alleles

In one cage with a low initial release frequency, the TARE drive increased in frequency only slowly, and in another cage with a lower release frequency, the drive frequency remained constant. These two cages experienced fluctuations in frequency (possibly due to differences in the health of initially released TARE and w1118 individuals) before quickly reaching drive carrier frequencies of 20% and 10% (the former of which then proceeded to slowly increase, while the latter remained stable), suggesting that the introduction threshold of the TARE drive may be approximately 5–15%, a value that depends on the efficiency and fitness cost of the drive [40, 42]. With the cleavage parameters calculated based on drive inheritance (Additional file 2: Data Set S1), the fitness cost of the drive in homozygotes was perhaps 10–15% based on previous models [40]. A maximum-likelihood based analysis [48] indicates that drive homozygotes had approximately a 13–15% fitness cost (fitness was 0.867 compared to wild-type fitness of 1, with a 95% confidence interval of 0.785–0.954) if costs impacted both sexes. The fitness cost was 25–27% if only one sex was impacted (with a 95% confidence interval of 0.622 - 0.906 for female fecundity alone), and models with a single fitness parameter for female fecundity, male mating success, or offspring viability all gave similar performance (Additional file 1: Table S2). All these models would result in an approximately 14% introduction frequency threshold with a 95% confidence interval of 3–21%. Based on the behavior of the two cages that were released close to the introduction threshold and where frequency trajectory should therefore be particularly sensitive to fitness effects, we favor the lower half of the fitness/introduction threshold range. The estimated effective population size was 3.6% of the census population size (with a 95% confidence interval of 2.4–5.3%), somewhat lower than similar cages analyzed by the same maximum likelihood method [26, 47,48,49]. This is perhaps indicative of a model that does not capture all relevant factors influencing genotype frequencies [48], and indeed, the rate of increase of the TARE drive was somewhat higher than predicted when it was at an intermediate frequency (40–70% carrier frequency). One possible explanation for this is underestimation of drive cut rates. Another is patchy distribution of eggs from drive and wild-type individuals, either between or within the eight food bottles of the cage experiment, thus creating areas with higher proportions of eggs from drive-carrying mothers. In these areas, some offspring from drive mothers would be nonviable, allowing their adjacent siblings access to more resources at the early larval stage than individuals in more crowded areas with a higher proportion of viable eggs.

In the two cages where the TARE drive reached 100% frequency, individuals homozygous for the TARE drive and heterozygous for one of the homing drives were released. The haplolethal homing drive was introduced at 6% frequency, and the homing suppression drive was introduced at 33% frequency (Fig. 7, Additional file 2: Data Set S4). The homing drives were allowed to spread for several generations, with all adults phenotyped for DsRed. The haplolethal homing drive eventually reached a frequency of 100%, while the homing suppression drive approached an apparent equilibrium frequency of somewhat over 50%. This latter result was due to low drive efficiency and a fitness cost, which prevented complete suppression and was seen in another study using the same drive [47]. The rate of increase of these homing drives was lower than in other studies that used these lines [26, 47], as expected from the somewhat lower drive conversion rate of these drives when combined with the TARE drive compared to the previously used nanos-Cas9 supporting element. The additional gRNAs in the TARE drive may also have somewhat reduced the cleavage activity at each site in the split homing drives due to saturation of Cas9, though this reduction was likely small [25].

Discussion

In this study, we performed an experimental trial of the “tethered gene drive” system. This system was designed to allow a potentially confined drive system, such as a TARE drive or an underdominance system, to support a homing drive, thus allowing costly alleles to be spread with the powerful homing drive, while also preventing it from spreading to populations beyond the supporting confined drive [45]. Our demonstration with TARE followed by a tethered population modification drive [26] was successful, though the tethered homing suppression drive did not perform as well. However, the equilibrium frequency of the suppression drive is determined by its fitness costs, drive conversion efficiency, and other characteristics of the drive itself [47], rather than by any feature of the tethered system, which would only be expected to slightly reduce drive performance due to saturation of Cas9 with additional gRNAs in the combined system.

Since the TARE drive utilized the same rescue element, gRNA promoter, and target sites as a previous split TARE system [42], its reduced performance as indicated by reduced drive inheritance was almost certainly due to the Cas9 element. Though both of these TARE drives had identical nanos-Cas9 elements (with the same orientation with respect to the 3xP3 promoter of the fluorescent element), their genomic locations were substantially different. In this study, the Cas9 gene was located within the TARE drive on chromosome 3L, rather than a separate genomic site on chromosome 2R. This led to reduced cleavage rates, particularly in the early embryo, a critical parameter that can influence drive inheritance in individual crosses. In this case, such reduced performance could potentially be useful to make the drive more confined to a target region. However, for the homing suppression element, high performance would likely be needed to achieve high enough drive conversion efficiency and genetic load for population suppression [25], particularly in complex natural environments [50, 51]. Higher expression rates of Cas9 within the TARE element could possibly be achieved by reorienting the Cas9 element within the drive allele. It would also be straightforward to choose a different site within the h gene, or more likely, a different target gene entirely. Essential but haplosufficient genes can likely be found in most genomes, and similar CRISPR toxin-antidote ClvR elements have already targeted other such genes with success [43, 52].

In general, the tethered system provides a highly promising strategy for confining a homing drive to a target population that is sufficiently isolated from other non-target populations. If an efficient homing drive can be made, then a suitable TARE system can likely be engineered as well. This is because a suitable homing system would have a high drive conversion rate, implying at least an equal rate of germline cleavage (which includes both the drive conversion rate and germline resistance allele formation rate for homing drives). A TARE system would be expected to still have good efficiency without embryo cutting if the germline cut rate was high (without embryo cutting, it would take only a few more generations to spread to the whole population [40]), allowing it to utilize the same promoter. While TARE would additionally need a rescue element, several studies have already shown that their engineering is feasible in flies [26, 42, 43, 52] and mosquitoes [36]. Furthermore, a TARE drive could be engineered to have an additional CRISPR nuclease with a different promoter and different gRNA specificity [53,54,55], allowing germline-restricted CRISPR cleavage for the homing element and both germline and embryo cleavage for the TARE element. This could increase the efficiency of the homing drive, which generally suffers from embryo cutting. However, modeling indicates that a homing drive with high drive conversion efficiency or lacking significant fitness costs would still be able to tolerate high embryo cut rates, as long as functional resistance alleles did not form [25, 56, 57].

For more stringent confinement than what can be provided by a standard TARE drive, other CRISPR toxin-antidote systems with higher introduction thresholds could likely be engineered with similar ease, such as 2-locus TARE systems (Additional file 1: Table S1) [41]. This additional confinement could provide a buffer against an evolved reduction in fitness costs or high variation in migration levels that could otherwise occasionally result in a drive exceeding its introduction threshold in a non-target population. However, in situations where an efficient homing drive could not be engineered due to difficulty in achieving a high enough rate of drive conversion by homology-directed repair after germline cleavage, then TADE suppression systems [40, 41] may present a useful alternative option for confined suppression if targeting of haplolethal genes could be successfully engineered.

Conclusions

This study experimentally demonstrated the feasibility of using tethered homing gene drives in cage populations based on CRISPR toxin-antidote drive systems. Such systems can allow for the power of homing drives to be limited to target populations. This could be useful and even necessary when dealing with invasive species outside their native range, when sociopolitical considerations demand a limited drive, or even during an initial testing phase of a homing gene drive system.

Methods

Plasmid construction

The starting plasmid EGDh2 was constructed previously [42] (see Supplementary Information). Restriction enzymes for plasmid digestion, Q5 Hot Start DNA Polymerase for PCR, and Assembly Master Mix for Gibson assembly were acquired from New England Biolabs. Oligonucleotides were obtained from Integrated DNA Technologies. 5-α competent Escherichia coli from New England Biolabs and ZymoPure Midiprep kit from Zymo Research were used to transform and purify plasmids. See Additional file 1 for a list of DNA fragments, plasmids, primers, and restriction enzymes used for cloning.

Generation of transgenic lines

Injections were conducted by Rainbow Transgenic Flies. The donor plasmid TAREhNU2G (507 ng/μL) was injected along with plasmid EGDhg2t (149 ng/μL), which provided gRNAs for transformation and was constructed previously [42], and pBS-Hsp70-Cas9 (442 ng/μL, Addgene plasmid #45945), which provided Cas9. A 10 mM Tris-HCl, 100 μM EDTA solution at pH 8.5 was used for the injection. Both lines containing split homing modification [26] and suppression [47] drives were generated in previous studies.

Genotypes and phenotypes

TARE drive carriers are indicated by expression of EGFP driven by the 3xP3 promoter, which is highly visible in the eyes of w1118 flies. Flies carrying either of the homing drives are similarly marked by 3xP3-DsRed, which can be easily distinguished from EGFP. For phenotyping, flies were frozen, and scored for red and green fluorescence in the eyes using the NIGHTSEA system.

Fly rearing

Flies were reared at 25 °C with a 14/10 h day/night cycle. Bloomington Standard medium was provided as food. The cage study used 30 × 30 × 30 cm (Bugdorm, BD43030D) enclosures. Flies homozygous for the TARE drive were allowed to lay eggs in one or two food bottles for one day, and a proportionately higher number of w1118 individuals were separately allowed to lay eggs in six or seven food bottles for one day. All adults were removed, and the eight bottles were placed in one cage. Eleven days later, the old bottles were replaced with fresh food bottles, and the adult flies were left in the cage and allowed to lay eggs for 1 day before being removed and frozen for phenotyping. Adults emerging from the original bottles were considered to be generation zero. This 12-day cycle was repeated until the TARE drive had nearly reached fixation. Flies heterozygous for each homing drive and homozygous for the TARE drive were separately allowed to lay eggs in food bottles for one day at the same time as the cage flies were laying their eggs. The adults were removed before the bottles were placed into the cages with TARE drive bottles so that they were under the same conditions and would hatch at approximately the same time. The same 12-day cycle was repeated, with all flies phenotyped shortly after being frozen.

To prevent accidental releases of gene drive flies, all live drive-carrying flies were quarantined at the Sarkaria Arthropod Research Laboratory at Cornell University under Arthropod Containment Level 2 protocols in accordance with USDA APHIS standards. In addition, the homing drives contained no Cas9 and were incapable of drive conversion in wild-type flies. All safety standards were approved by the Cornell University Institutional Biosafety Committee.

Phenotype data analysis

To calculate drive parameters, we pooled all offspring from the same type of cross together and used the combined counts to calculate rates. Because offspring had different parents and were raised in separate vials, this pooling approach could potentially distort rate and error estimates. To account for such batch effects, we performed an alternate analysis similar to the one used in previous studies [25, 26, 42]. This involved fitting a generalized linear mixed-effects model with a binomial distribution (fit by maximum likelihood, Adaptive Gauss-Hermite Quadrature, nAGQ = 25). In this model, offspring from a single vial were considered as a separate batch, even if they had the same parents as offspring from other vials. This approach allows for variance between batches and results in slight changes to rate and error estimates. The analysis was performed using the R statistical computing environment (3.6.1) with the packages lme4 (1.1-21, https://cran.r-project.org/web/packages/lme4/index.html) and emmeans (1.4.2, https://cran.r-project.org/web/packages/emmeans/index.html). The specific R script we used for this analysis is available on Github (https://github.com/MesserLab/Binomial-Analysis). The resulting rate estimates and errors were similar to the pooled analysis and are provided in Additional file 2: Data Sets S1-S3.

Genotyping

To extract DNA, flies were frozen and ground in 30 μL of 10 mM Tris-HCl pH 8, 1 mM EDTA, 25 mM NaCl, and 200 μg/mL recombinant proteinase K (Thermo Scientific). The mixture was incubated at 37 °C for 30 min, then at 95 °C for 5 min. The DNA was used as a template for PCR using Q5 Hot Start DNA Polymerase (New England Biolabs), following the manufacturer’s protocol. The region of interest containing gRNA target sites was amplified using DNA oligo primers hCut_S_F and hCut_S_R (see Additional file 1 for primer sequences). DNA fragments were isolated by gel electrophoresis, Sanger sequenced, and analyzed with ApE software (http://biologylabs.utah.edu/jorgensen/wayned/ape).

Estimation of drive fitness

We used a maximum likelihood approach to estimate the fitness parameters for the gene drive lines. Details of this approach and its application to cage population data were previously described [48]. The fitness parameters for both homing lines used here have been previously estimated [26, 47]. For the TARE line, we estimated the viability, female fecundity, and male mating success parameters and the effective population size using a model that assumes a co-dominant, multiplicative fitness effect of the TARE allele (heterozygotes were assigned fitness equal to the square root of homozygotes). The genomic site of TARE (the h locus) is haplosufficient, so heterozygotes that bear one wild-type and one disrupted h allele were assumed to have the same fitness as wildtype homozygotes. The germline and embryonic rates of cleavage of the wildtype allele by TARE were set at experimentally estimated values. Parameter values were estimated by maximizing the likelihood across all generation transitions from the four cages combined (i.e., a single estimate per parameter is generated for all four cages). The first two generational transitions for each cage were discounted because of apparent parental effects temporarily reducing the fitness of drive individuals (resulting in fewer drive carriers in the next generation than expected), though this could also be partly explained by inbreeding between drive individuals eclosing in the same food bottle. The effective population size parameter was estimated as a fraction of the census population size in each generation transition, using the average of both generations that were part of the transition. The code is available on GitHub (github.com/MesserLab/TetheredDrives) [58].

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its supplementary information files. The code used to estimate the drive fitness parameters is available at github.com/MesserLab/TetheredDrives [58], and the code used to fit the generalized linear fixed-effects model is available at https://github.com/MesserLab/Binomial-Analysis.

References

  1. Champer J, Buchman A, Akbari OS. Cheating evolution: engineering gene drives to manipulate the fate of wild populations. Nat Rev Genet. 2016;17:146–59.

    Article  CAS  Google Scholar 

  2. Burt A. Heritable strategies for controlling insect vectors of disease. Philos Trans R Soc L B Biol Sci. 2014;369:20130432.

    Article  Google Scholar 

  3. Alphey L. Genetic control of mosquitoes. Annu Rev Entomol. 2014;59:205–24.

    Article  CAS  Google Scholar 

  4. Esvelt KM, Smidler AL, Catteruccia F, Church GM. Concerning RNA-guided gene drives for the alteration of wild populations. Elife. 2014:e03401. https://doi.org/10.7554/eLife.03401.

  5. Wang G-H, Du J, Chu CY, Madhav M, Hughes GL, Champer J. Symbionts and gene drive: two strategies to combat vector-borne disease. Trends Genet. 2022;0. https://doi.org/10.1016/J.TIG.2022.02.013.

  6. Unckless RL, Messer PW, Connallon T, Clark AG. Modeling the manipulation of natural populations by the mutagenic chain reaction. Genetics. 2015;201:425–31. https://doi.org/10.1534/genetics.115.177592.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Noble C, Adlam B, Church GM, Esvelt KM, Nowak MA. Current CRISPR gene drive systems are likely to be highly invasive in wild populations. Elife. 2018;7:e33423. https://doi.org/10.7554/eLife.33423.

    Article  PubMed  PubMed Central  Google Scholar 

  8. DiCarlo JE, Chavez A, Dietz SL, Esvelt KM, Church GM. Safeguarding CRISPR-Cas9 gene drives in yeast. Nat Biotechnol. 2015;33:1250–5. https://doi.org/10.1038/nbt.3412.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Roggenkamp E, Giersch RM, Schrock MN, Turnquist E, Halloran M, Finnigan GC. Tuning CRISPR-Cas9 gene drives in Saccharomyces cerevisiae. G3. 2018;8:999–1018. https://doi.org/10.1534/g3.117.300557.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Basgall EM, Goetting SC, Goeckel ME, Giersch RM, Roggenkamp E, Schrock MN, et al. Gene drive inhibition by the anti-CRISPR proteins AcrIIA2 and AcrIIA4 in Saccharomyces cerevisiae. Microbiology. 2018;164:464–74. https://doi.org/10.1099/mic.0.000635.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Shapiro RS, Chavez A, Porter CBM, Hamblin M, Kaas CS, DiCarlo JE, et al. A CRISPR–Cas9-based gene drive platform for genetic interaction analysis in Candida albicans. Nat Microbiol. 2018;3:73–82. https://doi.org/10.1038/s41564-017-0043-0.

    Article  CAS  PubMed  Google Scholar 

  12. Oberhofer G, Ivy T, Hay BA. Behavior of homing endonuclease gene drives targeting genes required for viability or female fertility with multiplexed guide RNAs. Proc Natl Acad Sci. 2018;115:E9343–52. https://doi.org/10.1073/pnas.1805278115.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Carrami EM, Eckermann KN, Ahmed HMM, Sánchez CHM, Dippel S, Marshall JM, et al. Consequences of resistance evolution in a Cas9-based sex-conversion suppression gene drive for insect pest management. Proc Natl Acad Sci. 2018:201713825. https://doi.org/10.1073/pnas.1713825115.

  14. Chae D, Lee J, Lee N, Park K, Moon SJ, Kim HH. Chemical controllable gene drive in Drosophila. ACS Synth Biol. 2020. https://doi.org/10.1021/acssynbio.0c00117.

  15. López Del Amo V, Leger BS, Cox KJ, Gill S, Bishop AL, Scanlon GD, et al. Small-molecule control of super-Mendelian inheritance in gene drives. Cell Rep. 2020;31:107841.

    Article  Google Scholar 

  16. Kandul NP, Liu J, Buchman A, Gantz VM, Bier E, Akbari OS. Assessment of a split homing based gene drive for efficient knockout of multiple genes. G3. 2020;10:827–37. https://doi.org/10.1534/g3.119.400985.

    Article  CAS  PubMed  Google Scholar 

  17. Xu X-RS, Bulger EA, Gantz VM, Klanseck C, Heimler SR, Auradkar A, et al. Active genetic neutralizing elements for halting or deleting gene drives. Mol Cell. 2020. https://doi.org/10.1016/j.molcel.2020.09.003.

  18. López Del Amo V, Bishop AL, Sánchez CHM, Bennett JB, Feng X, Marshall JM, et al. A transcomplementing gene drive provides a flexible platform for laboratory investigation and potential field deployment. Nat Commun. 2020;11:1–12. https://doi.org/10.1038/s41467-019-13977-7.

    Article  CAS  Google Scholar 

  19. Gantz VM, Bier E. Genome editing. The mutagenic chain reaction: a method for converting heterozygous to homozygous mutations. Science (80- ). 2015;348:442–4.

    Article  CAS  Google Scholar 

  20. Champer J, Reeves R, Oh SY, Liu C, Liu J, Clark AG, et al. Novel CRISPR/Cas9 gene drive constructs reveal insights into mechanisms of resistance allele formation and drive efficiency in genetically diverse populations. PLoS Genet. 2017;13:e1006796. https://doi.org/10.1101/112011.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Champer J, Liu J, Oh SY, Reeves R, Luthra A, Oakes N, et al. Reducing resistance allele formation in CRISPR gene drive. Proc Natl Acad Sci. 2018;115:5522–7. https://doi.org/10.1073/pnas.1720354115.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Champer J, Wen Z, Luthra A, Reeves R, Chung J, Liu C, et al. CRISPR gene drive efficiency and resistance rate is highly heritable with no common genetic loci of large effect. Genetics. 2019. https://doi.org/10.1534/genetics.119.302037.

  23. Champer J, Chung J, Lee YL, Liu C, Yang E, Wen Z, et al. Molecular safeguarding of CRISPR gene drive experiments. Elife. 2019;8. https://doi.org/10.7554/eLife.41439.

  24. Guichard A, Haque T, Bobik M, Xu X-RS, Klanseck C, Kushwah RBS, et al. Efficient allelic-drive in Drosophila. Nat Commun. 2019;10:1640. https://doi.org/10.1038/s41467-019-09694-w.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Champer SE, Oh SY, Liu C, Wen Z, Clark AG, Messer PW, et al. Computational and experimental performance of CRISPR homing gene drive strategies with multiplexed gRNAs. Sci Adv. 2020;6:eaaz0525. https://doi.org/10.1126/sciadv.aaz0525.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Champer J, Yang E, Lee E, Liu J, Clark AG, Messer PW. A CRISPR homing gene drive targeting a haplolethal gene removes resistance alleles and successfully spreads through a cage population. Proc Natl Acad Sci. 2020;117:24377–83. https://doi.org/10.1073/PNAS.2004373117.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Hammond AM, Kyrou K, Bruttini M, North A, Galizi R, Karlsson X, et al. The creation and selection of mutations resistant to a gene drive over multiple generations in the malaria mosquito. PLoS Genet. 2017;13:e1007039. https://doi.org/10.1371/journal.pgen.1007039.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Hammond A, Galizi R, Kyrou K, Simoni A, Siniscalchi C, Katsanos D, et al. A CRISPR-Cas9 gene drive system targeting female reproduction in the malaria mosquito vector Anopheles gambiae. Nat Biotechnol. 2015;34:78–83. https://doi.org/10.1038/nbt.3439.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Gantz VM, Jasinskiene N, Tatarenkova O, Fazekas A, Macias VM, Bier E, et al. Highly efficient Cas9-mediated gene drive for population modification of the malaria vector mosquito Anopheles stephensi. Proc Natl Acad Sci U S A. 2015;112:E6736–43. https://doi.org/10.1073/pnas.1521077112.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Kyrou K, Hammond AM, Galizi R, Kranjc N, Burt A, Beaghton AK, et al. A CRISPR-Cas9 gene drive targeting doublesex causes complete population suppression in caged Anopheles gambiae mosquitoes. Nat Biotechnol. 2018. https://doi.org/10.1038/nbt.4245.

  31. Pham TB, Phong CH, Bennett JB, Hwang K, Jasinskiene N, Parker K, et al. Experimental population modification of the malaria vector mosquito, Anopheles stephensi. PLoS Genet. 2019;15:e1008440. https://doi.org/10.1371/journal.pgen.1008440.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Li M, Yang T, Kandul NP, Bui M, Gamez S, Raban R, et al. Development of a confinable gene drive system in the human disease vector Aedes aegypti. Elife. 2020;9.

  33. Simoni A, Hammond AM, Beaghton AK, Galizi R, Taxiarchi C, Kyrou K, et al. A male-biased sex-distorter gene drive for the human malaria vector Anopheles gambiae. Nat Biotechnol. 2020:1–7. https://doi.org/10.1038/s41587-020-0508-1.

  34. Carballar-Lejarazú R, Ogaugwu C, Tushar T, Kelsey A, Pham TB, Murphy J, et al. Next-generation gene drive for population modification of the malaria vector mosquito, Anopheles gambiae. Proc Natl Acad Sci U S A. 2020;117:22805–14. https://doi.org/10.1073/pnas.2010214117.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Adolfi A, Gantz VM, Jasinskiene N, Lee HF, Hwang K, Terradas G, et al. Efficient population modification gene-drive rescue system in the malaria mosquito Anopheles stephensi. Nat Commun. 2020;11:1–13. https://doi.org/10.1038/s41467-020-19426-0.

    Article  CAS  Google Scholar 

  37. Grunwald HA, Gantz VM, Poplawski G, Xu X-RS, Bier E, Cooper KL. Super-Mendelian inheritance mediated by CRISPR–Cas9 in the female mouse germline. Nature. 2019;566:105–9. https://doi.org/10.1038/s41586-019-0875-2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Hay BA, Oberhofer G, Guo M. Engineering the composition and fate of wild populations with gene drive. Annu Rev Entomol. 2021;66:annurev-ento-020117-043154. https://doi.org/10.1146/annurev-ento-020117-043154.

    Article  CAS  Google Scholar 

  39. Leftwich PT, Edgington MP, Harvey-Samuel T, Carabajal Paladino LZ, Norman VC, Alphey L. Recent advances in threshold-dependent gene drives for mosquitoes. Biochem Soc Trans. 2018;46:1203–12. https://doi.org/10.1042/BST20180076.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Champer J, Kim IK, Champer SE, Clark AG, Messer PW. Performance analysis of novel toxin-antidote CRISPR gene drive systems. BMC Biol. 2020;18:27. https://doi.org/10.1186/s12915-020-0761-2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

  42. Champer J, Lee E, Yang E, Liu C, Clark AG, Messer PW. A toxin-antidote CRISPR gene drive system for regional population modification. Nat Commun. 2020;11:1082. https://doi.org/10.1038/s41467-020-14960-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Oberhofer G, Ivy T, Hay BA. Cleave and Rescue, a novel selfish genetic element and general strategy for gene drive. Proc Natl Acad Sci. 2019:201816928. https://doi.org/10.1073/pnas.1816928116.

  44. Noble C, Min J, Olejarz J, Buchthal J, Chavez A, Smidler AL, et al. Daisy-chain gene drives for the alteration of local populations. Proc Natl Acad Sci. 2019:201716358. https://doi.org/10.1073/pnas.1716358116.

  45. Dhole S, Lloyd AL, Gould F. Tethered homing gene drives: a new design for spatially restricted population replacement and suppression. Evol Appl. 2019:eva.12827. https://doi.org/10.1111/eva.12827.

  46. Dhole S, Vella MR, Lloyd AL, Gould F. Invasion and migration of spatially self-limiting gene drives: a comparative analysis. Evol Appl. 2018;11:794–808. https://doi.org/10.1111/eva.12583.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Yang E, Metzloff M, Langmüller AM, Clark AG, Messer PW, Champer J. A homing suppression gene drive with multiplexed gRNAs maintains high drive conversion efficiency and avoids functional resistance alleles. G3. 2022. https://doi.org/10.1093/g3journal/jkac081.

  48. Liu J, Champer J, Langmüller AM, Liu C, Chung J, Reeves R, et al. Maximum likelihood estimation of fitness components in experimental evolution. Genetics. 2019;211:1005–17. https://doi.org/10.1534/genetics.118.301893.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Langmüller AM, Champer J, Lapinska S, Xie L, Metzloff M, Liu J, et al. Fitness effects of CRISPR endonucleases in Drosophila melanogaster populations. bioRxiv. 2021. https://doi.org/10.1101/2021.05.13.444039.

  50. North AR, Burt A, Godfray HCJ. Modelling the suppression of a malaria vector using a CRISPR-Cas9 gene drive to reduce female fertility. BMC Biol. 2020;18:98. https://doi.org/10.1186/s12915-020-00834-z.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Champer J, Kim IK, Champer SE, Clark AG, Messer PW. Suppression gene drive in continuous space can result in unstable persistence of both drive and wild-type alleles. Mol Ecol. 2021;30:1086–101. https://doi.org/10.1111/mec.15788.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Oberhofer G, Ivy T, Hay BA. Gene drive and resilience through renewal with next generation Cleave and Rescue selfish genetic elements. Proc Natl Acad Sci U S A. 2020;117:9013–21. https://doi.org/10.1073/pnas.1921698117.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Chatterjee P, Jakimo N, Jacobson JM. Minimal PAM specificity of a highly similar SpCas9 ortholog. Sci Adv. 2018;4:eaau0766. https://doi.org/10.1126/sciadv.aau0766.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Chatterjee P, Jakimo N, Lee J, Amrani N, Rodríguez T, Koseki SRT, et al. An engineered ScCas9 with broad PAM range and high specificity and activity. Nat Biotechnol. 2020;38:1154–8.

  55. Nishimasu H, Shi X, Ishiguro S, Gao L, Hirano S, Okazaki S, et al. Engineered CRISPR-Cas9 nuclease with expanded targeting space. Science (80- ). 2018;361:1259–62. https://doi.org/10.1126/science.aas9129.

    Article  CAS  Google Scholar 

  56. Beaghton AK, Hammond A, Nolan T, Crisanti A, Burt A. Gene drive for population genetic control: non-functional resistance and parental effects. Proceedings Biol Sci. 2019;286:20191586. https://doi.org/10.1098/rspb.2019.1586.

    Article  CAS  Google Scholar 

  57. Liu Y, Champer J. Modeling homing suppression gene drive in haplodiploid organisms. Proc Biol Sci. 2022. https://doi.org/10.1098/rspb.2022.0320.

  58. TetheredDrives Supplemental Files. 2022. https://doi.org/10.5281/ZENODO.6422591.

Download references

Acknowledgements

We thank Jim Bull and two anonymous reviewers for their suggestions, which helped improve the quality of our manuscript.

Funding

This study was supported by the National Institutes of Health awards R21AI130635 to JC, AGC, and PWM, award F32AI138476 to JC, and award R01GM127418 to PWM.

Author information

Authors and Affiliations

Authors

Contributions

JC designed the study. MM, EY, and JC performed the experimental investigation. SD performed the computational investigation. PWM, AGC, and JC provided supervision. MM and JC drafted the manuscript. All authors revised the manuscript and approved the final version.

Corresponding author

Correspondence to Jackson Champer.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

List of plasmid construction primers and methods. Figure S1. Tethered drive mechanism. Figure S2. Relations between genotypes in the TARE drive. Figure S3. Egg-to-pupae viability for haplolethal homing drive. Figure S4. Population sizes for cage study. Table S1. Comparison of different gene drive types. Table S2. Combined maximum likelihood parameter estimates from cage populations.

Additional file 2: Data Set S1.

TARE drive crosses. Data Set S2. Haplolethal homing drive crosses. Data Set S3. Homing suppression drive crosses. Data Set S4. Cage data.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Metzloff, M., Yang, E., Dhole, S. et al. Experimental demonstration of tethered gene drive systems for confined population modification or suppression. BMC Biol 20, 119 (2022). https://doi.org/10.1186/s12915-022-01292-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12915-022-01292-5

Keywords