Dynamics of gene introgression in the African malaria vector Anopheles gambiae

Anopheles gambiae is a major malaria vector in Africa and a popular model species for a variety of ecological, evolutionary, and genetic studies on vector control. Genetic manipulation of mosquito vectorial capacity is a promising new weapon for the control of malaria. However, the release of exotic transgenic mosquitoes will bring in novel alleles in addition to the parasite-inhibiting genes, which may have unknown effects on the local population. Therefore, it is necessary to develop methodologies that can be used to evaluate the spread rate of introduced genes in A. gambiae. In this study, the effects and dynamics of genetic introgression between two geographically distinct A. gambiae populations from western Kenya (Mbita) and eastern Tanzania (Ifakara) were investigated with amplified fragment length polymorphisms (AFLPs) and microsatellite markers. Microsatellites and polymorphic cDNA markers revealed a large genetic differentiation between the two populations (average F(ST) = 0.093, P < 0.001). When the two strains were crossed in random mating between the two populations, significant differences in the rate of genetic introgression were found in the mixed populations. Allele frequencies of 18 AFLP markers (64.3%) for Mbita and of 26 markers (92.9%) for Ifakara varied significantly from F5 to F20. This study provides basic information on how a mosquito release program would alter the genetic makeup of natural populations, which is critical for pilot field testing and ecological risk evaluation of transgenic mosquitoes.     

Revisiting the role of introgression vs shared ancestral polymorphisms as key processes shaping genetic diversity in the recently separated sibling species of the Anopheles gambiae complex

The role of interspecific hybridisation in the evolution of pest species is poorly understood. In mosquito disease vectors this is of particular importance due to the evolution of insecticide resistance and the proposed release of transgenic strains that are refractory to the malaria parasite. In this study, we apply population genetic methods in a novel manner to determine whether mitochondrial DNA sequences have introgressed between the closely related African malaria vectors Anopheles gambiae and A. arabiensis. Our results suggest that speciation was geologically recent and ancestral haplotypes at the ND5 locus are retained in both species. In addition, comparing haplotype frequencies in allopatric and sympatric populations, suggest locale specific unidirectional introgression of mitochondria from A. arabiensis into A. gambiae.     

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

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

The toxin and antidote puzzle: new ways to control insect pest populations through manipulating inheritance

Insects carry out essential ecological functions, such as pollination, but also cause extensive damage to agricultural crops, and transmit human diseases such as malaria and dengue fever. Advances in insect transgenesis are making it increasingly feasible to engineer genes conferring desirable phenotypes, and gene drive systems are required to spread these genes into wild populations. Medea provides one solution, being able to spread into a population from very low initial frequencies through the action of a maternally-expressed toxin linked to a zygotically-expressed antidote. Several other toxin-antidote combinations are imaginable that distort the offspring ratio in favor of a desired transgene, or drive the population towards an all-male crash. We explore two such systems--Semele, which is capable of spreading a desired transgene into an isolated population in a confined manner; and Merea, which is capable of inducing a local population crash when located on the Z chromosome of a Lepidopteron pest.     

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

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

Using underdominance to bi-stably transform local populations

Underdominance refers to natural selection against individuals with a heterozygous genotype. Here, we analyze a single-locus underdominant system of two large local populations that exchange individuals at a certain migration rate. The system can be characterized by fixed points in the joint allele frequency space. We address the conditions under which underdominance can be applied to transform a local population that is receiving wildtype immigrants from another population. In a single population, underdominance has the benefit of complete removal of genetically modified alleles (reversibility) and coexistence is not stable. The two population system that exchanges migrants can result in internal stable states, where coexistence is maintained, but with additional release of wildtype individuals the system can be reversed to a fully wildtype state. This property is critically controlled by the migration rate. We approximate the critical minimum frequency required to result in a stable population transformation. We also concentrate on the destabilizing effects of fitness and migration rate asymmetry. Practical implications of our results are discussed in the context of utilizing underdominance to genetically modify wild populations. This is of importance especially for genetic pest management strategies, where locally stable and potentially reversible transformations of populations of disease vector species are of interest.     

