Homing endonuclease genes (HEGs) are ‘selfish’ genetic elements that combine the capability to selectively disrupt specific gene sequences with the ability to rapidly spread from a few individuals to an entire population through homologous recombination repair events. Because of these properties, HEGs are regarded as promising candidates to transfer genetic modifications from engineered laboratory mosquitoes to wild-type populations including Anopheles gambiae the vector of human malaria. Here we show that I-SceI and I-PpoI homing endonucleases cleave their recognition sites with high efficiency in A. gambiae cells and embryos and we demonstrate HEG-induced homologous and non-homologous repair events in a variety of functional assays. We also propose a gene drive system for mosquitoes that is based on our finding that I-PpoI cuts genomic rDNA located on the X chromosome in A. gambiae , which could be used to selectively incapacitate X-carrying spermatozoa thereby imposing a severe male-biased sex ratio.
Mosquito species of the Anopheles gambiae complex represent the major vectors of human malaria and they pose an enormous burden on global health and economies. Every year 300–500 million people are infected by malaria and over a million people die as consequence of Plasmodium parasite infections ( 1 ). While many insect pests have long been successful targeted with population control measures such as insecticides or release of sterile males ( 2 ), for others, including A. gambiae , classical control measures have largely failed to deliver long-term solutions. Disease endemic countries often do not have the economic resources and the logistics to sustain control efforts like the massive and prolonged use of insecticides. New control strategies that are affordable, easy to implement and sustainable are desperately needed.
This global health problem has prompted an unprecedented effort aimed at generating new molecular tools and a better understanding of the biology and the genetics of Anopheline mosquitoes that culminated in the sequencing of the A. gambiae genome ( 3 ) and development of gene transfer technology for a series of vectors species ( 4 , 5 ). These molecular advances have made it possible to express genes that can block the transmission of Plasmodium in model systems ( 6–8 ) or express traits facilitating the implementation of sterile insect techniques for vector control ( 9 ).
The translation of these achievements in suitable control measures still represents a major scientific and technical challenge. Genetically modified mosquitoes carrying a desired trait such as malaria refractoriness would need to be released on a gigantic scale given the vast numbers of these insects and the wide areas that are inhabited by vectors of human tropical diseases. Therefore, a mechanism must be developed to spread the desired genetic modification from a few laboratory-reared mosquitoes to a large fraction of the wild-type vector population ( 10 , 11 ). Naturally occurring ‘selfish’ genetic elements that have non-mendelian inheritance mechanisms and spread through populations even when they provide no benefit to the host organism ( 12 ) have been proposed to transform wild-type mosquito populations.
Homing endonuclease genes (HEGs) are highly specific DNA endonucleases found in some viruses, bacteria and eukaryotes. The endonuclease promotes the movement of its encoding DNA from one allele to the other by creating a double-strand break (DSB) at a specific, long (15–40 bp) target site in an allele that lacks the HEG. Homologous DNA repair then copies the HEG to the cut chromosome in a process called ‘gene conversion’ ( 13 , 14 ).
The observation that HEGs can be engineered to cleave novel DNA sequences ( 15–18 ) offers a multitude of opportunities to utilize these elements for mosquito control. For example, HEGs could be used to disrupt genes regulating the ability of Anopheles mosquitoes to function as efficient vectors for Plasmodium parasites, or to drive recombinant refractory genes through a mosquito population, rendering them unable to transmit malaria. Alternatively, HEGs designed to target an essential mosquito gene or a gene required for female fertility could be utilized to introduce a genetic load on the population leading to population size reduction or collapse ( 19 ). More recently, it has been suggested that a harmful selfish element put under the control of a promoter which is active in individuals susceptible to Plasmodium infection but inactive in refractory individuals should drive alleles causing refractoriness through the population ( 20 ). Finally, HEGs could be used to bias the sex ratio towards males, using an endonuclease that targets X-linked sequences and is expressed during male spermatogenesis from the Y chromosome ( 19 ).
To investigate the feasibility of using HEGs as driving genetic element in mosquitoes, we have analysed the activity of ectopic HEGs in both A. gambiae cells and embryos, using experimental systems that are highly predictive of in vivo behaviour of mobile genetic elements ( 21 , 22 ). We determined the activity of two well-characterized HEGs, I-PpoI (a member of the His-Cys box family of endonucleases from the slime mold Physarum polycephalum ) ( 23–25 ) and I-SceI (a LAGLIDADG class endonuclease originally isolated from Saccharomyces cerevisiae mitochondria) ( 26–28 ), both of which have been used in a variety of organisms (including Drosophila ) to induce DNA DSBs ( 29–32 ). We systematically analysed the nature of HEG-mediated integration and recombination events in A. gambiae and the effect of expressing these endonucleases on cell proliferation.
MATERIALS AND METHODS
Plasmids and cloning
The target plasmids pBC/SacRB S1/S2/S3 and pBC/SacRB P1 were constructed as follows: pBC/SacRB ( 21 ) was cut with SalI and PstI, ends filled in with T4 DNA polymerase and religated to remove a redundant EcoRI site. Linkers were created by annealing 5′ phosphorylated oligonucleotides (S1: AATTCATTACCCTGTTATCCTAG and AATTCTAGGGATAACAGGGTAATG; S2: AATTCGATAGGGATAACAGGGTAATTG and AATTCAATTACCCTGTTATCCCTATCG; S3: AATTCAATTACCCTGTTATCCCTACCG and AATTCGGTAGGGATAACAGGGTAATTG; P1: AATTCCGCTACCTTAAGAGAGTCG and AATTCGACTCTCTTAAGGTAGCGG), which were cloned into the now unique EcoRI site in the SacRB CDS. Selection against SacRB was performed in LB agar lacking NaCl containing 15% sucrose (w/v) and chloramphenicol (25 μg/ml).
To create pSL-SacRB Tet , the minimal 1.5 kb tetracycline (Tet) resistance cassette from pBR322 was amplified with the primers TTCAAGAATTCTCATGTTTGACAG and ATGAATTCTGCTAACCAGTAAGGCAACC and cloned into pBC/SacRB with EcoRI. The Tet cassette replaces the I-SceI site and disrupts the SacRB gene. The 3.4 kb SacRB Tet cassette was moved to pSLfa1180fa ( 33 ) using XhoI/XbaI to create pSL-SacRB Tet .
To create pDR-CMV-GFP, the CMV promoter from pEGFP-Ppo was amplified with the primers AAAGGGCCCTAGTTATTAATAGTAATCAATTACGGGGTCATTAG and AAAGAATTCGATCTGACGGTTCACTAAACCAGCTC and cut with ApaI and EcoRI. The fragment was ligated into partially ApaI and EcorI digested pDR-GFP ( 34 ). This replaces the 3′ part of the chicken β-actin promoter with CMV. To remove the remaining sequences of the chicken β-actin promoter, the resulting vector was cut with SnaBI and religated.
PSL-Act-EGFP was constructed in pSLFa1180fa to contain the 2.5 kb Drosophila Actin 5C promoter driving EGFP (BamHI/XbaI fragment) and the Drosophila Hsp70 terminator.
Cell culture and transfections
Cells from the stable anchorage-dependent A. gambiae cell line, Suakoko 4 (Sua 4.0) (Müller,H.M. et al. , 1999) were cultured in Schneider's Drosophila medium (Invitrogen) supplemented with 10% FCS (Invitrogen) and 200 U/ml penicillin and 200 μg/ml streptomycin sulphate (Invitrogen) in a cooled incubator at 27°C.
In vivo HEG activity assays were performed by lipid-mediated transient transfections (Effecten, Qiagen) of 1–3 × 10 5 cells with 2 μg/ml (culture volume) recipient plasmid and 4 μg/ml donor plasmid. When necessary, cells were heat-shocked for 1 h at 41°C, 24 h post-transfection to induce expression of I-SceI from the Drosophila Hsp70 promoter on the pP[v+, 70 I-SceI] plasmid. Total DNA was extracted 48 h post-transfection (Promega Wizard Genomic DNA purification Kit) and re-suspended at 25–40 μl. This preparation was used to transform Escherichia coli DH5α strain.
Embryo microinjections and preparation of low molecular weight DNA
Anopheles gambiae adult females (G3 strain) were allowed to deposit their embryos 72 h after a blood meal on a moist filter paper. Embryos were injected 60–120 min after oviposition, essentially as described ( 35 ). Embryos were injected with a mixture of I-SceI donor plasmid pP[v+, 70 I-SceI] (400 μg/ml), Tet donor plasmid pSL-SacRB Tet (500 μg/ml) and HEG target plasmid pBC/SacRB S1 (500 μg/ml). Between 150 and 200 embryos were microinjected and 12 h later embryos were heat-shocked for 1 h at 41°C and left to recover for 2 h. Low molecular weight DNA was extracted ( 21 ), re-suspended in 20 μl and 150 ng of the recovered DNA was used to transform E. coli DH5α strain.
Southern blot and primer extension
Genomic DNA was digested with ClaI in the presence and absence of I-PpoI. As a probe we used a 2 kb rDNA PCR fragment amplified from genomic DNA using the primers GCCGAAGCAATTAGCCCTTAAAATGGATG and CACCAGTAGGGTAAAACTAACCTGTCTCACG. The probe was P 32 labelled using the High Prime DNA labelling kit (Roche) and purified with ProbeQuant™ G-50 columns (GE Healthcare). For primer extension, genomic DNA was digested with HincII (a.k.a. HindII). The reaction was performed essentially as described ( 36 ) using the 5′ P 32 -labelled primer rPrex GTTAATCCATTCATGCGCGTCACTAATTAG and vent (exo-) polymerase (New England Biolabs). The reaction products were resolved on a 6% denaturing polyacrylamide gel. Results for both experiments were visualized using a FUJIFILM-FLA-5000 Phosphoimager (Fuji Photo Film Co. Ltd, Stamford, CT, USA). For in vitro digestions, we used commercially available I-PpoI (Promega) and I-SceI (New England Biolabs) enzymes.
Cell proliferation assay
Sua 4.0 cells were transfected with either of the two endonuclease plasmids (4 μg/ml) together with pIB/V5-His (2 μg/ml), conferring resistance to blasticidin (Invitrogen). Forty-eight hours post-transfection, blasticidin was supplemented to complete medium at 50 μg/ml. Cells were incubated in blasticidin for an initial proliferation period of 5 days, at which point they were harvested and re-seeded at 1.5 × 10 5 cells/ml and grown for a further 5 days. Transfections were performed in triplicates and cell numbers for each were counted in four replicates.
GFP gene conversion assay
Sua 4.0 cells were transfected with pDR-CMV-GFP (2 μg/ml) in the presence or absence of donor plasmid pP[v+, 70 I-SceI] (4 μg/ml). Twenty-four hours post-transfection, cells were heat-shocked for 1 h at 41°C. Gene conversion was measured by fluorescence activated cell sorting (FACS) 72 h post-transfection using a Beckman FACS Calibur; 50 000 size-dependent gated events were analysed for GFP fluorescence. Data were analysed using the FlowJo software package. Transfection efficiency was assessed by parallel co-transfections with an Actin 5C-DsRed plasmid.
Anopheles gambiae Sua 4.0 cells were transfected with pEGFP-Ppo, pEGFP-Ppo H98A and the control pSL-Act-GFP, and 48 h later cells were fixed for 5 min in 4% paraformaldehyde and permiabilized for 10 min with 0.1% Triton X-100. Nuclei were stained with DAPI (2 ng/μl) and actin filaments with Alexa546-phalloidin (1 U/ml, Invitrogen). Cell micrographs were taken at ×40 magnification using a Zeiss widefield microscope.
Development of a reporter system for HEG activity
To assess the functionality of HEGs in Anopheles cells, we developed a rapid and efficient reporter system to score HEG-mediated activity (site specific cleavage and homologous recombination events) in both A. gambiae cells and embryos. This assay, based on the interplasmid transposition assay utilized in insect cells and embryos to assess the activity of transposable elements ( 21 ), utilizes two sets of plasmids ( Figure 1 a). A donor plasmid directs the production of either I-SceI [pP(v+,70I-SceI)] or I-PpoI (pEGFP-PpoI), and the target plasmid (pBC/SacRB) contains the bacterial suicidal gene SacRB , coding for levansucrase, engineered to contain either the I-SceI or the I-PpoI recognition sequences. SacRB catalyses the hydrolysis of sucrose and the synthesis of levans, high-molecular-weight fructose polymers that accumulate in the periplasmic space, and are toxic to gram negative bacteria ( 37 , 38 ). HEG cleavage of its recognition sequence results in the inactivation of SacRB genes, which is then detected in bacteria transformed with plasmid DNA recovered from transiently transfected cells and injected embryos. Negative selection against the functional SacRB gene in medium containing sucrose is accompanied by positive selection for the chloramphenicol (Cam) resistance marker, which is also present on pBC/ SacRB ( Figure 1 ).
The I-SceI or I-PpoI recognition sites were inserted in the SacRB sequence after Glu67, downstream of the signal peptide required for secretion ( 39 ), at a position predicted not to be essential ( 40 ). We assessed the effect of inserting several recognition site variants in different frames on SacRB activity. None of the inserted amino acid variants (pBC/SacRB S1: ITLLSL , pBC/SacRB S2: D RDNRV I, pBC/SacRB S3: N YPVIP I for I-SceI; pBC/SacRB P1: ATLRE for I-PpoI) interfered with SacRB function as inferred by the continued inability of bacteria transformed with these vectors to grow on Cam supplemented with 15% sucrose (data not shown). We therefore concluded that the SacRB gene is tolerant to amino acid insertions at this position and a variety of recognition sequences of natural and reengineered HEGs could be inserted and tested using this approach.
I-SceI and I-PpoI homing endonuclease activity in A. gambiae cell lines
We transfected the plasmids pP[v+,70I-SceI] and pEGFP-Ppo together with their corresponding target pBC/SacRB S1 and pBC/SacRB P1 into A. gambiae Sua4.0 cell lines ( Figure 1 a). The plasmid pP[v+,70I-SceI] ( 30 ) expresses the I-SceI ORF with an N-terminal SV40 NLS and hemagglutinin (HA) tag under the control of the Drosophila Hsp70 promoter, which directs a significant inducible expression of I-SceI in Sua4.0 cells ( Figure 1 d, lanes 2 and 3). The plasmid pEGFP-Ppo expresses the EGFP-I-PpoI fusion protein containing a SV40 NLS ( Figure 1 a and d) under the control of the CMV promoter. An inactive variant pEGFP-Ppo H98A (Raymond Monnat personal communication) was used as a control. Cells were heat-shocked for 1 h at 41°C 24 h post transfection and DNA extracted after an additional 24 h. Under these experimental conditions, we observed that the target plasmids were cut in vivo by the corresponding endonuclease ( Figure 1 b). Digestion of DNA extracted from transfected cells with purified endonucleases in vitro revealed the presence of a fraction of endonuclease resistant plasmids ( Figure 1 b, lanes 6–10, white arrow), possibly generated by non-homologous religation of HEG-mediated cleavage events.
To recover and analyse non-homologous end joining (NHEJ) events generated in mosquito cells, total DNA (which includes plasmid DNA) extracted from transfected cells was used to transform bacterial cells plated on Cam or Cam/Suc selective media. Compared to control experiments (target plasmid only), transfection of mosquito cells with I-SceI- or I-PpoI expressing plasmids increases the number of Cam/Suc resistant clones by 15- and 8-fold, respectively ( Figure 1 c). By comparing the colony numbers on Cam and Cam/Suc plates of bacteria transformed with recovered DNA from A. gambiae cells, we found that ∼0.5–2% of recovered plasmids allowed growth on Cam/Suc as compared to Cam alone. This indicates that although HEG cleavage is efficient ( Figure 1 b) only a small number of these plasmids are subsequently religated.
Plasmids from Cam/Suc-resistant bacterial cells were isolated and digested with BamHI and HindIII endonucleases. Only plasmids showing a 1.8 kb band, indicating the presence of an intact SacRB gene ( Figure 1 a), were analysed by sequencing ( Figure 2 ). This step was undertaken to ensure that sucrose resistance is due to HEG-mediated disruption of SacRB, as we occasionally observed the growth of Cam/Suc-resistant bacteria in the absence of a HEG expression vector ( Figure 1 c). Digestion of plasmids with BamHI/HindIII in the presence of I-SceI or I-PpoI showed that all recovered clones were also resistant to endonuclease cleavage in vitro (data not shown). The sequence of the regions surrounding the HEG cleavage sites was analysed in plasmids recovered from Cam/Suc-resistant bacteria ( Figure 2 ). All 27 sequenced clones showed deletion events of variable sizes, ranging from 1 to 80 bp, close to the predicted HEG cleavage site. We did not observe any nucleotide insertions, contrary to reports of NHEJ repair products in human cell lines ( 41 , 42 ).
I-SceI expression induces homologous repair events in A. gambiae cells
Copying of HEGs from one allelic site to the other requires that cleavage of the target site is followed by gene conversion events using the HEG-containing allele as a repair template. We therefore adapted our reporter system to test for the occurrence HEG-induced gene conversion events in mosquito cells. We used plasmid pDR-CMV-GFP ( 34 ) ( Figure 3 a), which contains the CMV promoter 5′ of a non-functional eGFP gene. In this gene, 11 bp of the GFP CDS have been deleted and replaced by an I-SceI site which also introduces two stop codons. The plasmid also carries a second promoterless and inactive GFP locus, which lacks 220 bp of the GFP C terminus, functioning as a repair template. Cells were transfected with pDR-CMV-GFP or co-transfected with pDR-CMV-GFP and pP[v+,70I-SceI]. As expected, transfection of mosquito cells with pDR-CMV-GFP alone did not result in a high rate of spontaneous repair of the GFP gene as shown by the low frequency (0.02%) of fluorescent cells ( Figure 3 b), whereas we found that co-transfection with pP[v+,70I-SceI] increased the number of GFP+ cells 50-fold (0.99%). Adjusted for transfection efficiency, GFP positive cells represent 0.15% (pDR-CMV-GFP) and 10.9% (pDR-CMV-GFP + pP[v+,70I-SceI]) of transfected cells. Co-transfection with the I-SceI expressing plasmid increases the average and peak intensity of fluorescence of GFP+ cells ( Figure 3 b, lower right panel and data not shown), thus suggesting the presence of a higher number of repaired plasmids per cell in these experiments.
I-SceI is expressed and active in A. gambiae embryos
We also investigated HEG activity in A. gambiae embryos. Briefly, I-SceI donor plasmid pP[v+,70I-SceI] and pBC SacRB S1 were co-injected into preblastoderm embryos. I-PpoI was not used in these experiments due to its toxic effect on insect cells (see subsequently). One day post-injection embryos were heat-shocked for 1 h at 41°C and low molecular weight DNA was extracted 2 h later. Selection for target plasmids modified by NHEJ was performed as described above for cells. As previously observed in mosquito cells, co-injection of the plasmid pP[v+,70I-SceI] increased the number of Cam/Suc-resistant bacteria colonies ∼10-fold: 0.8% of the recovered pBC/SacRB S1 pool were able to grow on Cam/Suc, a similar frequency to that obtained from cells. We analysed the sequence of the I-SceI recognition site in 20 plasmids recovered from injected embryos ( Figure 2 b). In contrast to what we observed in mosquito cells, both deletion and insertion events were found. One of the recovered plasmids contained a 66 bp insertion homologous to an A. gambiae genomic sequence.
I-SceI expression induces inter-plasmid gene conversion events in A. gambiae embryos
In order to test for HEG-induced gene conversion events in A. gambiae embryos, we used the plasmid pSL-SacRB Tet . This plasmid contains a Tet resistance cassette inserted into the SacRB gene on a pSL backbone ( Figures 1 a and 3 c). Therefore, this plasmid shares two regions of homology with the pBC target plasmids which flank the Tet resistance cassette but contains a different backbone ( Figure 3 c). As shown in Figure 3 c and d, the outcome of a perfect gene conversion event would generate a plasmid that contains both Cam and Tet resistance cassettes but a non-functional SacRB gene. Plasmids pBC/SacRB S1, pSL-SacRB Tet and pP[v+,70I- Sce I] were co-injected into pre-blastoderm stage embryos and DNA extracted as described above. As a control a second set of injections was performed excluding the pP[v+,70I-SceI] plasmid. In control experiments, where only the HEG target plasmid was injected we obtained six Cam/Suc bacterial colonies and none of these were able to grow when plated on Tet. To demonstrate that gene conversion cannot be achieved by the bacteria alone, bacteria were transformed with the injection mixture and plated on Cam, Cam/Suc, Tet as well as Cam/Suc/Tet. Whereas bacteria transformed with injection mix can grow on Cam or Tet, no colonies were obtained on Cam/Suc or Cam/Suc/Tet plates.
From injection experiments performed with all three plasmids (pBC/SacRB S1, pSL-SacRB Tet and pP[v+,70I- Sce I]) we obtained 98 bacterial colonies growing on Cam/Suc in the initial selection round. Of these 11 were also able to grow on Tet (11.2%). Further analysis revealed that two of them were unable to grow on ampicillin (Amp) (the Amp resistance gene is found in the pSL-SacRB Tet but not in the pBC/SacRB S1 backbone). We analysed these two colonies by restriction endonuclease digestion ( Figure 3 e) and sequencing. The sequence analysis revealed the occurrence of perfect gene conversion events, in which the Tet cassette had been inserted into pBC/SacRB S1 as predicted. The nine clones that were able to grow on Amp were resistant to cleavage with I-SceI but they appeared to be larger than the expected product. All nine plasmids were cleavable by AscI, a rare cutting enzyme that flanks the SacRB Tet cassette ( Figures 1 a and 3 c), indicating that larger parts of pSL-SacRB Tet had been transferred during the process of homologous recombination.
The A. gambiae 28S rDNA gene contains a conserved I-PpoI recognition site
I-PpoI mediates homing of intron 3 (Pp LSU 3) in the extrachromosomal nuclear rDNA of the acellular slime mold Physarum polycephalum ( 23 ). This region of the 28S rDNA is highly conserved in eukaryotes and I-PpoI has been shown to cleave human rDNA repeats in vivo ( 42 ). In order to determine whether an I-PpoI site is present in A. gambiae rDNA, the 28S gene of A. gambiae was assembled from the sequence of 31 cDNA sequences obtained from the AnoEST database. Sequence analysis indicated the presence of the full 29 bp I-PpoI recognition site. The assembled gene sequence was used to design a set of primers (rDfwd and rDrev, Figure 4 a) that amplify a 2 kb fragment of the A. gambiae 28S rDNA gene containing the I-PpoI site ( Figure 4 a). The PCR product from A. gambiae strains KWA and G3 was digested with I-PpoI to confirm the presence of the site in both strains (data not shown).
I-PpoI cleaves chromosomal rDNA repeats in vivo
To determine whether I-PpoI can bind and cleave chromosomal rDNA repeats in vivo , A. gambiae Sua 4.0 cells were transfected with pEGFP-Ppo or pEGFP-Ppo H98A plasmids that express a functional and an inactive version of I-PpoI respectively. Genomic DNA was extracted 24 and 48 h post-transfection and digested with ClaI or ClaI and I-PpoI ( Figure 4 a). ClaI is predicted to cut 1.5 kb up and 1.4 kb downstream of the I-PpoI site. Southern blot analysis using the 2 kb 28S rDNA PCR product as a probe indicated that chromosomal 28S rDNA was cut efficiently by I-PpoI but not I-PpoI H98A ( Figure 4 b). The probe hybridized to several fragments, the shortest of which corresponds to the expected 2.9 kb ClaI fragment. All but one of these fragments were cleaved by I-PpoI ( Figure 4 b). This seems to indicate some heterogeneity in the 28S genes and we obtained a similar result with another enzyme combination (data not shown). Presumably this is a result of retrotransposon insertions downstream of the I-PpoI site ( 43 ). The appearance of resistant bands, retrieved when using I-SceI in Anopheles cells ( Figure 1 b, lanes 7 and 10), was not observed when analysing genomic rDNA cleaved with I-PpoI in vitro ( Figure 4 b, lanes 7–10). To confirm the result obtained by Southern blotting, we performed a primer extension analysis using primer rPrex, which binds 21 bases downstream of the I-PpoI cleavage site ( Figure 4 a). Figure 4 c shows that I-PpoI expression causes premature stops of the primer extension at a position corresponding to the I-PpoI site. These data together indicate that I-PpoI is able to efficiently cut mosquito rDNA in transfected cells.
I-PpoI but not I-SceI expression causes cytotoxicity in A. gambiae cell lines
In several mosquito species including A. gambiae , the 28S rDNA genes are clustered as tandem repeats on the X chromosome. The use of I-PpoI offers therefore the possibility to selectively disrupt the X chromosome. While I-PpoI cleaves the essential 28S rDNA genes, I-SceI is not predicted to have a target site in the A. gambiae genome. To study the effect of expressing these two HEGs on the viability of mosquito cells, we transfected Sua 4.0 cells with pP[v+,70I-SceI], pEGFP-Ppo, pEGFP-Ppo H98A or pSL-Act-EGFP, together with the plasmid pIB/V5-His, which confers resistance to the translation inhibitor blasticidin S. Cell proliferation was assessed as a measure of growth in medium containing 50 μg/ml blasticidine, which will kill all cells that do contain the pIB/V5-His plasmid. The results of these experiments are shown in Figure 5 a. While neither I-SceI nor mutant I-PpoI H98A expression interferes with cell proliferation, I-PpoI expression leads to growth arrest. Expression of the I-PpoI-GFP fusion protein seemed to induce nuclear fragmentation ( Figure 1 d) as observed by fluorescent microscopy. It appears that the nucleoli of cells transfected with pEGFP-Ppo ( Figure 5 b) are disintegrating. This can be explained by the fact that ribosomal DNA repeats which are cut by I-PpoI form the nucleolus organizer regions of the cell.
Whereas natural HEGs have evolved under constant selection pressure towards reduced toxicity and are selected to cut only at their native homing sites, engineered HEGs might be less specific. Experiments with designed zinc finger nucleases have shown that expression of these enzymes is often accompanied by general toxicity ( 44–46 ). Our assay confirms that I-SceI, although highly active, is not toxic when expressed in cells. I-PpoI on the other hand cleaves the essential rDNA genes and induces proliferation arrest, as has been described for human cells ( 42 ). The HEG activity, GFP reporter and cell proliferation assays described can be applied to the initial assessments for the suitability of any engineered HEG candidates designed against Anopheles genes.
Our results demonstrate that two different HEGs, I-SceI and I-PpoI, when expressed in A. gambiae cells and embryos, retain their ability to recognize and cleave their cognate target sequences. Analysis of I-SceI and I-PpoI function using an interplasmid activity assay in cultured mosquito cells showed that, in the absence of a template for homologous recombination, recovered target sites contained deletions ranging from 1 to 80 bp created by HEG-induced cleavage and subsequent NHEJ repair. Sequences recovered from embryos injected with the I-SceI gene contained both insertions (1–66 bp) and deletions (2–62 bp). It is not yet known whether the observed bias towards deletions in cells reflects a feature of Anopheles biology or whether it is a consequence of our plasmid-based selection system. All targeted sites analysed were resistant to subsequent I-SceI and I-PpoI cleavage in vitro as nucleotide deletions or insertions in the recognition site significantly impair endonuclease activity ( 47 ).
Recently, it has been shown that HEGs can be engineered to confer new sequence specificities ( 15–18 ), thus offering the possibility of using these highly specific nucleases for gene disruption and gene therapy ( 48 , 49 ). Our findings indicate that engineered HEGs could be used to selectively disrupt mosquito genes either in the germ line or in a tissue-specific manner depending on the spatial and temporal expression pattern of the driving promoter.
It has been suggested that the ‘selfish’ genetic behaviour of HEGs could be exploited to rapidly spread a genetic modification affecting mosquito vectorial capacity from a few individuals to an entire population ( 11 , 19 ). This process would require both DNA cleavage and repair via homologous recombination at the targeted site using the HEG-carrying allele as template. To determine the feasibility of such an approach, we determined whether repair of HEG-mediated cleavage occurred by homologous recombination when appropriate template sequences are present. A functional study carried out in mosquito cells showed, that in the presence of I-SceI, a plasmid carrying a transcriptionally silent and incomplete eGFP coding sequence could function as template to repair (as inferred by the gain of fluorescence) another non-functional eGFP sequence disrupted by an I-SceI cleavage site. The frequency of homologous repair events in mosquito cells is comparable to that previously described in mammalian cell lines ( 34 ). To extend this observation, we adapted the SacRB-based interplasmid assay to score homologous recombination events. We used for this purpose a relatively large conversion cassette containing a Tet resistance transcription unit (1.5 kb), in an attempt to simulate the size of a complete HEG element. While no recombinant clones were observed in the absence of the HEG, the precise transfer of the Tet cassette between two plasmids was found in 2 out of 98 events. We also observed nine recombination events showing larger insertions encompassing regions of the pSL plasmid flanking the Tet cassette. We speculate that these events arose by imprecise extensions during homologous recombination. Results from these experiments suggest an overall frequency of homologous repair in embryos of ∼10%. Our Tet selection system and GFP reporter system only allow us to recover events that have resulted in the complete insertions of the Tet/GFP sequences. The actual rate of homologous repair might be higher as the occurrence of putative recombination events that resulted in the transfer of smaller fragments of Tet or GFP can be expected but were selected against. Also, in none of our experiments were the repair donor fragments flanked by HEG half sites, a fact that might reduce the rate of homologous repair. In our experiments, DNA DSBs are induced in vivo by the activity of HEGs. Efficient homologous repair in A. gambiae Ag55 cells has also been demonstrated by the transfection of linear recombination substrates ( 50 ).
Together, these findings indicate in two highly predictive experimental systems that HEG would retain their gene-driving activity in mosquitoes. Furthermore, HEG-mediated cleavage of a target gene in embryos followed by gene conversion using a repair template provided in trans could represent an alternative approach to traditional transposon or integrase-mediated germ line mosquito transformation.
We have shown that I-PpoI specifically recognizes and cleaves a vital target sequence in A. gambiae cells inducing cell proliferation arrest and presumably cell death. The I-PpoI recognition site is located in the 28S rDNA gene within domain V of the 28S rRNA, a region that forms the peptidyl transferase centre of the ribosome. This sequence is among the most conserved sequences in the entire eukaryotic kingdom. In some mosquito species, including at least two members of the A. gambiae complex ( 51–53 ), the rDNA repeats are exclusively located in the centromeric region of the X chromosome. This arrangement raises the possibility of using I-PpoI to selectively attack the X chromosome: if specifically expressed during male meiosis, I-PpoI could distort the sex ratio towards males by selectively incapacitating X-carrying spermatozoa ( 19 ). Natural driving Y chromosomes in Aedes aegypti and Culex pipiens have been described and can produce extreme sex ratios of more than 90% males ( 54 ). Although the molecular details of how these distorters act are unknown, cytological evidence suggests that they are associated with breaks in the X chromosome during male meiosis I ( 55 , 56 ). These observations are relevant for the development of malaria vector control measures. A male sex distorter combining the specificity of I-PpoI together with its potential genetic drive activity if placed on the Y chromosome could knock down a wild-type Anopheline population in a few tens of generations ( 57 ).
We are indebted to Raymond Monnat for his generosity in providing the I-PpoI expression vectors as well as pDR-GFP and for helpful comments. We would also like to thank Kalle Magnusson, Sofia Pinto, Michael Povelones and Rita Tewari for technical help. F.C. is currently sponsored by an MRC Career Development Award. Funded by a Grant from the Foundation for the National Institutes of Health through the Grand Challenges in Global Health initiative. Funding to pay the Open Access publication charges for this article was provided by the Foundation for the National Institutes of Health through the Grand Challenges in Global Health initiative.
Conflict of interest statement . None declared.