Gene Family Evolution in the Pea Aphid Based on Chromosome-Level Genome Assembly

DNA methylation machinery is involved in development and reproduction in the viviparous pea aphid (Acyrthosiphon pisum)

Epigenetic mechanisms, such as DNA methylation, have been proposed to mediate plastic responses in insects. The pea aphid (Acyrthosiphon pisum), like the majority of extant aphids, displays cyclical parthenogenesis - the ability of mothers to switch the reproductive mode of their offspring from reproducing parthenogenetically to sexually in response to environmental cues. The pea aphid genome encodes two paralogs of the de novo DNA methyltransferase gene, dnmt3a and dnmt3x. Here we show, using phylogenetic analysis, that this gene duplication event occurred at least 150 million years ago, likely after the divergence of the lineage leading to the Aphidomorpha (phylloxerans, adelgids and true aphids) from that leading to the scale insects (Coccomorpha) and that the two paralogs are maintained in the genomes of all aphids examined. We also show that the mRNA of both dnmt3 paralogs is maternally expressed in the viviparous aphid ovary. During development both paralogs are expressed in the germ cells of embryos beginning at stage 5 and persisting throughout development. Treatment with 5-azactyidine, a chemical that generally inhibits the DNA methylation machinery, leads to defects of oocytes and early-stage embryos and causes a proportion of later stage embryos to be born dead or die soon after birth. These phenotypes suggest a role for DNA methyltransferases in reproduction, consistent with that seen in other insects. Taking the vast evolutionary history of the dnmt3 paralogs, and the localisation of their mRNAs in the ovary, we suggest there is a role for dnmt3a and/or dnmt3x in early development, and a role for DNA methylation machinery in reproduction and development of the viviparous pea aphid.

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.

A pipeline contributes to efficient identification of salivary proteins in short-headed planthopper, Epeurysa nawaii

Saliva, an oral secretion primarily originating from salivary glands (SGs), exert critical roles in the ongoing evolutionary interaction between insects and plants. However, identifying insect salivary components poses challenges due to the tiny size of insects, low secretion amounts, and the propensity for degradation after secretion. In this study, we developed a transcriptome-based approach to comprehensively analyze the salivary proteins of the short-headed planthopper, Epeurysa nawaii, a species with unique feeding habits on bamboo. A total of 165 salivary proteins were identified, with 114 secretory genes highly and specifically expressed in SGs. Consistent with most phloem-feeding insects, digestive enzymes, calcium-binding proteins, oxidoreductases, and a few previously reported salivary effectors were ubiquitously distributed in E. nawaii saliva. However, we also identified a substantial portion of salivary proteins exhibiting taxonomy specificity, including 60 E. nawaii-specific and 62 Delphacidae-specific proteins. These taxonomy-restricted proteins potentially play a role in insect adaptation to specific host plants. Our study provides an efficient pipeline for salivary protein identification and serves as a valuable resource for the functional characterization of effectors.

Transcriptomic basis of sex loss in the pea aphid

BACKGROUND

Transitions from sexual to asexual reproduction are common in eukaryotes, but the underlying mechanisms remain poorly known. The pea aphid-Acyrthosiphon pisum-exhibits reproductive polymorphism, with cyclical parthenogenetic and obligate parthenogenetic lineages, offering an opportunity to decipher the genetic basis of sex loss. Previous work on this species identified a single 840 kb region controlling reproductive polymorphism and carrying 32 genes. With the aim of identifying the gene(s) responsible for sex loss and the resulting consequences on the genetic programs controlling sexual or asexual embryogenesis, we compared the transcriptomic response to photoperiod shortening-the main sex-inducing cue-of a sexual and an obligate asexual lineage of the pea aphid, focusing on heads (where the photoperiodic cue is detected) and embryos (the final target of the cue).

RESULTS

Our analyses revealed that four genes (one expressed in the head, and three in the embryos) of the region responded differently to photoperiod in the two lineages. We also found that the downstream genetic programs expressed during embryonic development of a future sexual female encompass ∼1600 genes, among which miRNAs, piRNAs and histone modification pathways are overrepresented. These genes mainly co-localize in two genomic regions enriched in transposable elements (TEs).

CONCLUSIONS

Our results suggest that the causal polymorphism(s) in the 840 kb region somehow impair downstream epigenetic and post-transcriptional regulations in obligate asexual lineages, thereby sustaining asexual reproduction.

Genomic and Transcriptomic Insights into the Genetic Basis of Foam Secretion in Rice Spittlebug Callitettix versicolor

Many animal species produce protective foams, the majority of which exhibit evolutionary adaptability. Although the function and composition of foams have been widely studied, the genetic basis of foam secretion remains unknown. Unlike most species that produce foam under specific situations, spittlebugs continuously secrete foams throughout all nymphal stages. Here, we capitalize on the rice spittlebug (Callitettix versicolor) to explore the genetic basis of foam secretion through genomic and transcriptomic approaches. Our comparative genomic analysis for C. versicolor and eight other insect species reveals 606 species-specific gene families and 66 expanded gene families, associated with carbohydrate and lipid metabolism. These functions are in accordance with the composition of foams secreted by spittlebugs. Transcriptomic analyses of malpighian tubules across developmental stages detected 3192 differentially expressed genes. Enrichment analysis of these genes highlights functions also revealed by our comparative genomic analysis and aligns with previous histochemical and morphological observations of foam secretion. This consistency suggests the important roles of these candidate genes in foam production. Our study not only provides novel insights into the genetic basis of foam secretion in rice spittlebugs but also contributes valuable knowledge for future evolutionary studies of spittlebugs and the development of pest control strategies for C. versicolor.

Chromosome-level genome assembly of the Asian citrus psyllid, Diaphorina citri

Diaphorina citri is a global citrus pest. As a vector insect, it can transmit the causative agents of citrus huanglongbing, causing irreversible losses to the citrus industry. The acquisition of genomic information can provide a molecular genetic basis for effective control of D. citri. Here, the DNBSEQ™ , Oxford Nanopore Technologies, and Hi-C technologies are applied to generate a high-quality chromosome-level genome of D. citri. The genome size of D. citri was 523.78 Mb with a scaffold N50 of 47.05 Mb distributed on 13 chromosomes. A total of 250.64 Mb (47.85%) repeat sequences and 24 048 protein-coding genes were predicted. Genome resequencing of female and male individuals indicated that the sex chromosome system of D. citri is XO. Phylogenetic analysis demonstrated that D. citri and Pachypsylla venusta, which separated from their most recent common ancestor about 336.62 million years ago, were the most closely related. Additionally, we identified genes potentially involved in detoxification metabolism, pathogen transmission, and honeydew secretion for further investigation. The high-quality genome provides an important reference for developing effective management strategies of D. citri.

Chromosome-Level Genome Assembly and Annotation of a Periodical Cicada Species: Magicicada septendecula

