Next Article in Journal
The Addition of Sterols and Cryoprotectants to Optimize a Diet Developed for Eldana saccharina Walker (Lepidoptera: Pyralidae) Using the Carcass Milling Technique
Next Article in Special Issue
Genome-Wide Screening of Transposable Elements in the Whitefly, Bemisia tabaci (Hemiptera: Aleyrodidae), Revealed Insertions with Potential Insecticide Resistance Implications
Previous Article in Journal
The Effect of Wind Speed on Male Potato Aphid, Macrosiphum euphorbiae, Responses to Primary Host Plant Volatiles and Female Sex Pheromone
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Miniature Inverted-Repeat Transposable Elements (MITEs) in the Two Lepidopteran Genomes of Helicoverpa armigera and Helicoverpa zea

1
Laboratory of Biochemistry and Biotechnology (LR01ES05), Faculty of Sciences of Tunis, University of Tunis El Manar, Tunis 1068, Tunisia
2
Laboratory of Biology of Organisms, Stress, Health, Environment (BiOSSE), Le Mans University, 72085 Le Mans, France
*
Authors to whom correspondence should be addressed.
Insects 2022, 13(4), 313; https://doi.org/10.3390/insects13040313
Submission received: 21 January 2022 / Revised: 10 March 2022 / Accepted: 20 March 2022 / Published: 23 March 2022
(This article belongs to the Special Issue Transposable Elements in Insects)

Abstract

:

Simple Summary

Miniature inverted-repeat transposable elements (MITEs) are non-autonomous transposable elements that play important roles in genome organization and evolution. Helicoverpa armigera and Helicoverpa zea shows a high number of reported cases of insecticide resistance worldwide, having evolved resistance against pyrethroids, organophosphates, carbamates, organochlorines, and recently to macrocyclic lactone spinosad and several Bacillus thuringiensis toxins. In the present study, we conducted a genome screening of MITEs in the H. armigera and H. zea genomes using bioinformatics approaches, and the results revealed a total of 3570 and 7405 MITE sequences in the H. armigera and H. zea genomes, respectively. Among these MITEs, we highlighted eleven MITE insertions in the H. armigera defensome genes and only one MITE insertion in those of H. zea.

Abstract

Miniature inverted-repeat transposable elements MITEs are ubiquitous, non-autonomous class II transposable elements. The moths, Helicoverpa armigera and Helicoverpa zea, are recognized as the two most serious pest species within the genus. Moreover, these pests have the ability to develop insecticide resistance. In the present study, we conducted a genome-wide analysis of MITEs present in H. armigera and H. zea genomes using the bioinformatics tool, MITE tracker. Overall, 3570 and 7405 MITE sequences were identified in H. armigera and H. zea genomes, respectively. Comparative analysis of identified MITE sequences in the two genomes led to the identification of 18 families, comprising 140 MITE members in H. armigera and 161 MITE members in H. zea. Based on target site duplication (TSD) sequences, the identified families were classified into three superfamilies (PIF/harbinger, Tc1/mariner and CACTA). Copy numbers varied from 6 to 469 for each MITE family. Finally, the analysis of MITE insertion sites in defensome genes showed intronic insertions of 11 MITEs in the cytochrome P450, ATP-binding cassette transporter (ABC) and esterase genes in H. armigera whereas for H. zea, only one MITE was retrieved in the ABC-C2 gene. These insertions could thus be involved in the insecticide resistance observed in these pests.

1. Introduction

Transposable elements (TEs) can play key innovating roles in their host genomes [1]. Indeed, their capacity to amplify in genomes and their ability to create new genetic variability by insertion/excision make them a rich source of genomic variants that can be selected through evolution [2,3]. Although most TE insertions are presumably deleterious or neutral, some insertions are expected to be beneficial to their carriers [4].
Miniature inverted-repeat transposable elements (MITEs) constitute a group of non-autonomous class II transposons, widespread and abundant in eukaryotic genomes [5]. MITEs are characterized by a small size (≤800 bp), a high copy number and a lack of any coding capacity [6]. They are further divided into sub-groups based on the similarity of their terminal inverted repeats (TIRs) and target site duplications (TSDs) to those of autonomous class II transposons [6]. MITEs were first described in plant genomes [7] and later found in a wide range of organisms, including insects [5,8]. It has been shown that MITEs are significant components in several insect genomes and that major MITE superfamilies are Tc1/mariner, PIF/harbinger, hAT, Mutator, and CACTA [9]. The small size of MITEs would allow them to escape the defense system of the host genome leading, therefore, to their accumulation [10]. Most MITE insertions, in multiple organisms, have been found in close proximity or even within genes [11]. Such a location entails the possibility of a functional impact on the nearby gene [12]. These results suggested that MITEs might have great effects on gene regulation and genome evolution. For example, the insertion of a MITE 200 bp upstream of the P450 gene CYP9M10 has been correlated to pyrethroid resistance in Culex quinquefasciatus [13].
The moths, Helicoverpa armigera and Helicoverpa zea, are both recognized as the most serious lepidopteran pests within the genus [14]. These two species diverged around 1.5 million years ago [15] and are considered sibling species due to their high morphological similarity, their genetic proximity, their emission of the same pheromone compounds but in different concentrations, and their capacity of interspecific crosses under natural and controlled environmental conditions [16].
H. armigera, the cotton bollworm, is native to the Old World (Asia, Europe, Africa, and Australasia) and is one of the most important pests worldwide. It is a polyphagous agricultural pest, and it was reported in more than 180 cultivated and wild plants, encompassing about 45 plant families [17].
H. zea, the corn earworm, has a wide distribution in the Americas. It is located from Canada to the south of Argentina. This species is polyphagous, and its larvae have been identified affecting leaves and fruits in more than 100 plant species [18].
H. armigera and H. zea have exhibited reduced susceptibility to groups of insecticides, including carbamates, organophosphates, pyrethroids, and Bacillus thuringiensis proteins [19]. Potential changes in the susceptibility of these species to conventional insecticides represent a major threat to agriculture in areas with established populations of these species and their potential hybrids. The H. armigera and H. zea genomes were sequenced in 2017 [15]. In the H. zea genome, there is a single strong signature of introgression around a region containing a novel, chimeric gene implicated in insecticide resistance (CYP337B3 gene), which was previously only documented in H. armigera and identified in the invasive population in Brazil [20]. Both H. armigera and H. zea genomes have been poorly investigated regarding their TE content. However, a previous study showed that 49% of all annotated TEs in the genome of H. armigera are MITEs [21]. Concerning H. zea, a genome walking study targeted on environment-adaptation genes showed MITEs present within or in close proximity to xenobiotic-metabolizing cytochrome P450 genes [22].
Here, we report the characterization of all the different groups of MITEs in the genomes of H. armigera and H. zea, using in silico approaches. The aim of this study was to provide an accurate description and comparison of MITEs in the two genomes, and to explore insecticide resistance genes in their proximity.