Conditions for the persistence of partial zygotic reproductive isolation in an island-continent model

A theoretical study has been made of the contact between two populations (island-continent model) which are partially isolated insofar as they are monomorphic for different alleles at one or more strictly underdominant loci with multiplicative heterozygote fitnesses. This genetic model of an isolating mechanism is the simplest model of zygotic isolation and is also a very good representation of underdominant chromosomal rearrangements. The particular problem tackled is how the threshold value of the migration rate (m_c), below which differentiation at the underdominant loci can persist, varies as a function of the number of underdominant loci involved (n), for a given level of selection against hybrids (α). It has been found that the stability of genetic differentiation causing a given level α of reproductive isolation decreases as the number of loci involved increases. This effect is particularly noticeable for low levels of isolation. These results have some relevance with regard to the question of whether the establishment of zygotic reproductive isolation between populations occurs more probably as a result of many mutations with small effects on fitness or of only a few mutations with significant effects.     

Genetically engineered underdominance for manipulation of pest populations: a deterministic model

We theoretically investigate the potential for introgressing a desired engineered gene into a pest population by linking the desired gene to DNA constructs that exhibit underdominance properties. Our deterministic model includes two independently segregating engineered constructs that both carry a lethal gene, but suppress each other. Only genotypes containing both or neither construct are viable. Both constructs also carry the desired gene with an independent regulatory mechanism. We examine the minimal number of individuals of an engineered strain that must be released into a natural population to successfully introgress the desired gene. We compare results for strains carrying single and multiple insertions of the constructs. When there are no fitness costs associated with the inserted constructs (when the lethal sequences are not expressed), the number of individuals that must be released decreases as the number of insertions in the genome of the released strain increases. As fitness costs increase, the number of individuals that must be released increases at a greater rate for release strains with more insertions. Under specific conditions this results in the strain with only a single insertion of each construct being the most efficient for introgressing the desired gene. We discuss practical implications of our findings.     

The effect of gene drive on containment of transgenic mosquitoes

Mosquito-borne diseases such as malaria and dengue fever continue to be a major health problem through much of the world. Several new potential approaches to disease control utilize gene drive to spread anti-pathogen genes into the mosquito population. Prior to a release, these projects will require trials in outdoor cages from which transgenic mosquitoes may escape, albeit in small numbers. Most genes introduced in small numbers are very likely to be lost from the environment; however, gene drive mechanisms enhance the invasiveness of introduced genes. Consequently, introduced transgenes may be more likely to persist than ordinary genes following an accidental release. Here, we develop stochastic models to analyze the loss probabilities for several gene drive mechanisms, including homing endonuclease genes, transposable elements, Medea elements, the intracellular bacterium Wolbachia, engineered underdominance genes, and meiotic drive. We find that Medea and Wolbachia present the best compromise between invasiveness and containment for the six gene drive systems currently being considered for the control of mosquito-borne disease.     

Double drives and private alleles for localised population genetic control

Synthetic gene drive constructs could, in principle, provide the basis for highly efficient interventions to control disease vectors and other pest species. This efficiency derives in part from leveraging natural processes of dispersal and gene flow to spread the construct and its impacts from one population to another. However, sometimes (for example, with invasive species) only specific populations are in need of control, and impacts on non-target populations would be undesirable. Many gene drive designs use nucleases that recognise and cleave specific genomic sequences, and one way to restrict their spread would be to exploit sequence differences between target and non-target populations. In this paper we propose and model a series of low threshold double drive designs for population suppression, each consisting of two constructs, one imposing a reproductive load on the population and the other inserted into a differentiated locus and controlling the drive of the first. Simple deterministic, discrete-generation computer simulations are used to assess the alternative designs. We find that the simplest double drive designs are significantly more robust to pre-existing cleavage resistance at the differentiated locus than single drive designs, and that more complex designs incorporating sex ratio distortion can be more efficient still, even allowing for successful control when the differentiated locus is neutral and there is up to 50% pre-existing resistance in the target population. Similar designs can also be used for population replacement, with similar benefits. A population genomic analysis of CRISPR PAM sites in island and mainland populations of the malaria mosquito Anopheles gambiae indicates that the differentiation needed for our methods to work can exist in nature. Double drives should be considered when efficient but localised population genetic control is needed and there is some genetic differentiation between target and non-target populations.     

