The genetic architecture of a complex ecological trait: host plant use in the specialist moth, Heliothis subflexa

Chromosome-level genome assembly and comparative genomics shed light on Helicoverpa assulta ecology and pest management

BACKGROUND

The Oriental tobacco budworm, Helicoverpa assulta, a specialist herbivorous insect that exclusively feeds on plants of the Solanaceae family, causes considerable damage to crops, such as tobacco and hot pepper. The absence of a genome sequence for this species hinders further research on its pest management and ecological adaptation.

RESULTS

Here, we present a high-quality chromosome-level genome of a Korean strain of H. assulta (Pyeongchang strain, K18). The total assembly spans 424.4 Mb with an N50 length of 14.54 Mb and 37% GC content. The assembled genome (ASM2961881v1) comprises 31 chromosomes, similar to other congeneric generalist species including H. armigera and H. zea. In terms of genomic assembly quality, the complete BUSCOs and repeat content accounted for 98.3% and 33.01% of the genome, respectively. Based on this assembly, 19 485 protein-coding genes were predicted in the genome annotation. A comparative analysis was conducted using the identified number of protein-coding genes in H. armigera (24154) and H. zea (23696). Out of the 19 485 predicted genes, 137 genes in 15 orthogroups were found to have expanded significantly in H. assulta, while 149 genes in 95 orthogroups contracted rapidly.

CONCLUSION

This study revealed specific gene expansions and contractions in H. assulta compared to those in its close relatives, indicating potential adaptations related to its specialized feeding habits. Also, the comparative genome analysis provides valuable insights for the integrated pest management of H. assulta and other globally significant pests in the Heliothinae subfamily. © 2024 Society of Chemical Industry.

Population Genomics for Insect Conservation

Insects constitute vital components of ecosystems. There is alarming evidence for global declines in insect species diversity, abundance, and biomass caused by anthropogenic drivers such as habitat degradation or loss, agricultural practices, climate change, and environmental pollution. This raises important concerns about human food security and ecosystem functionality and calls for more research to assess insect population trends and identify threatened species and the causes of declines to inform conservation strategies. Analysis of genetic diversity is a powerful tool to address these goals, but so far animal conservation genetics research has focused strongly on endangered vertebrates, devoting less attention to invertebrates, such as insects, that constitute most biodiversity. Insects' shorter generation times and larger population sizes likely necessitate different analytical methods and management strategies. The availability of high-quality reference genome assemblies enables population genomics to address several key issues. These include precise inference of past demographic fluctuations and recent declines, measurement of genetic load levels, delineation of evolutionarily significant units and cryptic species, and analysis of genetic adaptation to stressors. This enables identification of populations that are particularly vulnerable to future threats, considering their potential to adapt and evolve. We review the application of population genomics to insect conservation and the outlook for averting insect declines.

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.

Response to selection for parasitism of a suboptimal, low-preference host in an aphid parasitoid

Risks of postintroduction evolution in insects introduced to control invasive pests have been discussed for some time, but little is known about responses to selection or genetic architectures of host adaptation and thus about the likelihood or rapidity of evolutionary shifts. We report here results on the response to selection and genetic architecture of parasitism of a suboptimal, low-preference host species by an aphid parasitoid, Aphelinus rhamni, a candidate for introduction against the soy bean aphid, Aphis glycines. We selected A. rhamni for increased parasitism of Rhopalsiphum padi by rearing the parasitoid on this aphid for three generations. We measured parasitism of R. padi at generations 2 and 3, and at generation 3, we crossed and backcrossed parasitoids from the populations reared on R. padi with those from populations reared on Aphis glycines and compared parasitism of both R. padi and Aphis glycines among F 1 and backcross females. Aphelinus rhamni responded rapidly to selection for parasitism of R. padi. Selection for R. padi parasitism reduced parasitism of Aphis glycines, the original host of A. rhamni. However, parasitism of R. padi did not increase from generation 2 to generation 3 of selection, suggesting reduced variance available for selection, which was indeed found. We tested the associations between 184 single nucleotide polymorphisms (SNP) and increased parasitism of R. padi and found 28 SNP loci, some of which were associated with increased and others with decreased parasitism of R. padi. We assembled and annotated the A. rhamni genome, mapped all SNP loci to contigs and tested whether genes on contigs with SNP loci associated with parasitism were enriched for candidate genes or gene functions. We identified 80 genes on these contigs that mapped to 1.2 Mb of the 483 Mb genome of A. rhamni but found little enrichment of candidate genes or gene functions.

