Large-scale candidate gene scan reveals the role of chemoreceptor genes in host plant specialization and speciation in the pea aphid

Partial Correspondence between Host Plant-Related Differentiation and Symbiotic Bacterial Community in a Polyphagous Insect

Host plants play a vital role in insect population differentiation, while symbiotic associations between bacteria and insects are ubiquitous in nature. However, existing studies have given limited attention to the connection between host-related differentiation and symbiotic bacterial communities in phytophagous insects. In this study, we collected 58 samples of Aphis odinae from different host plants in southern China and constructed phylogenetic trees to investigate their differentiation in relation to host plants. We also selected aphid samples from the five most preferred host plants and analyzed their symbiotic bacterial composition using Illumina sequencing of the V3-V4 hypervariable region of the 16S rRNA gene. The phylogeny and symbiotic bacterial community structure of A. odinae populations on different host plants showed that samples from Triadica sebifera (Euphorbiaceae) had a consistent presence of Wolbachia as the predominant secondary symbiont and suggested the possibility of undergoing differentiation. Conversely, although differentiation was observed in samples from Rhus chinensis (Anacardiaceae), no consistent presence of predominant secondary symbionts was found. Additionally, the samples from Heptapleurum heptaphyllum (Araliaceae) consistently carried Serratia, but no host differentiation was evident. In summary, this study reveals a partial correspondence between symbiotic bacterial communities and host-related differentiation in A. odinae. The findings contribute to our understanding of the microevolutionary influencing the macroevolutionary relationships between bacterial symbionts and phytophagous insects. The identification of specific symbionts associated with host-related differentiation provides valuable insights into the intricate dynamics of insect-bacteria interactions.

The genome sequence of a spider mite, Tetranychus truncatus, provides insights into interspecific host range variation and the genetic basis of adaptation to a low-quality host plant

The phytophagous mite Tetranychus truncatus is a serious pest in East Asia but has a relatively narrower host range than the pest mite Tetranychus urticae, which can feed on over 1200 plant species. Here, we generated a high-quality chromosomal level genome of T. truncatus and compared it with that of T. urticae, with an emphasis on the genes related to detoxification and chemoreception, to explore the genomic basis underlying the evolution of host range. We also conducted population genetics analyses (in 86 females from 10 populations) and host transfer experiments (in 4 populations) to investigate transcription changes following transfer to a low-quality host (Solanum melongena, eggplant), and we established possible connections between fitness on eggplant and genes related to detoxification and chemoreception. We found that T. truncatus has fewer genes related to detoxification, transport, and chemoreception than T. urticae, with a particularly strong reduction in gustatory receptor (GR) genes. We also found widespread transcriptional variation among T. truncatus populations, which varied in fitness on eggplant. We characterized selection on detoxification-related genes through ω values and found a negative correlation between expression levels and ω values. Based on the transcription results, as well as the fitness and genetic differences among populations, we identified genes potentially involved in adaptation to eggplant in T. truncatus. Our work provides a genomic resource for this pest mite and new insights into mechanisms underlying the adaptation of herbivorous mites to host plants.

Molecular Mechanisms Underlying Host Plant Specificity in Aphids

Aphids are serious pests of agricultural and ornamental plants and important model systems for hemipteran-plant interactions. The long evolutionary history of aphids with their host plants has resulted in a variety of systems that provide insight into the different adaptation strategies of aphids to plants and vice versa. In the past, various plant-aphid interactions have been documented, but lack of functional tools has limited molecular studies on the mechanisms of plant-aphid interactions. Recent technological advances have begun to reveal plant-aphid interactions at the molecular level and to increase our knowledge of the mechanisms of aphid adaptation or specialization to different host plants. In this article, we compile and analyze available information on plant-aphid interactions, discuss the limitations of current knowledge, and argue for new research directions. We advocate for more work that takes advantage of natural systems and recently established molecular techniques to obtain a comprehensive view of plant-aphid interaction mechanisms.

Genome sequencing of a predominant clonal lineage of the grain aphid Sitobion avenae

The English grain aphid, Sitobion avenae, is a cosmopolitan pest that feeds on cereals, provoking substantial yield losses by injuring plant tissue and by vectoring plant viruses. Here we report a highly complete, de novo draft genome of the grain aphid using long-read sequencing. We generated an assembly of 2740 contigs with a N50 of 450 kb. We compared this draft genome with that of other aphid species, inspecting gene family evolution, genome-wide positive selection, and searched for horizontal gene transfer events. In addition, we described a recent copy number variant expansion of gene families involving aconitase, ABC transporter, and esterase genes that could be associated with resistance to insecticides and plant chemical defenses. This S. avenae genome obtained from a predominant invasive genotype can provide a framework for studying the spatial-temporal success of these clonal lineages in invaded agroecosystems.

Evolution of the insect PPK gene family

Insect pickpocket (PPK) receptors mediate diverse functions, among them the detection of mechano- and chemo-sensory stimuli. Notwithstanding their relevance, studies on their evolution only focused on Drosophila. We have analyzed the genomes of 26 species of 8 orders including holometabolous and hemimetabolous insects (Blattodea, Orthoptera, Hemiptera, Phthiraptera, Hymenoptera, Lepidoptera, Coleoptera, and Diptera), to characterize the evolution of this gene family. PPKs were detected in all genomes analyzed, with 578 genes distributed in 7 subfamilies. According to our phylogeny ppk17 is the most divergent member, composing the new subfamily VII. PPKs evolved under a gene birth-and-death model that generated lineage-specific expansions usually located in clusters, while purifying selection affected several orthogroups. Subfamily V was the largest, including a mosquito-specific expansion that can be considered a new target for pest control. PPKs present a high gene turnover generating considerable variation. On one hand, Musca domestica (59), Aedes albopictus (51), Culex quinquefasciatus (48), and Blattella germanica (41) presented the largest PPK repertoires. On the other hand, Pediculus humanus (only ppk17), bees and ants (6-9) had the smallest PPK sets. A subset of prevalent PPKs was identified, indicating very conserved functions for these receptors. Finally, at least twenty percent of the sequences presented calmodulin-binding motifs, suggesting that these PPKs may amplify sensory responses similarly as proposed for D. melanogaster ppk25. Overall, this work characterized the evolutionary history of these receptors revealing relevant unknown gene sequence features and clade-specific expansions.

