De novo comparative transcriptome analysis of a rare cicada, with identification of candidate genes related to adaptation to a novel host plant and drier habitats

The impact of geographic isolation and host shifts on population divergence of the rare cicada Subpsaltria yangi

The contributions of divergent selection and spatial isolation to population divergence are among the main focuses of evolutionary biology. Here we employed integrated methods to explore genomic divergence, demographic history and calling-song differentiation in the cicada Subpsaltria yangi, and compared the genotype and calling-song phenotype of different populations occurring in distinct habitats. Our results indicate that this species comprises four main lineages with unique sets of haplotypes and calling-song structure, which are distinctly associated with geographic isolation and habitats. The populations occurring on the Loess Plateau underwent substantial expansion at ∼0.130-0.115 Ma during the Last Interglacial. Geographic distance and host shift between pairs of populations predict genomic divergence, with geographic distance and acoustical signal together explaining > 60% of the divergence among populations. Differences in calling songs could reflect adaptation of populations to novel environments with different host plants, habitats and predators, which may have resulted from neutral divergence at the molecular level followed by natural selection. Geomorphic barriers and climate oscillations associated with Pleistocene glaciation may have been primary factors in shaping the population genetic structure of this species. Ultimately this may couple with a host shift in leading toward allopatric speciation in S. yangi, i.e., isolation by distance. Our findings improve understanding of divergence in allopatry of herbivorous insects, and may inform future studies on the molecular mechanisms underlying the association between genetic/phenotypic changes and adaptation of insects to novel niches and host plants.

Gene expression plasticity facilitates different host feeding in Ips sexdentatus (Coleoptera: Curculionidae: Scolytinae)

Host shift is ecologically advantageous and a crucial driver for herbivore insect speciation. Insects on the non-native host obtain enemy-free space and confront reduced competition, but they must adapt to survive. Such signatures of adaptations can often be detected at the gene expression level. It is astonishing how bark beetles cope with distinct chemical environments while feeding on various conifers. Hence, we aim to disentangle the six-toothed bark beetle (Ips sexdentatus) response against two different conifer defences upon host shift (Scots pine to Norway spruce). We conducted bioassay and metabolomic analysis followed by RNA-seq experiments to comprehend the beetle's ability to surpass two different terpene-based conifer defence systems. Beetle growth rate and fecundity were increased when reared exclusively on spruce logs (alternative host) compared to pine logs (native host). Comparative gene expression analysis identified differentially expressed genes (DEGs) related to digestion, detoxification, transporter activity, growth, signalling, and stress response in the spruce-feeding beetle gut. Transporter genes were highly abundant during spruce feeding, suggesting they could play a role in pumping a wide variety of endogenous and xenobiotic compounds or allelochemicals out. Trehalose transporter (TRET) is also up-regulated in the spruce-fed beetle gut to maintain homeostasis and stress tolerance. RT-qPCR and enzymatic assays further corroborated some of our findings. Taken together, the transcriptional plasticity of key physiological genes plays a crucial role after the host shift and provides vital clues for the adaptive potential of bark beetles on different conifer hosts.

Population Genomics for Insect Conservation

Insects constitute vital components of ecosystems. There is alarming evidence for global declines in insect species diversity, abundance, and biomass caused by anthropogenic drivers such as habitat degradation or loss, agricultural practices, climate change, and environmental pollution. This raises important concerns about human food security and ecosystem functionality and calls for more research to assess insect population trends and identify threatened species and the causes of declines to inform conservation strategies. Analysis of genetic diversity is a powerful tool to address these goals, but so far animal conservation genetics research has focused strongly on endangered vertebrates, devoting less attention to invertebrates, such as insects, that constitute most biodiversity. Insects' shorter generation times and larger population sizes likely necessitate different analytical methods and management strategies. The availability of high-quality reference genome assemblies enables population genomics to address several key issues. These include precise inference of past demographic fluctuations and recent declines, measurement of genetic load levels, delineation of evolutionarily significant units and cryptic species, and analysis of genetic adaptation to stressors. This enables identification of populations that are particularly vulnerable to future threats, considering their potential to adapt and evolve. We review the application of population genomics to insect conservation and the outlook for averting insect declines.

De novo transcriptomes for two species of snappers from the Western Atlantic (Lutjanus griseus and Lutjanus synagris)

