G-U base-paired hpRNA confers potent inhibition of small RNA function in plants.

MicroRNA (miRNA) target mimicry technologies, utilizing naturally occurring miRNA decoy molecules, represent a potent tool for analyzing miRNA function. In this study, we present a highly efficient small RNA (sRNA) target mimicry design based on G-U base-paired hairpin RNA (hpG:U), which allows for the simultaneous targeting of multiple sRNAs. The hpG:U constructs consistently generate high amounts of intact, polyadenylated stem-loop (SL) RNA outside the nuclei, in contrast to traditional hairpin RNA designs with canonical base pairing (hpWT), which were predominantly processed resulting in a loop. By incorporating a 460-bp G-U base-paired double-stranded stem and a 312-576 nt loop carrying multiple miRNA target mimicry sites (GUMIC), the hpG:U construct displayed effective repression of three Arabidopsis miRNAs, namely miR165/166, miR157, and miR160, both individually and in combination. Additionally, a GUMIC construct targeting a prominent cluster of siRNAs derived from cucumber mosaic virus (CMV) Y-satellite RNA (Y-Sat) effectively inhibited Y-Sat siRNA-directed silencing of the chlorophyll biosynthetic gene CHLI, thereby reducing the yellowing symptoms in infected Nicotiana plants. Therefore, the G-U base-paired hpRNA, characterized by differential processing compared to traditional hpRNA, acts as an efficient decoy for both miRNAs and siRNAs. This technology holds great potential for sRNA functional analysis and the management of sRNA-mediated diseases.

Historical museum samples reveal signals of selection and drift in response to changing insecticide use in an agricultural pest moth.

In response to environmental and human-imposed selective pressures, agroecosystem pests frequently undergo rapid evolution, with some species having a remarkable capacity to rapidly develop pesticide resistance. Temporal sampling of genomic data can comprehensively capture such adaptive changes over time, for example, by elucidating allele frequency shifts in pesticide resistance loci in response to different pesticides. Here, we leveraged museum specimens spanning over a century of collections to generate temporal contrasts between pre- and post-insecticide populations of an agricultural pest moth, Helicoverpa armigera. We used targeted exon sequencing of 254 samples collected across Australia from the pre-1950s (prior to insecticide introduction) to the 1990s, encompassing decades of changing insecticide use. Our sequencing approach focused on genes that are known to be involved in insecticide resistance, environmental sensation, and stress tolerance. We found an overall lack of spatial and temporal population structure change across Australia. In some decades (e.g., 1960s and 1970s), we found a moderate reduction of genetic diversity, implying stochasticity in evolutionary trajectories due to genetic drift. Temporal genome scans showed extensive evidence of selection following insecticide use, although the majority of selected variants were low impact, and alternating trajectories of allele frequency change were suggestive of potential antagonistic pleiotropy. Our results provide new insights into recent evolutionary responses in an agricultural pest and show how temporal contrasts using museum specimens can improve mechanistic understanding of rapid evolution.

Adaptive Introgression across Semipermeable Species Boundaries between Local Helicoverpa zea and Invasive Helicoverpa armigera Moths.

Hybridization between invasive and native species has raised global concern, given the dramatic increase in species range shifts and pest outbreaks due to anthropogenic dispersal. Nevertheless, secondary contact between sister lineages of local and invasive species provides a natural laboratory to understand the factors that determine introgression and the maintenance or loss of species barriers. Here, we characterize the early evolutionary outcomes following secondary contact between invasive Helicoverpa armigera and native H. zea in Brazil. We carried out whole-genome resequencing of Helicoverpa moths from Brazil in two temporal samples: during the outbreak of H. armigera in 2013 and 2017. There is evidence for a burst of hybridization and widespread introgression from local H. zea into invasive H. armigera coinciding with H. armigera expansion in 2013. However, in H. armigera, the admixture proportion and the length of introgressed blocks were significantly reduced between 2013 and 2017, suggesting selection against admixture. In contrast to the genome-wide pattern, there was striking evidence for adaptive introgression of a single region from the invasive H. armigera into local H. zea, including an insecticide resistance allele that increased in frequency over time. In summary, despite extensive gene flow after secondary contact, the species boundaries are largely maintained except for the single introgressed region containing the insecticide-resistant locus. We document the worst-case scenario for an invasive species, in which there are now two pest species instead of one, and the native species has acquired resistance to pyrethroid insecticides through introgression.