The draft genome of the specialist flea beetle Altica viridicyanea (Coleoptera: Chrysomelidae)

BACKGROUND

Altica (Coleoptera: Chrysomelidae) is a highly diverse and taxonomically challenging flea beetle genus that has been used to address questions related to host plant specialization, reproductive isolation, and ecological speciation. To further evolutionary studies in this interesting group, here we present a draft genome of a representative specialist, Altica viridicyanea, the first Alticinae genome reported thus far.

RESULTS

The genome is 864.8 Mb and consists of 4490 scaffolds with a N50 size of 557 kb, which covered 98.6% complete and 0.4% partial insect Benchmarking Universal Single-Copy Orthologs. Repetitive sequences accounted for 62.9% of the assembly, and a total of 17,730 protein-coding gene models and 2462 non-coding RNA models were predicted. To provide insight into host plant specialization of this monophagous species, we examined the key gene families involved in chemosensation, detoxification of plant secondary chemistry, and plant cell wall-degradation.

CONCLUSIONS

The genome assembled in this work provides an important resource for further studies on host plant adaptation and functionally affiliated genes. Moreover, this work also opens the way for comparative genomics studies among closely related Altica species, which may provide insight into the molecular evolutionary processes that occur during ecological speciation.

Speciation across the Tree of Life

Much of what we know about speciation comes from detailed studies of well-known model systems. Although there have been several important syntheses on speciation, few (if any) have explicitly compared speciation among major groups across the Tree of Life. Here, we synthesize and compare what is known about key aspects of speciation across taxa, including bacteria, protists, fungi, plants, and major animal groups. We focus on three main questions. Is allopatric speciation predominant across groups? How common is ecological divergence of sister species (a requirement for ecological speciation), and on what niche axes do species diverge in each group? What are the reproductive isolating barriers in each group? Our review suggests the following patterns. (i) Based on our survey and projected species numbers, the most frequent speciation process across the Tree of Life may be co-speciation between endosymbiotic bacteria and their insect hosts. (ii) Allopatric speciation appears to be present in all major groups, and may be the most common mode in both animals and plants, based on non-overlapping ranges of sister species. (iii) Full sympatry of sister species is also widespread, and may be more common in fungi than allopatry. (iv) Full sympatry of sister species is more common in some marine animals than in terrestrial and freshwater ones. (v) Ecological divergence of sister species is widespread in all groups, including ~70% of surveyed species pairs of plants and insects. (vi) Major axes of ecological divergence involve species interactions (e.g. host-switching) and habitat divergence. (vii) Prezygotic isolation appears to be generally more widespread and important than postzygotic isolation. (viii) Rates of diversification (and presumably speciation) are strikingly different across groups, with the fastest rates in plants, and successively slower rates in animals, fungi, and protists, with the slowest rates in prokaryotes. Overall, our study represents an initial step towards understanding general patterns in speciation across all organisms.

Comparative transcriptome analysis of three invasive leafminer flies provides insights into interspecific competition

Liriomyza spp. (Diptera: Agromyzidae) represent a group of economically-significant highly polyphagous pests of plants grown in field and greenhouse conditions. Liriomyza spp. share similar biological and morphological characteristics, and complex interspecific interactions have been documented among these species in various geographical regions. Where the displacement of one of these species by the other has been studied, no unique mechanisms have been identified as causing it. The impact of competitive factors (such as, insecticide tolerance, thermotolerance, and adaptability to cropping systems) may be unique to specific geographic regions of Liriomyza spp., but more research is needed to confirm these hypotheses. In this study, RNA-seq was used to determine the transcriptomes of three closely-related leafminers, e.g. L. sativae, L. trifolii, and L. huidobrensis. Over 20 Gb of clean reads were generated and assembled into unique transcriptomes, and 38,747 unigenes were annotated in different databases. In pairwise comparisons, L. trifolii and L. sativae had more up-regulated genes than L. huidobrensis. With respect to common differentially-expressed genes (Co-DEGs), the three leafminers exhibited distinct groups of highly-expressed gene clusters. When genes related to competitive factors were compared, expression patterns in L. trifolii and L. sativae were more closely related to each other than to L. huidobrensis. The data suggest that DEGs involved in competitive factors may play a key role in competition and displacement of leafminers. The divergent genes identified in this study will be valuable in revealing possible mechanisms of invasion, displacement and interspecific competition in Liriomyza spp.

Genomic content of chemosensory receptors in two sister blister beetles facilitates characterization of chemosensory evolution

Background

More than 2500 species belong to the Meloidae family (Coleoptera: Tenebrionoidea), members of which produce the potent defensive blistering agent cantharidin and are commonly known as blister beetles or Spanishflies. Cantharidin has recently been used for cancer therapy. Hycleus cichorii and Hycleus phaleratus have been used in traditional Chinese medicine for more than 2000 years due to their ability to biosynthesize cantharidin. To understand the role of the chemosensory system in beetle evolution, we comparatively analysed the chemosensory receptor families of both blister beetle species and compared them with those of other beetles.