Transcriptomes are appropriate resources for studying species that lack sequenced genomes, as they can serve as references for a broad suite of genetic applications, including: phylogenetic assessments, population genomics, and evaluate responses to environmental fluctuations. Here, we present the transcriptomes of two species of marine fishes of commercial and ecological relevance in the Western Atlantic: Lutjanus griseus and L. synagris. This project represents a step forward on developing genomic resources for important species of the Atlantic Ocean.

Role of Genes in Regulating Host Plants Expansion in Tephritid Fruit Flies (Diptera) and Potential for RNAi-Based Control

Host plant expansion is an important survival strategy for tephritids as they expand their range. Successful host expansion requires tephritids to adapt to the chemical and nonchemical properties of a novel host fruit, such as fruit color, phenology, and phytochemicals. These plant properties trigger a series of processes in tephritids, with each process having its own genetic basis, which means that various genes are involved in regulating host plant expansion by tephritids. This review summarizes current knowledge on the categories and roles of genes involved in host plant expansion in several important tephritid species, including genes related to chemoreception (olfactory and gustation), vision, digestion, detoxification, development, ribosomal and energy metabolism. Chemoreception- and detoxification- and digestion-related genes are stimulated by volatile chemicals and secondary chemicals of different hosts, respectively, which are involved in the regulation of nervous signal transduction that triggers behavioral, physical, and chemical responses to the novel host fruit. Vision-, nerve-, and development-related genes and metabolism-associated genes are activated in response to nonchemical stimuli from different hosts, such as color and phenology, to regulate a comprehensive adaptation of the extending host for tephritids. The chemical and nonchemical signals of hosts activate ribosomal and energy-related genes that result in the basic regulation of many processes of host expansion, including detoxification and development. These genes do not regulate novel host use individually, but multiple genes regulate multilevel adaptation to novel host fruits via multiple mechanisms. These genes may also be potential target genes for RNAi-based control of tephritid pests.

Transcriptional Regulation of Reproductive Diapause in the Convergent Lady Beetle, Hippodamia convergens

Simple Summary

Diapause is a dormant period typically controlled by daylength that ensures an insect’s survival through harsh environmental conditions. The convergent lady beetle, Hippodamia convergens, undergoes a reproductive diapause in winter, where female ovaries remain immature and no eggs are laid. This species is an important biological control agent, but during diapause, beetles are less likely to eat pest insects. Thus, knowledge of diapause mechanisms may facilitate manipulation thereof to improve biological control. Further, molecular studies of adult diapause and diapause in Coleoptera are relatively lacking. Here, we assembled and annotated a transcriptome for this species and quantified transcript expression changes during diapause. Female beetles were sampled at three times in diapause (early, mid, and late diapause), which allowed us to characterize the molecular processes occurring at distinct transitions throughout diapause. We found that transcripts involved in flight were consistently upregulated during diapause, which is consistent with dispersal flights at this stage, while transcripts involved in ovarian development were downregulated, which is consistent with the shutdown of reproduction in diapausing females. These findings identify key regulators of diapause in H. convergens and contribute to a growing body of literature on the molecular mechanisms of diapause across the insect phylogeny.

Abstract

Diapause is an alternate development program that synchronizes an insect’s life cycle with seasonally abundant resources and ensures survival in unfavorable conditions. The physiological basis of diapause has been well characterized, but the molecular mechanisms regulating it are still being elucidated. Here, we present a de novo transcriptome and quantify transcript expression during diapause in the convergent lady beetle Hippodamia convergens. H. convergens is used as an augmentative biocontrol agent, and adult females undergo reproductive diapause that is regulated by photoperiod. We sampled females at three stages (early, mid, and late diapause) and compared transcript expression to non-diapausing individuals. Based on principle component analysis, the transcriptomes of diapausing beetles were distinct from non-diapausing beetles, and the three diapausing points tended to cluster together. However, there were still classes of transcripts that differed in expression across distinct phases of diapause. In general, transcripts involved in muscle function and flight were upregulated during diapause, likely to support dispersal flights that occur during diapause, while transcripts involved in ovarian development were downregulated. This information could be used to improve biological control by manipulating diapause. Additionally, our data contribute to a growing understanding of the genetic regulation of diapause across diverse insects.

Genomic insight into the scale specialization of the biological control agent Novius pumilus (Weise, 1892)

BACKGROUND