2. Materials and Methods

2.1. Identification of MITEs in H. zea

Both H. zea and H. armigera genomes were sequenced by Pearce et al. (2017) using the same technologies, from pupae of laboratory colonies derived from material collected 20 years ago [15]. The genomes were assembled by the same team using AllPaths version. JAN-2012 method. The quality of the two assemblies was tested using the QUAST tool [23] (Figures S1 and S2).
The H. zea genome available in GenBank-NCBI (BioProject PRJNA378438) is 341,147 Mb. This genome is assembled in 27,984 contigs and 2975 scaffolds corresponding to 20,602 kb and 201,477 kb N50 in length, respectively [15]. The H. zea assembled genomic sequences were downloaded and submitted to the MITE tracker tool with default parameters [6]. To find MITE candidates, MITE Tracker first searches for valid inverted repeat sequences of a given length (between 50 and 800 bp). At this step, a nucleotide−nucleotide BLAST search is used to align each MITE candidate to its reverse complement sequence (MITE maximum length was designed at 800 bp and TIRs at 10 bp). In the second step, putative MITEs were aligned and clustered into families by Vsearch [24] based on TSD and TIR sequences. Vsearch is executed with -id 0.8 (a similarity of 80% for clustering). The analysis of the obtained features allows a classification of MITE sequences into different superfamilies. For each element, right- and left-flanking sequences (sequences surrounding the element outside the TSD, by default 50 nucleotides length) are retrieved and compared with the flanking sequences of all other elements from the same family using a local alignment algorithm. The family is conserved only if the number of different individuals is equal or above a user-defined minimum copy number threshold (i.e., 3 by default). Finally, identified MITEs in H. zea and MITE sequences of H. armigera from previous study [21] were searched against GenBank, Repbase and i-MITE databases to find homologous MITEs from other insect species [9,25,26].
The evolutionary dynamics of MITEs was analyzed using the “TE” package implemented in R and based on statistical models [27]. This package estimates the age distribution based on mismatch values reflecting the deletion rate and random mutations in the sequences.

2.2. MITE Sequences Comparison in the H. armigera and H. zea Genomes

All annotated MITEs in H. armigera and H. zea were compared using BLAST searches with thresholds of 80% alignment length and 80% identity. The obtained MITE sequences were classified into families and superfamilies, using the 80:80:80 similarity rule and according their TSD sequences [28].
The identified MITE sequences in the two genomes were compared to each other and clustered based on a sequence identity higher than 80% and a minimum alignment coverage of 80% to their longest sequence. The obtained lineages were visualized by Cytoscape 3.7.2 software [29]. Additionally, the clustered MITE sequences were aligned based on a MUSCLE sequence alignment and a phylogenetic tree was generated using the maximum-likelihood method with 500 replications of bootstrap in MEGAX software [30].

2.3. Search for MITEs in Defensome Genes

The MITE sequences from the H. armigera and the H. zea genomes have been extended by 50 kb both upstream and downstream of their nucleotide sequences. The nucleotide BLAST was used to find defensome genes in the extended regions.

3. Results

3.1. MITEs Identification in the H. armigera and H. zea Genomes