Results

We identified 89 odorant receptors (ORs), 86 gustatory receptors (GRs), and 45 ionotropic receptors (IRs) in H. phaleratus and 149 ORs, 102 GRs and 50 IRs in H. cichorii. Nine groups of beetle ORs were recovered, and a similar pattern of ORs in Coleoptera emerged. Two evident expanded clades in Hycleus (Groups 5A and 3) were reconstructed in the phylogenetic tree. Four of eight genes with evidence of positive selection were clustered in the expanded clades of Group 5A. Three, eight and three orthologous pairs of CO₂, sugar and fructose receptors, respectively, were identified in both blister beetles. Two evident expanded clades of putative bitter GRs in Hycleus were also found, and the GR in one clade had notably low divergence. Interestingly, IR41a was specifically expanded in blister beetles compared to other insects identified to date, and IR75 was also clearly expanded in both blister beetles based on our phylogenetic tree analysis. Moreover, evidence of positive selection was detected for eight ORs, three GRs and two IRs, half of which were from five duplicate clades.

Conclusions

We first annotated the chemosensory receptor families in a pair of sister beetle genomes (Meloidae: Hycleus), which facilitated evolutionary analysis of chemosensory receptors between sibling species in the Coleoptera group. Our analysis suggests that changes in chemosensory receptors have a possible role in chemical-based species evolution in blister beetles. Future studies should include more species to verify this correlation, which will help us understand the evolution of blister beetles.

Genomic signatures of host-associated divergence and adaptation in a coral-eating snail, Coralliophila violacea (Kiener, 1836)

The fluid nature of the ocean, combined with planktonic dispersal of marine larvae, lowers physical barriers to gene flow. However, divergence can still occur despite gene flow if strong selection acts on populations occupying different ecological niches. Here, we examined the population genomics of an ectoparasitic snail, Coralliophila violacea (Kiener 1836), that specializes on Porites corals in the Indo-Pacific. Previous genetic analyses revealed two sympatric lineages associated with different coral hosts. In this study, we examined the mechanisms promoting and maintaining the snails' adaptation to their coral hosts. Genome-wide single nucleotide polymorphism (SNP) data from type II restriction site-associated DNA (2b-RAD) sequencing revealed two differentiated clusters of C. violacea that were largely concordant with coral host, consistent with previous genetic results. However, the presence of some admixed genotypes indicates gene flow from one lineage to the other. Combined, these results suggest that differentiation between host-associated lineages of C. violacea is occurring in the face of ongoing gene flow, requiring strong selection. Indeed, 2.7% of all SNP loci were outlier loci (73/2,718), indicative of divergence with gene flow, driven by adaptation of each C. violacea lineage to their specific coral hosts.

Identification and comparative expression profiles of chemosensory genes in major chemoreception organs of a notorious pests, Laodelphax striatellus

The small brown planthopper, Laodelphax striatellus (Stål) (SBPH), is a notorious rice pest in East Asia and damages the host by feeding on the phloem and transmitting virus particles. Although SBPH relies on chemosensory perception for seeking the host, courtship, selecting oviposition sites and spreading virus particles, a systematic study of chemosensory genes in SBPH is lacking. In this study, we identified multi-gene chemosensory families from the transcriptome of SBPH olfactory organs and analyzed their expression patterns in male and female tissues. Among the chemosensory genes, 14 odorant-binding proteins (OBPs), 12 chemosensory proteins (CSPs), 7 sensory neuron membrane proteins (SNMPs) and 95 odorant receptors (ORs) were identified and annotated in SBPH olfactory organs. Based on expression profile and phylogenetic analysis, LstrOBP1, 2, 5, 6, 7, 10, LstrSNMP1, and most LstrORs showed an antennae-enriched expression pattern, which suggests an olfactory role for these genes. Relative expression of LstrOBPs was validated by quantitative real-time PCR. Our findings provide the genetic information for disrupting the feeding behavior of SBPH, which is essential for developing eco-friendly pest management technologies.

Genome-wide identification of chemosensory receptor genes in the small brown planthopper, Laodelphax striatellus

The small brown planthopper (SBPH), Laodelphax striatellus is one of the major insect pests of rice, but little is known about the molecular-level means by which it locates its hosts. SBPH host-seeking behavior heavily relies on chemosensory receptors (CRs). In this study, we utilized genome analysis of the SBPH to identify 169 CRs, including: 133 odorant receptors (ORs), 13 gustatory receptors (GRs) and 23 ionotropic receptors (IRs). The phylogenetic relationships of OR genes from three rice planthoppers and other insect species revealed that the odorant co-receptor (Orco) clade is the most conserved group. Among the candidate GRs, two sugar receptors and five fructose receptors have been identified but no carbon dioxide receptors investigated. Furthermore, we identified homologs of the three highly conserved IR co-receptors. The obtained results will provide us with precious information needed to better understand the interaction between insect pests and crop plants required for effective crop protection.

Comparing Pool-seq, Rapture, and GBS genotyping for inferring weak population structure: The American lobster (Homarus americanus) as a case study

