Heritable genome editing with CRISPR/Cas9 induces anosmia in a crop pest moth

Abstract

Lepidoptera suffer critical lack of genetic tools and heritable genome edition has been achieved only in a few model species. Here we demonstrate that the CRISPR/Cas9 system is highly efficient for genome editing in a non-model crop pest Lepidoptera, the noctuid moth Spodoptera littoralis. We knocked-out the olfactory receptor co-receptor Orco gene to investigate its function in Lepidoptera olfaction. We find that 89.6% of the injected individuals carried Orco mutations, 70% of which transmitted them to the next generation. CRISPR/Cas9-mediated Orco knockout caused defects in plant odor and sex pheromone olfactory detection in homozygous individuals. Our work genetically defines Orco as an essential OR partner for both host and mate detection in Lepidoptera, and demonstrates that CRISPR/Cas9 is a simple and highly efficient genome editing technique in noctuid pests opening new routes for gene function analysis and the development of novel pest control strategies.

Heritable genome editing with

Fotini A. ,*, Christelle ,*, Marie-Christine Franois, Anne de Cian, Corinne Royer, det

Spodoptera littoralis Orco Orco Orco

Lepidoptera represent more than 10% of the total described species of living organisms. They are studied in all elds of biological research, but suer critical lack of reverse genetic tools. RNA interference (RNAi) approaches are usually inefficient in these species1 and heritable mutagenesis has been established in only a limited number of model species. First developed in Bombyx mori using a transposon approach2, such germ line transformation in other Lepidoptera has been mostly inefficient. Recently, new genome editing tools, such as Zinc Finger Nucleases (ZFN), Transcription Activator-Like Eector Nucleases (TALEN) and the Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/Cas9 System3, open new routes to manipulate genes in a diversity of organisms. These methods use sequence-specic endonucleases to induce double-strand DNA breaks (DSBs) in a target gene. The non-homologous end-joining (NHEJ) DNA repair pathway is then attempting to repair the lesion, generating errors that eventually interrupt the open reading frame, inactivating the gene by incomplete protein translation3. ZFN-mediated mutagenesis has been recently reported for B. mori4 and Danaus plexippus5; TALEN has been developed in B. mori68; and CRISPR/Cas9 proved to be efficient in this species9,10 as well as in the model butteries Papilio xuthus11 and Danaus plexippus12. The latter approach is in full expansion13,14 and represents a groundbreaking milestone for non-model species. The present study aims at extending its use in non-model Lepidoptera.

Focusing on the noctuid crop pest Spodoptera littoralis, a highly polyphagous pest causing important damage on cotton and vegetable crops in Europe, Asia and Africa, we used CRISPR/Cas9 to target a gene that could interrupt the chemical communication of this species, the receptor co-receptor Orco. As demonstrated for its orthologue in the fruit y Drosophila melanogaster, Orco forms heterodimers with olfactory receptors (ORs)15,16 that

act as odor-gated ion channels at the membrane of olfactory sensory neurons1719. In D. melanogaster, the ORs are involved in host/partner odorant recognition while another family of receptors, the ionotropic receptors (IRs), are mostly involved in the recognition of fermentation products, such as acids and amines20,21. Orco disruption

via knock-out or knock-down in Diptera models16,22,23 and other insect species2426 induced impaired olfactory detection capacities. A noctuid Orco could rescue the olfactory abilities of a D. melanogaster Orco mutant27, but the essential nature of Orco in Lepidoptera has never been challenged yet. Thus, Orco represents a good target gene to investigate olfactory pathways in Lepidoptera and gives a unique opportunity for testing the development of genome engineering strategies in crop pest moths.

*

R A

P

SCIENTIFIC REPORTS

1

www.nature.com/scientificreports/

Figure 1. S. littoralis Orco (SlitOrco) gene structure and guide RNA positions. SlitOrco consists of 10 exons (grey boxes, E) and 9 introns. The number of amino acids encoded by each exon is indicated in the box. Two gRNAs, cr104 and cr105, were designed in exon 2 (E2) and one, cr106, in exon 4 (E4). Exon 2 and 4 sequences are detailed (underligned parts) showing positions of gRNAs (in blue ; PAM in italic and in box). Primers used for PCR amplication in the genotyping assay are in red italic.

Figure 2. Comparison of the activity of the dierent guide RNAs (gRNAs). Three gRNAs (cr104, cr105, cr106) were co-injected with the Cas9 protein in dierent egg batches (9 batches illustrated for cr104 and 3 batches illustrated for each of the two other guides). Batches of emerging larvae were genotyped by the T7EI assay following genomic DNA extraction and PCR amplication of an Orco fragment (for fragment sequences and primers, see Fig.1). DNA modication could be obtained only for cr104 and in all egg batches, as revealed by the gel pattern: the highest band corresponds to the wild type gDNA amplication and the two small bands (arrows) result from T7EI action. Water injected larvae (H2O) and no template () controls are shown for each primer combination. The ladder used is the Quick-load 2-Log DNA Ladder (0.110.0kb, BioLabs).