Long-Term Population Studies Uncover the Genome Structure and Genetic Basis of Xenobiotic and Host Plant Adaptation in the Herbivore Tetranychus urticae

Pesticide resistance arises rapidly in arthropod herbivores, as can host plant adaptation, and both are significant problems in agriculture. These traits have been challenging to study as both are often polygenic and many arthropods are genetically intractable. Here, we examined the genetic architecture of pesticide resistance and host plant adaptation in the two-spotted spider mite, Tetranychus urticae, a global agricultural pest. We show that the short generation time and high fecundity of T. urticae can be readily exploited in experimental evolution designs for high-resolution mapping of quantitative traits. As revealed by selection with spirodiclofen, an acetyl-CoA carboxylase inhibitor, in populations from a cross between a spirodiclofen-resistant and a spirodiclofen-susceptible strain, and which also differed in performance on tomato, we found that a limited number of loci could explain quantitative resistance to this compound. These were resolved to narrow genomic intervals, suggesting specific candidate genes, including acetyl-CoA carboxylase itself, clustered and copy variable cytochrome P450 genes, and NADPH cytochrome P450 reductase, which encodes a redox partner for cytochrome P450s. For performance on tomato, candidate genomic regions for response to selection were distinct from those responding to the synthetic compound and were consistent with a more polygenic architecture. In accomplishing this work, we exploited the continuous nature of allele frequency changes across experimental populations to resolve the existing fragmented T. urticae draft genome to pseudochromosomes. This improved assembly was indispensable for our analyses, as it will be for future research with this model herbivore that is exceptionally amenable to genetic studies.

The genetic architecture of ecological adaptation: intraspecific variation in host plant use by the lepidopteran crop pest Chloridea virescens

Intraspecific variation in ecologically important traits is a cornerstone of Darwin's theory of evolution by natural selection. The evolution and maintenance of this variation depends on genetic architecture, which in turn determines responses to natural selection. Some models suggest that traits with complex architectures are less likely to respond to selection than those with simple architectures, yet rapid divergence has been observed in such traits. The simultaneous evolutionary lability and genetic complexity of host plant use in the Lepidopteran subfamily Heliothinae suggest that architecture may not constrain ecological adaptation in this group. Here we investigate the response of Chloridea virescens, a generalist that feeds on diverse plant species, to selection for performance on a novel host, Physalis angulata (Solanaceae). P. angulata is the preferred host of Chloridea subflexa, a narrow specialist on the genus Physalis. In previous experiments, we found that the performance of C. subflexa on P. angulata depends on many loci of small effect distributed throughout the genome, but whether the same architecture would be involved in the generalist's adoption of P. angulata was unknown. Here we report a rapid response to selection in C. virescens for performance on P. angulata, and establish that the genetic architecture of intraspecific variation is quite similar to that of the interspecific differences in terms of the number, distribution, and effect sizes of the QTL involved. We discuss the impact of genetic architecture on the ability of Heliothine moths to respond to varying ecological selection pressures.

Contemporary evolution of a Lepidopteran species, Heliothis virescens, in response to modern agricultural practices

Adaptation to human-induced environmental change has the potential to profoundly influence the genomic architecture of affected species. This is particularly true in agricultural ecosystems, where anthropogenic selection pressure is strong. Heliothis virescens primarily feeds on cotton in its larval stages, and US populations have been declining since the widespread planting of transgenic cotton, which endogenously expresses proteins derived from Bacillus thuringiensis (Bt). No physiological adaptation to Bt toxin has been found in the field, so adaptation in this altered environment could involve (i) shifts in host plant selection mechanisms to avoid cotton, (ii) changes in detoxification mechanisms required for cotton-feeding vs. feeding on other hosts or (iii) loss of resistance to previously used management practices including insecticides. Here, we begin to address whether such changes occurred in H. virescens populations between 1997 and 2012, as Bt-cotton cultivation spread through the agricultural landscape. For our study, we produced an H. virescens genome assembly and used this in concert with a ddRAD-seq-enabled genome scan to identify loci with significant allele frequency changes over the 15-year period. Genetic changes at a previously described H. virescens insecticide target of selection were detectable in our genome scan and increased our confidence in this methodology. Additional loci were also detected as being under selection, and we quantified the selection strength required to elicit observed allele frequency changes at each locus. Potential contributions of genes near loci under selection to adaptive phenotypes in the H. virescens cotton system are discussed.

