Next Article in Journal
Evaluation of a Non-Chemical Compared to a Non-Chemical Plus Silica Gel Approach to Bed Bug Management
Next Article in Special Issue
Insights from Population Genomics to Enhance and Sustain Biological Control of Insect Pests
Previous Article in Journal
Oral Ingestion of Bacterially Expressed dsRNA Can Silence Genes and Cause Mortality in a Highly Invasive, Tree-Killing Pest, the Emerald Ash Borer
Previous Article in Special Issue
Phylogeography of Organophosphate Resistant ace Alleles in Spanish Olive Fruit Fly Populations: A Mediterranean Perspective in the Global Change Context
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Diversity in Invasive Populations of Argentine Stem Weevil Associated with Adaptation to Biocontrol

1
Genomics Aotearoa and Department of Biochemistry, University of Otago, Dunedin 9054, Aotearoa, New Zealand
2
AgResearch, Grasslands Research Centre, Palmerston North 4410, New Zealand
3
AgResearch, Invermay Agricultural Centre, Private Bag 50034, Mosgiel, New Zealand
4
AgResearch, Lincoln Science Centre, Private Bag 4749, Christchurch 8140, New Zealand
5
Bio-Protection Research Centre, PO Box 85084, Lincoln University, Lincoln 7647, New Zealand
*
Author to whom correspondence should be addressed.
Insects 2020, 11(7), 441; https://doi.org/10.3390/insects11070441
Submission received: 5 June 2020 / Revised: 21 June 2020 / Accepted: 4 July 2020 / Published: 14 July 2020
(This article belongs to the Special Issue Population Genetics of Insects)

Abstract

:
Modified, agricultural landscapes are susceptible to damage by insect pests. Biological control of pests is typically successful once a control agent has established, but this depends on the agent’s capacity to co-evolve with the host. Theoretical studies have shown that different levels of genetic variation between the host and the control agent will lead to rapid evolution of resistance in the host. Although this has been reported in one instance, the underlying genetics have not been studied. To address this, we measured the genetic variation in New Zealand populations of the pasture pest, Argentine stem weevil (Listronotus bonariensis), which is controlled with declining effectiveness by a parasitoid wasp, Microctonus hyperodae. We constructed a draft reference genome of the weevil, collected samples from a geographical survey of 10 sites around New Zealand, and genotyped them using a modified genotyping-by-sequencing approach. New Zealand populations of Argentine stem weevil have high levels of heterozygosity and low population structure, consistent with a large effective population size and frequent gene flow. This implies that Argentine stem weevils were able to evolve more rapidly than their biocontrol agent, which reproduces asexually. These findings show that monitoring genetic diversity in biocontrol agents and their targets is critical for long-term success of biological control.

1. Introduction

Biological control of pests via the release of specialist natural predators can provide continued, self-sustaining, non-polluting and inexpensive management. Although the estimated chance of establishment of biocontrol programs is low [1], biocontrol agents maintain their efficacy once established [2], partly because the agent can evolve adaptations to counter adaptations in the host [3,4].
A biocontrol system has been in use since the 1990s to manage a destructive, invasive pest of New Zealand pastures, the Argentine stem weevil (ASW; Listronotus bonariensis Kuschel) (Coleoptera: Curculionidae). New Zealand pastures are highly modified, based on a very low number of introduced Palearctic plant species, and are particularity susceptible to pest impacts [5]. This susceptibility is due to low plant and animal diversity, resulting in low biotic resistance to invasive species [5]. In New Zealand, adult ASW populations can reach densities of 700 adults m-2 and cause economic impacts of up to NZ$200M per annum [6,7,8]. Conventional chemical control of ASW is ineffective, environmentally damaging and uneconomical (reviewed in [9,10]), because the stem-mining larvae avoid direct contact with the pesticides [9]. To complement endophyte-based plant resistance [11,12], the solitary wasp Microctonus hyperodae Loan (Hymenoptera: Braconidae) was released for biological control of ASW in 1992. Within three years of its release, parasitism of ASW by M. hyperodae had reached 90% [13], reducing or eliminating damage to pasture [13,14,15].
Although ASW was initially suppressed by M. hyperodae, this began to fail after about 14 generations [16,17,18]. Loss of efficacy may be the result of adaptation in weevil populations resulting from selection pressure by the parasitoid [17,18]. Because ASW reproduces sexually, ASW populations may have greater capacity to evolve than populations of M. hyperodae, which reproduces parthenogenetically. Empirical modelling of the ASW–M. hyperodae interaction has indicated that resistance is inevitable when hosts have more genetic variation than their predators [19]. Despite this theoretical pathway for resistance, other examples of evolution of resistance to classical biological control have not been reported [20].
Population-level studies of genetic variation in host and parasitoid are required to explain the evolution of resistance in this case. We address this with a genotyping-by-sequencing study of a geographical survey of 10 Argentine stem weevil populations collected from across the North and South Islands of New Zealand. Our experiments revealed a repetitive genome with high heterozygosity and a high proportion of unstructured variation across populations. This is consistent with large effective population size and gene flow between populations. Genetic variation was found along a latitudinal cline, and was associated with signatures of selection in regions of the genome, indicating a level of local adaptation within populations, but at the resolution of this study we found no evidence of genetic adaptation in parasitised weevils compared to parasitoid-free weevils. Our results showed that the amount of genetic variation in New Zealand populations of ASW is far greater than detected by traditional molecular markers [21,22], implying that ASW populations have evolved resistance via weak selection acting on variants of minor effect that existed before the introduction of M. hyperodae.

2. Materials and Methods

2.1. Weevil Sampling

We collected ASW samples from commercially-farmed ryegrass (Lolium perenne L.)/white clover (Trifolium repens L.) pastures, at 10 sites across New Zealand, using a suction device to collect ground litter (Table 1). Weevils were extracted from the litter in the laboratory. The locations sampled are illustrated in Figure 1. The map was plotted with the ggmap package for ggplot2 [23], using map tiles by Stamen Design under CC BY 3.0, with data by OpenStreetMap under ODbL.
For the comparison between parasitised and unparasitised weevils, samples were collected from ryegrass/clover pasture at Ruakura and Lincoln (as in Table 1) in August 2017. These samples were dissected as described by Goldson and Emberson [24] to determine whether they were parasitised. After dissection, heads were removed and used for genotyping.

2.2. Genome Assembly

