Next Article in Journal
An Odorant Receptor from the Proboscis of the Cotton Bollworm Helicoverpa armigera (Lepidoptera: Noctuidae) Narrowly Tuned to Indole
Next Article in Special Issue
Functional Characterization of the Nuclear Receptor Gene SaE75 in the Grain Aphid, Sitobion avenae
Previous Article in Journal
New Fossil Xyelidae (Hymenoptera: Symphyta) from the Mesozoic of Northeastern China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

CRISPR/Cas9-Mediated Mutagenesis of Abdominal-A and Ultrabithorax in the Asian Corn Borer, Ostrinia furnacalis

1
College of Plant Protection, Shenyang Agricultural University, Shenyang 110866, China
2
State Key Laboratory of Cotton Biology, Key Laboratory of Plant Stress Biology, School of Life Sciences, Henan University, Kaifeng 475004, China
3
Department of Entomology, University of Kentucky, Lexington, KY 40546, USA
4
Key Laboratory of Insect Developmental and Evolutionary Biology, Institute of Plant Physiology and Ecology, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai 200032, China
*
Author to whom correspondence should be addressed.
Submission received: 14 March 2022 / Revised: 2 April 2022 / Accepted: 9 April 2022 / Published: 13 April 2022
(This article belongs to the Special Issue Opportunities and Challenges in Insect Functional Genomics)

Abstract

:

Simple Summary

Homeotic genes encode transcription factors that coordinated the anatomical structure formation during the early embryonic development of organisms. In this study, we functionally characterized two homeotic genes, Abdominal-A (Abd-A) and Ultrabithorax (Ubx), in the Asian corn borer, Ostrinia furnacalis (a maize pest that has devastated the Asia-Pacific region) by using a CRISPR/Cas9 genome editing system. Our results show that the mutagenesis of OfAbd-A and OfUbx led to severe morphological defects in O. furnacalis, which included fused segments and segmental twist during the larval stage, and hollowed and incision-like segments during the pupal stage in OfAbd-A mutants, as well as defects in the wing-pad development in pupal and adult OfUbx mutants. Overall, knocking out Abd-A and Ubx in O. furnacalis resulted in the embryonic lethality to, and pleiotropic impact on, other homeotic genes. This study not only confirms the conserved body planning functions in OfAbd-A and OfUbx, but it also strengthens the control implications of these homeotic genes for lepidopteran pests.

Abstract

(1) Background: Abdominal-A (Abd-A) and Ultrabithorax (Ubx) are homeotic genes that determine the identity and morphology of the thorax and abdomen in insects. The Asian corn borer, Ostrinia furnacalis (Guenée) (Lepidoptera: Pyralidae), is a devastating maize pest throughout Asia, the Western Pacific, and Australia. Building on previous knowledge, we hypothesized that the knockout of Abd-A and Ubx would disrupt the abdominal body planning in O. furnacalis. (2) Methods: CRISPR/Cas9-targeted mutagenesis was employed to decipher the functions of these homeotic genes. (3) Results: Knockout insects demonstrated classical homeotic transformations. Specifically, the mutagenesis of OfAbd-A resulted in: (1) Fused segments and segmental twist during the larval stage; (2) Embryonic lethality; and (3) The pleiotropic upregulation of other homeotic genes, including Lab, Pd, Dfd, Antp, and Abd-B. The mutagenesis of OfUbx led to: (1) Severe defects in the wing pads, which limited the ability of the adults to fly and mate; (2) Female sterility; and (3) The pleiotropic upregulation of other homeotic genes, including Dfd, Abd-B, and Wnt1. (4) Conclusions: These combined results not only support our hypothesis, but they also strengthen the potential of using homeotic genes as molecular targets for the genetic control of this global insect pest.

1. Introduction

Asian corn borer, Ostrinia furnacalis (Guenée) (Lepidoptera: Pyralidae), is a devastating maize pest throughout Asia, Western Pacific, and Australia [1,2]. Damage caused by O. furnacalis results in a loss of approximately 9 million tons of corn annually [3]. Traditional pest management strategies for O. furnacalis and European corn borer, Ostrinia nubilalis (Hübner), rely on chemical pesticides; biocontrol agents, such as the parasitoid wasp, Trichogramma ostriniae [4,5,6,7]; and biopesticides, such as Bt toxins [8,9,10,11]. Bt resistance, however, has developed in O. furnacalis in the laboratory [12], and in O. nubilalis in the field [13,14,15]. Multiple strategies have been developed to manage Bt resistance, including refuge, stacked traits, and gene editing techniques [16].
Hox genes, which are structurally conserved with identical Homeodomain motifs across taxa [17], dictate the anterior–posterior body identity and morphological features in insects [18,19]. Genes in the Hox family are subcategorized into two clusters: the Antennapedia complex (ANT-C) and the Bithorax complex (BX-C) [20]. There are a total of 8 Hox genes: Labial (lab), Antennapedia (Antp), proboscipedia (pb), Deformed (Dfd), and Sex combs reduced (Scr) belong to the ANT-C cluster, while Ultrabithorax (Ubx), Abdominal-A (Abd-A), and Abdominal-B (Abd-B) belong to the BX-C cluster [21]. Hox genes are typically expressed at embryonic stage to specify the arrangement of thorax and abdominal segments [17]. Mechanisms that underlie how Hox genes work, individually or collaboratively, are still not fully understood [22,23].
Abd-A is essential in the formation of embryonic and abdominal segments in insects [24]. In Drosophila, Abd-A plays integral roles in suppressing limb development and in specifying abdominal segments [25,26]. Mutation of Abd-A results in homeotic transformation of parasegments 5 and 6 during the embryonic stage, which eventually leads to neonatal mortality and defects in head and thorax [17]. Abd-A is also involved in the formation of female internal genitalia [27,28], nervous system, fat body [29], and midgut [30]; aorta and heart cardioblast differentiation [31]; and ectopic pigmentation [32].
Ubx acts as a genetic switch to modify specific morphological features in thoracic region in insects [17,33,34,35]. In Drosophila, Ubx promotes the formation of halteres through expression in the meso and metathorax [36,37]. The size and identity of the appendages is modulated by the expression of Ubx in larvae [33,35,36,37,38,39,40,41]. In brown planthopper, Nilaparvata lugens (Stål), Ubx is the key regulator of the transformation between short and long wing forms [22]. In lepidopterans, Ubx is expressed in the metathoracic segment in Bombyx mori and it regulates wing development [38]. Silencing of Ubx and Abd-A leads to suspended limb development in Drosophila, and in many lepidopteran insects, through the suppression of Distal-less (Dll) gene expression [27,39,40,41].
The overall goal of this study was to gain a better understanding of body planning and organ development in O. furnacalis. Building on previous knowledge and preliminary research, we hypothesized that Abd-A and Ubx are involved in body patterning of thoracic and abdominal segments in O. furnacalis. To examine this overarching hypothesis, we: (1) Established a CRISPR/Cas9 genome editing system in O. furnacalis; and (2) Functionally characterized Abd-A and Ubx in this devastating corn pest.

2. Materials and Methods

2.1. Phylogenetic Analysis of OfAbd-A and OfUbx Genes

2.1.1. Phylogenetic Tree