Members of the genus Novius Mulsant, 1846 (= Rodolia Mulsant, 1850) (Coleoptera, Coccinellidae), play important roles in the biological control of cotton cushion scale pests, especially those belonging to Icerya. Since the best-known species, the vedalia beetle Novius cardinalis (Mulsant, 1850) was introduced into California from Australia, more than a century of successful use in classical biological control, some species of Novius have begun to exhibit some field adaptations to novel but related prey species. Despite their economic importance, relatively little is known about the underlying genetic adaptations associated with their feeding habits. Knowledge of the genome sequence of Novius is a major step towards further understanding its biology and potential applications in pest control.

RESULTS

We report the first high-quality genome sequence for Novius pumilus (Weise, 1892), a representative specialist of Novius. Computational Analysis of gene Family Evolution (CAFE) analysis showed that several orthogroups encoding chemosensors, digestive, and immunity-related enzymes were significantly expanded (P < 0.05) in N. pumilus compared to the published genomes of other four ladybirds. Furthermore, some of these orthogroups were under significant positive selection pressure (P < 0.05). Notably, transcriptome profiling demonstrated that many genes among the significantly expanded and positively selected orthogroups, as well as genes related to detoxification were differentially expressed, when N. pumilus feeding on the nature prey Icerya compared with the no feeding set. We speculate that these genes are vital in the Icerya adaptation of Novius species.

CONCLUSIONS

We report the first Novius genome thus far. In addition, we provide comprehensive transcriptomic resources for N. pumilus. The results from this study may be helpful for understanding the association of the evolution of genes related to chemosensing, digestion, detoxification and immunity with the prey adaptation of insect predators. This will provide a reference for future research and utilization of Novius in biological control programs. Moreover, understanding the possible molecular mechanisms of prey adaptation also inform mass rearing of N. pumilus and other Novius, which may benefit pest control.

Gene expression profiles in Malpighian tubules of the vector leafhopper Psammotettix striatus (L.) revealed regional functional diversity and heterogeneity

Background

Many leafhoppers are known as pests and disease vectors of economically important plants. Previous studies of the physiological functions of vector leafhoppers have mainly focused on the salivary glands and the alimentary tract that are deemed to be associated with digestion, host defense and phytoplasma and/or virus transmission. By contrast, the significance of Malpighian tubules (MTs) is less studied. To clarify the physiological function of MTs of the vector leafhopper Psammotettix striatus that transmits phytoplasma triggering the wheat blue dwarf disease, we performed a transcriptome study on P. striatus MTs and compared gene expression profiles among different anatomical regions in the tubules (i.e., MT1+2, the anterior segment together with the sub-anterior segment; MT3, the inflated segment; and MT4, the distal segment).

Results

Transcriptome of P. striatus MTs generate a total of 42,815 high-quality unigenes, among which highly expressed unigenes are mainly involved in organic solute transport, detoxification and immunity in addition to osmoregulation. Region-specific comparative analyses reveal that all these MT regions have functions in osmoregulation, organic solute transport and detoxification, but each region targets different substrates. Differential expression and regional enrichment of immunity-related effector activities and molecules involved in phagocytosis and the biosynthesis of antimicrobial peptides among different regions indicate that MT1+2 and MT4 have the ability to eliminate the invading pathogens. However, in MT3 which secrets brochosomes to the integument and eggs as physical barriers, disulfide-isomerase, acidic ribosomal protein P and many other unigenes were highly expressed, which can be attractive candidate genes for future studies of the biosynthesis and the origin of brochosomes.

Conclusions

Psammotettix striatus MTs perform multiple physiological functions as versatile organs than just excretory organs with osmoregulatory function. Heterogeneity of physiological functions among different MT regions is related to organic solute transport, detoxification, immunity and brochosome biosynthesis in addition to osmoregulation, and each region targets different substrates. These functions may be helpful for P. striatus to resist pathogens from habitats and to utilize a wider range of host plants, which may assist the transmission and spread of phytoplasmas. The results provide potential molecular targets for the exploit of chemical and/or gene-silencing insecticides.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08300-6.

Behavioral and Physiological Plasticity Provides Insights into Molecular Based Adaptation Mechanism to Strain Shift in Spodoptera frugiperda

