Next Article in Journal
Silver: Forging almost Gold Standard Datasets
Previous Article in Journal
Characterization of a Missense Mutation in the Catalytic Domain and a Splicing Mutation of Coagulation Factor X Compound Heterozygous in a Chinese Pedigree
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 12
Issue 10
10.3390/genes12101522
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Knockouts Indicate That the ABCC2 Protein in the Bollworm Helicoverpa zea Is Not a Major Receptor for the Cry1Ac Insecticidal Protein

by
Omaththage P. Perera
1,*,
Nathan S. Little
1,
Heba Abdelgaffar
2,
Juan Luis Jurat-Fuentes
2 and
Gadi V. P. Reddy
1
1
Southern Insect Management Research Unit, USDA, Agricultural Research Service, Stoneville, MS 38776, USA
2
Department of Entomology and Plant Pathology, University of Tennessee, Knoxville, TN 37996, USA
*
Author to whom correspondence should be addressed.
Genes 2021, 12(10), 1522; https://doi.org/10.3390/genes12101522
Submission received: 1 September 2021 / Revised: 24 September 2021 / Accepted: 26 September 2021 / Published: 28 September 2021
(This article belongs to the Section Animal Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
Members of the insect ATP binding cassette transporter subfamily C2 (ABCC2) in several moth species are known as receptors for the Cry1Ac insecticidal protein from Bacillus thuringiensis (Bt). Mutations that abolish the functional domains of ABCC2 are known to cause resistance to Cry1Ac, although the reported levels of resistance vary widely depending on insect species. In this study, the function of the ABCC2 gene as a putative Cry1Ac receptor in Helicoverpa zea, a major pest of over 300 crops, was evaluated using CRISPR/Cas9 to progressively eliminate different functional ABCC2 domains. Results from bioassays with edited insect lines support that mutations in ABCC2 were associated with Cry1Ac resistance ratios (RR) ranging from 7.3- to 39.8-fold. No significant differences in susceptibility to Cry1Ac were detected between H. zea with partial or complete ABCC2 knockout, although the highest levels of tolerance were observed when knocking out half of ABCC2. Based on >500–1000-fold RRs reported in similar studies for closely related moth species, the low RRs observed in H. zea knockouts support that ABCC2 is not a major Cry1Ac receptor in this insect.

1. Introduction

The entomopathogenic bacterium Bacillus thuringiensis (Bt) produces Cry insecticidal proteins which are the main active components of sprays and are expressed by transgenic crop plants for environmentally sensitive pest control. The cry1Ac gene was first commercialized in the USA in 1996 as a plant incorporated protectant (PIP) in cotton [1], and since then has also been commercialized or authorized as PIP in corn (1997), eggplant in Bangladesh (2013) and the Philippines (2021), rice in China (2009) and the USA (2018), and soybean in several countries [2]. The main threat to the sustainable use of Cry1Ac as a PIP is the development of resistance in target pests. Practical resistance to Cry1Ac as a PIP in cotton has been reported in populations of bollworm (Helicoverpa zea) in the USA [3], and the pink bollworm (Pectinophora gossypiella) in India [4]. In P. gossypiella, this resistance was associated with alternative splicing which generated aberrant transcripts of cadherin, a receptor for Cry1Ac in that insect [5]. This cadherin that lacks a functional receptor is associated with reduced Cry1Ac binding and resistance [6]. Although the mechanistic description of resistance in H. zea remains elusive, mutations in a novel cadherin gene [7] and reduced Cry1Ac processing [8] have been proposed as putative causes.
Despite the existence of cadherin-based resistance to Cry1Ac, the most commonly reported mechanism for practical resistance to Cry protein PIPs in Lepidoptera involves alterations in ATP-binding cassette (ABC) transporter genes serving as receptors for these proteins on the midgut epithelium [9]. Among the ABC transporter superfamily, members of subfamily C2 (ABCC2) have been identified as functional receptors for Cry1Ac in Bombyx mori [10], Chloridea (Heliothis) virescens [11], Helicoverpa armigera [12], Plutella xylostella [13,14], Spodoptera frugiperda [15] and S. exigua [16,17]. Importantly, since these ABCC2 proteins are thought to facilitate oligomerization and insertion of the Cry1Ac protein into the plasma membrane [18], alterations in ABCC2 genes are linked with resistance to Cry1Ac [9]. The goal of this study was to use the CRISPR/Cas9 gene editing system in testing ABCC2 as a candidate Cry1Ac receptor in H. zea. In attaining this goal, we targeted different exons in the H. zea ABCC2 gene to progressively truncate the ABCC2 transporter protein and then tested the effects on susceptibility to Cry1Ac using diet bioassays.

2. Materials and Methods

2.1. ABCC2 Transporter Gene Sequence

The publicly available assembled H. zea genome [19] and a custom genome assembly from an H. zea female from the SIMRU laboratory colony were used to create a searchable database (BlastStation-Local 64 v1.3, TM Software, Inc., Arcadia, CA, USA). The custom genome was sequenced and assembled by a service provider (Hudson Alpha, Huntsville, AL, USA) using the Chromium Genome Solution and Supernova v1.1.5 genome assembly platforms (10× Genomics, Pleasanton, CA, USA). The database was interrogated with the H. armigera ABCC2 transporter cDNA sequence (GenBank accession number KF479231) using BlastStation-Local 64, and a genomic scaffold containing the ABCC2 gene was annotated and submitted to GenBank (accession number KY701524). The predicted ABCC2 gene transcript constructed from the exons identified in the genomic DNA sequence was aligned with ABCC2 mRNA sequences from GenBank to validate the annotation. The translated amino acid sequence was used to identify functional domains of the ABCC2 protein using the Simple Modular Architecture Research Tool (SMART: http://smart.embl-heidelberg.de, accessed on 16 August 2021 [20]). Transmembrane (TMD) and ATPase or nucleotide-binding (NBD) domains were identified and genomic exons corresponding to the TMD and NBD domains were selected as potential CRISPR guide RNA targets.

2.2. Guide RNA Design

All CRISPR RNAs (crRNA) used in the experiments were designed as described in Perera et al., [21] using the tools available at http://crispr.mit.edu/m (accessed on 15 October 2016) [22]. Briefly, target ABCC2 exon sequences were submitted to the design pipeline and the best crRNA design for each target location was purchased from IDT DNA Inc (Coralville, IA, USA) in the Alt-R™ crRNA format (Table 1). Companion reagents for the Alt-R™ crRNA, universal transactivating CRISPR RNA (Alt-R™ tracrRNA) and a modified Cas9 nuclease containing one N-terminal nuclear localization signal (NLS) and two C-terminal NLS (Cas9-3NLS), were also purchased from IDT DNA Inc.

2.3. Insect Colony Maintenance and Egg Collection

All methodology pertaining to H. zea colony maintenance, egg collection, crRNA design, and microinjection have been described in detail in Perera et al. [21]. The H. zea colony used is routinely maintained at the USDA-ARS Southern Insect Management Research Unit (SIMRU, Stoneville, MS, USA) following the methods described by Gore et al. [23]. Emerging adults were placed in 3.8 L paper containers with 3% sugar solution for 3–4 days to facilitate mating. Females were identified by examining the genitalia and were placed in a 5 × 5 mm wire mesh cage with a piece of wax paper wrapped around it. The wire cage was then placed inside a dark laboratory cabinet for 30 min until egg collection. The wax paper was removed, eggs were dislodged from the wax paper using a fine paint brush, and then mounted on 1 mm wide strips of double-sided tape (3M, St. Paul, MN, USA) attached to 25 × 25 mm plastic coverslips. The eggs were then dehydrated for 15 min by placing in a desiccator containing silica gel and approximately −100 KPa (−30.5 inches of mercury) vacuum before being used for microinjections.

2.4. Microinjections

In vitro assembled nucleoprotein complexes containing single guide RNA (sgRNA) and Cas9 were optimized previously in H. zea to obtain high rates of genome editing [21]. Using the same protocol, crRNA and tracrRNA were resuspended separately in the nuclease-free duplex buffer (30 mM Hepes pH 7.5, 100 mM potassium acetate) to yield a 100 µM solution. To form functional sgRNA, 1 µL (100 µM) each of reconstituted crRNA and tracrRNA and 8 µL of duplex buffer was mixed to yield a final concentration of 10 µM. When two or more crRNAs were used in a single injection mixture, the tracrRNA was increased to match the sum of all crRNA concentrations to yield a final concentration of 10 µM in a 10 µL solution. This mixture was heated to 95 °C for 5 min in a thermal cycler block and allowed to cool to 25 °C on the bench top. The reconstituted sgRNA mixture (4 µL) was transferred to a new tube where Cas9-3NLS nuclease was added to a concentration equivalent to the sum of the sgRNA concentrations, and the volume of the injection mixture was adjusted to 10 µL with 1 µL of 10× injection buffer (5 mM KCl; 0.1 M sodium phosphate, pH 6.8) and nuclease free water. The final concentration of the sgRNA-Cas9-3NLS nucleoprotein complex in the injection mixture was calculated to be 4 µM. If more than two sgRNAs were used together, the injection mixture was diluted to a final concentration of 2 to 3 µM using 1× injection buffer to prevent clogging of needles due to high concentration of sgRNA-Cas9-3NLS nucleoprotein complexes. Borosilicate microinjection needles prepared with Narishige model PC-10 needle puller (Narishige International USA, Inc., Amityville, NY, USA) and beveled with Sutter Instrument Company (Novato, CA, USA) Model BV-10 beveller were used in initial embryo injections. Microinjection needles were backfilled with approximately 3 µL of the injection mix and the eggs were injected using a Narishige micromanipulator model MMN-333 and an Olympus SZ stereo microscope (Olympus Corporation, Waltham, MA, USA) equipped with a mechanical stage. Each egg was injected with approximately 5 nL of injection mix. Control egg hatch rates were calculated using the eggs processed up to the desiccation stage, but without subjecting them to injection.
Initial microinjection experiments (Group 1) were performed with guide RNA designed to exons 21, 22 and 24 that targeted the second NBD domain (NBD2) in the ABCC2 protein. After evaluation of the success rate of the first ABCC2 gene editing effort, two additional rounds of injections were performed by multiplexing different combinations of sgRNA to target the first and second TM domains of the ABCC2 protein. The group 2A sgRNAs targeted exons 3, 8, 13, and 21, while the group 2B sgRNAs targeted exons 13, 19, and 21.
Coverslips containing injected eggs were placed in a 100 mm plastic petri dish layered with a moist filter paper, covered with a lid, sealed with plastic tape, and placed in a plastic box designated as secondary containment. Eggs were incubated at room temperature for approximately 4 hours prior to transferring to a designated incubator set to 26 °C and relative humidity over 75%. Eggs were observed for larval hatch and the neonates were transferred to diet cups containing Helicoverpa diet [23]. The total number of eggs injected, and the number of eggs hatched were recorded for each experiment.

2.5. Genotyping and DNA Analysis

Larvae collected from injected embryos (P0) were raised to fourth instar and chilled on ice for 30 min prior to obtaining 5–15 µL of hemolymph by cutting the tip of one of the medial (abdominal) prolegs with a pair of micro scissors. Droplets of hemolymph were collected using a pipette and transferred to a tube containing Tissue and Cell Lysis buffer from the MasterPure DNA purification kit (Epicentre Technologies, Madison, WI, USA) supplemented with 20 mg/mL RNAse. Genomic DNA was extracted from hemolymph following the manufacturer’s protocol and re-suspended in 20 µL of 10 mM Tris-HCl. Forward and reverse primers designed for exon 2 (3800Hz_ABCC2_Ex2F) and exon 24 (3809Hz_ABCC2_Ex24R) of the ABCC2 gene were used to amplify by PCR the ABCC2 gene using the LongAmp Taq polymerase mix (New England Biolabs, Ipswich, MA, USA) on a PTC-200 DNA Engine (Supplementary data Table S2: Primers). Amplifications were performed with a thermal cycling profile containing 30 second initial denaturation at 95 °C, followed by 10 cycles of 10 sec denaturation at 95 °C, 15 sec annealing at 52 °C, and 8 min 30 sec extension at 65 °C. At the end of the 11th cycle, 25 additional amplification cycles were performed by increasing the 65 °C extension by 10 sec in each successive cycle. The final extension time at the 36th cycle was 12 min for 40 sec. The PCR buffer contained 60 mM Tris-SO4, 20 mM ammonium sulfate, 3% glycerol, 0.06% IGEPAL® CA-630 (octylphenoxypolyethoxyethanol), 0.05% Tween® 20, and 2.5 mM MgSO4 (pH 9.0). Amplicons from each larva were cleaned by binding to AmPure XP paramagnetic beads (Beckman Coulter, Indianapolis, IN, USA) at DNA:Beads ratio of 1:1.8. Nucleotide sequences of the purified amplicons were obtained by direct sequencing with Sanger dideoxy sequencing using primers designed to flank target sites. Primers used for PCR amplification and sequencing of the ABCC2 gene are listed in Supplementary data Table S1.

2.6. Insect Husbandry

Insects were reared at 28 °C and 14:10 light: dark cycles in environmental chambers (Percival Scientific, Perry, IA, USA). Insects from P0 with mutations at the target sites were selected for mating as single pairs with wild type insects of the opposite sex from the SIMRU laboratory colony. Resulting heterozygotes (F1) from each single-pair mating were inbred to obtain F2 progeny. Fourth instar F2 larvae were genotyped by sequencing the DNA extracted from a droplet of hemolymph as described previously and the number of the homozygous mutant, heterozygotes, and homozygous wild type insects were determined. Group matings were set up among homozygous mutants, and as a recovery option in case the homozygous mutant insects did not reproduce, matings were also set between heterozygotes and homozygous mutants, and among heterozygotes. Homozygous insects from each knockout line were inbred for at least four generations, genotyped to verify mutations and used in bioassays to evaluate tolerance to Cry1Ac protein.

2.7. Validation of Gene Knockouts by mRNA and Genomic DNA Sequencing

Total RNA extracted from homozygous larvae from each knockout line was treated with DNAse I (NEB) at a concentration of 10 µg/mL for 30 min at 37 °C to remove any carryover genomic DNA. The DNAse was inactivated by heating to 90 °C for 10 min and the cDNA was synthesized using the SuperScript cDNA synthesis kit (Invitrogen, Carlsbad, CA, USA) with an anchored oligo dT primer (5’-GGTAATACGACTCACTATAGGGAGAAGAGGCGAGCACAGAATTAATACGACTTTTTTTTTTTTTTTTTTTTV-3’). Amplification by PCR of the ABCC2 cDNA was performed using primers 3800 Hz_ABCC2_Ex2F and 3809 Hz_ABCC2_Ex24R. Genomic DNA was amplified using primer pairs as outlined in the Table S1 and purified using the QiaX II gel extraction kit (Qiagen, Germantown, MD, USA). Amplicons were submitted to the USDA-ARS Genomics and Bioinformatics Research Unit (Stoneville, MS, USA) for Sanger dideoxy sequencing. Nucleotide sequence analyses were carried out using the Vector NTI Advance v11.5 suite (Invitrogen). Additional sequencing of cDNA and genomic DNA amplicons from wild type and knockout lines was performed using the Oxford Nanopore Minion system (Oxford Nanopore, Oxford, UK) to obtain long reads, and resulting sequences were assembled using the DNAStar NGen module of DNAStar Lasergene software (DNAStar, Madison, WI, USA). The final curated sequences were submitted to GenBank.

2.8. Off-Target Sequence Analysis

The 12-nucleotide seed sequence of each sgRNA (i.e. the 12 nucleotides immediately upstream of PAM) [24] was used to search a local database created from the published H. zea genome [19: GCA_002150865.1] using BlastStation-Local64 v1.3 to identify genomic scaffolds with potential off-target sites. Off-target sequences with 0 to 3 mismatches in the 12-nucleotide seed sequence and containing a 5’-NGG-3’ PAM sequence [22,24,25] were selected for further review. Primers were developed to amplify and sequence off-targets with 0 to 2 mismatches in the seed sequence. Amplicons were sequenced on an ABI 3730Xl instrument at the USDA ARS Genomics and Bioinformatics Research Unit (Stoneville, MS, USA). Nucleotide sequences with potential off-target sites were analyzed with the Vector NTI 11.5 Advance suite (Invitrogen).

2.9. Bioassays

The Bacillus thuringiensis subsp. kurstaki HD-73 strain obtained from the Bacillus Genetic Stock Center (Columbus, OH, USA) was used to produce full length Cry1Ac protoxin as described in Little et al. [26]. Bioassays were performed in 128-well plastic trays (B-D International, Franklin Lakes, NJ, USA) with Cry1Ac protein incorporated in meridic diet. Each replicate consisted of 16 bioassay wells containing 0.5 mL of diet supplemented with 1, 10, 100, 500 or 1000 ppm of purified full length Cry1Ac protoxin (approx. 120 kDa). In untreated controls Cry1Ac was replaced with water in diet. A single neonate was placed per well before sealing with a ventilated plastic cover (B-D international, Franklin Lakes, NJ, USA). Trays were incubated at 28 °C for seven days, when stunting and mortality were recorded. Three bioassay replicates were performed for each knockout line, and the SIMRU laboratory colony of H. zea was used as a control. Mortality data were evaluated using regression analyses with the PROC PROBIT procedure in SAS v9.4 [27,28] to estimate the concentration that killed 50% of the tested population (LC50) and the fiducial limits. Resistance ratios were calculated by dividing the LC50 of the knockout and control strains.

2.10. Brush Border Membrane Vesicle (BBMV) Preparation

Midguts of fourth-instar H. zea from the SIMRU wild type and the edited SPM-8 and SPM-B1 strains were dissected and used to prepare BBMV using a differential centrifugation method [29], with minor modifications [30]. Isolated BBMV proteins were quantified in a Qubit fluorometer (Invitrogen), and then kept at −80 °C until used (less than two months). Brush border enzyme enrichment in the BBMV preparations was tested measuring specific activity of aminopeptidase-N (APN) using leucine-ρ-nitroanilide as substrate, as described elsewhere [31]. APN activities in the final BBMV preparations were enriched 3 to 5-fold when compared to initial midgut homogenates.

2.11. Cry1Ac Radioiodination and Binding Assays

The Cry1Ac protoxin was activated by incubation with bovine trypsin and then purified using anion exchange chromatography, as described previously [32]. Activated Cry1Ac (25 μg) was radiolabeled with 0.5 mCi of NaI125 (Perkin Elmer, Boston, MA, USA) using chloramine T, as previously described [33]. The labeled Cry1Ac was purified from free iodine using a PD-10 desalting column (GE Healthcare Life Sciences, Piscataway, NJ, USA) equilibrated in column buffer (20 mM Tris-HCl, 150 mM NaCl, 0.1% BSA, pH 8.6). Fractions containing the radiolabeled Cry1Ac protein based on autoradiography of fractions resolved by SDS-8% PAGE, were pooled and used for binding assays. The specific activity of the labeled Cry1Ac protein was 5.12 µCi/µg, based on input protein.
Binding saturation assays were performed with a constant BBMV protein concentration (0.15 μg/µL) and increasing concentrations of 125I-Cry1Ac (from 0.1 to 5 nM) as ligand in a 0.1 mL final reaction volume of binding buffer (PBS, pH 7.4, 0.1% BSA). Reactions were incubated at room temperature for 1 h and then stopped by centrifugation (14,500× g for 10 min), and pellets washed with 0.5 mL of ice-cold binding buffer. After the second centrifugation, the radioactivity in the final pellets was measured in a Wizard2 gamma counter (Perkin Elmer). Non-specific binding was determined by including an excess (300-fold the amount of ligand) of unlabeled Cry1Ac in the binding reactions. Specific binding was calculated by subtracting non-specific from the total (in the absence of competitor) binding. Data from two replicated experiments performed in duplicate were fitted to a one binding site model as the best-fitting model in the SigmaPlot v.11.0 software (Systat Software, San Jose, CA, USA) to obtain the apparent dissociation constant (Kd) and concentration of binding sites (Bmax).

3. Results

3.1. Genomic Organization of ABCC2 Transporter Gene in Helicoverpa zea

Searches with the H. armigera ABCC2 mRNA sequence (GenBank accession KF479231) identified genomic scaffold KZ118297.1 (188,800 bp in length) in the public H. zea genome [19] and scaffold 4812 (190,673 bp) in a custom H. zea genome assembly. Gaps in scaffold KZ118297.1 spanned some exons, while the gaps in genomic scaffold 4812 were limited to introns. Therefore, genomic scaffold 4812 was annotated and submitted to GenBank (accession number KY701524). Genomic scaffold 4812 contained genes similar to ribosomal protein L8 of Bombyx mori (accession number NP_001037141.1), cryptochrome CRY2-like protein of H. armigera (AGA11662.1), golgin subfamily A member 4-like protein of Amyelois transitella (XP_013199561.1), decapping and exoribonuclease-like protein of Papilio machaon (XP_014360359.1), multidrug resistance-associated protein 4-like of H. armigera (XP_021196505.1), SMC domain protein of Heliothis subflexa (ADH16742.1), and ABCC3 protein of H. armigera (AHL68987.1) The organization of the ABCC2, ABCC3, and SMC domain protein genes was highly conserved between H. zea and H. subflexa (GenBank accession number GQ332573.1). The ABCC2 gene of H. zea consisted of 25 exons from the putative promoter to the polyadenylation site, spanning nucleotides 97,034 to 110,370 of scaffold 4812. However, only the first 24 exons coded for the ABCC2 transporter polypeptide. Functional domains of the H. zea ABCC2 protein predicted using online SMART are shown in Table S2.

3.2. Embryo Injections Egg Hatch and Mutation Rates

The initial round of embryo microinjections targeting exons 21, 22, and 24 of the ABCC2 gene produced a 42.5% hatch rate (74 larvae from 174 injected eggs). Considering that the highly inbred laboratory colony has an average hatch rate of 58.3%, the adjusted hatch rate [34] of injected embryos approximated 72.9%. Of the 74 larvae, 69 pupated and 61 adults emerged. Genetic analysis revealed that 56 adults had mutations in the ABCC2 gene. Most mutations in the ABCC2 gene identified in this round of injections were frame-shift deletions ranging from 1 to 19 bp at the target site in exon 21. Of the single pair matings with injected P0 adults, only 18 produced progeny. Screening of inbred F2 progeny from resulting families identified insects from several families homozygous for ABCC2 gene mutations in exons 21 and 22. One family (SPM8) with a 7 bp frame shift deletion (Supplementary, Figure S1) that truncated the NBD2 was selected to establish a colony for further evaluation. SPM-8 also had a 6 bp insertion at the sgRNA target site in exon 22. Another mutant family, SPM16, with a 426 bp deletion in the genomic DNA starting at the sgRNA target site in exon 21 that eliminated the exon 21 downstream of the target site along with a part of intron 21 and a 4 bp insertion at the exon 22 target site was also selected for further evaluation (Supplementary, Figure S2).
Subsequent embryo injections targeting TMD1 and TMD2 with multiplexed sgRNAs of groups 2A and 2B yielded 22.8% and 18.5% hatch rates, respectively, with a combined average hatch rate of 20.7%. The hatch rate of the uninjected control eggs collected from the same batch of moths was 40.6% (69 out of 170). The corrected average [34] hatch rate for this set of injected eggs was 50.9%. A total of 10 larvae from injected embryos (P0) died during pupation due to inability to complete transformation. In addition, 20 pupae also failed to develop into adults. Of the 82 adults that emerged, seven from each sex had severe wing and body deformities that precluded successful mating. Single pair matings of 37 females and 31 males of the P0 with wild type insects from the SIMRU laboratory colony gave rise to 30 and 20 F1 families, respectively. The results of embryo injections are summarized in Table S3. After evaluating all nucleotide sequence data, families SPM-A28C and SPM-B1 were selected to represent deletions of TMD1 and TMD2 domains, respectively. A large deletion of 4686 bp spanning from sgRNA target sites in exon 3 to exon 13 was identified in the knockout line SPM-A28C. This knockout line also had a 20 bp insertion at the exon 21 sgRNA target site (Supplementary, Figure S3). In knockout line SPM-B1, a deletion of 1896 bp from the sgRNA target site in exons 13 to intron 17 and a 10 bp deletion at the exon 21 target site were identified (Supplementary, Figure S4). Deletions and insertions caused by CRISPR/Cas9 at target sites of the knockout lines used in bioassays are shown in Figure 1. Sizes of genomic DNA amplicon from knockout lines and wild type insects are shown in Figure 2. True breeding lines for conducting bioassays were established by inbreeding the selected knockout lines.
Evaluation of cDNA sequences generated from mRNA isolated from ABCC2 gene knockout lines identified proteins that lack functional NBD2 (SPM8 and SPM-16), TMD2 (SPM-B1), or complete knockout of the ABCC2 protein (SPM-A28C). Alignment of cDNA sequences obtained from knockout lines with the reference full-length ABCC2 mRNA (KM360184) and putative amino acid sequences translated from the cDNA sequences are shown in Figures S5 and S6, respectively. Deletion of 7 bp in exon 21 in SPM-8 truncated the ABCC2 protein at amino acid 1128, eliminating most of the NBD2 and adding 17 random amino acids to the end of the truncated protein (Supplementary, Figure S6). Deletion of the latter half of exon 21 and a part of intron 21 in SPM-16 eliminated the splice donor site (GT) at the end of exon 21, leading to splicing of 49 nucleotides of intron 21 to exon 22. This deletion and the mis-splicing combination deleted 27 amino acids from the NBD2 of the ABCC2 protein (1127 to 1134 aa) while inserting 15 random amino acids into the protein. After the insertion of random amino acids, the reading frame of the exon 22 was restored but the 4 bp insertion at the exon 22 sgRNA target site caused a frameshift that inserted 19 random amino acids at the carboxy terminal end of the truncated ABCC2 protein, which rendered the NBD2 non-functional. The deletion in SMP-B1 between the sgRNA target site in exon 13 and intron 18 apparently led to use of an alternative splice donor site upstream of the exon 13 deletion to join 236 nucleotides of exon 13 to exon 19. This deletion of a part of exon 13 and all of exons 14 to 18 eliminated 255 amino acids forming the TMD2 (from 760 to 1013 aa) of the ABCC2 protein. The deletion of 10 bp at the exon 19 sgRNA target site also led to an alternative splicing event that eliminated 19 amino acids while inserting 10 random amino acids for a net deletion of nine amino acids. However, the amino acid sequence was restored to wild type sequence after amino acid position 1077 (Supplementary, Figure S6). In knockout line SPM-A28C, the ABCC2 protein was truncated after amino acid 127, leaving only the first transmembrane helix of the TMD1 intact (Supplementary, Figure S6 and Table S2). Nucleotide sequences from the genomic DNA and mRNA of the knockout lines were submitted to GenBank under accession numbers MZ913261 to MZ913265 and MZ954929 to MZ954932, respectively.

3.3. Off-Target Analysis

Searching the genome of H. zea with sgRNA sequences identified one potential off-target site for exon 13 sgRNA in scaffold KZ118710 and two and one potential off-target site for exon 19 sgRNA in genomic scaffolds KZ117297 and KZ117617, respectively (Table 2). No off-targets matching the search criteria were identified for the remaining sgRNA sequences. The exon 13 sgRNA was used for embryo injections generating lines SPM-A28C and SPM-B1, and the exon 19 sgRNA was used in generating knockout line SPM-B1. Therefore, genomic DNA from both SPM-A28C and SPM-B1 was amplified and sequenced using primers developed for the off-target sites in scaffold KZ118710 and for SPM-B1 also for off-target sites in scaffolds KZ117297 and KZ117617 (Supplementary, Table S1). Amplicons from wild type insects were used as the reference. Alignments of nucleotide sequences obtained from direct sequencing of amplicons from knockout lines SPM-A28C, SPM-B1, and wild type with reference genomic scaffold sequences indicated lack of mutations in the off-target sites evaluated (Supplementary, Figure S7).

3.4. Bioassays

Bioassays conducted with purified Cry1Ac protoxin indicated different levels of tolerance associated to distinct mutated ABCC2 domains (Table 3). SPM-8, which had the second ATPase domain (NBD2) inactivated, showed the lowest tolerance (RR = 7.3). The knockout line SPM-16 with NBD2 deletion and in-frame insertion was slightly more tolerant to Cry1Ac, yet the LC50 fiducial limits overlapped with SPM-8, supporting no significant differences. The SPM-B1 line, with both TMD 2 and NBD2 knocked out, had the highest resistance ratio (RR = 39.8). The LC50 fiducial limits for SPM-B1 did not overlap with those of SPM-8 but overlapped with the fiducial limits of SPM-16. Knockout line SPM-A28C, in which the ABCC2 transporter was almost completely knocked out, had an LC50 that was not significantly different from SPM-8 or SPM-16.

3.5. Saturation Cry1Ac Binding Assays

Based on results from bioassays, we selected the SPM-8 and SPM-B1 lines, as having the lowest and highest susceptibility among edited strains, for Cry1Ac-binding assays. In saturation assays with radiolabeled Cry1Ac and BBMV from SIMRU we detected high-affinity saturable binding (Figure 3 and Table 4). In contrast, specific Cry1Ac binding was dramatically reduced in BBMV from SPM-8 and SPM-B1 strains when compared to SIMRU (Figure 3). Estimation of binding parameters for SPM-8 and SPM-B1 detected a ≥10-fold reduction in the concentration of binding sites (Bmax) relative to SIMRU BBMV (Table 4). We were unable to obtain accurate estimates of affinity (Kd) for SPM-8 and SPM-B1 from the dataset (p in both cases >0.05).

4. Discussion

Members of the ABC transporter family C2 (ABCC2) are functional receptors for Cry1 proteins in lepidopteran larvae [35]. Mutations in ABCC2 genes have been genetically linked with high levels of resistance to Cry1Ac in the Heliothinae subfamily, including Heliothis virescens [11] and H. armigera [36]. However, not much is known about the role of ABCC2 as Cry1Ac receptor in the closely related H. zea. In addressing this knowledge gap, we obtained progressive deletions and a complete ABCC2 gene knockout using the CRISPR/Cas9 editing system, and then tested the putative Cry1Ac functional receptor role of ABCC2 protein in H. zea. Bioassays testing the edited H. zea lines detected moderately reduced susceptibility compared to the parental (unedited) line. Similarly, knockout of ABCC2 gene in the closely related species H. armigera resulted in 3.8-fold reduced susceptibility to Cry1Ac toxin, and high levels of resistance (>1500-fold) were only observed with additional knockout of ABCC3 [37]. Akin results were observed in P. xylostella strains with natural mutations or with CRISPR-edited single and double ABCC2 and ABCC3 knockouts [38]. More recently, higher levels of resistance (112-fold) to Cry1Ac protoxin were reported after ABCC2 gene knockout in H. armigera, and even higher levels of resistance (816-fold) were detected with an additional knockout of cadherin [39]. In contrast, knockout of ABCC2 gene in P. xylostella resulted in high resistance (724-fold) to Cry1Ac protoxin [14], suggesting that ABCC2 has distinct relevance to Cry1Ac intoxication depending on the species. In addition, we found that there was no association between the size of the ABCC2 deletion and susceptibility to Cry1Ac. For example, knockout of the NBD2 in lines SPM-8 and SPM-16 resulted in only marginal tolerance to Cry1Ac, which was not different from tolerance in the complete ABCC2 knockout line SPM-A28C. The SPM-B1 line, which had TM2 and NBD2 deleted, showed the highest tolerance to Cry1Ac (39.8-fold) compared to the wild type, but still was not significantly different from tolerance in SPM-16.
Previous reports identified residues in loop 1 of TMD1 as critical for Cry1Ac binding to ABCC2 in Spodoptera spp., and these residues are conserved in ABCC2 of H. zea. [40,41]. The only knockout line tested in the present study in which these residues would not be present in ABCC2 was the almost full-length knockout SPM-A28C. However, susceptibility to Cry1Ac in SPM-A28C was not different from susceptibility in the SPM-8 strain, which was edited to eliminate the NBD2. This observation suggests that ABCC2 in H. zea must be functionally active and that simple interactions with the extracellular loops of TMDs are not sufficient to support Cry1Ac toxicity. This hypothesis is not supported by results from mutant ABCC2 proteins lacking NBDs still functioning as Cry receptors [16,42]. Alternatively, it is also possible that truncated ABCC2 is not produced and present in the BBMV, as was reported for S. frugiperda with a mutational disruption of the sfABCC2 gene [43]. In support of this hypothesis, binding of Cry1Ac was dramatically reduced in the tested ABCC2 gene knockout strains, independently of the domain targeted or their level of susceptibility to Cry1Ac. This observation also suggests that the remaining low levels of Cry1Ac binding to alternative receptors other than ABCC2 in H. zea BBMV are still conducive to toxicity. An alternative previously reported Cry1Ac receptor in H. zea BBMV is the alkaline phosphatase ALP2 [44]. Since mutations in a novel cadherin gene are associated with field-evolved resistance of H. zea to Cry1Ac, this could also be an alternative receptor [7]. In fact, cadherin knockout in H. armigera resulted in higher resistance to Cry1Ac compared to ABCC2 knockout [37,39], which may be indicative of their relative importance for Cry1Ac toxicity in Helicoverpa spp.
In summary, we used synthetic sgRNAs and recombinant Cas9 nuclease to assemble nucleoproteins that were multiplexed to target multiple protospacer sites in the H. zea ABCC2 gene to efficiently recover desired knockout lines. Mutations that truncated the ABCC2 at TMD1, TMD2, or NBD2 resulted in low levels of tolerance to Cry1Ac, suggesting that by itself the ABCC2 protein is not a critical Cry1Ac receptor. Future work will focus on examining alternative Cry1Ac receptors and their putative interactions with ABCC2 during Cry1Ac intoxication, as described for cadherin in H. armigera [45].

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/genes12101522/s1, Figure S1: Alignment of Helicoverpa zea ATP binding cassette transporter subfamily C2 (ABCC2) genome sequence from 107,087 to 108,334 bp and amplicon sequence from knockout line SPM-8 with 7 bp deletion at the exon 21 sgRNA target site and 6 bp net insertion at the exon 22 sgRNA target site. Figure S2: Alignment of Helicoverpa zea ATP binding cassette transporter subfamily C2 (ABCC2) genome sequence from 107,089 to 108,333 bp and amplicon sequence from knockout line SPM-16 with a deletion of 426 bp and an insertion of four bp at the sgRNA target sites in exon 21 and 22, respectively. Figure S3: Alignment of Helicoverpa zea ATP binding cassette transporter subfamily C2 (ABCC2) reference genome sequence from 99,062 to 108,844 bp and amplicon sequences from wild type SIMRU and knockout line SPM-A28C. Figure S4: Alignment of Helicoverpa zea ATP binding cassette transporter subfamily C2 (ABCC2) reference genome sequence from 104,174 to 108,822 bp and amplicon sequence from the knockout line SPM-B1 with 1876 bp and 10 bp deletions at the sgRNA target sites of exon 13 and exon 19, respectively. Figure S5: Alignment of partial mRNA sequences of Helicoverpa zea ATP binding cassette transporter subfamily C2 (ABCC2) from knockout lines SPM-8, SPM-16, SPM-B1, and SPM-A28C with H. zea ABCC2 mRNA sequence KM360184 from GenBank. Figure S6: Alignment of putative amino acid sequences predicted from partial mRNA sequences of Helicoverpa zea ATP binding cassette transporter subfamily C2 (ABCC2) gene knockout lines SPM-8, SPM-16, SPM-B1, and SPM-A28C with the ABCC2 amino acid sequence AKH49600 from GenBank. Figure S7: Alignment of nucleotide sequences from potential off-target sites for exon 13 and exon 19 sgRNA in Helicoverpa zea ATP binding cassette transporter subfamily C2 (ABCC2) knockout lines SPM-A28C and SPM-B1. Table S1: Primer names and sequences used for PCR amplification and nucleotide sequencing of Helicoverpa zea ATP binding cassette transporter subfamily C2 (ABCC2) gene and off-target sites. Table S2: Functional domains predicted in the ATP binding cassette transporter subfamily C2 (ABCC2) protein of Helicoverpa zea predicted using online Simple Modular Architecture Research Tool (SMART: http://smart.embl-heidelberg.de accessed on 16 August 2021). Table S3: The number of eggs injected, percent larvae hatched, the number of males and females, and the total adults recovered from three rounds of embryo injections with sgRNA/Cas9 nucleoprotein complexes targeting ATP binding cassette transporter subfamily C2 (ABCC2) gene in Helicoverpa zea.

Author Contributions

Conceptualization, O.P.P.; methodology, O.P.P., N.S.L., H.A. and J.L.J.-F.; formal analysis, O.P.P., J.L.J.-F., and N.S.L.; investigation, O.P.P., N.S.L., H.A. and J.L.J.-F.; resources, O.P.P., and G.V.P.R.; data curation, O.P.P. and J.L.J.-F.; writing—original draft preparation, O.P.P.; writing—review and editing, O.P.P., N.S.L., J.L.J.-F., and G.V.P.R.; project administration, O.P.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the USDA Agricultural Research Service in house research project 6066-22000-091-00D. Partial support was also provided by the NC246 Multistate Hatch project and Agriculture and Food Research Initiative Foundational Program competitive grant No. 2018-67013-27820 from the USDA National Institute of Food and Agriculture (NIFA).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All pertinent nucleotide sequence data have been deposited in public databases and/or have been reported in the text, figures, or tables.

Acknowledgments

We thank Calvin Pierce (USDA ARS SIMRU) for assistance with genotyping, sequencing library construction, PCR, bioassays, and insect husbandry, Priya Chatakondi (SIMRU) for insect husbandry, Michelle Mullen (SIMRU) for bioassays and PROBIT analysis, Fanny Liu (USDA ARS Genomics and Bioinformatics Research Unit, Stoneville, MS, USA) for assistance with DNA sequencing, and Margaret Allen (USDA ARS Biological Control of Pests Research Unit, Stoneville, MS, USA) for providing equipment for preparing microinjection needles. The use or mention of a trademark or proprietary product does not constitute an endorsement, guarantee, or warranty of the product and does not imply its approval to the exclusion of other suitable products by the U.S. Department of Agriculture, an equal opportunity employer.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Perlak, F.J.; Oppenhuizen, M.; Gustafson, K.; Voth, R.; Sivasupramaniam, S.; Heering, D.; Carey, B.; Ihrig, R.A.; Roberts, J.K. Development and commercial use of Bollgard cotton in the USA—early promises versus today’s reality. Plant J. 2001, 27, 489–501. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. GM Approval Database [Internet]. 2021. Available online: http://www.isaaa.org/gmapprovaldatabase/ (accessed on 10 August 2021).
  3. Tabashnik, B.E.; Gassmann, A.J.; Crowder, D.W.; Carrière, Y. Insect resistance to Bt crops: Evidence versus theory. Nat. Biotechnol. 2008, 26, 199–202. [Google Scholar] [CrossRef]
  4. Dhurua, S.; Gujar, G.T. Field-evolved resistance to Bt toxin Cry1Ac in the pink bollworm, Pectinophora gossypiella (Saunders) (Lepidoptera: Gelechiidae), from India. Pest Manag. Sci. 2011, 67, 898–903. [Google Scholar] [CrossRef] [PubMed]
  5. Fabrick, J.A.; Ponnuraj, J.; Singh, A.; Tanwar, R.K.; Unnithan, G.C.; Yelich, A.J.; Li, X.; Carriere, Y.; Tabashnik, B.E. Alternative splicing and highly variable cadherin transcripts associated with field-evolved resistance of pink bollworm to Bt cotton in India. PLoS ONE 2014, 9, e97900. [Google Scholar] [CrossRef] [PubMed]
  6. Ojha, A.; Sree, K.S.; Sachdev, B.; Rashmi, M.; Ravi, K.; Suresh, P.; Mohan, K.S.; Bhatnagar, R.K. Analysis of resistance to Cry1Ac in field-collected pink bollworm, Pectinophora gossypiella (Lepidoptera: Gelechiidae), populations. GM Crops Food 2014, 5, 280–286. [Google Scholar] [CrossRef] [Green Version]
  7. Fritz, M.L.; Nunziata, S.O.; Guo, R.; Tabashnik, B.E.; Carrière, Y. Mutations in a novel cadherin gene associated with Bt resistance in Helicoverpa zea. G3 Genes Genomes Genet. 2020, 10, 1563–1574. [Google Scholar] [CrossRef] [Green Version]
  8. Zhang, M.; Wei, J.; Ni, X.; Zhang, J.; Jurat-Fuentes, J.L.; Fabrick, J.A.; Carriere, Y.; Tabashnik, B.E.; Li, X. Decreased Cry1Ac activation by midgut proteases associated with Cry1Ac resistance in Helicoverpa zea. Pest Manag. Sci. 2019, 75, 1099–1106. [Google Scholar] [CrossRef]
  9. Jurat-Fuentes, J.L.; Heckel, D.; Ferré, J. Mechanisms of resistance to insecticidal proteins from Bacillus thuringiensis. Ann. Rev. Entomol. 2021, 66, 121–140. [Google Scholar] [CrossRef]
  10. Tanaka, S.; Miyamoto, K.; Noda, H.; Jurat-Fuentes, J.L.; Yoshizawa, Y.; Endo, H.; Sato, R. The ATP-binding cassette transporter subfamily C member 2 in Bombyx mori larvae is a functional receptor for Cry toxins from Bacillus thuringiensis. FEBS J. 2013, 280, 1782–1794. [Google Scholar] [CrossRef]
  11. Gahan, L.J.; Pauchet, Y.; Vogel, H.; Heckel, D.G. An ABC transporter mutation is correlated with insect resistance to Bacillus thuringiensis Cry1Ac toxin. PLoS Genet. 2010, 6, e1001248. [Google Scholar] [CrossRef] [Green Version]
  12. Xiao, Y.; Zhang, T.; Liu, C.; Heckel, D.G.; Li, X.; Tabashnik, B.E.; Wu, K. Mis-splicing of the ABCC2 gene linked with Bt toxin resistance in Helicoverpa armigera. Sci. Rep. 2014, 4, 6184. [Google Scholar] [CrossRef]
  13. Guo, Z.; Kang, S.; Chen, D.; Wu, Q.; Wang, S.; Xie, W.; Zhu, X.; Baxter, S.W.; Zhou, X.; Jurat-Fuentes, J.L.; et al. MAPK signaling pathway alters expression of midgut ALP and ABCC genes and causes resistance to Bacillus thuringiensis Cry1Ac toxin in diamondback moth. PLoS Genet. 2015, 11, e1005124. [Google Scholar] [CrossRef] [Green Version]
  14. Guo, Z.; Sun, D.; Kang, S.; Zhou, J.; Gong, L.; Qin, J.; Guo, L.; Zhu, L.; Bai, Y.; Luo, L.; et al. CRISPR/Cas9-mediated knockout of both the PxABCC2 and PxABCC3 genes confers high-level resistance to Bacillus thuringiensis Cry1Ac toxin in the diamondback moth, Plutella xylostella (L.). Insect. Biochem. Mol. Biol. 2019, 107, 31–38. [Google Scholar] [CrossRef]
  15. Banerjee, R.; Hasler, J.; Meagher, R.; Nagoshi, R.; Hietala, L.; Huang, F.; Narva, K.; Jurat-Fuentes, J.L. Mechanism and DNA-based detection of field-evolved resistance to transgenic Bt corn in fall armyworm (Spodoptera frugiperda). Sci. Rep. 2017, 7, 10877. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Pinos, D.; Martínez-Solís, M.; Herrero, S.; Ferré, J.; Hernández-Martínez, P. The Spodoptera exigua ABCC2 acts as a Cry1A receptor independently of its nucleotide binding domain II. Toxins 2019, 11, 172. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Ren, X.L.; Jiang, W.L.; Ma, Y.J.; Hu, H.Y.; Ma, X.Y.; Ma, Y.; Li, G.Q. The Spodoptera exigua (Lepidoptera: Noctuidae) ABCC2 mediates Cry1Ac cytotoxicity and, in conjunction with cadherin, contributes to cnhance Cry1Ca toxicity in Sf9 cells. J. Econ. Entomol. 2016, 109, 2281–2289. [Google Scholar] [CrossRef] [PubMed]
  18. Ocelotl, J.; Sánchez, J.; Gómez, I.; Tabashnik, B.E.; Bravo, A.; Soberón, M. ABCC2 is associated with Bacillus thuringiensis Cry1Ac toxin oligomerization and membrane insertion in diamondback moth. Sci. Rep. 2017, 7, 2386. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Pearce, S.L.; Clarke, D.F.; East, P.D.; Elfekih, S.; Gordon, K.H.J.; Jermiin, L.S.; McGaughran, A.; Oakeshott, J.G.; Papanikolaou, A.; Perera, O.P.; et al. Genomic innovations, transcriptional plasticity and gene loss underlying the evolution and divergence of two highly polyphagous and invasive Helicoverpa pest species. BMC Biol. 2017, 15, 63. [Google Scholar] [CrossRef] [Green Version]
  20. Letunic, I.; Doerks, T.; Bork, P. SMART 7: Recent updates to the protein domain annotation resource. Nucleic Acids Res. 2012, 40, D302–D305. [Google Scholar] [CrossRef] [PubMed]
  21. Perera, O.P.; Little, N.S.; Pierce, C.A., III. CRISPR/Cas9 mediated high efficiency knockout of the eye color gene Vermillion in Helicoverpa zea (Boddie). PLoS ONE 2018, 13, e0197567. [Google Scholar] [CrossRef] [Green Version]
  22. Hsu, P.D.; Scott, D.A.; Weinstein, J.A.; Ran, F.A.; Konermann, S.; Agarwala, V.; Li, Y.; Fine, E.J.; Wu, X.; Shalem, O.; et al. DNA targeting specificity of RNA-guided Cas9 nucleases. Nat. Biotechnol. 2013, 31, 827–832. [Google Scholar] [CrossRef] [PubMed]
  23. Gore, J.; Adamczyk, J.J.; Blanco, C.A. Selective feeding of tobacco budworm and bollworm (Lepidoptera: Noctuidae) on meridic diet with different concentrations of Bacillus thuringiensis proteins. J. Econ. Entomol. 2005, 98, 88–94. [Google Scholar] [CrossRef] [PubMed]
  24. Jinek, M.; Chylinski, K.; Fonfara, I.; Hauer, M.; Doudna, J.A.; Charpentier, E. A programmable dual-RNA-guided DNA endonuclease in adaptive bacterial immunity. Science 2012, 337, 816–821. [Google Scholar] [CrossRef]
  25. Mali, P.; Aach, J.; Stranges, P.B.; Esvelt, K.M.; Moosburner, M.; Kosuri, S.; Yang, L.; Church, G.M. CAS9 transcriptional activators for target specificity screening and paired nickases for cooperative genome engineering. Nat. Biotechnol. 2013, 31, 833–838. [Google Scholar] [CrossRef]
  26. Little, N.S.; Elkins, B.H.; Mullen, R.M.; Perera, O.P.; Parys, K.A.; Allen, K.C.; Boykin, D.L. Differences between two populations of bollworm, Helicoverpa zea (Lepidoptera: Noctuidae), with variable measurements of laboratory susceptibilities to Bt toxins exposed to non-Bt and Bt cottons in large field cages. PLoS ONE 2019, 14, e0212567. [Google Scholar] [CrossRef] [PubMed]
  27. Finney, D. The adjustment for a natural response rate in probit analysis. Ann. Appl. Biol. 1949, 36, 187–195. [Google Scholar] [CrossRef] [PubMed]
  28. SAS-Institute. SAS Software, version 9.4; SAS Institute: Cary, NC, USA, 2013. [Google Scholar]
  29. Wolfersberger, M.; Luthy, P.; Maurer, A.; Parenti, P.; Sacchi, V.F.; Giordana, B.; Hanozet, G.M. Preparation and partial characterization of amino acid transporting brush border membrane vesicles from the larval midgut of the cabbage butterfly (Pieris brassicae). Comp. Biochem. Phys. B 1987, 86, 301–308. [Google Scholar] [CrossRef]
  30. Jurat-Fuentes, J.L.; Gould, F.L.; Adang, M.J. Altered Glycosylation of 63- and 68-kilodalton microvillar proteins in Heliothis virescens correlates with reduced Cry1 toxin binding, decreased pore formation, and increased resistance to Bacillus thuringiensis Cry1 toxins. Appl. Environ. Microb. 2002, 68, 5711–5717. [Google Scholar] [CrossRef] [Green Version]
  31. Jurat-Fuentes, J.L.; Karumbaiah, L.; Jakka, S.R.; Ning, C.; Liu, C.; Wu, K.; Jackson, J.; Gould, F.; Blanco, C.; Portilla, M.; et al. Reduced levels of membrane-bound alkaline phosphatase are common to lepidopteran strains resistant to Cry toxins from Bacillus thuringiensis. PLoS ONE 2011, 6, e17606. [Google Scholar] [CrossRef] [Green Version]
  32. Perera, O.P.; Willis, J.D.; Adang, M.J.; Jurat-Fuentes, J.L. Cloning and characterization of the Cry1Ac-binding alkaline phosphatase (HvALP) from Heliothis virescens. ICES J. Mar. Sci. J. Du Cons. 2009, 39, 294–302. [Google Scholar] [CrossRef]
  33. Gouffon, C.; Van Vliet, A.; Van Rie, J.; Jansens, S.; Jurat-Fuentes, J.L. Binding sites for Bacillus thuringiensis Cry2Ae toxin on heliothine brush border membrane vesicles are not shared with Cry1A, Cry1F, or Vip3A toxin. Appl. Environ. Microb. 2011, 77, 3182–3188. [Google Scholar] [CrossRef] [Green Version]
  34. Abbott, W.S. A method of computing the effectiveness of an insecticide. J. Econ. Entomol. 1925, 18, 265–267. [Google Scholar] [CrossRef]
  35. Heckel, D.G. The essential and enigmatic role of ABC transporters in Bt resistance of noctuids and other insect pests of agriculture. Insects 2021, 12, 389. [Google Scholar] [CrossRef]
  36. Zhou, Z.; Wang, Z.; Liu, Y.; Liang, G.; Shu, C.; Song, F.; Zhou, X.; Bravo, A.; Soberón, M.; Zhang, J. Identification of ABCC2 as a binding protein of Cry1Ac on brush border membrane vesicles from Helicoverpa armigera by an improved pull-down assay. MicrobiologyOpen 2016, 5, 659–669. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Wang, J.; Ma, H.; Zhao, S.; Huang, J.; Yang, Y.; Tabashnik, B.E.; Wu, Y. Functional redundancy of two ABC transporter proteins in mediating toxicity of Bacillus thuringiensis to cotton bollworm. PLoS Pathog. 2020, 16, e1008427. [Google Scholar] [CrossRef]
  38. Liu, Z.; Fu, S.; Ma, X.; Baxter, S.W.; Vasseur, L.; Xiong, L.; Huang, Y.; Yang, G.; You, S.; You, M. Resistance to Bacillus thuringiensis Cry1Ac toxin requires mutations in two Plutella xylostella ATP-binding cassette transporter paralogs. PLoS Pathog. 2020, 16, e1008697. [Google Scholar] [CrossRef] [PubMed]
  39. Zhang, D.; Jin, M.; Yang, Y.; Zhang, J.; Yang, Y.; Liu, K.; Soberón, M.; Bravo, A.; Xiao, Y.; Wu, K. Synergistic resistance of Helicoverpa armigera to Bt toxins linked to cadherin and ABC transporters mutations. Insect. Biochem. Mol. Biol. 2021, 137, 103635. [Google Scholar] [CrossRef]
  40. Liu, L.; Chen, Z.; Yang, Y.; Xiao, Y.; Liu, C.; Ma, Y.; Soberon, M.; Bravo, A.; Yang, Y.; Liu, K. A single amino acid polymorphism in ABCC2 loop 1 is responsible for differential toxicity of Bacillus thuringiensis Cry1Ac toxin in different Spodoptera (Noctuidae) species. Insect. Biochem. Mol. Biol. 2018, 100, 59–65. [Google Scholar] [CrossRef]
  41. Liu, Y.; Jin, M.; Wang, L.; Wang, H.; Xia, Z.; Yang, Y.; Bravo, A.; Soberón, M.; Xiao, Y.; Liu, K. SfABCC2 transporter extracellular loops 2 and 4 are responsible for the Cry1Fa insecticidal specificity against Spodoptera frugiperda. Insect. Biochem. Mol. Biol. 2021, 135, 103608. [Google Scholar] [CrossRef]
  42. Tanaka, S.; Endo, H.; Adegawa, S.; Iizuka, A.; Imamura, K.; Kikuta, S.; Sato, R. Bombyx mori ABC transporter C2 structures responsible for the receptor function of Bacillus thuringiensis Cry1Aa toxin. Insect. Biochem. Mol. Biol. 2017, 91, 44–54. [Google Scholar] [CrossRef]
  43. Flagel, L.; Lee, Y.; Wanjugi, H.; Swarup, S.; Brown, A.; Wang, J.; Kraft, E.; Greenplate, J.; Simmons, J.; Adams, N. Mutational disruption of the ABCC2 gene in fall armyworm, Spodoptera frugiperda, confers resistance to the Cry1Fa and Cry1A. 105 insecticidal proteins. Sci. Rep. 2018, 8, 7255. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Wei, J.; Zhang, M.; Liang, G.; Li, X. Alkaline phosphatase 2 is a functional receptor of Cry1Ac but not Cry2Ab in Helicoverpa zea. Insect. Mol. Biol. 2019, 28, 372–379. [Google Scholar] [CrossRef] [PubMed]
  45. Ma, Y.; Zhang, J.; Xiao, Y.; Yang, Y.; Liu, C.; Peng, R.; Yang, Y.; Bravo, A.; Soberón, M.; Liu, K. The cadherin Cry1Ac binding-region is necessary for the cooperative effect with ABCC2 transporter enhancing insecticidal activity of Bacillus thuringiensis Cry1Ac toxin. Toxins 2019, 11, 538. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Genes 12 01522 g001 550
Figure 1. Organization of Helicoverpa zea ATP binding cassette transporter subfamily C2 (ABCC2) transporter gene and CRISPR/Cas9 mediated deletions generated and tested in this study. Exons and introns, spaced to scale, are shown by blue hatched rectangles and solid black lines, respectively. Numbers below the predicted polypeptide represent amino acids. (A). Wild type ABCC2 gene containing 25 exons and predicted polypeptide with transmembrane (TMD) and nucleotide-binding (NBD) domains. Exons 1, 5, 10, 15, 20, and 25 are numbered above the gene diagram. Approximate location of the promoter and exons coding for TMD 1, TMD 2, NBD 1, and NBD 2 are shown by horizontal double arrows. Locations for target sites in exons 3, 8, 13, 19, 21, 22, and 24 for which single guide RNA (sgRNA) were designed are shown with red down arrows. (B). Knockout line SPM-A28C generated by sgRNA targeting exons 3 and 13 which truncated the ABCC2 protein at amino acid 127. (C). Knockout line SPM-B1 with a deletion from exon 13 to 19 recovered from the experiment targeting exons 13, 19, and 21 in a polypeptide that had a deletion of 255 amino acids at the beginning of the TMD 2 (starting at amino acid 760). Alternative splicing using a splice donor (GT) site in exon 13 to exon 19 restored the NBD2 domain of the ABCC2 protein. (D,E). Knockouts generated by sgRNA targeting exons 21, 22, and 24. In SPM-8 (D), a 7 bp deletion in exon 21 (from 93,539 to 93,545 bp in genomic scaffold) truncated the polypeptide at amino acid 1122 and added eight random amino acids, inactivating NBD2. In knockout line SPM-16 (E), a deletion of 426 bp from 93,540 to 93,965 (exon 21 to intron 21) created an open reading frame and extended the truncated NBD2 polypeptide sequence with 34 random amino acids. In this knockout line there was an additional 4 bp insertion at the exon 22 sgRNA target site (position 94,196 bp of the genomic scaffold) that did not affect the outcome due to truncation of the protein by the deletion in exon 21. Red arrows in (B–E) indicate the position of sgRNA that caused the mutation in each line.
Figure 1. Organization of Helicoverpa zea ATP binding cassette transporter subfamily C2 (ABCC2) transporter gene and CRISPR/Cas9 mediated deletions generated and tested in this study. Exons and introns, spaced to scale, are shown by blue hatched rectangles and solid black lines, respectively. Numbers below the predicted polypeptide represent amino acids. (A). Wild type ABCC2 gene containing 25 exons and predicted polypeptide with transmembrane (TMD) and nucleotide-binding (NBD) domains. Exons 1, 5, 10, 15, 20, and 25 are numbered above the gene diagram. Approximate location of the promoter and exons coding for TMD 1, TMD 2, NBD 1, and NBD 2 are shown by horizontal double arrows. Locations for target sites in exons 3, 8, 13, 19, 21, 22, and 24 for which single guide RNA (sgRNA) were designed are shown with red down arrows. (B). Knockout line SPM-A28C generated by sgRNA targeting exons 3 and 13 which truncated the ABCC2 protein at amino acid 127. (C). Knockout line SPM-B1 with a deletion from exon 13 to 19 recovered from the experiment targeting exons 13, 19, and 21 in a polypeptide that had a deletion of 255 amino acids at the beginning of the TMD 2 (starting at amino acid 760). Alternative splicing using a splice donor (GT) site in exon 13 to exon 19 restored the NBD2 domain of the ABCC2 protein. (D,E). Knockouts generated by sgRNA targeting exons 21, 22, and 24. In SPM-8 (D), a 7 bp deletion in exon 21 (from 93,539 to 93,545 bp in genomic scaffold) truncated the polypeptide at amino acid 1122 and added eight random amino acids, inactivating NBD2. In knockout line SPM-16 (E), a deletion of 426 bp from 93,540 to 93,965 (exon 21 to intron 21) created an open reading frame and extended the truncated NBD2 polypeptide sequence with 34 random amino acids. In this knockout line there was an additional 4 bp insertion at the exon 22 sgRNA target site (position 94,196 bp of the genomic scaffold) that did not affect the outcome due to truncation of the protein by the deletion in exon 21. Red arrows in (B–E) indicate the position of sgRNA that caused the mutation in each line.
Genes 12 01522 g001
Genes 12 01522 g002 550
Figure 2. Agarose gel electrophoresis of amplicons of ATP binding cassette transporter subfamily C2 (ABCC2) knockout lines and wild type insects of H. zea. Lanes 2, 4, 6, and 8: SPM-8, SPM-16, SPM-B1, and SPM-A28C, respectively. Lanes 3, 5, 7, and 9: amplicons from wild type insects. Lane 1 and 10: DNA ladder with molecular weights of the main bands (in Kbp) marked.
Figure 2. Agarose gel electrophoresis of amplicons of ATP binding cassette transporter subfamily C2 (ABCC2) knockout lines and wild type insects of H. zea. Lanes 2, 4, 6, and 8: SPM-8, SPM-16, SPM-B1, and SPM-A28C, respectively. Lanes 3, 5, 7, and 9: amplicons from wild type insects. Lane 1 and 10: DNA ladder with molecular weights of the main bands (in Kbp) marked.
Genes 12 01522 g002
Genes 12 01522 g003 550
Figure 3. Saturation Cry1Ac-binding assays with midgut brush border membrane vesicles (BBMV) from larvae of wild type (SIMRU) and gene-edited (SPM-8 and SPM-B1) strains of Helicoverpa zea. Data points shown are the specific binding estimated from total and non-specific binding (as detailed in Materials and Methods) from two independent experiments each performed in duplicate. Curves are the result from non-linear regression fitting of the data to a one-site binding model as the best fit.
Figure 3. Saturation Cry1Ac-binding assays with midgut brush border membrane vesicles (BBMV) from larvae of wild type (SIMRU) and gene-edited (SPM-8 and SPM-B1) strains of Helicoverpa zea. Data points shown are the specific binding estimated from total and non-specific binding (as detailed in Materials and Methods) from two independent experiments each performed in duplicate. Curves are the result from non-linear regression fitting of the data to a one-site binding model as the best fit.
Genes 12 01522 g003
Table
Table 1. CRISPR RNA (crRNA) sequence (red) and protospacer adjacent motif (blue) used for targeting different exons in the Helicoverpa zea ATP binding cassette transporter subfamily C2 (ABCC2) gene (GenBank KY701524). Target sites in the coding DNA strand are denoted by a plus (+) sign, and those on the reverse strand are denoted by a minus (−) sign.
Table 1. CRISPR RNA (crRNA) sequence (red) and protospacer adjacent motif (blue) used for targeting different exons in the Helicoverpa zea ATP binding cassette transporter subfamily C2 (ABCC2) gene (GenBank KY701524). Target sites in the coding DNA strand are denoted by a plus (+) sign, and those on the reverse strand are denoted by a minus (−) sign.
crRNA NamecrRNA SequenceStrandTarget
Hz_Ex3GCTACTGTCGTACTGGTCGGTGG+Exon 3
Hz_Ex8TTAACAAAGTAAGCGCATCGTGG+Exon 8
Hz_Ex13CTGCCGACTATTGGTTGAGTTGG−Exon 13
Hz_Ex19TTGCTCCATATTGGTCTCCGTGG−Exon 19
Hz_Ex21CGGCAAGTCATCGCTCATCGCGG+Exon 21
Hz_Ex22ACAGCGACGACGATATTTGGAGG+Exon 22
Hz_Ex24GGTCATGGACCAGGGCGAAGTGG+Exon 24
Table
Table 2. Potential off-target sequences for single guide RNA (sgRNA) targeting exon 13 and 19 in the genome of Helicoverpa zea. Scaffold number, nucleotide position, off-target sequence, and the strand of DNA are given. Nucleotides of the protospacer adjacent motif are shown in red bold text and 12-nucleotide seed sequences are shown in bold Italics. Nucleotides identical to each sgRNA are shown by a period.
Table 2. Potential off-target sequences for single guide RNA (sgRNA) targeting exon 13 and 19 in the genome of Helicoverpa zea. Scaffold number, nucleotide position, off-target sequence, and the strand of DNA are given. Nucleotides of the protospacer adjacent motif are shown in red bold text and 12-nucleotide seed sequences are shown in bold Italics. Nucleotides identical to each sgRNA are shown by a period.
sgRNA/ScaffoldPositionSequenceStrand
sgRNA Exon 13 ACTCAACCAATAGTCGGCAGTGG
KZ118710.1289298-289320CAAT...T............TGG+
sgRNA Ex19TTGCTCCATATTGGTCTCCGTGG
KZ117297.1190397-190417G..TGATTC...........TGG−
204522-204544GCTAGTA..T...T......TGG-
KZ117617.170340-70362C.TACT.T.A..........TGG-
Table
Table 3. Bioassay data for control (SIMRU) and ATP binding cassette transporter subfamily C2 (ABCC2) knockout lines of Helicoverpa zea to estimate the concentrations Cry1Ac that killed 50 and 90% larvae (LC50 and LC90, respectively) and their fiducial limits (Lower-Upper). Slope of the regression line, standard error of slope, Chi-square value of the slope and the p-value are also presented. Significantly different LC50 values (based on non-overlapping fiducial limits) are denoted by different symbols (†, ‡, and §). Cry1Ac concentrations (µg/mL) used in bioassays: a 0, 1, 10, 50,100; b 0, 1, 10, 100, 500; c 0, 1, 10, 100, 1000; d 0, 1, 10, 100; e 0, 0.1, 1, 10, 100. RR = LC50 ratio between edited and control strain.
Table 3. Bioassay data for control (SIMRU) and ATP binding cassette transporter subfamily C2 (ABCC2) knockout lines of Helicoverpa zea to estimate the concentrations Cry1Ac that killed 50 and 90% larvae (LC50 and LC90, respectively) and their fiducial limits (Lower-Upper). Slope of the regression line, standard error of slope, Chi-square value of the slope and the p-value are also presented. Significantly different LC50 values (based on non-overlapping fiducial limits) are denoted by different symbols (†, ‡, and §). Cry1Ac concentrations (µg/mL) used in bioassays: a 0, 1, 10, 50,100; b 0, 1, 10, 100, 500; c 0, 1, 10, 100, 1000; d 0, 1, 10, 100; e 0, 0.1, 1, 10, 100. RR = LC50 ratio between edited and control strain.
Insect LineABCC2 Domain KnockoutLC50 µg/mL
(Lower-Upper)		RRLC90 µg/mL
(Lower-Upper)		Slope ± SEΧ2 Slope
SPM-A28C aComplete21.87 †
(13.92–33.77)		9.2251.65
(141.04–576.35)		1.22 ± 0.1474.00
p < 0.0001
SPM-B1 bTM2 to NBD294.32 ‡
(65.75–137.50)		39.8455.69
(282.81–946.91)		1.87 ± 0.2554.58
p < 0.0001
SPM-16 cNBD2 (insertion)42.58 †
(24.71–74.62)		18.01118.00
(491.32–3876.00)		0.90 ± 0.1165.79
p < 0.0001
SPM-8 dNBD217.35 †,‡
(10.82–29.13)		7.3355.42
(156.54–1357.00)		0.98 ± 0.1352.79
p < 0.0001
Control eNone2.37 §
(1.75–3.20)		122.31
(14.85–37.31)		1.32 ± 0.10162.63
p < 0.0001
Table
Table 4. Binding parameters from saturation Cry1Ac binding assays and goodness of fit criteria for wild type (SIMRU) and gene edited (SPM-8 and SPM-B1) strains of Helicoverpa zea. Kd = dissociation constant (in nM units); Bmax = concentration of receptors (in ng/mg BBMV protein units); SE = standard error; P = statistic for the parameter estimate (significant estimates (p < 0.05) are underlined); R2 = coefficient of determination for the regression model used for fitting (one binding site model).
Table 4. Binding parameters from saturation Cry1Ac binding assays and goodness of fit criteria for wild type (SIMRU) and gene edited (SPM-8 and SPM-B1) strains of Helicoverpa zea. Kd = dissociation constant (in nM units); Bmax = concentration of receptors (in ng/mg BBMV protein units); SE = standard error; P = statistic for the parameter estimate (significant estimates (p < 0.05) are underlined); R2 = coefficient of determination for the regression model used for fitting (one binding site model).
StrainKd ± SEPBmax ± SEPR2
SIMRU1.55 ± 0.370.000249.61 ± 8.080.00240.9978
SPM-80.16 ± 0.090.10945.20 ± 0.840.00020.9883
SPM-B10.04 ± 0.050.42443.07 ± 0.500.00020.7669
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Perera, O.P.; Little, N.S.; Abdelgaffar, H.; Jurat-Fuentes, J.L.; Reddy, G.V.P. Genetic Knockouts Indicate That the ABCC2 Protein in the Bollworm Helicoverpa zea Is Not a Major Receptor for the Cry1Ac Insecticidal Protein. Genes 2021, 12, 1522. https://doi.org/10.3390/genes12101522

AMA Style

Perera OP, Little NS, Abdelgaffar H, Jurat-Fuentes JL, Reddy GVP. Genetic Knockouts Indicate That the ABCC2 Protein in the Bollworm Helicoverpa zea Is Not a Major Receptor for the Cry1Ac Insecticidal Protein. Genes. 2021; 12(10):1522. https://doi.org/10.3390/genes12101522

Chicago/Turabian Style

Perera, Omaththage P., Nathan S. Little, Heba Abdelgaffar, Juan Luis Jurat-Fuentes, and Gadi V. P. Reddy. 2021. "Genetic Knockouts Indicate That the ABCC2 Protein in the Bollworm Helicoverpa zea Is Not a Major Receptor for the Cry1Ac Insecticidal Protein" Genes 12, no. 10: 1522. https://doi.org/10.3390/genes12101522

APA Style

Perera, O. P., Little, N. S., Abdelgaffar, H., Jurat-Fuentes, J. L., & Reddy, G. V. P. (2021). Genetic Knockouts Indicate That the ABCC2 Protein in the Bollworm Helicoverpa zea Is Not a Major Receptor for the Cry1Ac Insecticidal Protein. Genes, 12(10), 1522. https://doi.org/10.3390/genes12101522

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop