Targeted re-sequencing confirms the importance of chemosensory genes in aphid host race differentiation

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.

Insights into chemosensory genes of Pagiophloeus tsushimanus adults using transcriptome and qRT-PCR analysis

Pagiophloeus tsushimanus is a new, destructive, and monophagous weevil pest that thrives on Cinnamomum camphora, found in Shanghai. The functions of chemosensory genes involved in the host location and intraspecific communication of P. tsushimanus remain unknown. The male-female transcriptomes of P. tsushimanus adults were assembled using Illumina sequencing, and we focused on all chemosensory genes in transcriptomes. In general, 58,088 unigenes with a mean length of 1018.19 bp were obtained. In total, 39 odorant binding proteins (OBPs), 10 chemosensory proteins (CSPs), 22 olfactory receptors (ORs), 16 gustatory receptors (GRs), eight ionotropic receptors (IRs), and five sensory neuron membrane proteins (SNMPs) were identified. PtsuOBPs comprised four subfamilies (20 Minus-C, one Plus-C, two Dimer, and 15 Classic). Both PtsuOBPs and PtsuCSPs contained a highly conserved sequence motif of cysteine residues. PtsuORs including one olfactory receptor co-receptors (Ptsu/Orco) comprised seven predicted transmembrane domains. Phylogenetic analysis revealed that PtsuOBPs, PtsuCSPs, and PtsuORs in P. tsushimanus exhibited low homology compared to other insect species. The results of tissue- and sex-specific expression patterns indicated that PtsuOBPs and PtsuORs were highly abundant in the antennae; whereas, PtsuCSPs were not only highly abundant in antennae, but also abdominal apexes, wings, and legs. In conclusion, these results enrich the gene database of P. tsushimanus, which may serve as a basis for identifying novel targets to disrupt olfactory key genes and may provide a reverse validation method to identify attractants for formulating potential eco-friendly control strategies for this pest.

Alarm Pheromone Responses Depend on Genotype, but Not on the Presence of Facultative Endosymbionts in the Pea Aphid Acyrthosiphon pisum

Aphids use an alarm pheromone, E-β farnesene (EBF), to warn conspecifics of potential danger. The antennal sensitivity and behavioural escape responses to EBF can be influenced by different factors. In the pea aphid, Acyrthosiphon pisum, different biotypes are adapted to different legume species, and within each biotype, different genotypes exist, which can carry or not Hamiltonella defensa, a bacterial symbiont that can confer protection against natural enemies. We investigate here the influence of the aphid genotype and symbiotic status on the escape behaviour using a four-way olfactometer and antennal sensitivity for EBF using electroantennograms (EAGs). Whereas the investigated three genotypes from two biotypes showed significantly different escape and locomotor behaviours in the presence of certain EBF doses, the infection with H. defensa did not significantly modify the escape behaviour and only marginally influenced the locomotor behaviour at high doses of EBF. Dose-response curves of EAG amplitudes after stimulation with EBF differed significantly between aphid genotypes in correlation with behavioural differences, whereas antennal sensitivity to EBF did not change significantly as a function of the symbiotic status. The protective symbiont H. defensa does thus not modify the olfactory sensitivity to the alarm pheromone. How EBF sensitivity is modified between genotypes or biotypes remains to be investigated.

Catching the complexity of salmon-louse interactions

The study of host-parasite relationships is an integral part of the immunology of aquatic species, where the complexity of both organisms has to be overlayed with the lifecycle stages of the parasite and immunological status of the host. A deep understanding of how the parasite survives in its host and how they display molecular mechanisms to face the immune system can be applied for novel parasite control strategies. This review highlights current knowledge about salmon and sea louse, two key aquatic animals for aquaculture research worldwide. With the aim to catch the complexity of the salmon-louse interactions, molecular information gleaned through genomic studies are presented. The host recognition system and the chemosensory receptors found in sea lice reveal complex molecular components, that in turn, can be disrupted through specific molecules such as non-coding RNAs.

The chemoreceptors and odorant binding proteins of the soybean and pea aphids