To produce the short read dataset, an Illumina TruSeq PCR-free 350 bp insert library was generated from DNA extracted from a single, male Argentine stem weevil collected from endophyte-free hybrid ryegrass (L. perenne × Lolium multiflorum) at Lincoln, New Zealand. Library preparation and sequencing were performed by Macrogen Inc. (Seoul, Korea). A total of 158 Gb of 100 b and 150 b paired-end reads were generated from the TruSeq PCR-free library. After removing common sequencing contaminants and trimming adaptor sequences using BBTools [25], the short-read-only genome was assembled with Meraculous 2.2.6 [26,27,28]. Reproducible code for assembling the short-read dataset and assessing the assemblies is hosted at github.com/tomharrop/asw-nopcr.
To produce long reads from a single individual, we produced high molecular weight DNA from a single, male ASW collected from Ruakura, New Zealand, using a modified QIAGEN Genomic-tip 20/G extraction protocol [29]. We amplified the DNA using Φ29 multiple displacement amplification (QIAGEN REPLI-g Midi Kit) and debranched the amplified DNA using T7 Endonuclease I (New England Biolabs) according to the Oxford Nanopore Technologies Premium whole genome amplification protocol version WGA_kit9_v1. Debranching reduced the raw read N50 length to 9.0 kb. Amplified DNA was sequenced on 6 R9.4.1 flowcells using a MinION Mk1B sequencer (Oxford Nanopore Technologies). We also extracted high molecular weight DNA from three pools, each of 20 unsexed individuals collected from Ruakura, New Zealand. We sequenced this pooled DNA on 5 R9.4.1 flowcells, following the Genomic DNA by Ligation protocol (SQK-LSK109; Oxford Nanopore Technologies). We basecalled raw Nanopore data with Guppy 3.4.1 (Oxford Nanopore Technologies). We removed adaptor sequences from the long reads with Porechop 0.2.4 (github.com/rrwick/Porechop) and assembled with Flye 2.6 [30].
All genome assemblies were assessed by size and contiguity statistics and BUSCO analysis [31]. Redundant contigs were removed from the combined, long read assembly with Purge Haplotigs 0b9afdf [32] using a low, mid and high cutoff of 60, 120 and 190, respectively.
We were not able to estimate repeat content in the full genomes, because RepeatModeler 2.0.1 [33] identified >500 M High-scoring Segment Pairs (HSPs) and did not finish after running for 6 weeks with ~200 GB of physical RAM (results not shown). We estimated repeat content by subsetting the assemblies using the leave-one-out alignment method implemented in Funannotate clean 1.7.4 [34]. We then used RepeatModeler 2.0.1 and RepeatClassifier 2.0.1 [33] and RepeatMasker 4.1.0 [35] from the Dfam TE Tools Container v1.1 (github.com/Dfam-consortium/TETools) to estimate the repeat content of the subset assemblies. We identified less than 1 M HSPs in the subset assemblies, indicating that the repeat content of the subset assemblies is an underestimate of the repeat content in the full assemblies.
Reproducible code for assembling and assessing the long-read ASW genomes is hosted at github.com/TomHarrop/asw-flye-withpool.
We annotated the final, draft assembly with Funannotate 1.7.4 [34], using five RNA sequencing libraries generated from abdomens and heads of unparasitised adult ASW collected from Ruakura. Reproducible code for annotating the draft ASW genome is hosted at github.com/TomHarrop/asw-annotate.

2.3. Reduced-Representation Genome Sequencing, Processing and Analysis

DNA extraction and double digest RADseq (genotyping-by-sequencing, GBS; [36]) were performed by AgResearch, Invermay, New Zealand. DNA was extracted from individual weevil heads using the ZR-96 Tissue & Insect DNA Kit (Zymo Research, CA, USA). The DNA was digested with ApeKI and MspI and barcoded based on the Elshire method [37] with modifications [38]. Pooled libraries were size selected on a BluePippin (Sage Science, MA, USA) with a window size of 150–500 bp. 100 nt single-end reads were generated from libraries an Illumina HiSeq 2500 instrument.
We used a strict processing pipeline to prepare the raw GBS reads for locus assembly. Samples were demultiplexed with zero allowed barcode mismatches to 91–93 b reads, depending on barcode length. Reads were trimmed by searching for adaptors with a minimum match of 11 b. Reads shorter than 80 b after trimming were discarded. All remaining reads were truncated to 80 b to account for unmatched adaptor sequence < 11 b that may have been present at the end of reads. To remove overamplified samples, we calculated the GC content for each library and discarded samples with median read GC > 45%. We assembled loci against our draft genome using gstacks 2.53 [39].
For analysis, we used BCFtools 1.10 to remove sites with more than 2 alleles, minor allele frequency <0.05, or missing genotypes in more than 20% of individuals. After filtering loci, we also removed individuals that had missing genotypes at more than 20% of loci. We ran the Stacks 2.53 populations module [39] to calculate inbreeding (F) and heterozygosity statistics. We used PLINK 1.9 [40] to prune sites in linkage disequilibrium for principal components analysis and discriminant analysis of principal components with the adegenet 2.1.2 package for R [41,42], using the first four principal components. We used PGDSpider 2.1.1.5 [43] to convert the un-pruned dataset for detection of loci under selection with BayeScan 2.1 [44]. We analysed cross-population extended haplotype homozygosity with the R package rehh 3.1.0 [45]. For demographic analysis, we converted the pruned dataset to minor allele (folded) site frequency spectra using easySFS commit c2b26c5 from github.com/isaacovercast/easySFS. We estimated likelihood for each demographic model ten times using fastsimcoal2 2.6 [46] with 1 million simulations and 60 optimisation cycles per run. We compared model runs using delta likelihood (maximum observed likelihood - maximum estimated likelihood) and Akaike information criteria [47].
All the code we used to process the raw reads, assemble loci and run downstream analyses is hosted at github.com/TomHarrop/stacks-asw, including the parameters and software containers for each step.

2.4. Reproducibility and Data Availability

Raw sequence data for the ASW genome assembly and annotation and raw GBS reads are hosted at the National Center for Biotechnology Information Sequence Read Archive (NCBI SRA) under accession PRJNA640511. We used Snakemake [48] to arrange analysis steps into workflows and monitor dependencies, and Singularity [49] to capture the computing environment. Using the code repositories listed in each methods section, the final results can be reproduced from the raw data with a single command using Snakemake and Singularity. The source for this manuscript is hosted at github.com/TomHarrop/asw-gbs-genome-paper.

3. Results

3.1. The Argentine Stem Weevil Genome Is Repetitive and Polymorphic

