Whole genome sequence of the soybean aphid, Aphis glycines

Genome assembly and transcriptomic analysis to elucidate the ability of Nasonovia ribisnigri to break host plant resistance

Aphid genomic resources enable the study of complex life history traits and provide information on vector biology, host adaption and speciation. The currant-lettuce aphid (Nasonovia ribisnigri (Hemiptera: Aphididae) (Mosley)) is a cosmopolitan pest of outdoor lettuce (Lactuca sativa (Asterales: Asteraceae) (Linnaeus)). Until recently, the use of resistant cultivars was an effective method for managing N. ribisnigri. A resistant cultivar containing a single gene (Nr-locus), introduced in the 1980s, conferred complete resistance to feeding. Overreliance of this Nr-locus in lettuce resulted in N. ribisnigri's ability to break resistance mechanism, with first reports during 2003. Our work attempts to understand which candidate gene(s) are associated with this resistance-breaking mechanism. We present two de novo draft assembles for N. ribisnigri genomes, corresponding to both avirulent (Nr-locus susceptible) and virulent (Nr-locus resistant) biotypes. Changes in gene expression of the two N. ribisnigri biotypes were investigated using transcriptomic analyses of RNA-sequencing (RNA-seq) data to understand the potential mechanisms of resistance to the Nr-locus in lettuce. The draft genome assemblies were 94.2% and 91.4% complete for the avirulent and virulent biotypes, respectively. Out of the 18,872 differentially expressed genes, a single gene/locus was identified in N. ribisnigri that was shared between two resistant-breaking biotypes. This locus was further explored and validated in Real-Time Quantitative Reverse Transcription PCR (qRT-PCR) experiments and has predicted localisations in both the cytoplasm and nucleus. This is the first study to provide evidence that a single gene/locus is likely responsible for the ability of N. ribisnigri to overcome the Nr-locus resistance in the lettuce host.

Design, Synthesis, and Insecticidal Activity of Pyridino[1,2-a]pyrimidines Containing Indole Moeites at the 1-Position

Research on mesoionic structures in pesticide design has gained significant attention in recent years. However, the 1-position of pyridino[1,2-a]pyrimidine is usually designed with 2-chlorothiazole, 2-chloropyridine, or cyano moieties commonly found in neonicotinoid insecticides. In order to enrich the available pharmacophore library, here, we disclose a series of new pyridino[1,2-a]pyrimidine mesoionics bearing indole-containing substituents at the 1-position. Most of these target compounds are confirmed to have good insecticidal activity against aphids through bioevaluation. In addition, a three-dimensional structure-activity relationship model is established to allow access to optimal compound F45 with an LC50 value of 2.97 mg/L. This value is comparable to the property achieved by the positive control triflumezopyrim (LC50 = 2.94 mg/L). Proteomics and molecular docking analysis suggest that compound F45 has the potential to modulate the functioning of the aphid nervous system through its interaction with neuronal nicotinic acetylcholine receptors. This study expands the existing pharmacophore library for the future development of new mesoionic insecticides based on 1-position modifications of the pyridino[1,2-a]pyrimidine scaffold.

Elucidating the ecophysiology of soybean pod-sucking stinkbug Riptortus pedestris (Hemiptera: Alydidae) based on de novo genome assembly and transcriptome analysis

Food security is important for the ever-growing global population. Soybean, Glycine max (L.) Merr., is cultivated worldwide providing a key source of food, protein and oil. Hence, it is imperative to maintain or to increase its yield under different conditions including challenges caused by abiotic and biotic stresses. In recent years, the soybean pod-sucking stinkbug Riptortus pedestris has emerged as an important agricultural insect pest in East, South and Southeast Asia. Here, we present a genomics resource for R. pedestris including its genome assembly, messenger RNA (mRNA) and microRNA (miRNA) transcriptomes at different developmental stages and from different organs. As insect hormone biosynthesis genes (genes involved in metamorphosis) and their regulators such as miRNAs are potential targets for pest control, we analyzed the sesquiterpenoid (juvenile) and ecdysteroid (molting) hormone biosynthesis pathway genes including their miRNAs and relevant neuropeptides. Temporal gene expression changes of these insect hormone biosynthesis pathways were observed at different developmental stages. Similarly, a diet-specific response in gene expression was also observed in both head and salivary glands. Furthermore, we observed that microRNAs (bantam, miR-14, miR-316, and miR-263) of R. pedestris fed with different types of soybeans were differentially expressed in the salivary glands indicating a diet-specific response. Interestingly, the opposite arms of miR-281 (-5p and -3p), a miRNA involved in regulating development, were predicted to target Hmgs genes of R. pedestris and soybean, respectively. These observations among others highlight stinkbug's responses as a function of its interaction with soybean. In brief, the results of this study not only present salient findings that could be of potential use in pest management and mitigation but also provide an invaluable resource for R. pedestris as an insect model to facilitate studies on plant-pest interactions.

A high-quality genome assembly of the waterlily aphid Rhopalosiphum nymphaeae

Waterlily aphid, Rhopalosiphum nymphaeae (Linnaeus), is a host-alternating aphid known to feed on both terrestrial and aquatic hosts. It causes damage through direct herbivory and acting as a vector for plant viruses, impacting worldwide Prunus spp. fruits and aquatic plants. Interestingly, R. nymphaeae's ability to thrive in both aquatic and terrestrial conditions sets it apart from other aphids, offering a unique perspective on adaptation. We present the first high-quality R. nymphaeae genome assembly with a size of 324.4 Mb using PacBio long-read sequencing. The resulting assembly is highly contiguous with a contig N50 reached 12.7 Mb. The BUSCO evaluation suggested a 97.5% completeness. The R. nymphaeae genome consists of 16.9% repetitive elements and 16,834 predicted protein-coding genes. Phylogenetic analysis positioned R. nymphaeae within the Aphidini tribe, showing close relations to R. maidis and R. padi. The high-quality reference genome R. nymphaeae provides a unique resource for understanding genome evolution in aphids and paves the foundation for understanding host plant adaptation mechanisms and developing pest control strategies.

Untangling an insect's virome from its endogenous viral elements

BACKGROUND

Insects are an important reservoir of viral biodiversity, but the vast majority of viruses associated with insects have not been discovered. Recent studies have employed high-throughput RNA sequencing, which has led to rapid advances in our understanding of insect viral diversity. However, insect genomes frequently contain transcribed endogenous viral elements (EVEs) with significant homology to exogenous viruses, complicating the use of RNAseq for viral discovery.

METHODS

In this study, we used a multi-pronged sequencing approach to study the virome of an important agricultural pest and prolific vector of plant pathogens, the potato aphid Macrosiphum euphorbiae. We first used rRNA-depleted RNAseq to characterize the microbes found in individual insects. We then used PCR screening to measure the frequency of two heritable viruses in a local aphid population. Lastly, we generated a quality draft genome assembly for M. euphorbiae using Illumina-corrected Nanopore sequencing to identify transcriptionally active EVEs in the host genome.

RESULTS

We found reads from two insect-specific viruses (a Flavivirus and an Ambidensovirus) in our RNAseq data, as well as a parasitoid virus (Bracovirus), a plant pathogenic virus (Tombusvirus), and two phages (Acinetobacter and APSE). However, our genome assembly showed that part of the 'virome' of this insect can be attributed to EVEs in the host genome.

CONCLUSION

Our work shows that EVEs have led to the misidentification of aphid viruses from RNAseq data, and we argue that this is a widespread challenge for the study of viral diversity in insects.

Patterns of asexual reproduction of the soybean aphid, Aphis glycines (Matsumura), with and without the secondary symbionts Wolbachia and Arsenophonus, on susceptible and resistant soybean genotypes

