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

Arsenophonus Interacts with Buchnera to Improve Growth Performance of Aphids under Amino Acid Stress

Amino acids play a crucial role in the growth and development of insects. Aphids cannot ingest enough amino acids in plant phloem to meet their requirements, and therefore, they are mainly dependent on the obligate symbiont Buchnera aphidicola to synthesize essential amino acids. Besides Buchnera, aphids may harbor another facultative symbiont, Arsenophonus, which alters the requirement of the cotton-melon aphid Aphis gossypii for amino acid. However, it is unclear how Arsenophonus regulates the requirement. Here, we found that Arsenophonus ameliorated growth performance of A. gossypii on an amino acid-deficient diet. A deficiency in lysine (Lys) or methionine (Met) led to changes in the abundance of Arsenophonus. Arsenophonus suppressed the abundance of Buchnera when aphids were fed a normal amino acid diet, but this suppression was eliminated or reversed when aphids were on a Lys- or Met-deficient diet. The relative abundance of Arsenophonus was positively correlated with that of Buchnera, but neither of them was correlated with the body weight of aphids. The relative expression levels of Lys and Met synthase genes of Buchnera were affected by the interaction between Arsenophonus infections and Buchnera abundance, especially in aphids reared on a Lys- or Met-deficient diet. Arsenophonus coexisted with Buchnera in bacteriocytes, which strengthens the interaction. IMPORTANCE The obligate symbiont Buchnera can synthesize amino acids for aphids. In this study, we found that a facultative symbiont, Arsenophonus, can help improve aphids' growth performance under amino acid deficiency stress by changing the relative abundance of Buchnera and the expression levels of amino acid synthase genes. This study highlights the interaction between Arsenophonus and Buchnera to ameliorate aphid growth under amino acid stress.

Buchnera breaks the specialization of the cotton-specialized aphid (Aphis gossypii) by providing nutrition through zucchini

The cotton aphid, Aphis gossypii Glover, is a species of polyphagous aphid with many biotypes, and its host transfer has always been the focus of research on the control of cotton aphid. An important factor affecting aphid specialization is the nutritional association with microbial symbionts that provide the host with nutrients lacking in the diet. We analyzed the microbial composition and biodiversity of reared on zucchini for 10 generations (T1-T10) and cotton as a control (CK), by high-throughput Illumina sequencing of 16S ribosomal RNA genes. The findings showed that the change in plant hosts decreased the richness and variety of microbial species. Regardless of whether the plant host is altered or not, Proteobacteria and Firmicutes are the predominate phyla in cotton-specialized aphid. Additionally, cotton-specialized aphids that live in zucchini had considerably lower relative abundances of non-dominant phyla (Bacteroidetes) than cotton hosts. At the genus level the dominant communities were Buchnera, Acinetobacter, and Arsenophonus. The relative abundance of Buchnera was significantly higher in aphids reared on zucchini than those on cotton, whereas the opposite was observed for Acinetobacter, as well as for some non-dominant communities (Stenotrophomonas, Pseudomons, Flavobacterium, Novosphingobium). Collectively, this study clarifies the dynamic changes of symbiotic bacteria in cotton-specialized aphids reared on zucchini for multiple generations. Among them, Buchnera is crucial for the cotton-specialized aphid to get nutrients during the transfer of the host and has a favorable impact on the colonization of cotton-specialized aphid populations on zucchini hosts. It not only enriches our understanding of the relationship between the bacterial microbiota of aphids and their adaptability to new hosts, zucchini, but also expands the current body of research on the mechanisms underlying the host shifting ability of cotton-specialized aphids.

Silencing of Cytochrome P450 Gene AgoCYP6CY19 Reduces the Tolerance to Host Plant in Cotton- and Cucumber-Specialized Aphids, Aphis gossypii

Cytochrome P450 monooxygenases play important roles in insect metabolism and detoxification of toxic plant substances. However, the function of CYP6 family genes in degrading plant toxicants in Aphis gossypii has yet to be elucidated. In this study, AgoCYP6CY19, an A. gossypii CYP gene that differentially expresses in cotton- and cucumber-specialized aphids, was characterized. Spatiotemporal expression profiling revealed that AgoCYP6CY19 expression was higher in second instar nymph and 7 day old adults than in other developmental stages. Although the expression of AgoCYP6CY19 was significantly higher in cotton-specialized aphids, AgoCYP6CY19 silencing significantly increased larval and adult mortality and reduced total fecundity in both cotton- and cucumber-specialized aphids. What is more, the expression of AgoCYP6CY19 was significantly induced after the cotton-specialized and cucumber-specialized aphids fed on epigallocatechin gallate (EGCG) and cucurbitacin B (CucB), respectively. These findings demonstrate that AgoCYP6CY19 plays a pivotal role in toxic plant substance detoxification and metabolism. Functional knowledge about plant toxicity tolerance genes in this major pest can provide new insights into insect detoxification of toxic plant substances and insecticides and offer new targets for agricultural pest control strategies.