In a previous study, we have reported on 3570 MITE sequences corresponding to seven known superfamilies and divided into 333 families in the H. armigera genome [21]. In the present study, MITE annotation in the genome of H. zea led to the identification of 7405 MITE sequences grouped in seven superfamilies and 435 families (Table 1).
Each MITE family was characterized by a short length, ranging from 50 to 800 bp and from 51 to 796 bp in H. armigera and H. zea, respectively. The TIR sequences were 10 to 32 bp long for H. armigera and 10 to 84 pb for H. zea. The Tc1/mariner was the main and the most diversified superfamily representing 1817 sequences belonging to 142 families in H. armigera. For H. zea, the main superfamily was PIF/harbinger (4587 sequences) and the most diversified family was TC1/mariner superfamily (188 families). CACTA, PiggyBac, hAT, Transib and Mavercik were the less represented superfamilies in the two analysed genomes.
Blast analysis of all MITEs identified in the two lepidopteran genomes against the Repbase and iMITE databases showed that H. armigera carries only one MITE which has a high similarity (83%) with a PIF/harbinger MITE from Bombyx mori. For H. zea, also only one MITE having an identity of 85% with a PIF/harbinger MITE from Spodoptera frugiperda was identified.
The MITE dynamic was evaluated, and the age distribution of these MITEs showed an expansion of MITE sequences at ∼90 million years ago (Mya) in H. armigera and a peak at 30 Mya in H. zea (Figure S3).

3.2. Characterization of MITE Families in H. armigera and H. zea Genomes

Comparative analysis of MITE sequences revealed a high identity between 140 MITEs (1501 copies) in H. armigera and 161 MITEs (1832 copies) in H. zea (Table 2, File S1). Based on their TIRs and TSD, we successfully identified 18 families belonging to three known superfamilies: Tc1/mariner, PIF/harbinger and CACTA. We named the PIF/harbinger families as Helicoverpa PIF/harbinger (HelPIF) 1 to 9, the Tc1/mariner families as Helicoverpa Tc1/mariner (HelTc1mar) one to six and the CACTA families as Helicoverpa CACTA (HelCac) one to three (Table 2). The clustering analysis of these MITEs using the cytoscape tool showed the separation of the 18 groups. The MITE sequences from H. armigera and H. zea belonging to the same family were grouped into the same cluster (Figure S4). These clusters were retrieved in the phylogenetic analysis with some difference for some families (Figure 1). For the PIF/harbinger superfamily, the results showed a clear separation of MITE members belonging to seven families, except for HelPIF-4 and HelPIF-7, showing a high sequence variation between MITE members. The MITEs of some TC1/mariner families appeared dispersed, in particular, the HelTc1mar-1 family split into five clades and the HelTc1mar-3 family clustered with the HelCac-3 family. MITEs belonging to the HelCac-1 and HelCac-2 are divided in four and two clades, respectively. However, all the 18 families showed a clear and a strong link between H. armigera and H. zea MITEs.
Eleven MITE families (HelPIF-1, -2, -5, -7, -9, HelTc1mar-1, -3, -4 and HelCac-1, -2, -3) were similar in copy numbers between both Helicoverpa genomes but only 7 MITE families showed 2- to 6-fold differences in total sequence number (underlined in red in Table 2). The HelPIF-2 family displays the highest copy numbers with 469 and 324 copies in H. armigera and H. zea, respectively, whereas the HelTc1mar-6 and HelTc1mar-4 families have the fewest copies corresponding to 6 and 10 copies in H. armigera and H. zea, respectively (Table 2).

3.3. MITE Insertions Analysis in Relation with Defensome Genes

BLAST searches led to the identification of 11 MITE insertions in defensome genes in the H. armigera genome. The TSD classification of the involved MITEs revealed that these sequences are members of Tc1/mariner and CACTA superfamilies (Table 3). Further analysis of MITE insertions has shown that six MITEs were inserted in cytochrome P450 genes, four elements were retrieved in ATP binding cassette transporter subfamily G gene and five were hosted by esterase genes. All inserted MITEs have intronic insertion sites.
For H. zea, we successfully identified 20 defensome genes by similarity with H. armigera (Table S1). We identified only one MITE sequence inserted in an ATP binding cassette transporter subfamily C2 (ABCC2) gene. This MITE was 137 bp long and belonged to the PIF/harbinger superfamily (Table 3).
Three MITE sequences were previously identified in the genome of H. zea and were inserted in the cytochrome P450 gene CYP9A14 (GenBank accession DQ788840.1) [22]. Blast searches with the CYP9A14 gene sequence identified the genomic scaffold KZ117493.1 (35,331 bp in length) in the public H. zea genome assembly (80% identity and 65% query cover). The gene position in this scaffold was between 13,600 bp and 16,000 bp, and only one MITE sequence, located from 17,076 bp to 17,206 bp, was identified upstream the gene. This MITE was 130 bp and belonged to the TC1/mariner superfamily (Table 3).

4. Discussion