To construct a reference for genotyping populations of Argentine stem weevils, we produced a draft assembly of the ASW genome. We initially attempted assembly from a single individual using PCR-free, short read sequencing. This resulted in a fragmented assembly with low BUSCO scores (Table 2). k-mer analysis on the raw short reads suggested 2.1 polymorphisms per 100 bp and a genomic repeat content of 28–48% in the individual we sequenced (Supplementary Materials Figure S1). We then attempted to produce a long-read genome assembly using whole-genome amplification (WGA) of high molecular weight (HMW) DNA from a single individual, followed by sequencing on the Oxford Nanopore Technologies (ONT) MinION sequencer. We produced 29.8 Gb of quality-filtered reads with an N50 length of 9.0 kb. Assembling the single individual, long read genome resulted in improved contiguity and BUSCO scores compared to the short-read assembly (Table 2). Consistent with the raw short read data, the single individual, long read genome was at least 70% repetitive (Table 2). To improve assembly across long repeats, we produced a second ONT dataset with longer reads from HMW DNA from three pools of 20 individuals each, without amplification. Sequencing these samples on the MinION sequencer produced a total of 12.0 Gb of quality-filtered reads with an N50 length of 19.5 kb. Assembling the longer reads generated from the pooled sample alone resulted in a more contiguous genome, but with lower BUSCO scores (Table 2). We constructed a combined, long-read genome using the pooled, long-read dataset for contig construction, and the single-individual, long-read dataset for assembly polishing. This improved the BUSCO scores, but produced a large number of redundant contigs (Table 2), presumably because of the high rate of heterozygosity in the pooled, long-read dataset. We then used the PCR-free, short read sequencing data from a single individual with the Purge Haplotigs pipeline to remove redundant contigs from the combined long read assembly [32]. This resulted in a final draft assembly of 1.1 Gb with an N50 length of 122.3 kb and a BUSCO completeness of 83.9%. The final draft assembly had a repeat content of at least 70% (Table 2; Supplementary Table S2), with a maximum repeat size of 30.4 kb and a repeat N50 length of 494 bp. The majority of the repeats were unclassified when compared against the Dfam 3.1 database [50], with 9.2% of the genome detected as retroelements and 7.5% as DNA transposons (Supplementary Materials Table S2) The non-repetitive regions had an N50 length of 1066 bp.

3.2. Genetic Variation Is Associated with Geography in NZ Populations of Argentine Stem Weevil

To measure genetic variation in invasive New Zealand populations of ASW, we collected individuals from 10 sites across the North and South Islands of New Zealand (Figure 1A). We genotyped 183 individuals with a modified genotyping-by-sequencing (GBS) protocol [37,38]. After strict trimming and filtering of the raw GBS data, we mapped reads from each individual against our draft genome and used gstacks to assemble loci [39]. For analysis, we removed loci with more than two alleles, minor allele frequency less than 0.05, or missing genotypes in more than 20% of individuals. We also removed individuals missing genotypes at more than 20% of loci. The filtered dataset comprised 7–15 individuals per location (total 116), genotyped at 52,051 biallelic SNPs. The mean observed heterozygosity ranged from 0.18–0.21 across populations (Figure 1B), and pairwise FST values between populations ranged from 0.024–0.051 (Figure 1C). For principal components analysis (PCA), we pruned the dataset to 18,715 biallelic SNPs that were not in linkage disequilibrium, using a correlation threshold of 0.1. PCA of genotypes at these sites revealed overlapping populations of ASW, with 9.2% of total variance explained by the first two components (Figure 1D). These populations of ASW are highly heterozygous, but the low proportion of total variance explained by the major principal components suggests that variation is not highly structured between populations. This is consistent with a large effective population size and gene flow between populations. To find variance between populations, we used discriminant analysis of principal components (DAPC) on the same set of pruned SNPs [41]. The major linear discriminant, which explains 96.7% of between-population variation, separates populations from North and South of the Main Divide (Figure 1E), although we found evidence of mixing in all populations except Lincoln (Figure 1F). Although the PCA suggests that the majority of the total variance is not structured, the DAPC indicates a degree of genetic isolation between populations from North and South of the Main Divide. This suggests that the Main Divide, which runs along the Southern Alps and divides the South Island, is the main geographic barrier to ASW populations in New Zealand.

3.3. Genetic Variation Is not Associated with Parasitism by M. hyperodae

To detect large-effect variants associated with susceptibility to parasitism by M. hyperodae, we genotyped weevils that had also been tested for the presence of a parasitoid larva. We used a total of 200 individuals, collected from Lincoln and Ruakura (Table 3), because of the extent of historical declines in parasitism rates recorded at these locations [18]. The weevils were examined for a parasitoid larva and genotyped at the same loci used for the geographical diversity survey. After filtering and pruning sites in linkage disequilibrium, we used 19,482 SNPs for PCA and DAPC in 95 parasitised inviduals and 84 individuals where a parasitoid was not detected (Table 3). We did not detect any genetic differentiation associated with the presence of a parasitoid, either within populations or between populations, or any evidence of skewed allele frequencies in these groups using BayeScan (lowest Q-value 0.97).

3.4. Genetic Differentiation between ASW Populations North and South of the Main Divide

Although we did not detect variation associated with presence of a parasitoid, parasitism rates vary across sites in NZ [18]. Regional genetic differences could be related to selection acting on different loci North and South of the Main Divide and/or genetic drift acting on isolated populations. To investigate the genetic differentiation between regions, we grouped individuals that were collected from North and South of the Main Divide (Figure 1). The two groups had a mean FST of 0.013. We detected 47 SNPs with skewed allele frequencies across 24 contigs in the draft genome with BayeScan (Figure 2). The contigs containing these SNPs had a total of 3–36 SNPs, and all 47 of the detected SNPs had positive α values, suggesting diversifying selection (Supplementary Materials Table S2). The SNPs identified by BayeScan were an average of 11.5 kb from the nearest genes. None of the closest genes were homologous to genes with characterized functions in insects. Using an additional method, 3 SNPs on another contig had outlying cross-population extended haplotype homozygosity (XPEHH) scores (Supplementary Materials Table S2; [51,52]). No common regions were identified by both methods.

3.5. Separate Incursions of ASW into New Zealand

The genetic differentiation between weevils from North and South of the Main Divide suggests the possibility of either multiple routes of entry, or incursion via a single route of entry followed by isolation and diversification. The level of heterozygosity we measured across populations also suggests that incursions were large and/or repeated. To test these different possibilities, we simulated site frequency spectra (SFS) under 10 different models of demographic history and compared them to the observed SFS. Our models covered single and multiple introductions, with either moderate or strong reduction in effective population size during introduction, and multiple introductions from different source populations (Figure 3A). The models that best matched the observed SFS had the North and South populations separated before bottlenecks (Figure 3A, models ii, iii and v), and support for these models was better when migration between populations was included. The model with the lowest mean delta likelihood supports separate routes of entry into New Zealand, with a bottleneck in each population prior to entry and migration between the North and South populations (Figure 3B).