Plant breeding is used to develop crops with host resistance to aphids, however, virulent biotypes often develop that overcome host resistance genes. We tested whether the symbionts, Arsenophonus (A) and Wolbachia (W), affect virulence and fecundity in soybean aphid biotypes Bt1 and Bt3 cultured on whole plants and detached leaves of three resistant, Rag1, Rag2 and Rag1 + 2, and one susceptible, W82, soybean genotypes. Whole plants and individual aphid experiments of A. glycines with and without Arsenophonus and Wolbachia did not show differences in overall fecundity. Differences were observed in peak fecundity, first day of deposition, and day of maximum nymph deposition of individual aphids on detached leaves. Bt3 had higher fecundity than Bt1 on detached leaves of all plant genotypes regardless of bacterial profile. Symbionts did not affect peak fecundity of Bt1 but increased it in Bt3 (A+W+) and all Bt3 strains began to deposit nymphs earlier than the Bt1 (A+W-). Arsenophonus in Bt1 delayed the first day of nymph deposition in comparison to aposymbiotic Bt1 except when reared on Rag1 + 2. For the Bt1 and Bt3 strains, symbionts did not result in a significant difference in the day they deposited the maximum number of nymphs nor was there a difference in survival or variability in number of nymphs deposited. Variability of number of aphids deposited was higher in aphids feeding on resistant plant genotypes. The impact of Arsenophonus on soybean aphid patterns of fecundity was dependent on the aphid biotype and plant genotype. Wolbachia alone had no detectable impact but may have contributed to the increased fecundity of Bt3 (A+W+). An individual based model, using data from the detached leaves experiment and with intraspecific competition removed, found patterns similar to those observed in the greenhouse and growth chamber experiments including a significant interaction between soybean genotype and aphid strain. Combining individual data with the individual based model of population growth isolated the impact of fecundity and host resistance from intraspecific competition and host health. Changes to patterns of fecundity, influenced by the composition and concentration of symbionts, may contribute to competitive interactions among aphid genotypes and influence selection on virulent aphid populations.

Hybridisation has shaped a recent radiation of grass-feeding aphids

BACKGROUND

Aphids are common crop pests. These insects reproduce by facultative parthenogenesis involving several rounds of clonal reproduction interspersed with an occasional sexual cycle. Furthermore, clonal aphids give birth to live young that are already pregnant. These qualities enable rapid population growth and have facilitated the colonisation of crops globally. In several cases, so-called "super clones" have come to dominate agricultural systems. However, the extent to which the sexual stage of the aphid life cycle has shaped global pest populations has remained unclear, as have the origins of successful lineages. Here, we used chromosome-scale genome assemblies to disentangle the evolution of two global pests of cereals-the English (Sitobion avenae) and Indian (Sitobion miscanthi) grain aphids.

RESULTS

Genome-wide divergence between S. avenae and S. miscanthi is low. Moreover, comparison of haplotype-resolved assemblies revealed that the S. miscanthi isolate used for genome sequencing is likely a hybrid, with one of its diploid genome copies closely related to S. avenae (~ 0.5% divergence) and the other substantially more divergent (> 1%). Population genomics analyses of UK and China grain aphids showed that S. avenae and S. miscanthi are part of a cryptic species complex with many highly differentiated lineages that predate the origins of agriculture. The complex consists of hybrid lineages that display a tangled history of hybridisation and genetic introgression.

CONCLUSIONS

Our analyses reveal that hybridisation has substantially contributed to grain aphid diversity, and hence, to the evolutionary potential of this important pest species. Furthermore, we propose that aphids are particularly well placed to exploit hybridisation events via the rapid propagation of live-born "frozen hybrids" via asexual reproduction, increasing the likelihood of hybrid lineage formation.

Chromosome-level genome assembly of the spotted alfalfa aphid Therioaphis trifolii

The spotted alfalfa aphid (SAA, Therioaphis trifolii) (Hemiptera: Aphididae) is a destructive pest of cultivated alfalfa (Medicago sativa L.) that leads to large financial losses in the livestock industry around the world. Here, we present a chromosome-scale genome assembly of T. trifolii, the first genome assembly for the aphid subfamily Calaphidinae. Using PacBio long-read sequencing, Illumina sequencing, and Hi-C scaffolding techniques, a 541.26 Mb genome was generated, with 90.01% of the assembly anchored into eight scaffolds, and the contig and scaffold N50 are 2.54 Mb and 44.77 Mb, respectively. BUSCO assessment showed a completeness score of 96.6%. A total of 13,684 protein-coding genes were predicted. The high-quality genome assembly of T. trifolii not only provides a genomic resource for the more complete analysis of aphid evolution, but also provides insights into the ecological adaptation and insecticide resistance of T. trifolii.

Tracking Adaptive Pathways of Invasive Insects: Novel Insight from Genomics

Despite the huge human and economic costs of invasive insects, which are the main group of invasive species, their environmental impacts through various mechanisms remain inadequately explained in databases and much of the invasion biology literature. High-throughput sequencing technology, especially whole-genome sequencing, has been used as a powerful method to study the mechanisms through which insects achieve invasion. In this study, we reviewed whole-genome sequencing-based advances in revealing several important invasion mechanisms of invasive insects, including (1) the rapid genetic variation and evolution of invasive populations, (2) invasion history and dispersal paths, (3) rapid adaptation to different host plant ranges, (4) strong environmental adaptation, (5) the development of insecticide resistance, and (6) the synergistic damage caused by invasive insects and endosymbiotic bacteria. We also discussed prevention and control technologies based on whole-genome sequencing and their prospects.

Functional Characterization of the Nuclear Receptor Gene SaE75 in the Grain Aphid, Sitobion avenae

Ecdysteroid hormones are key regulators of insect development and metamorphosis. Ecdysone-inducible E75, a major component of insect ecdysone signaling pathway, has been well characterized in holometabolous insects, however, barely in hemimetabolous species. In this study, a total of four full-length E75 cDNAs from the English grain aphid, Sitobion avenae, were identified, cloned, and characterized. The four SaE75 cDNAs contained 3048, 2625, 2505, and 2179 bp open reading frames (ORF), encoding 1015, 874, 856, and 835 amino acids, respectively. Temporal expression profiles showed that SaE75 expression was low in adult stages, while high in pseudo embryo and nymphal stages. SaE75 was differentially expressed between winged and wingless morphs. RNAi-mediated suppression of SaE75 led to substantial biological impacts, including mortality and molting defects. As for the pleiotropic effects on downstream ecdysone pathway genes, SaHr3 (hormone receptor like in 46) was significantly up-regulated, while Sabr-c (broad-complex core protein gene) and Saftz-f1 (transcription factor 1) were significantly down-regulated. These combined results not only shed light on the regulatory role of E75 in the ecdysone signaling pathway, but also provide a potential novel target for the long-term sustainable management of S. avenae, a devastating global grain pest.

NLRs are highly relevant resistance genes for aphid pests

Since the 20th century, when plant resistance to aphids was available, it has been widely used by farmers and the inheritance of plant resistance has been understood for several crops. However, it is only when the plant-aphid relationship was compared with that of microbial pathogens, that aphid resistance has begun to be understood and integrated into the plant immune network. Three of the four genes identified for plant resistance to aphid encode nucleotide-binding site leucine-rich repeat receptor (NLR) proteins responsible for aphid-effector triggered immunity, and NLRs are serious candidates for aphid resistance in four other plant species. Aphids are vectors for plant viruses, and aphid-effectors triggering immunity when they pierce plant cells are expected to trigger resistance to the viruses transmitted to the plant with effectors, as has been shown for aphid resistance in melon. This dual phenotype increases the interest of NLRs in the control of aphids.

Genome-Wide Identification and Expression Pattern of Cytochrome P450 Genes in the Social Aphid Pseudoregma bambucicola

Cytochrome P450 monooxygenases (P450s) have a variety of functions, including involvement in the metabolism of exogenous substances and the synthesis and degradation of endogenous substances, which are important for the growth and development of insects. Pseudoregma bambucicola is a social aphid that produces genetically identical but morphologically and behaviorally distinct first-instar soldiers and normal nymphs within colonies. In this study, we identified 43 P450 genes based on P. bambucicola genome data. Phylogenetic analysis showed that these genes were classified into 4 clans, 13 families, and 23 subfamilies. The CYP3 and CYP4 clans had a somewhat decreased number of genes. In addition, differential gene expression analysis based on transcriptome data showed that several P450 genes, including CYP18A1, CYP4G332, and CYP4G333, showed higher expression levels in soldiers compared to normal nymphs and adult aphids. These genes may be candidates for causing epidermal hardening and developmental arrest in soldiers. This study provides valuable data and lays the foundation for the study of functions of P450 genes in the social aphid P. bambucicola.