We examined the genome of the soybean aphid, Aphis glycines, and an updated genome assembly of the pea aphid, Acyrthosiphon pisum, for members of the three major families of chemoreceptors, the Odorant Receptors (ORs), Gustatory Receptors (GRs) and Ionotropic Receptors (IRs), as well as the Odorant Binding Proteins (OBPs). The soybean aphid has 47 ORs, 61 GRs, 19 IRs, and 10 OBPs, compared with 87 ORs, 78 Grs, 19 IRs, and 18 OBPs in the pea aphid, with variable numbers of pseudogenes in the OR and GR families. Phylogenetic analysis reveals that while all of the IRs are simple orthologs between these two species, the OR, GR, and OBP families in the pea aphid have experienced major expansions of particular gene lineages and fewer losses of gene lineages. This imbalance in birth-and-death of chemosensory genes has led to the larger pea aphid gene repertoire, which might be related to the broader host range of pea aphids versus the specialization of soybean aphids on a single summer host plant. Examination of the expression levels of these chemosensory genes in parthenogenetic and sexual females and males of pea aphids revealed multiple genes that are differentially expressed in sexual females or males and might be involved in reproductive biology. Examination of the soybean aphid genes in parthenogenetic females under multiple stressors revealed multiple genes whose expression levels changed with heat or starvation stress, the latter potentially important in finding new food sources.

Genomic transitions during host race and species formation

Darwin recognized species as discontinuous, yet considered them to be formed by an incremental process of natural selection. Recent theoretical work on 'genome-wide congealing' is bridging this gap between the gradualism of divergent selection and rapid genome-wide divergence, particularly during ecological speciation-with-gene-flow. Host races and species of phytophagous insects, displaying a spectrum of divergence and gene flow among member taxa, provide model systems for testing predicted non-linear transitions from 'genic' divergence at a few uncoupled loci to 'genomic' divergence with genome-wide coupling of selected loci and strong reproductive isolation. Integrating across natural history, genomics, and evolutionary theory, emerging research suggests a tipping point from 'genic' to 'genomic' divergence between host races and species, during both sympatric speciation and secondary contact.

Intrinsic pre-zygotic reproductive isolation of distantly related pea aphid host races

Human activities may weaken or destroy reproductive isolation between young taxa, leading to their fusion with consequences for population and community ecology. Pea aphid host races are adapted to different legume taxa, providing a degree of pre-mating isolation mediated by habitat choice. Yet, all races can feed and reproduce on the broad bean (Vicia faba), a major crop which represents a 'universal host plant', which can promote hybridization between races. Here, we ask if pea aphid host races have reproductive barriers which prevent or reduce gene flow when they co-occur on the universal host plant. We observed mating behaviour, female survival, number of eggs and egg fertilization rates for three types of crosses: among individuals of the same host race, between closely related host races and between distantly related host races. We did not find significant differences in mating behaviour and female survival among the three types of crosses. However, we observed a drastic reduction in the number of eggs laid, and in the number of fertilized eggs, in distant crosses. We conclude that widespread broad bean cultivation in agriculture may predispose closely related-but not distantly related-host races to hybridize, disrupting reproductive isolation between incipient species.

Identification and Comparison of Chemosensory Genes in the Antennal Transcriptomes of Eucryptorrhynchus scrobiculatus and E. brandti Fed on Ailanthus altissima

The key to the coexistence of two or more species on the same host is ecological niche separation. Adult Eucryptorrhynchus scrobiculatus and E. brandti both feed on the tree of heaven, Ailanthus altissima, but on different sections of the plant. Olfaction plays a vital role in foraging for food resources. Chemosensory genes on the antennae, the main organ for insect olfaction, might explain their feeding differentiation. In the present study, we identified 130 and 129 putative chemosensory genes in E. scrobiculatus and E. brandti, respectively, by antennal transcriptome sequencing, including 31 odorant-binding proteins (OBPs), 11 chemosensory proteins (CSPs), 49 odorant receptors (ORs), 17 ionotropic receptors (IRs), 19 gustatory receptors (GRs), and three sensory neuron membrane proteins (SNMPs) in E. scrobiculatus and 28 OBPs, 11 CSPs, 45 ORs, 25 IRs, 17 GRs, and three SNMPs in E. brandti. We inferred that EscrOBP8 (EscrPBP1), EscrOBP24 (EscrPBP2) and EbraOBP8 (EbraPBP1), EbraOBP24 (EbraPBP2) were putative PBPs by the phylogenetic analysis. We identified species-specific OR transcripts (10 EscrORs and 8 EbraORs) with potential roles in the recognition of specific volatiles of A. altissima. In addition to conserved “antennal IRs,” we also found several “divergent IRs” orthologues in E. scrobiculatus and E. brandti, such as EscrIR16, EbraIR19, and EbraIR20. Compared with other chemosensory genes, GRs between E. scrobiculatus and E. brandti shared lower amino acid identities, which could explain the different feeding habits of the species. We examined OBP expression patterns in various tissues and sexes. Although amino acid sequence similarities were high between EscrOBPs and EbraOBPs, the homologous OBPs showed different tissue expression pattern between two weevils. Our systematic comparison of chemosensory genes in E. scrobiculatus and E. brandti provides a foundation for studies of olfaction and olfactory differentiation in the two weevils as well as a theoretical basis for studying species differentiation.