In this study, we performed a systematic analysis to identify MITE families in the H. armigera and H. zea genomes published in year 2017 and publicly available. These genomes were both obtained from the pupae of laboratory colonies, by the same team, with the same technologies of sequencing and assembling, with a scaffold assembly level and comparable quality. The identification of MITE sequences was performed using a computational tool, MITE-tracker based on TIRs and TSDs features. In a previous study, a total of 3570 MITE sequences were identified in the genome of H. armigera [21]. In the present analysis, results showed a 2-fold difference in total MITE sequences in H. zea genome (7405 sequences) versus H. armigera genome. From all identified MITEs in H. armigera, only few MITEs were retrieved in the new genome assembly of H. armigera larvae available on NCBI from March 2021 (ASM1716586v1). This difference can be due to the molecular (difference in sequencing technology), biological (larvae vs. pupae) and bioinformatics difference between the two sequenced genomes.
The evolutionary dynamics and age distribution of identified MITEs in the H. armigera and H. zea genomes revealed a more recent and important burst of these sequences in H. zea genome at 30 Mya. However, it should be noted that the evolutionary dynamics and transposition timing may not be informative for very short and degenerated sequences such MITEs. This may partly explain the difference between the present expansion times and the divergence time of the two species, estimated to be 1.5 Mya according to the mitochondrial genome analysis performed by Behere et al. (2017) [31]. In addition, this difference may also be explained by a reactivation or invasion of active elements (by horizontal transfer) whose intact protein would interact with TIRs of MITEs and trigger their expansion. However, a more detailed analysis would be necessary to confirm this hypothesis.
The TSD classification of annotated MITEs showed that in H. armigera the main MITE superfamily was the Tc1/mariner superfamily, which is one of the predominant transposable element superfamilies in H. armigera genome [21] but in H. zea, the main MITE superfamily was the PIF/harbinger superfamily suggesting important genomic differences between the two sibling species.
In previous study, two MITE sequences (HzMITE2-1 and HzMITE3) were identified in a laboratory strain and a midgut cell line RP-HzGUT-AW1 of H. zea [22] which were retrieved in the present study and in the sequenced genome with less than 40% identity and cover.
The total MITE sequences identified in the two Helicoverpa genomes is much higher than the data from Han et al. (2016) in which the MITE hunter tool identified only 726 MITE sequences in nine lepidopteran species. This confirms the interest of the MITE tracker tool to increase the detection of MITEs in insect genomes, as it appears to be more efficient in detecting MITEs compared to other currently available tools.
Several studies tried to identify MITEs in insect genomes. Recently, 84 MITE families covering 1.3 Mb were identified in the genome of the diptera Mayetiola destructor [32]. In the hemidipteran Bemisia tabaci genome, 71 MITE sequences were identified as mariner-like element derivatives [33]. Eight novel families of MITEs were discovered in the African malaria mosquito, Anopheles gambiae [5]. In Bombyx mori, 17 MITE families with a total of 5785 members were identified [8]. However, MITE sequences available in databases are very restricted. The database of transposable elements, Repbase, contains only seven insect MITE sequences from Aedes aegypti [25] and the iMITE database include 6012 MITE sequences from 98 insect genomes [9]. Blast analysis of MITEs against Repbase and iMITE databases have shown a high similarity for only two MITEs of H. armigera and H. zea with two known MITEs in Bombyx mori and Spodoptera frugiperda, suggesting that the two lepidopteran genomes contain a high number of new or unknown MITE sequences. This may also be due to the lack of significant sequence conservation between MITEs of divergent species.
The identified MITE sequences were divided into 333 and 435 families in the H. armigera and H. zea genome, respectively, suggesting a tremendous diversity in the H. armigera and H. zea MITEs. The TIRs and TSDs similarity of MITE sequences inside different families suggests that they may arise from the amplification of a few progenitor copies. After blast searches and manual editing, a total of 18 MITE families were identified in H. armigera and H. zea genomes. Cytoscape and phylogenetic analysis showed that the MITE sequences in both genomes were closely related and showed a high identity level between MITEs of the two genomes. MITE members from certain family showed a dispersion into different clades, particularly the HelTc1mar-1 family. This dispersion may be due to the small size of the MITEs and great care must be taken in interpreting the phylogenetic analysis of such sequences without their original TEs.
The numbers of MITEs in 12 families were similar, suggesting that the major sequences of these families evolved before the divergence of the two species around 1.5 Mya. However, seven MITE families (HelPIF-3, 4, 6, 8, HelTc1mar-2, 5, 6) displayed large variation in copy numbers between the two species. The MITE families examined herein were classified into three superfamilies based on their TSDs, which were the main MITE superfamilies in the H. armigera and H. zea genomes.
After searching for defensome genes, eleven MITEs were found inserted in cytochrome P450, ABC transporter and esterase defensome genes of H. armigera. All inserted MITEs are located in intron. They could result in exon skipping, alternative splicing, or even in alterations in expression profiles if the corresponding introns contain regulatory sequences as exemplified by the Mu insertion into an intron of the knotted locus in maize [34].
In P450 genes, six MITEs were identified in H. armigera. The overexpression of P450 genes has been reported to increase the ability to metabolize insecticides in insect pests with agricultural importance and public health [35]. In Aedes aegypti, the P450 4C1 was highly upregulated in a resistant strain [36]. Three MITEs were inserted in carboxylesterase genes in H. armigera. Temephos resistance has been associated with the up-regulation, through gene amplification, of two carboxylesterase (CCE) genes closely linked on the genome of Aedes albopictus [37]. In the ABCG member-20 gene of H. armigera, a single MITE was identified in the position 89,853–90,098 bp. In a previous study, a hAT transposon in the position 1486–2436 bp and a Tc1/mariner transposon from 91,998 bp to 93,071 bp were inserted in this ABCG20 gene [21] indicating that several types of TE can be inserted into this ABCG20 gene. It is likely that each type of TEs may induce different genomic changes upon transposition. Consequently, this ABC gene should be gnomically highly variable.
For H. zea, only one MITE belonging to the PIF/harbinger superfamily was inserted in an ABCC2 gene. The function of this ABCC2 gene was recently evaluated in H. zea, using CRISPR/Cas9 and results observed support that this gene is not a major Cry1Ac receptor in this insect [38].
Thus, the present study will deepen our knowledge of MITEs and their position near defensome genes, which are fundamental to understanding the evolution and adaptation of these pest genomes and insecticides resistance of these species.