Mobilome of the Rhus Gall Aphid Schlechtendalia chinensis Provides Insight into TE Insertion-Related Inactivation of Functional Genes

Transposable elements (TEs) comprise a considerable proportion of insect genomic DNA; how they contribute to genome structure and organization is still poorly understood. Here, we present an analysis of the TE repertoire in the chromosome-level genome assembly of Rhus gall aphid Schlechtendalia chinensis. The TE fractions are composed of at least 32 different superfamilies and many TEs from different families were transcriptionally active in the S. chinensis genome. Furthermore, different types of transposase-derived proteins were also found in the S. chinensis genome. We also provide insight into the TEs related insertional inactivation, and exogenization of TEs in functional genes. We considered that the presence of TE fragments in the introns of functional genes could impact the activity of functional genes, and a large number of TE fragments in introns could lead to the indirect inactivation of functional genes. The present study will be beneficial in understanding the role and impact of TEs in genomic evolution of their hosts.

Selection of Reference Genes for Normalization of qRT‒PCR Analysis in the Soybean Aphid Aphis glycines Matsumura (Hemiptera: Aphididae)

The soybean aphid Aphis glycines Matsumura is a predominant insect pest in Asia and North America and causes great losses to soybean. The release of genome data for A. glycines will facilitate gene function research in the future. However, suitable reference genes for A. glycines under various experimental conditions are scarce. To search for appropriate reference genes for A. glycines, nine candidate reference genes, including Act, α-Tub, β-Tub, RPS12, RPS18, RPL5, RPL27, EF1α, and Fer, were tested under six experimental conditions to evaluate their suitability for use in the normalization of qRT‒PCR data. Results showed that EF1α and RPS12 were optimal for the developmental stages of A. glycines, RPS18 and RPS12 were appropriate for wing dimorphism, β-Tub and RPS18 were suitable for different tissues and RPL5, and α-Tub could be used for normalization at different temperatures. β-Tub and EF1α could be proposed as reference genes for insecticide treatment, and RPL5 and RPS12 were found to be the most stable reference genes in different photoperiods. The results provide appropriate reference genes for analyzing gene expression in A. glycines and contribute to future research on the molecular physiology and biochemistry of A. glycines.

Octopamine receptor genes are involved in the starvation response of Rhopalosiphum padi (Hemiptera: Aphididae)

Insect octopamine (OA) receptors are G-protein coupled receptors (GPCRs) that play essential roles in physiological and behavioural processes. However, there is little information about the function of OA receptors in the aphids' response to stress. From the genome sequence of Rhopalosiphum padi genome sequence, a cosmopolitan cereal pest, we identified six OA receptor genes RpOAMB, RpOctR, RpOctβ1R, RpOctβ2R, RpOctβ3R, RpOctR-like with two, one, one, four, four, seven exons, respectively. All the OA receptors contain seven transmembrane domains, which were the signature of GPCRs. Our results showed that (1) the contents of OA increased significantly after food starvation, (2) the transcription levels of RpOAMB, RpOctR, RpOctβ2R and RpOctβ3R increased after starvation and were restored after re-feeding, and (3) the expression levels of these four genes decreased significantly 48 h post-injection of dsRNA that targeted the respective genes. Knockdown of RpOctR, RpOctβ2R or RpOctβ3R genes significantly increased aphid mortality under 24 h starvation conditions. Mortality of R. padi injected with dsRpOctR or dsRpOctβ2R was significantly higher than control under 48 h starvation treatments. This is the first report on the role of OA receptors in the starvation response of aphids. The current study provides knowledge for a better understanding the physiological roles of insect OA receptors.

Influence on Soybean Aphid by the Tripartite Interaction between Soybean, a Rhizobium Bacterium, and an Arbuscular Mycorrhizal Fungus

The inoculation of arbuscular mycorrhizal (AM) fungi and rhizobia in legumes has been proven to increase plant growth and yield. To date, studies of the effects of these interactions on phytophagous insects have shown them to be context-dependent depending on the inoculant strain, the plant, and the insect species. Here, we document how a symbiosis involving an AM fungus, Rhizophagus irregularis; a rhizobium, Bradyrhizobium japonicum; and soybean, Glycine max, influences the soybean aphid, Aphis glycines. Soybean co-inoculated with the AM fungus–rhizobium pair increased the plant’s biomass, nodulation, mycorrhizal colonization, nitrogen, and carbon concentrations, but decreased phosphorus concentration. Similar effects were observed with rhizobium alone, with the exception that root biomass was unaffected. With AM fungus alone, we only observed an increase in mycorrhizal colonization and phosphorus concentration. The aphids experienced an increased reproductive rate with the double inoculation, followed by rhizobium alone, whereas no effect was observed with the AM fungus. The size of individual aphids was not affected. Furthermore, we found positive correlation between nitrogen concentration and aphid population density. Our results confirm that co-inoculation of two symbionts can enhance both plant and phytophagous insect performance beyond what either symbiont can contribute alone.

Chromosome‐level genome assembly for the horned‐gall aphid provides insights into interactions between gall‐making insect and its host plant

The aphid Schlechtendalia chinensis is an economically important insect that can induce horned galls, which are valuable for the medicinal and chemical industries. Up to now, more than twenty aphid genomes have been reported. Most of the sequenced genomes are derived from free‐living aphids. Here, we generated a high‐quality genome assembly from a galling aphid. The final genome assembly is 271.52 Mb, representing one of the smallest sequenced genomes of aphids. The genome assembly is based on contig and scaffold N50 values of the genome sequence are 3.77 Mb and 20.41 Mb, respectively. Nine‐seven percent of the assembled sequences was anchored onto 13 chromosomes. Based on BUSCO analysis, the assembly involved 96.9% of conserved arthropod and 98.5% of the conserved Hemiptera single‐copy orthologous genes. A total of 14,089 protein‐coding genes were predicted. Phylogenetic analysis revealed that S. chinensis diverged from the common ancestor of Eriosoma lanigerum approximately 57 million years ago (MYA). In addition, 35 genes encoding salivary gland proteins showed differentially when S. chinensis forms a gall, suggesting they have potential roles in gall formation and plant defense suppression. Taken together, this high‐quality S. chinensis genome assembly and annotation provide a solid genetic foundation for future research to reveal the mechanism of gall formation and to explore the interaction between aphids and their host plants.

Aphid hologenomics: current status and future challenges

Aphids are important model organisms in ecological, developmental, and evolutionary studies of, for example, symbiosis, insect-plant interactions, pest management, and developmental polyphenism. Here, we review the recent progress made in the genomics of aphids and their symbionts: hologenomics. The reference genome of Acyrthosiphon pisum has been greatly improved, and chromosome-level assembly is now available. The genomes of over 20 aphid species have been sequenced, and comparative genomic analyses have revealed pervasive gene duplication and dynamic chromosomal rearrangements. Over 120 symbiont genomes (both obligate and facultative) have been sequenced, and modern deep-sequencing technologies have identified novel symbionts. The advances in hologenomics have helped to elucidate the dynamic evolution of facultative and co-obligate symbionts with the ancient obligate symbiont Buchnera.

Genome sequence of the English grain aphid, Sitobion avenae and its endosymbiont Buchnera aphidicola

The English grain aphid, Sitobion avenae, is a major agricultural pest of wheat, barley and oats, and one of the principal vectors of barley yellow dwarf virus leading to significant reductions in grain yield, annually. Emerging resistance to and increasing regulation of insecticides has resulted in limited options for their control. Using PacBio HiFi data, we have produced a high-quality draft assembly of the S. avenae genome; generating a primary assembly with a total assembly size of 475.7 Mb, and an alternate assembly with a total assembly size of 430.8 Mb. Our primary assembly was highly contiguous with only 326 contigs and a contig N50 of 15.95 Mb. Assembly completeness was estimated at 97.7% using BUSCO analysis and 31,007 and 29,037 protein-coding genes were predicted from the primary and alternate assemblies, respectively. This assembly, which is to our knowledge the first for an insecticide resistant clonal lineage of English grain aphid, will provide novel insight into the molecular and mechanistic determinants of resistance and will facilitate future research into mechanisms of viral transmission and aphid behavior.