Abd-A protein sequences that were used in the phylogenetic analysis included Apis mellifera (GenBank accession number: XP_016772643.1), Aedes aegypti (NP_001345961.1), Tribolium castaneum (NP_001034518.1), Bombyx mori (NP_001166808.1), Drosophila melanogaster (NP_001247145.1), Danaus plexippus plexippus (OWR52832.1), Pieris rapae (XP_022113611.1), Bicyclus anynana (XP_023954095.1), Vanessa tameamea (XP_026485495.1), Papilio xuthus (XP_013173869.1), Dendrolimus punctatus (AQM32554.1), Plutella xylostella (XP_011569267.1), Spodoptera litura (XP_022825998.1), Trichoplusia ni (XP_026731481.1), Myrmica rubra (AAK06846.2), Papilio polytes (XP_013139462.1), and Galleria mellonella (XP_026759670.1). Ubx protein sequences included Drosophila melanogaster (NP_536752.1), Bombyx mori (NP_001107632.1), Apis mellifera (NP_001162171.1), Ceratitis capitata (XP_004524337.1), Spodoptera litura (XP_022837384.1), Trichoplusia ni (XP_026735287.1), Bactrocera dorsalis (XP_011203163.1), Plutella xylostella (NP_001303599.1), Tribolium castaneum (NP_001034497.1), Anopheles gambiae (AAC31942.1), Sogatella furcifera (ATW63192.1), Papilio xuthus (XP_013173873.1), Pieris rapae (XP_022113620.1), Helicoverpa armigera (XP_021195138.1), Dendrolimus punctatus (AQM32553.1), Orchesella cincta (CDI44538.1), Papilio machaon (XP_014359749.1), Bicyclus anynana (XP_023954092.1), Biston betularia (ADO33070.2), and Galleria mellonella (XP_026759674.1). Sequence alignment was constructed by using the maximum likelihood method based on CLUSTAL W2 [42,43]. All of the ambiguous positions were removed for each sequence pair. The neighbor-joining method was used to create the tree from 18 available Abd-A sequences and 21 Ubx sequences. The percentages of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) are shown next to the branches. The evolutionary distances were computed using the Poisson correction method and are displayed as the number of amino acid substitutions per site. The evolutionary history was inferred by using the neighbor-joining method and MEGA-7 with 1000 bootstrap replicates [44].

2.1.2. Multiple Alignment

Amino acid sequences encoded by Abd-A from Drosophila melanogaster (NP_001247145), Plutella xylostella (XP_011569267), Bombyx mori (ACD10794), Tribolium castaneum (NP_001034518), and Apis mellifera (XP_016772643), as well as the predicted amino acid sequence of Abd-A from Ostrinia furnacalis, and similarly, amino acid sequences of Ubx derived from Drosophila melanogaster (NP_536752), Plutella xylostella (NP_001303599), Bombyx mori (NP_001107632), Tribolium castaneum (NP_001034497), Apis mellifera (NP_001162171), and O. furnacalis were analyzed by Clustal Omega, which is a multiple sequence alignment program that uses seeded guide trees and HMM profile–profile techniques to generate alignments between three or more sequences [45].

2.2. Temporal and Spatial Expression Profiles of OfAbd-A and OfUbx

2.2.1. Sample Collection

To investigate the spatial expressions of OfAbd-A and OfUbx genes, total RNA was isolated from eggs, from larvae during the first day of each instar, from prepupae (PP), from pupae (P), and from adults (A), by using Trizol reagent (Invitrogen, Carlsbad, CA, USA), and it was treated with RNase-free DNase I (Ambion, Austin, TX, USA), according to the manufacturer’s protocols. For tissue-specific analysis, fifth-instar larvae were dissected three days after ecdysis to obtain tissue from head, midgut, foregut, fat body, epidermis, testes, and ovaries, which were kept in liquid nitrogen for sample collection and total RNA extraction. To address the question of how Hox proteins are regulated by the mutagenesis of OfAbd-A and OfUbx, we collected pupal O. furnacalis for the total RNA extraction to analyze the relative transcript levels of the genes in the Hox family in OfAbd-A and OfUbx mutants.

2.2.2. Reverse Transcription Quantitative Real-Time PCR (RT-qPCR) Analysis

cDNAs were synthesized using 1 μg total RNA from each developmental stage as templates by using an Omniscript reverse transcriptase kit (Qiagen, Hilden, Germany) in a 20 μL reaction mixture. The primers used for the RT-qPCR are listed in Table S1. RT-qPCR was performed on an Eppendorf Real-time PCR System using the following conditions: a 2 min denaturing cycle at 95 °C, 35 cycles of 1 min at 95 °C, 30 s at 55 °C, and 30 s at 72 °C, followed by a final extension at 72 °C for 10 min.

2.3. CRISPR/Cas9-Targeted Mutagenesis in O. furnacalis

2.3.1. Molecular Cloning and Target Selection

Total RNA was isolated from fifth-instar larvae using Trizol Reagent (Invitrogen, Carlsbad, USA) and was treated with RNase-free DNase I (Ambion, Austin, TX, USA), according to the manufacturer’s protocol. cDNA was synthesized with Omniscript reverse transcriptase kit (Qiagen, Hilden, Germany), using manufacturer’s instructions. Putative OfAbd-A and OfUbx genes were identified using the NCBI Blast system. OfAbd-A cDNA fragments were amplified by PCR with the following pair of primers: forward, 5′-ATGGCAGCGGCTGCCCAGTT-3′; and reverse, 5′-TTACGTGGGGACTTTGTTCA-3′. OfUbx cDNA fragments were amplified by PCR with the following pair of primers: forward, 5′-ATGAACTCCTACTTTGAGCAGGGT-3′; and reverse, 5′-TTACGTGGGGACTTTGTTCA -3′. PCR was carried out using KOD -Plus- polymerase (TOYOBO, Osaka, Japan) under the following conditions: 98 °C for 2 min, followed by 30 cycles at 98 °C for 30 s, 55 °C for 30 s, 68 °C for 1 min, and an elongation phase at 68 °C for 10 min. Amplified products were sequenced after cloning into a PJET1.2-T vector (Fermentas, Burlington, ON, Canada). The primers are listed in Table S1.

2.3.2. Synthesis of Cas9 mRNA and sgRNAs

We selected one 23 bp sgRNA to target the OfAbd-A genome locus on the third exon of the common region in the four spliced variants, and two 23 bp sgRNAs targeting the OfUbx genome locus in exon 1 of the common region for the two spliced variants. The sgRNA was subcloned into a 500 bp linearized CloneJet PJET1.2-T vector (Thermo Fisher, Waltham, MA, USA), upstream of the protospacer adjacent motif (PAM) sequence to allow sgRNA expression under the control of the T7 promoter. The sgRNA was synthesized in vitro with a MEGAScript T7 kit (Ambion, Austin, TX, USA), according to the manufacturer’s instructions. Cas9 mRNA was synthesized in vitro using the mMESSAGE T7 Kit (Ambion, Austin, TX, USA) and a PTD1-T7-Cas9 vector as the template, according to the manufacturer’s instructions.

2.3.3. Embryo Microinjection

A laboratory strain (Shanghai) of O. furnacalis larvae was reared with an artificial diet [46]. Insects were kept at 25 °C with 80% relative humidity and a 16:8 light:dark photoperiod [47]. Adults of O. furnacalis were maintained in transparent plastic bags and were fed using cotton balls soaked in sugar water (10% honey in distilled water) to allow them to lay eggs. Eggs were collected from the plastic bags, which were cut into pieces for easier collection, and were arranged in a row on sterilized cover slips, as described previously [48]. Eggs were injected on the lateral side with 10 nL of a mixture containing 300 ng/μL Cas9 mRNA and 150 ng/μL sgRNA, within 1 h of oviposition. After injection, eggs were incubated in a humidified chamber at 25 °C for 4 days until hatching.

2.4. Genomic DNA Extraction to Identify Successful Mutants

The genomic DNA of O. furnacalis larvae was extracted by using phenol:chloroform and an isopropanol precipitation extraction. Specifically, newly hatched larvae were collected and incubated with proteinase K, followed by DNA purification and RNaseA treatment. PCR was carried out to identify OfAbd-A and OfUbx mutant alleles by using primers spanning the target sites in OfAbd-A and OfUbx (Table S1). The PCR conditions were as follows: 98 °C for 2 min, followed by 35 cycles of 94 °C for 10 s, 55 °C for 30 s, and 72 °C for 1 min, followed by a final extension period of 72 °C for 10 min. The PCR products were cloned into pJET1.2-T vectors (Fermentas, Burlington, ON, Canada) and were sent for sequencing. The mutants were photographed with a digital stereoscope (Nikon AZ100).
To detect OfAbd-A and OfUbx mutants, RT-qPCR was carried out using gene-specific primers for Abd-A (forward, 5′-CGGCAAACTTACACGAGGTT-3′; and reverse, 5′-TCCTGCTCCTCTCTCTCTCG-3′) to amplify a 221 bp fragment. RT-qPCR reactions were carried out by using gene-specific primers for the Ubx gene (forward, 5′-CCACACCTTCTACCCTTGGA-3′; and reverse, 5′- TCATCCTCCGATTCTGGAAC-3′) to amplify a 221 bp fragment. Ofactin was used as an internal control for the RT-qPCR [47]. The primers are listed in Table S1.