We report here successful CRISPR/Cas9 genome editing in a non-model pest crop moth. This work provides new perspectives for functional genomics in such emerging-species, and oers opportunities for new pest control strategies.

Orco The genomic sequence of the S. littoralis Orco gene was identified via PCR, according to the B. mori Orco sequence (SilkDB: http://silkworm.genomics.org.cn/) assuming that intron positions would be similar. The gene structure consisted of 10 exons and 9 introns, comparable in size and positions with its B. mori orthologue (Fig.1).

As in most non-model species, S. littoralis embryogenesis timing is not known. We optimized early events by both selecting freshly laid eggs (within the rst hour aer oviposition) and injecting Cas9 protein instead of RNA, avoiding translation delay. As reasoned in Merlin et al.5,

the objective was to target the beginning of the mitotic stages of the eggs, assuming that egg fertilization occurs during oviposition when sperm from the spermatheca passes through the micropyl28.

Three gRNAs (cr104, cr105 and cr106) were designed (Fig.1), injected in eggs and tested by genotyping pools of newly hatched larvae, to select the more efficient gRNA. PCR amplication of an Orco fragment encompassing the targeted sequence was followed by a T7 endonuclease I (T7EI) assay that recognized and cleaved heteroduplex DNA. By this assay, batches of larvae without any mutation revealed only one band corresponding to the wild type PCR amplied Orco fragment and batches of larvae containing mutations revealed two additional bands that resulted from T7EI action (Fig.2). Out of the three gRNAs, only one was active at inducing mutations (cr104), and mutations were observed in all egg batches injected (Fig.2). This gRNA was used further for individual studies by injection with Cas9 in 2569 eggs. In S. littoralis, eggs are laid in batches and percentages of injection

SCIENTIFIC REPORTS

2

www.nature.com/scientificreports/

Number Mutated sequence number (see Fig 3) Total % Injection injected eggs 2569hatched larvae 694 27 G0 G0 T7EI genotyped larvae 58

G0 mutated larvae (T7EI gel) 52 89.6 G0 sequenced larvae 33G0 larvae with multiple mutations

(uncharacterized sequences) 17

G0 larvae with multiple mutations

(characterized sequences) 7 6,7,9,10,11,13,14

G0 larvae with unique mutation

(characterized sequence) 9 15,8,12,15,16

G0 sequence mutations provoking

ORF frameshis 11 27,9,11,14,15,16

Founders G0 mutant lines tested for founder

eects (mutation transmission) 10 16,9,11,15,16

G0 founders 8 1,36,9,15,16 80 G1 genotyping Mutated sequences (as in Fig.3) 1 2 3 4 5 6 9 11 15 16

G1 genotyped larvae (T7EI) 30 30 30 30 30 30 30 30 30 30 300 G1 heterozygote larvae (T7EI gel) 10 10 5 6 10 5 11 11 9 13 33 G1 heterozygote larvae

(sequencing) 10 0 3 2 7 5 11 0 6 13 57 Germline

mutation rate

survival (694 larvae hatched, 27%) ranged from 2.5% to 46% depending on the injected egg batch (Table1). These percentages were similar to those observed for control eggs injected with water or just punctured (1.5% to 45%, 24% on average), suggesting that the protein/gRNA injected mix was not toxic. Aer individual genotyping, positive larvae were further reared until adults.

All data described here are summarized in Table1. Among the larvae that emerged from cr104/Cas9-injected eggs, 58 were individually genotyped, and 89.6% (52 larvae) of those presented mismatch at the Orco target region (T7EI gel assay). Thirty-three were sequenced, out of which 24 harboured multiple mutations. Due to sequence superposition in chromatograms, 17 of them could not be characterised. Seven (sequences 6, 7, 9, 10, 11, 13 and 14, Fig.3) consisted of only two or three superposed sequences and could be easily characterised by visual curation of chromatograms. Further sequencing of G1 individuals that all carried a single mutated sequence conrmed the predicted sequences. All of the seven G0 superposed sequences were carrying the same background sequence, consisting of a 6 base pair (bp) deletion (sequence 1, Fig.3). The IUPAC code was used when the superposing bases could not be clearly identied, typically when three sequences were superposed (Fig.3). For nine larvae (out of 33, 27%), unique mutated Orco sequences were detected (sequences 15, 8, 12, 15 and 16 in Fig.3). Sequence 1 corresponded to the 6bp deletion. This frequent 6bp deletion mutation (50%: 8 out of 16 characterized mutated larvae) does not induce stop codon and has possibly no consequence on protein function. Eleven mutated sequences resulted in frameshi in the open reading frame generating a premature stop codon (sequences 27, 9, 11, 14, 15 and 16 in Fig.3), and are thus expected to generate non-functional Orco proteins.

Ten injected individuals harbouring a unique or two dierent mutations (sequences 16, 9, 11, 15 and 16, Fig.3) were tested for their founder (mutation transmission) capacities at the G1 generation by cross with wild-type individuals. From each cross, thirty G1 larvae randomly chosen were genotyped. Out of the ten putative founders, 8 (80%, sequences 1, 36, 9, 15 and 16 in Fig.3) produced heterozygous mutant progeny at percentages ranging from 6.6 to 43.3% (Table1), revealing mutation inheritance. Three G1 lines with mutations causing stop codons (sequences 5 , 15 and 16 in Fig.3) as well as the line with the 6 bp deletion were mated with their siblings to give second-generation mutants (G2 generation). In G2, around sixty larvae were genotyped per line. Contrary to the expected Mendelian inheritance ratio (% 25/50/25) of homozygous/heterozygous mutants/wild-types, we got an average ratio of % 5/63/32 (Table1). Interestingly, we could not obtain any progeny aer crossing homozygous mutated males and females; females laid scarce and dispersed eggs lacking scales protection (Supplementary Fig. S1), which dried and never hatched.

% heterozygote larvae

(sequencing) 33.3 0 10 6.6 23.3 16.6 36.6 0 20 43.3 19G2 G2 genotyped larvae (T7E1) 61 60 59 64 244 100

G2 +/+ larvae (T7EI) 18 23 19 18 78 32 G2 /+ larvae (T7EI) 40 37 37 40 154 63 G2 / larvae (T7EI, veried by

sequencing) (see Table2 for details) 3 0 3 6 12 5

Table 1. Summary of the experiments from G0 to G2.

SCIENTIFIC REPORTS

3

www.nature.com/scientificreports/

Figure 3. Sequences of the CRISPR/Cas9 induced mutations in the S. littoralis Orco gene obtained at G0.

The PAM motif is highlighted as a grey box and the 20bp cr104 guide RNA (gRNA) is shown with blue letters. The sequence modications () are highlighted in red letters (insertions) and in red dashes (deletions). The sequences that were obtained from samples containing also sequence 1 are numbered with green numbers. The sequences with mutations that produced stop codons within the Orco coding sequence are designated with asterisks (*).

Sexe mutation number / +/ +/+

Males 1 3 12 10

Females 1 0 20 5Males 5 0 12 2Females 5 0 13 7Males 15 1 9 5Females 15 2 10 6Males 16 2 14 8Females 16 4 14 4

Table 2. Details of G2 genotyping by Orco sequencing. Numbers of homozygote mutants (/), heterozygote mutants (+/) and wild-type individuals (+/+) obtained at the G2 generation per Orco mutation (see Fig.3 for mutation numbers).

Orco We used electroantennography (EAG) to evaluate the eect of Orco knock-out (KO) on the olfactory abilities of 163 G2 individuals obtained from four mutant lines (1, 5, 15 and 16 mutant lines, Fig.3), including / individuals (no mutation characterized), /+ individuals (heterozygote mutants) and +/+ individuals (homozygote mutants), all genotypes being conrmed by sequencing (Table2). EAG was also conducted on 14 individuals from our rearing colony and whose parents did not encounter any egg injection (further designed as wild type). One antenna from each individual was consecutively stimulated with dierent plant-odorants known to be detected by S. littoralis antennae (E-ocimene, Z3-hexenyl-acetate, E2-hexenol, benzyl alcohol and phenyl acetaldehyde)29, the main sex pheromone component Z9,E11-tetradecenyl acetate (Z9,E11-14:Ac)30 and one acid (propionic acid) as well as their specic solvants (paraffin oil, hexane and water respectively). Plant-odorants and sex pheromones are supposed to be transduced via the OR-Orco pathway, while acids are supposed to be transduced via the IR pathway16,20,21,23,31. The use of these stimuli allowed us to test the eects of Orco KO on both pathways.

Non-mutated G2 (/) presented similar EAG responses as wild type adults when stimulated with plant-odorants and the main pheromone component32. Furthermore, we evidenced propionic acid detection in both wild type and / G2 males and females. The heterozygotes responded to all odorants with no statistically signicant dierence from the wild-types (p> 0,001, Fig.4a,b), except for E-ocimen that induced very low EAG responses.

Whatever the sex, homozygous mutants were anosmic to the plant odorants and Z9,E11-14:Ac; the responses were significantly different from wild-type, p < 0,001, and close to zero. However, homozygous mutants responded as wild-type to the propionic acid (Fig.4a,b). In all experiments, no dierence was found between male and female, at the exception of the pheromone responses. Heterozygous and wild-type male responses to the pheromone were at least twice as high as the female responses, as previously reported32. The responses of individuals carrying the homozygous 6bp mutation were not statistically dierent (p> 0,001) from the responses of heterozygotes and wild-types, although a weak anosmic phenotype was observed against the Z9,E11-14:Ac (Fig.4c). This suggests that Orco was still functional for this mutation but that the two lacking amino acids (positions 61 and 62) may be important for complete functioning of the protein. Although these amino acid positions have yet never been assigned a functional role in Orco, it is known that single amino acid modication may aect Orco functioning33.

To test whether pheromone anosmia in Orco KO mutants would be responsible for the unsuccessful matings, we investigated alteration of mating performance in antennectomised wild type animals. Antennal ablation has been indeed reported to induce efficient anosmia in insects34,35. Mating capacity was not altered when the females were deprived of their antennae (100% mating and fertilized eggs, as for intact insect pairs), but was totally disrupted when we used antennectomised males (0% of mating observation and no larvae hatching).

SCIENTIFIC REPORTS

4

www.nature.com/scientificreports/

Figure 4. Electrophysiological impact of Orco KO in S. littoralis antennae. Electroantennogram responses (EAG, in mV SEM, the response to solvent was subtracted) of S. littoralis antennae isolated from wild-type (WT, non-injected parents) and CRISPR/Cas9 G2 individuals toward plant odorants (10g), the main pheromone component (Z9,E11-14:Ac, 1g) and propionic acid (10g). (a) Female responses: WT (black), +/+ (grey, CRISPR/Cas9 G2 without any mutation), +/ (pink, CRISPR/Cas9 G2 with heterozygous mutation) and / (red, CRISPR/Cas9 G2 with homozygous mutation) (sequences 5, 15 and 16 as in Fig.3, causing a truncated Orco protein); (b) Male responses as in (a); (c) Responses from males generated from parents carrying Orco sequence 1 (2 amino acid loss compared to the Orco wild-type sequence, no stop codon produced). Dierent letters above each odorant response indicate signicant dierences (t-test; p< 0.001); a is dierent from band ab is not dierent from either a or b. PAA: 2-phenyl acetaldehyde; Ac: acetate ; OH: alcohol; n: number of individuals tested for each genotype.

Our experiments demonstrate that CRISPR technology is an eective method for generating targeted knock-outs in a non-model crop pest noctuid. We took advantage of the large number of eggs that pests usually lay in a short period, while the low number of eggs is oen the major limiting parameter of induced mutagenesis/transgenesis success in other Lepidoptera. We favoured early NHEJ events by injecting gRNA/Cas9 protein within the rst hour aer oviposition, as close as possible to the one-nucleus-stage or the pronuclei fusion stage and surely before the blastoderm cellularisation. However, 72.7% (24 out of 33 sequenced G0 larvae) of the eggs showing eective mutation were certainly injected aer the one-nucleus-stage because they carried several mutations in the Orco sequence. Anyhow, this strategy appears successful since genotyping and sequencing analyses performed on G1 and G2 conrmed the germline transformation, demonstrating a stable transmission of the mutations. As usually observed with CRISPR/Cas9, we obtained both deletion and insertion events. Insertions usually correspond to

SCIENTIFIC REPORTS

5

www.nature.com/scientificreports/

plasmid or genome sequence integrations. Since we did not inject any plasmid constructs, it is probable that the long fragments inserted corresponded to genome fragment, but this could not be veried since no genome is yet sequenced for S. littoralis.

The CRISPR/Cas9 system induced Orco mutations at very high efficiency (89.6%) in adults derived from injected embryos. 80% of mutated individuals transmitted mutation to their progeny with up to 43.3% germline transmission efficiency. Comparable efficiencies were obtained in B. mori using CRISPR/Cas9 to target the BmBLOS2 gene10: 95.6% of mosaic phenotype could be obtained at G0 and germline transmission efficiency was 35.6%. Using TALEN and ZFN mutagenesis methods, 46% and 72% of G0 mutants, respectively, were obtained from B. mori injected embryos4,6. However, the TALEN mutagenesis gave 31% founders and a germline transmission efficiency of 61%, while the ZFNs showed very low efficiency in transmitting mutations to the G1 (9%). In monarch butteries, 50% germline mutations has been achieved using ZFNs5.

In our experiment, one specic mutation (a 6bp deletion) was found in 50% of the G0 genotyped individuals. Sequence inspection revealed that this deletion could have resulted from joining of 2 small homology stretches of sequence GTAT that ank the Cas9 cleavage site. The prevalence of this mutation may have been favoured by microhomology-mediated end-joining. This is reminiscent of studies in human cells that reported a high frequency of mutations involving stretches of sequence microhomology36 and suggests that when designing guide RNA for gene inactivation in S. littoralis, special care should be taken to avoid generating in this manner a majority of in-frame deletions as found here. As shown previously36, this can be done using appropriate bioinformatic tools, for example available at crispor.tefor.net or rgenome.net.

Targeted mutagenesis allowed us to knock-out the Orco gene in S. littoralis. The function of this gene in Lepidoptera has not yet been deeply investigated. As a noctuid Orco could rescue the detection capacities to some odorants in a Drosophila Orco mutant27, one could suggest that Orco has the same function in Lepidoptera. Here, we validate this hypothesis since the Orco KO homozygous mutants we obtained were not able to detect the tested host plant odorants. No eect was observed in Orco KO heterozygous mutants, as observed in mosquitoes22,

suggesting that both mutated alleles are required for efficient KO. Moreover, adults were also anosmic to the main sex pheromone component. Orco coupling with the sex pheromone receptors (PRs), a well-dened subclass of Lepidoptera ORs, has been previously debated. In B. mori, the two PRs and Orco do not co-express in the same olfactory receptor neurons37, whereas co-expression would be expected if they form heterodimers. Our results suggest that lepidopteran PRs function as ORs, via interaction with Orco.

Orco KO homozygous mutants still detected the acid as did the wild-types. Acid sensing in Lepidoptera has not been investigated before, although they are known to express as many IRs as D. melanogaster, if not more38,39.

Our study suggests that S. littoralis adults detect acids via an Orco independent pathway, probably the IR pathway.

Apart from disabling olfaction, Orco KO had dramatic consequences in homozygous moths. First, the homozygous ratio at G2 was lower than the expected Mendelian ratio, including in the 6 bp mutants. We suggested upper that the two lacking amino acids of these 6bp mutants may be located in an Orco domain important for the detection of some molecules since the detection of the pheromone was slightly modied. It is possible that this Orco modication may alter the detection of some other key odorants involved in food detection, leading to an increase in larval death. Second, we could not maintain homozygotes since couples gave only few sparse eggs that never hatched. It is possible that Orco KO impaired sperm activity, as demonstrated in mosquitoes40.

We cannot exclude impaired oogenesis and spermatogenesis, which were not investigated in the homozygotes, but we suspected that Orco KO-induced default in pheromone detection would impact mate detection. We could conrm this latter hypothesis since removing antennae in wild-type males was sufficient to impair mating, as already demonstrated in other moth species34,35. We cannot exclude that default in homozygote mating is due to o-target eects that the CRISPR/Cas9 system frequently induces4143. Scanning the genome searching for possible additional sequences matching the gRNA was not possible since S. littoralis is a non-model species with no genome yet available. However, since the 6bp deletion mutation (that did not modify the Orco reading frame) did not resume the olfactory responses in homozygous mutants and did not impair their reproduction, we are condent that the induced anosmia in the Orco KO mutants is specic.

The technological resource for noctuids described here unlocks their potential as future genetic model systems. This study appears as a proof of concept that would be potentially applicable with little adjustment to a wide variety of related species. With the development of sequencing technologies, an increased number of Lepidoptera genomes will be soon available, including those from crop pest species, thus such reverse genetic tools are highly anticipated. Modication of genomic fragments driven by pairs of double-strand breaks can be used to test the function of the targeted gene, as demonstrated here. The efficiency of this technology opens new routes since it not only allows targeting specic exons, but also putative cis-regulatory elements. This technic development may be also used for knock-in, allowing transgene replacement or even gene drive44, that may lead to the long term development of novel pest control technologies45.

Methods

SlitOrco Genomic DNA was prepared from S. littoralis larvae using the Wizard

Genomic DNA Purication Kit (Promega, Madison WI, USA) and used as template for PCRs conducted with primers designed to amplify putative intron-exon boundaries (Supplementary Table S1)

Three RNA guides were designed against exons 2 and 4 of the Orco gene (Fig.1) using the CRISPOR gRNA design tool cripsor.tefor.net and the SlitOrco genomic DNA sequence as target: cr104 (5-GGCCATGTTGATACCCATAC-3), cr105 (5-GGTGTGAGTGAAGAAGAGGA-3) and cr106 (5-GGCCTTCAGAGGCAGTCGAG-3). Guide sequences were subcloned in DR274 (http://www.addgene. org/42250) derived vector. Plasmids were digested by DraI, puried and transcribed using Hiscribe T7 high yield

SCIENTIFIC REPORTS

6

www.nature.com/scientificreports/

transcription kit (New England Biolabs). Reactions were puried using EZNA microelute RNA clean-up kit (OMEGA Biotek). Streptococcus pyogenes Cas9 protein, bearing 3 nuclear localization sequences, was produced in Escherichia coli and puried as previously described46.

S. littoralis was reared in the laboratory at 24C, 70% relative humidity on semi-articial diet47 under light:dark 8:16 photophase. To get eggs, boxes of two to three days-old females and two males were prepared and followed up every half hour. Freshly laid egg batches were collected, scales were removed and the mix Cas9 (12.3 M in 20 mM Hepes-NaOH pH 7.5, 150 mM KCl, 1mM TCEP, 10% Glycerol) and gRNA (38.8M in water) was injected without any further dilution through the chorion using a pulled borosilicate glass capillary attached to an Eppendorf - Transjector 5246, within one hour aer egg-laying. Eggs were let to develop and young larvae were reared as described above.

Genomic DNA was isolated from batches of young larvae (gRNA efficiency evaluation) or from one pseudopod cut from fourth instar larvae (individual non-invasive genotyping) with the Wizard Genomic

DNA Purication kit (Promega) and used as template for PCR using specic primers (for cr104 and cr105 gRNA: forward primer TGCCCAATTGATAAGCTCCT and reverse primer AATGCAACTCACCCCAGTTC; for cr106 gRNA: forward primer TGCGATCCAGTTCACTTTGA, reverse primer GCACTTGTAACGCCTTTGGT) amplifying a DNA fragment encompassing the target Orco sequence (Fig.1). Mutagenic events were detected with the T7EI (New England Biolabs, Ipswich, MA USA) assay as previously described48 (Fig.2). Mutated samples were sequenced at all generations (Biodal, Vaulx-en-Velin, France). Individuals carrying KO mutations were backcrossed with wild-type adults. The G1 progeny was screened for targeted mutations as in G0, using both T7EI assay and sequencing. Male and female G1 heterozygotes carrying the same mutation were then coupled together to generate homozygous Orco G2 mutants. G2 were genotyped as G1 except that a wild-type Orco fragment was included in the PCR template to reveal homozygous mutations, and sequenced for verication.

EAG was performed on one isolated antenna from each G2 wild-type, heterozygous or homozygous adults, to evaluate the Orco mutation eect on the antennal capacity to detect specic odor-ants. Mounted between two glass electrodes containing Ringer solution49 and continuously humidified by charcoal-ltered airow (70L/h), antennae were stimulated, with a one minute interval, with 10g of E-ocimene, Z3-hexenyl-acetate, E2-hexenol, benzyl alcohol and 2-phenyl acetaldehyde diluted in paraffin oil, 10g of propionic acid in water (all purchased from Sigma-Aldrich) and 1g of the main sex pheromone component Z9,E11-14:Ac (Gi from Martine Lettere, Versailles, France) in hexane (Carlo Erba) (pus of 500 ms, 10 L/h). EAG amplitudes were calculated by subtracting the solvent response. Analyses were done with Clampt 10 soware (Molecular Devices).

Single pair matings were organized as follows: 1) male and female having both their antennae (positive control), 2) both male and female with no antennae, 3) male with no antennae and intact female and 4) intact male and female with no antennae. Six pairs of each of the three categories were kept with sugar water under the same conditions as the regular rearing and observed every two hours during three scotophases for mating. Egg laying was recorded aer the 3 days and eggs were kept until larvae hatching, if any.