The phylogenetic and evolutionary analyses of detoxification gene families in Aphidinae species

Detoxification enzymes play significant roles in the interactions between insects and host plants, wherein detoxification-related genes make great contributions. As herbivorous pests, aphids reproduce rapidly due to parthenogenesis. They are good biological materials for studying the mechanisms that allow insect adaptation to host plants. Insect detoxification gene families are associated with insect adaptation to host plants. The Aphidinae is the largest subfamily in the Aphididae with at least 2483 species in 256 genera in 2 tribes: the Macrosiphini (with 3/4 of the species) and the Aphidini. Most aphid pests on crops and ornamental plants are Aphidinae. Members of the Aphidinae occur in nearly every region of the world. The body shape and colour vary significantly. To research the role that detoxification gene families played in the process of aphid adaptation to host evolution, we analyzed the phylogeny and evolution of these detoxification gene families in Aphidinae. In general, the P450/GST/CCE gene families contract, whereas the ABC/UGT families are conserved in Aphidinae species compared to these families in other herbivorous insects. Genus-specific expansions of P450 CYP4, and GST Delta have occurred in the genus Acyrthosiphon. In addition, the evolutionary rates of five detoxification gene families in the evolution process of Aphidinae are different. The comparison of five detoxification gene families among nine Aphidinae species and the estimated relative evolutionary rates provided herein support an understanding of the interaction between and the co-evolution of Aphidinae and plants.

A scaffold-level genome assembly of a minute pirate bug, Orius laevigatus (Hemiptera: Anthocoridae), and a comparative analysis of insecticide resistance-related gene families with hemipteran crop pests

BACKGROUND

Orius laevigatus, a minute pirate bug, is a highly effective beneficial predator of crop pests including aphids, spider mites and thrips in integrated pest management (IPM) programmes. No genomic information is currently available for O. laevigatus, as is the case for the majority of beneficial predators which feed on crop pests. In contrast, genomic information for crop pests is far more readily available. The lack of publicly available genomes for beneficial predators to date has limited our ability to perform comparative analyses of genes encoding potential insecticide resistance mechanisms between crop pests and their predators. These mechanisms include several gene/protein families including cytochrome P450s (P450s), ATP binding cassette transporters (ABCs), glutathione S-transferases (GSTs), UDP-glucosyltransferases (UGTs) and carboxyl/cholinesterases (CCEs).

METHODS AND FINDINGS

In this study, a high-quality scaffold level de novo genome assembly for O. laevigatus has been generated using a hybrid approach with PacBio long-read and Illumina short-read data. The final assembly achieved a scaffold N50 of 125,649 bp and a total genome size of 150.98 Mb. The genome assembly achieved a level of completeness of 93.6% using a set of 1658 core insect genes present as full-length genes. Genome annotation identified 15,102 protein-coding genes - 87% of which were assigned a putative function. Comparative analyses revealed gene expansions of sigma class GSTs and CYP3 P450s. Conversely the UGT gene family showed limited expansion. Differences were seen in the distributions of resistance-associated gene families at the subfamily level between O. laevigatus and some of its targeted crop pests. A target site mutation in ryanodine receptors (I4790M, PxRyR) which has strong links to diamide resistance in crop pests and had previously only been identified in lepidopteran species was found to also be present in hemipteran species, including O. laevigatus.

CONCLUSION AND SIGNIFICANCE

This assembly is the first published genome for the Anthocoridae family and will serve as a useful resource for further research into target-site selectivity issues and potential resistance mechanisms in beneficial predators. Furthermore, the expansion of gene families often linked to insecticide resistance may be an indicator of the capacity of this predator to detoxify selective insecticides. These findings could be exploited by targeted pesticide screens and functional studies to increase effectiveness of IPM strategies, which aim to increase crop yields by sustainably, environmentally-friendly and effectively control pests without impacting beneficial predator populations.

Chromosome-level genome assemblies of two cotton-melon aphid Aphis gossypii biotypes unveil mechanisms of host adaption

The cotton-melon aphid Aphis gossypii is a sap-sucking insect that is considered a serious global pest. The species is distributed over a large geographical range and uses a wide variety of hosts, with some populations being specialized to attack different plant species. Here, we provide de novo chromosome-level genome assemblies of a cotton specialist population (Hap1) and a cucurbit specialist population (Hap3). We achieved this by using a combination of third-generation sequencing platforms, namely Illumina and Hi-C sequencing technologies. We were able to anchor a total of 334.89 Mb (scaffold N50 of 89.13 Mb) and 359.95 Mb (scaffold N50 of 68.88 Mb) to four chromosomes for Hap1 and Hap3, respectively. Moreover, our results showed that the X-chromosome of Hap3 (113.01 Mb) was significantly longer than that of Hap1 (100.26 Mb), with a high level of sequence conservation between the aphid species. We also report variation in the number of protein-coding genes and repeat sequences between Hap1 and Hap3. In particular, olfactory and gustatory receptor genes underwent a high level of gene duplication and expansion events in A. gossypii, including between Hap1 and Hap3. Moreover, we identified two glutathione S-transferase genes which underwent single gene duplications in Hap3, and tandem duplication and inversion events affecting the cytochrome P450 monooxygenase between Hap1 and Hap3, all of which include the CYP3 family. Our results illustrate the variance in the genomic composition of two specialized A. gossypii populations and provide a helpful resource for the study of aphid population evolution, host adaption and insecticide resistance.

Involvement of chemosensory proteins in host plant searching in the bird cherry-oat aphid

Chemosensory systems are considered to play an important role in host plant selection in herbivorous insects. However, few studies have focused on chemosensory proteins (CSPs) for aphid host-location mechanisms. The roles of CSPs in searching for different Poaceae species (wheat, barley, triticale, maize and sorghum) were tested in Rhopalosiphum padi, an important cereal pest. The olfactometer assays showed that R. padi responds to plant odors. Seven R. padi CSP genes were identified. Influence of aphid morph, tissue and starvation state on expression patterns of CSPs was evaluated. Expression levels of CSP1, CSP4, CSP5 and CSP6 in winged aphids were significantly higher than those in wingless ones. Transcription levels of four genes (CSP1, CSP4, CSP5 and CSP6) were relatively higher in the head with antennae, and the four genes tended to be upregulated following starvation. Silencing of three CSPs (CSP4, CSP5 and CSP6) altered aphid host-location behavior in response to the five different host plants tested. Three volatile compounds of host plants (octanal, [E]-2-hexenol and linalool) have significant attraction to winged R. padi according to the four-arm olfactometer tests. Molecular docking predicted hydrogen bonding sites which played key roles in the binding of CSP4, CSP5 and CSP6 with volatile compounds. Knockdown of CSP4 or CSP5 significantly decreased the staying time of R. padi in the arms with octanal. However, knockdown of CSP6 could not affect the response of R. padi to octanal. These results bring evidence for the involvement of three CSPs in R. padi host-location behavior.

Next-Generation Sequencing and Its Impacts on Entomological Research in Ecology and Evolution

The advent of NGS-based methods has been profoundly transforming entomological research. Through continual development and improvement of different methods and sequencing platforms, NGS has promoted mass elucidation of partial or whole genetic materials associated with beneficial insects, pests (of agriculture, forestry and animal, and human health), and species of conservation concern, helping to unravel ecological and evolutionary mechanisms and characterizing survival, trophic interactions, and dispersal. It is shifting the scale of biodiversity and environmental analyses from individuals and biodiversity indicator species to the large-scale study of communities and ecosystems using bulk samples of species or a mixed "soup" of environmental DNA. As the NGS-based methods have become more affordable, complexity demystified, and specificity and sensitivity proven, their use in entomological research has spread widely. This article presents several examples on how NGS-based methods have been used in entomology to provide incentives to apply them when appropriate and to open our minds to the expected advances in entomology that are yet to come.