How herbivorous insects adapt to host plants is a key question in ecological and evolutionary biology. The fall armyworm, (FAW) Spodoptera frugiperda (J.E. Smith), although polyphagous and a major pest on various crops, has been reported to have a rice and corn (maize) feeding strain in its native range in the Americas. The species is highly invasive and has recently established in China. We compared behavioral changes in larvae and adults of a corn population (Corn) when selected on rice (Rice) and the molecular basis of these adaptational changes in midgut and antennae based on a comparative transcriptome analysis. Larvae of S. frugiperda reared on rice plants continuously for 20 generations exhibited strong feeding preference for with higher larval performance and pupal weight on rice than on maize plants. Similarly, females from the rice selected population laid significantly more eggs on rice as compared to females from maize population. The most highly expressed DEGs were shown in the midgut of Rice vs. Corn. A total of 6430 DEGs were identified between the populations mostly in genes related to digestion and detoxification. These results suggest that potential adaptations for feeding on rice crops, may contribute to the current rapid spread of fall armyworm on rice crops in China and potentially elsewhere. Consistently, highly expressed DEGs were also shown in antennae; a total of 5125 differentially expressed genes (DEGs) s were identified related to the expansions of major chemosensory genes family in Rice compared to the Corn feeding population. These results not only provide valuable insight into the molecular mechanisms in host plants adaptation of S. frugiperda but may provide new gene targets for the management of this pest.

Females prefer males producing a high-rate song with shorter timbal–stridulatory sound intervals in a cicada species

Uncovering mate choice and factors that lead to the choice are very important to understanding sexual selection in evolutionary change. Cicadas are known for their loud sounds produced by males using the timbals. However, males in certain cicada species emit 2 kinds of sounds using respectively timbals and stridulatory organs, and females may produce their own sounds to respond to males. What has never been considered is the mate choice in such cicada species. Here, we investigate the sexual selection and potential impact of predation pressure on mate choice in the cicada Subpsaltria yangi Chen. It possesses stridulatory sound-producing organs in both sexes in addition to the timbals in males. Results show that males producing calling songs with shorter timbal–stridulatory sound intervals and a higher call rate achieved greater mating success. No morphological traits were found to be correlated with mating success in both sexes, suggesting neither males nor females display mate preference for the opposite sex based on morphological traits. Males do not discriminate among responding females during mate searching, which may be due to the high energy costs associated with their unusual mate-seeking activity and the male-biased predation pressure. Females generally mate once but a minority of them re-mated after oviposition which, combined with the desirable acoustic traits of males, suggest females may maximize their reproductive success by choosing a high-quality male in the first place. This study contributes to our understanding mechanisms of sexual selection in cicadas and other insects suffering selective pressure from predators.

Comparative transcriptome analysis of the newly discovered insect vector of the pine wood nematode in China, revealing putative genes related to host plant adaptation

BACKGROUND

In many insect species, the larvae/nymphs are unable to disperse far from the oviposition site selected by adults. The Sakhalin pine sawyer Monochamus saltuarius (Gebler) is the newly discovered insect vector of the pine wood nematode (Bursaphelenchus xylophilus) in China. Adult M. saltuarius prefers to oviposit on the host plant Pinus koraiensis, rather than P. tabuliformis. However, the genetic basis of adaptation of the larvae of M. saltuarius with weaken dispersal ability to host environments selected by the adult is not well understood.

RESULTS

In this study, the free amino and fatty acid composition and content of the host plants of M. saltuarius larvae, i.e., P. koraiensis and P. tabuliformis were investigated. Compared with P. koraiensis, P. tabuliformis had a substantially higher content of various free amino acids, while the opposite trend was detected for fatty acid content. The transcriptional profiles of larval populations feeding on P. koraiensis and P. tabuliformis were compared using PacBio Sequel II sequencing combined with Illumina sequencing. The results showed that genes relating to digestion, fatty acid synthesis, detoxification, oxidation-reduction, and stress response, as well as nutrients and energy sensing ability, were differentially expressed, possibly reflecting adaptive changes of M. saltuarius in response to different host diets. Additionally, genes coding for cuticle structure were differentially expressed, indicating that cuticle may be a potential target for plant defense. Differential regulation of genes related to the antibacterial and immune response were also observed, suggesting that larvae of M. saltuarius may have evolved adaptations to cope with bacterial challenges in their host environments.

CONCLUSIONS

The present study provides comprehensive transcriptome resource of M. saltuarius relating to host plant adaptation. Results from this study help to illustrate the fundamental relationship between transcriptional plasticity and adaptation mechanisms of insect herbivores to host plants.

Mechanisms of Novel Host Use by Bactrocera tau (Tephritid: Diptera) Revealed by RNA Transcriptomes