2.5. Analysis of Hatch Rate after Mutagenesis

We analyzed the hatching rate of eggs injected with different concentrations of OfAbd-A sgRNA, and the hatching rate of eggs produced by adult OfUbx mutants. When OfUbx mutants reach adult stage, single-paired adults of different combinations, including mutant male with mutant female, mutant female with WT male, mutant male with WT female, and GFP male with GFP female, as control, were collected to examine hatching rate of G1. Each combination was set in five replicates. All adults were placed in a transparent plastic bag supplemented with cotton balls soaked in 10% honey water to allow them to lay eggs. After five days, we collected eggs from each bag and allowed them to hatch in order to analyze hatching rate of each combination.

2.6. Statistical Analysis

RT-qPCR results were collected and presented in figures by using GraphPad Prism 7.0 software (GraphPad, San Diego, CA, USA). Data were analyzed using an unpaired Student’s t test (SPSS Statistics 25.0 software, IBM, Armonk, NY, USA). Probability values of less than 0.05 were considered significant. Data are presented as means with SEMs.

3. Results

3.1. Phylogenetic Analysis of OfAbd-A and OfUbx

Hox genes in both Drosophila and other invertebrates perform the same overall function of body organization along the anterior–posterior axis [49]. In Drosophila, bithorax complex (BX-C) includes Ubx, Abd-A, and Abd-B genes, which implies that Ubx and Abd-A have a close relationship [50]. We analyzed the genomic sequence of the Abd-A and Ubx genes in O. furnacalis and performed multiple alignments to carry out a sequence analysis. The genomic sequence of OfAbd-A is 41,431 base pairs (bps) in length, and it has three exons with four alternative splicing variants. OfAbd-A mRNA contains four ORFs that are 1050, 1050, 1035, and 1023 bp in length, which encode proteins of lengths of 349, 349, 344, and 340 amino acids (aa), respectively. Comparably, length of Abd-A in Drosophila melanogaster is 990 bp, which encodes a 330aa protein, which exhibits a highly conserved domain with OfAbd-A. Among other moth species, open reading frame (ORF) of Abd-A in Bombyx mori (NP_001166808.1) is 1056 bp, which encodes a protein of 351aa, while, in Plutella xylostella (XP_011569267), 1062 bp Abd-A encodes a protein of 352aa. Both genes are slightly larger than the O. furnacalis Abd-A gene. The genomic sequence of OfUbx is 123,000 bp in length and it has two exons, with two ORFs of 774 and 762 bp, which encode 257 and 254aa proteins, respectively. Ubx encodes a 389aa protein in D. melanogaster (NP_536752.1) and a 254aa protein in B. mori (NP_001107632.1). On the basis of our phylogenetic analysis, OfAbd-A is clustered with homologs in other Lepidopteran species, such as Plutella xylostella, Bombyx mori, and Papilio Xuthus (Figure 1A), while OfUbx is grouped with homologs in Dendrolimus punctatus and Biston betularia (Figure 1B). Our sequence alignment analysis shows that OfAbd-A and OfUbx both contain a Homeodomain and the respective Abdominal and Ultrabithorax domains, which demonstrates a structural conservation with their homologues in other insects (Figure 2A,B).

3.2. Temporal–Spatial Distribution of OfAbd-A and OfUbx

To investigate the transcript changes of OfAbd-A and OfUbx during different developmental stages in Ostrinia furnacalis, we collected eggs, larval instars (first–fifth), from hatching to wandering stages, pupae, and adults of both females and males for the total RNA extraction and RT-qPCR analysis. The results of the RT-qPCR show that the OfAbd-A gene is highly expressed in the egg, wandering, pupal, and female adult stages (Figure 3A), whereas the OfUbx gene is highly expressed in the pupal and adult stages of both females and males (Figure 3B).
To address the question of whether the transcript levels of OfAbd-A and OfUbx are tissue specific, we investigated the expression levels of these two genes from different tissue regions in O. furnacalis larvae by using RT-qPCR. Tissue from head, foregut (FG), midgut (MG), fat body (FB), epidermis (EPI), testes in male larvae (TE), and ovaries in female larvae (OV) were dissected using three-day-old fifth-instar larvae (L5D3). Results from RT-qPCR analysis show that OfAbd-A is highly expressed in the epidermis and ovaries (Figure 3C), while OfUbx is highly expressed in fat body, epidermis, and ovaries (Figure 3D).

3.3. CRISPR/Cas9-Mediated Mutagenesis in O. furnacalis

We used CRISPR/Cas9 genome editing system to knock out Abd-A and Ubx genes in O. furnacalis [32]. There are four isoforms of OfAbd-A, and two isoforms of OfUbx (Figure 4A and Figure 5A). To identify the mutated alleles, we extracted and sequenced the genomic DNA from larvae with mutant phenotypes for OfAbd-A, and from pupae with mutant phenotypes for OfUbx. The genome sequencing showed the successful deletion of sequences within a single target site in the OfAbd-A gene (Figure 4C), and within two target sites in the OfUbx gene (Figure 5C).

3.4. Phenotypic Impacts of OfAbd-A and OfUbx Mutagenesis

3.4.1. Morphological Impacts

The OfAbd-A knockout resulted in twisted abdominal segments in the larvae via the abnormal combination of adjacent segments from A2 to A7 (Figure 6A). In the pupal stage, mutagenesis of OfAbd-A led to hollowed and incision-like segments at abdominal segments, A2 to A7 (Figure 6B). OfUbx gene mutagenesis caused abnormal folding of wing in pupal stage, which hinders the complete covering of the thoracic region (Figure 7A). Mutagenesis of OfUbx also led to the disruption of wing development in pupae. Moreover, severe mutagenesis prevented the adults from spreading their wings to fly, which thus limited their ability to mate and reproduce (Figure 7B).

3.4.2. Physiological Impacts

To investigate the physiological changes induced by OfAbd-A mutation in O. furnacalis, we analyzed the hatching rate of larvae post OfAbd-A injection by using GFP as a control (Table 1). The hatching rate of the GFP control larvae was approximately 62%, while that of the embryos injected with OfAbd-A Cas9/sgRNA was approximately 30%, which suggests that the decrease in hatching rate was specifically caused by OfAbd-A mutagenesis. To further clarify the relationship between the hatching rate and OfAbd-A mutagenesis, O. furnacalis eggs were injected with different doses of Cas9 and Abd-A sgRNA. Results show that, when eggs were injected with Cas9/sgRNA at a concentration of 300 ng/μL, the hatching rate of the larvae decreased to 21% with a mutation rate of approximately 56% (Table 1). When the injection concentration of Cas9/sgRNA was decreased to 150 ng/μL, the hatching rate of the larvae reached 33%, and the mutation rate decreased to approximately 38%. Pupation rate of the treatment group was also lower than that of the GFP control group (Table 1). Our results demonstrate that mutagenesis of the OfAbd-A gene leads to embryonic lethality and affects larval development (Figure 8).
Similarly, higher concentrations of OfUbx sgRNA and Cas9 mRNA induced a higher mutation rate. Only 50% of the pupal OfUbx mutants were able to emerge to the adult stage (Table 1), compared to the emergence rate of GFP Cas9/sgRNA-treated pupae, which was 71%. We also observed that the adult male OfUbx mutants were unable to mate with either wild-type or mutant females, and that female OfUbx mutants could not mate with males from either group. On the basis of our observations, OfUbx mutation can induce sterility in both sexes of adult O. furnacalis (Figure 8).

3.4.3. Pleiotropic Impacts

To address the question of how Hox proteins are regulated by the mutagenesis of OfAbd-A and OfUbx, we analyzed the relative expression levels of the other Hox genes within the OfAbd-A and OfUbx mutants in the pupal stage. OfAbd-A mutagenesis gave rise to the pleiotropic upregulation of Labial (Lab), Proboscipedia (Pb), Deformed (Dfd), Antennapedia (Antp), and Abd-B (Figure S1A), while OfUbx mutagenesis resulted in the upregulation of Dfd, Abd-B, and Wingless Integrated family member 1 (Wnt1) (Figure S1B). The relative transcript level of Dfd increased dramatically following OfUbx and OfAbd-A mutagenesis. Our results show no significant transcript variation in Scr in response to OfUbx and OfAbd-A mutagenesis, while Antp was highly upregulated after OfAbd-A was knocked out (Figure S1).