We present a high-quality assembly and annotation of the periodical cicada species, Magicicada septendecula (Hemiptera: Auchenorrhyncha: Cicadidae). Periodical cicadas have a significant ecological impact, serving as a food source for many mammals, reptiles, and birds. Magicicada are well known for their massive emergences of 1 to 3 species that appear in different locations in the eastern United States nearly every year. These year classes ("broods") emerge dependably every 13 or 17 yr in a given location. Recently, it has become clear that 4-yr early or late emergences of a sizeable portion of a population are an important part of the history of brood formation; however, the biological mechanisms by which they track the passage of time remain a mystery. Using PacBio HiFi reads in conjunction with Hi-C proximity ligation data, we have assembled and annotated the first whole genome for a periodical cicada, an important resource for future phylogenetic and comparative genomic analysis. This also represents the first quality genome assembly and annotation for the Hemipteran superfamily Cicadoidea. With a scaffold N50 of 518.9 Mb and a complete BUSCO score of 96.7%, we are confident that this assembly will serve as a vital resource toward uncovering the genomic basis of periodical cicadas' long, synchronized life cycles and will provide a robust framework for further investigations into these insects.

Whole-genome sequence of the Cooley spruce gall adelgid, Adelges cooleyi (Hemiptera: Sternorrhyncha: Adelgidae)

The adelgids (Adelgidae) are a small family of sap-feeding insects, which, together with true aphids (Aphididae) and phylloxerans (Phylloxeridae), make up the infraorder Aphidomorpha. Some adelgid species are highly destructive to forest ecosystems such as Adelges tsugae, Adelges piceae, Adelges laricis, Pineus pini, and Pineus boerneri. Despite this, there are no high-quality genomic resources for adelgids, hindering advanced genomic analyses within Adelgidae and among Aphidomorpha. Here, we used PacBio continuous long-read and Illumina RNA-sequencing to construct a high-quality draft genome assembly for the Cooley spruce gall adelgid, Adelges cooleyi (Gillette), a gall-forming species endemic to North America. The assembled genome is 270.2 Mb in total size and has scaffold and contig N50 statistics of 14.87 and 7.18 Mb, respectively. There are 24,967 predicted coding sequences, and the assembly completeness is estimated at 98.1 and 99.6% with core BUSCO gene sets of Arthropoda and Hemiptera, respectively. Phylogenomic analysis using the A. cooleyi genome, 3 publicly available adelgid transcriptomes, 4 phylloxera transcriptomes, the Daktulosphaira vitifoliae (grape phylloxera) genome, 4 aphid genomes, and 2 outgroup coccoid genomes fully resolves adelgids and phylloxerans as sister taxa. The mitochondrial genome is 24 kb, among the largest in insects sampled to date, with 39.4% composed of noncoding regions. This genome assembly is currently the only genome-scale, annotated assembly for adelgids and will be a valuable resource for understanding the ecology and evolution of Aphidomorpha.

The nutritional dimension of facultative bacterial symbiosis in aphids: Current status and methodological considerations for future research

Aphids are valuable models for studying the functional diversity of bacterial symbiosis in insects. In addition to their ancestral obligate nutritional symbiont Buchnera aphidicola, these insects can host a myriad of so-called facultative symbionts. The diversity of these heritable bacterial associates is now well known, and some of the ecologically important traits associated with them have been well documented. Some twenty years ago, it was suggested that facultative symbionts could play an important role in aphid nutrition, notably by improving feeding performance on specific host plants, thus influencing the adaptation of these insects to host plants. However, the underlying mechanisms have never been elucidated, and the nutritional role that facultative symbionts might perform in aphids remains enigmatic. In this opinion piece, I put forward a series of arguments in support of the hypothesis that facultative symbionts play a central role in aphid nutrition and emphasize methodological considerations for testing this hypothesis in future work. In particular, I hypothesize that the metabolic capacities of B. aphidicola alone may not always be able to counterbalance the nutritional deficiencies of phloem sap. The association with one or several facultative symbionts with extensive metabolic capabilities would then be necessary to buffer the insect from host plant-derived nutrient deficiencies, thus enabling it to gain access to certain host plants.

Phylloxera and Aphids Show Distinct Features of Genome Evolution Despite Similar Reproductive Modes

Genomes of aphids (family Aphididae) show several unusual evolutionary patterns. In particular, within the XO sex determination system of aphids, the X chromosome exhibits a lower rate of interchromosomal rearrangements, fewer highly expressed genes, and faster evolution at nonsynonymous sites compared with the autosomes. In contrast, other hemipteran lineages have similar rates of interchromosomal rearrangement for autosomes and X chromosomes. One possible explanation for these differences is the aphid's life cycle of cyclical parthenogenesis, where multiple asexual generations alternate with 1 sexual generation. If true, we should see similar features in the genomes of Phylloxeridae, an outgroup of aphids which also undergoes cyclical parthenogenesis. To investigate this, we generated a chromosome-level assembly for the grape phylloxera, an agriculturally important species of Phylloxeridae, and identified its single X chromosome. We then performed synteny analysis using the phylloxerid genome and 30 high-quality genomes of aphids and other hemipteran species. Unexpectedly, we found that the phylloxera does not share aphids' patterns of chromosome evolution. By estimating interchromosomal rearrangement rates on an absolute time scale, we found that rates are elevated for aphid autosomes compared with their X chromosomes, but this pattern does not extend to the phylloxera branch. Potentially, the conservation of X chromosome gene content is due to selection on XO males that appear in the sexual generation. We also examined gene duplication patterns across Hemiptera and uncovered horizontal gene transfer events contributing to phylloxera evolution.

Horizontally Transferred Salivary Protein Promotes Insect Feeding by Suppressing Ferredoxin-Mediated Plant Defenses

Herbivorous insects such as whiteflies, planthoppers, and aphids secrete abundant orphan proteins to facilitate feeding. Yet, how these genes are recruited and evolve to mediate plant-insect interaction remains unknown. In this study, we report a horizontal gene transfer (HGT) event from fungi to an ancestor of Aleyrodidae insects approximately 42 to 190 million years ago. BtFTSP1 is a salivary protein that is secreted into host plants during Bemisia tabaci feeding. It targets a defensive ferredoxin 1 in Nicotiana tabacum (NtFD1) and disrupts the NtFD1-NtFD1 interaction in plant cytosol, leading to the degradation of NtFD1 in a ubiquitin-dependent manner. Silencing BtFTSP1 has negative effects on B. tabaci feeding while overexpressing BtFTSP1 in N. tabacum benefits insects and rescues the adverse effect caused by NtFD1 overexpression. The association between BtFTSP1 and NtFD1 is newly evolved after HGT, with the homologous FTSP in its fungal donor failing to interact and destabilize NtFD1. Our study illustrates the important roles of horizontally transferred genes in plant-insect interactions and suggests the potential origin of orphan salivary genes.

Parallel Losses of Blue Opsin Correlate with Compensatory Neofunctionalization of UV-Opsin Gene Duplicates in Aphids and Planthoppers