Application of a dense genetic map for assessment of genomic responses to selection and inbreeding in Heliothis virescens

Adaptation of pest species to laboratory conditions and selection for resistance to toxins in the laboratory are expected to cause inbreeding and genetic bottlenecks that reduce genetic variation. Heliothis virescens, a major cotton pest, has been colonized in the laboratory many times, and a few laboratory colonies have been selected for Bacillus thuringiensis (Bt) resistance. We developed 350-bp double-digest restriction-site associated DNA-sequencing (ddRAD-seq) molecular markers to examine and compare changes in genetic variation associated with laboratory adaptation, artificial selection and inbreeding in this nonmodel insect species. We found that allelic and nucleotide diversity declined dramatically in laboratory-reared H. virescens as compared with field-collected populations. The declines were primarily a result of the loss of low frequency alleles present in field-collected H. virescens. A further, albeit modest decline in genetic diversity was observed in a Bt-selected population. The greatest decline was seen in H. virescens that were sib-mated for 10 generations, in which more than 80% of loci were fixed for a single allele. To determine which regions of the genome were resistant to fixation in our sib-mated line, we generated a dense intraspecific linkage map containing three PCR-based and 659 ddRAD-seq markers. Markers that retained polymorphism were observed in small clusters spread over multiple linkage groups, but this clustering was not statistically significant. Overall, we have confirmed and extended the general expectations for reduced genetic diversity in laboratory colonies, provided tools for further genomic analyses and produced highly homozygous genomic DNA for future whole genome sequencing of H. virescens.

Factors influencing the effect size distribution of adaptive substitutions

The distribution of effect sizes of adaptive substitutions has been central to evolutionary biology since the modern synthesis. Early theory proposed that because large-effect mutations have negative pleiotropic consequences, only small-effect mutations contribute to adaptation. More recent theory suggested instead that large-effect mutations could be favoured when populations are far from their adaptive peak. Here we suggest that the distributions of effect sizes are expected to differ among study systems, reflecting the wide variation in evolutionary forces and ecological conditions experienced in nature. These include selection, mutation, genetic drift, gene flow, and other factors such as the degree of pleiotropy, the distance to the phenotypic optimum, whether the optimum is stable or moving, and whether new mutation or standing genetic variation provides the source of adaptive alleles. Our goal is to review how these factors might affect the distribution of effect sizes and to identify new research directions. Until more theory and empirical work is available, we feel that it is premature to make broad generalizations about the effect size distribution of adaptive substitutions important in nature.

Genetic architecture underlying host choice differentiation in the sympatric host races of Lochmaea capreae leaf beetles

Speciation in herbivorous insects has received considerable attention during the last few decades. Much of this group's diversity originates from adaptive population divergence onto different host plants, which often involves the evolution of specialized patterns of host choice behaviour. Differences in host choice often translates directly into divergence in mating sites, and therefore positive assortative mating will be created which will act as a strong barrier to gene flow. In this study, we first explored whether host choice is a genetically determined trait in the sympatric willow and birch host races of the leaf feeding beetle Lochmaea capreae, or whether larval experience influences adult host choice. Once we had established that host choice is a genetically based trait we determined its genetic architecture. To achieve this, we employed a reciprocal transplant design in which offspring from pure willow and birch cross-types, F1, F2 and backcrosses were raised on each host plant and their preference was determined upon reaching adulthood. We then applied joint-scaling analysis to uncover the genetic architecture of host preference. Our results suggest that rearing host does not have a pronounced effect on adult's host choice; rather the segregation pattern implies the existence of genetic loci affecting host choice in these host races. The joint-scaling analysis revealed that population differences in host choice are mainly influenced by the contribution of additive genetic effects and also maternally inherited cytoplasmic effects. We explore the implications of our findings for evolutionary dynamics of sympatric host race formation and speciation.