A versatile contribution of both aminopeptidases N and ABC transporters to Bt Cry1Ac toxicity in the diamondback moth

Background Biopesticides and transgenic crops based on Bacillus thuringiensis (Bt) toxins are extensively used to control insect pests, but the rapid evolution of insect resistance seriously threatens their effectiveness. Bt resistance is often polygenic and complex. Mutations that confer resistance occur in midgut proteins that act as cell surface receptors for the toxin, and it is thought they facilitate its assembly as a membrane-damaging pore. However, the mechanistic details of the action of Bt toxins remain controversial. Results We have examined the contribution of two paralogous ABC transporters and two aminopeptidases N to Bt Cry1Ac toxicity in the diamondback moth, Plutella xylostella, using CRISPR/Cas9 to generate a series of homozygous polygenic knockout strains. A double-gene knockout strain, in which the two paralogous ABC transporters ABCC2 and ABCC3 were deleted, exhibited 4482-fold resistance to Cry1A toxin, significantly greater than that previously reported for single-gene knockouts and confirming the mutual functional redundancy of these ABC transporters in acting as toxin receptors in P. xylostella. A double-gene knockout strain in which APN1 and APN3a were deleted exhibited 1425-fold resistance to Cry1Ac toxin, providing the most direct evidence to date for these APN proteins acting as Cry1Ac toxin receptors, while also indicating their functional redundancy. Genetic crosses of the two double-gene knockouts yielded a hybrid strain in which all four receptor genes were deleted and this resulted in a > 34,000-fold resistance, indicating that while both types of receptor need to be present for the toxin to be fully effective, there is a level of functional redundancy between them. The highly resistant quadruple knockout strain was less fit than wild-type moths, but no fitness cost was detected in the double knockout strains. Conclusion Our results provide direct evidence that APN1 and APN3a are important for Cry1Ac toxicity. They support our overarching hypothesis of a versatile mode of action of Bt toxins, which can compensate for the absence of individual receptors, and are consistent with an interplay among diverse midgut receptors in the toxins’ mechanism of action in a super pest. Supplementary Information The online version contains supplementary material available at 10.1186/s12915-022-01226-1.


Background
The entomopathogenic bacterium Bacillus thuringiensis (Bt) produces different types of insecticidal crystal proteins which are used to control insect pests both in sprayable formulations and transgenic crops [1]. Btbased biotechnology products are the most successful alternatives to chemical pesticides for insect pest control [2,3]. However, the economic, environmental, and social benefits of these Bt products are threatened by the evolution of insect resistance [4,5]. To date, field-evolved Bt resistance has been observed in at least nine different lepidopteran and one coleopteran species [6]. Although a range of countermeasures have been introduced to prevent the development of resistance in diverse insect populations [7][8][9][10][11], an in-depth unraveling of the molecular mechanisms underlying Bt mechanism of action and resistance is particularly important for the improvement of Bt technology and its sustainability [12].
Bt Cry toxins are toxic upon ingestion by the host insect and target the midgut epithelium, the interactions of Cry toxins with their midgut-specific receptors trigger toxin oligomerization and pore formation [13,14]. The best established Bt resistance mechanism in insects is via alterations of these midgut receptors such as cadherin (CAD), aminopeptidases-N (APNs), alkaline phosphatases (ALPs), and ABC transporters, that interrupt toxin binding [15,16]. Among these receptors, structural mutations or expression alterations of ABC transporters (especially ABCC2) have been widely reported to be associated with high levels of resistance in lepidopteran insects [17][18][19][20][21][22][23][24][25][26][27][28]. In addition, as the first identified functional receptors of Bt Cry toxins, GPI-anchored APN proteins also play a crucial role [27,[29][30][31][32]. The alteration of a single receptor gene can however affect only one pathway of toxicity for Cry toxins [33]. Numerous instances of pest resistance to chemical insecticides are polygenic, and the eventual levels of resistance are decided by the interactions of different resistance alleles [34,35]. Previous studies suggested that the interactions of ABC transporters with cadherin could enhance cell toxicity, and co-operation between these two receptors was demonstrated by in vitro ectopic expression in different cell lines [36][37][38][39].
With the emergence of complete whole genome information for insects, the clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) system has been used to explore gene function. The creation of multiple or large fragment gene deletions by the CRISPR method makes it possible to decipher the relationship between functional genes in host insects [40][41][42]. Besides the easily achieved single-gene knockouts [25,27], double-gene knockout of paralogous midgut receptor genes with the CRISPR/Cas9 system has recently been conducted to confirm their in vivo functions with respect to Bt toxin receptors [43][44][45]. The molecular mechanism underlying insect resistance to Bt toxins is often perplexing and multifaceted, with the resistance trait associated with different types of midgut receptors [46]. Resistance to Cry1A toxins in the strains of Plutella xylostella that we are studying is associated with the downregulation of multiple receptors including ABCC2, ABCC3, ABCB1, ABCG1, APN1, APN3a and ALP [25,27] but is not associated with any mutation in, or change in expression of, cadherin [47,48].
In this study, CRISPR/Cas9-mediated double knockouts of two pairs of the same types of midgut receptor genes (PxAPN1/PxAPN3a and PxABCC2/PxABCC3) resulted in small and moderate levels of resistance enhancement, respectively, when compared to knocking out the individual receptors [25,27]. However, when all four receptors were knocked out, an extremely high level of resistance was observed. This result provides a multiple pathway model for understanding how different types of midgut receptors are involved and interplay in Bt pathogenicity and resistance.