Unraveling genetic population structure is challenging in species potentially characterized by large population size and high dispersal rates, often resulting in weak genetic differentiation. Genotyping a large number of samples can improve the detection of subtle genetic structure, but this may substantially increase sequencing cost and downstream bioinformatics computational time. To overcome this challenge, alternative, cost-effective sequencing approaches, namely Pool-seq and Rapture, have been developed. We empirically measured the power of resolution and congruence of these two methods in documenting weak population structure in nonmodel species with high gene flow comparatively to a conventional genotyping-by-sequencing (GBS) approach. For this, we used the American lobster (Homarus americanus) as a case study. First, we found that GBS, Rapture, and Pool-seq approaches gave similar allele frequency estimates (i.e., correlation coefficient over 0.90) and all three revealed the same weak pattern of population structure. Yet, Pool-seq data showed F ST estimates three to five times higher than GBS and Rapture, while the latter two methods returned similar F ST estimates, indicating that individual-based approaches provided more congruent results than Pool-seq. We conclude that despite higher costs, GBS and Rapture are more convenient approaches to use in the case of species exhibiting very weak differentiation. While both GBS and Rapture approaches provided similar results with regard to estimates of population genetic parameters, GBS remains more cost-effective in project involving a relatively small numbers of genotyped individuals (e.g., <1,000). Overall, this study illustrates the complexity of estimating genetic differentiation and other summary statistics in complex biological systems characterized by large population size and migration rates.

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.

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.

Positive selection at sites of chemosensory genes is associated with the recent divergence and local ecological adaptation in cactophilic Drosophila

Background

Adaptation to new hosts in phytophagous insects often involves mechanisms of host recognition by genes of sensory pathways. Most often the molecular evolution of sensory genes has been explained in the context of the birth-and-death model. The role of positive selection is less understood, especially associated with host adaptation and specialization. Here we aim to contribute evidence for this latter hypothesis by considering the case of Drosophila mojavensis, a species with an evolutionary history shaped by multiple host shifts in a relatively short time scale, and its generalist sister species, D. arizonae.

Results

We used a phylogenetic and population genetic analysis framework to test for positive selection in a subset of four chemoreceptor genes, one gustatory receptor (Gr) and three odorant receptors (Or), for which their expression has been previously associated with host shifts. We found strong evidence of positive selection at several amino acid sites in all genes investigated, most of which exhibited changes predicted to cause functional effects in these transmembrane proteins. A significant portion of the sites identified as evolving positively were largely found in the cytoplasmic region, although a few were also present in the extracellular domains.

Conclusions

The pattern of substitution observed suggests that some of these changes likely had an effect on signal transduction as well as odorant recognition and protein-protein interactions. These findings support the role of positive selection in shaping the pattern of variation at chemosensory receptors, both during the specialization onto one or a few related hosts, but as well as during the evolution and adaptation of generalist species into utilizing several hosts.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1250-x) contains supplementary material, which is available to authorized users.

Measuring Genetic Differentiation from Pool-seq Data

The advent of high throughput sequencing and genotyping technologies enables the comparison of patterns of polymorphisms at a very large number of markers. While the characterization of genetic structure from individual sequencing data remains expensive for many nonmodel species, it has been shown that sequencing pools of individual DNAs (Pool-seq) represents an attractive and cost-effective alternative. However, analyzing sequence read counts from a DNA pool instead of individual genotypes raises statistical challenges in deriving correct estimates of genetic differentiation. In this article, we provide a method-of-moments estimator of [Formula: see text] for Pool-seq data, based on an analysis-of-variance framework. We show, by means of simulations, that this new estimator is unbiased and outperforms previously proposed estimators. We evaluate the robustness of our estimator to model misspecification, such as sequencing errors and uneven contributions of individual DNAs to the pools. Finally, by reanalyzing published Pool-seq data of different ecotypes of the prickly sculpin Cottus asper, we show how the use of an unbiased [Formula: see text] estimator may question the interpretation of population structure inferred from previous analyses.

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.

The genomic landscape at a late stage of stickleback speciation: High genomic divergence interspersed by small localized regions of introgression

Speciation is a continuous process and analysis of species pairs at different stages of divergence provides insight into how it unfolds. Previous genomic studies on young species pairs have revealed peaks of divergence and heterogeneous genomic differentiation. Yet less known is how localised peaks of differentiation progress to genome-wide divergence during the later stages of speciation in the presence of persistent gene flow. Spanning the speciation continuum, stickleback species pairs are ideal for investigating how genomic divergence builds up during speciation. However, attention has largely focused on young postglacial species pairs, with little knowledge of the genomic signatures of divergence and introgression in older stickleback systems. The Japanese stickleback species pair, composed of the Pacific Ocean three-spined stickleback (Gasterosteus aculeatus) and the Japan Sea stickleback (G. nipponicus), which co-occur in the Japanese islands, is at a late stage of speciation. Divergence likely started well before the end of the last glacial period and crosses between Japan Sea females and Pacific Ocean males result in hybrid male sterility. Here we use coalescent analyses and Approximate Bayesian Computation to show that the two species split approximately 0.68-1 million years ago but that they have continued to exchange genes at a low rate throughout divergence. Population genomic data revealed that, despite gene flow, a high level of genomic differentiation is maintained across the majority of the genome. However, we identified multiple, small regions of introgression, occurring mainly in areas of low recombination rate. Our results demonstrate that a high level of genome-wide divergence can establish in the face of persistent introgression and that gene flow can be localized to small genomic regions at the later stages of speciation with gene flow.

Molecular characterization and evolution of a chemosensory receptor gene family in three notorious rice planthoppers, Nilaparvata lugens, Sogatella furcifera and Laodelphax striatellus, based on genome and transcriptome analyses

BACKGROUND

The white-backed planthopper (WBPH) Sogatella furcifera, the brown planthopper (BPH) Nilaparvata lugens, and the small brown planthopper (SBPH) Laodelphax striatellus (Hemiptera: Delphacidae) are rice pests that damage rice plants by sap-sucking and by transmitting viruses. Host-seeking behavior involves chemosensory receptor genes that include odorant receptors (ORs), ionotropic receptors (IRs) and gustatory receptors (GRs).

RESULTS