Use of novel plant hosts can facilitate the establishment and range expansion of herbivorous invasive species. However, the inherent mechanisms of novel host use are still unclear in many herbivorous species. Here, we examine mechanisms of novel host use in the invasive tephritid fruit fly Bactrocera tau (Walker)(Diptera: Tephritidae) by documenting changes in the RNA transcriptomes associated with a novel host. RNA transcripts of B. tau were obtained with high-throughput sequencing from samples continuously reared on two traditional Cucurbitaceae hosts and a novel host (banana). We found transcriptome variation was strongly associated with feeding on banana. Moreover, B. tau feeding on banana contained more differentially expressed genes (DEGs) and more annotated categories of DEGs in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database with 1,595 DEGs and 21 major annotated pathways. The annotated categories of DEGs in individuals reared on banana differed with from those individuals feeding on other hosts and were enriched in oxidative phosphorylation, citrate cycle pathway, and four other carbohydrate pathways. For B. tau feeding on banana, the predominant numbers of upregulated genes in the mitochondrial NADH (56 on average) and a relatively higher numbers of upregulated genes (13 on average) were found in oxidative phosphorylation and the TCA pathway, respectively. Changes in RNA transcriptomes associated with novel host use, especially for genes related to energy and carbohydrate metabolism, help to explain how B. tau can be successful in use of novel hosts and may be useful in developing novel strategies for control of tephritid flies.

Transcriptome Profiling Revealed Potentially Critical Roles for Digestion and Defense-Related Genes in Insects' Use of Resistant Host Plants: A Case Study with Sitobion Avenae

Using host plant resistance (HPR) in management of insect pests is often environmentally friendly and suitable for sustainable development of agricultural industries. However, this strategy can be limited by rapid evolution of insect populations that overcome HPR, for which the underlying molecular factors and mechanisms are not well understood. To address this issue, we analyzed transcriptomes of two distinct biotypes of the grain aphid, Sitobion avenae (Fabricius), on wheat and barley. This analysis revealed a large number of differentially expressed genes (DEGs) between biotypes 1 and 3 on wheat and barley. The majority of them were common DEGs occurring on both wheat and barley. GO and KEGG enrichment analyses for these common DEGs demonstrated significant expression divergence between both biotypes in genes associated with digestion and defense. Top defense-related common DEGs with the most significant expression changes included three peroxidases, two UGTs (UDP-glycosyltransferase), two cuticle proteins, one glutathione S-transferases (GST), one superoxide dismutase, and one esterase, suggesting their potentially critical roles in the divergence of S. avenae biotypes. A relatively high number of specific DEGs on wheat were identified for peroxidases (9) and P450s (8), indicating that phenolic compounds and hydroxamic acids may play key roles in resistance of wheat against S. avenae. Enrichment of specific DEGs on barley for P450s and ABC transporters suggested their key roles in this aphid's detoxification against secondary metabolites (e.g., alkaloids) in barley. Our results can provide insights into the molecular factors and functions that explain biotype adaptation in insects and their use of resistant plants. This study also has significant implications for developing new resistant cultivars, developing strategies that limit rapid development of insect biotypes, and extending resistant crop cultivars' durability and sustainability in integrated management programs.

Candidate detoxification-related genes in brown planthopper, Nilaparvata lugens, in response to β-asarone based on transcriptomic analysis

Nilaparvata lugens(Stål) is a serious pest of rice and has evolved different levels of resistance against most chemical pesticides. β-asarone is the main bioactive insecticidal compound of Acorus calamus L. that shows strong insecticidal activity against pests. In this study, we conducted a bioassay experiment to determine the contact toxicity of β-asarone to N. lugens nymphs. The LD₃₀ sublethal dose was 0.106 μg per nymph, with 95% confidence limits of 0.070-0.140 μg. We applied the LD₃₀ concentration of β-asarone to nymphs for 24 h or 72 h and then performed a transcriptome sequence analysis by referencing the N. lugens genome to characterize the variation. The transcriptomic analysis showed that several GO terms and KEGG pathways presented significant changes. Individually, 126 differentially expressed genes (DEGs), including 72 upregulated and 54 downregulated genes, were identified at 24 h, and 1771 DEGs, including 882 upregulated and 889 downregulated genes, were identified at 72 h. From the DEGs, we identified a total of 40 detoxification-related genes, including eighteen Cytochrome P450 monooxygenase genes (P450s), three Glutathione S-transferase genes, one Carboxylesterase gene, twelve UDP-glucosyltransferases and six ATP-binding cassette genes. We selected the eighteen P450s for subsequent verification by quantitative PCR. These findings indicated that β-asarone presented strong contact toxicity to N. lugens nymphs and induced obvious variation of detoxification-related genes that may be involved in the response to β-asarone.