Hybridization and gene flow in the mega-pest lineage of moth, Helicoverpa.

Within the mega-pest lineage of heliothine moths are a number of polyphagous, highly mobile species for which the exchange of adaptive traits through hybridization would affect their properties as pests. The recent invasion of South America by one of the most significant agricultural pests, Helicoverpa armigera, raises concerns for the formation of novel combinations of adaptive genes following hybridization with the closely related Helicoverpa zea To investigate the propensity for hybridization within the genus Helicoverpa, we carried out whole-genome resequencing of samples from six species, focusing in particular upon H. armigera population structure and its relationship with H. zea We show that both H. armigera subspecies have greater genetic diversity and effective population sizes than do the other species. We find no signals for gene flow among the six species, other than between H. armigera and H. zea, with nine Brazilian individuals proving to be hybrids of those two species. Eight had largely H. armigera genomes with some introgressed DNA from H. zea scattered throughout. The ninth resembled an F1 hybrid but with stretches of homozygosity for each parental species that reflect previous hybridization. Regions homozygous for H. armigera-derived DNA in this individual included one containing a gustatory receptor and esterase genes previously associated with host range, while another encoded a cytochrome P450 that confers insecticide resistance. Our data point toward the emergence of novel hybrid ecotypes and highlight the importance of monitoring H. armigera genotypes as they spread through the Americas.

Genomic insights into a population of introduced European rabbits Oryctolagus cuniculus in Australia and the development of genetic resistance to rabbit hemorrhagic disease virus.

The European rabbit (Oryctolagus cuniculus) is one of the most devastating invasive species in Australia. Since the 1950s, myxoma virus (MYXV) and rabbit haemorrhagic disease virus (RHDV) have been used to manage overabundant rabbit populations. Resistance to MYXV was observed within a few years of the release. More recently, resistance to lethal RHDV infection has also been reported, undermining the efficiency of landscape-scale rabbit control. Previous studies suggest that genetic resistance to lethal RHDV infection may differ locally between populations, yet the mechanisms of genetic resistance remain poorly understood. Here, we used genotyping by sequencing (GBS) data representing a reduced representation of the genome, to investigate Australian rabbit populations. Our aims were to understand the relationship between populations and identify possible genomic signatures of selection for RHDV resistance. One population we investigated had previously been reported to show levels of resistance to lethal RHDV infection. This population was compared to three other populations with lower or no previously reported RHDV resistance. We identified a set of novel candidate genes that could be involved in host-pathogen interactions such as virus binding and infection processes. These genes did not overlap with previous studies on RHDV resistance carried out in different rabbit populations, suggesting that multiple mechanisms are feasible. These findings provide useful insights into the different potential mechanisms of genetic resistance to RHDV virus which will inform future functional studies in this area.

Global population genomic signature of Spodoptera frugiperda (fall armyworm) supports complex introduction events across the Old World.

Native to the Americas, the invasive Spodoptera frugiperda (fall armyworm; FAW) was reported in West Africa in 2016, followed by its chronological detection across the Old World and the hypothesis of an eastward Asia expansion. We explored population genomic signatures of American and Old World FAW and identified 12 maternal mitochondrial DNA genome lineages across the invasive range. 870 high-quality nuclear single nucleotide polymorphic DNA markers identified five distinct New World population clusters, broadly reflecting FAW native geographical ranges and the absence of host-plant preferences. We identified unique admixed Old World populations, and admixed and non-admixed Asian FAW individuals, all of which suggested multiple introductions underpinning the pest's global spread. Directional gene flow from the East into eastern Africa was also detected, in contrast to the west-to-east spread hypothesis. Our study demonstrated the potential of population genomic approaches via international partnership to address global emerging pest threats and biosecurity challenges.

Whole-genome sequencing to detect mutations associated with resistance to insecticides and Bt proteins in Spodoptera frugiperda.