Expanding on previous efforts to survey the visual opsin repertoires of the Hemiptera, this study confirms that homologs of the UV- and LW-opsin subfamilies are conserved in all Hemiptera, while the B-opsin subfamily is missing from the Heteroptera and subgroups of the Sternorrhyncha and Auchenorrhyncha, i.e., aphids (Aphidoidea) and planthoppers (Fulgoroidea), respectively. Unlike in the Heteroptera, which are characterized by multiple independent expansions of the LW-opsin subfamily, the lack of B-opsin correlates with the presence of tandem-duplicated UV-opsins in aphids and planthoppers. Available data on organismal wavelength sensitivities and retinal gene expression patterns lead to the conclusion that, in both groups, one UV-opsin paralog shifted from ancestral UV peak sensitivity to derived blue sensitivity, likely compensating for the lost B-opsin. Two parallel bona fide tuning site substitutions compare to 18 non-corresponding amino acid replacements in the blue-shifted UV-opsin paralogs of aphids and planthoppers. Most notably, while the aphid blue-shifted UV-opsin clade is characterized by a replacement substitution at one of the best-documented UV/blue tuning sites (Rhodopsin site 90), the planthopper blue-shifted UV-opsin paralogs retained the ancestral lysine at this position. Combined, the new findings identify aphid and planthopper UV-opsins as a new valuable data sample for studying adaptive opsin evolution.

Contrasting Evolutionary Patterns Between Sexual and Asexual Lineages in a Genomic Region Linked to Reproductive Mode Variation in the pea aphid

Although asexual lineages evolved from sexual lineages in many different taxa, the genetics of sex loss remains poorly understood. We addressed this issue in the pea aphid Acyrthosiphon pisum, whose natural populations encompass lineages performing cyclical parthenogenesis (CP) and producing one sexual generation per year, as well as obligate parthenogenetic (OP) lineages that can no longer produce sexual females but can still produce males. An SNP-based, whole-genome scan of CP and OP populations sequenced in pools (103 individuals from 6 populations) revealed that an X-linked region is associated with the variation in reproductive mode. This 840-kb region is highly divergent between CP and OP populations (FST = 34.9%), with >2,000 SNPs or short Indels showing a high degree of association with the phenotypic trait. In OP populations specifically, this region also shows reduced diversity and Tajima's D, consistent with the OP phenotype being a derived trait in aphids. Interestingly, the low genetic differentiation between CP and OP populations at the rest of the genome (FST = 2.5%) suggests gene flow between them. Males from OP lineages thus likely transmit their op allele to new genomic backgrounds. These genetic exchanges, combined with the selection of the OP and CP reproductive modes under different climates, probably contribute to the long-term persistence of the cp and op alleles.

Phylloxera and aphids show distinct features of genome evolution despite similar reproductive modes

Genomes of aphids (family Aphididae) show several unusual evolutionary patterns. In particular, within the XO sex determination system of aphids, the X chromosome exhibits a lower rate of interchromosomal rearrangements, fewer highly expressed genes, and faster evolution at nonsynonymous sites compared to the autosomes. In contrast, other hemipteran lineages have similar rates of interchromosomal rearrangement for autosomes and X chromosomes. One possible explanation for these differences is the aphid's life cycle of cyclical parthenogenesis, where multiple asexual generations alternate with one sexual generation. If true, we should see similar features in the genomes of Phylloxeridae, an outgroup of aphids which also undergoes cyclical parthenogenesis. To investigate this, we generated a chromosome-level assembly for the grape phylloxera, an agriculturally important species of Phylloxeridae, and identified its single X chromosome. We then performed synteny analysis using the phylloxerid genome and 30 high-quality genomes of aphids and other hemipteran species. Unexpectedly, we found that the phylloxera does not share aphids' patterns of chromosome evolution. By estimating interchromosomal rearrangement rates on an absolute time scale, we found that rates are elevated for aphid autosomes compared to their X chromosomes, but this pattern does not extend to the phylloxera branch. Potentially, the conservation of X chromosome gene content is due to selection on XO males that appear in the sexual generation. We also examined gene duplication patterns across Hemiptera and uncovered horizontal gene transfer events contributing to phylloxera evolution.

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.

High-quality chromosome-level scaffolds of the plant bug Pachypeltis micranthus provide insights into the availability of Mikania micrantha control

BACKGROUND

The plant bug, Pachypeltis micranthus Mu et Liu (Hemiptera: Miridae), is an effective potential biological control agent for Mikania micrantha H.B.K. (Asteraceae; one of the most notorious invasive weeds worldwide). However, limited knowledge about this species hindered its practical application and research. Accordingly, sequencing the genome of this mirid bug holds great significance in controlling M. micrantha.

RESULTS

Here, 712.72 Mb high-quality chromosome-level scaffolds of P. micranthus were generated, of which 707.51 Mb (99.27%) of assembled sequences were anchored onto 15 chromosome-level scaffolds with contig N50 of 16.84 Mb. The P. micranthus genome had the highest GC content (42.43%) and the second highest proportion of repetitive sequences (375.82 Mb, 52.73%) than the three other mirid bugs (i.e., Apolygus lucorum, Cyrtorhinus lividipennis, and Nesidiocoris tenuis). Phylogenetic analysis showed that P. micranthus clustered with other mirid bugs and diverged from the common ancestor approximately 200 million years ago. Gene family expansion and/or contraction were analyzed, and significantly expanded gene families associated with P. micranthus feeding and adaptation to M. micrantha were manually identified. Compared with the whole body, transcriptome analysis of the salivary gland revealed that most of the upregulated genes were significantly associated with metabolism pathways and peptidase activity, particularly among cysteine peptidase, serine peptidase, and polygalacturonase; this could be one of the reasons for precisely and highly efficient feeding by the oligophagous bug P. micranthus on M. micrantha.

CONCLUSION

Collectively, this work provides a crucial chromosome-level scaffolds resource to study the evolutionary adaptation between mirid bug and their host. It is also helpful in searching for novel environment-friendly biological strategies to control M. micrantha.

Pigment-dispersing factor is present in circadian clock neurons of pea aphids and may mediate photoperiodic signalling to insulin-producing cells

The neuropeptide pigment-dispersing factor (PDF) plays a pivotal role in the circadian clock of most Ecdysozoa and is additionally involved in the timing of seasonal responses of several photoperiodic species. The pea aphid, Acyrthosiphon pisum, is a paradigmatic photoperiodic species with an annual life cycle tightly coupled to the seasonal changes in day length. Nevertheless, PDF could not be identified in A. pisum so far. In the present study, we identified a PDF-coding gene that has undergone significant changes in the otherwise highly conserved insect C-terminal amino acid sequence. A newly generated aphid-specific PDF antibody stained four neurons in each hemisphere of the aphid brain that co-express the clock protein Period and have projections to the pars lateralis that are highly plastic and change their appearance in a daily and seasonal manner, resembling those of the fruit fly PDF neurons. Most intriguingly, the PDF terminals overlap with dendrites of the insulin-like peptide (ILP) positive neurosecretory cells in the pars intercerebralis and with putative terminals of Cryptochrome (CRY) positive clock neurons. Since ILP has been previously shown to be crucial for seasonal adaptations and CRY might serve as a circadian photoreceptor vital for measuring day length, our results suggest that PDF plays a critical role in aphid seasonal timing.

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.

Comparative transcriptomics of aphid species that diverged > 22 MYA reveals genes that are important for the maintenance of their symbiosis

Most plant-sap feeding insects have obligate relationships with maternally transmitted bacteria. Aphids require their nutritional endosymbiont, Buchnera aphidicola, for the production of essential amino acids. Such endosymbionts are harbored inside of specialized insect cells called bacteriocytes. Here, we use comparative transcriptomics of bacteriocytes between two recently diverged aphid species, Myzus persicae and Acyrthosiphon pisum, to identify key genes that are important for the maintenance of their nutritional mutualism. The majority of genes with conserved expression profiles in M. persicae and A. pisum are for orthologs previously identified in A. pisum to be important for the symbiosis. However, asparaginase which produces aspartate from asparagine was significantly up-regulated only in A. pisum bacteriocytes, potentially because Buchnera of M. persicae encodes its own asparaginase enzyme unlike Buchnera of A. pisum resulting in Buchnera of A. pisum to be dependent on its aphid host for aspartate. One-to-one orthologs that explained the most amount of variation for bacteriocyte specific mRNA expression for both species includes a collaborative gene for methionine biosynthesis, multiple transporters, a horizontally transmitted gene, and secreted proteins. Finally, we highlight species-specific gene clusters which may contribute to host adaptations and/or accommodations in gene regulation to changes in the symbiont or the symbiosis.

Planthopper salivary sheath protein LsSP1 contributes to manipulation of rice plant defenses

Salivary elicitors secreted by herbivorous insects can be perceived by host plants to trigger plant immunity. However, how insects secrete other salivary components to subsequently attenuate the elicitor-induced plant immunity remains poorly understood. Here, we study the small brown planthopper, Laodelphax striatellus salivary sheath protein LsSP1. Using Y2H, BiFC and LUC assays, we show that LsSP1 is secreted into host plants and binds to salivary sheath via mucin-like protein (LsMLP). Rice plants pre-infested with dsLsSP1-treated L. striatellus are less attractive to L. striatellus nymphs than those pre-infected with dsGFP-treated controls. Transgenic rice plants with LsSP1 overexpression rescue the insect feeding defects caused by a deficiency of LsSP1 secretion, consistent with the potential role of LsSP1 in manipulating plant defenses. Our results illustrate the importance of salivary sheath proteins in mediating the interactions between plants and herbivorous insects.

Differential regulation of the Tor gene homolog drives the red/green pigmentation phenotype in the aphid Myzuspersicae

In some aphid species, intraspecific variation in body colour is caused by differential carotenoid content: whilst green aphids contain only yellow carotenoids (β-, γ-, and β,γ-carotenes), red aphids additionally possess red carotenoids (torulene and 3,4-didehydrolycopene). Unusually, within animals who typically obtain carotenoids from their diet, ancestral horizontal gene transfer of carotenoid biosynthetic genes from fungi (followed by gene duplication), have imbued aphids with the intrinsic gene repertoire necessary to biosynthesise carotenoids. In the pea aphid, Acyrthosiphon pisum a lycopene (phytoene) desaturase gene (Tor) underpins the red/green phenotype, with this locus present in heterozygous form in red individuals but absent in green aphids, resulting in them being unable to convert lycopene into the red compounds 3,4-didehydrolycopene and torulene. The green peach aphid, Myzus persicae, separated from the pea aphid for ≈45MY also exists as distinct colour variable morphs, with both red and green individuals present. Here, we examined genomic data for both red and green morphs of M. persicae and identified an enlarged (compared to A. pisum) repertoire of 16 carotenoid biosynthetic genes (11 carotenoid desaturases and five carotenoid cyclase/synthase genes). From these, we identify the homolog of A. pisum Tor (here called carotene desaturase 2 or CDE-2) and show through 3D modelling that this homolog can accommodate the torulene precursor lycopene and, through RNA knockdown feeding experiments, demonstrate that disabling CDE-2 expression in red M. persicae clones results in green-coloured offspring. Unlike in A. pisum, we show that functional CDE-2 is present in the genomes of both red and green aphids. However, expression differences between the two colour morphs (350-700 fold CDE-2 overexpression in red clones), potentially driven by variants identified in upstream putative regulatory elements, underpin this phenotype. Thus, whilst aphids have a common origin of their carotenoid biosynthetic pathway, two aphid species separated for over 40MY have evolved very different drivers of intraspecific colour variation.

Addressing the challenges of symbiont-mediated RNAi in aphids

Because aphids are global agricultural pests and models for bacterial endosymbiosis, there is a need for reliable methods to study and control their gene function. However, current methods available for aphid gene knockout and knockdown of gene expression are often unreliable and time consuming. Techniques like CRISPR-Cas genome editing can take several months to achieve a single gene knockout because they rely on aphids going through a cycle of sexual reproduction, and aphids often lack strong, consistent levels of knockdown when fed or injected with molecules that induce an RNA interference (RNAi) response. In the hopes of addressing these challenges, we attempted to adapt a new method called symbiont-mediated RNAi (smRNAi) for use in aphids. smRNAi involves engineering a bacterial symbiont of the insect to continuously supply double-stranded RNA (dsRNA) inside the insect body. This approach has been successful in thrips, kissing bugs, and honeybees. We engineered the laboratory Escherichia coli strain HT115 and the native aphid symbiont Serratia symbiotica CWBI-2.3^T to produce dsRNA inside the gut of the pea aphid (Acyrthosiphon pisum) targeting salivary effector protein (C002) or ecdysone receptor genes. For C002 assays, we also tested co-knockdown with an aphid nuclease (Nuc1) to reduce RNA degradation. However, we found that smRNAi was not a reliable method for aphid gene knockdown under our conditions. We were unable to consistently achieve the expected phenotypic changes with either target. However, we did see indications that elements of the RNAi pathway were modestly upregulated, and expression of some targeted genes appeared to be somewhat reduced in some trials. We conclude with a discussion of the possible avenues through which smRNAi, and aphid RNAi in general, could be improved in the future.

Genomic Variation in the Tea Leafhopper Reveals the Basis of Adaptive Evolution

Tea green leafhopper (TGL), Empoasca onukii, is of biological and economic interest. Despite numerous studies, the mechanisms underlying its adaptation and evolution remain enigmatic. Here, we use previously untapped genome and population genetics approaches to examine how the pest adapted to different environmental variables and thus has expanded geographically. We complete a chromosome-level assembly and annotation of the E. onukii genome, showing notable expansions of gene families associated with adaptation to chemoreception and detoxification. Genomic signals indicating balancing selection highlight metabolic pathways involved in adaptation to a wide range of tea varieties grown across ecologically diverse regions. Patterns of genetic variations among 54 E. onukii samples unveil the population structure and evolutionary history across different tea-growing regions in China. Our results demonstrate that the genomic changes in key pathways, including those linked to metabolism, circadian rhythms, and immune system functions, may underlie the successful spread and adaptation of E. onukii. This work highlights the genetic and molecular bases underlying the evolutionary success of a species with broad economic impacts, and provides insights into insect adaptation to host plants, which will ultimately facilitate more sustainable pest management.

The genomic basis of evolutionary novelties in a leafhopper

Evolutionary innovations generate phenotypic and species diversity. Elucidating the genomic processes underlying such innovations is central to understanding biodiversity. In this study, we addressed the genomic basis of evolutionary novelties in the glassy-winged sharpshooter (Homalodisca vitripennis, GWSS), an agricultural pest. Prominent evolutionary innovations in leafhoppers include brochosomes, proteinaceous structures that are excreted and used to coat the body, and obligate symbiotic associations with two bacterial types that reside within cytoplasm of distinctive cell types. Using PacBio long-read sequencing and Dovetail Omni-C technology, we generated a chromosome-level genome assembly for the GWSS and then validated the assembly using flow cytometry and karyotyping. Additional transcriptomic and proteomic data were used to identify novel genes that underlie brochosome production. We found that brochosome-associated genes include novel gene families that have diversified through tandem duplications. We also identified the locations of genes involved in interactions with bacterial symbionts. Ancestors of the GWSS acquired bacterial genes through horizontal gene transfer (HGT), and these genes appear to contribute to symbiont support. Using a phylogenomics approach, we inferred HGT sources and timing. We found that some HGT events date to the common ancestor of the hemipteran suborder Auchenorrhyncha, representing some of the oldest known examples of HGT in animals. Overall, we show that evolutionary novelties in leafhoppers are generated by the combination of acquiring novel genes, produced both de novo and through tandem duplication, acquiring new symbiotic associations that enable use of novel diets and niches, and recruiting foreign genes to support symbionts and enhance herbivory.

Chromosomal-level assembly of Bactericera cockerelli reveals rampant gene family expansions impacting genome structure, function and insect-microbe-plant-interactions

Lineage specific expansions and gene duplications are some of the most important sources of evolutionary novelty in eukaryotes. Although not as prevalent in eukaryotes compared to bacteria, horizontal gene transfer events can also result in key adaptations for insects, especially for those involved in insect-microbe interactions. In this study we assemble the first chromosomal assembly of the psyllid Bactericera cockerelli and reveal that the B. cockerelli genome has experienced significantly more gene expansion events compared to other Hemipteran representatives with fully sequenced genomes. We also reveal that B. cockerelli's genome is the largest psyllid genome (567 Mb) sequenced to date and is ~15% larger than the other two psyllid species genomes sequenced (Pachypsylla venusta and Diaphorina citri). Structurally, B. cockerelli appears to have an additional chromosome compared to the distantly related psyllid species P. venusta due to a previous chromosomal fission or fusion event. The increase in genome size and dynamic nature of the B. cockerelli genome may largely be contributed to the widespread expansion of type I and type II repeat elements that are rampant across all of B. cockerelli's. chromosomes. These repeat elements are distributed near equally in both euchromatic and heterochromatic regions. Furthermore, significant gene family expansions and gene duplications were uncovered for genes that are expected to be important in its adaptation to insect-plant and microbe interactions, which include transcription factors, proteases, odorant receptors, and horizontally transferred genes that are involved in the nutritional symbioses with their long-term nutritional endosymbiont Carsonella.

Gene Editing and Genetic Control of Hemipteran Pests: Progress, Challenges and Perspectives

The origin of the order Hemiptera can be traced to the late Permian Period more than 230 MYA, well before the origin of flowering plants 100 MY later in during the Cretaceous period. Hemipteran species consume their liquid diets using a sucking proboscis; for phytophagous hemipterans their mouthparts (stylets) are elegant structures that enable voracious feeding from plant xylem or phloem. This adaptation has resulted in some hemipteran species becoming globally significant pests of agriculture resulting in significant annual crop losses. Due to the reliance on chemical insecticides for the control of insect pests in agricultural settings, many hemipteran pests have evolved resistance to insecticides resulting in an urgent need to develop new, species-specific and environmentally friendly methods of pest control. The rapid advances in CRISPR/Cas9 technologies in model insects such as Drosophila melanogaster, Tribolium castaneum, Bombyx mori, and Aedes aegypti has spurred a new round of innovative genetic control strategies in the Diptera and Lepidoptera and an increased interest in assessing genetic control technologies for the Hemiptera. Genetic control approaches in the Hemiptera have, to date, been largely overlooked due to the problems of introducing genetic material into the germline of these insects. The high frequency of CRISPR-mediated mutagenesis in model insect species suggest that, if the delivery problem for Hemiptera could be solved, then gene editing in the Hemiptera might be quickly achieved. Significant advances in CRISPR/Cas9 editing have been realized in nine species of Hemiptera over the past 4 years. Here we review progress in the Hemiptera and discuss the challenges and opportunities for extending contemporary genetic control strategies into species in this agriculturally important insect orderr.

Ploidy dynamics in aphid host cells harboring bacterial symbionts

Aphids have evolved bacteriocytes or symbiotic host cells that harbor the obligate mutualistic bacterium Buchnera aphidicola. Because of the large cell size (approximately 100 μm in diameter) of bacteriocytes and their pivotal role in nutritional symbiosis, researchers have considered that these cells are highly polyploid and assumed that bacteriocyte polyploidy may be essential for the symbiotic relationship between the aphid and the bacterium. However, little is known about the ploidy levels and dynamics of aphid bacteriocytes. Here, we quantitatively analyzed the ploidy levels in the bacteriocytes of the pea-aphid Acyrthosiphon pisum. Image-based fluorometry revealed the hyper polyploidy of the bacteriocytes ranging from 16- to 256-ploidy throughout the lifecycle. Bacteriocytes of adult parthenogenetic viviparous females were ranged between 64 and 128C DNA levels, while those of sexual morphs (oviparous females and males) were comprised of 64C, and 32–64C cells, respectively. During post-embryonic development of viviparous females, the ploidy level of bacteriocytes increased substantially, from 16 to 32C at birth to 128–256C in actively reproducing adults. These results suggest that the ploidy levels are dynamically regulated among phenotypes and during development. Our comprehensive and quantitative data provides a foundation for future studies to understand the functional roles and biological significance of the polyploidy of insect bacteriocytes.

Characterization of carotenoid biosynthetic pathway genes in the pea aphid (Acyrthosiphon pisum) revealed by heterologous complementation and RNA interference assays

Carotenoids are involved in many essential physiological functions and are produced from geranylgeranyl pyrophosphate through synthase, desaturase, and cyclase activities. In the pea aphid (Acyrthosiphon pisum), the duplication of carotenoid biosynthetic genes, including carotenoid synthases/cyclases (ApCscA-C) and desaturases (ApCdeA-D), through horizontal gene transfer from fungi has been detected, and ApCdeB has known dehydrogenation functions. However, whether other genes contribute to aphid carotenoid biosynthesis, and its specific regulatory pathway, remains unclear. In the current study, functional analyses of seven genes were performed using heterologous complementation and RNA interference assays. The bifunctional enzymes ApCscA-C were responsible for the synthase of phytoene, and ApCscC may also have a cyclase activity. ApCdeA, ApCdeC, and ApCdeD had diverse dehydrogenation functions. ApCdeA catalyzed the enzymatic conversion of phytoene to neurosporene (three-step product), ApCdeC catalyzed the enzymatic conversion of phytoene to ζ-carotene (two-step product), and ApCdeD catalyzed the enzymatic conversion of phytoene to lycopene (four-step product). Silencing of ApCscs reduced the expression levels of ApCdes, and silencing these carotenoid biosynthetic genes reduced the α-, β-, and γ-carotene levels, as well as the total carotenoid level. The results suggest that these genes were activated and led to carotenoid biosynthesis in the pea aphid.

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.

Annotation of Hox cluster and Hox cofactor genes in the Asian citrus psyllid, Diaphorina citri, reveals novel features