5. Conclusions

Genomic analysis of H. zea identified 7405 MITE sequences belonging to the Tc1/mariner, PIF harbinger, CACTA, PiggyBac, hAT, Transib and Maverick superfamilies. The genome-wide comparative study of around 300 MITE members belonging to 18 families in H. armigera and H. zea genomes highlighted important genomic similarities between the two sibling species. At all, 11 MITEs were identified in defensome genes in H. armigera and only one MITE in the genome of H. zea. This disparity is mainly due to the poor annotation of H. zea genome, for which, in particular, a very small number of defensome genes have been identified. The identification and characterization of MITEs in the two lepidopteran genomes provide a basis for further studies on their genomic impact. It will also facilitate the use of MITE insertions to identify genetic factors associated with defensome genes of H. armigera and H zea.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/insects13040313/s1. File S1. MITE sequences of 18 families in H. armigera and H. zea genomes. Table S1. Defensome gene identification in H. zea genome by similarity with H. armigera. Table S2. Mismatch calculation of identified MITEs in H. armigera. Table S3. Mismatch calculation of identified MITEs in H. zea, Figure S1. Cumulative length in H.armigera and H. zea assemblies. Figure S2. GC content in the H. armigera and H. zea genomes. Figure S3. Age distribution of MITEs in H. armigera (A) and H. zea (B). Figure S4. Cytoscape visualisation of MITE families belonging to PIF-harbinger (A), Tc1/mariner (B) and CACTA (C) superfamily in H. armigera and H. zea genomes.

Author Contributions

