Next Article in Journal
Seasonal Flight Pattern of the Kissing Bugs Triatoma rubida and T. protracta (Hemiptera: Reduviidae: Triatominae) in Southern Arizona, United States
Previous Article in Journal
Modeling of Control Efforts against Rhipicephalus sanguineus, the Vector of Rocky Mountain Spotted Fever in Sonora Mexico
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Genomics and Population Genomics of the Light Brown Apple Moth, Epiphyas postvittana, an Invasive Tortricid Pest of Horticulture

1
The New Zealand Institute of Plant and Food Research Limited, Auckland 1025, New Zealand
2
The New Zealand Institute of Plant and Food Research Limited, Palmerston North 4410, New Zealand
3
School of Biological Sciences, University of Auckland, Auckland 1010, New Zealand
4
Dipartimento di Biologia, Università degli Studi di Padova, 35131 Padova, Italy
*
Author to whom correspondence should be addressed.
Insects 2022, 13(3), 264; https://doi.org/10.3390/insects13030264
Submission received: 17 January 2022 / Revised: 1 March 2022 / Accepted: 4 March 2022 / Published: 7 March 2022
(This article belongs to the Section Insect Molecular Biology and Genomics)

Abstract

:

Simple Summary

In this study, we produced a genomic resource for the light brown apple moth, Epiphyas postvittana, to understand the biological basis of adaptation to a high number of hosts (polyphagy) and the invasive nature of this and other lepidopteran pests. The light brown apple moth is an invasive pest of horticultural plants, with over 500 recorded plant hosts. With origins in Australia, the pest has subsequently spread to New Zealand, Hawaii, California and Europe, causing significant economic losses for fruit producers. Comparative genomic analyses with other lepidopteran genomes indicate that a high proportion of the genome is made up of repetitive sequences, with the majority of the known elements being DNA transposable elements and retrotransposons. Twenty gene families show significant expansions, including some likely to have a role in its pest status. Finally, population genomics, investigated by a RAD-tag approach, indicated likely patterns of invasion and admixture, with Californian moths most probably being derived from Australia.

Abstract

The light brown apple moth, Epiphyas postvittana is an invasive, polyphagous pest of horticultural systems around the world. With origins in Australia, the pest has subsequently spread to New Zealand, Hawaii, California and Europe, where it has been found on over 500 plants, including many horticultural crops. We have produced a genomic resource, to understand the biological basis of the polyphagous and invasive nature of this and other lepidopteran pests. The assembled genome sequence encompassed 598 Mb and has an N50 of 301.17 kb, with a BUSCO completion rate of 97.9%. Epiphyas postvittana has 34% of its assembled genome represented as repetitive sequences, with the majority of the known elements made up of longer DNA transposable elements (14.07 Mb) and retrotransposons (LINE 17.83 Mb). Of the 31,389 predicted genes, 28,714 (91.5%) were assigned to 11,438 orthogroups across the Lepidoptera, of which 945 were specific to E. postvittana. Twenty gene families showed significant expansions in E. postvittana, including some likely to have a role in its pest status, such as cytochrome p450s, glutathione-S-transferases and UDP-glucuronosyltransferases. Finally, using a RAD-tag approach, we investigated the population genomics of this pest, looking at its likely patterns of invasion.

1. Introduction

The light brown apple moth, Epiphyas postvittana, is a polyphagous, invasive tortricid moth [1]. Larvae have been recorded from more than 500 host plant species in 363 genera, from 121 families across vascular plants [2]. The ability to tolerate the toxicity associated with many of the chemical components of leaves from different plants has been attributed to the diversity of detoxification enzymes expressed in the larval gut [3,4]. In Australia and New Zealand, E. postvittana is an important horticultural pest, mostly of pome fruit, citrus and grapes, and its pest status is a major concern in the USA [5]. In Australasia, the impacts on horticulture are through both direct damage to crops and the market access of exported products. Epiphyas postvittana is endemic to mainland Australia, where it is distributed throughout south eastern parts of the landmass, but it has been introduced into Tasmania and countries, such as New Zealand, the USA (Hawaii and California) and northern Europe [1]. The timings of these successful invasions, outside Australia, are thought to vary from the late nineteenth century for New Zealand, early twentieth century for Hawaii, to the late twentieth century for Europe and the early twenty-first century for California. Efforts to determine the origin of the Californian invasion using mitochondrial sequence analysis found evidence for two separate incursions and ruled out a UK or Hawaiian origin for either [6,7].
Due to its pest status and recurrent incursions, there are ongoing efforts to control E. postvittana [1]. Insecticides have been an important tool, although resistance, withdrawal of insecticides and the development of alternative measures and Integrated Pest Management programmes have reduced their use. Biopesticides, such as Bacillus thuringiensis, Nuclear Polyhedrosis Virus (NPV) and spinosad are effective but not widely used. A number of classical biological control introductions (137) have been made into New Zealand and the Sterile Insect Technique has been considered for controlling the pest [8]. Pheromone-based tactics have been developed for both monitoring and control, using methods such as mating disruption, mass trapping and lure-and-kill. The advent of new gene technologies provides opportunities to develop novel control tactics, such as RNAi [9].
Members of the Lepidoptera (moths and butterflies) have become some of the most serious pests of primary production globally. Their destructive nature is often associated with their ability to complete their life cycle on many different host plants, invade new territories and to be resistant to methods to control them [10,11]. Understanding the biology associated with polyphagy, invasiveness and xenobiotic resistance is now being informed by genomic sequencing of pest species, revealing the genetic basis for these important traits [12]. Indeed, it may well become possible to predict the potential species likely to invade new areas or become pests based on their genetics [13]. More and more genomes of lepidopteran pests are becoming available [14,15,16,17,18], with many being highly polyphagous and invasive. Comparisons of these genomes have identified gene families that are likely to underpin polyphagy, enable invasiveness and confer insecticide resistance [16,18,19]. For example, both polyphagy and invasiveness have been associated with expansion of gene families involved in chemosensation and detoxification [14,16,19,20].
The Tortricidae is one of the largest families of microlepidopteras, which comprise more than 10,000 described species worldwide. The phylogenetic relationships of the Tortricidae, also in relation to the other Lepidoptera, have been studied by Regier and colleagues [21,22]. The family includes a number of important pests of horticulture. Genomic sequences have been deposited in Genbank for eight species, of which six from 2021 are yet to have an accompanying publication. The two genomes with publications are both from pests with relatively narrow host ranges, with codling moth (Cydia pomonella) being a pest of apple and related Rosaceae and the smaller tea tortrix (Adoxophyes honmai) being a pest of tea. Within the codling moth genome, Wan et al. [17] identified a duplication in an odorant receptor involved in pheromone and host plant recognition, and an association between insecticide resistance and a particular cytochrome P450 (CYP450). From the smaller tea tortrix genome, Uchibori-Asano et al. [23] identified target site resistance and enhanced detoxification by carboxylesterases and CYP450s underlying tebufenozide resistance. To date, however, no highly polyphagous pest within the tortricid family has been genome sequenced.
Here, we describe a genome assembly for E. postvittana and compare it with other moth genomes, use it to investigate gene families potentially involved in key pest traits and to inform a population genetic analysis of its invasion patterns.

2. Materials and Methods

2.1. Genome Sequencing