4. Discussion

The purpose of this work was to investigate genetic variation in New Zealand populations of ASW and its possible relationship to resistance to M. hyperodae. Previous reports using randomly amplified polymorphic DNA (RAPD) markers and cytochrome C oxidase subunit I (COI) sequencing suggested a high degree of genetic similarity and identified a single COI haplotype in New Zealand populations [21,22]. In contrast, our results from a genome-wide genotype-by-sequencing (GBS) approach reveal a high level of genetic diversity within and between populations. We suggest that this standing variation provides an evolutionary advantage to ASW populations in comparison to their biocontrol agent, M. hyperodae. We expect variation to be limited by asexual reproduction in M. hyperodae (e.g., [19]). This lack of variation, and the inability of M. hyperodae to switch hosts [53], would limit the capacity of M. hyperodae to co-evolve with ASW. This indicates that genetic variation in both host and biocontrol agent need to be monitored with high-resolution genotyping to maintain success of biological control. More work will be required to describe the genetic mechanism of resistance and its prevalence in weevil populations, and to measure the amount of variation and population structure in M. hyperodae.
ASW was thought to have arrived in New Zealand in the early 20th century, probably via trade in pasture seeds or hay used for feed during stock transit [54]. The earlier reports of low genetic diversity, based on traditional molecular markers [21,22], suggested a limited incursion followed by dispersal and expansion. Our results provide three main pieces of evidence to update the proposed history of ASW incursions in New Zealand. The high heterozygosity across populations could be explained by a large initial incursion, repeated introductions, and/or an unusually high mutation rate. The genetic differentiation between populations from North and South of the Main Divide points to low migration rates between these regions. Our demographic modelling suggests that the populations expanded to their current effective sizes after already being separated into North and South populations. The most likely scenario is separate introductions from the same source population to North and South of the Main Divide, with some migration between the two populations. The power to resolve the possible evolutionary histories that led to the current population structure of New Zealand weevils was provided by the increased resolution of genome-wide genotyping.
Despite the increased resolution of GBS compared to traditional markers, we did not detect regions of the genome associated with parasitism by M. hyperodae. Possible reasons for this include one or more of the following: i. resistance to biocontrol may not be genetic; ii. resistance may be encoded by part of the genome not captured in our assembly; iii. microscopic detection of the parasitoid may not be a strong enough phenotype to separate resistant and susceptible individuals, because individuals without a detectable parasitoid are not necessarily resistant, e.g., if they had not been exposed to the parasitoid before collection from the field; or iv. resistance is encoded by multiple regions of small effect, which we were unable to detect in our study. In model organisms, adaptive evolution in response to selective agents acting within the survivability distribution of a population takes the form of polygenic responses on standing variation [55,56]. The highest reported parasitism rate of ASW by M. hyperodae is 90% [13], implying that some individuals in a population survive predation. In other words, selection by M. hyperodae acts within the survivability distribution of ASW populations. Because we detected a large amount of standing variation in our survey of ASW populations, which may encode phenotypic variation in parasitism survivability, we suggest that a polygenic response is the most probable scenario. The number of markers yielded by legacy genotyping-by-sequencing approaches provides low power to detect polygenic responses resulting from weak selection on standing variation. Higher-resolution, genome-wide association studies using whole-genome resequencing with more individuals and a stronger resistance phenotype may allow detection of regions of the genome associated with resistance of the weevils to biocontrol.
Two draft weevil (Coleoptera: Curculionidae) genomes constructed from short reads have been deposited in the NCBI Genome database. The coffee berry borer, Hypothenemus hampei, has a draft genome size of 163 MB [57], and the mountain pine beetle, Dendroctonus ponderosae, has a draft genome size of 202 MB in males and 213 MB in females [58]. Draft genomes that incorporate long reads have been deposited for the red palm weevil (Rhynchophorus ferrugineus; GCA_012979105.1) and the rice weevil (Sitophilus oryzae; GCF_002938485.1). These assemblies are 782 MB and 771 MB, respectively. Assemblies using long reads capture more of the genome, presumably because larger repeat regions can be assembled. Our ASW genome of 1.1 GB is larger than other available weevil genomes, and has a high proportion of repetitive sequences. The contiguity statistics and BUSCO scores indicate draft quality, and we expect gaps in the assembly at larger repeat regions that were not sufficiently covered by long reads. Our attempt at short-read assembly of the Argentine stem weevil genome was not effective because of the extreme repeat content. The heterozygosity in weevil populations and lack of an inbred, laboratory strain made pooling individuals for sequencing undesirable. This is highlighed by the number of multiple-copy genes in the combined, long read assembly. Our strategy to assemble the ASW genome included contig construction with the longest reads, followed by assembly polishing with long reads from a single individual, and then redundant contig removal with PCR-free short reads from another single individual. This allowed us to optimise the contiguity and completeness of the genome whilst managing the number of redundant contigs (Table 2).

5. Conclusions

Our results show that New Zealand populations of ASW have a large amount of heterozygosity, and we suggest that this allowed them to evolve resistance to their biological control agent. This highlights the need for monitoring biological control systems by genome-wide genotyping.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-4450/11/7/441/s1, Figure S1: k-mer analysis of the short read sequencing dataset at k = 31. A k-mer distribution in the short read dataset before and after normalisation. The peak at 84–216× coverage (highlighted in green) suggests a diploid genome with a haploid size of 967 Mb. The haploid peak (22–56× coverage) indicates a high number of polymorphic sites in the single, male individual used to generate this dataset. BBTools [25] estimated 2.1 polymorphisms per 100 bp in this dataset. B Cumulative percentage of all k-mers vs. frequency. 72% of all k-mers are at 216× or lower coverage, indicating that the remaining 28% of k-mers are repetitive. The earlier inflexion point at 48% suggests that k-mers in the region from 48–72% may also be from repetitive regions. Table S1: Repeat content of the subset draft assembly. We used RepeatModeler to detect repeats and the Dfam_3.1 Combined Database with RepeatClassifier to classify them [33,50], then used RepeatMasker to locate repeats in the assembly [35]. 1: Most repeats fragmented by insertions or deletions have been counted as one element. Table S2: Number of SNPs under selection using BayeScan [44] (Q < 0.01) or cross-population extended haplotype homozygosity (XPEHH) analysis [51,52] (−log10p > 4). α is BayeScan’s locus-specific component of FST coefficient [44]. Positive values suggest diversifying selection. Positive XPEHH scores suggest selection in the South group, and negative scores suggest selection in the North group.

Author Contributions