4. Discussion

In insects, homeotic genes, such as Abd-A and Ubx, specify the distinct identities of body segments [51,52]. We examined the transcription levels of OfAbd-A and OfUbx across O. furnacalis developmental stages to generate a temporal profile of expression for these two genes. OfAbd-A is highly expressed in eggs, wandering-stage larvae, pupae, and female adults (Figure 3A). In other insect species, Abd-A has been reported to initiate its function of segmental identity determination during embryogenesis [53]. OfUbx is highly expressed in pupal and adult stages (Figure 3B). This result is in accordance with characterized function of Ubx in rice planthoppers, where Ubx is expressed in both forewing and hindwing, and it acts as the key regulator for the switch between long and short wing forms in adult stages of both sexes [22]. In addition, Ubx was reported to regulate the development of A1 and T3 legs during adult stage in Oncopeltus fasciatus [52].
Spatial expression pattern of Hox family genes is relatively well understood in Drosophila [19,54] and other insects [55,56], especially in the thoracic and abdominal segments. However, studies on tissue-specific characterization in whole larvae are comparably lacking. We investigated the expression of OfAbd-A and OfUbx in seven larval tissue types in O. furnacalis (Figure 3C,D). Our results show that the relative expression of OfAbd-A is dramatically higher in epidermis, and is comparably high in ovary tissue, compared to other body regions (Figure 3C), which is in accordance with previous reports that show its expression in the epidermal and neural cells in Drosophila [19,26]. The OfUbx gene is highly expressed within epidermis, fat body, and ovary tissues (Figure 3D), which is consistent with previous studies that show that Ubx is involved in gut and muscle development in Drosophila and that it has distinct identities in segmental and structural identification [26,57]. The extensive expression of Ubx was examined in Drosophila in a contiguous region from T3 to A6/A7 in abdomen and gonads, alongside the expression in the fat body [29]. The expression of Ubx was also reported to extend to cover anterior and posterior body regions, and to suppress limb and/or wing development [58,59,60].
CRISPR/Cas9 is a genome editing tool that has recently been exploited as an alternative for pest management [61,62,63]. In this study, we hypothesized that OfAbd-A and OfUbx genes modulate the morphological identity and development of thorax and abdomen in O. furnacalis. The first gene editing system in O. furnacalis was established by You et al. (2018) to characterize the functional role of Ago1 in pigmentation [48], which was then followed by precise gene manipulation targeting yellow [46]. In addition to Ago1 and yellow, research regarding the functional role of Bt toxin receptors, such as ABCC2 and cadherin in O. furnacalis, has been conducted by using gene editing systems [47,64]. By using the CRISPR/Cas9 genome editing system, we successfully induced deletion mutations of 59 and 250 bp into OfAbd-A, respectively (Figure 4), and deletion mutations of 9 and 82 bp into OfUbx, respectively (Figure 5), in O. furnacalis. After knocking out OfAbd-A, the abnormal development of the abdominal segments was observed in larvae (Figure 6), while the abnormal folding of the wings in pupae and adults was observed after OfUbx mutagenesis (Figure 7A,B). In the latter case, adults with severely malformed wings were unable to fly, which restricted their ability to mate and reproduce.
Abd-A regulates wing disc cuticle protein genes that control larval-to-pupal metamorphosis in Bombyx mori, and silencing of Abd-A has demonstrated its role in patterning the third to sixth abdominal segments during embryonic development [41]. In Spodoptera litura, Abd-A is essential for larval segmentation, and mutation can lead to ectopic pigmentation during embryonic development [32]. Phenotypic impacts of Abd-A mutants observed in O. furnacalis (Figure 6) are consistent with McGinnis et al. (1992), which are that Abd-A modulates the development of posterior abdominal segment 1 to anterior abdominal segment 7 in Drosophila [19].
Ubx carries similar functions for wing development and scale morphology in different insect species [38,65]. It modulates the transformation of wings to halteres in third thoracic segment in Dipterans through the repression of wing-specific genes [65]. Ubx, in beetle, Tribolium, regulates the development of elytra in larvae, while RNAi on both Ubx and Abd-A results in the formation of elytron primordia in all abdominal segments [66]. The short-and-long-wing switch in rice planthoppers is also regulated by Ubx gene, and nutritional status of the brown planthopper affects the expression of Ubx, thereby influencing the development between short and long wing forms [22]. This observation is in accordance with our result that shows that knockdown of Ubx disrupts the wing development in O. furnacalis (Figure 7). In our study, knockout of OfUbx gene induced disabled mating behavior and decreased the fecundity of adults (Figure 7B). Our research on OfUbx mutagenesis verifies the functional role of Ubx gene in mediating wing development, from pupal to the adult stages, in O. furnacalis (Figure 7A).
Hox genes play important roles in thoracic and abdominal segmentation, but the genetic relationship among genes in the Hox family is still controversial [1]. Hox genes encode functional transcription factors that modulate distinct identities of thoracic and segmental structures through interactions with a series of downstream genes [67]. To further characterize the role of OfAbd-A and OfUbx in the Hox family, we examined the relative transcript level of Hox genes after mutagenesis (Figure S1). Abd-A control the expression of Lab, which is localized in gut endothelial cells in the visceral mesoderm [20]. Given that the relative transcript level of Lab dramatically increased after OfAbd-A mutagenesis in the pupal stage, the role of Lab may be to compensate for the lost expression of OfAbd-A, and to replicate its function (Figure S1). In addition, lab was reported to specify the intercalary and posterior head segments, while pb, together with Scr, regulates the development of proboscis [68]. Dfd and Scr act downstream of Abd-A to define antennal identity in head segments [69]. Dfd is responsible for the development of eye-antennal discs [70]. In our study, increased expression pattern of Dfd in the OfAbd-A mutants agrees with the report that states that Dfd is involved in the homeotic transformation of thoracic epidermis formation, in addition to being responsible for the ventral and labial epidermis formation on head [49]. As Antp is involved in the development of thoracic segments in B. mori [71], we also examined its relative expression level in OfAbd-A and OfUbx mutants (Figure S1). Our study shows that Antp was highly expressed after OfAbd-A mutation, which suggests that it may act as an essential factor in controlling thoracic and abdominal development (Figure S1A). Abd-A could activate the wingless gene to induce the extra formation of the abdominal segment [72]. In our study, the expression of Abd-B was dramatically increased, while the expression of Wnt1 was significantly decreased in the OfAbd-A mutant (Figure S1A), which further confirms the role of Abd-B, which suppresses the wingless gene through suppression of Abd-A in the A7 region for the elimination of abdomen in Drosophila [72].
Analysis of the downstream gene in OfUbx mutant found that Dfd was upregulated after OfUbx was knocked out (Figure S1), which suggests that Dfd could be suppressed by Ubx. A previous study on Drosophila and B. mori shows that Ubx regulates a series of downstream genes, such as Scr and wingless [38]. Our result found that the expression of Abd-A was significantly downregulated, while the expressions of Abd-B and Wnt1 were highly elevated in the OfUbx mutants (Figure S1B).
Ostrinia furnacalis is one of the most destructive pests of maize, especially in China and northeast Asia [73]. The extensive use of chemical pesticides has led to adverse environmental and non-target impacts, and the development of resistance has also become a serious concern for O. furnacalis management [73]. Therefore, novel pest management approaches, such as RNA-based control alternatives, are urgently needed [74,75].
Our study demonstrates that disruption of OfAbd-A can lead to embryonic lethality, and that OfUbx mutation induces adult sterility (Figure 8), which suggest their potential as targets of genetic tools for pest management. Since the transgenic line in O. furnacalis has been constructed for genetic control by using a piggybac transposon [76], the transgenic insect technique can be carried out by targeting Hox genes, such as OfAbd-A and OfUbx, in order to exert female-specific interruption (Figure S2). By releasing transgenic adult males that carry a female-specific promoter that initiates the Cas9 protein and the U6 promoter to drive the targeted gene-specific sgRNA expression, males that mate with females in the field can induce female-specific lethality and sterility in the next generation as a form of pest control (Figure S2). Future research regarding how to target OfAbd-A and OfUbx to induce the abnormal development of pupae and adults in the field is highly warranted. The two Hox genes that are identified and characterized in this study provide fundamental knowledge regarding the development of O. furnacalis and show potential as targets for genetic pest control that would be beneficial to maize growers, and that could also extend to other lepidopteran pests (Figure S2).