A Detailed Spatial Expression Analysis of Wing Phenotypes Reveals Novel Patterns of Odorant Binding Proteins in the Soybean Aphid, Aphis glycines

The wide range of insect niches has led to a rapid expansion of chemosensory gene families as well as their relatively independent evolution and a high variation. Previous studies have revealed some functions for odorant-binding proteins (OBPs) in processes beyond olfaction, such as gustation and reproduction. In this study, a comparative transcriptomic analysis strategy was applied for the soybean aphid, Aphis glycines, focusing on various functional tissues and organs of winged aphids, including the antenna, head, leg, wing, thorax, cauda, and cornicle. Detailed spatial OBP expression patterns in winged and wingless parthenogenetic aphids were detected by RT-qPCR. Twelve OBPs were identified, and three new OBPs in A. glycines are first reported. All OBPs showed comparatively higher expression in sensory organs and tissues, such as the antenna, head, or leg. Additionally, we found some novel expression patterns for aphid OBPs (Beckendorf et al., 2008). Five OBPs exhibited high-expression levels in the cauda and four in the cornicle (Biasio et al., 2015). Three genes (OBP2/3/15) were highly expressed in the wing (Calvello et al., 2003). Two (OBP3/15) were significantly more highly expressed in the wingless thorax than in the winged thorax with the wings removed, and these transcripts were significantly enriched in the removed wings. More details regarding OBP spatial expression were revealed under our strategy. These findings supported the existence of carrier transport functions other than for foreign chemicals and therefore broader ligand ranges of aphid OBPs. It is important for understanding how insect OBPs function in chemical perception as well as their other potential physiological functions.

Triple RNA-Seq characterizes aphid gene expression in response to infection with unequally virulent strains of the endosymbiont Hamiltonella defensa

Background

Secondary endosymbionts of aphids provide benefits to their hosts, but also impose costs such as reduced lifespan and reproductive output. The aphid Aphis fabae is host to different strains of the secondary endosymbiont Hamiltonella defensa, which encode different putative toxins. These strains have very different phenotypes: They reach different densities in the host, and the costs and benefits (protection against parasitoid wasps) they confer to the host vary strongly.

Results

We used RNA-Seq to generate hypotheses on why four of these strains inflict such different costs to A. fabae. We found different H. defensa strains to cause strain-specific changes in aphid gene expression, but little effect of H. defensa on gene expression of the primary endosymbiont, Buchnera aphidicola. The highly costly and over-replicating H. defensa strain H85 was associated with strongly reduced aphid expression of hemocytin, a marker of hemocytes in Drosophila. The closely related strain H15 was associated with downregulation of ubiquitin-related modifier 1, which is related to nutrient-sensing and oxidative stress in other organisms. Strain H402 was associated with strong differential regulation of a set of hypothetical proteins, the majority of which were only differentially regulated in presence of H402.

Conclusions

Overall, our results suggest that costs of different strains of H. defensa are likely caused by different mechanisms, and that these costs are imposed by interacting with the host rather than the host’s obligatory endosymbiont B. aphidicola.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07742-8.

Aphid resistance is the future for soybean production, and has been since 2004: efforts towards a wider use of host plant resistance in soybean

The soybean aphid (Aphis glycines) is an important pest of soybeans in the Midwestern US. The first aphid resistance genes were identified in the early 21st century and resistant varieties have been commercially available for 10 years, but have been very underutilized. Major seed companies have avoided commercializing aphid resistant soybean varieties for conventional farmers (i.e., not organic), in part because of the discovery of virulent biotypes in North America. The emergence of soybean aphid populations resistant to insecticides creates a greater incentive for the use of host plant resistance. New research on aphid genetics and markers, plant gene expression and in-plant refuges, suggest important avenues for insect resistance management (IRM) which may encourage more widescale commercialization of this valuable pest management tool.

Chemosensory proteins participate in insecticide susceptibility in Rhopalosiphum padi, a serious pest on wheat crops

Rhopalosiphum padi is a worldwide agricultural pest. Chemosensory proteins (CSPs) are considered to be a type of transporters which can bind chemicals from external environments. Previous research showed that the expression of some insect CSPs were significantly increased after exposure to insecticides, and CSPs were involved in insecticide resistance or susceptibility. However, the role of CSPs in the susceptibility and response of R. padi to insecticides is still unknown. In this study, we identified eight CSP (RpCSP) from R. padi by genome-wide investigation. Seven RpCSP genes had two exons, while RpCSP7 had three exons. qPCR analyses showed that the mRNA levels of the eight RpCSP genes were significantly affected by imidacloprid and beta-cypermethrin in different post-treatment periods. Molecular docking predicted that there were hydrogen bonding sites which played key roles in binding of RpCSP4, RpCSP5, RpCSP6, RpCSP7 and RpCSP10 with imidacloprid and beta-cypermethrin. Knockdown of RpCSP4, RpCSP5, RpCSP6 and RpCSP10 by RNA interference significantly increased the aphid mortality under two sublethal concentrations of imidacloprid. Mortalities under two sublethal concentrations of beta-cypermethrin conditions were significantly higher after injection of R. padi with dsCSP4 and dsCSP6. The results indicate that some RpCSP genes are involved in the insecticide susceptibility of R. padi.

Chromosome-Scale Genome Assemblies of Aphids Reveal Extensively Rearranged Autosomes and Long-Term Conservation of the X Chromosome

Chromosome rearrangements are arguably the most dramatic type of mutations, often leading to rapid evolution and speciation. However, chromosome dynamics have only been studied at the sequence level in a small number of model systems. In insects, Diptera and Lepidoptera have conserved genome structure at the scale of whole chromosomes or chromosome arms. Whether this reflects the diversity of insect genome evolution is questionable given that many species exhibit rapid karyotype evolution. Here, we investigate chromosome evolution in aphids-an important group of hemipteran plant pests-using newly generated chromosome-scale genome assemblies of the green peach aphid (Myzus persicae) and the pea aphid (Acyrthosiphon pisum), and a previously published assembly of the corn-leaf aphid (Rhopalosiphum maidis). We find that aphid autosomes have undergone dramatic reorganization over the last 30 My, to the extent that chromosome homology cannot be determined between aphids from the tribes Macrosiphini (Myzus persicae and Acyrthosiphon pisum) and Aphidini (Rhopalosiphum maidis). In contrast, gene content of the aphid sex (X) chromosome remained unchanged despite rapid sequence evolution, low gene expression, and high transposable element load. To test whether rapid evolution of genome structure is a hallmark of Hemiptera, we compared our aphid assemblies with chromosome-scale assemblies of two blood-feeding Hemiptera (Rhodnius prolixus and Triatoma rubrofasciata). Despite being more diverged, the blood-feeding hemipterans have conserved synteny. The exceptional rate of structural evolution of aphid autosomes renders them an important emerging model system for studying the role of large-scale genome rearrangements in evolution.

Genome Sequence of the Banana Aphid, Pentalonia nigronervosa Coquerel (Hemiptera: Aphididae) and Its Symbionts

The banana aphid, Pentalonia nigronervosa Coquerel (Hemiptera: Aphididae), is a major pest of cultivated bananas (Musa spp., order Zingiberales), primarily due to its role as a vector of Banana bunchy top virus (BBTV), the most severe viral disease of banana worldwide. Here, we generated a highly complete genome assembly of P. nigronervosa using a single PCR-free Illumina sequencing library. Using the same sequence data, we also generated complete genome assemblies of the P. nigronervosa symbiotic bacteria Buchnera aphidicola and Wolbachia. To improve our initial assembly of P. nigronervosa we developed a k-mer based deduplication pipeline to remove genomic scaffolds derived from the assembly of haplotigs (allelic variants assembled as separate scaffolds). To demonstrate the usefulness of this pipeline, we applied it to the recently generated assembly of the aphid Myzus cerasi, reducing the duplication of conserved BUSCO genes by 25%. Phylogenomic analysis of P. nigronervosa, our improved M. cerasi assembly, and seven previously published aphid genomes, spanning three aphid tribes and two subfamilies, reveals that P. nigronervosa falls within the tribe Macrosiphini, but is an outgroup to other Macrosiphini sequenced so far. As such, the genomic resources reported here will be useful for understanding both the evolution of Macrosphini and for the study of P. nigronervosa. Furthermore, our approach using low cost, high-quality, Illumina short-reads to generate complete genome assemblies of understudied aphid species will help to fill in genomic black spots in the diverse aphid tree of life.

