Estimation of long terminal repeat element content in the Helicoverpa zea genome from high-throughput sequencing of bacterial artificial chromosome pools

Validation of reference-assisted assembly using existing and novel Heliothine genomes

Chloridea subflexa and Chloridea virescens are a pair of closely related noctuid species exhibiting pheromone-based sexual isolation and divergent host plant preferences. We produced a novel Illumina short read C. subflexa genome assembly and an improved C. virescens genome assembly, which offer opportunities to study the genomic basis for evolutionarily important traits in this lepidopteran family with few genomic resources. We then examined the feasibility of reference-assisted assembly, an approach that leverages existing high quality genomic resources for genome improvement in closely related taxa and applied it to our Heliothine genomes. Our work demonstrates that reference-assisted assembly has the potential to enhance contiguity and completeness of existing insect genomic resources with minimal additional laboratory costs. We conclude by discussing both the potential and pitfalls of reference-assisted assembly according to the intended downstream assembly application.

Novel genetic basis of resistance to Bt toxin Cry1Ac in Helicoverpa zea

Crops genetically engineered to produce insecticidal proteins from the bacterium Bacillus thuringiensis have advanced pest management, but their benefits are diminished when pests evolve resistance. Elucidating the genetic basis of pest resistance to Bacillus thuringiensis toxins can improve resistance monitoring, resistance management, and the design of new insecticides. Here, we investigated the genetic basis of resistance to Bacillus thuringiensis toxin Cry1Ac in the lepidopteran Helicoverpa zea, one of the most damaging crop pests in the United States. To facilitate this research, we built the first chromosome-level genome assembly for this species, which has 31 chromosomes containing 375 Mb and 15,482 predicted proteins. Using a genome-wide association study, fine-scale mapping, and RNA-seq, we identified a 250-kb quantitative trait locus on chromosome 13 that was strongly associated with resistance in a strain of Helicoverpa zea that had been selected for resistance in the field and lab. The mutation in this quantitative trait locus contributed to but was not sufficient for resistance, which implies alleles in more than one gene contributed to resistance. This quantitative trait locus contains no genes with a previously reported role in resistance or susceptibility to Bacillus thuringiensis toxins. However, in resistant insects, this quantitative trait locus has a premature stop codon in a kinesin gene, which is a primary candidate as a mutation contributing to resistance. We found no changes in gene sequence or expression consistently associated with resistance for 11 genes previously implicated in lepidopteran resistance to Cry1Ac. Thus, the results reveal a novel and polygenic basis of resistance.

Genome Size Reversely Correlates With Host Plant Range in Helicoverpa Species

In organisms with very low percentages of transposable elements (TEs), genome size may positively or negatively correlate with host range, depending on whether host adaptation or host modification is the main route to host generalism. To test if this holds true for insect herbivores with greater percentages of TEs, we conducted flow cytometry to measure the endopolyploidy levels and C-values of the host modification (salivary gland and mandibular gland in head), host adaptation (midgut), and host use-independent tissues (male gonad, hemolymph, and Malpighian tubules) of the generalist Helicoverpa armigera and the head of its older specialist sister H. assulta. Larval salivary gland displayed a consecutive chain of endopolyploidy particles from 8Cx to higher than 32Cx and larval head and midgut had endopolyploidy nuclei clusters of 16Cx and 32Cx, whereas larval male gonad, hemolymph, and Malpighian tubules possessed no endopolyploidy nuclei of higher than 8Cx. The estimated genome size of the Solanaceae plant specialist H. assulta is 430 Mb, significantly larger than that of its older generalist sister Heliothis virescens (408 Mb) and those of its two generalist descendants H. armigera (394 Mb) and H. zea (363 Mb). These data not only reveal a negative correlation between host plant range and genome size in this terminal lineage, but also imply that Helicoverpa species appear to depend more on host modification than on host adaptation to achieve polyphagy.

Genomic innovations, transcriptional plasticity and gene loss underlying the evolution and divergence of two highly polyphagous and invasive Helicoverpa pest species

BACKGROUND

Helicoverpa armigera and Helicoverpa zea are major caterpillar pests of Old and New World agriculture, respectively. Both, particularly H. armigera, are extremely polyphagous, and H. armigera has developed resistance to many insecticides. Here we use comparative genomics, transcriptomics and resequencing to elucidate the genetic basis for their properties as pests.

RESULTS

We find that, prior to their divergence about 1.5 Mya, the H. armigera/H. zea lineage had accumulated up to more than 100 more members of specific detoxification and digestion gene families and more than 100 extra gustatory receptor genes, compared to other lepidopterans with narrower host ranges. The two genomes remain very similar in gene content and order, but H. armigera is more polymorphic overall, and H. zea has lost several detoxification genes, as well as about 50 gustatory receptor genes. It also lacks certain genes and alleles conferring insecticide resistance found in H. armigera. Non-synonymous sites in the expanded gene families above are rapidly diverging, both between paralogues and between orthologues in the two species. Whole genome transcriptomic analyses of H. armigera larvae show widely divergent responses to different host plants, including responses among many of the duplicated detoxification and digestion genes.

CONCLUSIONS

The extreme polyphagy of the two heliothines is associated with extensive amplification and neofunctionalisation of genes involved in host finding and use, coupled with versatile transcriptional responses on different hosts. H. armigera's invasion of the Americas in recent years means that hybridisation could generate populations that are both locally adapted and insecticide resistant.