We used genome and transcriptome data to identify 141 ORs, 28 GRs and 25 IRs in BPH; 135 ORs, 18 GRs and 16 IRs in WBPH; and 37 ORs, 14 GRs and 6 IRs in SBPH. A phylogenetic analysis identified several specific OR clades of rice planthoppers, the results indicating that these OR members might be used to respond to specific host volatiles. OR co-receptor (Orco) is the most conserved and essential OR gene among these species and RNA interference (RNAi) can decrease their mRNA expression level to <50%. RNAi knockdown rice planthoppers were anosmia and were unable to seek or locate rice plants in behavioral tests.

CONCLUSION

The results demonstrate the importance of the planthopper Orco genes in locating rice plants. This information may aid in the development of RNAi-based transgenic rice and other pest management technologies. © 2018 Society of Chemical Industry.

From population connectivity to the art of striping Russian dolls: the lessons from Pocillopora corals

Here, we examined the genetic variability in the coral genus Pocillopora, in particular within the Primary Species Hypothesis PSH09, identified by Gélin, Postaire, Fauvelot and Magalon (2017) using species delimitation methods [also named Pocillopora eydouxi/meandrina complex sensu, Schmidt-Roach, Miller, Lundgren, & Andreakis (2014)] and which was found to split into three secondary species hypotheses (SSH09a, SSH09b, and SSH09c) according to assignment tests using multi-locus genotypes (13 microsatellites). From a large sampling (2,507 colonies) achieved in three marine provinces [Western Indian Ocean (WIO), Tropical Southwestern Pacific (TSP), and Southeast Polynesia (SEP)], genetic structuring analysis conducted with two clustering analyses (structure and DAPC) using 13 microsatellites revealed that SSH09a was restricted to the WIO while SSH09b and SSH09c were almost exclusively in the TSP and SEP. More surprisingly, each SSH split into two to three genetically differentiated clusters, found in sympatry at the reef scale, leading to a pattern of nested hierarchical levels (PSH > SSH > cluster), each level hiding highly differentiated genetic groups. Thus, rather than structured populations within a single species, these three SSHs, and even the eight clusters, likely represent distinct genetic lineages engaged in a speciation process or real species. The issue is now to understand which hierarchical level (SSH, cluster, or even below) corresponds to the species one. Several hypotheses are discussed on the processes leading to this pattern of mixed clusters in sympatry, evoking formation of reproductive barriers, either by allopatric speciation or habitat selection.

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

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

The chemical signatures underlying host plant discrimination by aphids

The diversity of phytophagous insects is largely attributable to speciation involving shifts between host plants. These shifts are mediated by the close interaction between insects and plant metabolites. However, there has been limited progress in understanding the chemical signatures that underlie host preferences. We use the pea aphid (Acyrthosiphon pisum) to address this problem. Host-associated races of pea aphid discriminate between plant species in race-specific ways. We combined metabolomic profiling of multiple plant species with behavioural tests on two A. pisum races, to identify metabolites that explain variation in either acceptance or discrimination. Candidate compounds were identified using tandem mass spectrometry. Our results reveal a small number of compounds that explain a large proportion of variation in the differential acceptability of plants to A. pisum races. Two of these were identified as L-phenylalanine and L-tyrosine but it may be that metabolically-related compounds directly influence insect behaviour. The compounds implicated in differential acceptability were not related to the set correlated with general acceptability of plants to aphids, regardless of host race. Small changes in response to common metabolites may underlie host shifts. This study opens new opportunities for understanding the mechanistic basis of host discrimination and host shifts in insects.

Genomic and phenotypic analysis of Vavilov's historic landraces reveals the impact of environment and genomic islands of agronomic traits

The Vavilov Institute of Plant Genetic Resources (VIR), in St. Petersburg, Russia, houses a unique genebank, with historical collections of landraces. When they were collected, the geographical distribution and genetic diversity of most crops closely reflected their historical patterns of cultivation established over the preceding millennia. We employed a combination of genomics, computational biology and phenotyping to characterize VIR's 147 chickpea accessions from Turkey and Ethiopia, representing chickpea's center of origin and a major location of secondary diversity. Genotyping by sequencing identified 14,059 segregating polymorphisms and genome-wide association studies revealed 28 GWAS hits in potential candidate genes likely to affect traits of agricultural importance. The proportion of polymorphisms shared among accessions is a strong predictor of phenotypic resemblance, and of environmental similarity between historical sampling sites. We found that 20 out of 28 polymorphisms, associated with multiple traits, including days to maturity, plant phenology, and yield-related traits such as pod number, localized to chromosome 4. We hypothesize that selection and introgression via inadvertent hybridization between more and less advanced morphotypes might have resulted in agricultural improvement genes being aggregated to genomic 'agro islands', and in genotype-to-phenotype relationships resembling widespread pleiotropy.

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.

Genetic divergence with ongoing gene flow is maintained by the use of different hosts in phytophagous ladybird beetles genus Henosepilachna

Adaptation to different environments can promote population divergence via natural selection even in the presence of gene flow - a phenomenon that typically occurs during ecological speciation. To elucidate how natural selection promotes and maintains population divergence during speciation, we investigated the population genetic structure, degree of gene flow and heterogeneous genomic divergence in three closely related Japanese phytophagous ladybird beetles: Henosepilachna pustulosa, H. niponica and H. yasutomii. These species act as a generalist, a wild thistle (Cirsium spp.) specialist and a blue cohosh (Caulophyllum robustum) specialist, respectively, and their ranges differ accordingly. The two specialist species widely co-occur but are reproductively isolated solely due to their high specialization to a particular host plant. Genomewide amplified fragment-length polymorphism (AFLP) markers and mitochondrial cytochrome c oxidase subunit I (COI) gene sequences demonstrated obvious genomewide divergence associated with both geographic distance and ecological divergence. However, a hybridization assessment for both AFLP loci and the mitochondrial sequences revealed a certain degree of unidirectional gene flow between the two sympatric specialist species. Principal coordinates analysis (PCoA) based on all of the variable AFLP loci demonstrated that there are genetic similarities between populations from adjacent localities irrespective of the species (i.e. host range). However, a further comparative genome scan identified a few fractions of loci representing approximately 1% of all loci as different host-associated outliers. These results suggest that these three species had a complex origin, which could be obscured by current gene flow, and that ecological divergence can be maintained with only a small fraction of the genome is related to different host use even when there is a certain degree of gene flow between sympatric species pairs.