Introgression and persistence of cultivar alleles in wild carrot (Daucus carota) populations in the United States

Premise

Cultivated species and their wild relatives often hybridize in the wild, and the hybrids can survive and reproduce in some environments. However, it is unclear whether cultivar alleles are permanently incorporated into the wild genomes or whether they are purged by natural selection. This question is key to accurately assessing the risk of escape and spread of cultivar genes into wild populations.

Methods

We used genomic data and population genomic methods to study hybridization and introgression between cultivated and wild carrot (Daucus carota) in the United States. We used single nucleotide polymorphisms (SNPs) obtained via genotyping by sequencing for 450 wild individuals from 29 wild georeferenced populations in seven states and 144 cultivars from the United States, Europe, and Asia.

Results

Cultivated and wild carrot formed two genetically differentiated groups, and evidence of crop–wild admixture was detected in several but not all wild carrot populations in the United States. Two regions were identified where cultivar alleles were present in wild carrots: California and Nantucket Island (Massachusetts). Surprisingly, there was no evidence of introgression in some populations with a long-known history of sympatry with the crop, suggesting that post-hybridization barriers might prevent introgression in some areas.

Conclusions

Our results provide support for the introgression and long-term persistence of cultivar alleles in wild carrots populations. We thus anticipate that the release of genetically engineered (GE) cultivars would lead to the introduction and spread of GE alleles in wild carrot populations.     

The impact of dissociation on transposon-mediated disease control strategies

Vector-borne diseases such as malaria and dengue fever continue to be a major health concern through much of the world. The emergence of chloroquine-resistant strains of malaria and insecticide-resistant mosquitoes emphasize the need for novel methods of disease control. Recently, there has been much interest in the use of transposable elements to drive resistance genes into vector populations as a means of disease control. One concern that must be addressed before a release is performed is the potential loss of linkage between a transposable element and a resistance gene. Transposable elements such as P and hobo have been shown to produce internal deletion derivatives at a significant rate, and there is concern that a similar process could lead to loss of the resistance gene from the drive system following a transgenic release. Additionally, transposable elements such as Himar1 have been shown to transpose significantly more frequently when free of exogenous DNA. Here, we show that any transposon-mediated gene drive strategy must have an exceptionally low rate of dissociation if it is to be effective. Additionally, the resistance gene must confer a large selective advantage to the vector to surmount the effects of a moderate dissociation rate and transpositional handicap.     

Medusa: A Novel Gene Drive System for Confined Suppression of Insect Populations

Gene drive systems provide novel opportunities for insect population suppression by driving genes that confer a fitness cost into pest or disease vector populations; however regulatory issues arise when genes are capable of spreading across international borders. Gene drive systems displaying threshold properties provide a solution since they can be confined to local populations and eliminated through dilution with wild-types. We propose a novel, threshold-dependent gene drive system, Medusa, capable of inducing a local and reversible population crash. Medusa consists of four components - two on the X chromosome, and two on the Y chromosome. A maternally-expressed, X-linked toxin and a zygotically-expressed, Y-linked antidote results in suppression of the female population and selection for the presence of the transgene-bearing Y because only male offspring of Medusa-bearing females are protected from the effects of the toxin. At the same time, the combination of a zygotically-expressed, Y-linked toxin and a zygotically-expressed, X-linked antidote selects for the transgene-bearing X in the presence of the transgene-bearing Y. Together these chromosomes create a balanced lethal system that spreads while selecting against females when present above a certain threshold frequency. Simple population dynamic models show that an all-male release of Medusa males, carried out over six generations, is expected to induce a population crash within 12 generations for modest release sizes on the order of the wild population size. Re-invasion of non-transgenic insects into a suppressed population can result in a population rebound; however this can be prevented through regular releases of modest numbers of Medusa males. Finally, we outline how Medusa could be engineered with currently available molecular tools.     