Gene Expression and Diet Breadth in Plant-Feeding Insects: Summarizing Trends

Transcriptomic studies lend insights into the role of transcriptional plasticity in adaptation and specialization. Recently, there has been growing interest in understanding the relationship between variation in herbivorous insect gene expression and the evolution of diet breadth. We review the studies that have emerged on insect gene expression and host plant use, and outline the questions and approaches in the field. Many candidate genes underlying herbivory and specialization have been identified, and a few key studies demonstrate increased transcriptional plasticity associated with generalist compared with specialist species. Addressing the roles that transcriptional variation plays in insect diet breadth will have important implications for our understanding of the evolution of specialization and the genetic and environmental factors that govern insect-plant interactions.

Bacterial communities in digestive and excretory organs of cicadas

Bacteriocyte-associated symbionts are essential for the health of many sap-sucking insects, such as cicadas, leafhoppers and treehoppers, etc., but little is known about the bacterial community in the gut and other related organs in these insects. We characterized the bacterial communities in the salivary glands, alimentary canal and the Malpighian tubules of two populations of the cicada Subpsaltria yangi occurring in different habitats and feeding on different hosts. A high degree of similarity of core microbiota was revealed between the two populations, both with the top three bacteria belonging to Meiothermus, Candidatus Sulcia and Halomonas. The bacterial communities in various organs clustered moderately by populations possibly reflect adaptive changes in the microbiota of related S. yangi populations, which provide a better understanding of the speciation and adaptive mechanism of this species to different diets and habitats. When compared with two phylogenetically distant cicada species, Hyalessa maculaticollis and Meimuna mongolica, the core microbiota in S. yangi was significantly different to that of these species. In addition, our results confirm that Ca. Sulcia distributes in the digestive and excretory organs besides the bacteriomes and gonads, which provide potential important information onto the trophic functions of this obligate endosymbiont to the host insects.

Off-target capture data, endosymbiont genes and morphology reveal a relict lineage that is sister to all other singing cicadas

Phylogenetic asymmetry is common throughout the tree of life and results from contrasting patterns of speciation and extinction in the paired descendant lineages of ancestral nodes. On the depauperate side of a node, we find extant ‘relict’ taxa that sit atop long, unbranched lineages. Here, we show that a tiny, pale green, inconspicuous and poorly known cicada in the genus Derotettix, endemic to degraded salt-plain habitats in arid regions of central Argentina, is a relict lineage that is sister to all other modern cicadas. Nuclear and mitochondrial phylogenies of cicadas inferred from probe-based genomic hybrid capture data of both target and non-target loci and a morphological cladogram support this hypothesis. We strengthen this conclusion with genomic data from one of the cicada nutritional bacterial endosymbionts, Sulcia, an ancient and obligate endosymbiont of the larger plant-sucking bugs (Auchenorrhyncha) and an important source of maternally inherited phylogenetic data. We establish Derotettiginae subfam. nov. as a new, monogeneric, fifth cicada subfamily, and compile existing and new data on the distribution, ecology and diet of Derotettix. Our consideration of the palaeoenvironmental literature and host-plant phylogenetics allows us to predict what might have led to the relict status of Derotettix over 100 Myr of habitat change in South America.

Comparative sialotranscriptome analysis of the rare Chinese cicada Subpsaltria yangi, with identification of candidate genes related to host-plant adaptation

Little is known about the mechanisms underlying the relationship between changes in salivary proteins and the adaptation of insects to different host-plants. To address this knowledge gap, the transcriptional profiles of salivary glands were compared among three populations of the rare cicada Subpsaltria yangi, in which two populations specialize on Zizyphus jujuba var. spinosa, but the population occurring in the Helan (HL) Mountains is locally specialized on the endemic plant Ephedra lepidosperma. The comparisons indicate that genes related to digestion and detoxification are differentially regulated in populations feeding on different plants, possibly reflecting adaptative changes in salivary proteins of S. yangi in response to different host chemistries. In detail, 38 differentially expressed genes and 21 up-regulated genes related to digestion and detoxification were identified respectively in two pairwise comparisons among the populations using different hosts, with some genes exclusively expressed in the HL population. Our results are consistent with the hypothesis that the host plant shift in the HL population was facilitated by differential regulation of genes related to digestion and detoxification. This study provides new information for elucidating the molecular mechanisms underlying the relationship between changed salivary proteins and the adaptability of plant-feeding insects to novel host plants.