A t-test was used for statistical comparison of the EAG results.

1. Terenius, O. et al. RNA interference in Lepidoptera: An overview of successful and unsuccessful studies and implications for experimental design. Journal of Insect Physiology 57, 231245 (2011).

2. Tamura, T. et al. Germline transformation of the silkworm Bombyx mori L-using a piggyBac transposon-derived vector. Nature Biotechnology 18, 8184 (2000).

3. Carroll, D. In Annual Review of Biochemistry Vol. 83 (ed Kornberg, R. D.) 409439 (2014).4. Takasu, Y. et al. Targeted mutagenesis in the silkworm Bombyx mori using zinc nger nuclease mRNA injection. Insect biochemistry and molecular biology 40, 759765 (2010).

5. Merlin, C., Beaver, L. E., Taylor, O. R., Wolfe, S. A. & Reppert, S. M. Efficient targeted mutagenesis in the monarch buttery using zinc-nger nucleases. Genome Research 23, 159168 (2013).

6. Ma, S. et al. Highly efficient and specic genome editing in silkworm using custom TALENs. Plos one 7, e45035 (2012).7. Sajwan, S. et al. Efficient disruption of endogenous Bombyx gene by TAL eector nucleases. Insect Biochem Mol Biol 43, 1723 (2013).

8. Sakurai, T. et al. Targeted disruption of a single sex pheromone receptor gene completely abolishes in vivo pheromone response in the silkmoth. Scientic reports 5, 11001 (2015).