Conceptualization, K.K., M.M.K., A.C., B.C. and N.C.; methodology, K.K., M.Z.; investigation, data curation, K.K.; writing—original draft preparation, K.K., M.M.K., B.C. and N.C.; writing—review and editing, M.M.K., B.C., N.C., A.C. and F.D.; supervision, M.M.K., B.C. and N.C.; project administration, M.M.K. and N.C.; funding acquisition, M.M.K. and N.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by PHC-Utique 19G0901 “Hubert Curien Franco-Tunisian partnership” (41703ZJ) and the Tunisian Ministry of Higher Education and Scientific Research and the University of Tunis El Manar.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank Johann CONFAIS for his useful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Roquis, D.; Robertson, M.; Yu, L.; Thieme, M.; Julkowska, M.; Bucher, E. Genomic impact of stress-induced transposable element mobility in Arabidopsis. Nucleic Acids Res. 2021, 49, 10431–10447. [Google Scholar] [CrossRef] [PubMed]
  2. Chénais, B.; Caruso, A.; Hiard, S.; Casse, N. The impact of transposable elements on eukaryotic genomes: From genome size increase to genetic adaptation to stressful environments. Gene 2012, 509, 7–15. [Google Scholar] [CrossRef]
  3. Castanera, R.; Vendrell-Mir, P.; Bardil, A.; Carpentier, M.C.; Panaud, O.; Casacuberta, J.M. Amplification dynamics of miniature inverted-repeat transposable elements and their impact on rice trait variability. Plant J. 2021, 107, 118–135. [Google Scholar] [CrossRef] [PubMed]
  4. Gilbert, C.; Peccoud, J.; Cordaux, R. Transposable elements and the evolution of insects. Annu. Rev. Entomol. 2021, 66, 355–372. [Google Scholar] [CrossRef] [PubMed]
  5. Tu, Z. Eight novel families of miniature inverted repeat transposable elements in the African malaria mosquito, Anopheles gambiae. Proc. Natl. Acad. Sci. USA 2001, 98, 1699–1704. [Google Scholar] [CrossRef]
  6. Crescente, J.M.; Zavallo, D.; Helguera, M.; Vanzetti, L.S. MITE Tracker: An accurate approach to identify miniature inverted-repeat transposable elements in large genomes. BMC Bioinform. 2018, 19, 348. [Google Scholar] [CrossRef]
  7. Bennetzen, J.L. Transposable element contributions to plant gene and genome evolution. Plant Mol. Biol. 2000, 42, 251–269. [Google Scholar] [CrossRef]
  8. Han, Y.; Wessler, S.R. MITE-Hunter: A program for discovering miniature inverted-repeat transposable elements from genomic sequences. Nucleic Acids Res. 2010, 38, e199. [Google Scholar] [CrossRef] [Green Version]
  9. Han, M.J.; Zhou, Q.Z.; Zhang, H.H.; Tong, X.; Lu, C.; Zhang, Z.; Dai, F. iMITEdb: The genome-wide landscape of miniature inverted-repeat transposable elements in insects. Database 2016, 2016, baw148. [Google Scholar] [CrossRef] [Green Version]
  10. Sabot, F.; Schulman, A.H. Parasitism and the retrotransposon life cycle in plants: A hitchhiker’s guide to the genome. Heredity 2006, 97, 381–388. [Google Scholar] [CrossRef] [Green Version]
  11. Fattash, I.; Rooke, R.; Wong, A.; Hui, C.; Luu, T.; Bhardwaj, P.; Yang, G. Miniature inverted-repeat transposable elements: Discovery, distribution, and activity. Genome 2013, 56, 475–486. [Google Scholar] [CrossRef] [PubMed]
  12. Drongitis, D.; Aniello, F.; Fucci, L.; Donizetti, A. Roles of transposable elements in the different layers of gene expression regulation. Int. J. Mol. Sci. 2019, 20, 5755. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Itokawa, K.; Komagata, O.; Kasai, S.; Okamura, Y.; Masada, M.; Tomita, T. Genomic structures of Cyp9m10 in pyrethroid resistant and susceptible strains of Culex quinquefasciatus. Insect Biochem. Mol. Biol. 2010, 40, 631–640. [Google Scholar] [CrossRef] [PubMed]
  14. Amano, T.; Nomura, M. Re-examination of morphological variations in the female internal genitalia of Helicoverpa armigera and Helicoverpa zea (Lepidoptera: Noctuidae) for identification and pest management. Fla. Entomol. 2021, 104, 218–221. [Google Scholar] [CrossRef]
  15. Pearce, S.L.; Clarke, D.F.; East, P.D.; Elfekih, S.; Gordon, K.H.J.; Jermiin, L.S.; McGaughran, A.; Oakeshott, J.G.; Papanikolaou, A.; Perera, O.P.; et al. Genomic innovations, transcriptional plasticity and gene loss underlying the evolution and divergence of two highly polyphagous and invasive Helicoverpa pest species. BMC Biol. 2017, 15, 63. [Google Scholar] [CrossRef] [Green Version]
  16. Cordeiro, E.M.; Pantoja-Gomez, L.M.; de Paiva, J.B.; Nascimento, A.R.; Omoto, C.; Michel, A.P.; Correa, A.S. Hybridization and introgression between Helicoverpa armigera and H. zea: An adaptational bridge. BMC Evol. Biol. 2020, 20, 61. [Google Scholar] [CrossRef]
  17. Murùa, M.G.; Fogliata, S.V.; Herrero, M.I.; Vera, M.A.; Casmuz, A.S.; Sosa-Gomez, D.R. Biological and reproductive parameters of Helicoverpa armigera and Helicoverpa zea reared on artificial diet in Argentina. Bull. Insectol. 2021, 74, 55–64. [Google Scholar]
  18. Capinera, J.L. Corn earworm, Helicoverpa zea (Boddie)(Lepidoptera: Noctuidae). EDIS 2002, 7, 145. [Google Scholar] [CrossRef]
  19. da Silva, F.R.; Trujillo, D.; Bernardi, O.; Verle Rodrigues, J.C.; Bailey, W.D.; Gilligan, T.M.; Carrillo, D. Comparative toxicity of Helicoverpa armigera and Helicoverpa zea (Lepidoptera: Noctuidae) to selected insecticides. Insects 2020, 11, 431. [Google Scholar] [CrossRef]
  20. Valencia-Montoya, W.A.; Elfekih, S.; North, H.L.; Meier, J.I.; Warren, I.A.; Tay, W.T.; Gordon, K.H.J.; Specht, A.; Paula-Moraes, S.V.; Rane, R.; et al. Adaptive introgression across semipermeable species boundaries between local Helicoverpa zea and invasive Helicoverpa armigera moths. Mol. Biol. Evol. 2020, 37, 2568–2583. [Google Scholar] [CrossRef]
  21. Klai, K.; Chenais, B.; Zidi, M.; Djebbi, S.; Caruso, A.; Denis, F.; Confais, J.; Badawi, M.; Casse, N.; Mezghani Khemakhem, M. Screening of Helicoverpa armigera mobilome revealed transposable element insertions in insecticide resistance genes. Insects 2020, 11, 879. [Google Scholar] [CrossRef] [PubMed]
  22. Chen, S.; Li, X. Transposable elements are enriched within or in close proximity to xenobiotic-metabolizing cytochrome P450 genes. BMC Evol. Biol. 2007, 7, 46. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Gurevich, A.; Saveliev, V.; Vyahhi, N.; Tesler, G. QUAST: Quality assessment tool for genome assemblies. Bioinformatics 2013, 29, 1072–1075. [Google Scholar] [CrossRef] [PubMed]
  24. Rognes, T.; Flouri, T.; Nichols, B.; Quince, C.; Mahé, F. VSEARCH: A versatile open-source tool for metagenomics. PeerJ 2016, 4, e2584. [Google Scholar] [CrossRef] [PubMed]
  25. Benson, D.A.; Cavanaugh, M.; Clark, K.; Karsch-Mizrachi, I.; Lipman, D.J.; Ostell, J.; Sayers, E.W. GenBank. Nucleic Acids Res. 2012, 41, D36–D42. [Google Scholar] [CrossRef] [Green Version]
  26. Kapitonov, V.V.; Jurka, J. A universal classification of eukaryotic transposable elements implemented in Repbase. Nat. Rev. Genet. 2008, 9, 411–412. [Google Scholar] [CrossRef]
  27. Dai, X.; Wang, H.; Dvorak, J.; Bennetzen, J.; Mueller, H.G.; Dai, M.X. Package ‘TE’. 2018. Available online: https://cran.r-project.org/ (accessed on 19 March 2022).
  28. Wicker, T.; Sabot, F.; Hua-Van, A.; Bennetzen, J.L.; Capy, P.; Chalhoub, B.; Flavell, A.; Leroy, P.; Morgante, M.; Panaud, O.; et al. A unified classification system for eukaryotic transposable elements. Nat. Rev. Genet. 2007, 8, 973–982. [Google Scholar] [CrossRef]
  29. Demchak, B.; Hull, T.; Reich, M.; Liefeld, T.; Smoot, M.; Ideker, T.; Mesirov, J.P. Cytoscape: The network visualization tool for Genome Space. F1000Research 2014, 3, 151. [Google Scholar] [CrossRef]
  30. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547. [Google Scholar] [CrossRef]
  31. Behere, G.T.; Tay, W.T.; Russell, D.A.; Heckel, D.G.; Appleton, B.R.; Kranthi, K.R.; Batterham, P. Mitochondrial DNA analysis of field populations of Helicoverpa armigera (Lepidoptera: Noctuidae) and of its relationship to H. zea. BMC Evol. Biol. 2007, 7, 117. [Google Scholar] [CrossRef] [Green Version]
  32. Ben Amara, W.; Quesneville, H.; Khemakhem, M.M. A Genomic Survey of Mayetiola destructor Mobilome Provides New Insights into the Evolutionary History of Transposable Elements in the Cecidomyiid Midges. PLoS ONE 2021, 16, e0257996. [Google Scholar] [CrossRef] [PubMed]
  33. Zidi, M.; Denis, F.; Klai, K.; Chénais, B.; Caruso, A.; Djebbi, S.; Mezghani, M.; Casse, N. Genome-wide characterization of Mariner-like transposons and their derived MITEs in the Whitefly Bemisia tabaci (Hemiptera: Aleyrodidae). G3 2021, 11, jkab287. [Google Scholar] [CrossRef] [PubMed]
  34. Greene, B.; Walko, R.; Hake, S. Mutator insertions in an intron of the maize knotted1 gene result in dominant suppressible mutations. Genetics 1994, 138, 1275–1285. [Google Scholar] [CrossRef] [PubMed]
  35. Wang, Y.; Tian, J.; Han, Q.; Zhang, Y.; Liu, Z. Genomic organization and expression pattern of cytochrome P450 genes in the wolf spider Pardosa pseudoannulata. Comp. Biochem. Physiol. Part C Toxicol. Pharmacol. 2021, 248, 109118. [Google Scholar] [CrossRef] [PubMed]
  36. Lien, N.T.K.; Ngoc, N.T.H.; Lan, N.N.; Hien, N.T.; Van Tung, N.; Ngan, N.T.T.; Hoang, N.H.; Binh, N.T.H. Transcriptome sequencing and analysis of changes associated with insecticide resistance in the dengue mosquito (Aedes aegypti) in Vietnam. Am. J. Trop. Med. Hyg. 2019, 100, 1240. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Grigoraki, L.; Pipini, D.; Labbe, P.; Chaskopoulou, A.; Weill, M.; Vontas, J. Carboxylesterase gene amplifications associated with insecticide resistance in Aedes albopictus: Geographical distribution and evolutionary origin. PLoS Negl. Trop. Dis. 2017, 11, e0005533. [Google Scholar] [CrossRef] [Green Version]
  38. Perera, O.P.; Little, N.S.; Abdelgaffar, H.; Jurat-Fuentes, J.L.; Reddy, G.V. Genetic knockouts indicate that the ABCC2 protein in the bollworm Helicoverpa zea is not a major receptor for the Cry1Ac insecticidal protein. Genes 2021, 12, 1522. [Google Scholar] [CrossRef]