Conceptualization, T.W.R.H., J.M.E.J., S.L.G. and P.K.D.; Formal analysis, T.W.R.H., M.F.L.L., R.J., S.G. and R.L.A.; Funding acquisition, J.M.E.J., S.L.G. and P.K.D.; Investigation, T.W.R.H., S.E.T., S.N.I., T.v.S., H.H., J.S. and S.L.G.; Resources, S.E.T., J.M.E.J. and S.L.G.; Writing—original draft, T.W.R.H., M.F.L.L., S.L.G. and P.K.D.; Writing—review & editing, J.M.E.J. All authors have read and agreed to the published version of the manuscript.

Funding

This project (2) was funded by the New Zealand Bio-Protection Research Centre to SLG, supported by the Tertiary Education Commission of New Zealand, and by the New Zealand Ministry of Business, Innovation and Employment via its funding of the ‘Genomics for Production & Security in a Biological Economy’ programme (C10X1306) to AgResearch.

Acknowledgments

We thank the AgResearch Animal Genomics team (Invermay), Nathan Brandt (USDA), Shola Olaniyan (NZMPI) and Malvika Bana for their expert dissection of the ASW samples used in this study. We are grateful to Diane Barton and Colin Ferguson (AgResearch, Invermay) for help with weevil retrieval and identification, and Joseph Guhlin (University of Otago & Genomics Aotearoa) for help with annotation of the ASW genome. We also thank the participants, organisers and instructors of the Genomics Aotearoa ONT MinION workshop in Dunedin in April 2018, where we generated the Nanopore reads for the pooled weevil samples. The workshop was funded by Genomics Aotearoa, supported by the New Zealand Ministry of Business Innovation and employment.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gurr, G.M.; Barlow, N.D.; Memmott, J.; Wratten, S.D.; Greathead, D.J. A History of Methodological, Theoretical and Empirical Approaches to Biological Control. In Biological Control: Measures of Success; Gurr, G., Wratten, S., Eds.; Springer: Dordrecht, The Netherlands, 2000; pp. 3–37. [Google Scholar] [CrossRef]
  2. Holt, R.D.; Hochberg, M.E. When Is Biological Control Evolutionarily Stable (or Is It)? Ecology 1997, 78, 1673–1683. [Google Scholar] [CrossRef]
  3. Pastoret, P.-P. Biological control of vertebrate pests: The history of myxomatosis—An experiment in evolution. Vet. J. 2000, 159, 219. [Google Scholar] [CrossRef] [Green Version]
  4. Lively, C.M.; Dybdahl, M.F. Parasite adaptation to locally common host genotypes. Nature 2000, 405, 679–681. [Google Scholar] [CrossRef] [PubMed]
  5. Goldson, S.L.; Barker, G.M.; Chapman, H.M.; Popay, A.J.; Stewart, A.V.; Caradus, J.R.; Barratt, B.I.P. Severe Insect Pest Impacts on New Zealand Pasture: The Plight of an Ecological Outlier. J. Insect Sci 2020, 20, 1–17. [Google Scholar] [CrossRef] [PubMed]
  6. Barker, G.; Addison, P.J. Argentine stem weevil populations and damage in ryegrass swards of contrasting Acremonium infection. In Proceedings of the 6th Australasian Conference on Grassland Invertebrate Ecology; Prestidge, R.A., Ed.; AgResearch: Hamilton, New Zealand, 1993. [Google Scholar]
  7. Prestidge, R.A.; Barker, G.M.; Pottinger, R.P. The economic cost of Argentine stem weevil in pastures in New Zealand. In Proceedings of the 44th New Zealand Weed and Pest Control Conference, New Zealand, 8 January 1991; pp. 165–170. [Google Scholar]
  8. Ferguson, C.M.; Barratt, B.I.P.; Bell, N.; Goldson, S.L.; Hardwick, S.; Jackson, M.; Jackson, T.A.; Phillips, C.B.; Popay, A.J.; Rennie, G.; et al. Quantifying the economic cost of invertebrate pests to New Zealand’s pastoral industry. N. Z. J. Agric. Res. 2019, 62, 255–315. [Google Scholar] [CrossRef] [Green Version]
  9. Goldson, S.L.; McNeill, M.R.; Stufkens, M.W.; Proffitt, J.R.; Pottinger, R.P.; Farrell, J.A. Importation and quarantine of Microctonus hyperodae a South American parasitoid of Argentine stem weevil. In Proceedings of the Forty Third New Zealand Weed and Pest Control Conference, Palmerston North, New Zealand, 8 January 1990; pp. 334–338. [Google Scholar] [CrossRef] [Green Version]
  10. Barlow, N.D.; Goldson, S.L. Alien invertebrates in New Zealand. In Biological Invasions: Economic and Environmental Costs of Alien Plant, Animal, and Microbe Species; Pimentel, D., Ed.; CRC Press: Boca Raton, FL, USA, 2002; pp. 195–216. [Google Scholar] [CrossRef]
  11. Johnson, L.J.; de Bonth, A.C.M.; Briggs, L.R.; Caradus, J.R.; Finch, S.C.; Fleetwood, D.J.; Fletcher, L.R.; Hume, D.E.; Johnson, R.D.; Popay, A.J.; et al. The exploitation of epichloae endophytes for agricultural benefit. Fungal Divers. 2013, 60, 171–188. [Google Scholar] [CrossRef]
  12. Kauppinen, M.; Saikkonen, K.; Helander, M.; Pirttilä, A.M.; Wäli, P.R. Epichloë grass endophytes in sustainable agriculture. Nat. Plants 2016, 2, 1–7. [Google Scholar] [CrossRef]
  13. Barker, G.M.; Addison, P.J. Early Impact of Endoparasitoid Microctonus hyperodae (Hymenoptera: Braconidae) After Its Establishment in Listronotus bonariensis (Coleoptera: Curculionidae) Populations of Northern New Zealand Pastures. J. Econ. Entomol. 2006, 99, 273–287. [Google Scholar] [CrossRef]
  14. Goldson, S.L.; Barron, M.C.; Kean, J.M.; van Koten, C. Argentine stem weevil (Listronotus bonariensis, Coleoptera: Curculionidae) population dynamics in Canterbury, New Zealand dryland pasture. Bull. Entomol. Res. 2011, 101, 295–303. [Google Scholar] [CrossRef]
  15. Barker, G.M. Biology of the Introduced Biocontrol Agent Microctonus hyperodae (Hymenoptera: Braconidae) and Its Host Listronotus bonariensis (Coleoptera: Curculionidae) in Northern New Zealand. Environ. Entomol 2013, 42, 902–914. [Google Scholar] [CrossRef] [Green Version]
  16. Popay, A.J.; McNeill, M.R.; Goldson, S.L.; Ferguson, C.M. The current status of Argentine stem weevil (Listronotus bonariensis) as a pest in the North Island of New Zealand. N. Z. Plant. Prot. 2011, 64, 55–62. [Google Scholar] [CrossRef] [Green Version]
  17. Goldson, S.L.; Tomasetto, F. Apparent Acquired Resistance by a Weevil to Its Parasitoid Is Influenced by Host Plant. Front. Plant. Sci. 2016, 7, 1259. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Tomasetto, F.; Tylianakis, J.M.; Reale, M.; Wratten, S.; Goldson, S.L. Intensified agriculture favors evolved resistance to biological control. Proc. Natl. Acad. Sci. USA 2017, 201618416. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Casanovas, P.; Goldson, S.L.; Tylianakis, J.M. Asymmetry in reproduction strategies drives evolution of resistance in biological control systems. PLoS ONE 2018, 13, e0207610. [Google Scholar] [CrossRef]
  20. Pennisi, E. In a first, natural selection defeats a biocontrol insect. Science 2017, 356, 570. [Google Scholar] [CrossRef]
  21. Williams, C.L.; Goldson, S.L.; Baird, D.B.; Bullock, D.W. Geographical origin of an introduced insect pest, Listronotus bonariensis (Kuschel), determined by RAPD analysis. Heredity 1994, 72, 412–419. [Google Scholar] [CrossRef] [Green Version]
  22. Vink, C.J.; Kean, J.M. PCR gut analysis reveals that Tenuiphantes tenuis (Araneae: Linyphiidae) is a potentially significant predator of Argentine stem weevil, Listronotus bonariensis (Coleoptera: Curculionidae), in New Zealand pastures. N. Z. J. Zool. 2013, 40, 304–313. [Google Scholar] [CrossRef]
  23. Kahle, D.; Wickham, H. ggmap: Spatial Visualization with ggplot2. R J. 2013, 5, 144–161. [Google Scholar] [CrossRef] [Green Version]
  24. Goldson, S.L.; Emberson, R.M. Reproductive morphology of the Argentine stem weevil, Hyperodes bonariensis (Coleoptera: Curculionidae). N. Z. J. Zool. 1981, 8, 67–77. [Google Scholar] [CrossRef] [Green Version]
  25. Bushnell, B. BBMap: A Fast, Accurate, Splice-Aware Aligner. Available online: https://www.osti.gov/biblio/1241166-bbmap-fast-accurate-splice-aware-aligner (accessed on 6 March 2020).
  26. Chapman, J.A.; Ho, I.; Sunkara, S.; Luo, S.; Schroth, G.P.; Rokhsar, D.S. Meraculous: De Novo Genome Assembly with Short Paired-End Reads. PLoS ONE 2011, 6, e23501. [Google Scholar] [CrossRef]
  27. Chapman, J.A.; Ho, I.Y.; Goltsman, E.; Rokhsar, D.S. Meraculous2: Fast accurate short-read assembly of large polymorphic genomes. arXiv 2017, arXiv:1608.01031. Available online: http://arxiv.org/abs/1608.01031 (accessed on 27 June 2017).
  28. Goltsman, E.; Ho, I.; Rokhsar, D. Meraculous-2D: Haplotype-sensitive Assembly of Highly Heterozygous genomes. arXiv 2017, arXiv:1703.09852. Available online: http://arxiv.org/abs/1703.09852 (accessed on 27 June 2017).
  29. Harrop, T. HMW DNA Extraction for Insects. Available online: https://www.protocols.io/view/hmw-dna-extraction-for-insects-pnwdmfe/metrics (accessed on 16 June 2020).
  30. Kolmogorov, M.; Yuan, J.; Lin, Y.; Pevzner, P.A. Assembly of long, error-prone reads using repeat graphs. Nat. Biotechnol. 2019, 37, 540–546. [Google Scholar] [CrossRef]
  31. Simão, F.A.; Waterhouse, R.M.; Ioannidis, P.; Kriventseva, E.V.; Zdobnov, E.M. BUSCO: Assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 2015, 31, 3210–3212. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Roach, M.J.; Schmidt, S.A.; Borneman, A.R. Purge Haplotigs: Allelic contig reassignment for third-gen diploid genome assemblies. BMC Bioinform. 2018, 19, 460. [Google Scholar] [CrossRef] [PubMed]
  33. Smit, A.F.A.; Hubley, R. RepeatModeler Open-1.0. Available online: http://www.repeatmasker.org (accessed on 6 March 2020).
  34. Love, J.; Palmer, J.; Stajich, J.; Esser, T.; Kastman, E.; Bogema, D.; Winter, D. Nextgenusfs/Funannotate: Funannotate v1.7.4. Available online: https://zenodo.org/record/3679386#.Xwvdxud5tPY (accessed on 6 March 2020).
  35. Smit, A.F.A.; Hubley, R.; Green, P. RepeatMasker Open-4.0. Available online: http://www.repeatmasker.org (accessed on 6 March 2020).
  36. Peterson, B.K.; Weber, J.N.; Kay, E.H.; Fisher, H.S.; Hoekstra, H.E. Double Digest RADseq: An Inexpensive Method for De Novo SNP Discovery and Genotyping in Model and Non-Model Species. PLoS ONE 2012, 7, e37135. [Google Scholar] [CrossRef] [Green Version]
  37. Elshire, R.J.; Glaubitz, J.C.; Sun, Q.; Poland, J.A.; Kawamoto, K.; Buckler, E.S.; Mitchell, S.E. A Robust, Simple Genotyping-by-Sequencing (GBS) Approach for High Diversity Species. PLoS ONE 2011, 6, e19379. [Google Scholar] [CrossRef] [Green Version]
  38. Dodds, K.G.; McEwan, J.C.; Brauning, R.; Anderson, R.M.; van Stijn, T.C.; Kristjánsson, T.; Clarke, S.M. Construction of relatedness matrices using genotyping-by-sequencing data. BMC Genom. 2015, 16, 1047. [Google Scholar] [CrossRef] [Green Version]
  39. Catchen, J.; Hohenlohe, P.A.; Bassham, S.; Amores, A.; Cresko, W.A. Stacks: An analysis tool set for population genomics. Mol. Ecol. 2013, 22, 3124–3140. [Google Scholar] [CrossRef] [Green Version]
  40. Chang, C.C.; Chow, C.C.; Tellier, L.C.; Vattikuti, S.; Purcell, S.M.; Lee, J.J. Second-generation PLINK: Rising to the challenge of larger and richer datasets. Gigascience 2015, 4. [Google Scholar] [CrossRef]
  41. Jombart, T.; Devillard, S.; Balloux, F. Discriminant analysis of principal components: A new method for the analysis of genetically structured populations. BMC Genet. 2010, 11, 94. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. R Core Team R: A Language and Environment for Statistical Computing. Available online: http://www.R-project.org (accessed on 8 December 2016).
  43. Lischer, H.E.L.; Excoffier, L. PGDSpider: An automated data conversion tool for connecting population genetics and genomics programs. Bioinformatics 2012, 28, 298–299. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Foll, M.; Gaggiotti, O. A Genome-Scan Method to Identify Selected Loci Appropriate for Both Dominant and Codominant Markers: A Bayesian Perspective. Genetics 2008, 180, 977–993. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Gautier, M.; Klassmann, A.; Vitalis, R. rehh 2.0: A reimplementation of the R package rehh to detect positive selection from haplotype structure. Mol. Ecol. Resour. 2017, 17, 78–90. [Google Scholar] [CrossRef]
  46. Excoffier, L.; Dupanloup, I.; Huerta-Sánchez, E.; Sousa, V.C.; Foll, M. Robust Demographic Inference from Genomic and SNP Data. PLoS Genet. 2013, 9, e1003905. [Google Scholar] [CrossRef] [Green Version]
  47. Akaike, H. A new look at the statistical model identification. IEEE Trans. Autom. Control. 1974, 19, 716–723. [Google Scholar] [CrossRef]
  48. Köster, J.; Rahmann, S. Snakemake—A scalable bioinformatics workflow engine. Bioinformatics 2012, 28, 2520–2522. [Google Scholar] [CrossRef] [Green Version]
  49. Kurtzer, G.M.; Sochat, V.; Bauer, M.W. Singularity: Scientific containers for mobility of compute. PLoS ONE 2017, 12, e0177459. [Google Scholar] [CrossRef]
  50. Hubley, R.; Finn, R.D.; Clements, J.; Eddy, S.R.; Jones, T.A.; Bao, W.; Smit, A.F.A.; Wheeler, T.J. The Dfam database of repetitive DNA families. Nucleic Acids Res. 2016, 44, D81–D89. [Google Scholar] [CrossRef] [Green Version]
  51. Sabeti, P.C.; Varilly, P.; Fry, B.; Lohmueller, J.; Hostetter, E.; Cotsapas, C.; Xie, X.; Byrne, E.H.; McCarroll, S.A.; Gaudet, R.; et al. Genome-wide detection and characterization of positive selection in human populations. Nature 2007, 449, 913–918. [Google Scholar] [CrossRef]
  52. Gautier, M.; Vitalis, R. rehh: An R package to detect footprints of selection in genome-wide SNP data from haplotype structure. Bioinformatics 2012, 28, 1176–1177. [Google Scholar] [CrossRef]
  53. Goldson, S.L.; McNeill, M.R.; Phillips, C.B.; Proffitt, J.R. Host specificity testing and suitability of the parasitoid Microctonus hyperodae (Hym.: Braconidae, Euphorinae) as a biological control agent of Listronotus bonariensis (Col.: Curculionidae) in New Zealand. Entomophaga 1992, 37, 483–498. [Google Scholar] [CrossRef]
  54. Brooking, T.; Pawson, E. Silences of Grass: Retrieving the Role of Pasture Plants in the Development of New Zealand and the British Empire. J. Imp. Commonw. Hist. 2007, 35, 417–435. [Google Scholar] [CrossRef]
  55. McKenzie, J.A.; Batterham, P. The genetic, molecular and phenotypic consequences of selection for insecticide resistance. Trends Ecol. Evol. 1994, 9, 166–169. [Google Scholar] [CrossRef]
  56. Green, L.; Battlay, P.; Fournier-Level, A.; Good, R.T.; Robin, C. Cis- and trans-acting variants contribute to survivorship in a naïve Drosophila melanogaster population exposed to ryanoid insecticides. Proc. Natl. Acad. Sci. USA 2019, 116, 10424–10429. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Vega, F.E.; Brown, S.M.; Chen, H.; Shen, E.; Nair, M.B.; Ceja-Navarro, J.A.; Brodie, E.L.; Infante, F.; Dowd, P.F.; Pain, A. Draft genome of the most devastating insect pest of coffee worldwide: The coffee berry borer, Hypothenemus hampei. Sci. Rep. 2015, 5, 12525. [Google Scholar] [CrossRef] [Green Version]
  58. Keeling, C.I.; Yuen, M.M.; Liao, N.Y.; Roderick Docking, T.; Chan, S.K.; Taylor, G.A.; Palmquist, D.L.; Jackman, S.D.; Nguyen, A.; Li, M.; et al. Draft genome of the mountain pine beetle, Dendroctonus ponderosae Hopkins, a major forest pest. Genome Biol. 2013, 14, R27. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Genetic diversity in NZ populations of Argentine stem weevil. (A) Weevil sampling locations. We collected Argentine stem weevils from 4 locations in the North Island and 6 locations in the South Island of New Zealand. Greymouth is in the South Island, but North of the Main Divide, which runs along the Southern Alps and partitions the South Island. The number of weevils genotyped (after filtering) from each location is shown on the map. Map tiles by Stamen Design under CC BY 3.0, with data by OpenStreetMap under ODbL. (B) Mean observed heterozygosity for each population. (C) Pairwise FST values between populations. (D) Principal components analysis (PCA) describing total variability and (E) discriminant analysis of principle components (DAPC) describing between-population variability of 116 individuals genotyped at 18,715 biallelic sites. In the PCA (D), populations overlap on the first two principal components (PC1 and PC2), but weevils sampled from higher latitudes have lower scores on PC1. PC1 and PC2 together explain 9.2% of variance in the dataset, indicating a high level of unstructured genetic variation in weevil populations. In the DAPC (E), the major linear discriminant (LD1) explains 96.7% of between-group variability. LD1 splits individuals from North and South of the Main Divide. LD2 separates Lincoln individuals from other individuals. (F) Posterior probability of group assignment for each individual. All populations contain individuals with high posterior probabilities of membership to other populations. This is consistent with gene flow between populations. We did not detect evidence of gene flow between populations from North and South of the Main Divide. Individuals sampled from Lincoln had the lowest posterior probabilities of membership to other populations.