9. Ma, S. Y. et al. CRISPR/Cas9 mediated multiplex genome editing and heritable mutagenesis of BmKu70 in Bombyx mori. Scientic reports 4, 4489 (2014).

10. Wang, Y. Q. et al. The CRISPR/Cas System mediates efficient genome engineering in Bombyx mori. Cell Research 23, 14141416 (2013).

11. Li, X. et al. Outbred genome sequencing and CRISPR/Cas9 gene editing in butteries. Nature communications 6, 82128212 (2015).12. Markert, M. J. et al. Genomic Access to Monarch Migration Using TALEN and CRISPR/Cas9-Mediated Targeted Mutagenesis. G3 (Bethesda), doi: 10.1534/g3.116.027029 (2016).

13. Doudna, J. A. & Charpentier, E. The new frontier of genome engineering with CRISPR-Cas9. Science 346, doi: 125809610.1126/ science.1258096 (2014).

14. Hsu, P. D., Lander, E. S. & Zhang, F. Development and applications of CRISPR-Cas9 for genome engineering. Cell 157, 12621278 (2014).

15. Benton, R., Sachse, S., Michnick, S. W. & Vosshall, L. B. Atypical membrane topology and heteromeric function of Drosophila odorant receptors in vivo. PLos Biology 4, 240257 (2006).