5. Conclusions

In this study, we hypothesized that OfAbd-A and OfUbx modulate morphological identity and development of thorax and abdomen in O. furnacalis. By using a newly developed CRISPR/Cas9 genome editing system in O. furnacalis, we knocked out OfAbd-A and OfUbx respectively, and subsequently observed substantial phenotypic impacts during insect development. Specifically, after OfAbd-A mutagenesis, abnormal development of abdominal segments was observed in larvae (Figure 6). The abnormal folding of the wings in pupae and adults was observed after OfUbx mutagenesis, which prevented the wings from completely covering the thoracic region of the pupae (Figure 7A), or limited the ability of adults to fly for the purposes of mating or reproducing (Figure 7B). Interestingly, the disruption of OfAbd-A led to embryonic lethality, while OfUbx mutation induced adult sterility. These combined results not only support our hypothesis, but they also provide a potential molecular target in homeotic genes for the genetic control of this global insect pest.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/insects13040384/s1, Figure S1: Relative transcript levels of downstream genes in OfAbd-A and OfUbx mutants. Figure S2: Prospective CRISPR/Cas9 mediated female specific transgenic pest control schematic. Table S1: Primers used in this study.

Author Contributions

Q.Z. designed the experiments; H.B. and X.L. performed the experiments; H.B. and J.G. analyzed the data; H.B. and Q.Z. drafted the manuscript; A.M. revised the manuscript; and X.Z. oversaw the entire project and revised and edited manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by grants from the Scientific and Technological Project of Henan Province (222102110108) to B.H.L., and the National Natural Science Foundation of China (31601687 and 32072482) to Z.Q.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available within the article and Supplementary Materials.

Acknowledgments

We are grateful to Yongping Huang from Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, for his support of experimental design and platform to initiate this project. In addition, the authors are grateful for the suggestions and comments by anonymous reviewers.

Conflicts of Interest

We declare that we have no competing interests.