A chromosome-level genome assembly of the woolly apple aphid, Eriosoma lanigerum Hausmann (Hemiptera: Aphididae)

Woolly apple aphid (WAA, Eriosoma lanigerum Hausmann) (Hemiptera: Aphididae) is a major pest of apple trees (Malus domestica, order Rosales) and is critical to the economics of the apple industry in most parts of the world. Here, we generated a chromosome-level genome assembly of WAA-representing the first genome sequence from the aphid subfamily Eriosomatinae-using a combination of 10X Genomics linked-reads and in vivo Hi-C data. The final genome assembly is 327 Mb, with 91% of the assembled sequences anchored into six chromosomes. The contig and scaffold N50 values are 158 kb and 71 Mb, respectively, and we predicted a total of 28,186 protein-coding genes. The assembly is highly complete, including 97% of conserved arthropod single-copy orthologues based on Benchmarking Universal Single-Copy Orthologs (busco) analysis. Phylogenomic analysis of WAA and nine previously published aphid genomes, spanning four aphid tribes and three subfamilies, reveals that the tribe Eriosomatini (represented by WAA) is recovered as a sister group to Aphidini + Macrosiphini (subfamily Aphidinae). We identified syntenic blocks of genes between our WAA assembly and the genomes of other aphid species and find that two WAA chromosomes (El5 and El6) map to the conserved Macrosiphini and Aphidini X chromosome. Our high-quality WAA genome assembly and annotation provides a valuable resource for research in a broad range of areas such as comparative and population genomics, insect-plant interactions and pest resistance management.

Apolygus lucorum genome provides insights into omnivorousness and mesophyll feeding

Apolygus lucorum (Miridae) is an omnivorous pest that occurs worldwide and is notorious for the serious damage it causes to various crops and substantial economic losses. Although some studies have examined the biological characteristics of the mirid bug, no reference genome is available in Miridae, limiting in-depth studies of this pest. Here, we present a chromosome-scale reference genome of A. lucorum, the first sequenced Miridae species. The assembled genome size was 1.02 Gb with a contig N50 of 785 kb. With Hi-C scaffolding, 1,016 Mb contig sequences were clustered, ordered and assembled into 17 large scaffolds with scaffold N50 length 68 Mb, each corresponding to a natural chromosome. Numerous transposable elements occur in this genome and contribute to the large genome size. Expansions of genes associated with omnivorousness and mesophyll feeding such as those related to digestion, chemosensory perception, and detoxification were observed in A. lucorum, suggesting that gene expansion contributed to its strong environmental adaptability and severe harm to crops. We clarified that a salivary enzyme polygalacturonase is unique in mirid bugs and has significantly expanded in A. lucorum, which may contribute to leaf damage from this pest. The reference genome of A. lucorum not only facilitates biological studies of Hemiptera as well as an understanding of the damage mechanism of mesophyll feeding, but also provides a basis on which to develop efficient control technologies for mirid bugs.

Big Genes, Small Effectors: Pea Aphid Cassette Effector Families Composed From Miniature Exons

Aphids secrete proteins from their stylets that evidence indicates function similar to pathogen effectors for virulence. Here, we describe two small candidate effector gene families of the pea aphid, Acyrthosiphon pisum, that share highly conserved secretory signal peptide coding regions and divergent non-secretory coding sequences derived from miniature exons. The KQY candidate effector family contains eleven members with additional isoforms, generated by alternative splicing. Pairwise comparisons indicate possible four unique KQY families based on coding regions without the secretory signal region. KQY1a, a representative of the family, is encoded by a 968 bp mRNA and a gene that spans 45.7 kbp of the genome. The locus consists of 37 exons, 33 of which are 15 bp or smaller. Additional KQY members, as well as members of the KHI family, share similar features. Differential expression analyses indicate that the genes are expressed preferentially in salivary glands. Proteomic analysis on salivary glands and saliva revealed 11 KQY members in salivary proteins, and KQY1a was detected in an artificial diet solution after aphid feeding. A single KQY locus and two KHI loci were identified in Myzus persicae, the peach aphid. Of the genes that can be anchored to chromosomes, loci are mostly scattered throughout the genome, except a two-gene region (KQY4/KQY6). We propose that the KQY family expanded in A. pisum through combinatorial assemblies of a common secretory signal cassette and novel coding regions, followed by classical gene duplication and divergence.

The genome sequence of the grape phylloxera provides insights into the evolution, adaptation, and invasion routes of an iconic pest

BACKGROUND

Although native to North America, the invasion of the aphid-like grape phylloxera Daktulosphaira vitifoliae across the globe altered the course of grape cultivation. For the past 150 years, viticulture relied on grafting-resistant North American Vitis species as rootstocks, thereby limiting genetic stocks tolerant to other stressors such as pathogens and climate change. Limited understanding of the insect genetics resulted in successive outbreaks across the globe when rootstocks failed. Here we report the 294-Mb genome of D. vitifoliae as a basic tool to understand host plant manipulation, nutritional endosymbiosis, and enhance global viticulture.

RESULTS

Using a combination of genome, RNA, and population resequencing, we found grape phylloxera showed high duplication rates since its common ancestor with aphids, but similarity in most metabolic genes, despite lacking obligate nutritional symbioses and feeding from parenchyma. Similarly, no enrichment occurred in development genes in relation to viviparity. However, phylloxera evolved > 2700 unique genes that resemble putative effectors and are active during feeding. Population sequencing revealed the global invasion began from the upper Mississippi River in North America, spread to Europe and from there to the rest of the world.

CONCLUSIONS

The grape phylloxera genome reveals genetic architecture relative to the evolution of nutritional endosymbiosis, viviparity, and herbivory. The extraordinary expansion in effector genes also suggests novel adaptations to plant feeding and how insects induce complex plant phenotypes, for instance galls. Finally, our understanding of the origin of this invasive species and its genome provide genetics resources to alleviate rootstock bottlenecks restricting the advancement of viticulture.

Influence of cell wall polymers and their modifying enzymes during plant-aphid interactions

Aphids are a major issue for commercial crops. These pests drain phloem nutrients and transmit ~50% of the known insect-borne viral diseases. During aphid feeding, trophic structures called stylets advance toward the phloem intercellularly, disrupting cell wall polymers. It is thought that cell wall-modifying enzymes (CWMEs) present in aphid saliva facilitate stylet penetration through this intercellular polymer network. Additionally, different studies have demonstrated that host settling preference, feeding behavior, and colony performance of aphids are influenced by modulating the CWME expression levels in host plants. CWMEs have been described as critical defensive elements for plants, but also as a key virulence factor for plant pathogens. However, whether CWMEs are elements of the plant defense mechanisms or the aphid infestation process remains unclear. Therefore, in order to better consider the function of CWMEs and cell wall-derived damage-associated molecular patterns (DAMPs) during plant-aphid interactions, the present review integrates different hypotheses, perspectives, and experimental evidence in the field of plant-aphid interactions and discusses similarities to other well-characterized models such as the fungi-plant pathosystems from the host and the attacker perspectives.

AcDCXR Is a Cowpea Aphid Effector With Putative Roles in Altering Host Immunity and Physiology

Cowpea, Vigna unguiculata, is a crop that is essential to semiarid areas of the world like Sub-Sahara Africa. Cowpea is highly susceptible to cowpea aphid, Aphis craccivora, infestation that can lead to major yield losses. Aphids feed on their host plant by inserting their hypodermal needlelike flexible stylets into the plant to reach the phloem sap. During feeding, aphids secrete saliva, containing effector proteins, into the plant to disrupt plant immune responses and alter the physiology of the plant to their own advantage. Liquid chromatography tandem mass spectrometry (LC-MS/MS) was used to identify the salivary proteome of the cowpea aphid. About 150 candidate proteins were identified including diacetyl/L-xylulose reductase (DCXR), a novel enzyme previously unidentified in aphid saliva. DCXR is a member of short-chain dehydrogenases/reductases with dual enzymatic functions in carbohydrate and dicarbonyl metabolism. To assess whether cowpea aphid DCXR (AcDCXR) has similar functions, recombinant AcDCXR was purified and assayed enzymatically. For carbohydrate metabolism, the oxidation of xylitol to xylulose was tested. The dicarbonyl reaction involved the reduction of methylglyoxal, an α-β-dicarbonyl ketoaldehyde, known as an abiotic and biotic stress response molecule causing cytotoxicity at high concentrations. To assess whether cowpea aphids induce methylglyoxal in plants, we measured methylglyoxal levels in both cowpea and pea (Pisum sativum) plants and found them elevated transiently after aphid infestation. Agrobacterium-mediated transient overexpression of AcDCXR in pea resulted in an increase of cowpea aphid fecundity. Taken together, our results indicate that AcDCXR is an effector with a putative ability to generate additional sources of energy to the aphid and to alter plant defense responses. In addition, this work identified methylglyoxal as a potential novel aphid defense metabolite adding to the known repertoire of plant defenses against aphid pests.

Brown marmorated stink bug, Halyomorpha halys (Stål), genome: putative underpinnings of polyphagy, insecticide resistance potential and biology of a top worldwide pest

BACKGROUND

Halyomorpha halys (Stål), the brown marmorated stink bug, is a highly invasive insect species due in part to its exceptionally high levels of polyphagy. This species is also a nuisance due to overwintering in human-made structures. It has caused significant agricultural losses in recent years along the Atlantic seaboard of North America and in continental Europe. Genomic resources will assist with determining the molecular basis for this species' feeding and habitat traits, defining potential targets for pest management strategies.

RESULTS

Analysis of the 1.15-Gb draft genome assembly has identified a wide variety of genetic elements underpinning the biological characteristics of this formidable pest species, encompassing the roles of sensory functions, digestion, immunity, detoxification and development, all of which likely support H. halys' capacity for invasiveness. Many of the genes identified herein have potential for biomolecular pesticide applications.

CONCLUSIONS

Availability of the H. halys genome sequence will be useful for the development of environmentally friendly biomolecular pesticides to be applied in concert with more traditional, synthetic chemical-based controls.

Cuticular Structure Proteomics in the Pea Aphid Acyrthosiphon pisum Reveals New Plant Virus Receptor Candidates at the Tip of Maxillary Stylets

Aphids are phloem-feeding insects known as major pests in agriculture that are able to transmit hundreds of plant viruses. The majority of these viruses, classified as noncirculative, are retained and transported on the inner surface of the cuticle of the needle-like mouthparts while the aphids move from plant to plant. Identification of receptors of viruses within insect vectors is a key challenge because they are promising targets for alternative control strategies. The acrostyle, an organ discovered earlier within the common food/salivary canal at the tip of aphid maxillary stylets, displays proteins at the cuticle-fluid interface, some of which are receptors of noncirculative viruses. To assess the presence of stylet- and acrostyle-specific proteins and identify putative receptors, we have developed a comprehensive comparative analysis of the proteomes of four cuticular anatomical structures of the pea aphid, stylets, antennae, legs, and wings. In addition, we performed systematic immunolabeling detection of the cuticular proteins identified by mass spectrometry in dissected stylets. We thereby establish the first proteome of stylets of an insect and determine the minimal repertoire of the cuticular proteins composing the acrostyle. Most importantly, we propose a short list of plant virus receptor candidates, among which RR-1 proteins are remarkably predominant. The data are available via ProteomeXchange (PXD016517).

Improved Genome Assembly and Annotation of the Soybean Aphid (Aphis glycines Matsumura)

Aphids are an economically important insect group due to their role as plant disease vectors. Despite this economic impact, genomic resources have only been generated for a small number of aphid species. The soybean aphid (Aphis glycines Matsumura) was the third aphid species to have its genome sequenced and the first to use long-read sequence data. However, version 1 of the soybean aphid genome assembly has low contiguity (contig N50 = 57 Kb, scaffold N50 = 174 Kb), poor representation of conserved genes and the presence of genomic scaffolds likely derived from parasitoid wasp contamination. Here, I use recently developed methods to reassemble the soybean aphid genome. The version 2 genome assembly is highly contiguous, containing half of the genome in only 40 scaffolds (contig N50 = 2.00 Mb, scaffold N50 = 2.51 Mb) and contains 11% more conserved single-copy arthropod genes than version 1. To demonstrate the utility of this improved assembly, I identify a region of conserved synteny between aphids and Drosophila containing members of the Osiris gene family that was split over multiple scaffolds in the original assembly. The improved genome assembly and annotation of A. glycines demonstrates the benefit of applying new methods to old data sets and will provide a useful resource for future comparative genome analysis of aphids.

Phylogenomics Identifies an Ancestral Burst of Gene Duplications Predating the Diversification of Aphidomorpha

Aphids (Aphidoidea) are a diverse group of hemipteran insects that feed on plant phloem sap. A common finding in studies of aphid genomes is the presence of a large number of duplicated genes. However, when these duplications occurred remains unclear, partly due to the high relatedness of sequenced species. To better understand the origin of aphid duplications we sequenced and assembled the genome of Cinara cedri, an early branching lineage (Lachninae) of the Aphididae family. We performed a phylogenomic comparison of this genome with 20 other sequenced genomes, including the available genomes of five other aphids, along with the transcriptomes of two species belonging to Adelgidae (a closely related clade to the aphids) and Coccoidea. We found that gene duplication has been pervasive throughout the evolution of aphids, including many parallel waves of recent, species-specific duplications. Most notably, we identified a consistent set of very ancestral duplications, originating from a large-scale gene duplication predating the diversification of Aphidomorpha (comprising aphids, phylloxerids, and adelgids). Genes duplicated in this ancestral wave are enriched in functions related to traits shared by Aphidomorpha, such as association with endosymbionts, and adaptation to plant defenses and phloem-sap-based diet. The ancestral nature of this duplication wave (106-227 Ma) and the lack of sufficiently conserved synteny make it difficult to conclude whether it originated from a whole-genome duplication event or, alternatively, from a burst of large-scale segmental duplications. Genome sequencing of other aphid species belonging to different Aphidomorpha and related lineages may clarify these findings.

Genome sequence of the corn leaf aphid (Rhopalosiphum maidis Fitch)

Background

The corn leaf aphid (Rhopalosiphum maidis Fitch) is the most economically damaging aphid pest on maize (Zea mays), one of the world's most important grain crops. In addition to causing direct damage by removing photoassimilates, R. maidis transmits several destructive maize viruses, including maize yellow dwarf virus, barley yellow dwarf virus, sugarcane mosaic virus, and cucumber mosaic virus.

Findings

The genome of a parthenogenetically reproducing R. maidis clone was assembled with a combination of Pacific Biosciences (207-fold coverage) and Illumina (83-fold coverage) sequencing. The 689 assembled contigs, which have an N50 size of 9.0 megabases (Mb) and a low level of heterozygosity, were clustered using Phase Genomics Hi-C interaction maps. Consistent with the commonly observed 2n = 8 karyotype of R. maidis, most of the contigs (473 spanning 321 Mb) were successfully oriented into 4 scaffolds. The genome assembly captured the full length of 95.8% of the core eukaryotic genes, indicating that it is highly complete. Repetitive sequences accounted for 21.2% of the assembly, and a total of 17,629 protein-coding genes were predicted with integrated evidence from ab initio and homology-based gene predictions and transcriptome sequences generated with both Pacific Biosciences and Illumina. An analysis of likely horizontally transferred genes identified 2 from bacteria, 7 from fungi, 2 from protozoa, and 9 from algae. Repeat elements, transposons, and genes encoding likely detoxification enzymes (cytochrome P450s, glutathione S-transferases, carboxylesterases, uridine diphosphate–glucosyltransferases, and ABC transporters) were identified in the genome sequence. Other than Buchnera aphidicola (642,929 base pairs, 602 genes), no endosymbiont bacteria were found in R. maidis.