Figure 1. Genetic diversity in NZ populations of Argentine stem weevil. (A) Weevil sampling locations. We collected Argentine stem weevils from 4 locations in the North Island and 6 locations in the South Island of New Zealand. Greymouth is in the South Island, but North of the Main Divide, which runs along the Southern Alps and partitions the South Island. The number of weevils genotyped (after filtering) from each location is shown on the map. Map tiles by Stamen Design under CC BY 3.0, with data by OpenStreetMap under ODbL. (B) Mean observed heterozygosity for each population. (C) Pairwise FST values between populations. (D) Principal components analysis (PCA) describing total variability and (E) discriminant analysis of principle components (DAPC) describing between-population variability of 116 individuals genotyped at 18,715 biallelic sites. In the PCA (D), populations overlap on the first two principal components (PC1 and PC2), but weevils sampled from higher latitudes have lower scores on PC1. PC1 and PC2 together explain 9.2% of variance in the dataset, indicating a high level of unstructured genetic variation in weevil populations. In the DAPC (E), the major linear discriminant (LD1) explains 96.7% of between-group variability. LD1 splits individuals from North and South of the Main Divide. LD2 separates Lincoln individuals from other individuals. (F) Posterior probability of group assignment for each individual. All populations contain individuals with high posterior probabilities of membership to other populations. This is consistent with gene flow between populations. We did not detect evidence of gene flow between populations from North and South of the Main Divide. Individuals sampled from Lincoln had the lowest posterior probabilities of membership to other populations.
Insects 11 00441 g001
Figure 2. Distribution of SNPs that are associated with altered allele frequencies between populations from North and South of the Main Divide over 24 contigs. 47 SNPs have altered allele frequencies, using the arbitrary Q-value cutoff of 0.01. α: BayeScan’s locus-specific component of FST coefficient [44]. Positive values suggest diversifying selection.
Figure 2. Distribution of SNPs that are associated with altered allele frequencies between populations from North and South of the Main Divide over 24 contigs. 47 SNPs have altered allele frequencies, using the arbitrary Q-value cutoff of 0.01. α: BayeScan’s locus-specific component of FST coefficient [44]. Positive values suggest diversifying selection.
Insects 11 00441 g002
Figure 3. Models used to simulate site frequency spectra (SFS) for comparison against the observed SFS. (A) Models i and ii represent introduction via a single route. In model i the population contracts on introduction and then begins expanding before dispersal. In model ii, dispersal occurs while the effective population size remains small. Models iii and iv represent multiple routes of entry, with or without reduction in population size compared to the current populations in New Zealand. Model v represents introduction to different sites from source populations that were isolated prior to introduction. All models were tested with and without migration between Northern and Southern populations. (B) Likelihood estimations from ten runs of each model. We used 1M simulations and 60 optimisation cycles per run [46]. Runs that ended with a non-finite delta likehood after 1M simulations are shown as an asterisk. The models that include a bottleneck after separation of the two populations had the lowest delta likelihoods, and model iii (with migration) had the lowest mean delta likelihood across runs.
Figure 3. Models used to simulate site frequency spectra (SFS) for comparison against the observed SFS. (A) Models i and ii represent introduction via a single route. In model i the population contracts on introduction and then begins expanding before dispersal. In model ii, dispersal occurs while the effective population size remains small. Models iii and iv represent multiple routes of entry, with or without reduction in population size compared to the current populations in New Zealand. Model v represents introduction to different sites from source populations that were isolated prior to introduction. All models were tested with and without migration between Northern and Southern populations. (B) Likelihood estimations from ten runs of each model. We used 1M simulations and 60 optimisation cycles per run [46]. Runs that ended with a non-finite delta likehood after 1M simulations are shown as an asterisk. The models that include a bottleneck after separation of the two populations had the lowest delta likelihoods, and model iii (with migration) had the lowest mean delta likelihood across runs.
Insects 11 00441 g003
Table 1. Weevil collection locations (see also Figure 1).
Table 1. Weevil collection locations (see also Figure 1).
LocationGPS Co-Ordinates (lat, lon)Date CollectedNumber Genotyped
Coromandel−37.20194, 175.59417June 201516
Ruakura−37.76750, 175.32361June 201516
Taranaki−39.61500, 174.30278July 201515
Wellington−41.13647, 175.35163July 201516
Greymouth−42.89506, 172.71926September 201616
Lincoln−43.64397, 172.44292July 201415
Ophir−45.10955, 169.58753August 201715
Mararoa Downs−45.50672, 167.97596May 201616
Mossburn−45.66966, 168.23884January 201616
Fortrose−46.57064, 168.79993November 201616
Table 2. Assembly statistics for the final draft genome and intermediate assemblies. 1. We estimated repeat content from a subset of contigs in the assembly (see Methods). n.d.: not determined.
Table 2. Assembly statistics for the final draft genome and intermediate assemblies. 1. We estimated repeat content from a subset of contigs in the assembly (see Methods). n.d.: not determined.
Short ReadSingle Individual, Long ReadPooled, Long ReadCombined, Long ReadFinal Draft
Assembly length (Gb)1.31.21.21.71.1
N5053,0464523295852812681
N50 length (kb)7.174.4112.686.4122.3
Gaps (%)3.50000
GC content (%)30.631.331.431.431.3
Complete, single-copy BUSCOs (%)32.772.27169.278.8
Complete, multiple-copy BUSCOs (%)17.27.55.917.45.1
Minimum 1 repeat content (%)n.d.7171.471.471.3
Table 3. Number of parasitised and unparasitised weevils from Ruakura and Lincoln.
Table 3. Number of parasitised and unparasitised weevils from Ruakura and Lincoln.
LocationParasitoidNumber GenotypedNumber after Filtering
RuakuraPresent5046
RuakuraNot detected5040
LincolnPresent5049
LincolnNot detected5044