The fall armyworm (FAW), Spodoptera frugiperda, is a major pest native to the Americas that has recently invaded the Old World. Point mutations in the target-site proteins acetylcholinesterase-1 (ace-1), voltage-gated sodium channel (VGSC) and ryanodine receptor (RyR) have been identified in S. frugiperda as major resistance mechanisms to organophosphate, pyrethroid and diamide insecticides respectively. Mutations in the adenosine triphosphate-binding cassette transporter C2 gene (ABCC2) have also been identified to confer resistance to Cry1F protein. In this study, we applied a whole-genome sequencing (WGS) approach to identify point mutations in the target-site genes in 150 FAW individuals collected from China, Malawi, Uganda and Brazil. This approach revealed three amino acid substitutions (A201S, G227A and F290V) of S. frugiperda ace-1, which are known to be associated with organophosphate resistance. The Brazilian population had all three ace-1 point mutations and the 227A allele (mean frequency = 0.54) was the most common. Populations from China, Malawi and Uganda harbored two of the three ace-1 point mutations (A201S and F290V) with the 290V allele (0.47-0.58) as the dominant allele. Point mutations in VGSC (T929I, L932F and L1014F) and RyR (I4790M and G4946E) were not detected in any of the 150 individuals. A novel 12-bp insertion mutation in exon 15 of the ABCC2 gene was identified in some of the Brazilian individuals but absent in the invasive populations. Our results not only demonstrate robustness of the WGS-based genomic approach for detection of resistance mutations, but also provide insights for improvement of resistance management tactics in S. frugiperda.

Oxidative stress delays development and alters gene expression in the agricultural pest moth, Helicoverpa armigera.

Stress is a widespread phenomenon that all organisms must endure. Common in nature is oxidative stress, which can interrupt cell homeostasis to cause cell damage and may be derived from respiration or from environmental exposure through diet. As a result of the routine exposure from respiration, many organisms can mitigate the effects of oxidative stress, but less is known about responses to oxidative stress from other sources. Helicoverpa armigera is a major agricultural pest moth that causes significant damage to crops worldwide. Here, we examined the effects of oxidative stress on H. armigera by chronically exposing individuals to paraquat-a free radical producer-and measuring changes in development (weight, developmental rate, lifespan), and gene expression. We found that oxidative stress strongly affected development in H. armigera, with stressed samples spending more time as caterpillars than control samples (>24 vs. ~15 days, respectively) and therefore living longer overall. We found 1,618 up- and 761 down-regulated genes, respectively, in stressed versus control samples. In the up-regulated gene set, was an over-representation of biological processes related to cuticle and chitin development, glycine metabolism, and oxidation-reduction. Oxidative stress clearly impacts physiology and biochemistry in H. armigera and the interesting finding of an extended lifespan in stressed individuals could demonstrate hormesis, the phenomenon whereby toxic compounds can actually be beneficial at low doses. Collectively, our findings provide new insights into physiological and gene expression responses to oxidative stress in invertebrates.

Mitochondrial DNA genomes of five major Helicoverpa pest species from the Old and New Worlds (Lepidoptera: Noctuidae).

Five species of noctuid moths, Helicoverpa armigera, H. punctigera, H. assulta, H. zea, and H. gelotopoeon, are major agricultural pests inhabiting various and often overlapping global distributions. Visual identification of these species requires a great deal of expertise and misidentification can have repercussions for pest management and agricultural biosecurity. Here, we report on the complete mitochondrial genomes of H. assulta assulta and H. assulta afra, H. gelotopoeon, H. punctigera, H. zea, and H. armigera armigera and H. armigera conferta' assembled from high-throughput sequencing data. This study significantly increases the mitogenome resources for these five agricultural pests with sequences assembled from across different continents, including an H. armigera individual collected from an invasive population in Brazil. We infer the phylogenetic relationships of these five Helicoverpa species based on the 13 mitochondrial DNA protein-coding genes (PCG's) and show that two publicly available mitogenomes of H. assulta (KP015198 and KR149448) have been misidentified or incorrectly assembled. We further consolidate existing PCR-RFLP methods to cover all five Helicoverpa pest species, providing an updated method that will contribute to species differentiation and to future monitoring efforts of Helicoverpa pest species across different continents. We discuss the value of Helicoverpa mitogenomes to assist with species identification in view of the context of the rapid spread of H. armigera in the New World. With this work, we provide the molecular resources necessary for future studies of the evolutionary history and ecology of these species.