Fast Evolution and Lineage-Specific Gene Family Expansions of Aphid Salivary Effectors Driven by Interactions with Host-Plants

Effector proteins play crucial roles in plant-parasite interactions by suppressing plant defenses and hijacking plant physiological responses to facilitate parasite invasion and propagation. Although effector proteins have been characterized in many microbial plant pathogens, their nature and role in adaptation to host plants are largely unknown in insect herbivores. Aphids rely on salivary effector proteins injected into the host plants to promote phloem sap uptake. Therefore, gaining insight into the repertoire and evolution of aphid effectors is key to unveiling the mechanisms responsible for aphid virulence and host plant specialization. With this aim in mind, we assembled catalogues of putative effectors in the legume specialist aphid, Acyrthosiphon pisum, using transcriptomics and proteomics approaches. We identified 3,603 candidate effector genes predicted to be expressed in A. pisum salivary glands (SGs), and 740 of which displayed up-regulated expression in SGs in comparison to the alimentary tract. A search for orthologs in 17 arthropod genomes revealed that SG-up-regulated effector candidates of A. pisum are enriched in aphid-specific genes and tend to evolve faster compared with the whole gene set. We also found that a large fraction of proteins detected in the A. pisum saliva belonged to three gene families, of which certain members show evidence consistent with positive selection. Overall, this comprehensive analysis suggests that the large repertoire of effector candidates in A. pisum constitutes a source of novelties promoting plant adaptation to legumes.

The Atlantic salmon (Salmo salar) antimicrobial peptide cathelicidin-2 is a molecular host-associated cue for the salmon louse (Lepeophtheirus salmonis)

Chemical signals are a key element of host-parasite interactions. In marine ecosystems, obligate ectoparasites, such as sea lice, use chemical cues and other sensory signals to increase the probability of encountering a host and to identify appropriate hosts on which they depend to complete their life cycle. The chemical compounds that underlie host identification by the sea lice are not fully described or characterized. Here, we report a novel compound - the Atlantic salmon (Salmo salar) antimicrobial peptide cathelicidin-2 (Cath-2) – that acts as an activation cue for the marine parasitic copepod Lepeophtheirus salmonis. L. salmonis were exposed to 0, 7, 70 and 700 ppb of Cath-2 and neural activity, swimming behaviour and gene expression profiles of animals in response to the peptide were evaluated. The neurophysiological, behavioural and transcriptomic results were consistent: L. salmonis detects Cath-2 as a water-soluble peptide released from the skin of salmon, triggering chemosensory neural activity associated with altered swimming behaviour of copepodids exposed to the peptide, and chemosensory-related genes were up-regulated in copepodids exposed to the peptide. L. salmonis are activated by Cath-2, indicating a tight link between this peptide and the salmon louse chemosensory system.

Why Do Herbivorous Mites Suppress Plant Defenses?

Plants have evolved numerous defensive traits that enable them to resist herbivores. In turn, this resistance has selected for herbivores that can cope with defenses by either avoiding, resisting or suppressing them. Several species of herbivorous mites, such as the spider mites Tetranychus urticae and Tetranychus evansi, were found to maximize their performance by suppressing inducible plant defenses. At first glimpse it seems obvious why such a trait will be favored by natural selection. However, defense suppression appeared to readily backfire since mites that do so also make their host plant more suitable for competitors and their offspring more attractive for natural enemies. This, together with the fact that spider mites are infamous for their ability to resist (plant) toxins directly, justifies the question as to why traits that allow mites to suppress defenses nonetheless seem to be relatively common? We argue that this trait may facilitate generalist herbivores, like T. urticae, to colonize new host species. While specific detoxification mechanisms may, on average, be suitable only on a narrow range of similar hosts, defense suppression may be more broadly effective, provided it operates by targeting conserved plant signaling components. If so, resistance and suppression may be under frequency-dependent selection and be maintained as a polymorphism in generalist mite populations. In that case, the defense suppression trait may be under rapid positive selection in subpopulations that have recently colonized a new host but may erode in relatively isolated populations in which host-specific detoxification mechanisms emerge. Although there is empirical evidence to support these scenarios, it contradicts the observation that several of the mite species found to suppress plant defenses actually are relatively specialized. We argue that in these cases buffering traits may enable such mites to mitigate the negative side effects of suppression in natural communities and thus shield this trait from natural selection.