Share and Cite

MDPI and ACS Style

Harrop, T.W.R.; Le Lec, M.F.; Jauregui, R.; Taylor, S.E.; Inwood, S.N.; van Stijn, T.; Henry, H.; Skelly, J.; Ganesh, S.; Ashby, R.L.; et al. Genetic Diversity in Invasive Populations of Argentine Stem Weevil Associated with Adaptation to Biocontrol. Insects 2020, 11, 441. https://doi.org/10.3390/insects11070441

AMA Style

Harrop TWR, Le Lec MF, Jauregui R, Taylor SE, Inwood SN, van Stijn T, Henry H, Skelly J, Ganesh S, Ashby RL, et al. Genetic Diversity in Invasive Populations of Argentine Stem Weevil Associated with Adaptation to Biocontrol. Insects. 2020; 11(7):441. https://doi.org/10.3390/insects11070441

Chicago/Turabian Style

Harrop, Thomas W. R., Marissa F. Le Lec, Ruy Jauregui, Shannon E. Taylor, Sarah N. Inwood, Tracey van Stijn, Hannah Henry, John Skelly, Siva Ganesh, Rachael L. Ashby, and et al. 2020. "Genetic Diversity in Invasive Populations of Argentine Stem Weevil Associated with Adaptation to Biocontrol" Insects 11, no. 7: 441. https://doi.org/10.3390/insects11070441

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