Detoxifying enzyme complements and host use phenotypes in 160 insect species.

We use the genomes of 160 insect species to test the hypothesis that the size of detoxifying enzyme families is greater in species using more chemically diverse food resources. Phylogenetically appropriate contrasts in subsamples of the data generally support the hypothesis. We find relatively high numbers of cytochrome P450, glutathione S-transferase and carboxyl/choline esterase genes in omnivores and herbivores feeding on chemically complex tissues and relatively low numbers of these genes in specialists on relatively simple diets, including plant sap, nectar and pollen, and blood. Among Lepidoptera feeding on green plant tissue and Condylognatha feeding on sap we also find more of these genes in highly polyphagous species, many of which are major agricultural pests. These genomic signatures of food resource use are consistent with the hypothesis that some taxa are preadapted for insecticide resistance evolution.

Protein-carbohydrate regulation in Helicoverpa amigera and H. punctigera and how diet protein-carbohydrate content affects insect susceptibility to Bt toxins.

Many animals, including insects, demonstrate a remarkable ability to regulate their intake of key macronutrients (e.g., soluble protein and digestible carbohydrates), which allows them to optimize fitness and performance. Additionally, regulating the intake of these two macronutrients enhances an animal's ability to defend itself against pathogens, mitigate the effects of secondary plant metabolites, and decrease susceptibility to toxins. In this study, we first compared how Bt-resistant and -susceptible lines of Helicoverpa armigera and Helicoverpa punctigera regulate their intake of protein (p) and digestible carbohydrates (c). We found that there was no difference in the self-selected protein-carbohydrate intake target between resistant and susceptible genotypes of either species. We then explored the extent to which food protein-carbohydrate content altered the susceptibility of these species to three Bt toxins: Cry1Ac, Cry2Ab, and Vip3Aa. We found that H. armigera on diets that had protein-carbohydrate profiles that matched their self-selected protein-carbohydrate intake target were significantly less susceptible to Cry1Ac. In contrast, diet protein-carbohydrate content did not affect H. punctigera susceptibility to Cry1Ac. For both H. armigera and H. punctigera, susceptibility to Cry2Ab and Vip3Aa toxins did not change as a function of diet protein-carbohydrate profile. These results, when combined with earlier work on H. zea, suggest food protein-carbohydrate content can modify susceptibility to some Bt toxins, but not others. An increased understanding of how the nutritional environment can modify susceptibility to different Bt toxins could help improve pest management and resistance management practices.

Are feeding preferences and insecticide resistance associated with the size of detoxifying enzyme families in insect herbivores?

The size of gene families associated with xenobiotic detoxification in insects may be associated with the complexity of their diets and their propensities to develop insecticide resistance. We test these hypotheses by collating the annotations of cytochrome P450, carboxyl/cholinesterase and glutathione S-transferase genes in 65 insect species with data on their host use and history of insecticide resistance. We find 2-4 fold variation across the species in the numbers of these genes and, in some orders, especially the Hymenoptera, there is a clear relationship between the numbers of genes and feeding preferences. However in other orders, in particular the Lepidoptera, no such relationship is apparent. The size of these three gene families also tend to correlate with insecticide resistance propensity but this may not be an independent effect because species with broader host ranges are more likely to be pests that are heavily sprayed with insecticides.

Complete mitochondrial genome of the soybean stem fly Melanagromyza sojae (Diptera: Agromyzidae).

We report the complete mitochondrial DNA genome of the soybean stem fly (SSF) Melanagromyza sojae from Brazil Santa Catarina state based on Illumina MiSeq sequence data. The estimated mitogenome is 15 475 base pairs (bp) (KT597923), with 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, and 22 tRNAs, and an estimated 579 bp AT-rich control region. Similar to other insects, the SSF mitogenome is A-T bias with 40.9% A, 36.7% T, 13.6% C, and 8.8% G. Molecular characterization of SSF mitogenome will facilitate the development of effective molecular markers, and robust and rapid identification of suspected biosecurity incursions and field infestations of this insect pest.