References

  1. Nafus, D.M.; Schreiner, I. Review of the biology and control of the Asian corn borer, Ostrinia furnacalis (Lep: Pyralidae). Trop. Pest Manag. 1991, 37, 41–56. [Google Scholar] [CrossRef]
  2. Afidchao, M.M.; Musters, C.J.; de Snoo, G.R. Asian corn borer (ACB) and non-ACB pests in GM corn (Zea mays L.) in the Philippines. Pest Manag. Sci. 2013, 69, 792–801. [Google Scholar] [CrossRef] [PubMed]
  3. Wang, Z.Y.; Wu, Y.; He, K.L.; Bai, S.X. Effects of transgenic Bt maize pollen on longevity and fecundity of Trichogramma ostriniae in laboratory conditions. Bull. Insectol. 2007, 60, 49–55. [Google Scholar]
  4. Wang, C.L.; Wang, H.X.; Gui, C.M.; Lu, H. Studies on the control of the Asian corn borer, Ostrinia furnacalis (Guenee), with Trichogramma ostriniae. In US. National Acedemy of Science Joint Symposium on Biological Control of Insects; Adkisson, P.L., Ma, S.J., Eds.; Science Press: Beijing, China, 1984; pp. 268–273. [Google Scholar]
  5. Hassan, S.A.; Guo, M.F. Selection of effective strains of egg parasites of the genus Trichogramma (Hym., Trichogrammatidae) to control the European corn borer Ostrinia nubilalis Hb. (Lep., Pyralidae). J. Appl. Entomol. 1991, 111, 335–341. [Google Scholar] [CrossRef]
  6. Hoffmann, M.P.; Wright, M.G.; Pitcher, S.A.; Gardner, J. Inoculative releases of Trichogramma ostriniae for suppression of Ostrinia nublilalis (European corn borer) in sweet corn: Field biology and population dynamics. Biol. Control 2002, 25, 249–258. [Google Scholar] [CrossRef]
  7. Wright, M.G.; Kuhar, T.P.; Hoffmann, M.P.; Chenus, S.A. Effect of inoculative releases of Trichogramma Ostriniae on populations of Ostrinia nubilalis and damage to sweet corn and field corn. Biol. Control 2002, 23, 149–155. [Google Scholar] [CrossRef]
  8. Munkvold, G.P.; Hellmich, R.L.; Rice, L.G. Comparison of Fumonisin Concentrations in Kernels of Transgenic Bt Maize Hybrids and Nontransgenic Hybrids. Plant Dis. 1999, 83, 130–138. [Google Scholar] [CrossRef] [Green Version]
  9. Hurley, T.M.; Secchi, S.; Babcock, B.A.; Hellmichet, R.L. Managing the risk of European corn borer resistance to Bt corn. Environ. Resour. Econ. 2002, 22, 537–558. [Google Scholar] [CrossRef]
  10. Ferry, N.; Edwards, M.G.; Gatehouse, J.; Capell, T.; Christou, P.; Gatehouse, A.M. Transgenic plants for insect pest control: A forward looking scientific perspective. Transgenic Res. 2006, 15, 13–19. [Google Scholar] [CrossRef]
  11. Thompson, G.D.; Dalmacia, S.C.; Criador, A.R.; Alvarez, E.R.; Hechanova, R.F. Field Performance of TC1507 Transgenic Corn Hybrids Against Asian Corn Borer in the Philippines. Philipp. Agric. Sci. 2010, 93, 375–383. [Google Scholar]
  12. Xu, L.N.; Wang, Y.Q.; Wang, Z.Y.; Hu, B.J.; Ling, Y.H.; He, K.L. Transcriptome differences between Cry1Ab resistant and susceptible strains of Asian corn borer. BMC Genom. 2015, 16, 173. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Crespo, A.L.S.; Spencer, T.A.; Alves, A.P.; Hellmich, R.L.; Blankenship, E.E.; Magalhaes, L.C.; Siegfried, B.D. On-plant survival and inheritance of resistance to Cry1Ab toxin from Bacillus thuringiensis in a field-derived strain of European corn borer, Ostrinia nubilalis. Pest Manag. Sci. 2009, 65, 1071–1081. [Google Scholar] [CrossRef] [PubMed]
  14. Crespo, A.L.; Rodrigo-Simón, A.; Siqueira, H.A.; Pereira, E.J.; Ferre, J.; Siegfried, B.D. Cross-resistance and mechanism of resistance to Cry1Ab toxin from Bacillus thuringiensis in a field-derived strain of European corn borer, Ostrinia nubilalis. J. Invert. Pathol. 2011, 107, 185–192. [Google Scholar] [CrossRef]
  15. Smith, L.J.; Farhan, Y.; Schaafsma, W.A. Practical resistance of Ostrinia nubilalis (Lepidoptera: Crambidae) to Cry1F Bacillus thuringiensis maize discovered in Nova Scotia, Canada. Sci. Rep. 2019, 9, 18247. [Google Scholar] [CrossRef] [Green Version]
  16. Carroll, S.B. Endless forms: The evolution of gene regulation and morphological diversity. Cell 2000, 101, 577–580. [Google Scholar] [CrossRef] [Green Version]
  17. Lewis, E.B. A gene complex controlling segmentation in Drosophila. Nature 1978, 276, 565–570. [Google Scholar] [CrossRef]
  18. Krumlauf, R. Hox genes in vertebrate development. Cell 1994, 78, 191–201. [Google Scholar] [CrossRef]
  19. McGinnis, W.; Krumlauf, R. Homeobox genes and axial patterning. Cell 1992, 68, 283–302. [Google Scholar] [CrossRef]
  20. Reuter, R.; Scott, M.P. Expression and function of the homoeotic genes Antennapedia and Sex combs reduced in the embryonic midgut of Drosophila. Development 1990, 109, 289–303. [Google Scholar] [CrossRef]
  21. Maeda, R.K.; Karch, F. The bithorax complex of Drosophila an exceptional Hox cluster. Curr. Top. Dev. Biol. 2009, 88, 1–33. [Google Scholar]
  22. Liu, F.; Li, X.; Zhao, M.; Guo, M.; Han, K.; Dong, X.; Zhao, J.; Cai, W.; Zhang, Q.; Hua, H. Ultrabithorax is a key regulator for the dimorphism of wings, a main cause for the outbreak of planthoppers in rice. Natl. Sci. Rev. 2020, 7, 1181–1189. [Google Scholar] [CrossRef]
  23. Lawrence, P.A.; Morata, G. Homeobox genes: Their function in Drosophila segmentation and pattern formation. Cell 1994, 78, 181–189. [Google Scholar] [CrossRef]
  24. Ponzielli, R.; Astier, M.; Chartier, A.; Gallet, A.; Sémériva, M. Heart tube patterning in Drosophila requires integration of axial and segmental information provided by the Bithorax Complex genes and hedgehog signaling. Development 2002, 129, 4509–4521. [Google Scholar] [CrossRef]
  25. Vachon, G.; Cohen, B.; Pfeifle, C.; McGuffin, M.E.; Botas, J.; Cohen, S.M. Homeotic genes of the Bithorax complex repress limb development in the abdomen of the Drosophila embryo through the target gene Distal-less. Cell 1992, 71, 437–450. [Google Scholar] [CrossRef]
  26. Harding, K.; Wedeen, C.; McGinnis, W.; Levine, M. Spatially regulated expression of homeotic genes in Drosophila. Science 1985, 229, 1236–1242. [Google Scholar] [CrossRef]
  27. Cumberledge, S.U.; Szabad, J.A.; Sakonju, S.H. Gonad formation and development requires the abd-A domain of the bithorax complex in Drosophila melanogaster. Development 1992, 115, 395–402. [Google Scholar] [CrossRef]
  28. Foronda, D. Requirement of Abdominal-A and Abdominal-B in the developing genitalia of Drosophila breaks the posterior downregulation rule. Development 2006, 133, 117–127. [Google Scholar] [CrossRef] [Green Version]
  29. Marchetti, M.F. Differential expression of the Drosophila BX-C in polytene chromosomes in cells of larval fat bodies: A cytological approach to identifying in vivo targets of the homeotic Ubx, Abd-A and Abd-B proteins. Development 2003, 130, 3683–3689. [Google Scholar] [CrossRef] [Green Version]
  30. Mathies, L.D.; Kerridge, S.; Scott, M.P. Role of the teashirt gene in Drosophila midgut morphogenesis: Secreted proteins mediate the action of homeotic genes. Development 1994, 120, 2799–2809. [Google Scholar] [CrossRef]
  31. Perrin, L.; Monier, B.; Ponzielli, R.; Astier, M.; Semeriva, M. Drosophila cardiac tube organogenesis requires multiple phases of Hox activity. Dev. Biol. 2004, 272, 419–431. [Google Scholar] [CrossRef] [Green Version]
  32. Bi, H.L.; Xu, J.; Tan, A.J.; Huang, Y.P. CRISPR/Cas9-mediated targeted gene mutagenesis in Spodoptera litura. Insect Sci. 2016, 23, 469–477. [Google Scholar] [CrossRef]
  33. Roch, F.; Akam, M. Ultrabithorax and the control of cell morphology in Drosophila halteres. Development 2000, 127, 97–107. [Google Scholar] [CrossRef]
  34. Averof, M.; Patel, N.H. Crustacean appendage evolution associated with changes in Hox gene expression. Nature 1997, 388, 682–686. [Google Scholar] [CrossRef]
  35. Kaufman, T.C.; Lewis, R.; Wakimoto, B. Cytogenetic Analysis of Chromosome 3 in Drosophila melanogaster: The Homoeotic Gene Complex in Polytene Chromosome Interval 84a-B. Genetics 1980, 94, 115–133. [Google Scholar] [CrossRef]
  36. Garcia-Bellido, A. Genetic control of wing disc development in Drosophila. Ciba Found. Symp. 1975, 29, 161–182. [Google Scholar]
  37. Struhl, G. Genes controlling segmental specification in the Drosophila thorax. Proc. Natl. Acad. Sci. USA 1982, 79, 7380–7384. [Google Scholar] [CrossRef] [Green Version]
  38. Prasad, N.; Tarikere, S.; Khanale, D.; Habib, F.; Shashidhara, L.S. A comparative genomic analysis of targets of Hox protein Ultrabithorax amongst distant insect species. Sci. Rep. 2016, 6, 27885. [Google Scholar] [CrossRef] [Green Version]
  39. Gebelein, B.; Culi, J.; Ryoo, H.D.; Zhang, W.; Mann, R.S. Specificity of Distalless repression and limb primordia development by abdominal Hox proteins. Dev. Cell. 2002, 3, 487–498. [Google Scholar] [CrossRef] [Green Version]
  40. Ronshaugen, M.; McGinnis, N.; McGinnis, W. Hox protein mutation and macroevolution of the insect body plan. Nature 2002, 415, 914–917. [Google Scholar] [CrossRef]
  41. Pan, M.H.; Wang, X.Y.; Chai, C.L.; Zhang, C.D.; Lu, C.; Xiang, Z.H. Identification and function of Abdominal-A in the silkworm, Bombyx mori. Insect Mol. Biol. 2009, 18, 155–160. [Google Scholar] [CrossRef]
  42. Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G. MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Larkin, M.A.; Blackshields, G.; Brown, N.P.; Chenna, R.; McGettigan, P.A.; McWilliam, H.; Valentin, F.; Wallace, I.M.; Wilm, A.; Lopez, R.; et al. Clustal W and Clustal X version 2.0. Bioinformatics 2007, 23, 2947–2948. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Saitou, N.; Nei, M. The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 1987, 4, 406–425. [Google Scholar]
  45. Madeira, F.; Park, Y.M.; Lee, J.; Buso, N.; Gur, T.; Madhusoodanan, N.; Basutkar, P.; Tivey, A.R.; Potter, S.C.; Finn, R.D.; et al. The EMBL-EBI search and sequence analysis tools APIs in 2019. Nucleic Acids Res. 2019, 47, W636–W641. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Fang, G.Q.; Zhang, Q.; Chen, X.E.; Cao, Y.H.; Wang, Y.H.; Qi, M.M.; Wu, N.N.; Qian, L.S.; Zhu, C.X.; Huang, Y.P.; et al. The draft genome of the Asian corn borer yields insights into ecological adaptation of a devastating maize pest. Insect Biochem. Mol. Biol. 2021, 138, 103638. [Google Scholar] [CrossRef] [PubMed]
  47. Wang, X.; Xu, Y.; Huang, J.; Jin, W.; Yang, Y.; Wu, Y. CRISPR-Mediated Knockout of the ABCC2 Gene in Ostrinia furnacalis Confers High-Level Resistance to the Bacillus thuringiensis Cry1Fa Toxin. Toxins 2020, 12, 246. [Google Scholar] [CrossRef]
  48. You, L.; Bi, H.L.; Wang, Y.H.; Li, X.W.; Chen, X.E.; Li, Z.Q. CRISPR/Cas9-based mutation reveals Argonaute 1 is essential for pigmentation in Ostrinia furnacalis. Insect Sci. 2019, 26, 1020–1028. [Google Scholar] [CrossRef]
  49. Merrill, V.K.; Turner, F.R.; Kaufman, T.C. A genetic and developmental analysis of mutations in the Deformed locus in Drosophila melanogaster. Dev. Biol. 1987, 122, 379–395. [Google Scholar] [CrossRef]
  50. Balavoine, G.; de Rosa, R.; Adoutte, A. Hox clusters and bilaterian phylogeny. Mol. Phylogenet. Evol. 2002, 24, 366–373. [Google Scholar] [CrossRef]
  51. Duncan, I. The bithorax complex. Annu. Rev. Genet. 1987, 21, 285–319. [Google Scholar] [CrossRef]
  52. Angelini, D.R.; Liu, P.Z.; Hughes, C.L.; Kaufman, T.C. Hox gene function and interaction in the milkweed bug Oncopeltus fasciatus (Hemiptera). Dev. Biol. 2005, 287, 440–455. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Hughes, C.L.; Kaufman, T.C. Hox genes and the evolution of the arthropod body plan. Evol. Dev. 2002, 4, 459–499. [Google Scholar] [CrossRef] [PubMed]
  54. Pearson, J.C.; Lemons, D.; McGinnis, W. Modulating Hox gene functions during animal body patterning. Nat. Rev. Genet. 2005, 6, 893–904. [Google Scholar] [CrossRef] [PubMed]
  55. Konopova, B.; Akam, M. The Hox genes Ultrabithorax and abdominal-A specify three different types of abdominal appendage in the springtail Orchesella cincta (Collembola). EvoDevo 2014, 5, 2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Robertson, L.K.; Mahaffey, J.W. Insect homeotic complex genes and development, lessons from Drosophila and beyond. In Comprehensive Insect Science; Elsevier: London, UK, 2005; Volume 1, pp. 247–303. [Google Scholar]
  57. Chauvet, S.; Maurel-Zaffran, C.; Miassod, R.; Jullien, N.; Pradel, J.; Aragnol, D. dlarp, a new candidate Hox target in Drosophila whose orthologue in mouse is expressed at sites of epithelium/mesenchymal interactions. Dev. Dyn. 2000, 218, 401–413. [Google Scholar] [CrossRef]
  58. Khila, A.; Abouheif, E.; Rowe, L. Evolution of a novel appendage ground plan in water striders is driven by changes in the Hox gene Ultrabithorax. PLoS Genet. 2009, 5, e1000583. [Google Scholar] [CrossRef] [Green Version]
  59. Castelli-Gair, J.; Akam, M. How the Hox gene Ultrabithorax specifies two different segments: The significance of spatial and temporal regulation within metameres. Development 1995, 121, 2973–2982. [Google Scholar] [CrossRef]
  60. Mahfooz, N.S.; Li, H.; Popadić, A. Differential expression patterns of the hox gene are associated with differential growth of insect hind legs. Proc. Natl. Acad. Sci. USA. 2004, 101, 4877–4882. [Google Scholar] [CrossRef] [Green Version]
  61. Chen, W.; Yang, F.Y.; Xu, X.J.; Kumar, U.; He, W.Y.; You, M.S. Genetic control of Plutella xylostella in omics era. Arch. Insect Biochem. Physiol. 2019, 102, e21621. [Google Scholar] [CrossRef]
  62. Paulo, F.D.; Williamson, M.E.; Arp, A.P.; Li, F.; Sagel, A.; Skoda, S.R.; Sanchez-Gallego, J.; Vasquez, M.; Quintero, G.; Pérez de León, A.A.; et al. Specific Gene Disruption in the Major Livestock Pests Cochliomyia hominivorax and Lucilia cuprina Using CRISPR/Cas9. G3-Genes Genomes Genet. 2019, 9, 3045–3055. [Google Scholar] [CrossRef] [Green Version]
  63. Bi, H.L.; Xu, X.; Li, X.W.; Zhang, Y.; Huang, Y.P.; Li, K.; Xu, J. CRISPR Disruption of BmOvo Resulted in the Failure of Emergence and Affected the Wing and Gonad Development in the Silkworm Bombyx mori. Insects 2019, 10, 254. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Jin, W.Z.; Zhai, Y.Q.; Yang, Y.H.; Wu, Y.D.; Wang, X.L. Cadherin Protein Is Involved in the Action of Bacillus thuringiensis Cry1Ac Toxin in Ostrinia furnacalis. Toxins 2021, 13, 658. [Google Scholar] [CrossRef] [PubMed]
  65. Weatherbee, S.D.; Nijhout, H.F.; Grunert, L.W.; Halder, G.; Galant, R.; Selegue, J.; Carroll, S. Ultrabithorax function in butterfly wings and the evolution of insect wing patterns. Curr. Biol. 1999, 9, 109–115. [Google Scholar] [CrossRef] [Green Version]
  66. Deutsch, J. Hox and Wings. Bioessays 2005, 27, 673–675. [Google Scholar] [CrossRef]
  67. Pavlopoulos, A.; Akam, M. Hox gene Ultrabithorax regulates distinct sets of target genes at successive stages of Drosophila haltere morphogenesis. Proc. Natl. Acad. Sci. USA 2011, 108, 2855–2860. [Google Scholar] [CrossRef] [Green Version]
  68. Abzhanov, A.; Holtzman, S.; Kaufman, T.C. The Drosophila proboscis is specified by two Hox genes, proboscipedia and Sex combs reduced, via repression of leg and antennal appendage genes. Development 2001, 128, 2803–2814. [Google Scholar] [CrossRef]
  69. Martin, A.; Serano, J.M.; Jarvis, E.; Bruce, H.S.; Wang, J.; Ray, S.; Barker, C.A.; O’Connell, L.C.; Patel, N.H. CRISPR/Cas9 Mutagenesis Reveals Versatile Roles of Hox Genes in Crustacean Limb Specification and Evolution. Curr. Biol. 2015, 26, 14–26. [Google Scholar] [CrossRef] [Green Version]
  70. Martinez-Arias, A.; Ingham, P.W.; Scott, M.P.; Akam, M.E. The spatial and temporal deployment of Dfd and Scr transcripts throughout development of Drosophila. Development 1987, 100, 673–683. [Google Scholar] [CrossRef]
  71. Chen, P.; Tong, X.L.; Li, D.D.; Fu, M.Y.; He, S.Z.; Hu, H.; Xiang, Z.H.; Lu, C.; Dai, F.Y. Antennapedia is involved in the development of thoracic legs and segmentation in the silkworm, Bombyx mori. Heredity 2013, 111, 182–188. [Google Scholar] [CrossRef] [Green Version]
  72. Singh, N.P.; Mishra, R.K. Role of Abd-A and Abd-B in development of abdominal epithelia breaks posterior prevalence rule. PLoS Genet. 2014, 10, e1004717. [Google Scholar] [CrossRef] [Green Version]
  73. Liu, Q.; Hallerman, E.M.; Peng, Y.; Li, Y. Development of Bt Rice and Bt Maize in China and Their Efficacy in Target Pest Control. Int. J. Mol. Sci. 2016, 17, 1561. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Xu, L.; Zeng, B.; Noland Jeffery, E.; Zhou, X. The coming of RNA-based pest controls. J. Plant. Protection. 2015, 42, 673–690. [Google Scholar]
  75. Enserink, M.; Hines, J.H.; Vignieri, S.N.; Wigginton, N.S. Smarter pest control. The pesticide paradox. Introduction. Science 2013, 341, 728–729. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Liu, D.; Yan, S.C.; Huang, Y.P.; Tan, A.J.; Stanley, D.W.; Song, Q.S. Genetic transformation mediated by piggyBac in the Asian corn borer, Ostrinia furnacalis (Lepidoptera: Crambidae). Arch. Insect Biochem. Physiol. 2012, 80, 140–150. [Google Scholar] [CrossRef]