Comparative Transcriptional Analysis of the Host-Specialized Aphids Aphis gossypii (Hemiptera: Aphididae)

Host specialization is an ubiquitous character in aphid populations. Many polyphagous aphid populations usually consist of several subpopulations that have strong fidelity to a specific host or a subset of host range. Host specialization is an evolutional result of food habit of insects. However, genetic basis and molecular mechanism of host specialization are still unclear. In this study, we presented a comparative analysis on global gene expression profiles of three lineages of Aphis gossypii Glover: cotton-specialized (CO), cucurbit-specialized (CU), and CU reared on cowpea (CU-cowpea), using RNA-Seq method. More than 157 million clean reads and 38,398 different unigenes were generated from transcriptomes of these three aphid lineages. The 1,106 down- and 2,835 up-regulated genes were found between CO and CU, and 812 down- and 14,492 up-regulated genes between CU-cowpea and CU. Differentially expressed genes between CO and CU were enriched in sugar metabolism, immune system process, pathogen infection or symbiosis, and salivary secretion. Genes associated with cytochrome P450, major facilitator superfamily, and salivary effector were differentially expressed between CO and CU, which might be involved in determining host specialization. UDP-glycosyltransferases genes were sensitive to host shift. Carboxylesterases and digestion-related protease genes were related to both the host specialization and host shift of aphids. Expression levels of 22 out of 24 genes of CO and CU measured by RT-qPCR method were as similar as the results from RNA-seq method. This study provides a road map for future study on molecular mechanism of host specialization in aphids.

Population Structure, Diversity and Reproductive Mode of the Grape Phylloxera (Daktulosphaira vitifoliae) across Its Native Range

Grape Phylloxera, Daktulosphaira vitifoliae, is a gall-forming insect that feeds on the leaves and roots of many Vitis species. The roots of the cultivated V. vinifera cultivars and hybrids are highly susceptible to grape phylloxera feeding damage. The native range of this insect covers most of North America, and it is particularly abundant in the eastern and central United States. Phylloxera was introduced from North America to almost all grape-growing regions across five of the temperate zone continents. It devastated vineyards in each of these regions causing large-scale disruptions to grape growers, wine makers and national economies. In order to understand the population diversity of grape phylloxera in its native range, more than 500 samples from 19 States and 34 samples from the introduced range (northern California, Europe and South America) were genotyped with 32 simple sequence repeat markers. STRUCTURE, a model based clustering method identified five populations within these samples. The five populations were confirmed by a neighbor-joining tree and principal coordinate analysis (PCoA). These populations were distinguished by their Vitis species hosts and their geographic locations. Samples collected from California, Europe and South America traced back to phylloxera sampled in the northeastern United States on V. riparia, with some influence from phylloxera collected along the Atlantic Coast and Central Plains on V. vulpina. Reproductive statistics conclusively confirmed that sexual reproduction is common in the native range and is combined with cyclical parthenogenesis. Native grape phylloxera populations were identified to be under Hardy-Weinberg equilibrium. The identification of admixed samples between many of these populations indicates that shared environments facilitate sexual reproduction between different host associated populations to create new genotypes of phylloxera. This study also found that assortative mating might occur across the sympatric range of the V. vulpina west and V. cinerea populations.

Comparative transcriptome analysis of chemosensory genes in two sister leaf beetles provides insights into chemosensory speciation

Divergence in chemosensory traits has been posited as an important component of chemosensory speciation in insects. In particular, chemosensory genes expressed in the peripheral sensory neurons are likely to influence insect behaviors such as preference for food, oviposition sites, and mates. Despite their key role in insect behavior and potentially speciation, the underlying genetic basis for divergence in chemosensory traits remains largely unexplored. One way to ascertain the role of chemosensory genes in speciation is to make comparisons of these genes across closely related species to detect the genetic signatures of divergence. Here, we used high throughput transcriptome analysis to compare chemosensory genes of the sister leaf beetles species Pyrrhalta maculicollis and P. aenescens, whose sexual isolation and host plant preference are mediated by divergent chemical signals. Although there was low overall divergence between transcriptome profiles, there were a number of genes that were differentially expressed between the species. Furthermore, we also detected two chemosensory genes under positive selection, one of which that was also differentially expressed between the species, suggesting a possible role for these genes in chemical-based premating reproductive isolation and host use. Combined with the available chemical and ecological work in this system, further studies of the divergent chemosensory genes presented here will provide insight into the process of chemosensory speciation among Pyrrhalta beetles.

Identification and expression analysis of an olfactory receptor gene family in green plant bug Apolygus lucorum (Meyer-Dür)