Hox genes and their cofactors are essential developmental genes specifying regional identity in animals. Hox genes have a conserved arrangement in clusters in the same order in which they specify identity along the anterior-posterior axis. A few insect species have breaks in the cluster, but these are exceptions. We annotated the 10 Hox genes of the Asian citrus psyllid Diaphorina citri, and found a split in its Hox cluster between the Deformed and Sex combs reduced genes - the first time a break at this position has been observed in an insect Hox cluster. We also annotated D. citri orthologs of the Hox cofactor genes homothorax, PKNOX and extradenticle and found an additional copy of extradenticle in D. citri that appears to be a retrogene. Expression data and sequence conservation suggest that the extradenticle retrogene may have retained the original extradenticle function and allowed divergence of the parental extradenticle gene.

Chromosome-level genome of the globe skimmer dragonfly (Pantala flavescens)

BACKGROUND

The globe skimmer dragonfly (Pantala flavescens) is a notable Odonata insect distributed in nature fields and farmlands worldwide, and it is commonly recognized as a natural enemy because it preys on agricultural pests and health pests. As one of the sister groups of winged insects, odonatan species are key to understanding the evolution of insect wings.

FINDINGS

We present a high-quality reference genome of P. flavescens, which is the first chromosome-level genome in the Palaeoptera (Odonata and Ephemeroptera). The assembled genome size was 662 Mb, with a contig N50 of 16.2 Mb. Via Hi-C scaffolding, 648 Mb (97.9%) of contig sequences were clustered, ordered, and assembled into 12 large scaffolds, each corresponding to a natural chromosome. The X chromosome was identified by sequence coverage depth. The repetitive sequences and gene density of the X chromosome are similar to those of autosomal sequences, but the X chromosome shows a much lower degree of heterozygosity. Our analysis shows that the effective population size experienced 3 declining events, which may have been caused by climate change and environmental pollution.

CONCLUSIONS

The genome of P. flavescens provides more information on the biology and evolution of insects and will help for the use of this species in pest control.

Genomic identification, annotation, and comparative analysis of Vacuolar-type ATP synthase subunits in Diaphorina citri

The hemipteran insect Diaphorina citri, or Asian citrus psyllid, is a vector for Candidatus Liberibacter asiaticus (CLas), the bacterium causing citrus greening disease, or Huanglongbing (HLB). Millions of citrus trees have been destroyed, and every grove in Florida, USA, has been directly affected by this disease. In eukaryotes, vacuolar-type ATP synthase (V-ATPase) is an abundant heterodimeric enzyme that serves the cell with essential compartment acidification through the active processes that transport protons across the membrane. Fifteen putative V-ATPase genes in the D. citri genome were manually curated. Comparative genomic analysis revealed that D. citri V-ATPase subunits share domains and motifs with other insects, including the V-ATPase-A superfamily domain. Phylogenetic analysis separates D. citri V-ATPase subunits into expected clades with orthologous sequences. Annotation of the D. citri genome is a critical step towards developing directed pest management strategies to reduce the spread of HLB throughout the citrus industry.

Evolutionary diversification of insulin-related peptides (IRPs) in aphids and spatiotemporal distribution in Acyrthosiphon pisum

Members of the insulin superfamily activate the evolutionarily highly conserved insulin/insulin-like growth factor signaling pathway, involved in regulation of growth, energy homeostasis, and longevity. In the current study we focus on aphids to gain more insight into the evolution of the IRPs and how they may contribute to regulation of the insulin-signaling pathway. Using the latest annotation of the pea aphid (Acyrthosiphon pisum) genome, and combining sequence alignments and phylogenetic analyses, we identified seven putative IRP encoding-genes, with IRP1-IRP4 resembling the classical insulin and insulin-like protein structures, and IRP5 and IRP6 bearing insulin-like growth factor (IGF) features. We also identified IRP11 as a new and structurally divergent IRP present in at least eight aphid genomes. Globally the ten aphid genomes analyzed in this work contain four to 15 IRPs, while only three IRPs were found in the genome of the grape phylloxera, a hemipteran insect representing an earlier evolutionary branch of the aphid group. Expression analyses revealed spatial and temporal variation in the expression patterns of the different A. pisum IRPs. IRP1 and IRP4 are expressed throughout all developmental stages and morphs in neuroendocrine cells of the brain, while IRP5 and IRP6 are expressed in the fat body. IRP2 is expressed in specific cells of the gut in aphids in non-crowded conditions and in the head of aphids under crowded conditions, IRP3 in salivary glands, and both IRP2 and IRP3 in the male morph. IRP11 expression is enriched in the carcass. This complex spatiotemporal expression pattern suggests functional diversification of the IRPs.

OUP accepted manuscript

Homology of highly divergent genes often cannot be determined from sequence similarity alone. For example, we recently identified in the aphid Hormaphis cornu a family of rapidly evolving bicycle genes, which encode novel proteins implicated as plant gall effectors, and sequence similarity search methods yielded few putative bicycle homologs in other species. Coding sequence-independent features of genes, such as intron-exon boundaries, often evolve more slowly than coding sequences, however, and can provide complementary evidence for homology. We found that a linear logistic regression classifier using only structural features of bicycle genes identified many putative bicycle homologs in other species. Independent evidence from sequence features and intron locations supported homology assignments. To test the potential roles of bicycle genes in other aphids, we sequenced the genome of a second gall-forming aphid, Tetraneura nigriabdominalis and found that many bicycle genes are strongly expressed in the salivary glands of the gall forming foundress. In addition, bicycle genes are strongly overexpressed in the salivary glands of a non-gall forming aphid, Acyrthosiphon pisum, and in the non-gall forming generations of H. cornu. These observations suggest that Bicycle proteins may be used by multiple aphid species to manipulate plants in diverse ways. Incorporation of gene structural features into sequence search algorithms may aid identification of deeply divergent homologs, especially of rapidly evolving genes involved in host-parasite interactions.

BUSCO: Assessing Genomic Data Quality and Beyond

Evaluation of the quality of genomic "data products" such as genome assemblies or gene sets is of critical importance in order to recognize possible issues and correct them during the generation of new data. It is equally essential to guide subsequent or comparative analyses with existing data, as the correct interpretation of the results necessarily requires knowledge about the quality level and reliability of the inputs. Using datasets of near universal single-copy orthologs derived from OrthoDB, BUSCO can estimate the completeness and redundancy of genomic data by providing biologically meaningful metrics based on expected gene content. These can complement technical metrics such as contiguity measures (e.g., number of contigs/scaffolds, and N50 values). Here, we describe the use of the BUSCO tool suite to assess different data types that can range from genome assemblies of single isolates and assembled transcriptomes and annotated gene sets to metagenome-assembled genomes where the taxonomic origin of the species is unknown. BUSCO is the only tool capable of assessing all these types of sequences from both eukaryotic and prokaryotic species. The protocols detail the various BUSCO running modes and the novel workflows introduced in versions 4 and 5, including the batch analysis on multiple inputs, the auto-lineage workflow to run assessments without specifying a dataset, and a workflow for the evaluation of (large) eukaryotic genomes. The protocols further cover the BUSCO setup, guidelines to interpret the results, and BUSCO "plugin" workflows for performing common operations in genomics using BUSCO results, such as building phylogenomic trees and visualizing syntenies. © 2021 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Assessing an input sequence with a BUSCO dataset specified manually Basic Protocol 2: Assessing an input sequence with a dataset automatically selected by BUSCO Basic Protocol 3: Assessing multiple inputs Alternate Protocol: Decreasing analysis runtime when assessing a large number of small genomes with BUSCO auto-lineage workflow and Snakemake Support Protocol 1: BUSCO setup Support Protocol 2: Visualizing BUSCO results Support Protocol 3: Building phylogenomic trees.

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.