First Steps towards Underdominant Genetic Transformation of Insect Populations

The idea of introducing genetic modifications into wild populations of insects to stop them from spreading diseases is more than 40 years old. Synthetic disease refractory genes have been successfully generated for mosquito vectors of dengue fever and human malaria. Equally important is the development of population transformation systems to drive and maintain disease refractory genes at high frequency in populations. We demonstrate an underdominant population transformation system in Drosophila melanogaster that has the property of being both spatially self-limiting and reversible to the original genetic state. Both population transformation and its reversal can be largely achieved within as few as 5 generations. The described genetic construct {Ud} is composed of two genes; (1) a UAS-RpL14.dsRNA targeting RNAi to a haploinsufficient gene RpL14 and (2) an RNAi insensitive RpL14 rescue. In this proof-of-principle system the UAS-RpL14.dsRNA knock-down gene is placed under the control of an Actin5c-GAL4 driver located on a different chromosome to the {Ud} insert. This configuration would not be effective in wild populations without incorporating the Actin5c-GAL4 driver as part of the {Ud} construct (or replacing the UAS promoter with an appropriate direct promoter). It is however anticipated that the approach that underlies this underdominant system could potentially be applied to a number of species.     

High-throughput sorting of mosquito larvae for laboratory studies and for future vector control interventions

ABSTRACT: BACKGROUND: Mosquito transgenesis offers new promises for the genetic control of vector-borne infectious diseases such as malaria and dengue fever. Genetic control strategies require the release of large number of male mosquitoes into field populations, whether they are based on the use of sterile males (sterile insect technique, SIT) or on introducing genetic traits conferring refractoriness to disease transmission (population replacement). However, the current absence of high-throughput techniques for sorting different mosquito populations impairs the application of these control measures. METHODS: A method was developed to generate large mosquito populations of the desired sex and genotype. This method combines flow cytometry and the use of Anopheles gambiae transgenic lines that differentially express fluorescent markers in males and females. RESULTS: Fluorescence-assisted sorting allowed single-step isolation of homozygous transgenic mosquitoes from a mixed population. This method was also used to select wild-type males only with high efficiency and accuracy, a highly desirable tool for genetic control strategies where the release of transgenic individuals may be problematic. Importantly, sorted males showed normal mating ability compared to their unsorted brothers. CONCLUSIONS: The developed method will greatly facilitate both laboratory studies of mosquito vectorial capacity requiring high-throughput approaches and future field interventions in the fight against infectious disease vectors.     

Marker-assisted introgression using non-unique marker alleles II: selection on probability of presence of the introgressed allele

If marker alleles that identify a gene for introgression are not completely unique to the different base populations, the trait allele can be lost quickly during the process of backcrossing. This study considers ways to deal with incompletely informative markers in order to retain the desired allele. Selection was based on the probability of the presence of the desired (introgressed) trait allele, which was calculated for each marker genotype, using a single marker or a diallelic or triallelic marker bracket. The percentage of individuals retaining the introgressed allele was calculated over five generations of backcrossing, for selected fractions between 0 and 1, for marker alleles that could occur in both base populations. The best results were obtained with a rather large selected fraction, when all individuals, heterozygous and homozygous for the most desirable allele at the marker loci, were selected. Additional selection against marker homozygotes (which might have the highest probability of carrying the desired-trait allele, but produce uninformative gametes) altered the optimum selected fraction, making the selected fraction more consistently inversely related to a better retention of the desired-trait allele. A marker bracket was found to give a better retention of the desired-trait allele than a single marker and triallelic markers were better than diallelic markers, giving a retention of almost 50%. The earlier that preselection of parents (on informativeness) took place the better the overall result; preselection should occur preferably in the base populations. Preselection could make marker alleles unique to alternative base populations and markers would effectively become fully informative. Selection in the base populations might not be possible or not desirable, for example, because of the available number of individuals. This is unlikely to be a problem when parents are paired up to exclude any common marker alleles.     

Balancing selection,genetic drift,and human‐mediated introgression interplay to shape MHC (functional) diversity in Mediterranean brown trout