Construction of double knockout strains by CRISPR/Cas9
Based on the fact that the PxABCC2 and PxABCC3 receptor genes are adjacent on the P. xylostella genome, as are the PxAPN1 and PxAPN3a receptor genes, a dual sgRNA CRISPR/Cas9 approach was used to knockout each pair of genes. Fresh eggs from the susceptible DBM1Ac-S strain were collected and co-injected with a mixture of Cas9 protein and two sgRNAs targeting either PxABCC3 and PxABCC2 (group 1) or PxAPN3a and PxAPN1 genes (group 2) to obtain G0 progeny (Figs. 1a, b and 2a).
Subsequently, single-pair reciprocal crosses between G0 and DBM1Ac-S were performed to generate G1 progeny (Table 1). Single pairs which produced fertile progeny were screened and ten second-instar P. xylostella larvae from each of these pairs were pooled to extract gDNA samples. PCR amplification with primer pair 3-2-F/2-3-R exhibited a 516-bp fragment (indicating a successful deletion) in only a single pair in group 1 (Fig. 1c). This PCR product was then directly sequenced using 3-2-F as primer. The sequencing chromatogram revealed successful mutagenesis of PxABCC2 and PxABCC3 genes (Fig. 1d) and the remaining larvae were reared to pupation. In group 2, PCR amplification with primer pair 3a-1-F/1-13-R revealed a 716 bp fragment in two single pairs (Fig. 2b), PCR products were directly sequenced in these two pairs using 3a-1-F to confirm the mutagenesis within the PxAPN3a and PxAPN1 locus (Fig. 2c). The pair group with the most larvae was chosen for further rearing. A nondestructive method for extracting the gDNA samples was used to screen pupae for the presence of the mutation (Figs. 1d and 2c) and heterozygotes containing the mutation were sib-mated to produce G2 progeny. Finally, homozygous individuals identified by PCR and DNA sequencing were further sib-crossed to establish two stable homozygous mutant strains designated C2-3KO (ABC transporter) and N1-3aKO (APN) (Additional file 1: Fig. S1).