SCIENTIFIC REPORTS

7

www.nature.com/scientificreports/

16. Larsson, M. C. et al. Or83b encodes a broadly expressed odorant receptor essential for Drosophila olfaction. Neuron 43, 703714 (2004).

17. Sato, K. et al. Insect olfactory receptors are heteromeric ligand-gated ion channels. Nature 452, 1002U1009 (2008).18. Smart, R. et al. Drosophila odorant receptors are novel seven transmembrane domain proteins that can signal independently of heterotrimeric G proteins. Insect biochemistry and molecular biology 38, 770780 (2008).

19. Wicher, D. et al. Drosophila odorant receptors are both ligand-gated and cyclic-nucleotide-activated cation channels. Nature 452, 10071011 (2008).

20. Abuin, L. et al. Functional architecture of olfactory ionotropic glutamate receptors. Neuron 69, 4460 (2011).21. Ai, M. et al. Acid sensing by the Drosophila olfactory system. Nature 468, 691U112 (2010).22. DeGennaro, M. et al. orco mutant mosquitoes lose strong preference for humans and are not repelled by volatile DEET. Nature 498, 487491 (2013).

23. Liu, C. et al. Distinct olfactory signaling mechanisms in the malaria vector mosquito Anopheles gambiae. PLos Biology 8, 117 (2010).

24. Franco, T. A., Oliveira, D. S., Moreira, M. F., Leal, W. S. & Melo, A. C. A. Silencing the odorant receptor co-receptor RproOrco aects the physiology and behavior of the Chagas disease vector Rhodnius prolixus. Insect biochemistry and molecular biology, 19 (2015).

25. Zhou, Y.-L. et al. Silencing in Apolygus lucorum of the olfactory coreceptor Orco gene by RNA interference induces EAG response declining to two putative semiochemicals. Journal of Insect Physiology 60, 3139 (2014).

26. Zhao, Y. Y., Liu, F., Yang, G. & You, M. S. PsOr1, a potential target for RNA interference-based pest management. Insect Molecular Biology 20, 97104 (2011).

27. Jones, W. D., Nguyen, T. A. T., Kloss, B., Lee, K. J. & Vosshall, L. B. Functional conservation of an insect odorant receptor gene across 250 million years of evolution. Current Biology 15, R119R121 (2005).

28. Kobayashi, Y., Tanaka, M. & Ando, H. In Lepidoptera, moths and buttereflies: Morphology, physiology, and development (ed Kristensen, N. P.) Ch. 19, (De Gruyter, 2003).

29. Binyameen, M. et al. Identication of plant semiochemicals and characterization of new olfactory sensory neuron types in a polyphagous pest moth, Spodoptera littoralis. Chemical Senses 39, 719733 (2014).