Conclusions

A high-quality R. maidis genome was assembled at the chromosome level. This genome sequence will enable further research related to ecological interactions, virus transmission, pesticide resistance, and other aspects of R. maidis biology. It also serves as a valuable resource for comparative investigation of other aphid species.

Topical dsRNA delivery induces gene silencing and mortality in the pea aphid

BACKGROUND

With the growing number of available aphid genomes and transcriptomes, an efficient and easy-to-adapt tool for gene function study is urgently required. RNA interference (RNAi), as a post-transcriptional gene silencing mechanism, is important as a research tool for determining gene functions and has potential as a novel insect control strategy. However, these applications have been hampered by the lack of effective dsRNA delivery approaches in aphids.

RESULTS

Here, we developed a convenient and efficient dsRNA delivery method, topical RNAi, in aphids. An investigation of its dose and time-dependent RNAi efficiencies revealed that with as little as 60 ng dsRNA per adult pea aphid (Acyrthosiphon pisum), the indicator gene, Aphunchback, could be significantly silenced within 2 h of exposure. The method was further validated by successfully silencing other different genes, and it was also efficient toward two other aphid species, Aphis citricidus and Myzus persicae. Furthermore, a noticeable mortality was also observed in pea aphids using topical RNAi-mediated gene silencing, within 4 days post-dsRNA application for four out of seven tested genes.

CONCLUSION

Compared with the currently used dsRNA delivery methods in aphids, microinjection and ingestion, topical RNAi is time- and cost-effective, which could greatly influence RNAi-based gene functional studies and potential candidate gene selection for developing RNAi-based aphid control strategies in the future. © 2019 Society of Chemical Industry.

Comparative Analysis of Intra- and Inter-Specific Genomic Variability in the Peach Potato Aphid, Myzus persicae

The availability of genomic data in the last decade relating to different aphid species has allowed the analysis of the genomic variability occurring among such species, whereas intra-specific variability has hitherto very largely been neglected. In order to analyse the intra-genomic variability in the peach potato aphid, Myzus persicae, comparative analyses were performed revealing several clone-specific gene duplications, together with numerous deletions/rearrangements. Our comparative approach also allowed us to evaluate the synteny existing between the two M. persicae clones tested and between the peach potato aphid and the pea aphid, Acyrthosiphon pisum. Even if part of the observed rearrangements are related to a low quality of some assembled contigs and/or to the high number of contigs present in these aphid genomes, our evidence reveals that aphid clones are genetically more different than expected. These results suggest that the choice of performing genomes sequencing combining different biotypes/populations, as revealed in the case of the soybean aphid, Aphis glycines, is unlikely to be very informative in aphids. Interestingly, it is possible that the holocentric nature of aphid chromosomes favours genome rearrangements that can be successively inherited transgenerationally via the aphid's apomictic (parthenogenetic) mode of reproduction. Lastly, we evaluated the structure of the cluster of genes coding for the five histones (H1, H2A, H2B, H3 and H4) in order to better understand the quality of the two M. persicae genomes and thereby to improve our knowledge of this functionally important gene family.

Molecular Basis of Soybean Resistance to Soybean Aphids and Soybean Cyst Nematodes

Soybean aphid (SBA; Aphis glycines Matsumura) and soybean cyst nematode (SCN; Heterodera glycines Ichninohe) are major pests of the soybean (Glycine max [L.] Merr.). Substantial progress has been made in identifying the genetic basis of limiting these pests in both model and non-model plant systems. Classical linkage mapping and genome-wide association studies (GWAS) have identified major and minor quantitative trait loci (QTLs) in soybean. Studies on interactions of SBA and SCN effectors with host proteins have identified molecular cues in various signaling pathways, including those involved in plant disease resistance and phytohormone regulations. In this paper, we review the molecular basis of soybean resistance to SBA and SCN, and we provide a synthesis of recent studies of soybean QTLs/genes that could mitigate the effects of virulent SBA and SCN populations. We also review relevant studies of aphid–nematode interactions, particularly in the soybean–SBA–SCN system.

A chromosome-level draft genome of the grain aphid Sitobion miscanthi

BACKGROUND

Sitobion miscanthi is an ideal model for studying host plant specificity, parthenogenesis-based phenotypic plasticity, and interactions between insects and other species of various trophic levels, such as viruses, bacteria, plants, and natural enemies. However, the genome information for this species has not yet to be sequenced and published. Here, we analyzed the entire genome of a parthenogenetic female aphid colony using Pacific Biosciences long-read sequencing and Hi-C data to generate chromosome-length scaffolds and a highly contiguous genome assembly.

RESULTS

The final draft genome assembly from 33.88 Gb of raw data was ∼397.90 Mb in size, with a 2.05 Mb contig N50. Nine chromosomes were further assembled based on Hi-C data to a 377.19 Mb final size with a 36.26 Mb scaffold N50. The identified repeat sequences accounted for 26.41% of the genome, and 16,006 protein-coding genes were annotated. According to the phylogenetic analysis, S. miscanthi is closely related to Acyrthosiphon pisum, with S. miscanthi diverging from their common ancestor ∼25.0-44.9 million years ago.

CONCLUSIONS

We generated a high-quality draft of the S. miscanthi genome. This genome assembly should help promote research on the lifestyle and feeding specificity of aphids and their interactions with each other and species at other trophic levels. It can serve as a resource for accelerating genome-assisted improvements in insecticide-resistant management and environmentally safe aphid management.

A polymer/detergent formulation improves dsRNA penetration through the body wall and RNAi-induced mortality in the soybean aphid Aphis glycines

BACKGROUND

It is difficult to efficiently silence gene expression in some insects, probably because of the degradation of dsRNA by enzymes present in the gut and hemolymph post-oral feeding or body injecting of dsRNA. We previously developed a nanocarrier delivery system that can systemically deliver dsRNA into chewing mouthpart insects by oral feeding and efficiently silence gene expression. For the purpose of pest control in the field, there is a great demand to develop a spray method to apply dsRNA formulation.

RESULTS

A formulation of the nanocarrier/dsRNA/detergent was developed and could be easily applied just by dropping it on the notum of the aphid. The formulation penetrated the body wall into the hemocoel and then spread into various tissues within 1 h. The delivered dsRNA efficiently silenced the target gene expression at a high knockdown effect (95.4%) and the aphid population was largely suppressed (80.5%).

CONCLUSION

A novel dsRNA formulation was developed with the help of a nanocarrier and detergent that can quickly penetrate the insect body wall and efficiently silence gene expression. The formulation may provide a fast and easy tool for gene silence in some tough insects and for pest control in the field. © 2019 Society of Chemical Industry.

A Laterally Transferred Viral Gene Modifies Aphid Wing Plasticity

Organisms often respond to changing environments by altering development of particular traits. These plastic traits exhibit genetic variation; i.e., genotypes respond differently to the same environmental cues. Theoretical studies have demonstrated the importance of this variation, which is targeted by natural selection, in adapting plastic responses to maximize fitness [1, 2]. However, little is known about the underlying genetic mechanisms. We identify two laterally transferred genes that contribute to variation in a classic example of phenotypic plasticity: the pea aphid's ability to produce winged offspring in response to crowding. We discovered that aphid genotypes vary extensively for this trait and that aphid genes of viral origin are upregulated in response to crowding solely in highly inducible genotypes. We knocked down expression of these genes to demonstrate their functional role in wing plasticity. Through phylogenetic analysis, we found that these genes likely originated from a virus that infects rosy apple aphids and causes their hosts to produce winged offspring [3]. The function of these genes has therefore been retained following transfer to pea aphids. Our results uncover a novel role for co-opted viral genes, demonstrating that they are used to modulate ecologically relevant, plastic phenotypes. Our findings also address a critical question about the evolution of environmentally sensitive traits: whether the genes that control the expression of plastic traits also underlie variation in plasticity. The genes we identify originated from outside aphids themselves, and thus, our work shows that genes formerly unrelated to plasticity can fine-tune the strength of plastic responses to the environment.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.