Combinational inactivation of ABC transporters and GPIanchored proteins
To produce a combined mutant strain containing both of the above deletions in the ABC transporter and APN genes, 30 single-pair reciprocal crosses were performed between C2-3KO and N1-3aKO (Additional file 1: Fig.   Fig. 1 CRISPR/Cas9-mediated double gene mutation of PxABCC2 and PxABCC3 in P. xylostella. a The contiguous genomic structure of PxABCC2 and PxABCC3. Exons are shown as boxes except for the first and last exons which are represented by triangles and arrows respectively. The sgRNA target sequences of PxABCC2 and PxABCC3 are highlighted by green text, the PAM sites are highlighted in red and underlined. The remaining sequences after cleavage are shown in blue. b Protein structures of PxABCC2 (right) and PxABCC3 (left), as well as the cleavage position following editing. c Genotyping of individual P. xylostella for deletion of PxABCC2 and PxABCC3 according to the banding profile of PCR products from the four primer pairs. d Representative chromatograms of direct sequencing of PCR products with the primer pair 3-2-F/2-3-R, the arrow indicates the junction of PxABCC2 and PxABCC3 after editing. In the nomenclature of these primers, the first number denotes the gene, the second number indicates the exon/intron. M: DNA marker III   [25,27].
Bioassays were subsequently performed to test the susceptibility changes to Cry1A protoxins among the newly-built strains in this study with the susceptible DBM1Ac-S used as a control ( Table 2). The results indicated that the LC 50 values of the double knockout strains C2-3KO and N1-3aKO to Cry1Ac protoxin were approximately 4482 and 1425-fold higher than that of the susceptible DBM1Ac-S strain, respectively. The quadruple knockout strain C-NKO showed > 34,000-fold resistance to Bt Cry1Ac protoxin. The three polygenic knockout strains also manifested high levels of cross-resistance to Cry1Ab protoxin (C2-3KO: 591-fold; N1-3aKO: 132-fold; C-NKO: > 4200-fold) but no cross-resistance to Cry1Aa protoxin was detected (C2-3KO: 1.2-fold; N1-3aKO: 1.3-fold; C-NKO: 1.5-fold). These data demonstrated the involvement of PxABCC2 and PxABCC3, as well as PxAPN1 and PxAPN3a genes in the mechanism of action of Cry1Ac and Cry1Ab in P. xylostella.
Inheritance of resistance to Cry1Ac in the engineered strains To analyze the inheritance of Cry1Ac resistance in the polygenic knockout strains, F1 progeny were mass crossed between pairs of strains from DBM1Ac-S, C2-3KO, N1-3aKO, and C-NKO, and resulting larvae were screened with 10 mg/L Cry1Ac protoxin to remove all susceptible individuals. As expected, the mortality was 100% for the susceptible DBM1Ac-S strain while the F1 progeny of the knockout strains had a high survival rate (96-100%) (Additional file 1: Table S2). F1 progeny from crosses between either double mutant strain C2-3KO or N1-3aKO or the quadruple knockout C-NKO strain and DBM1Ac-S showed relatively low survival rates (8-12%) as did progeny from the C2-3KO vs N1-3aKO cross (4%). In contrast, progeny from crosses between C-NKO and the two double mutant strains showed high survival rates (92-96%) (Fig. 3). These results confirm that the C2-3KO and N1-3aKO strains share different, essentially recessive, resistance loci, and that allelic complementation restored susceptibility to Cry1Ac protoxin and showed a wild-type phenotype. Table 2 Resistance to Cry1 toxin in larvae from the susceptible (DBM1Ac-S) and polygenic knockout strains Both Cry1Ac and Cry1Ab toxins produced < 15% mortality in larvae from the C-NKO strain when examined with the Cry1Ac concentration of 2.5× 10 4 mg/L and 3× 10 3 mg/L

Genetic linkage of polygenic mutations with Cry1Ac resistance
We further carried out genetic linkage analysis to test the cosegregation of different multiple gene inactivations with resistance to Cry1Ac in our polygenic knockout strains. Single-pair crosses between a male from the polygenic knockout strains (C2-3KO, N1-3aKO, or C-NKO) and a female from the susceptible DBM1Ac-S were performed to generate F1 progeny. After five reciprocal crosses between the F1 progeny and the respective polygenic knockout moths were performed to produce ten F2 backcross families (named backcross a and backcross b, as labeled in Fig. 4a), individuals were or were not treated with a diagnostic dose of Cry1Ac protoxin ( Fig. 4a, b). The genotypes of individuals were assessed by PCR analysis (Fig. 4c-e). Among 300 untreated progeny from the C2-3KO cross, almost half of the individuals were heterozygous (152) and the rest (148) were homozygous for the PxABCC2 and PxABCC3 mutations (Fig. 4c). The genotype ratio (rr:rs) between the number of moths in backcross a and backcross b were 73:77 and 75:75, respectively, which was consistent with the 1:1 random assortment ratio (p > 0.5 or p = 1.0, Fisher's exact test) (Fig. 4b, top right). When treated with the diagnostic dose of Cry1Ac (20 mg/L), only homozygous individuals were detected and accounted for about 50% (51% for backcross a, 49% for backcross b), which followed the expected Mendelian inheritance of a recessive resistance trait (p < 0.001, Fisher's exact test) (Fig. 4b, top left). A similar result was obtained in N1-3aKO with the genotype ratio (rr:rs) between the number of moths in backcross a and backcross b being 80:70 and 76:74 (p > 0.5 or p = 1.0, Fisher's exact test) (Fig. 4b, middle right). Rearing the larvae of each backcross with the diagnostic dose of Cry1Ac (20 mg/L) caused approximately 50% mortality (51% for backcross a, 49% for backcross b) (p < 0.001, Fisher's exact test) (Fig. 4b, middle left). These results suggest that the mutations of PxABCC2/PxABCC3 or PxAPN1/PxAPN3a genes are closely linked with resistance to Cry1Ac in the C2-3KO and N1-3aKO strains, respectively.
For the C-NKO strain, the genotype ratio (rr:rs) between the number of moths in backcross a and backcross b were 71:79 and 73:77 (Fig. 4b, bottom right). Rearing the larvae of each backcross with the diagnostic dose of Cry1Ac (4000 mg/L) caused 50% mortality for both backcrosses a and b (p < 0.001, Fisher's exact test). This further demonstrated that the combined mutations of ABC transporters and GPI-anchored proteins are also tightly linked to Cry1Ac resistance in the C-NKO strain.

Fitness cost analysis
Fitness costs related to life-history traits are frequently associated with insect resistance to Bt pathogens thereby influencing the evolution of resistance [49]. To confirm Fig. 3 Interstrain allelic complementation test with a diagnostic dose of Cry1Ac protoxin. F1 progeny were generated by six interstrain crosses between DBM1Ac-S, C2-3KO, N1-3aKO, and C-NKO. The mortality of the resulting strains (50 larvae from each strain) and their F1 progeny (100 larvae from each F1 offspring) was examined with a diagnostic dose (10 mg/L) of Cry1Ac protoxin, which should kill 100% heterozygous F1 larvae. Significant differences in larval mortality were analyzed by one-way ANOVA. The p values from Holm-Sidak's test, *p < 0.05, **p < 0.01, ***p < 0.001, ns, not significant whether individuals of the different knockout strains had growth penalties, a series of biological parameters including pupation rate, pupal weight, pupal duration, and hatching percentage were assessed. We observed that the double mutant C2-3KO and N1-3aKO strains had no significant differences when compared to the susceptible DBM1Ac-S strain, whereas the quadruple knockout strain C-NKO displayed some fitness disadvantages in all of the measured parameters ( Fig. 5a-d). Previously, we had found no fitness costs associated with Cry1Ac resistance in a strain in which down-regulation of the functional receptors PxABCC2/3 and PxAPN1/3a was accompanied by the up-regulation of the non-receptor paralogs PxABCC1 and PxAPN5/6. We speculated that this compensation mechanism minimized fitness costs while maintaining a high level of resistance [22,27,50]. Given the fitness costs associated with the quadruple knockout C-NKO (but not the two double knockout strains) we investigated the expression of the nonreceptor paralogs in these three strains. We observed no differences in expression of PxAPN5/6 or PxABCC1 in any of the knockout strains compared to the parent (Fig.  5e). This observation was consistent with fitness costs in the quadruple knockout strain being associated with the Fig. 4 Genetic linkage analysis of Cry1Ac resistance in the polygenic knockout strains of P. xylostella. a The crossing strategy for analysis of the linkage between polygenic knockout and Cry1Ac resistance. b To examine the genetic linkage between multiple resistance alleles and Cry1Ac resistance, 10 backcross families from the single-pair cross between C2-3KO (rr) and their F1 progeny (rs, C2-3KO×DBM1Ac-S) (top pictures), N1-3aKO (rr), and their F1 progeny (rs, N1-3aKO×DBM1Ac-S) (middle pictures), as well as C-NKO (rr) and their F1 progeny (rs, C-NKO×DBM1Ac-S) (bottom pictures) were used. Each picture contains 10 backcross families from bottom to top, the first lines 1-5 represent backcross a, and the following lines 6-10 showed backcross b. Bioassays were performed on 30 larvae from each backcross group, and a total of 300 larvae were treated with the diagnostic doses of Cry1Ac protoxin, another 300 untreated larvae were also collected as controls. The genotypes of individuals that were or were not treated with Cry1Ac were detected by the band size of the PCR products with the four primer pairs in PxABCC2 and PxABCC3 or PxAPN1 and PxAPN3a locus (shown in Additional file 1: Table S1). The number of individuals with the genotype of rr is listed on the left side of these figures, and the significant difference of genotype ratio in Cry1Ac treated or not backcross groups were calculated by the Fisher's exact test, and the p values are shown on the right of figures, *p < 0.05, **p < 0.01, ***p < 0.001. c-e Genotype detection of F2 backcross families from crossing the polygenic knockout strains (C2-3KO, N1-3aKO, and C-NKO, respectively) and their corresponding F1 progeny. The gDNA samples of all the individuals from F2 backcross families, that were or were not exposed to Cry1Ac toxin, were extracted for PCR amplification using the four primer pairs in PxABCC2/PxABCC3 or PxAPN1/PxAPN3a loci (shown in Figs. 1c and 2b). PCR products were purified and resolved by 1.5% agarose gel electrophoresis

Discussion
The interactions between insect midgut receptors and the Bt pesticidal proteins influence both the pathogenicity of the bacterium and the ability of the host to evolve resistance. Interactions between resistance alleles are known to play a vital role in the evolution or maintenance of polygenic insecticide resistance in the field [34]. Resistance alleles based on mutations to midgut receptors are well known to affect the pathogenicity of Bt and more recently reports of interactions between these alleles have been noted. An interaction between ABCC2 and cadherin BtR175 in Bt Cry1Ab resistance was proposed in Bombyx mori [37] and similar results were obtained with Cry1Ac resistance in Heliothis virescens and Trichoplusia ni [17,33,36]. ABCC2 cooperated with its paralog ABCC3 triggering > 15,000 and > 8000-fold resistance to Bt Cry1Ac toxin in Helicoverpa armigera and P. xylostella, respectively, when both were knocked out, however, mutations in either of these single genes had little or no effect [43,44]. Other reports observed levels of resistance when either ABCC2 or ABCC3 were knocked out individually [25,28] but much higher levels e The expression level of non-receptor genes in the DBM1Ac-S, C2-3KO, N1-3aKO, C-NKO, and NIL-R strains. Data are presented as mean value ± SEM, n = 5 biologically independent samples with ten larvae per replicate for a-d, and n = 3 biologically independent samples for e, *p < 0.05, **p < 0.01, ***p < 0.001, ns, not significant. One-way ANOVA with Holm-Sidak's test was used for comparison were reported for the double knockouts. In this study, a CRISPR/Cas9-mediated double knockout of ABCC2 and ABCC3 in P. xylostella resulted in a significantly greater level of resistance than was previously observed with individual knockouts [25,27]. Although the individual resistance ratios of the single and double knockouts varied between studies [43,44], there was a consistent pattern of the double knockout resulting in a greatly increased resistance (Fig. 6), indicating functional redundancy between the two ABCC receptors. In addition, we showed that knocking down the two APN receptors also resulted in a significant increase in resistance to Cry1Ac, suggesting that these APN receptors have a redundant function (Fig. 6). The level of Cry1Ac resistance in the quadruple knockout strain was so high that it was impossible to calculate an LC 50 value.
A commonly cited model for Bt toxin mechanism of action is the "sequential binding" model in which toxin monomers are proposed to the first bind to a membrane-anchored receptor increasing the concentration of toxins at the membrane surface. Subsequently, the toxin interacts with a second receptor which induces oligomerization of the protein following a proteolytic cleavage. Insertion of the oligomer into the membrane to form a pore may be facilitated by it further binding to a lipid raft-associated protein [13,15]. In a model in which a defined and obligatory sequential movement of the toxin between different receptors exists, it logically follows that blocking any one of the steps would lead to similar levels of resistance. For example, if movement between an APN and an ABCC receptor is a required step in the mechanism, then removing either the APN Fig. 6 Diagramatic representation of the response of P. xylostella larvae from multiple strains after treatment with Cry1Ac toxin. The dashed lines to which the arrow points represent the pathways with the corresponding receptors changed. Different P. xylostella strains: susceptible DBM1Ac-S strain; near-isogenic Cry1Ac-resistant NIL-R strain; CRISPR/Cas9-mediated single mutant strain with PxAPN1/PxAPN3a/PxABCC2/PxABCC3 deleted, respectively; CRISPR/Cas9-mediated double mutant strain with PxAPN1/PxAPN3a or PxABCC2/PxABCC3 deleted; the hybrid strain with both deletions. RR: Resistance ratio receptor(s) or the ABCC receptor(s) would lead to the same level of resistance. However, if redundant receptors with the same (receptor) function exist in the same insect, this would require multiple receptors to be knocked out in order to have high levels of resistance. In P. xylostella any redundancy between APN and ABCC proteins is most likely to exist in one of two ways. Firstly the two APN proteins substitute for each other at one stage of the sequential pathway and the two ABCC proteins at another stage. In such a scenario knocking out both the APNs, or both the ABCC proteins, would block the pathway and give similarly high levels of resistance. Secondly, all four proteins could independently act as receptors in a non-sequential pathway in which case all four would have to be knocked out to give high levels of resistance. The data presented here support the latter idea that there is full redundancy between the PxAPN1/3a and PxABCC2/3 receptors. There is good evidence from ectopic expression studies that APN and ABCC proteins are alone capable of facilitating Cry protein toxicity [27,36,[51][52][53][54]. Thus, although our data suggest that these four receptors are functionally redundant and independent, they are also consistent with the possibility of some cooperation between them to enhance the activity of the Cry protein.
There is convincing evidence in the literature for synergistic interactions between different receptors within the Bt toxin mechanism of action. For example, Bretschneider et al. [36] demonstrated that ectopic expression of the H. virescens ABCC2 receptor in a non-susceptible cell line could render that line susceptible to Cry1Ab. Ectopic expression of the cadherin receptor from the same insect resulted in little or no change in susceptibility, yet when both receptors were expressed in the same cell, the susceptibility was much greater than the additive effect of the two alone. Further studies using ABCC2 and cadherin receptors from H. armigera, Spodoptera litura, and Spodoptera exigua expressed in cell lines also observed synergistic reactions [55,56] as did studies in which the receptors were expressed in Xenopusoocytes [39] or in Drosophila melanogaster [57]. Work involving both naturally occurring and artificially created mutations in H. armigera ABCC2 and cadherin also showed evidence of synergistic interaction [58]. The mechanistic detail of this synergistic effect remains unclear although Heckel [59] has outlined a number of cisand trans-acting models. A common theme of these synergism models is that some receptors may be better at inducing oligomer formation whereas others may be better at facilitating insertion into the membrane. ABCC2 has been shown to be capable of enabling both of these processes, inducing oligomer formation and its insertion into the membrane [60], the overall effect though could be enhanced via cooperation with a second receptor such as an APN.
The pyramid strategy of transgenic Bt crops effectively controls insect pests by producing two or more toxins that kill the same pest [9]. However, the method is often ineffective in delaying or overcoming resistance because resistance to one toxin is often accompanied by crossresistance to others [61]. In our laboratory bioassays of P. xylostella, we discovered that the knockout strains C2-3KO, N1-3aKO, and C-NKO showed resistance to Cry1Ac and Cry1Ab but no obvious cross-resistance to Cry1Aa (Table 2). Various models describing multiple Cry protein binding sites in P. xylostella have been proposed [62], although these were derived based on resistant strains to which Cry1Aa bound but was not toxic. Although some resistant P. xylostella strains show crossresistance to Cry1Aa, resistance ratios to this toxin are often less than those for Cry1Ab or Cry1Ac [63,64]. Our data are consistent with the primary functional receptor for Cry1Aa being something other than APN1/3a or ABCC2/3, perhaps ALP which is known to be associated with resistance to Cry1Aa expressing Bt strains in other P. xylostella strains [22,27].
In nature, the evolutionary outcome of insect resistance to insecticides is determined by the fitness costs of each resistance allele, and the overall fitness disadvantages are traded off against the advantages of the resistance phenotype [34,65,66]. Since many Cry toxin receptors are physiologically important proteins in the insect, resistance-inducing mutations or loss of expression can alter digestive processes. However, the upregulation of paralogs with no receptor function and transcriptional plasticity could compensate for the functional impairment [67,68]. Our previous studies have shown that the MAPK signaling pathway can trans-regulate the differential expression of multiple midgut genes, including PxABCC1-3 and four PxAPN genes (PxAPN1, PxAPN3a, PxAPN5, and PxAPN6), thereby countering the virulence effect of Bt toxin in P. xylostella without occurring significance fitness costs [22,27,50]. Among them, the up-regulation of PxABCC1, PxAPN5, and PxAPN6 were not associated with Cry1Ac resistance. This up-regulation is speculated to make up for the loss of function of PxABCC2, PxABCC3, PxAPN1, and PxAPN3a and diminish the fitness cost of Cry1Ac resistance. A similar compensation mechanism of TnAPN1 and TnAPN6 was reported in Cry1Ac-resistance T. ni [30]. In this work, we found that the double mutant C2-3KO and N1-3aKO strains exhibited no statistically significant differences in the measured developmental or physiological features, although such a difference was observed in the quadruple knockout strain C-NKO (Fig.  6). One possible explanation for the lack of large fitness costs in the double knockouts is that the insect compensated for the loss of the receptor proteins through the upregulation of receptor paralogs, perhaps via hormonal signaling pathways that have previously been associated with resistance and insect fitness [69,70]. However, when we measured the levels of the non-receptor paralogs PxAPN5/6 and PxABCC1 we found no change in any of the knockout strains. This would indicate that no simple feedback system exists to compensate for the loss of those midgut proteins. It remains unclear why the two double knockout strains had little effect on fitness whereas the quadruple one did, although it could just be that individual knockouts incur small fitness penalties that are amplified when combined. Related work recently reported for H. armigera [58] found that in a strain where resistance had been linked to a mutation in HaABCC2, expression of HaABCC3 was also downregulated, whereas this was not the case in an artificially constructed HaABCC2 knockout. This result and our previous characterization of resistant P. xylostella strains indicate that, unlike lab-induced gene knockouts, naturally evolved resistant strains often present secondary effects that both increase the level of resistance while mitigating the physiological load on the insect.

Conclusion
In summary, our data support the hypothesis that functional redundancy goes beyond ABC transporters to include other, more diverse, midgut receptors. Specifically, knockouts of individuals or pairs of a suite of Bt receptors, including APN1, APN3a, ABCC2, and ABCC3, result in some degree of resistance, while high levels of resistance are only seen when all four proteins are eliminated. Previous studies [25,27,57] have demonstrated that individual members of this group can act as receptors in their own right and so it appears that there is no requirement for multiple receptors to be present. Further research is warranted to address whether the effect of multiple receptors on susceptibility is purely additive or if there are synergistic interactions. The functional redundancy that we have identified among diverse midgut receptors not only supports our overarching hypothesis of a compensatory mode of action of Bt toxins in P. xylostella, but also provides a plausible explanation as to why the resistance mechanisms of this super pest involve a coordinated transcriptional programing of multiple receptor proteins [22,27,50].

Insect strains
The susceptible DBM1Ac-S strain and its near-isogenic Cry1Ac-resistant NIL-R strain of P. xylostella used in this study have been described in detail elsewhere [25,71]. CRISPR/Cas9-mediated double mutant strains (C2-3KO and N1-3aKO) were constructed from the susceptible DBM1Ac-S strain by introducing 37 kb or 17 kb deletions into the PxABCC2/PxABCC3 or PxAPN1/ PxAPN3a gene loci resulting in loss of function of both encoded proteins. A hybrid strain (C-NKO) that contained mutations in all the above four genes was established by crossing the moths from C2-3KO and N1-3aKO strains. All P. xylostella strains were reared at 25°C with 65% relative humidity (RH) and a photoperiod 16:8 (light: dark) on Jing Feng No. 1 cabbage (Brassica oleracea var. capitata), with a 10% honey/ water solution for adults.

Preparation of sgRNA/Cas9 protein mixtures for microinjection
Freshly preblastoderm-stage eggs were laid on dry microscope slides (24 × 50 mm) precoated with fresh cabbage leaf juice to induce the spawning behavior of female adults. The Streptococcus pyogenes Cas9 protein was purchased from Thermo Fisher Scientific (Shanghai, China). About 1-nl mixture of sgRNAs and Cas9 protein were microinjected into individual eggs. The instruments and injection procedures were in line with our previous report [25]. For the double gene knockout, the two sgRNAs (sgRNA1 and 2 or sgRNA3 and 4) and Cas9 protein were simultaneously injected into individual eggs, and the final concentration of each sgRNA and Cas9 protein was both 100 ng/μl. The microinjection process was performed within 2 h, and the injected eggs were immediately incubated at 25°C with 65% RH for hatching.

Identification of double mutant moths mediated by CRISPR/Cas9
According to the arrangement of the PxABCC3/ PxABCC2 and PxAPN3a/PxAPN1 genes in the DBM genome database, two primer pairs were designed to detect the mutagenesis of the two double-gene regions. 3-2-F is located in intron 2 of the PxABCC3 gene as the forward primer and 2-3-R located in exon 3 of the PxABCC2 gene as the reverse primer, a 516 bp genome DNA fragment was amplified with the primers 3-2-F/2-3-R if PxABCC3 and PxABCC2 genes were simultaneously deleted. 3a-1-F is located in exon 1 of the PxAPN3a gene as the forward primer, and 1-13-R located in exon 13 of the PxAPN1 gene as a reverse primer, a 716 bp gDNA fragment was amplified with the primers 3a-12-F/1-13-R if the PxAPN3a and PxAPN1 genes were deleted. To identify homozygous or heterozygous individuals, three other PCR primer pairs were designed between the two contiguous genes. For the identification of the PxABCC3/ PxABCC2 knockout, the PCR primers 3-18-F/3-19-R are respectively located in exons 18 and 19 of the PxABCC3 gene giving a PCR product of 441 bp, the PCR primers 2-12-F/2-13-R are located in exons 12 and 13 of the PxABCC2 gene giving a PCR product of 334 bp, and the PCR primers 3-25-F/ 2-26-R are located at the junction of the PxABCC3 and PxABCC2 genes-the last exon 25 of the PxABCC3 and the last exon 26 of the PxABCC2 genes giving a PCR product of 891 bp. The mutations of the PxAPN3a and PxAPN1 were similarly identified, PCR primers 3a-13-F/ 3a-14-R are located in exons 13 and exon 14 of PxAPN3a giving a 522 bp gDNA fragment, and PCR primers 1-10-F/1-11-R are located in exons 10 and 11 of PxAPN1 gene giving a 332 bp PCR product fragment, an approximately 1381 bp PCR product between the PxAPN3a and PxAPN1 genes was amplified using PCR primers 3a-14-F/1-2-R, which were respectively located in exon 14 of the PxAPN3a gene and exon 2 of the PxAPN1 gene (Additional file 1: Table S1) Finally, the genotypes of the double mutant individuals (PxABCC3/ PxABCC2 and PxAPN3a/PxAPN1) were validated by the PCR profile of the resultant amplicons and DNA sequencing.

Construction and identification of genetic hybridization strain
To generate a moth line with PxABCC3, PxABCC2, PxAPN3a, and PxAPN1 gene mutations, a genetic cross strategy was employed. Taking into account that the PxABCC3/PxABCC2 and PxAPN3a/PxAPN1 genes in P.
xylostella are orthologous to genes located in chromosome 15 or 9 of B. mori, respectively [18], we performed reciprocal crosses between lines C2-3KO (BBaa) and N1-3aKO (bbAA) that produce a heterozygous genotype (AaBb) (the aa represents the homozygous alleles for PxABCC2 and PxABCC3 deletions, and bb shows the homozygous alleles for PxAPN1 and PxAPN3a mutations) (Additional file 1: Fig. S2). Subsequently, G1 heterozygous individuals were sib-crossed to produce G2 progeny and screened with a diagnostic dose (10 mg/L) of Cry1Ac protoxin which could kill all the susceptible larvae. The remaining surviving larvae were reared to pupation, and the gDNA samples of tiny exuviates were isolated from the final fourth-instar larvae to avoid damaging the pupae. Then, the genotype of individuals was detected by the profile of PCR amplicons using the corresponding four primer pairs between the PxABCC3/ PxABCC2 and PxAPN3a/PxAPN1 loci (Additional file 1: Table S1). Finally, G2 individuals containing PxABCC3, PxABCC2, PxAPN3a, and PxAPN1 mutations were selected and sib-mated to establishing a stable homozygous mutant strain designated C-NKO in G3.

Toxins and bioassays
The Bt protoxins were extracted and purified from Btk strain HD-73 via the isoelectric point precipitation method as described previously [72]. All the Bt protoxins and trypsin-activated Bt toxins were quantified by Bradford's method using bovine serum albumin (BSA) as a standard.
A 72 h leaf-disc-dip bioassay was performed to test the response of the DBM1Ac-S, C2-3KO, N1-3aKO, and C-NKO strains to three Bt Cry1A toxins (Cry1Aa, Cry1Ab, Cry1Ac) as reported elsewhere [71]. Ten third-instar larvae were respectively inoculated in seven gradient concentrations of each toxin, and each treatment was repeated three times. Larval mortality was recorded in each strain and the control mortality limit was set to 5%. After bioassays, the LC 50 values (the concentration that killed 50% tested larval) and 95% CL (95% fiducial limits of LC 50 values) values were calculated by Probit analysis using POLO Plus 2.0 software (LeOra Software, Berkeley, CA, USA). We examined the resistance ratio for each toxin as the LC 50 values of each knockout strain divided by the LC 50 of the susceptible DBM1Ac-S strain. The LC 50 values of pairwise comparison were perceived as significantly different if their 95% CL did not overlap.

Inheritance and allelic complementation analysis
To ascertain the inheritance of Cry1Ac resistance in both polygenic knockout strains, we used the same interstrain crossing strategy between two out of the DBM1Ac-S, C2-3KO, N1-3aKO, and C-NKO strains, the mortality of these four strains (50 larvae from each strain), and their F1 progeny (100 larvae from each F1 offspring) were examined with a diagnostic dose (10 mg/ L) of Cry1Ac protoxin, which could kill 100% heterozygous F1 larvae. Then, the dominance degree (h) of Cry1Ac resistance was calculated with the formula reported elsewhere [73], with the value of h ranging from 0 (completely recessive) to 1 (completely dominant).
To confirm if the locus responsible for resistance to Bt Cry1Ac toxin changed among different strains, interstrain allelic complementation tests were performed. Through bioassays, we examined the offspring of singlepair crossed between the DBM1Ac-S, C2-3KO, N1-3aKO, and C-NKO strains. If the moths from two resistant strains were crossed and the recessive resistance alleles were located in different loci, their progeny will exhibit susceptibility to Cry1Ac toxin, and if the recessive resistance alleles are located in the same locus, their progeny will be resistant to Cry1Ac toxin.

Genetic linkage analysis
For genetic linkage analysis of double gene knockout strains (C2-3KO, N1-3aKO) and hybrid strain (C-NKO) with Cry1Ac resistance, a single-pair cross between a male from the knockout strains (C2-3KO, N1-3aKO, C-NKO) and a female from the susceptible DBM1Ac-S strain was performed to produce F1 progeny. Subsequently, five reciprocal crosses between an F1 progeny and a knockout moth were performed to generate ten F2 backcross families. A diagnostic dose of Cry1Ac protoxin that could kill all the F1 heterozygous larvae (20 mg/L for C2-3KO and N1-3aKO strains or 4000 mg/L for C-NKO strain) was used to confirm susceptibility by bioassay. For F2 progenies from each backcross family, 30 third-instar larvae were respectively reared on control cabbage or treated cabbage (20 mg/L or 4000 mg/L Cry1Ac protoxin smeared on the leaves) for 72 h. The mortality was recorded and the gDNA sample from each surviving larva was extracted for genotyping detection by PCR profile using the primers mentioned above for detecting the heterozygous or homozygous mutations in PxABCC3/PxABCC2 or PxAPN3a/PxAPN1 (Additional file 1: Table S1).

Fitness cost analysis
We compared a series of physiological parameters of P. xylostella polygenic knockout strains to analyze the fitness cost caused by CRISPR/Cas9 and the genetic hybridization method. Biological parameters including pupation rate, pupation duration, pupal weight and hatching percentage were recorded, and larvae from the susceptible DBM1Ac-S were used as a control. Ten second-instar larvae from each strain were reared on fresh cabbage leaves without exposure to any Bt Cry toxin and the experiments were repeated in quintuplicate. The statistical significance difference in biological parameters between the control and polygenic knockout strains was determined by one-way ANOVA with Holm-Sidak's test (overall significance level > 0.05).

Quantitative PCR (qPCR) analysis
The transcript levels of non-receptor genes were quantified by real-time qPCR performed in a QuantStudio 3 Real-Time PCR System (Applied Biosystems). Genespecific primers for PxAPN5, PxAPN6, and PxABCC1 were used (Additional file 1: Table S1) in qPCR reactions with 2.5×SYBR Green MasterMix Kit (TIANGEN) as described in detail elsewhere [22,27]. The qPCR program consisted of an original denaturation step for 6 min at 95°C, subsequently, 40 cycles of denaturation for 30 s at 95°C, annealing for 30 s at 53°C for PxAPN6, 55°C for PxABCC1 and PxAPN5, followed by an extension for 35 s at 72°C. Relative quantification was calculated by utilizing the 2 −ΔΔCt method and standardized to the ribosomal protein L32 gene (GenBank accession no. AB180441). Three biological repetitions and four technical replicates were performed for each sample.
Additional file 1: Figure S1. Diagram of the crossing strategy to obtain homozygous double knockout strains in P. xylostella. Figure S2. Diagram of the crossing strategy used for generation of the genetic hybrid strain. Table S1. Primers used in this study. Table S2. Toxicity to Cry1Ac toxin in larvae from the polygenic knockout strains and their F1 progeny.
Additional file 2. The raw data of all the figures and statistical analyses in this study.