The extraordinary polymorphism of major histocompatibility complex (MHC) genes is considered a paradigm of pathogen‐mediated balancing selection, although empirical evidence is still scarce. Furthermore, the relative contribution of balancing selection to shape MHC population structure and diversity, compared to that of neutral forces, as well as its interaction with other evolutionary processes such as hybridization, remains largely unclear. To investigate these issues, we analyzed adaptive (MHC‐DAB gene) and neutral (11 microsatellite loci) variation in 156 brown trout (Salmo trutta complex) from six wild populations in central Italy exposed to introgression from domestic hatchery lineages (assessed with the LDH gene). MHC diversity and structuring correlated with those at microsatellites, indicating the substantial role of neutral forces. However, individuals carrying locally rare MHC alleles/supertypes were in better body condition (a proxy of individual fitness/parasite load) regardless of the zygosity status and degree of sequence dissimilarity of MHC, hence supporting balancing selection under rare allele advantage, but not heterozygote advantage or divergent allele advantage. The association between specific MHC supertypes and body condition confirmed in part this finding. Across populations, MHC allelic richness increased with increasing admixture between native and domestic lineages, indicating introgression as a source of MHC variation. Furthermore, introgression across populations appeared more pronounced for MHC than microsatellites, possibly because initially rare MHC variants are expected to introgress more readily under rare allele advantage. Providing evidence for the complex interplay among neutral evolutionary forces, balancing selection, and human‐mediated introgression in shaping the pattern of MHC (functional) variation, our findings contribute to a deeper understanding of the evolution of MHC genes in wild populations exposed to anthropogenic disturbance.     

Effects of stocking on the genetic diversity of brown trout populations of the Adriatic and Danubian drainages in Switzerland

This study focuses on genetic variation of brown trout Salmo trutta populations of the Adriatic and Danubian drainages in Switzerland. The allozyme and other protein loci data show a major replacement of native stocks from the Adriatic drainages by introduced hatchery trout of Atlantic basin origin. In most samples, diagnostic alleles for the Adriatic form of Salmo trutta f. fario and for the marbled trout Salmo trutta marmoratus are found at very low frequencies (f<0.15). Taking into account previous genetic studies on brown trout of this basin, the Danubian samples are not heavily contaminated with foreign alleles. The results are consistent with records of local stocking activities which account in part for the high introgression rates of Atlantic alleles into local populations of the Adriatic drainages. In addition, introgression is enhanced by a decrease of natural reproduction which is caused by a deterioration of trout habitats through human activities. Furthermore, a third mechanism is proposed that may contribute to the high introgression rates observed: if Atlantic trout are introduced, the reproductive barriers between the two native forms, marbled trout and Adriatic fario respectively, break down. Atlantic trout apparently hybridize with both native forms and generate gene flow between them. In some parts of Adriatic drainages in Switzerland, the patterns of introgression and hybridization are further complicated by introduction of trout from the Danubian system. Alleles of the marbled trout are also found in the samples of the Danubian drainage system. These are due to stocking activities across the watershed.     

Gene drive systems in mosquitoes: rules of the road

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

General principles of single-construct chromosomal gene drive

Gene drive systems are genetic elements capable of spreading into a population even if they confer a fitness cost to their host. We consider a class of drive systems consisting of a chromosomally located, linked cluster of genes, the presence of which renders specific classes of offspring arising from specific parental crosses unviable. Under permissive conditions, a number of these elements are capable of distorting the offspring ratio in their favor. We use a population genetic framework to derive conditions under which these elements spread to fixation in a population or induce a population crash. Many of these systems can be engineered using combinations of toxin and antidote genes, analogous to Medea, which consists of a maternal toxin and zygotic antidote. The majority of toxin-antidote drive systems require a critical frequency to be exceeded before they spread into a population. Of particular interest, a Z-linked Medea construct with a recessive antidote is expected to induce an all-male population crash for release frequencies above 50%. We suggest molecular tools that may be used to build these systems, and discuss their relevance to the control of a variety of insect pest species, including mosquito vectors of diseases such as malaria and dengue fever.