Figure 1. Phylogenetic analysis of MITE families from H. armigera and H. zea. The analysis was performed using the identified MITE family members having 80:80 coverage. The 18 MITE families were separated into different clades by a maximum likelihood phylogenetic analysis (model HKY85) with 500 bootstrap replications (only bootstrap value ≥ 0.5 are showed). The families belonging to the three MITE superfamilies, TC1/mariner, PIF harbinger and CACTA, are indicated in different colors and MITE members from H. armigera H. zea are indicated in black and red, respectively.
Figure 1. Phylogenetic analysis of MITE families from H. armigera and H. zea. The analysis was performed using the identified MITE family members having 80:80 coverage. The 18 MITE families were separated into different clades by a maximum likelihood phylogenetic analysis (model HKY85) with 500 bootstrap replications (only bootstrap value ≥ 0.5 are showed). The families belonging to the three MITE superfamilies, TC1/mariner, PIF harbinger and CACTA, are indicated in different colors and MITE members from H. armigera H. zea are indicated in black and red, respectively.
Insects 13 00313 g001
Table 1. MITE superfamilies detected in H. armigera and H. zea genomes.
Table 1. MITE superfamilies detected in H. armigera and H. zea genomes.
H. armigera *H. zea
SuperfamiliesTSDMITEs Length (pb)TIR Length (pb)MITE SequencesFamiliesMITEs Length (pb)TIR Length (pb)MITE SequencesFamilies
Tc1/marinerTA50–36010–21181714285–79410–712328188
PIF/HarbingerTWA, AT or AWT55–68515–32136811160–72013–844587154
CACTA2–4 bp78–77510–262507051–79617–4945086
PiggyBacTTAA50–80015–31932256–77521–69254
hAT8 bp56–26017–29203120–55028–63151
Transib5 pb83–38613–2716275–56436–5751
Maverick6 pb50–80010–2463306–65418–4251
Total 3570333 7415435
(*) Values from our previous study [21].
Table 2. Characteristics of the 18 MITE families in the H. armigera and H. zea genomes.
Table 2. Characteristics of the 18 MITE families in the H. armigera and H. zea genomes.
H. armigeraH. zea
FamiliesMITE Length (bp)TIR Length (bp)Sequences NumberCopies NumberMITE SequencesCopies Number
1HelPIF-1155–23425–652521013174
2HelPIF-2121–16022–381746911324
3 HelPIF-3 213–22958–72916121323
4 HelPIF-4 374–43930–432271153
5HelPIF-5111–14035–44573581
6 HelPIF-6 198–25653–69215531
7HelPIF-7163–17439–47447239
8 HelPIF-8 367–38526–3617521
9HelPIF-9290–33922–37221215
Total PIF/Harbinger1114392272671030751061
10HelTc1mar-1129–16421–572311427173
11 HelTc1mar-2 98–12528–381410320205
12HelTc1mar-3255–29516–24555276
13HelTc1mar-4195–21221–34318210
14 HelTc1mar-5 260–28337–5918448
15 HelTc1mar-6 457–67910–2916336
Total Tc1/mariner9867910594730458548
16HelCac-1348–46916–2513761381
17HelCac-2352–48616–46740867
18HelCac-3260–27415–20651775
Total CACTA26048615462616728223
Total MITEs 14015011611832
Table 3. MITE insertions in defensome genes of H. armigera and H. zea genomes.
Table 3. MITE insertions in defensome genes of H. armigera and H. zea genomes.
Gene FamiliesGene NameGene LengthMITE Inserted NameMITE Length (bp)TIR Length (bp)Insertion PositionExon/Intron Position
Helicoverpa armigeraCytochrome P450CYP450 4V2-like (LOC110375933)14,709MITE_armg_527 390189793–10,182Intron 11
CYP450 4V2-like (LOC110375407)8843MITE_armg_2987 189578301–8489Intron 11
CYP450 4C1-like (LOC110375947)10,768MITE_armg_4229 202563478–3651Intron 4
CYP450 6B5-like (LOC110371743)15,792MITE_armg_5115 129273727–3855Intron 1
Probable CYP450 49a1 (LOC110372238)38,869MITE_armg_6331 107178631–8737Intron 2
Probable CYP450 6d2 (LOC110383081)4170MITE_armg_6757 96261309–1404Intron 1
ABCABCG 49-like (LOC110377844)48,302MITE_armg_7483 4862125,885–26,370Intron 1
ABCG 49-like (LOC110374586)55,698MITE_armg_5600 1362130,248–30,383Intron 23
ABCG 20 (LOC110376033)96,146MITE_armg_1805 2461089,853–90,098Intron 12
ABCG 23 (LOC110373590)60,734MITE_armg_575812432855–978Intron 1
EsteraseEsterase FE4 like (LOC110380254)9333MITE_armg_1387 287133323–3609Intron 4
Esterase FE4 like (LOC110384365)8458MITE_armg_6918 641172829–3469Intron 1
Carboxylesterase 1C (LOC110375169)6830MITE_armg_54881401446127–6266Intron 10
Carboxylesterase 1E (LOC110379202)27,073MITE_armg_1116 3345018,270–18,603Intron 3
Venom carboxylesterase-6 (LOC110373494)27,998MITE_armg_335 393409211–9603Intron 1
Helicoverpa zeaABC transporterABC C2 (KY701524.1)11560MITE_zea_8125137152663–2799Intron 5
Cytochrome P450CYP9A14 (KZ117493.1)2400MITE_zea_510013010-Upstream
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Klai, K.; Zidi, M.; Chénais, B.; Denis, F.; Caruso, A.; Casse, N.; Mezghani Khemakhem, M. Miniature Inverted-Repeat Transposable Elements (MITEs) in the Two Lepidopteran Genomes of Helicoverpa armigera and Helicoverpa zea. Insects 2022, 13, 313. https://doi.org/10.3390/insects13040313

AMA Style

Klai K, Zidi M, Chénais B, Denis F, Caruso A, Casse N, Mezghani Khemakhem M. Miniature Inverted-Repeat Transposable Elements (MITEs) in the Two Lepidopteran Genomes of Helicoverpa armigera and Helicoverpa zea. Insects. 2022; 13(4):313. https://doi.org/10.3390/insects13040313

Chicago/Turabian Style

Klai, Khouloud, Marwa Zidi, Benoît Chénais, Françoise Denis, Aurore Caruso, Nathalie Casse, and Maha Mezghani Khemakhem. 2022. "Miniature Inverted-Repeat Transposable Elements (MITEs) in the Two Lepidopteran Genomes of Helicoverpa armigera and Helicoverpa zea" Insects 13, no. 4: 313. https://doi.org/10.3390/insects13040313

APA Style

Klai, K., Zidi, M., Chénais, B., Denis, F., Caruso, A., Casse, N., & Mezghani Khemakhem, M. (2022). Miniature Inverted-Repeat Transposable Elements (MITEs) in the Two Lepidopteran Genomes of Helicoverpa armigera and Helicoverpa zea. Insects, 13(4), 313. https://doi.org/10.3390/insects13040313

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