The boom and bust of the aphid's essential amino acid metabolism across nymphal development

Within long-term symbioses, animals integrate their physiology and development with their symbiont. In a model nutritional mutualism, aphids harbor the endosymbiont, Buchnera, within specialized bacteriocyte cells. Buchnera synthesizes essential amino acids (EAAs) and vitamins for their host, which are lacking from the aphid's plant sap diet. It is unclear if the aphid host differentially expresses aphid EAA metabolism pathways and genes that collaborate with Buchnera for the production of EAA and vitamins throughout nymphal development when feeding on plants. It is also unclear if aphid bacteriocytes are differentially methylated throughout aphid development as DNA methylation may play a role in gene regulation. By analyzing aphid gene expression, we determined that the bacteriocyte is metabolically more active in metabolizing Buchnera's EAAs and vitamins early in nymphal development compared to intermediate or later immature and adult lifestages. The largest changes in aphid bacteriocyte gene expression, especially for aphid genes that collaborate with Buchnera, occurred during the 3rd to 4th instar transition. During this transition, there is a huge shift in the bacteriocyte from a high energy "nutrient-consuming state" to a "recovery and growth state" where patterning and signaling genes and pathways are upregulated and differentially methylated, and de novo methylation is reduced as evidenced by homogenous DNA methylation profiles after the 2nd instar. Moreover, bacteriocyte number increased and Buchnera's titer decreased throughout aphid nymphal development. These data suggest in combination that bacteriocytes of older nymphal and adult lifestages depend less on the nutritional symbiosis compared to early nymphal lifestages.

Genetic innovations in animal-microbe symbioses

Animal hosts have initiated myriad symbiotic associations with microorganisms and often have maintained these symbioses for millions of years, spanning drastic changes in ecological conditions and lifestyles. The establishment and persistence of these relationships require genetic innovations on the parts of both symbionts and hosts. The nature of symbiont innovations depends on their genetic population structure, categorized here as open, closed or mixed. These categories reflect modes of inter-host transmission that result in distinct genomic features, or genomic syndromes, in symbionts. Although less studied, hosts also innovate in order to preserve and control symbiotic partnerships. New capabilities to sequence host-associated microbial communities and to experimentally manipulate both hosts and symbionts are providing unprecedented insights into how genetic innovations arise under different symbiont population structures and how these innovations function to support symbiotic relationships.

Global patterns in genomic diversity underpinning the evolution of insecticide resistance in the aphid crop pest Myzus persicae

The aphid Myzus persicae is a destructive agricultural pest that displays an exceptional ability to develop resistance to both natural and synthetic insecticides. To investigate the evolution of resistance in this species we generated a chromosome-scale genome assembly and living panel of >110 fully sequenced globally sampled clonal lines. Our analyses reveal a remarkable diversity of resistance mutations segregating in global populations of M. persicae. We show that the emergence and spread of these mechanisms is influenced by host-plant associations, uncovering the widespread co-option of a host-plant adaptation that also offers resistance against synthetic insecticides. We identify both the repeated evolution of independent resistance mutations at the same locus, and multiple instances of the evolution of novel resistance mechanisms against key insecticides. Our findings provide fundamental insights into the genomic responses of global insect populations to strong selective forces, and hold practical relevance for the control of pests and parasites.

Male-Biased microRNA Discovery in the Pea Aphid

Epigenetic mechanisms modulate gene expression levels during development, shaping how a single genome produces a diversity of phenotypes. Here, we begin to explore the epigenetic regulation of sexual dimorphism in pea aphids (Acyrthosiphon pisum) by focusing on microRNAs. Previous analyses of microRNAs in aphids have focused solely on females, so we performed deep sequencing of a sample containing early-stage males. We used this sample, plus samples from Genbank, to find 207 novel pea aphid microRNA coding loci. We localized microRNA loci to a chromosome-level assembly of the pea aphid genome and found that those on the X chromosome have lower overall expression compared to those on autosomes. We then identified a set of 19 putative male-biased microRNAs and found them enriched on the X chromosome. Finally, we performed protein-coding RNA-Seq of first instar female and male pea aphids to identify genes with lower expression in males. 10 of these genes were predicted targets of the 19 male-biased microRNAs. Our study provides the most complete set of microRNAs in the pea aphid to date and serves as foundational work for future studies on the epigenetic control of sexual dimorphism.

Chromosome-level genome assembly of the bean bug Riptortus pedestris

The bean bug (Riptortus pedestris) causes great economic losses of soybeans by piercing and sucking pods and seeds. Although R. pedestris has become the focus of numerous studies associated with insect-microbe interactions, plant-insect interactions, and pesticide resistance, a lack of genomic resources has limited deeper insights. Here, we report the first R. pedestris genome at the chromosomal level using PacBio, Illumina, and Hi-C technologies. The assembled genome was 1.080 Gb in size with a contig N50 of 2.882 Mb. More than 96.3% of the total genome bases were successfully anchored to six unique chromosomes. Genome resequencing of male and female individuals and chromosomic staining demonstrated that the sex chromosome system of R. pedestris is XO, and the shortest chromosome is the X chromosome. In total, 19,026 protein-coding genes were predicted, 18,745 of which were validated as being expressed. Temporospatial expression of R. pedestris genes in six tissues and 37 development stages revealed 4,657 and 7,793 genes mainly expressed in gonads and egg periods, respectively. Evolutionary analysis demonstrated that R. pedestris and Oncopeltus fasciatus formed a sister group and split ∼80 million years ago (Mya). Additionally, a 5.04 Mb complete genome of symbiotic Serratia marcescens Rip1 was assembled, and the virulence factors that account for successful colonization in the host midgut were identified. The high-quality R. pedestris genome provides a valuable resource for further research, as well as for the pest management of bug pests.

Selection following Gene Duplication Shapes Recent Genome Evolution in the Pea Aphid Acyrthosiphon pisum

Ecology of insects is as wide as their diversity, which reflects their high capacity of adaptation in most of the environments of our planet. Aphids, with over 4,000 species, have developed a series of adaptations including a high phenotypic plasticity and the ability to feed on the phloem sap of plants, which is enriched in sugars derived from photosynthesis. Recent analyses of aphid genomes have indicated a high level of shared ancestral gene duplications that might represent a basis for genetic innovation and broad adaptations. In addition, there are a large number of recent, species-specific gene duplications whose role in adaptation remains poorly understood. Here, we tested whether duplicates specific to the pea aphid Acyrthosiphon pisum are related to genomic innovation by combining comparative genomics, transcriptomics, and chromatin accessibility analyses. Consistent with large levels of neofunctionalization, we found that most of the recent pairs of gene duplicates evolved asymmetrically, showing divergent patterns of positive selection and gene expression. Genes under selection involved a plethora of biological functions, suggesting that neofunctionalization and tissue specificity, among other evolutionary mechanisms, have orchestrated the evolution of recent paralogs in the pea aphid and may have facilitated host-symbiont cooperation. Our comprehensive phylogenomics analysis allowed us to tackle the history of duplicated genes to pave the road toward understanding the role of gene duplication in ecological adaptation.