Identifying genomic hotspots of differentiation and candidate genes involved in the adaptive divergence of pea aphid host races

Identifying the genomic bases of adaptation to novel environments is a long-term objective in evolutionary biology. Because genetic differentiation is expected to increase between locally adapted populations at the genes targeted by selection, scanning the genome for elevated levels of differentiation is a first step towards deciphering the genomic architecture underlying adaptive divergence. The pea aphid Acyrthosiphon pisum is a model of choice to address this question, as it forms a large complex of plant-specialized races and cryptic species, resulting from recent adaptive radiation. Here, we characterized genomewide polymorphisms in three pea aphid races specialized on alfalfa, clover and pea crops, respectively, which we sequenced in pools (poolseq). Using a model-based approach that explicitly accounts for selection, we identified 392 genomic hotspots of differentiation spanning 47.3 Mb and 2,484 genes (respectively, 9.12% of the genome size and 8.10% of its genes). Most of these highly differentiated regions were located on the autosomes, and overall differentiation was weaker on the X chromosome. Within these hotspots, high levels of absolute divergence between races suggest that these regions experienced less gene flow than the rest of the genome, most likely by contributing to reproductive isolation. Moreover, population-specific analyses showed evidence of selection in every host race, depending on the hotspot considered. These hotspots were significantly enriched for candidate gene categories that control host-plant selection and use. These genes encode 48 salivary proteins, 14 gustatory receptors, 10 odorant receptors, five P450 cytochromes and one chemosensory protein, which represent promising candidates for the genetic basis of host-plant specialization and ecological isolation in the pea aphid complex. Altogether, our findings open new research directions towards functional studies, for validating the role of these genes on adaptive phenotypes.

The Evolutionary Dynamics of the Odorant Receptor Gene Family in Corbiculate Bees

Insects rely on chemical information to locate food, choose mates, and detect potential predators. It has been hypothesized that adaptive changes in the olfactory system facilitated the diversification of numerous insect lineages. For instance, evolutionary changes of Odorant Receptor (OR) genes often occur in parallel with modifications in life history strategies. Corbiculate bees display a diverse array of behaviors that are controlled through olfaction, including varying degrees of social organization, and manifold associations with floral resources. Here we investigated the molecular mechanisms driving the evolution of the OR gene family in corbiculate bees in comparison to other chemosensory gene families. Our results indicate that the genomic organization of the OR gene family has remained highly conserved for ∼80 Myr, despite exhibiting major changes in repertoire size among bee lineages. Moreover, the evolution of OR genes appears to be driven mostly by lineage-specific gene duplications in few genomic regions that harbor large numbers of OR genes. A selection analysis revealed that OR genes evolve under positive selection, with the strongest signals detected in recently duplicated copies. Our results indicate that chromosomal translocations had a minimal impact on OR evolution, and instead local molecular mechanisms appear to be main drivers of OR repertoire size. Our results provide empirical support to the longstanding hypothesis that positive selection shaped the diversification of the OR gene family. Together, our results shed new light on the molecular mechanisms underlying the evolution of olfaction in insects.

Rapid evolution of aphid pests in agricultural environments

Aphids constitute a major group of crop pests that inflict serious damages to plants, both directly by ingesting phloem and indirectly as vectors of numerous diseases. In response to intense and repeated human-induced pressures, such as insecticide treatments, the use of resistant plants and biological agents, aphids have developed a series of evolutionary responses relying on adaptation and phenotypic plasticity. In this review, we highlight some remarkable evolutionary responses to anthropogenic pressures in agroecosystems and discuss the mechanisms underlying the ecological and evolutionary success of aphids. We outline the peculiar mode of reproduction, the polyphenism for biologically important traits and the diverse and flexible associations with microbial symbionts as key determinants of adaptive potential and pest status of aphids.