Olfactory receptors are believed to play a central role in insects host-seeking, mating, and ovipositing. On the basis of male and female antennal transcriptome of adult Apolygus lucorum, a total of 110 candidate A. lucorum odorant receptors (AlucOR) were identified in this study including five previously annotated AlucORs. All the sequences were validated by cloning and sequencing. Tissue expression profiles analysis by RT-PCR indicated most AlucORs were antennal highly expressed genes. The qPCR measurements further revealed 40 AlucORs were significantly higher in the antennae. One AlucOR was primarily expressed in the female antennae, while nine AlucORs exhibited male-biased expression patterns. Additionally, both the RPKM value and RT-qPCR analysis showed AlucOR83 and AlucOR21 were much higher abundant in male antennae than in female antennae, suggesting their different roles in chemoreception of gender. Phylogenetic analysis of ORs from several Hemipteran species demonstrated that most AlucORs had orthologous genes, and five AlucOR-specific clades were defined. In addition, a sub-clade of potential male-based sex pheromone receptors were also identified in the phylogenetic tree of AlucORs. Our results will facilitate the functional studies of AlucORs, and thereby provide a foundation for novel pest management approaches based on these genes.

Gustatory receptors in Lepidoptera: chemosensation and beyond

Lepidoptera is one of the most widespread insect orders and includes several agriculturally important insect species. Ecological success of the lepidopteran insects partly depends on their adaptive chemoreception tactics, which play an important role in the selection of hosts, egg-laying sites and mates. Members of the G-protein coupled receptor family, gustatory receptors (GRs), are an integral part of the Lepidoptera chemosensory machinery. They are expressed in chemosensory neurones and are known to detect different environmental stimuli. Here, we discuss various aspects of the lepidopteran GRs with an emphasis on their roles in different processes such as chemosensation, host selection and adaptation. Phylogenetic analyses have shown that the large diversity of GR genes may have been generated through gene duplication and positive selection events, which also show lineage- and tissue-specific expression. Moreover, lepidopteran GR proteins are diverse and demonstrate broad ligand selectivity for several molecules including sugars, deterrents, salts and CO2 . Binding of ligands to GRs generates multiple downstream changes at the cellular level, which are followed by changes in behaviour. GRs play a critical role in chemosensation and influence the insect's behaviour. Overall, insect GRs are potential targets in the design of effective insect control strategies.

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

Host‐associated races of phytophagous insects provide a model for understanding how adaptation to a new environment can lead to reproductive isolation and speciation, ultimately enabling us to connect barriers to gene flow to adaptive causes of divergence. The pea aphid (Acyrthosiphon pisum) comprises host races specializing on legume species and provides a unique system for examining the early stages of diversification along a gradient of genetic and associated adaptive divergence. As host choice produces assortative mating, understanding the underlying mechanisms of choice will contribute directly to understanding of speciation. As host choice in the pea aphid is likely mediated by smell and taste, we use capture sequencing and SNP genotyping to test for the role of chemosensory genes in the divergence between eight host plant species across the continuum of differentiation and sampled at multiple locations across western Europe. We show high differentiation of chemosensory loci relative to control loci in a broad set of pea aphid races and localities, using a model‐free approach based on principal component analysis. Olfactory and gustatory receptors form the majority of highly differentiated genes and include loci that were already identified as outliers in a previous study focusing on the three most closely related host races. Consistent indications that chemosensory genes may be good candidates for local adaptation and barriers to gene flow in the pea aphid open the way to further investigations aiming to understand their impact on gene flow and to determine their precise functions in response to host plant metabolites.

Differential gene expression according to race and host plant in the pea aphid

Host-race formation in phytophagous insects is thought to provide the opportunity for local adaptation and subsequent ecological speciation. Studying gene expression differences amongst host races may help to identify phenotypes under (or resulting from) divergent selection and their genetic, molecular and physiological bases. The pea aphid (Acyrthosiphon pisum) comprises host races specializing on numerous plants in the Fabaceae and provides a unique system for examining the early stages of diversification along a gradient of genetic and associated adaptive divergence. In this study, we examine transcriptome-wide gene expression both in response to environment and across pea aphid races selected to cover the range of genetic divergence reported in this species complex. We identify changes in expression in response to host plant, indicating the importance of gene expression in aphid-plant interactions. Races can be distinguished on the basis of gene expression, and higher numbers of differentially expressed genes are apparent between more divergent races; these expression differences between host races may result from genetic drift and reproductive isolation and possibly divergent selection. Expression differences related to plant adaptation include a subset of chemosensory and salivary genes. Genes showing expression changes in response to host plant do not make up a large portion of between-race expression differences, providing confirmation of previous studies' findings that genes involved in expression differences between diverging populations or species are not necessarily those showing initial plasticity in the face of environmental change.

Identification and tissue expression profile of genes from three chemoreceptor families in an urban pest, Periplaneta americana

Periplaneta americana is a notorious urban pest prevalent in human habitats; very little is known about its chemosensory mechanism. Employing the advanced next-generation sequencing technique, in the present study, we conducted transcriptome sequencing and analysis of the antennae of the adult males and females as well as their mouthparts using an Illumina platform. This resulted in the discovery of a huge number of the members of all major known chemosensory receptor families in P. americana, including 96 odorant receptors (ORs), 53 ionotropic receptors (IRs), and 33 gustatory receptors (GRs). Tissue expression profiles showed most of them mainly expressed in antennae and phylogenetic analysis demonstrated the expansion in the clade distinguishing them from other functionally well-known Lepidoptera species. A high percentage of chemosensory receptor genes (ORs in particular) showing female antenna bias in mRNA expression was observed. Our results provide a basis for further investigations on how P. americana coordinates its chemosensory receptor genes in chemical communication with environments, and for development of novel pest management approaches.

Targeted resequencing reveals geographical patterns of differentiation for loci implicated in parallel evolution