The Aphid X Chromosome Is a Dangerous Place for Functionally Important Genes: Diverse Evolution of Hemipteran Genomes Based on Chromosome-Level Assemblies

Different evolutionary forces shape gene content and sequence evolution on autosomes versus sex chromosomes. Location on a sex chromosome can favor male-beneficial or female-beneficial mutations depending on the sex determination system and selective pressure on different sexual morphs. An X0 sex determination can lead to autosomal enrichment of male-biased genes, as observed in some hemipteran insect species. Aphids share X0 sex determination; however, models predict the opposite pattern, due to their unusual life cycles, which alternate between all-female asexual generations and a single sexual generation. Predictions include enrichment of female-biased genes on autosomes and of male-biased genes on the X, in contrast to expectations for obligately sexual species. Robust tests of these models require chromosome-level genome assemblies for aphids and related hemipterans with X0 sex determination and obligate sexual reproduction. In this study, we built the first chromosome-level assembly of a psyllid, an aphid relative with X0 sex determination and obligate sexuality, and compared it with recently resolved chromosome-level assemblies of aphid genomes. Aphid and psyllid X chromosomes differ strikingly. In aphids, female-biased genes are strongly enriched on autosomes and male-biased genes are enriched on the X. In psyllids, male-biased genes are enriched on autosomes. Furthermore, functionally important gene categories of aphids are enriched on autosomes. Aphid X-linked genes and male-biased genes are under relaxed purifying selection, but gene content and order on the X is highly conserved, possibly reflecting constraints imposed by unique chromosomal mechanisms associated with the unusual aphid life cycle.

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.

Chromosome-level genome assembly of the aphid parasitoid Aphidius gifuensis using Oxford Nanopore sequencing and Hi-C technology

Aphidius gifuensis is a parasitoid wasp that has been commercially bred and released in large scale as a biocontrol agent for the management of aphid pests. As a highly efficient endoparasitoid, it is also an important model for exploring mechanisms of parasitism. Currently, artificially bred populations of this wasp are facing rapid decline with undetermined cause, and mechanisms underlying its parasitoid strategy remain poorly understood. Exploring the mechanism behind its population decline and the host–parasitoid relationship is impeded partly due to the lack of a comprehensive genome data for this species. In this study, we constructed a high‐quality reference genome of A. gifuensis using Oxford Nanopore sequencing and Hi‐C (proximity ligation chromatin conformation capture) technology. The final genomic assembly was 156.9 Mb, with a contig N50 length of 3.93 Mb, the longest contig length of 10.4 Mb and 28.89% repetitive sequences. 99.8% of genome sequences were anchored onto six linkage groups. A total of 11,535 genes were predicted, of which 90.53% were functionally annotated. Benchmarking Universal Single‐Copy Orthologs (BUSCO) analysis showed the completeness of assembled genome is 98.3%. We found significantly expanded gene families involved in metabolic processes, transmembrane transport, cell signal communication and oxidoreductase activity, in particular ATP‐binding cassette (ABC) transporter, Cytochrome P450 and venom proteins. The olfactory receptors (ORs) showed significant contraction, which may be associated with the decrease in host recognition. Our study provides a solid foundation for future studies on the molecular mechanisms of population decline as well as host–parasitoid relationship for parasitoid wasps.

Characterization of Rhopalosiphum padi takeout-like genes and their role in insecticide susceptibility

Due to the extensive use of chemical insecticides, the field populations of Rhopalosiphum padi, a serious wheat pest worldwide, have developed resistance to insecticides. Therefore, deep understanding of the mechanisms of the aphid's physiological response to insecticides would be of importance for the management of insecticide resistance in pests. Takeout belongs to a protein superfamily found exclusively in insects. Previous research showed that the takeout gene had various functions in insect physiology and behavior. However, few studies have explored the functions of takeout in insect insecticide susceptibility. The susceptibility of R. padi to imidacloprid and beta-cypermethrin was tested. Thirteen takeout-like genes were identified based on the genome database of R. padi. The number of exons was variable in these takeout-like genes, and nine highly conserved amino acids (two Cysteine, two Proline, four Glycine and one Aspartic acid) were identified. Expression levels of takeout-like-2, takeout-like-3, takeout-like-5, takeout-like-8, takeout-like-10 and takeout-like-11 were significantly increased after imidacloprid treatment; seven genes (takeout-like-1, takeout-like-2, takeout-like-5, takeout-like-6, takeout-like-7, takeout-like-8 and takeout-like-11) tended to be upregulated after beta-cypermethrin treatment. RNA interference results showed that the mortalities of R. padi injected with dsTOL-2, dsTOL-5, dsTOL-8, dsTOL-10 and dsTOL-11 were significantly increased after exposure to imidacloprid in comparison with control (injection of dsGFP). Under two sublethal concentrations of beta-cypermethrin, the silencing of takeout-like-2, takeout-like-5 and takeout-like-11 significantly increased the mortalities of R. padi. These results provide evidence for the involvement of takeout-like genes in insecticide susceptibility of R. padi, which improves our understanding the determinant of insecticide susceptibility.

Evolutionary novelty in the apoptotic pathway of aphids

Apoptotic processes play an important role in the development and physiology of almost all metazoan clades. In the highly diverse group of insects, apoptotic pathways have been characterized in only a few dipteran and lepidopteran species, which may not be representative of all insect species. Here, we report the first complete annotation of the apoptotic pathway in a hemipteran insect, the pea aphid Acyrthosiphon pisum. We showed that its apoptotic pathway is rewired compared to other insects, with a significant increase in the number of inhibitors of apoptosis (IAPs) and evidence for functional diversification and structural modularity of this protein family. These novelties are widespread in the aphid lineage, suggesting a yet not understood novel aphid-specific function of IAPs.

Apoptosis, a conserved form of programmed cell death, shows interspecies differences that may reflect evolutionary diversification and adaptation, a notion that remains largely untested. Among insects, the most speciose animal group, the apoptotic pathway has only been fully characterized in Drosophila melanogaster, and apoptosis-related proteins have been studied in a few other dipteran and lepidopteran species. Here, we studied the apoptotic pathway in the aphid Acyrthosiphon pisum, an insect pest belonging to the Hemiptera, an earlier-diverging and distantly related order. We combined phylogenetic analyses and conserved domain identification to annotate the apoptotic pathway in A. pisum and found low caspase diversity and a large expansion of its inhibitory part, with 28 inhibitors of apoptosis (IAPs). We analyzed the spatiotemporal expression of a selected set of pea aphid IAPs and showed that they are differentially expressed in different life stages and tissues, suggesting functional diversification. Five IAPs are specifically induced in bacteriocytes, the specialized cells housing symbiotic bacteria, during their cell death. We demonstrated the antiapoptotic role of these five IAPs using heterologous expression in a tractable in vivo model, the Drosophila melanogaster developing eye. Interestingly, IAPs with the strongest antiapoptotic potential contain two BIR and two RING domains, a domain association that has not been observed in any other species. We finally analyzed all available aphid genomes and found that they all show large IAP expansion, with new combinations of protein domains, suggestive of evolutionarily novel aphid-specific functions.

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.