Epiphyas postvittana were from a laboratory colony that has been maintained at Plant & Food Research and its predecessors since 1967. The colony was established from a wild collection from Nelson, New Zealand, with an addition of wild males at generation 128. Larvae were reared on a general all-purpose diet [24]. Pupae were kept at room temperature and separated by sex before eclosion. Adults were provided with sugar and water.
Genomic DNA was extracted from purified nuclei produced using a method based on that of Naim et al. [25] starting from three grams of female pupae (about 50 individuals). We chose only females since in Lepidoptera, females are the usually heterogametic sex allowing a better representation of all chromosomes. Up to 0.2 g of purified nuclei were resuspended in 14 mL of lysis buffer (0.5% SDS, 5 mM EDTA, 150 mM Tris-borate pH 7.4), treated with RNase, DNA purified by chloroform:isoamyl alcohol (24:1) extraction, precipitated with isopropanol and washed with 70% ethanol. The detailed protocol is available at (https://dx.doi.org/10.17504/protocols.io.rncd5aw, accessed on 10 January 2022). DNA quality was assessed at OD 260/280 by nanodrop (Thermo Scientific, Waltham, MA, USA). Mitochondrial DNA contamination was assessed by PCR using primers specific for the mitochondrial cytochrome oxidase I gene relative to that of the nuclear gene Takeout 3.
Paired-end (PE) and both mate pair (MPE) and Nextera MPE libraries were constructed where the PE insert size varied between 169 bp and 264 bp, and MPE insert size ranged between 2.5 kb and 8.5 kb. Thereafter sequencing was carried out at the Allan Wilson Centre for Genome Sequencing (Palmerston North, New Zealand) and by Macrogen (Seoul, Korea) according to Illumina instructions. Raw paired-end and mate pair reads are available in the SRA archive, with a bioproject ID of PRJNA754242.

2.2. Genome Assembly and Assessment

The quality of resulting reads was assessed using FastQC (version 0.9.1). PE data were subject to quality and adapter trimming using the fastq-mcf tool of the ea-util suite (version 1.04.807) [26]. Further quality trimming was carried out using an in-house Perl script to remove reads with Ns and mononucleotides. MPE libraries were quality trimmed using the fastq-mcf tool of the ea-util suite (version 1.04.807), with read length set to 36 bp and PE contamination was removed using Bowtie (version 1.1.1) [27]. Removal of duplicates in both PE and MPE libraries was carried out using an in-house Perl script (available at https://github.com/rosscrowhurst/Perl/blob/fd49a32bd1ef14eca4f085f136afc3b656977547/Fastq/remove_duplicate_read_pairs_orig.pl, accessed on 10 January 2022). MPE data from the Nextera protocol were quality trimmed and duplicates removed using the NextClip tool (version 20140619) [28].
The characteristics of the genome were explored using GenomeScope (version 1) [29], where short-read PE libraries were used as input. Assembly of the genome was carried out using the quality trimmed reads with two assemblers, ALLPATHS-LG (version 50191) [30] and Platanus (version 1.2.4) [31]. The Platanus assembly was further scaffolded with SSPACE (version 2.0) [32] and a round of gap closing using the PE data was performed with GapCloser for SOAPdenovo (version 1.12) [33]. Thereafter both assemblies were checked for plant and microbial contaminants using Kraken (version 0.10.5) [34] where a custom database was made using the invertebrate, plant and microbial RefSeq databases. Subsequently, the NCBI ‘nt’ database was used to further classify the unclassified scaffolds from the first iteration. Based on the Kraken and the NCBI blast classifications, the scaffolds classified as ‘arthropoda’ and ‘unclassified’ were retained and the remainder were discarded as potential contaminants. Finally, Metassembler (version 1.3) [35] was used to merge the free-of-contaminants Platanus and Allpaths assemblies together, with the Platanus assembly set as the primary assembly for the merging.
Completeness of the E. postvittana genome was assessed with BUSCO (version 3.0.2) [36], using the ‘arthropoda_odb9′ database, which contains 1066 BUSCO markers and was run in the ‘long’ mode which turns on the Augustus optimization mode for self-training. To further test for completeness of the genome and to investigate the high number of gene models compared with those of other lepidopteran genomes, we mapped the light brown apple moth antennal transcriptomes [37,38] on the genome with Magic-BLAST [39] using default parameters. In addition, genomic reads were mapped back to the genome using Bowtie 2 (v2.3.4.3) [40].
The pattern of coverage and GC% were explored using BlobTools v1.0 [41]. Assembly contigs and scaffolds were blasted against NCBI nt database (E-value 1 × 10−10).

2.3. Genome Annotation

RepeatModeler (v 1.0.8) (http://www.repeatmasker.org; access date 4 March 2022) was used to predict and classify repetitive elements. This tool employs two de novo repeat finding programs, RECON and RepeatScout, to identify repeat element boundaries, followed by assigning repeats to the repeat classes based on the sequence feature. The resulting repeat models were searched against the GenBank non-redundant (nr) protein database (downloaded on 26 November 2015; E-value < 1 × 10−5) using Blastx (v2.2.30) to exclude potential protein-coding genes. The abundances of all predicted repeats were estimated in the genome assembly with RepeatMasker (v4-0-5).
Genome annotation was carried out using AUGUSTUS (version 3.3) [42], where in-house expressed sequence tags (ESTs) [4,37] and repeats were used as hints for the annotation with the BUSCO-trained model generated from the light brown apple moth genome fed in as “--species”. The protein sequences from the Augustus-predicted models were blasted against a protein database containing Drosophila melanogaster, Tribolium castaneum and the following lepidopteran species: Amyelois transitella, Bicyclus anyana, Bombyx mori, Chilo supressalis, Helicoverpa armigera, Papilio machaon, Papilio polytes, Papilio xuthus, Phoebis sennae, Pieris rapae, Plutella xylostella and Spodoptera litura. Blastp results were fed to OmicBox where mapping and annotation were performed with default parameters.

2.4. Comparative Genomics

To identify orthologous groups among E. postvittana and other insects we used OrthoFinder [43]. We retrieved protein sequences from the genome of 11 Lepidoptera including one from the superfamily Yponomeutoidea (P. xylostella, family: Plutellidae), two from Bombycoidea (B. mori, Bombycidae and Manduca sexta Sphingidae), three Noctuoidea (H. armigera, S. litura and T. ni, all from the family Noctuidae), three Papilionoidea (P. xuthus, Danaus plexippus, Papilionidae and two from the family Nymphalidae: Heliconius melpomene, P. rapae) and a Tortricoidea (C. pomonella, Tortricidae). We retrieved the protein data from the genome of a further 10 other insects (two Diptera: Anopheles gambiae and Drosophila melanogaster; two Coleoptera: Leptinotarsa decemlineata and T. castaneum; two Hymenoptera: Nasonia vitripennis and Apis mellifera, two Hemiptera: Diaphorina citri and Bemisia tabaci; two Blattodea: Blatella germanica and Cryptotermes secundus and one Collembola, (Folsomia candida, Entomobryomorpha). The choice of species was mostly based on the availability of annotated genomes at the time we preformed the analysis, therefore the most recently deposited six Tortricid genomes in Genbank from 2021 were not included. Protein sequences were extracted from genome sequences of each species and cleaned from multiple isoforms. Orthofinder was run, specifying the multiple sequence alignment (-M msa).
Single-copy orthologues shared by all species were aligned with MAFFT [44] and the resulting alignments were trimmed with trimAl v1.2 [45] to remove poorly aligned regions using the parameter “-automated1”. Single-copy aligned orthologues were merged in a multi-aligned supergene in Geneious. A maximum likelihood tree was obtained with IQ-TREE [46] using the LG+F+R6 substitution model calculated by ProtTest in IQ-TREE. Folsomia candida was used as the outgroup.
To identify gene families that significantly changed as a result of a stochastic birth and death process in each branch of the tree we applied the maximum likelihood method implemented in CAFE ver. 5.0 [47]. Using the number of genes per family obtained with OrthoFinder and the IQ-TREE tree as input files, we calculated the average gene expansion and contraction and identified the gene families with significant (p < 0.05) size changes (gain and loss of genes) accounting for the species phylogenetic history.

2.5. Population Genomics

Ten individuals, both male and female, were sampled from each of ten populations of E. postvittana. Two populations were from the North Island of New Zealand, Te Puke (TPU) and Hawke’s Bay (HB); three from the South Island of New Zealand, Motueka (MOT), Awatere (AWA) and Clyde (CLY); four populations from Australia, Wagga Wagga in New South Wales (NWS), Dalkeith in Western Australia (WA), Gumeracha in South Australia (GUM), one in Tasmania (TAS); and one population was from California, USA (CA).
DNA was extracted from each individual and pooled according to population in equal amounts. Each mixture of DNA from 10 moths was cut with Sbf1 for the RAD-tag library to subsample the genome and then sequenced on an Illumina Hiseq.
Genetic diversity measures for each population, such as nucleotide diversity π, Watterson’s θ and Tajima’s D, were calculated from the pooled samples of each population of E. postvittana using the package PoPoolation [48] following the walkthrough available with the software. Briefly, the cleaned reads of each population were mapped to the indexed reference genome (LBAM_v1.0) using bwa [49] with the following parameters (bwa aln -n 0.01 -l 100 -o 1 -d 12 -e 12 -t 8). The resulting sai files were converted to sam files with the function samse of bwa. Ambiguously mapped reads were removed with samtools [50] using the parameters (view -q 20 -bS). The population genetic measures were calculated using a sliding window approach along the genome scaffolds with min-coverage = 10 and min-quality = 20.
Differentiation between populations, expressed as Fst values [51], was calculated for the whole genome as well as for genes and intergenic regions using the package PoPoolation2 [52]. Significance of Fst values weas estimated per gene and intergenic regions by Fisher’s exact tests using PoPoolation2. The parameters used in all Fst calculations were min-count 6, min-coverage 20, max-coverage 2%, pool size 10. Population relationships were also investigated by inferring a maximum likelihood phylogenetic tree with TreeMix [53], grouping SNPs in windows of size 10, since SNPs were very sparse, migration to 3 estimated using the R package OptM [54], and with 100 bootstrap replicates obtained with the R package BITE [55]. The tree was rooted with the Western Australia population of Dalkeith.

3. Results

3.1. Genome Sequencing and Analyses

Genome characteristics using GenomeScope suggested a genome size of 360–437 Mb, with a heterozygosity > 1% (Figure 1).
Furthermore, the k-mer plots suggested the presence of a high repeat content in the genome, 20–25%. The initial Allpaths assembly size was 539 Mb and over 23,000 scaffolds, with a N50 of 174 kb, while the Platanus assembly was 660 Mb and over 34,000 scaffolds, with a N50 of 248 kb. A contamination check with Kraken, using RefSeq databases, revealed that a high portion of scaffolds were unclassified (~64% in both assemblies) and about 30% of the scaffolds were of arthropod origin, with Bombyx mori containing the highest portion (Supplementary Figure S1). Although it appeared that a high portion were categorised as unclassified in relation to assembly size, 493 Mb of the Allpaths and 594 Mb of the Platanus assembly were of arthropod origin, suggesting over 90% of the assembly comprised insects, with a low presence of contaminants. Subsequently, a further 9–11% of the unclassified scaffolds were classified as arthropod using the ‘nt’ database. After contaminant removal, the Allpaths assembly size was 521.6 Mb, with over 20,000 scaffolds, and a N50 of 196 kb, while the Platanus assembly was 638 Mb, with over 32,000 scaffolds, and a N50 of 265 kb. The final merged light brown apple moth genome assembly from Platanus and Allpaths had a total length of over 598 Mb, consisting of over 14,000 scaffolds (Table 1), with a contig N50 of 26.2 kb and a scaffold N50 of 301.17 kb. Assembly length was consistent with the genome size determined by flow cytometry, where the 1C-value was determined to be 0.62 pg, corresponding to an estimated haploid genome size of 606 Mb. Genome completeness level, estimated by BUSCO analysis, was 97.9%, with a duplication rate of 18.3% and 0.6% fragmented, indicating that good representation of the E. postvittana genome has been achieved. Furthermore, a high proportion of reads (84–93.5%) were observed to map back to the assembly using Bowtie 2 (v2.3.4.3).
Epiphyas postvittana has 34% of its assembled genome represented as repetitive sequences, with the majority of the known elements made up of longer DNA transposable elements (14.07 Mb) and retrotransposons (LINE 17.83 Mb) (Table 2). The genome assemblies of E. postvittana and P. xylostella have a similar percentage of the genome represented as repeats, with H. armigera having substantially less of the genome represented by repetitive sequences. The LINE family is the most abundant of the modelled element classes in both species. LTRs are more abundant in P. xylostella, while SINEs are more abundant in E. postvittana. There is a high percentage of unclassified repeats, in both E. postvittana and P. xylostella, compared with H. armigera, indicating that novel repeat classes are yet to be discovered in these species.

3.2. Genome Annotation

AUGUSTUS gene prediction resulted in 31,389 de novo predicted genes. This is high compared with those predicted in most other Lepidoptera and more similar to that found in the German cockroach, Blatella germanica (more than 28,000 protein-coding genes) [58], in the parthenogenetic springtail F. candida [59] and in the lepidopteran Heliconius melpomene [60], all with more than 20,000 protein-coding genes. To further investigate this result, we mapped the 23,591 non-overlapping codifying sequences from the antennal transcriptome of E. postvittana male and female [38]. As a result, 90.5% (21,350) of the coding sequences mapped to the genome, with 1406 (6.58%) of them mapping more than once.
The number of gene models with a blast hit against the NCBI protein database was 27,793, of which 27,003 (86%) contained an InterPro domain and 12,276 (38.7%) were assigned a GO term (Supplementary Table S1). No blast hits were recovered for 3596 of the gene models.

3.3. Comparative Genomics

Orthologue analysis with OrthoFinder identified 207 single-copy genes among the 22 insect species and F. candida. The phylogenetic tree from the concatenation of the 207 groups of aligned sequences (80,652 amino acids) is shown in Figure 2.
Of the 31,389 genes predicted in the E. postvittana genome, 28,714 (91.5%) were assigned to 11,438 orthogroups, of which 945 were specific to E. postvittana, with 3700 genes predicted to be unique to this species (Figure 2). Over one-third (14,376) of all E. postvittana genes share orthologues with the other lepidopteran species examined and 661 of these, in 275 orthologous groups, were unique to members of the Lepidoptera. A full list of the orthologous groups is provided in Supplementary Table S2. Statistical analysis of the evolution of the size of the gene families, using Cafe, given the ML phylogenetic tree and families sizes obtained with Orthofinder, predicted an increase in gene number in E. postvittana, across 2397 gene families and a decrease in 1033 from the common tortricid ancestor with C. pomonella (Figure 2). A further 903 gene families showed an increase in gene number in the tortricid cluster, which remained constant in E. postvittana, with only 41 showing a decrease (Figure 2).
Twenty gene families were significantly expanded in E. postvittana, including genes for larval cuticle proteins, genes involved in detoxification, such as cytochrome P450s and glutathione S-transferases (GST), genes involved in transcriptional regulation, such as zinc finger proteins, and genes involved in defence, such as serine protease and macrophage mannose receptors. Among the gene families involved in detoxification, we found almost 300 P450 gene models in 45 orthogroups in the E. postvittana genome, with two orthologous groups that showed significant expansion, the CYP 6B and CYP 9E subclasses. Among phase II metabolizing enzymes, more than 50 cytosolic GSTs in 14 orthogroups were detected in the E. postvittana genome, with one expanded group of 13 genes, belonging to the sigma class, and more than 100 UDP-glucuronosyltransferase in 16 orthogroups, of which two were expanded in E. postvittana.
Only three gene families showed a significant contraction in E. postvittana, of which two are associated with transposable elements. The full list is shown in Supplementary Table S2.

3.4. Population Genomics

Resequencing of one hundred individuals from ten populations revealed 18,767 SNPs, representing 0.0031% of the 598 Mb assembled genome—the variome. The majority of the SNPs were biallelic, found in genes (57.56%) (Table 3).
The number of SNPs per population and population genetic parameters are shown in Table 4.
The lowest number of SNPs, as well as the lowest nucleotide diversity (π), were found in Te Puke (3669 SNPs and 7.12 × 10−4, respectively), whereas the highest number of SNPs were found in California (7934), and the highest π in Gumeracha (1.24 × 10−3). The lowest population mutation rate (θ) was from Clyde (8.14 × 10−4), whereas the highest was in Gumeracha (1.48 × 10−3). Tajima’s D values were only slightly negative for all populations. There were no significant differences in any population genetic parameters among countries (Krustall–Wallis test: π: χ2 = 4.55, df = 2, p > 0.05, θ: χ2 = 3, df = 2, p > 0.05, D: χ2 = 1.59 df = 2, p > 0.05).
Average Fst among all populations was 0.146. Genetic differentiation was similar among populations, with the lowest Fst value being between the two North Island populations of New Zealand, Hawke’s Bay and Te Puke (0.086) (Table 5).
Of the 31,389 predicted genes from the assembled E. postvittana genome, 1025 (4.5%) contained SNPs and were informative for estimating the differentiation among populations (Figure 3). However, population differentiation was similar when calculated using SNPs within genes versus within intergenic regions (0.147 and 0.144, respectively) (Supplementary Table S3), and their pairwise population values were highly correlated (r2 = 0.914). Of these 1025 genes, 51 were in the top 5% of the Fst distribution (Fst > 0.37) and 50 of these genes were among the annotated genes (Supplementary Table S4).
No GO enrichment was detected among the outlier genes. Genetic differentiation among populations was correlated with their geographic distances, whether considering all SNPs (Corr. = 0.970321, p = 0.05) or either genic or intergenic regions (Corr. = 0.966436, p = 0.05; Corr. = 0.974291, p = 0.05, respectively).
In a phylogenetic analysis (Figure 4), using TreeMix, the New Zealand populations formed two clusters, with the two most northern populations from the South Island (MOT and AWA), sister to the Tasmanian population. The most southern population from the South Island (CLY) was sister to the two North Island populations (TPU and HBY). The Californian population, clustered together with the eastern Australian populations (GUM, TAS and NSW) and two of the South Island populations of New Zealand (AWA and MOT). Migration events were predicted between the South Island populations of New Zealand and between Wagga Wagga in New South Wales and Gumeracha in South Australia.

4. Discussion

In this paper we present the first genomic resources for the light brown apple moth, E. postvittana, an invasive and polyphagous lepidopteran pest of horticultural systems. We further used this genomic resource to investigate the expansion of gene families putatively linked to the evolution of polyphagy and the population genomics of this invasive species.
The predicted size from the genome assembly of 598 Mbp is similar to the genome size determined by flow cytometry of 606 Mbp, and is within the range of those of other Lepidoptera sequenced to date. The size of the E. postvittana genome does not seem due to an expansion of repetitive elements, which make up 33.8% of the genome (Table 2), a proportion lower than in B. mori and P. xylostella, which both have smaller genomes than E. postvittana. Additionally, the percentage of different types of repeat elements was quite similar to that found in other lepidopteran genomes, although there were a high number of unclassified repeats (26%), as was reported for the P. xyolstella genome (28%) [15]. The K-mer profiles (Figure 1) suggested high heterozygosity and high repetitiveness, probably due to the high number of individuals used to construct the sequencing libraries, suggesting difficulty in assembly, although this genome assembly is only based on Illumina reads. However, high and similar levels of heterozygosity have been found from the genome sequencing of a single individual of the inbred clonal honey bee, Apis mellifera capensis [61], as well as other insects we have sequenced (unpublished data). Nevertheless, additional resolution of the E. postvittana genome would probably benefit from the newer long-read technologies, as well as using high molecular weight genomic DNA extracted from a single insect.
In comparison with genome size, the gene number is reasonably high compared with those of other lepidopteran genomes (Table 1), which may well be related to genome assembly only being based on Illumina reads. However, the number of genes that are expressed in the antennae of both sexes [38] was also higher than in the gene models of other lepidopteran genomes, and most of these CDS (21,350; 90.5%) mapped to our genome assembly. This evidence may suggest that the E. postvittana genome could contain a higher number of genes than those of other sequenced members of the Lepidoptera.
The evolution of highly polyphagous species has been linked to the expansion of gene families involved in detoxification pathways, from active metabolism of toxic plant compounds to enhanced excretion [14,16,19,20]. These same gene families have also been implicated in the development of resistance to pesticides in many insect species, suggesting a tight link between polyphagy and ability to develop resistance to xenobiotics [62,63]. We found several gene families that significantly changed in copy number, in the lineage leading to E. postvittana, with expansions in gene families involved in detoxification, defence and cuticle structure. Among the families involved in detoxification, we found expansion in phase I metabolism genes, including two CYP450 families, and phase II metabolism genes, including members of the sigma class of GSTs and two UDP-glucuronosyltransferase (UGT) families. While it is difficult to determine whether these expansions are biologically real or artefacts of the assembly process, given the source DNA was from multiple individuals, the observed expansions are in gene families known to be expanded in other Lepidoptera [18,19,56,64]. Furthermore, the antennal transcripts that mapped multiple times on the E. postivittana assembled genome did not involve the GST family, and in only one case, the P450 6B and UDP orthogroups, out of 38 and 23 changes in gene copy number in E. postvittana, respectively.
For CYP450s, there is good evidence that new copies appear when herbivores colonize new host plants [65,66,67], supporting the hypothesis that an arms race between herbivorous insects and their host plants drives P450 diversification. Of the two expanded groups of CYP450 in E. postvittana, members of the 6B subfamily are able to metabolize linear furanocoumarin (xanthotoxin and bergapten) in the lepidopterans Papilio polyxenes and P. glaucus [68], whereas members of the A and E subfamilies are highly upregulated in insecticide-resistant strains of Aedes aegypti and Culex pipiens [69,70].
GSTs have an important physiological role in lipid peroxidation by eliminating by-products of oxidative stress, as well as roles in insecticide resistance and in the breakdown of plant protective compounds [71,72]. Sigma GSTs are ubiquitous among metazoans and are one of the most abundant and conserved GST classes in insects [73]. The expansion of the sigma class in the light brown apple moth is a new report amongst the Lepidoptera, where typically expansions of the epsilon and delta classes have been observed [20]. Both the epsilon and delta classes have been implicated in detoxification and resistance to insecticides [19,71]. These two classes contained 20 and 7 gene copies, respectively, in E. postvittana.
UGTs are membrane-bound phase II metabolizing enzymes that, in insects, conjugate glucose to various endogenous and exogenous substrates, playing an important role in the detoxification of plant allelochemicals by promoting their excretion [74,75]. Most UDPs are highly expressed at larval stages, indicating their importance in the protection against plant secondary metabolites [76]. However, UGTs have also been implicated in insecticide resistance because of their activity on xenobiotics [77,78] and, in addition, they play roles in several other processes, including cuticle formation and olfaction [74]. The two UDP families that are expanded in E. postvittana, UDP33 and 40, are the two largest UDP families in both B. mori and P. xylostella, with the expression of most of their members being affected by insecticide treatments [78].
Mechanisms to minimize the effects of ingestion of plant proteinase inhibitors in insects often involve the upregulation of serine peptidases of the chymotrypsin and trypsin family [79]. In Spodopera frugiperda, chronic exposure to proteinase inhibitors results in the upregulation of 16 chymotrypsin and trypsin genes. The responsive trypsin genes belong to a separate clade from unresponsive genes, supporting a possible mechanism of resistance involving gene number expansion [80]. This same gene family has also been implicated in insecticide resistance [81].
In other highly polyphagous pests, such as H. armigera and S. litura, an expansion of gustatory receptors has been observed, which is thought to allow the detection of the many host defence response compounds [19]. This was not seen in E. postvittana, however, we did notice that manual efforts to annotate this gene family resulted in higher estimates in the number of gustatory receptor genes in E. postvittana (unpublished data), in line with other polyphagous Lepidoptera [19]. The annotation of odorant receptors genes and other genes involved in odorant reception has been covered separately [37].
The light brown apple moth is a global pest. From Australia, it invaded New Zealand in the late 19th century [82], and soon after, colonized Hawaii, between 1900 and 1925 [83], then England, where it arrived in 1936 [84]. More recently it has invaded Ireland (1968), Sweden and the Netherlands, and in 2006, California [1,85,86]. The resequencing of pooled samples from ten E. postvittana populations found more than 18,000 SNPs, of which more than 50% were within gene models. Although the number of pooled individuals was low for a very accurate estimation of variant frequency, genomic population results were similar to those obtained from a survey of mtDNA variation [7]. All populations were similarly diverse, including those in Tasmania, which showed higher pairwise Fst than observed for mtDNA [7]. Similar to the mtDNA results, California showed the same degree of diversity as Australia and New Zealand, despite its more recent origin. On the basis of the mtDNA haplotypes of the 2006/2007 samples, Tooman et al. [7] suggested two probable invasions in California. The high genetic diversity in the 2007 samples also indicated either a source that was highly diverse genetically, or multiple sources, for which both Australia and New Zealand fit the profile. Our genomic data identify Australia as a probable source of the California invasion, as suggested by population relationships on the ML tree, in which California clusters with the population of Wagga Wagga and Gumeracha in South Eastern Australia (Figure 3). New Zealand was likely invaded twice, as shown by the two distinct clusters in Figure 3, and with signs of admixture of populations from the South Island to the North Island. Signs of admixture were also evident between the two closest populations in Australia, Wagga Wagga and Gumeracha. Nevertheless, this result should be taken with caution, since both the sample size within populations and the number of sampled populations are low.
In the USA, the light brown apple moth is classified as a Class A pest. Recent models predict its potential invasive range reaching most of the coastal areas of the globe, USA and central Mexico, South America in Chile, Perù and Argentina, most of western Europe, and all the Mediterranean basin, the south of the African continent and Southeast Asia, including inland areas [87]. The highly polyphagous nature of E. postvittana, its high rates of genetic diversity, and predictions of its potential range, including many horticulturally important growing regions, raises the degree of concern about this species. As such, the genomic resource presented here should be useful in developing new management tools for this pest and for understanding the biological basis for this and other species’ ability to invade new regions.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/insects13030264/s1, Table S1: Blast description, GO and InterPro annotation of gene models obtained from Augustus gene prediction; Table S2: List of all orthologous groups recovered by OrthoFinder with a representative gene model, where present, of Epiphyas postvittana and its blast description, GO and InterPro annotation. For each orthologous group the number of gene changes in the gene tree are indicated; Table S3: Average Fst calculated on genic SNPs (above the diagonal) and on the intergenic SNPs (below the diagonal); Table S4: List of 50 genes in the highest 5% of the Fst distribution and their relative blast hit, description and annotation; Figure S1: Krona plot illustrating the taxonomic classification of the scaffolds using the RefSeq sequences as the reference database.

Author Contributions

Conceptualization, R.D.N. and A.G.; data curation, A.H.T.; formal analysis, J.T.C., R.M.S., L.K.T., M.D.J. and R.C.; funding acquisition, R.D.N.; investigation, A.H.T., C.W. and A.G.; methodology, E.H. and D.B.; writing—original draft, A.H.T., C.W., R.D.N. and A.G.; writing—review and editing, R.M.S., D.B., R.D.N. and A.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by The Foundation for Research, Science, and Technology of New Zealand (CO6X0701, CO6X0207, and C06X0301).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

All raw sequencing data have been deposited in the SRAS archive, with bioproject ID PRJNA754242 in the NCBI database.

Acknowledgments

We thank Anne Barrington for insect rearing, Floragenex for conducting the RAD library construction and sequencing and Ting-Hsuan Chen and Cecilia Deng for useful comments on the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Suckling, D.M.; Brockerhoff, E.G. Invasion biology, ecology, and management of the light brown apple moth (Tortricidae). Annu. Rev. Entomol. 2010, 55, 285–306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Brockerhoff, E.G.; Suckling, D.M.; Ecroyd, C.E.; Wagstaff, S.J.; Raabe, M.C.; Dowell, R.V.; Wearing, C.H. Worldwide host plants of the highly polyphagous, invasive Epiphyas postvittana (Lepidoptera: Tortricidae). J. Econ. Entomol. 2011, 104, 1514–1524. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Christeller, J.T.; Poulton, J.; Markwick, N.M.; Simpson, R.M. The effect of diet on the expression of lipase genes in the midgut of the lightbrown apple moth (Epiphyas postvittana Walker; Tortricidae). Insect Mol. Biol. 2010, 19, 9–25. [Google Scholar] [CrossRef] [PubMed]
  4. Simpson, R.M.; Newcomb, R.D.; Gatehouse, H.S.; Crowhurst, R.N.; Chagné, D.; Gatehouse, L.N.; Markwick, N.P.; Beuning, L.L.; Murray, C.; Marshall, S.D.; et al. Expressed sequence tags from the midgut of Epiphyas postvittana (Walker) (Lepidoptera: Tortricidae). Insect Mol. Biol. 2007, 16, 675–690. [Google Scholar] [CrossRef]
  5. Brown, J.W.; Epstein, M.E.; Gilligan, T.M.; Passoa, S.C.; Powell, J.A. Biology, identification, and history of the light brown apple moth, Epiphyas postvittana (Walker) (Lepidoptera: Tortricidae: Archipini) in California: An example of the importance of local faunal surveys to document the establishment of exotic insects. Am. Entomol. 2010, 56, 34–43. [Google Scholar] [CrossRef] [Green Version]
  6. Rubinoff, D.; Holland, B.S.; Jose, M.S.; Powell, J.A. Geographic proximity not a prerequisite for invasion: Hawaii not the source of california invasion by light brown apple moth (Epiphyas postvittana). PLoS ONE 2011, 6, e16361. [Google Scholar] [CrossRef] [Green Version]
  7. Tooman, L.K.; Rose, C.J.; Carraher, C.; Suckling, D.M.; Paquette, S.R.; Ledezma, L.A.; Gilligan, T.M.; Epstein, M.; Barr, N.B.; Newcomb, R.D. Patterns of mitochondrial haplotype diversity in the invasive pest Epiphyas postvittana (Lepidoptera: Tortricidae). J. Econ. Entomol. 2011, 104, 920–932. [Google Scholar] [CrossRef]
  8. Kean, J.M.; Suckling, D.M.; Stringer, L.D.; Woods, B. Modeling the sterile insect technique for suppression of light brown apple moth (Lepidoptera: Tortricidae). J. Econ. Entomol. 2011, 104, 1462–1475. [Google Scholar] [CrossRef]
  9. Turner, C.T.; Davy, M.W.; MacDiarmid, R.M.; Plummer, K.M.; Birch, N.P.; Newcomb, R.D. RNA interference in the light brown apple moth, Epiphyas postvittana (Walker) induced by double-stranded RNA feeding. Insect Mol. Biol. 2006, 15, 383–391. [Google Scholar] [CrossRef]
  10. Jurat-Fuentes, J.L.; Heckel, D.G.; Ferré, J. Mechanisms of resistance to insecticidal proteins from Bacillus thuringiensis. Annu. Rev. Entomol. 2021, 66, 121–140. [Google Scholar] [CrossRef]
  11. Richardson, E.B.; Troczka, B.J.; Gutbrod, O.; Davies, T.G.E.; Nauen, R. Diamide resistance: 10 years of lessons from lepidopteran pests. J. Pest Sci. 2020, 93, 911–928. [Google Scholar] [CrossRef] [Green Version]
  12. Tay, W.T.; Gordon, K.H.J. Going global—genomic insights into insect invasions. Curr. Opin. Insect Sci. 2019, 31, 123–130. [Google Scholar] [CrossRef]
  13. Huang, C.; Wang, S.; Fan, X.; Pian, C.; Luo, J.; Li, X.; Lang, K.; Xing, L.; Jiang, M.; Liu, W.; et al. Predicting insect invasiveness with whole-genome sequencing data. BMC Genom. 2020. underreview. [Google Scholar] [CrossRef]
  14. Chen, W.; Yang, X.; Tetreau, G.; Song, X.; Coutu, C.; Hegedus, D.; Blissard, G.; Fei, Z.; Wang, P. A high-quality chromosome-level genome assembly of a generalist herbivore, Trichoplusia ni. Mol. Ecol. Resour. 2019, 19, 485–496. [Google Scholar] [CrossRef]
  15. Jouraku, A.; Yamamoto, K.; Kuwazaki, S.; Urio, M.; Suetsugu, Y.; Narukawa, J.; Miyamoto, K.; Kurita, K.; Kanamori, H.; Katayose, Y.; et al. KONAGAbase: A genomic and transcriptomic database for the diamondback moth, Plutella xylostella. BMC Genom. 2013, 14, 464. [Google Scholar] [CrossRef] [Green Version]
  16. 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]
  17. Wan, F.; Yin, C.; Tang, R.; Chen, M.; Wu, Q.; Huang, C.; Qian, W.; Rota-Stabelli, O.; Yang, N.; Wang, S.; et al. A chromosome-level genome assembly of Cydia pomonella provides insights into chemical ecology and insecticide resistance. Nat. Commun. 2019, 10, 4237. [Google Scholar] [CrossRef] [Green Version]
  18. Xiao, H.; Ye, X.; Xu, H.; Mei, Y.; Yang, Y.; Chen, X.; Yang, Y.; Liu, T.; Yu, Y.; Yang, W.; et al. The genetic adaptations of fall armyworm Spodoptera frugiperda facilitated its rapid global dispersal and invasion. Mol. Ecol. Resour. 2020, 20, 1050–1068. [Google Scholar] [CrossRef]
  19. Cheng, T.; Wu, J.; Wu, Y.; Chilukuri, R.V.; Huang, L.; Yamamoto, K.; Feng, L.; Li, W.; Chen, Z.; Guo, H.; et al. Genomic adaptation to polyphagy and insecticides in a major East Asian noctuid pest. Nat. Ecol. Evol. 2017, 1, 1747–1756. [Google Scholar] [CrossRef]
  20. Gouin, A.; Bretaudeau, A.; Nam, K.; Gimenez, S.; Aury, J.-M.; Duvic, B.; Hilliou, F.; Durand, N.; Montagné, N.; Darboux, I.; et al. Two genomes of highly polyphagous lepidopteran pests (Spodoptera frugiperda, Noctuidae) with different host-plant ranges. Sci. Rep. 2017, 7, 11816. [Google Scholar] [CrossRef]
  21. Regier, J.C.; Brown, J.W.; Mitter, C.; Baixeras, J.; Cho, S.; Cummings, M.P.; Zwick, A. A Molecular Phylogeny for the Leaf-Roller Moths (Lepidoptera: Tortricidae) and Its Implications for Classification and Life History Evolution. PLoS ONE 2012, 7, e35574. [Google Scholar] [CrossRef]
  22. Regier, J.C.; Mitter, C.; Zwick, A.; Bazinet, A.L.; Cummings, M.P.; Kawahara, A.Y.; Sohn, J.-C.; Zwickl, D.J.; Cho, S.; Davis, D.R.; et al. A Large-Scale, Higher-Level, Molecular Phylogenetic Study of the Insect Order Lepidoptera (Moths and Butterflies). PLoS ONE 2013, 8, e58568. [Google Scholar] [CrossRef]
  23. Uchibori-Asano, M.; Jouraku, A.; Uchiyama, T.; Yokoi, K.; Akiduki, G.; Suetsugu, Y.; Kobayashi, T.; Ozawa, A.; Minami, S.; Ishizuka, C.; et al. Genome-wide identification of tebufenozide resistant genes in the smaller tea tortrix, Adoxophyes honmai (Lepidoptera: Tortricidae). Sci. Rep. 2019, 9, 4203. [Google Scholar] [CrossRef]
  24. Singh, P. A general purpose laboratory diet mixture for rearing insects. Int. J. Trop. Insect Sci. 1983, 4, 357–362. [Google Scholar] [CrossRef]
  25. Naim, F.; Nakasugi, K.; Crowhurst, R.N.; Hilario, E.; Zwart, A.B.; Hellens, R.P.; Taylor, J.M.; Waterhouse, P.M.; Wood, C.C. Advanced engineering of lipid metabolism in Nicotiana benthamiana using a draft genome and the V2 viral silencing-suppressor protein. PLoS ONE 2012, 7, e52717. [Google Scholar] [CrossRef]
  26. Aronesty, E. Comparison of sequencing utility programs. Open Bioinform. J. 2013, 7, 1–8. [Google Scholar] [CrossRef]
  27. Langmead, B.; Trapnell, C.; Pop, M.; Salzberg, S.L. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10, R25. [Google Scholar] [CrossRef] [Green Version]
  28. Leggett, R.M.; Clavijo, B.J.; Clissold, L.; Clark, M.D.; Caccamo, M. NextClip: An analysis and read preparation tool for Nextera Long Mate Pair libraries. Bioinformatics 2014, 30, 566–568. [Google Scholar] [CrossRef] [Green Version]
  29. Vurture, G.W.; Sedlazeck, F.J.; Nattestad, M.; Underwood, C.J.; Fang, H.; Gurtowski, J.; Schatz, M.C. GenomeScope: Fast reference-free genome profiling from short reads. Bioinformatics 2017, 33, 2202–2204. [Google Scholar] [CrossRef] [Green Version]
  30. Gnerre, S.; MacCallum, I.; Przybylski, D.; Ribeiro, F.J.; Burton, J.N.; Walker, B.J.; Sharpe, T.; Hall, G.; Shea, T.P.; Sykes, S.; et al. High-quality draft assemblies of mammalian genomes from massively parallel sequence data. Proc. Natl. Acad. Sci. USA 2011, 108, 1513–1518. [Google Scholar] [CrossRef] [Green Version]
  31. Kajitani, R.; Toshimoto, K.; Noguchi, H.; Toyoda, A.; Ogura, Y.; Okuno, M.; Yabana, M.; Harada, M.; Nagayasu, E.; Maruyama, H.; et al. Efficient de novo assembly of highly heterozygous genomes from whole-genome shotgun short reads. Genome Res. 2014, 24, 1384–1395. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Boetzer, M.; Henkel, C.V.; Jansen, H.J.; Butler, D.; Pirovano, W. Scaffolding pre-assembled contigs using SSPACE. Bioinformatics 2011, 27, 578–579. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Li, R.; Zhu, H.; Ruan, J.; Qian, W.; Fang, X.; Shi, Z.; Li, Y.; Li, S.; Shan, G.; Kristiansen, K.; et al. De novo assembly of human genomes with massively parallel short read sequencing. Genome Res. 2010, 20, 265–272. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Wood, D.E.; Salzberg, S.L. Kraken: Ultrafast metagenomic sequence classification using exact alignments. Genome Biol. 2014, 15, R46. [Google Scholar] [CrossRef] [Green Version]
  35. Wences, A.H.; Schatz, M.C. Metassembler: Merging and optimizing de novo genome assemblies. Genome Biol. 2015, 16, 207. [Google Scholar] [CrossRef] [Green Version]
  36. Waterhouse, R.M.; Seppey, M.; Simão, F.A.; Zdobnov, E.M. Using BUSCO to assess insect genomic resources. In Insect Genomics; Brown, S., Pfrender, M., Eds.; Methods in Molecular Biology; Humana Press: New York, NY, USA, 2019; Volume 1858, pp. 59–74. [Google Scholar]
  37. Corcoran, J.A.; Jordan, M.D.; Thrimawithana, A.H.; Crowhurst, R.N.; Newcomb, R.D. The peripheral olfactory repertoire of the lightbrown apple moth, Epiphyas postvittana. PLoS ONE 2015, 10, e0128596. [Google Scholar] [CrossRef] [Green Version]
  38. Grapputo, A.; Thrimawithana, A.H.; Steinwender, B.; Newcomb, R.D. Differential gene expression in the evolution of sex pheromone communication in New Zealand’s endemic leafroller moths of the genera Ctenopseustis and Planotortrix. BMC Genom. 2018, 19, 94. [Google Scholar] [CrossRef]
  39. Boratyn, G.M.; Thierry-Mieg, J.; Thierry-Mieg, D.; Busby, B.; Madden, T.L. Magic-BLAST, an accurate RNA-seq aligner for long and short reads. BMC Bioinform. 2019, 20, 405. [Google Scholar] [CrossRef]
  40. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [Green Version]
  41. Laetsch, D.R.; Blaxter, M.L. BlobTools: Interrogation of genome assemblies. F1000Research 2017, 6, 1287. [Google Scholar] [CrossRef]
  42. Stanke, M.; Keller, O.; Gunduz, I.; Hayes, A.; Waack, S.; Morgenstern, B. AUGUSTUS: Ab initio prediction of alternative transcripts. Nucleic Acids Res. 2006, 34, W435–W439. [Google Scholar] [CrossRef] [Green Version]
  43. Emms, D.M.; Kelly, S. OrthoFinder: Solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy. Genome Biol. 2015, 16, 157. [Google Scholar] [CrossRef] [Green Version]
  44. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [Green Version]
  45. Capella-Gutiérrez, S.; Silla-Martínez, J.M.; Gabaldón, T. trimAl: A tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 2009, 25, 1972–1973. [Google Scholar] [CrossRef]
  46. Nguyen, L.-T.; Schmidt, H.A.; von Haeseler, A.; Minh, B.Q. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef]
  47. De Bie, T.; Cristianini, N.; Demuth, J.P.; Hahn, M.W. CAFE: A computational tool for the study of gene family evolution. Bioinformatics 2006, 22, 1269–1271. [Google Scholar] [CrossRef] [Green Version]
  48. Kofler, R.; Orozco-terWengel, P.; Maio, N.D.; Pandey, R.V.; Nolte, V.; Futschik, A.; Kosiol, C.; Schlötterer, C. PoPoolation: A toolbox for population genetic analysis of next generation sequencing data from pooled individuals. PLoS ONE 2011, 6, e15925. [Google Scholar] [CrossRef] [Green Version]
  49. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [Green Version]
  50. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. The sequence alignment/map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [Green Version]
  51. Weir, B.S.; Cockerham, C.C. Estimating F-statistics for the analysis of population structure. Evolution 1984, 38, 1358–1370. [Google Scholar]
  52. Kofler, R.; Pandey, R.V.; Schlötterer, C. PoPoolation2: Identifying differentiation between populations using sequencing of pooled DNA samples (Pool-Seq). Bioinformatics 2011, 27, 3435–3436. [Google Scholar] [CrossRef] [Green Version]
  53. Pickrell, J.K.; Pritchard, J.K. Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet. 2012, 8, e1002967. [Google Scholar] [CrossRef] [Green Version]
  54. Fitak, R.R. OptM: Estimating the optimal number of migration edges on population trees using Treemix. Biol. Methods Protoc. 2021, 6, bpab017. [Google Scholar] [CrossRef]
  55. Milanesi, M.; Capomaccio, S.; Vajana, E.; Bomba, L.; Garcia, J.F.; Ajmone-Marsan, P.; Colli, L. BITE: An R package for biodiversity analyses. BioRxiv 2017, 181610. [Google Scholar] [CrossRef] [Green Version]
  56. You, M.; Yue, Z.; He, W.; Yang, X.; Yang, G.; Xie, M.; Zhan, D.; Baxter, S.W.; Vasseur, L.; Gurr, G.M.; et al. A heterozygous moth genome provides insights into herbivory and detoxification. Nat. Genet. 2013, 45, 220–225. [Google Scholar] [CrossRef] [Green Version]
  57. The International Silkworm Genome Consortium. The genome of a lepidopteran model insect, the silkworm Bombyx mori. Insect Biochem. Mol. Biol. 2008, 38, 1036–1045. [Google Scholar] [CrossRef]
  58. Harrison, M.C.; Jongepier, E.; Robertson, H.M.; Arning, N.; Bitard-Feildel, T.; Chao, H.; Childers, C.P.; Dinh, H.; Doddapaneni, H.; Dugan, S.; et al. Hemimetabolous genomes reveal molecular basis of termite eusociality. Nat. Ecol. Evol. 2018, 2, 557–566. [Google Scholar] [CrossRef] [Green Version]
  59. Faddeeva-Vakhrusheva, A.; Kraaijeveld, K.; Derks, M.F.L.; Anvar, S.Y.; Agamennone, V.; Suring, W.; Kampfraath, A.A.; Ellers, J.; Le Ngoc, G.; van Gestel, C.A.M.; et al. Coping with living in the soil: The genome of the parthenogenetic springtail Folsomia candida. BMC Genom. 2017, 18, 493. [Google Scholar] [CrossRef] [Green Version]
  60. Davey, J.W.; Chouteau, M.; Barker, S.L.; Maroja, L.; Baxter, S.W.; Simpson, F.; Merrill, R.M.; Joron, M.; Mallet, J.; Dasmahapatra, K.K.; et al. Major improvements to the Heliconius melpomene genome assembly used to confirm 10 chromosome fusion events in 6 million years of butterfly evolution. G3 Genes Genomes Genet. 2016, 6, 695–708. [Google Scholar] [CrossRef] [Green Version]
  61. Smith, N.M.A.; Wade, C.; Allsopp, M.H.; Harpur, B.A.; Zayed, A.; Rose, S.A.; Engelstädter, J.; Chapman, N.C.; Yagound, B.; Oldroyd, B.P. Strikingly high levels of heterozygosity despite 20 years of inbreeding in a clonal honey bee. J. Evol. Biol. 2019, 32, 144–152. [Google Scholar] [CrossRef]
  62. Dermauw, W.; Wybouw, N.; Rombauts, S.; Menten, B.; Vontas, J.; Grbić, M.; Clark, R.M.; Feyereisen, R.; Van Leeuwen, T. A link between host plant adaptation and pesticide resistance in the polyphagous spider mite Tetranychus urticae. Proc. Natl. Acad. Sci. USA 2013, 110, E113–E122. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Després, L.; David, J.-P.; Gallet, C. The evolutionary ecology of insect resistance to plant chemicals. Trends Ecol. Evol. 2007, 22, 298–307. [Google Scholar] [CrossRef] [PubMed]
  64. Yu, Q.; Lu, C.; Li, B.; Fang, S.; Zuo, W.; Dai, F.; Zhang, Z.; Xiang, Z. Identification, genomic organization and expression pattern of glutathione S-transferase in the silkworm, Bombyx mori. Insect Biochem. Mol. Biol. 2008, 38, 1158–1164. [Google Scholar] [CrossRef] [PubMed]
  65. Berenbaum, M.R.; Favret, C.; Schuler, M.A. On defining “Key Innovations” in an adaptive radiation: Cytochrome P450s and Papilionidae. Am. Nat. 1996, 148, S139–S155. [Google Scholar] [CrossRef]
  66. Cohen, M.B.; Schuler, M.A.; Berenbaum, M.R. A host-inducible cytochrome P-450 from a host-specific caterpillar: Molecular cloning and evolution. Proc. Natl. Acad. Sci. USA 1992, 89, 10920–10924. [Google Scholar] [CrossRef] [Green Version]
  67. Li, W.; Schuler, M.A.; Berenbaum, M.R. Diversification of furanocoumarin-metabolizing cytochrome P450 monooxygenases in two papilionids: Specificity and substrate encounter rate. Proc. Natl. Acad. Sci. USA 2003, 100, 14593–14598. [Google Scholar] [CrossRef] [Green Version]
  68. Scott, J.G.; Wen, Z. Cytochromes P450 of insects: The tip of the iceberg. Pest Manag. Sci. 2001, 57, 958–967. [Google Scholar] [CrossRef]
  69. Lien, N.T.K.; Ngoc, N.T.H.; Lan, N.N.; Hien, N.T.; Tung, N.V.; 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–1248. [Google Scholar] [CrossRef] [Green Version]
  70. Xu, N.; Sun, X.-H.; Liu, Z.-H.; Xu, Y.; Sun, Y.; Zhou, D.; Shen, B.; Zhu, C.-L. Identification and classification of differentially expressed genes in pyrethroid-resistant Culex pipiens pallens. Mol. Genet. Genom. 2019, 294, 861–873. [Google Scholar] [CrossRef]
  71. Enayati, A.A.; Ranson, H.; Hemingway, J. Insect glutathione transferases and insecticide resistance. Insect Mol. Biol. 2005, 14, 3–8. [Google Scholar] [CrossRef] [Green Version]
  72. Sheehan, D.; Meade, G.; Foley, V.M.; Dowd, C.A. Structure, function and evolution of glutathione transferases: Implications for classification of non-mammalian members of an ancient enzyme superfamily. Biochem. J. 2001, 360, 1–16. [Google Scholar] [CrossRef]
  73. Gawande, N.D.; Subashini, S.; Murugan, M.; Subbarayalu, M. Molecular screening of insecticides with sigma glutathione S-transferases (GST) in cotton aphid Aphis gossypii using docking. Bioinformation 2014, 10, 679–683. [Google Scholar] [CrossRef] [Green Version]
  74. Ahn, S.-J.; Vogel, H.; Heckel, D.G. Comparative analysis of the UDP-glycosyltransferase multigene family in insects. Insect Biochem. Mol. Biol. 2012, 42, 133–147. [Google Scholar] [CrossRef]
  75. Bock, K.W. Vertebrate UDP-glucuronosyltransferases: Functional and evolutionary aspects. Biochem. Pharmacol. 2003, 66, 691–696. [Google Scholar] [CrossRef]
  76. Li, X.; Schuler, M.A.; Berenbaum, M.R. Molecular mechanisms of metabolic resistance to synthetic and natural xenobiotics. Annu. Rev. Entomol. 2007, 52, 231–253. [Google Scholar] [CrossRef]
  77. Cui, X.; Wang, C.; Wang, X.; Li, G.; Liu, Z.; Wang, H.; Guo, X.; Xu, B. Molecular mechanism of the UDP-glucuronosyltransferase 2B20-like gene (AccUGT2B20-like) in pesticide resistance of Apis cerana cerana. Front. Genet. 2020, 11. [Google Scholar] [CrossRef]
  78. Li, X.; Shi, H.; Gao, X.; Liang, P. Characterization of UDP-glucuronosyltransferase genes and their possible roles in multi-insecticide resistance in Plutella xylostella (L.). Pest Manag. Sci. 2018, 74, 695–704. [Google Scholar] [CrossRef]
  79. Dunse, K.M.; Kaas, Q.; Guarino, R.F.; Barton, P.A.; Craik, D.J.; Anderson, M.A. Molecular basis for the resistance of an insect chymotrypsin to a potato type II proteinase inhibitor. Proc. Natl. Acad. Sci. USA 2010, 107, 15016–15021. [Google Scholar] [CrossRef] [Green Version]
  80. Souza, T.P.; Dias, R.O.; Castelhano, E.C.; Brandão, M.M.; Moura, D.S.; Silva-Filho, M.C. Comparative analysis of expression profiling of the trypsin and chymotrypsin genes from Lepidoptera species with different levels of sensitivity to soybean peptidase inhibitors. Comp. Biochem. Physiol. Part B Biochem. Mol. Biol. 2016, 196–197, 67–73. [Google Scholar] [CrossRef]
  81. Hou, M.-Z.; Shen, G.-M.; Wei, D.; Li, Y.-L.; Dou, W.; Wang, J.-J. Characterization of Bactrocera dorsalis serine proteases and evidence for their indirect role in insecticide tolerance. Int. J. Mol. Sci. 2014, 15, 3272–3286. [Google Scholar] [CrossRef] [Green Version]
  82. Dugdale, J.S. Lepidoptera—Annotated Catalogue, and Keys to Family-Group Taxa; Fauna of New Zealand; Manaaki Whenua Press: Lincoln, New Zealand, 1988; ISBN 978-0-477-02518-8. [Google Scholar]
  83. Zimmerman, E.C. Insects of Hawaii. Volume 9, Microlepidoptera; University of Hawai’i Press: Honolulu, HI, USA, 1978. [Google Scholar]
  84. Danthanarayana, W. The bionomics, distribution and host range of the light brown apple moth, Epiphyas postvittana (Walk.) (Tortricidae). Aust. J. Zool. 1975, 23, 419–437. [Google Scholar] [CrossRef]
  85. Svensson, I. Anmärkningsvärda fynd av småfjärilar (Microlepidoptera) i Sverige 2008, Remarkable records of Microlepidoptera in Sweden during 2008. Entomol. Tidskr. 2009, 130, 61–72. [Google Scholar]
  86. Varela, L.G.; Johnson, M.W.; Strand, L.; Wilen, C.A.; Pickel, C. Light brown apple moth’s arrival in California worries commodity groups. Calif. Agric. 2008, 62, 57–61. [Google Scholar] [CrossRef]
  87. Lozier, J.D.; Mills, N.J. Predicting the potential invasive range of light brown apple moth (Epiphyas postvittana) using biologically informed and correlative species distribution models. Biol. Invasions 2011, 13, 2409. [Google Scholar] [CrossRef]
Figure 1. GenomeScope k-mer profile plots for short-read PE libraries illustrating the distribution of 21-mers in (A) 169 bp insert size PE data and (B) 260 bp insert size PE data.
Figure 1. GenomeScope k-mer profile plots for short-read PE libraries illustrating the distribution of 21-mers in (A) 169 bp insert size PE data and (B) 260 bp insert size PE data.
Insects 13 00264 g001
Figure 2. A phylogenetic tree of 22 insect species with Folsomia candida as an outgroup, obtained using the aligned protein sequences of 207 single-copy orthologs present in all 23 species, constructed using the ML method implemented in IQ-TREE. Branch lengths represent substitutions per site. All nodes have a bootstrap support of 100%. Numbers at the nodes and after each species name represent the number of gene families with an increase/decrease in gene number, respectively, estimated by CAFE. On the right, each bar shows the total gene counts per species grouped in the following categories: single-copy genes in all 23 species and in all except one species, gene duplicated in all or except one species, lineage-specific ortholog genes (in Lepidoptera, Diptera, Coleoptera, Hymenoptera, Hemiptera, Blattodea and Entomobryomorpha), species-specific ortholog genes, genes with other ortholog relationships and unassigned genes.
Figure 2. A phylogenetic tree of 22 insect species with Folsomia candida as an outgroup, obtained using the aligned protein sequences of 207 single-copy orthologs present in all 23 species, constructed using the ML method implemented in IQ-TREE. Branch lengths represent substitutions per site. All nodes have a bootstrap support of 100%. Numbers at the nodes and after each species name represent the number of gene families with an increase/decrease in gene number, respectively, estimated by CAFE. On the right, each bar shows the total gene counts per species grouped in the following categories: single-copy genes in all 23 species and in all except one species, gene duplicated in all or except one species, lineage-specific ortholog genes (in Lepidoptera, Diptera, Coleoptera, Hymenoptera, Hemiptera, Blattodea and Entomobryomorpha), species-specific ortholog genes, genes with other ortholog relationships and unassigned genes.
Insects 13 00264 g002
Figure 3. Distribution of average Fst values among populations of Epiphyas postvittana calculated from the SNPs within gene models in the genome.
Figure 3. Distribution of average Fst values among populations of Epiphyas postvittana calculated from the SNPs within gene models in the genome.
Insects 13 00264 g003
Figure 4. Historical relationships and inferred migration events (indicated by the coloured arrows) among Epiphyas postvittana populations obtained with TreeMix.
Figure 4. Historical relationships and inferred migration events (indicated by the coloured arrows) among Epiphyas postvittana populations obtained with TreeMix.
Insects 13 00264 g004
Table 1. Summary statistics of the assembly, features and quality assessment of the Epiphyas postvittana genome in comparison with those of other lepidopteran species.
Table 1. Summary statistics of the assembly, features and quality assessment of the Epiphyas postvittana genome in comparison with those of other lepidopteran species.
SpeciesEpiphyas postvittanaCydia pomonella [17]Trichoplusia ni [14]Helicoverpa armigeraPlutella xylostella [56]Spodoptera litura [19]Bombyx mori [57]
Genome assembly:
Assembly size (Mb)598.1772.9333.0337.0394.1438.3431.7
Number of scaffolds14,077171719169971819359743,622
Max scaffold length (Mb)3.7034.608.926.153.944.4016.12
N50 scaffold size (kb)301.178915.454648.101000.40737.18915.403717.00
Number of contigs46,1802221788524,22814,35713,63688,842
N50 contig length (kb)26.2862.5140.018.349.468.415.5
Genomic_features:
Protein coding gene31,38917,18414,38417,08618,07115,31714,623
Repeats (%)33.842.916.714.634.031.843.6
GC (%)38.337.435.536.138.436.638.8
Mean CDS length (bp)1204.01461.01512.31330.91385.01565.21212.7
Mean exon length (bp)230.0256.9316.7309.2214.0235.8222.2
Mean intron length (bp)1254.01205.31603.31379.41224.01938.41084.1
Mean number of exons per gene5.205.687.046.206.476.645.45
Quality_assessment:
BUSCO % present * (complete)
Genome98.5 (97.9)98.5 (97.8)98.8 (98.4)98.5 (97.1)92.7 (89.9)99.0 (98.4)94.5 (92.1)
Protein (OGS)95.7 (92.0)93.8 (90.6)97.8 (97.3)98.4 (96.0)80.1 (75.0)99.0 (97.9)96.3 (91.5)
* Present = complete + fragmented.
Table 2. Repeat classes in the Epiphyas postvittana genome compared with those of other Lepidoptera.
Table 2. Repeat classes in the Epiphyas postvittana genome compared with those of other Lepidoptera.
SpeciesEpiphyas postvittanaCydia pomonellaTrichoplusia niHelicoverpa armigeraPlutella xylostellaSpodoptera lituraBombyx mori
Repeat TypesLength (Mb)P%Length (Mb)P%Length (Mb)P%Length (Mb)P%Length (Mb)P% Length (Mb)P% Length (Mb)P%
DNA elements14.072.3526.953.4921.926.580.220.457.501.908.41.9669.0710.44
LINE17.832.9868.518.866.662.000.130.2720.625.2346.2710.81101.2515.31
LTR4.590.7711.391.471.510.450.060.129.862.52.410.564.550.69
SINE6.121.0221.032.725.471.64--2.010.518.511.998.621.30
Simple repeat5.110.855.080.663.191.92----4.521.066.260.95
Other0.510.091.811.393.191.92--0.00005700.670.160.930.14
Unclassified153.7425.71--16.795.040.080.16111.2528.2373.8917.2654.558.25
Total201.9733.77341.5442.8755.5416.6849.0014.60151.2438.37136.1731.8245.2337.33
Table 3. Number of SNP types, percentage and distribution in genic and intergenic regions of the Epiphyas postvittana genome.
Table 3. Number of SNP types, percentage and distribution in genic and intergenic regions of the Epiphyas postvittana genome.
SNP TypeGenicIntergenicTotal
Biallelic10,803 (57.56)5986 (31.9)16,789 (89.46)
Triallelic1261 (6.72)585 (3.12)1846 (9.84)
Tetraallelic96 (0.51)36 (0.19)132 (0.70)
Total12,160 (64.79)6607 (35.21)18,767 (100.0)
Among brackets % calculated on the total number of SNPs.
Table 4. Epiphyas postvittana population genetic parameters per population.
Table 4. Epiphyas postvittana population genetic parameters per population.
CountryPopulationNpiThetaD
AUSGumeracha (GUM)70520.00120.0015−0.0325
Wagga Wagga (NSW)59980.00120.0014−0.0423
Tasmania (TAS)50060.00090.0010−0.0291
Dalkeith (WAA)58900.00090.0012−0.0352
USACalifornia (CA)79340.00090.0011−0.0282
NZAwatere (AWA)43250.00090.0010−0.0339
Clyde (CLY)40050.00070.0008−0.0259
Hawke’s Bay (HB)45230.00090.0011−0.0373
Motueka (MOT)56250.00090.0011−0.0299
Te Puke (TPU)36690.00080.0009−0.0422
Table 5. Average FST across the whole genome of Epiphyas postvittana (calculated as average of all gene and intergenic FSTs). For population abbreviation, see Table 4.
Table 5. Average FST across the whole genome of Epiphyas postvittana (calculated as average of all gene and intergenic FSTs). For population abbreviation, see Table 4.
GUMNSWTASWAACAAWACLYHBMOTTPU
GUM-0.106110.11990.140670.121270.108160.150750.119490.115230.14398
NSW -0.145480.13590.135670.124610.160920.122510.133780.14405
TAS -0.18850.160990.128620.181550.156950.132250.19199
WAA -0.171580.161460.19730.16350.17030.18416
CA -0.148440.192140.15660.150680.18032
AWA -0.145050.113110.098550.14918
CLY -0.126740.148360.15415
HB -0.131350.08647
MOT -0.16317
TPU -
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Thrimawithana, A.H.; Wu, C.; Christeller, J.T.; Simpson, R.M.; Hilario, E.; Tooman, L.K.; Begum, D.; Jordan, M.D.; Crowhurst, R.; Newcomb, R.D.; et al. The Genomics and Population Genomics of the Light Brown Apple Moth, Epiphyas postvittana, an Invasive Tortricid Pest of Horticulture. Insects 2022, 13, 264. https://doi.org/10.3390/insects13030264

AMA Style

Thrimawithana AH, Wu C, Christeller JT, Simpson RM, Hilario E, Tooman LK, Begum D, Jordan MD, Crowhurst R, Newcomb RD, et al. The Genomics and Population Genomics of the Light Brown Apple Moth, Epiphyas postvittana, an Invasive Tortricid Pest of Horticulture. Insects. 2022; 13(3):264. https://doi.org/10.3390/insects13030264

Chicago/Turabian Style

Thrimawithana, Amali H., Chen Wu, John T. Christeller, Robert M. Simpson, Elena Hilario, Leah K. Tooman, Doreen Begum, Melissa D. Jordan, Ross Crowhurst, Richard D. Newcomb, and et al. 2022. "The Genomics and Population Genomics of the Light Brown Apple Moth, Epiphyas postvittana, an Invasive Tortricid Pest of Horticulture" Insects 13, no. 3: 264. https://doi.org/10.3390/insects13030264

APA Style

Thrimawithana, A. H., Wu, C., Christeller, J. T., Simpson, R. M., Hilario, E., Tooman, L. K., Begum, D., Jordan, M. D., Crowhurst, R., Newcomb, R. D., & Grapputo, A. (2022). The Genomics and Population Genomics of the Light Brown Apple Moth, Epiphyas postvittana, an Invasive Tortricid Pest of Horticulture. Insects, 13(3), 264. https://doi.org/10.3390/insects13030264

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