Parallel divergence and speciation provide evidence for the role of divergent selection in generating biological diversity. Recent studies indicate that parallel phenotypic divergence may not have the same genetic basis in different geographical locations - 'outlier loci' (loci potentially affected by divergent selection) are often not shared among parallel instances of phenotypic divergence. However, limited sharing may be due, in part, to technical issues if false-positive outliers occur. Here, we test this idea in the marine snail Littorina saxatilis, which has evolved two partly isolated ecotypes (adapted to crab predation vs. wave action) in multiple locations independently. We argue that if the low extent of sharing observed in earlier studies in this system is due to sampling effects, we expect outliers not to show elevated FST when sequenced in new samples from the original locations and also not to follow predictable geographical patterns of elevated FST . Following a hierarchical sampling design (within vs. between country), we applied capture sequencing, targeting outliers from earlier studies and control loci. We found that outliers again showed elevated levels of FST in their original location, suggesting they were not generated by sampling effects. Outliers were also likely to show increased FST in geographically close locations, which may be explained by higher levels of gene flow or shared ancestral genetic variation compared with more distant locations. However, in contrast to earlier findings, we also found some outlier types to show elevated FST in geographically distant locations. We discuss possible explanations for this unexpected result.

Molecular Characterization and Sex Distribution of Chemosensory Receptor Gene Family Based on Transcriptome Analysis of Scaeva pyrastri

Chemosensory receptors play key roles in insect behavior. Thus, genes encoding these receptors have great potential for use in integrated pest management. The hover fly Scaeva pyrastri (L.) is an important pollinating insect and a natural enemy of aphids, mainly distributed in the Palearctic and Nearctic regions. However, a systematic identification of their chemosensory receptor genes in the antennae has not been reported. In the present study, we assembled the antennal transcriptome of S. pyrastri by using Illumina sequencing technology. Analysis of the transcriptome data identified 60 candidate chemosensory genes, including 38 for odorant receptors (ORs), 16 for ionotropic receptors (IRs), and 6 for gustatory receptors (GRs). The numbers are similar to those of other Diptera species, suggesting that we were able to successfully identify S. pyrastri chemosensory genes. We analyzed the expression patterns of all genes by using reverse transcriptase PCR (RT-PCR), and found that some genes exhibited sex-biased or sex-specific expression. These candidate chemosensory genes and their tissue expression profiles provide information for further studies aimed at fully understanding the molecular basis behind chemoreception-related behaviors in S. pyrastri.

Variation in life history traits and transcriptome associated with adaptation to diet shifts in the ladybird Cryptolaemus montrouzieri

BACKGROUND

Despite the broad diet range of many predatory ladybirds, the mechanisms involved in their adaptation to diet shifts are not completely understood. Here, we explored how a primarily coccidophagous ladybird Cryptolaemus montrouzieri adapts to feeding on aphids.

RESULTS

Based on the lower survival rate, longer developmental time, and lower adult body weight and reproduction rate of the predator, the aphid Megoura japonica proved being less suitable to support C. montrouzieri as compared with the citrus mealybug Planococcus citri. The results indicated up-regulation of genes related to ribosome and translation in fourth instars, which may be related to their suboptimal development. Also, several genes related to biochemical transport and metabolism, and detoxification were up-regulated as a result of adaptation to the changes in nutritional and non-nutritional (toxic) components of the prey.

CONCLUSION

Our results indicated that C. montrouzieri succeeded in feeding on aphids by regulation of genes related to development, digestion and detoxification. Thus, we argue that these candidate genes are valuable for further studies of the functional evolution of ladybirds led by diet shifts.

Targeted capture in evolutionary and ecological genomics

The rapid expansion of next-generation sequencing has yielded a powerful array of tools to address fundamental biological questions at a scale that was inconceivable just a few years ago. Various genome-partitioning strategies to sequence select subsets of the genome have emerged as powerful alternatives to whole-genome sequencing in ecological and evolutionary genomic studies. High-throughput targeted capture is one such strategy that involves the parallel enrichment of preselected genomic regions of interest. The growing use of targeted capture demonstrates its potential power to address a range of research questions, yet these approaches have yet to expand broadly across laboratories focused on evolutionary and ecological genomics. In part, the use of targeted capture has been hindered by the logistics of capture design and implementation in species without established reference genomes. Here we aim to (i) increase the accessibility of targeted capture to researchers working in nonmodel taxa by discussing capture methods that circumvent the need of a reference genome, (ii) highlight the evolutionary and ecological applications where this approach is emerging as a powerful sequencing strategy and (iii) discuss the future of targeted capture and other genome-partitioning approaches in the light of the increasing accessibility of whole-genome sequencing. Given the practical advantages and increasing feasibility of high-throughput targeted capture, we anticipate an ongoing expansion of capture-based approaches in evolutionary and ecological research, synergistic with an expansion of whole-genome sequencing.

Molecular Characterization and Differential Expression of an Olfactory Receptor Gene Family in the White-Backed Planthopper Sogatella furcifera Based on Transcriptome Analysis

The white-backed planthopper, Sogatella furcifera, a notorious rice pest in Asia, employs host plant volatiles as cues for host location. In insects, odor detection is mediated by two types of olfactory receptors: odorant receptors (ORs) and ionotropic receptors (IRs). In this study, we identified 63 SfurORs and 14 SfurIRs in S. furcifera based on sequences obtained from the head transcriptome and bioinformatics analysis. The motif-pattern of 130 hemiptera ORs indicated an apparent differentiation in this order. Phylogenetic trees of the ORs and IRs were constructed using neighbor-joining estimates. Most of the ORs had orthologous genes, but a specific OR clade was identified in S. furcifera, which suggests that these ORs may have specific olfactory functions in this species. Our results provide a basis for further investigations of how S. furcifera coordinates its olfactory receptor genes with its plant hosts, thereby providing a foundation for novel pest management approaches based on these genes.