30. Malo, E. A., Renou, M. & Guerrero, A. Analytical studies of Spodoptera littoralis sex pheromone components by electroantennography and coupled gas chromatography-electroantennographic detection. Talanta 52, 525532 (2000).

31. Nakagawa, T., Sakurai, T., Nishioka, T. & Touhara, K. Insect sex-pheromone signals mediated by specic combinations of olfactory receptors. Science 307, 16381642 (2005).

32. Ljungberg, H., Anderson, P. & Hansson, B. S. Physiology and morphology of pheromone-specic sensilla on the antennae of male and female Spodoptera littoralis (Lepidoptera: Noctuidae). Journal of Insect Physiology 39, 253260 (1993).

33. Hopf, T. A. et al. Amino acid coevolution reveals three-dimensional structure and functional domains of insect odorant receptors. Nat Commun 6, 6077 (2015).

34. Charlton, R. E. & Carde, R. T. Factors mediating copulatory behavior and close-range mate recognition in the male gypsy moth, Lymantria dispar (L.). Canadian Journal of Zoology 68, 19952004 (1990).

35. Pathak, P. H. & Krishna, S. S. Variation in the reproductive capacity of Earias vittella (F.) (Lepidoptera: Noctuidae) following antennectomy or alactomy in males or wing loss in both sexes. Proceedings of the Indian Academy of Sciences-Animal Sciences 95, 613616 (1986).