Figure 1. Phylogenetic analysis of OfAbd-A and OfUbx genes. The evolutionary histories of Abd-A and Ubx were inferred using the neighbor-joining method. The percentages of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) are shown next to the branches. Evolutionary distances were computed using the Poisson correction method and are listed as the number of amino acid substitutions per site. (A) The phylogenetic analysis of Abd-A involved 18 amino acid sequences. OfAbd-A was denoted by a square. (B) The phylogenetic analysis of Ubx involved 21 amino acid sequences. OfUbx was denoted by a square.
Figure 1. Phylogenetic analysis of OfAbd-A and OfUbx genes. The evolutionary histories of Abd-A and Ubx were inferred using the neighbor-joining method. The percentages of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) are shown next to the branches. Evolutionary distances were computed using the Poisson correction method and are listed as the number of amino acid substitutions per site. (A) The phylogenetic analysis of Abd-A involved 18 amino acid sequences. OfAbd-A was denoted by a square. (B) The phylogenetic analysis of Ubx involved 21 amino acid sequences. OfUbx was denoted by a square.
Insects 13 00384 g001
Figure 2. Sequence alignment of Abd-A and Ubx proteins. (A) Abd-A and (B) Ubx proteins were aligned from six insects, including previously reported sequences from D. melanogaster, A. mellifera, T. castaneum, P. xylostella, B. mori, and our newly identified sequences from O. furnacalis. Alignments of Abd-A and Ubx have common Homeodomain and specific conserved Abdominal-A and Ultrabithorax domains.
Figure 2. Sequence alignment of Abd-A and Ubx proteins. (A) Abd-A and (B) Ubx proteins were aligned from six insects, including previously reported sequences from D. melanogaster, A. mellifera, T. castaneum, P. xylostella, B. mori, and our newly identified sequences from O. furnacalis. Alignments of Abd-A and Ubx have common Homeodomain and specific conserved Abdominal-A and Ultrabithorax domains.
Insects 13 00384 g002
Figure 3. Spatial and temporal expression patterns of OfAbd-A and Ubx. The relative expressions of (A) OfAbd-A and (B) OfUbx in eggs, from the first day of first larval instar (L1D1) to the first day of fifth larval instar (L5D1), and at wandering stage, pupal stage, female adult stage (FA), and male adult stage (MA). The relative expressions of (C) OfAbd-A and (D) OfUbx in head, midgut (MG), foregut (FG), fat body (FB), epidermis (EPI), testes (TE), and ovaries (OV) at the third day of fifth instar (L5D3). Means labeled with different letters indicate significant difference at p < 0.05 (n = 3).
Figure 3. Spatial and temporal expression patterns of OfAbd-A and Ubx. The relative expressions of (A) OfAbd-A and (B) OfUbx in eggs, from the first day of first larval instar (L1D1) to the first day of fifth larval instar (L5D1), and at wandering stage, pupal stage, female adult stage (FA), and male adult stage (MA). The relative expressions of (C) OfAbd-A and (D) OfUbx in head, midgut (MG), foregut (FG), fat body (FB), epidermis (EPI), testes (TE), and ovaries (OV) at the third day of fifth instar (L5D3). Means labeled with different letters indicate significant difference at p < 0.05 (n = 3).
Insects 13 00384 g003
Figure 4. CRISPR/Cas9-mediated mutationswithin OfAbd-A target sites. (A) Target site of OfAbd-A genome locus focused on the third exon of the common region in the four splice variants. (B) Sequencing chromatogram of OfAbd-A mutants. The red wedge indicates position of cleavage by the CRISPR/Cas9 genome editing system. M1 (Mutant 1); M5 (Mutant 5). (C) Mutations detected by sequencing. The PAM sequence is in red. The black line represents the target site.
Figure 4. CRISPR/Cas9-mediated mutationswithin OfAbd-A target sites. (A) Target site of OfAbd-A genome locus focused on the third exon of the common region in the four splice variants. (B) Sequencing chromatogram of OfAbd-A mutants. The red wedge indicates position of cleavage by the CRISPR/Cas9 genome editing system. M1 (Mutant 1); M5 (Mutant 5). (C) Mutations detected by sequencing. The PAM sequence is in red. The black line represents the target site.
Insects 13 00384 g004
Figure 5. CRISPR/Cas9-mediated mutations within OfUbx target sites. (A) The target site of OfUbx genome locus in exon 1 of the common region for the two splice variants. (B) Sequencing chromatogram of OfUbx mutants. The red wedge indicates the position of cleavage by the CRISPR/Cas9 genome editing system. M1 (Mutant1); M2 (Mutant2). (C) Genotype detection of OfUbx genome sequence. The PAM sequence is in red. The black line represents the target site. The two target sites are separated by 151 bp.
Figure 5. CRISPR/Cas9-mediated mutations within OfUbx target sites. (A) The target site of OfUbx genome locus in exon 1 of the common region for the two splice variants. (B) Sequencing chromatogram of OfUbx mutants. The red wedge indicates the position of cleavage by the CRISPR/Cas9 genome editing system. M1 (Mutant1); M2 (Mutant2). (C) Genotype detection of OfUbx genome sequence. The PAM sequence is in red. The black line represents the target site. The two target sites are separated by 151 bp.
Insects 13 00384 g005
Figure 6. Segmentation malformation in OfAbd-A mutantsat larval and pupal stages. (A) The phenotypes of OfAbd-A mutants in the larval stage. Arrows show abnormal segments. Bar = 0.3 mm. (B) The phenotypes of OfAbd-A mutants in the pupal stage. Arrows show abnormal segments. Bar = 0.2 mm.
Figure 6. Segmentation malformation in OfAbd-A mutantsat larval and pupal stages. (A) The phenotypes of OfAbd-A mutants in the larval stage. Arrows show abnormal segments. Bar = 0.3 mm. (B) The phenotypes of OfAbd-A mutants in the pupal stage. Arrows show abnormal segments. Bar = 0.2 mm.
Insects 13 00384 g006
Figure 7. Wingdeformation in OfUbx mutantsat pupal and adult stages. (A) The phenotypes of OfUbx mutants in the larval stage. Arrows show abnormal and deficient wing discs. Bar = 0.3 mm. (B) The phenotypes of OfUbx mutants in the adult stage. Arrows show folded and deficient wings. Bar = 0.2 mm.
Figure 7. Wingdeformation in OfUbx mutantsat pupal and adult stages. (A) The phenotypes of OfUbx mutants in the larval stage. Arrows show abnormal and deficient wing discs. Bar = 0.3 mm. (B) The phenotypes of OfUbx mutants in the adult stage. Arrows show folded and deficient wings. Bar = 0.2 mm.
Insects 13 00384 g007
Figure 8. Phenotypic impacts of OfAbd-A and OfUbx mutagenesis on O. furnacalis hatching rate. Knocking out of OfAbd-A resulted in embryonic lethality, while OfUbx mutation led to adult sterility. From left to right, bars depict the hatching rates of eggs injected with different concentrations of OfAbd-A sgRNA, and the hatching rates of eggs produced by adult OfUbx mutants mated with each other. The asterisks (***) indicate significant differences (p < 0.01 or p < 0.001), compared to results from wild-type adults with a two-tailed t-test.
Figure 8. Phenotypic impacts of OfAbd-A and OfUbx mutagenesis on O. furnacalis hatching rate. Knocking out of OfAbd-A resulted in embryonic lethality, while OfUbx mutation led to adult sterility. From left to right, bars depict the hatching rates of eggs injected with different concentrations of OfAbd-A sgRNA, and the hatching rates of eggs produced by adult OfUbx mutants mated with each other. The asterisks (***) indicate significant differences (p < 0.01 or p < 0.001), compared to results from wild-type adults with a two-tailed t-test.
Insects 13 00384 g008
Table 1. Mutation frequency of CRISPR/Cas9-mediated mutagenesis in O. furnacalis.
Table 1. Mutation frequency of CRISPR/Cas9-mediated mutagenesis in O. furnacalis.
sgRNA Conc (ng/μL)Injected 1Hatched 2L Mutant 3Pupation 3P Mutant 4Adult 5A Mutant 6
OfAbd-A300785165 (21.0)92 (55.8)87 (52.7)29 (33.3)--
150596197 (33.1)75 (38.1)128 (65.1)18 (24.0)--
Ubx300486179 (36.8)-118 (65.9)53 (44.9)62 (52.5)18 (29.3)
150568234 (41.2)-160 (68.4)58 (36.2)95 (59.4)24 (25.3)
GFP300245152 (62.0)-108 (71.1)-70 (64.8)-
1 Number of injected individuals; 2 number and percent (%) of hatched individuals; 3 number and percent (%) of larvae that entered pupal stage; 4 number and percent (%) of pupal mutants; 5 number and percent (%) of pupae that entered adult stage; 6 number and percent (%) of adult mutants.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bi, H.; Merchant, A.; Gu, J.; Li, X.; Zhou, X.; Zhang, Q. CRISPR/Cas9-Mediated Mutagenesis of Abdominal-A and Ultrabithorax in the Asian Corn Borer, Ostrinia furnacalis. Insects 2022, 13, 384. https://0-doi-org.brum.beds.ac.uk/10.3390/insects13040384

AMA Style

Bi H, Merchant A, Gu J, Li X, Zhou X, Zhang Q. CRISPR/Cas9-Mediated Mutagenesis of Abdominal-A and Ultrabithorax in the Asian Corn Borer, Ostrinia furnacalis. Insects. 2022; 13(4):384. https://0-doi-org.brum.beds.ac.uk/10.3390/insects13040384

Chicago/Turabian Style

Bi, Honglun, Austin Merchant, Junwen Gu, Xiaowei Li, Xuguo Zhou, and Qi Zhang. 2022. "CRISPR/Cas9-Mediated Mutagenesis of Abdominal-A and Ultrabithorax in the Asian Corn Borer, Ostrinia furnacalis" Insects 13, no. 4: 384. https://0-doi-org.brum.beds.ac.uk/10.3390/insects13040384

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