The molecular mechanisms that determine different degrees of polyphagy in the Bemisia tabaci species complex

The whitefly Bemisia tabaci is a closely related group of >35 cryptic species that feed on the phloem sap of a broad range of host plants. Species in the complex differ in their host-range breadth, but the mechanisms involved remain poorly understood. We investigated, therefore, how six different B. tabaci species cope with the environmental unpredictability presented by a set of four common and novel host plants. Behavioral studies indicated large differences in performances on the four hosts and putative specialization of one of the species to cassava plants. Transcriptomic analyses revealed two main insights. First, a large set of genes involved in metabolism (>85%) showed differences in expression between the six species, and each species could be characterized by its own unique expression pattern of metabolic genes. However, within species, these genes were constitutively expressed, with a low level of environmental responsiveness (i.e., to host change). Second, within each species, sets of genes mainly associated with the super-pathways "environmental information processing" and "organismal systems" responded to the host switching events. These included genes encoding for proteins involved in sugar homeostasis, signal transduction, membrane transport, and immune, endocrine, sensory and digestive responses. Our findings suggested that the six B. tabaci species can be divided into four performance/transcriptomic "Types" and that polyphagy can be achieved in multiple ways. However, polyphagy level is determined by the specific identity of the metabolic genes/pathways that are enriched and overexpressed in each species (the species' individual metabolic "tool kit").

Elevated rates of positive selection drive the evolution of pestiferousness in the Colorado potato beetle (Leptinotarsa decemlineata, Say)

Contextualizing evolutionary history and identifying genomic features of an insect that might contribute to its pest status is important in developing early detection and control tactics. In order to understand the evolution of pestiferousness, which we define as the accumulation of traits that contribute to an insect population's success in an agroecosystem, we tested the importance of known genomic properties associated with rapid adaptation in the Colorado potato beetle (CPB), Leptinotarsa decemlineata Say. Within the leaf beetle genus Leptinotarsa, only CPB, and a few populations therein, has risen to pest status on cultivated nightshades, Solanum. Using whole genomes from ten closely related Leptinotarsa species native to the United States, we reconstructed a high-quality species tree and used this phylogenetic framework to assess evolutionary patterns in four genomic features of rapid adaptation: standing genetic variation, gene family expansion and contraction, transposable element abundance and location, and positive selection at protein-coding genes. Throughout approximately 20 million years of history, Leptinotarsa species show little evidence of gene family turnover and transposable element variation. However, there is a clear pattern of CPB experiencing higher rates of positive selection on protein-coding genes. We determine that these rates are associated with greater standing genetic variation due to larger effective population size, which supports the theory that the demographic history contributes to rates of protein evolution. Furthermore, we identify a suite of coding genes under positive selection that are putatively associated with pestiferousness in the Colorado potato beetle lineage. They are involved in the biological processes of xenobiotic detoxification, chemosensation and hormone function.

Transcriptome profiling revealed potentially important roles of defensive gene expression in the divergence of insect biotypes: a case study with the cereal aphid Sitobion avenae

Background

Many insects can develop differential biotypes on variable host plants, but the underlying molecular factors and mechanisms are not well understood. To address this issue, transcriptome profiling analyses were conducted for two biotypes of the cereal aphid, Sitobion avenae (Fabricius), on both original and alternative plants.

Results

Comparisons between both biotypes generated 4174 differentially expressed unigenes (DEGs). In their response to host plant shift, 39 DEGs were shared by both biotypes, whereas 126 and 861 DEGs occurred only in biotypes 1 and 3, respectively. MMC (modulated modularity clustering) analyses showed that specific DEGs of biotypes 1 and 3 clustered into five and nine transcriptional modules, respectively. Among these DEGs, defense-related genes underwent intensive expression restructuring in both biotypes. However, biotype 3 was found to have relatively lower gene transcriptional plasticity than biotype 1. Gene enrichment analyses of the abovementioned modules showed functional divergence in defensive DEGs for the two biotypes in response to host transfer. The expression plasticity for some defense related genes was showed to be directly related to fecundity of S. avenae biotypes on both original and alternative plants, suggesting that expression plasticity of key defensive genes could have significant impacts on the adaptive potential and differentiation of S. avenae biotypes on different plants.

Conclusions

The divergence patterns of transcriptional plasticity in defense related genes may play important roles in the phenotypic evolution and differentiation of S. avenae biotypes. Our results can provide insights into the role of gene expression plasticity in the divergence of insect biotypes and adaptive evolution of insect populations.

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.

Transcriptional Responses in Defense-Related Genes of Sitobion avenae (Hemiptera: Aphididae) Feeding on Wheat and Barley

Many aphids can adapt on plants of differential resistance levels, but molecular interactions underlying host plant utilization and shift of aphids are still not well understood. Here, we sequenced the transcriptome and compared global gene expression profiles of Sitobion avenae (Fabricius) feeding on wheat (i.e., the susceptible plant) and barley (i.e., the resistant plant). In total, 65,200 high-quality unigenes were identified from the merged transcriptomes, and 861 of them were differentially expressed. Relative to their expression on barley, all differentially expressed P450 (cytochrome 450 monooxygenase) genes (e.g., SavCYP6A13, SavCYP4C1, and SavCYP4G15) in S. avenae on wheat were upregulated, suggesting the significance of P450s in xenobiotic metabolism for this aphid on wheat. For S. avenae on barley, some genes encoding for ATP-binding cassette (ABC) transporters (e.g., ABCG1, ABCG4, ABCB7, and ABCA5) and UDP-glucuronosyltransferases (UGTs; e.g., UGT2B1 and UGT2C1) showed a dramatic increase in expression, suggesting that ABC transporters and UGTs could be critical for detoxification metabolism in S. avenae on barley. In addition, the expression for nearly all differentially expressed proteases was upregulated for S. avenae on barley, reflecting xenobiotic challenges facing S. avenae on resistant barley. Thus, various detoxification and other defense-related genes (e.g., proteases and oxidases) showed highly inducible transcript regulation, suggesting great adaptive potential for S. avenae on different plants. All the abovementioned genes will be prime candidates for further studies of molecular interactions underlying host plant use and specialization in this aphid. Our results provide insights into aphids' defenses against variable phytochemicals, and the molecular basis underlying aphid-plant interactions as well.

Draft genome of the cotton aphid Aphis gossypii

The cotton aphid Aphis gossypii Glover is a worldwide agricultural pest that feeds on cotton, melon, and other landscape plants, causing a high level of economic loss. In addition to the common characteristics shared with other aphids, Ap. gossypii has evolved multiple biotypes that present substantial differences in host adaption. These intriguing biological features are of interest from both a fundamental and applied perspective. However, the molecular studies of Ap. gossypii have been restrained by the lack of a reference genome. Furthermore, in order to establish a platform for the development of novel and sustainable control methods, it is necessary to generate genomic resources for Ap. gossypii. Here, we present a 294 Mb draft genome sequence of Ap. gossypii, which consists of 4,724 scaffolds with an N50 size of 438 kb. Compared to other aphid species with published genomes, Ap. gossypii presents the most compact genome size. A total of 14,694 protein-coding genes were predicted and annotated in the consensus gene set, 98.03% of CEGMA genes and 93.5% of BUSCO genes were captured respectively. Genome-wide selection analyses revealed that significantly evolving pathways in the genus Aphis are related to biological processes of detoxification, steroid biosynthesis, and ethylbenzene degradation. The acquisition of the genome of Ap. gossypii makes it possible to understand the molecular mechanism of intricate biological traits of this species, and will further facilitate the study of aphid evolution.

Deep Transcriptomic Analysis of Black Rockfish (Sebastes schlegelii) Provides New Insights on Responses to Acute Temperature Stress

In the present study, we conducted an RNA-Seq analysis to characterize the genes and pathways involved in acute thermal and cold stress responses in the liver of black rockfish, a viviparous teleost that has the ability to cope with a wide range of temperature changes. A total of 584 annotated differentially expressed genes (DEGs) were identified in all three comparisons (HT vs NT, HT vs LT and LT vs NT). Based on an enrichment analysis, DEGs with a potential role in stress accommodation were classified into several categories, including protein folding, metabolism, immune response, signal transduction, molecule transport, membrane, and cell proliferation/apoptosis. Considering that thermal stress has a greater effect than cold stress in black rockfish, 24 shared DEGs in the intersection of the HT vs LT and HT vs NT groups were enriched in 2 oxidation-related gene ontology (GO) terms. Nine important heat-stress-reducing pathways were significantly identified and classified into 3 classes: immune and infectious diseases, organismal immune system and endocrine system. Eight DEGs (early growth response protein 1, bile salt export pump, abcb11, hsp70a, rtp3, 1,25-dihydroxyvitamin d(3) 24-hydroxylase, apoa4, transcription factor jun-b-like and an uncharacterized gene) were observed among all three comparisons, strongly implying their potentially important roles in temperature stress responses.