36. Bae, S., Kweon, J., Kim, H. S. & Kim, J.-S. Microhomology-based choice of Cas9 nuclease target sites. Nature Methods 11, 705706 (2014).

37. Krieger, J., Grosse-Wilde, E., Gohl, T. & Breer, H. Candidate pheromone receptors of the silkmoth Bombyx mori. European Journal of Neuroscience 21, 21672176 (2005).

38. Croset, V. et al. Ancient protostome origin of chemosensory ionotropic glutamate receptors and the evolution of insect taste and olfaction. PLos Genetics 6, e1001064 (2010).

39. Olivier, V., Monsempes, C., Francois, M. C., Poivet, E. & Jacquin-Joly, E. Candidate chemosensory ionotropic receptors in a Lepidoptera. Insect Molecular Biology 20, 189199 (2011).

40. Pitts, R. J., Liu, C., Zhou, X. F., Malpartida, J. C. & Zwiebel, L. J. Odorant receptor-mediated sperm activation in disease vector mosquitoes. Proceedings of the National Academy of Sciences of the United States of America 111, 25662571 (2014).

41. Daimon, T., Kiuchi, T. & Takasu, Y. Recent progress in genome engineering techniques in the silkworm, Bombyx mori. Development Growth & Dierentiation 56, 1425 (2014).

42. Lin, S. C., Chang, Y. Y. & Chan, C. C. Strategies for gene disruption in Drosophila. Cell and Bioscience 4, 63, doi: 10.1186/2045-3701-4-63 (2014).

43. Pattanayak, V. et al. High-throughput proling of o-target DNA cleavage reveals RNA-programmed Cas9 nuclease specicity. Nature Biotechnology 31, 839843 (2013).

44. Gantz, V. M. & Bier, E. The mutagenic chain reaction: A method for converting heterozygous to homozygous mutations. Science 348, 442444 (2015).

45. Webber, B. L., Raghu, S. & Edwards, O. R. Opinion: Is CRISPR-based gene drive a biocontrol silver bullet or global conservation threat? Proceedings of the National Academy of Sciences of the United States of America 112, 1056510567 (2015).

46. Menoret, S. et al. Homology-directed repair in rodent zygotes using Cas9 and TALEN engineered proteins. Scientic reports 5, 14410 (2015).

47. Poitout, S., Bues, R. & Lerumeur, C. Rearing on a simple articial medium 2 Noctuid parasites of cottonEarias insulana and Spodoptera littoralis. Entomologia Experimentalis et Applicata 15, 341350 (1972).

48. Auer, T. O., Duroure, K., De Cian, A., Concordet, J.-P. & Del Bene, F. Highly efficient CRISPR/Cas9-mediated knock-in in zebrash by homology-independent DNA repair. Genome Research 24, 142153 (2014).

49. Kaissling, K. E. & Thorson, J. In Receptor for neurotransmitters, hormones and pheromones in insects (eds Sattelle, D. B., Hall, L. M. & Hildebrand, J. G.) 261282 (Elsevier/North-Holland Biomedical Press, 1980).

We thank Christine Merlin (Texas A&M, USA) for precious advices in egg injection and Aurlie Janvier (iEES-Paris, France) for help in genotyping. This work was supported by an INRA grant (Plant Health and Environment department) and an INRA fellowship to FK. We also thank the EFOR network and TEFOR (Investissement d'avenir-ANR-II-INBS-0014) for promoting our model and innovative transgenesis tools.

F.A.K. carried out egg injections, participated in genotyping, in electrophysiology experiments and in draing the manuscript. C.M. participated in genotyping, electrophysiology experiments and statistical analysis. M.-C.F. participated in insect rearing and genotyping. A.d.C. carried out RNA guide and Cas9 synthesis and participated in the study conception. C.R. participated in egg injection and helped to dra the manuscript. J.-P.C. participated in the study conception, carried out RNA guide design and helped to dra manuscript. E.J.-J. conceived and

SCIENTIFIC REPORTS

8

www.nature.com/scientificreports/

coordinated the study and participated in draing the manuscript. All authors read and approved the nal manuscript.

Supplementary information accompanies this paper at http://www.nature.com/srep

Competing nancial interests: The authors declare no competing nancial interests.

How to cite this article: Koutroumpa, F. A. et al. Heritable genome editing with CRISPR/Cas9 induces anosmia in a crop pest moth. Sci. Rep. 6, 29620; doi: 10.1038/srep29620 (2016).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

SCIENTIFIC REPORTS

9