Intrinsic incompatibilities evolving as a by-product of divergent ecological selection: Considering them in empirical studies on divergence with gene flow

The possibility of intrinsic barriers to gene flow is often neglected in empirical research on local adaptation and speciation with gene flow, for example when interpreting patterns observed in genome scans. However, we draw attention to the fact that, even with gene flow, divergent ecological selection may generate intrinsic barriers involving both ecologically selected and other interacting loci. Mechanistically, the link between the two types of barriers may be generated by genes that have multiple functions (i.e., pleiotropy), and/or by gene interaction networks. Because most genes function in complex networks, and their evolution is not independent of other genes, changes evolving in response to ecological selection can generate intrinsic barriers as a by-product. A crucial question is to what extent such by-product barriers contribute to divergence and speciation-that is whether they stably reduce gene flow. We discuss under which conditions by-product barriers may increase isolation. However, we also highlight that, depending on the conditions (e.g., the amount of gene flow and the strength of selection acting on the intrinsic vs. the ecological barrier component), the intrinsic incompatibility may actually destabilize barriers to gene flow. In practice, intrinsic barriers generated as a by-product of divergent ecological selection may generate peaks in genome scans that cannot easily be interpreted. We argue that empirical studies on divergence with gene flow should consider the possibility of both ecological and intrinsic barriers. Future progress will likely come from work combining population genomic studies, experiments quantifying fitness and molecular studies on protein function and interactions.

A nonrandom subset of olfactory genes is associated with host preference in the fruit fly Drosophila orena

Specialization onto different host plants has been hypothesized to be a major driver of diversification in insects, and traits controlling olfaction have been shown to play a fundamental role in host preferences. A diverse set of olfactory genes control olfactory traits in insects, and it remains unclear whether specialization onto different hosts is likely to involve a nonrandom subset of these genes. Here, we test the role of olfactory genes in a novel case of specialization in Drosophila orena. We report the first population‐level sample of D. orena on the West African island of Bioko, since its initial collection in Cameroon in 1975, and use field experiments and behavioral assays to show that D. orena has evolved a strong preference for waterberry (Syzygium staudtii). We then show that a nonrandom subset of genes controlling olfaction‐–those controlling odorant‐binding and chemosensory proteins–‐have an enriched signature of positive selection relative to the rest of the D. orena genome. By comparing signatures of positive selection on olfactory genes between D. orena and its sister species, D. erecta we show that odorant‐binding and chemosensory have evidence of positive selection in both species; however, overlap in the specific genes with evidence of selection in these two classes is not greater than expected by chance. Finally, we use quantitative complementation tests to confirm a role for seven olfactory loci in D. orena’s preference for waterberry fruit. Together, our results suggest that D. orena and D. erecta have specialized onto different host plants through convergent evolution at the level of olfactory gene family, but not at specific olfactory genes.

Molecular ecology studies of species radiations: current research gaps, opportunities and challenges

Understanding the drivers and limits of species radiations is a crucial goal of evolutionary genetics and molecular ecology, yet research on this topic has been hampered by the notorious difficulty of connecting micro- and macroevolutionary approaches to studying the drivers of diversification. To chart the current research gaps, opportunities and challenges of molecular ecology approaches to studying radiations, we examine the literature in the journal Molecular Ecology and revisit recent high-profile examples of evolutionary genomic research on radiations. We find that available studies of radiations are highly unevenly distributed among taxa, with many ecologically important and species-rich organismal groups remaining severely understudied, including arthropods, plants and fungi. Most studies employed molecular methods suitable over either short or long evolutionary time scales, such as microsatellites or restriction site-associated DNA sequencing (RAD-seq) in the former case and conventional amplicon sequencing of organellar DNA in the latter. The potential of molecular ecology studies to address and resolve patterns and processes around the species level in radiating groups of taxa is currently limited primarily by sample size and a dearth of information on radiating nuclear genomes as opposed to organellar ones. Based on our literature survey and personal experience, we suggest possible ways forward in the coming years. We touch on the potential and current limitations of whole-genome sequencing (WGS) in studies of radiations. We suggest that WGS and targeted ('capture') resequencing emerge as the methods of choice for scaling up the sampling of populations, species and genomes, including currently understudied organismal groups and the genes or regulatory elements expected to matter most to species radiations.