Deep sequencing of the transcriptomes of soybean aphid and associated endosymbionts

Host-specific co-evolution likely driven by diet in Buchnera aphidicola

BACKGROUND

Russian wheat aphid (Diuraphis noxia Kurd.) is a severe pest to wheat, and even though resistance varieties are available to curb this pest, they are becoming obsolete with the development of new virulent aphid populations. Unlike many other aphids, D noxia only harbours a single endosymbiont, Buchnera aphidicola. Considering the importance of Buchnera, this study aimed to elucidate commonalities and dissimilarities between various hosts, to better understand its distinctiveness within its symbiotic relationship with D. noxia. To do so, the genome of the D. noxia's Buchnera was assembled and compared to those of other aphid species that feed on diverse host species.

RESULTS

The overall importance of several features such as gene length and percentage GC content was found to be critical for the maintenance of Buchnera genes when compared to their closest free-living relative, Escherichia coli. Buchnera protein coding genes were found to have percentage GC contents that tended towards a mean of ~ 26% which had strong correlation to their identity to their E. coli homologs. Several SNPs were identified between different aphid populations and multiple isolates of Buchnera were confirmed in single aphids.

CONCLUSIONS

Establishing the strong correlation of percentage GC content of protein coding genes and gene identity will allow for identifying which genes will be lost in the continually shrinking Buchnera genome. This is also the first report of a parthenogenically reproducing aphid that hosts multiple Buchnera strains in a single aphid, raising questions regarding the benefits of maintaining multiple strains. We also found preliminary evidence for post-transcriptional regulation of Buchnera genes in the form of polyadenylation.

Patterns of asexual reproduction of the soybean aphid, Aphis glycines (Matsumura), with and without the secondary symbionts Wolbachia and Arsenophonus, on susceptible and resistant soybean genotypes

Plant breeding is used to develop crops with host resistance to aphids, however, virulent biotypes often develop that overcome host resistance genes. We tested whether the symbionts, Arsenophonus (A) and Wolbachia (W), affect virulence and fecundity in soybean aphid biotypes Bt1 and Bt3 cultured on whole plants and detached leaves of three resistant, Rag1, Rag2 and Rag1 + 2, and one susceptible, W82, soybean genotypes. Whole plants and individual aphid experiments of A. glycines with and without Arsenophonus and Wolbachia did not show differences in overall fecundity. Differences were observed in peak fecundity, first day of deposition, and day of maximum nymph deposition of individual aphids on detached leaves. Bt3 had higher fecundity than Bt1 on detached leaves of all plant genotypes regardless of bacterial profile. Symbionts did not affect peak fecundity of Bt1 but increased it in Bt3 (A+W+) and all Bt3 strains began to deposit nymphs earlier than the Bt1 (A+W-). Arsenophonus in Bt1 delayed the first day of nymph deposition in comparison to aposymbiotic Bt1 except when reared on Rag1 + 2. For the Bt1 and Bt3 strains, symbionts did not result in a significant difference in the day they deposited the maximum number of nymphs nor was there a difference in survival or variability in number of nymphs deposited. Variability of number of aphids deposited was higher in aphids feeding on resistant plant genotypes. The impact of Arsenophonus on soybean aphid patterns of fecundity was dependent on the aphid biotype and plant genotype. Wolbachia alone had no detectable impact but may have contributed to the increased fecundity of Bt3 (A+W+). An individual based model, using data from the detached leaves experiment and with intraspecific competition removed, found patterns similar to those observed in the greenhouse and growth chamber experiments including a significant interaction between soybean genotype and aphid strain. Combining individual data with the individual based model of population growth isolated the impact of fecundity and host resistance from intraspecific competition and host health. Changes to patterns of fecundity, influenced by the composition and concentration of symbionts, may contribute to competitive interactions among aphid genotypes and influence selection on virulent aphid populations.

Which Is Stronger? A Continuing Battle Between Cry Toxins and Insects

In this article, we review the latest works on the insecticidal mechanisms of Bacillus thuringiensis Cry toxins and the resistance mechanisms of insects against Cry toxins. Currently, there are two models of insecticidal mechanisms for Cry toxins, namely, the sequential binding model and the signaling pathway model. In the sequential binding model, Cry toxins are activated to bind to their cognate receptors in the mid-intestinal epithelial cell membrane, such as the glycophosphatidylinositol (GPI)-anchored aminopeptidases-N (APNs), alkaline phosphatases (ALPs), cadherins, and ABC transporters, to form pores that elicit cell lysis, while in the signaling pathway model, the activated Cry toxins first bind to the cadherin receptor, triggering an extensive cell signaling cascade to induce cell apoptosis. However, these two models cannot seem to fully describe the complexity of the insecticidal process of Cry toxins, and new models are required. Regarding the resistance mechanism against Cry toxins, the main method insects employed is to reduce the effective binding of Cry toxins to their cognate cell membrane receptors by gene mutations, or to reduce the expression levels of the corresponding receptors by trans-regulation. Moreover, the epigenetic mechanisms, host intestinal microbiota, and detoxification enzymes also play significant roles in the insects' resistance against Cry toxins. Today, high-throughput sequencing technologies like transcriptomics, proteomics, and metagenomics are powerful weapons for studying the insecticidal mechanisms of Cry toxins and the resistance mechanisms of insects. We believe that this review shall shed some light on the interactions between Cry toxins and insects, which can further facilitate the development and utilization of Cry toxins.

Microbiome diversity of cotton aphids (Aphis gossypii) is associated with host alternation

Aphids are infected by a series of bacteria that can help them survive on specific host plants. However, the associations between aphids and these bacteria are not clear, and the bacterial communities in many aphid species are poorly characterized. Here, we investigated the bacterial communities of cotton aphids (Aphis gossypii) on 2 representative winter host plants and transferred to 3 summer host plants by 16S rDNA sequencing using the Illumina MiSeq platform. Our results revealed that the bacterial communities varied among cotton aphids on hibiscus, cotton aphids on pomegranate, cotton aphids on cotton transferred from hibiscus, cotton aphids on muskmelon transferred from hibiscus, cotton aphids on cucumber transferred from hibiscus,. The diversity and richness of the bacterial communities were significantly higher in aphids on muskmelon and aphids on cucumber than in the other treatments. There were two main factors influencing the distribution of internal bacterial OTUs revealed by principal component analysis, including the differences among Punicaceae, Malvaceae and Cucurbitaceae. There were 28 bacterial communities with significant differences between two arbitrary treatments, which could be grouped into 6 main clusters depending on relative abundance. Moreover, our results indicated that in addition to the obligate endosymbiont Buchnera, with a dominant position (> 52%), A. gossypii also harbored 3 facultative endosymbiotic bacteria (Serratia, Arsenophonus, and Wolbachia) and 3 possibly symbiotic bacteria (Acinetobacter, Pantoea, and Flavobacterium). There were several correspondences between the symbiotic bacteria in cotton aphids and the specific host plants of the aphids. This study provides a better understanding of the interactions among symbiotic bacteria, aphids and host plants, suggesting that the selection pressure on aphid bacterial communities is likely to be exerted by the species of host plants.

Molecular cloning, characterization, and expression profiling analysis of Cry toxin receptor genes from sugarcane shoot borer Chilo infuscatellus (Snellen)

The sugarcane shoot borer Chilo infuscatellus (Snellen) is known for causing severe damage to sugarcane yield in China. Methods have been developed to control this pest, including Cry toxin pesticide and transgenic Bt plants. In order to investigate the molecular mechanism of the Cry toxin binding process and provide a basis for understanding the insect's resistance mechanism, we used a high throughput sequencing platform to perform a de novo transcriptome assembly across different larval developmental stages and analyzed Cry toxin receptors based on our assembled transcripts. We cloned twelve Cry toxin receptor genes including 1 cadherin (Cad), 7 aminopeptidase-Ns (APNs), 3 alkaline phosphatases (ALPs), and 1 ATP-binding cassette transporter subfamily C2 (ABCC2), and three of them with full length. The sublethal dosage of Cry1Ac toxin was applied to sugarcane shoot borer and identified some Cry toxin receptor genes that were significantly induced after 48 h of exposure. Furthermore, quantitative RT-PCR was conducted to detect the expression profiles of these genes. Our transcriptome sequence data provided a valuable molecular resource for further study and the identified Cry toxin receptor data gave insights for improved research into the mechanism of Bt resistance.

Variation and diversification of the microbiome of Schlechtendalia chinensis on two alternate host plants

Schlechtendalia chinensis, a gall-inducing aphid, has two host plants in its life cycle. Its wintering host is a moss (typically Plagiomnium maximoviczii) and its main host is Rhus chinensis (Sumac), on which it forms galls during the summer. This study investigated bacteria associated with S. chinensis living on the two different host plants by sequencing 16S rRNAs. A total of 183 Operational Taxonomic Units (OTUs) from 50 genera were identified from aphids living on moss, whereas 182 OTUs from 49 genera were found from aphids living in Sumac galls. The most abundant bacterial genus among identified OTUs from aphids feeding on both hosts was Buchnera. Despite similar numbers of OTUs, the composition of bacterial taxa showed significant differences between aphids living on moss and those living on R. chinensis. Specifically, there were 12 OTUs from 5 genera (family) unique to aphids living on moss, and 11 OTUs from 4 genera (family) unique to aphids feeding in galls on R. chinensis. Principal Coordinate Analysis (PCoA) also revealed that bacteria from moss-residing aphids clustered differently from aphids collected from galls. Our results provide a foundation for future analyses on the roles of symbiotic bacteria in plant-aphid interactions in general, and how gall-specific symbionts differ in this respect.

Sequence analysis of the potato aphid Macrosiphum euphorbiae transcriptome identified two new viruses

The potato aphid, Macrosiphum euphorbiae, is an important agricultural pest that causes economic losses to potato and tomato production. To establish the transcriptome for this aphid, RNA-Seq libraries constructed from aphids maintained on tomato plants were used in Illumina sequencing generating 52.6 million 75-105 bp paired-end reads. The reads were assembled using Velvet/Oases software with SEED preprocessing resulting in 22,137 contigs with an N50 value of 2,003bp. After removal of contigs from tomato host origin, 20,254 contigs were annotated using BLASTx searches against the non-redundant protein database from the National Center for Biotechnology Information (NCBI) as well as IntereProScan. This identified matches for 74% of the potato aphid contigs. The highest ranking hits for over 12,700 contigs were against the related pea aphid, Acyrthosiphon pisum. Gene Ontology (GO) was used to classify the identified M. euphorbiae contigs into biological process, cellular component and molecular function. Among the contigs, sequences of microbial origin were identified. Sixty five contigs were from the aphid bacterial obligate endosymbiont Buchnera aphidicola origin and two contigs had amino acid similarities to viruses. The latter two were named Macrosiphum euphorbiae virus 2 (MeV-2) and Macrosiphum euphorbiae virus 3 (MeV-3). The highest sequence identity to MeV-2 had the Dysaphis plantaginea densovirus, while to MeV-3 is the Hubei sobemo-like virus 49. Characterization of MeV-2 and MeV-3 indicated that both are transmitted vertically from adult aphids to nymphs. MeV-2 peptides were detected in the aphid saliva and only MeV-2 and not MeV-3 nucleic acids were detected inside tomato leaves exposed to virus-infected aphids. However, MeV-2 nucleic acids did not persist in tomato leaf tissues, after clearing the plants from aphids, indicating that MeV-2 is likely an aphid virus.

Analysis of Homologs of Cry-toxin Receptor-Related Proteins in the Midgut of a Non-Bt Target, Nilaparvata lugens (Stål) (Hemiptera: Delphacidae)

The brown planthopper (BPH) Nilaparvata lugens is one of the most destructive insect pests in the rice fields of Asia. Like other hemipteran insects, BPH is not susceptible to Cry toxins of Bacillus thuringiensis (Bt) or transgenic rice carrying Bt cry genes. Lack of Cry receptors in the midgut is one of the main reasons that BPH is not susceptible to the Cry toxins. The main Cry-binding proteins (CBPs) of the susceptible insects are cadherin, aminopeptidase N (APN), and alkaline phosphatase (ALP). In this study, we analyzed and validated de novo assembled transcripts from transcriptome sequencing data of BPH to identify and characterize homologs of cadherin, APN, and ALP. We then compared the cadherin-, APN-, and ALP-like proteins of BPH to previously reported CBPs to identify their homologs in BPH. The sequence analysis revealed that at least one cadherin, one APN, and two ALPs of BPH contained homologous functional domains identified from the Cry-binding cadherin, APN, and ALP, respectively. Quantitative real-time polymerase chain reaction used to verify the expression level of each putative Cry receptor homolog in the BPH midgut indicated that the CBPs homologous APN and ALP were expressed at high or medium-high levels while the cadherin was expressed at a low level. These results suggest that homologs of CBPs exist in the midgut of BPH. However, differences in key motifs of CBPs, which are functional in interacting with Cry toxins, may be responsible for insusceptibility of BPH to Cry toxins.

Host Plants Indirectly Influence Plant Virus Transmission by Altering Gut Cysteine Protease Activity of Aphid Vectors

The green peach aphid, Myzus persicae, is a vector of the Potato leafroll virus (PLRV, Luteoviridae), transmitted exclusively by aphids in a circulative manner. PLRV transmission efficiency was significantly reduced when a clonal lineage of M. persicae was reared on turnip as compared with the weed physalis, and this was a transient effect caused by a host-switch response. A trend of higher PLRV titer in physalis-reared aphids as compared with turnip-reared aphids was observed at 24 h and 72 h after virus acquisition. The major difference in the proteomes of these aphids was the up-regulation of predicted lysosomal enzymes, in particular the cysteine protease cathepsin B (cathB), in aphids reared on turnip. The aphid midgut is the site of PLRV acquisition, and cathB and PLRV localization were starkly different in midguts of the aphids reared on the two host plants. In viruliferous aphids that were reared on turnip, there was near complete colocalization of cathB and PLRV at the cell membranes, which was not observed in physalis-reared aphids. Chemical inhibition of cathB restored the ability of aphids reared on turnip to transmit PLRV in a dose-dependent manner, showing that the increased activity of cathB and other cysteine proteases at the cell membrane indirectly decreased virus transmission by aphids. Understanding how the host plant influences virus transmission by aphids is critical for growers to manage the spread of virus among field crops.

Aphis Glycines Virus 2, a Novel Insect Virus with a Unique Genome Structure

The invasive soybean aphid, Aphis glycines, is a major pest in soybeans, resulting in substantial economic loss. We analyzed the A. glycines transcriptome to identify sequences derived from viruses of A. glycines. We identified sequences derived from a novel virus named Aphis glycines virus 2 (ApGlV2). The assembled virus genome sequence was confirmed by reverse transcription polymerase chain reaction (RT-PCR) and Sanger sequencing, conserved domains were characterized, and distribution, and transmission examined. This virus has a positive sense, single-stranded RNA genome of ~4850 nt that encodes three proteins. The RNA-dependent RNA polymerase (RdRp) of ApGlV2 is a permuted RdRp similar to those of some tetraviruses, while the capsid protein is structurally similar to the capsid proteins of plant sobemoviruses. ApGlV2 also encodes a larger minor capsid protein, which is translated by a readthrough mechanism. ApGlV2 appears to be widespread in A. glycines populations and to persistently infect aphids with a 100% vertical transmission rate. ApGlV2 is susceptible to the antiviral RNA interference (RNAi) pathway. This virus, with its unique genome structure with both plant- and insect-virus characteristics, is of particular interest from an evolutionary standpoint.

Heritable symbionts in a world of varying temperature

Heritable microbes represent an important component of the biology, ecology and evolution of many plants, animals and fungi, acting as both parasites and partners. In this review, we examine how heritable symbiont-host interactions may alter host thermal tolerance, and how the dynamics of these interactions may more generally be altered by thermal environment. Obligate symbionts, those required by their host, are considered to represent a thermally sensitive weak point for their host, associated with accumulation of deleterious mutations. As such, these symbionts may represent an important determinant of host thermal envelope and spatial distribution. We then examine the varied relationship between thermal environment and the frequency of facultative symbionts that provide ecologically contingent benefits or act as parasites. We note that some facultative symbionts directly alter host thermotolerance. We outline how thermal environment will alter the benefits/costs of infection more widely, and additionally modulate vertical transmission efficiency. Multiple patterns are observed, with symbionts being cold sensitive in some species and heat sensitive in others, with varying and non-coincident thresholds at which phenotype and transmission are ablated. Nevertheless, it is clear that studies aiming to predict ecological and evolutionary dynamics of symbiont-host interactions need to examine the interaction across a range of thermal environments. Finally, we discuss the importance of thermal sensitivity in predicting the success/failure of symbionts to spread into novel species following natural/engineered introduction.

Transcriptome analysis provides insight into venom evolution in a seed-parasitic wasp, Megastigmus spermotrophus

One of the most striking host range transitions is the evolution of plant parasitism from animal parasitism. Parasitoid wasps that have secondarily evolved to attack plants (ie gall wasps and seed-feeders) demonstrate intimate associations with their hosts, yet the mechanism of plant-host manipulation is currently not known. There is, however, emerging evidence suggesting that ovipositional secretions play a role in plant manipulation. To investigate whether parasites have modified pre-existing adaptations to facilitate dramatic host shifts we aimed to characterize the expression of venom proteins in a plant parasite using a collection of parasitoid venom sequences as a guide. The transcriptome of a seed-feeding wasp, Megastigmus spermotrophus, was assembled de novo and three putative venoms were found to be highly expressed in adult females. One of these putative venoms, aspartylglucosaminidase, has been previously identified as a major venom component in two distantly related parasitoid wasps (Asobara tabida and Leptopilina heterotoma) and may have originated via gene duplication within the Hymenoptera. Our study shows that M. spermotrophus, a specialized plant parasite, expresses putative venom transcripts that share homology to venoms identified in Nasonia vitripennis (both superfamily Chalcidoidea), which suggests that M. spermotrophus may have co-opted pre-existing machinery to develop as a plant parasite.

Bacterial communities of the cotton aphid Aphis gossypii associated with Bt cotton in northern China

Aphids are infected with a wide variety of endosymbionts that can confer ecologically relevant traits. However, the bacterial communities of most aphid species are still poorly characterized. This study investigated the bacterial diversity of the cotton aphid Aphis gossypii associated with Bt cotton in northern China by targeting the V4 region of the 16S rDNA using the Illumina MiSeq platform. Our sequencing data revealed that bacterial communities of A. gossypii were generally dominated by the primary symbiont Buchnera, together with the facultative symbionts Arsenophonus and Hamiltonella. To our knowledge, this is the first report documenting the facultative symbiont Hamiltonella in A. gossypii. Moreover, the bacterial community structure was similar within aphids from the same province, but distinct among those from different provinces. The taxonomic diversity of the bacterial community is greater in Hebei Province compared with in samples from Henan and Shandong Provinces. The selection pressure exerted by the different geographical locations could explain the differences found among the various provinces. These findings broaden our understanding of the interactions among aphids, endosymbionts and their environments, and provide clues to develop potential biocontrol techniques against this cotton aphid.

Comparative transcriptomics and proteomics of three different aphid species identifies core and diverse effector sets

BACKGROUND

Aphids are phloem-feeding insects that cause significant economic losses to agriculture worldwide. While feeding and probing these insects deliver molecules, called effectors, inside their host to enable infestation. The identification and characterization of these effectors from different species that vary in their host range is an important step in understanding the infestation success of aphids and aphid host range variation. This study employs a multi-disciplinary approach based on transcriptome sequencing and proteomics to identify and compare effector candidates from the broad host range aphid Myzus persicae (green peach aphid) (genotypes O, J and F), and narrow host range aphids Myzus cerasi (black cherry aphid) and Rhopalosiphum padi (bird-cherry oat aphid).

RESULTS

Using a combination of aphid transcriptome sequencing on libraries derived from head versus body tissues as well as saliva proteomics we were able to predict candidate effectors repertoires from the different aphid species and genotypes. Among the identified conserved or core effector sets, we identified a significant number of previously identified aphid candidate effectors indicating these proteins may be involved in general infestation strategies. Moreover, we identified aphid candidate effector sequences that were specific to one species, which are interesting candidates for further validation and characterization with regards to species-specific functions during infestation. We assessed our candidate effector repertoires for evidence of positive selection, and identified 49 candidates with DN/DS ratios >1. We noted higher rates of DN/DS ratios in predicted aphid effectors than non-effectors. Whether this reflects positive selection due to co-evolution with host plants, or increased neofunctionalization upon gene duplication remains to be investigated.

CONCLUSION

Our work provides a comprehensive overview of the candidate effector repertoires from three different aphid species with varying host ranges. Comparative analyses revealed candidate effectors that are most likely are involved in general aspects of infestation, whereas others, that are highly divergent, may be involved in specific processes important for certain aphid species. Insights into the overlap and differences in aphid effector repertoires are important in understanding how different species successfully infest different ranges of plant species.

Loop replacements with gut-binding peptides in Cry1Ab domain II enhanced toxicity against the brown planthopper, Nilaparvata lugens (Stål)

Bacillus thuringiensis (Bt) Cry toxins have been used widely in pest managements. However, Cry toxins are not effective against sap-sucking insects (Hemiptera), which limits the application of Bt for pest management. In order to extend the insecticidal spectrum of Bt toxins to the rice brown planthopper (BPH), Nilaparvata lugens, we modified Cry1Ab putative receptor binding domains with selected BPH gut-binding peptides (GBPs). Three surface exposed loops in the domain II of Cry1Ab were replaced with two GBPs (P2S and P1Z) respectively. Bioassay results showed that toxicity of modified toxin L2-P2S increased significantly (~9 folds) against BPH nymphs. In addition, damage of midgut cells was observed from the nymphs fed with L2-P2S. Our results indicate that modifying Cry toxins based on the toxin-gut interactions can broaden the insecticidal spectrum of Bt toxin. This method provides another approach for the development of transgenic crops with novel insecticidal activity against hemipteran insects and insect populations resistant to current Bt transgenic crops.

Stress-induced changes in abundance differ among obligate and facultative endosymbionts of the soybean aphid

Bacterial endosymbionts can drive evolutionary novelty by conferring adaptive benefits under adverse environmental conditions. Among aphid species there is growing evidence that symbionts influence tolerance to various forms of stress. However, the extent to which stress inflicted on the aphid host has cascading effects on symbiont community dynamics remains poorly understood. Here we simultaneously quantified the effect of host-plant induced and xenobiotic stress on soybean aphid (Aphis glycines) fitness and relative abundance of its three bacterial symbionts. Exposure to soybean defensive stress (Rag1 gene) and a neurotoxic insecticide (thiamethoxam) substantially reduced aphid composite fitness (survival × reproduction) by 74 ± 10% and 92 ± 2%, respectively, which in turn induced distinctive changes in the endosymbiont microbiota. When challenged by host-plant defenses a 1.4-fold reduction in abundance of the obligate symbiont Buchnera was observed across four aphid clonal lines. Among facultative symbionts of Rag1-stressed aphids, Wolbachia abundance increased twofold and Arsenophonus decreased 1.5-fold. A similar pattern was observed under xenobiotic stress, with Buchnera and Arsenophonus titers decreasing (1.3-fold) and Wolbachia increasing (1.5-fold). Furthermore, variation in aphid virulence to Rag1 was positively correlated with changes in Arsenophonus titers, but not Wolbachia or Buchnera. A single Arsenophonus multi-locus genotype was found among aphid clonal lines, indicating strain diversity is not primarily responsible for correlated host-symbiont stress levels. Overall, our results demonstrate the nature of aphid symbioses can significantly affect the outcome of interactions under stress and suggests general changes in the microbiome can occur across multiple stress types.

Whole Genome Sequence of the Soybean Aphid Endosymbiont Buchnera aphidicola and Genetic Differentiation among Biotype-Specific Strains

Endosymbiosis with microorganisms is common in insects, with more than 10% of species requiring the metabolic capabilities of intracellular bacteria for their nutrient acquisition. Aphids harbor an obligate mutualism with the vertically transferred endosymbiont, Buchnera aphidicola, which produces key nutrients lacking in the aphid's phloem-based diet that are necessary for normal development and reproduction. It is thought that, in some groups of insects, bacterial symbionts may play key roles in biotype evolution against host-plant resistance. The genome of Buchnera has been sequenced in several aphid strains but little genomic data is currently available for the soybean aphid (Aphis glycines), one of the most important pests of soybean in North America. In this study, DNA sequencing was used to assemble and annotate the genome sequence of the Buchnera A. glycines strain and to reconstruct phylogenetic relationships among different strains. In addition, we identified several fixed Buchnera SNPs between Aphis glycines biotypes that were avirulent or virulent to a soybean aphid resistance gene (Rag1). The results of this study describe the genetic and evolutionary relationships of the Buchnera A. glycines strain, and begin to define the roles of an aphid symbiont in host-plant resistance.

Shifts in Buchnera aphidicola density in soybean aphids (Aphis glycines) feeding on virus-infected soybean

Vertically transmitted bacterial symbionts are common in arthropods. Aphids undergo an obligate symbiosis with Buchnera aphidicola, which provides essential amino acids to its host and contributes directly to nymph growth and reproduction. We previously found that newly adult Aphis glycines feeding on soybean infected with the beetle-transmitted Bean pod mottle virus (BPMV) had significantly reduced fecundity. We hypothesized that the reduced fecundity was attributable to detrimental impacts of the virus on the aphid microbiome, namely Buchnera. To test this, mRNA sequencing and quantitative real-time PCR were used to assay Buchnera transcript abundance and titre in A. glycines feeding on Soybean mosaic virus-infected, BPMV-infected, and healthy soybean for up to 14 days. Our results indicated that Buchnera density was lower and ultimately suppressed in aphids feeding on virus-infected soybean. While the decreased Buchnera titre may be associated with reduced aphid fecundity, additional mechanisms are probably involved. The present report begins to describe how interactions among insects, plants, and plant pathogens influence endosymbiont population dynamics.

Whitefly genome expression reveals host-symbiont interaction in amino acid biosynthesis

Background

Whitefly (Bemisia tabaci) complex is a serious insect pest of several crop plants worldwide. It comprises several morphologically indistinguishable species, however very little is known about their genetic divergence and biosynthetic pathways. In the present study, we performed transcriptome sequencing of Asia 1 species of B. tabaci complex and analyzed the interaction of host-symbiont genes in amino acid biosynthetic pathways.

Methodology/Principal Findings

We obtained about 83 million reads using Illumina sequencing that assembled into 72716 unitigs. A total of 21129 unitigs were annotated at stringent parameters. Annotated unitigs were mapped to 52847 gene ontology (GO) terms and 131 Kyoto encyclopedia of genes and genomes (KEGG) pathways. Expression analysis of the genes involved in amino acid biosynthesis pathways revealed the complementation between whitefly and its symbiont partner Candidatus Portiera aleyrodidarum. Most of the non-essential amino acids and intermediates of essential amino acid pathways were supplied by the host insect to its symbiont. The symbiont expressed the pathways for the essential amino acids arginine, threonine and tryptophan and the immediate precursors of valine, leucine, isoleucine and phenyl-alanine. High level expression of the amino acid transporters in the whitefly suggested the molecular mechanisms for the exchange of amino acids between the host and the symbiont.

Conclusions/Significance

Our study provides a comprehensive transcriptome data for Asia 1 species of B. tabaci complex that focusses light on integration of host and symbiont genes in amino acid biosynthesis pathways.

RNA virus discovery in insects

The advent of next generation sequencing (NGS) technology has allowed for significant advances in insect RNA virus discovery, particularly for identification of covert viruses and for identification of previously undescribed virus groups. Recent research has revealed the complexity and diversity both of viral populations within an insect, and of different virus groups infecting insects. While appropriate verification of the ability of a given virus to infect a putative host is essential, this is frequently challenging for newly discovered viruses. In this chapter, we describe recently discovered RNA viruses of insects, optimization of NGS data analysis for virus discovery, and challenges associated with virus nomenclature.

Emerging strategies for RNA interference (RNAi) applications in insects

RNA interference (RNAi) in insects is a gene regulatory process that also plays a vital role in the maintenance and in the regulation of host defenses against invading viruses. Small RNAs determine the specificity of the RNAi through precise recognition of their targets. These small RNAs in insects comprise small interfering RNAs (siRNAs), micro RNAs (miRNAs) and Piwi interacting RNAs (piRNAs) of various lengths. In this review, we have explored different forms of the RNAi inducers that are presently in use, and their applications for an effective and efficient fundamental and practical RNAi research with insects. Further, we reviewed trends in next generation sequencing (NGS) technologies and their importance for insect RNAi, including the identification of novel insect targets as well as insect viruses. Here we also describe a rapidly emerging trend of using plant viruses to deliver the RNAi inducer molecules into insects for an efficient RNAi response.

Abiotic and biotic stressors causing equivalent mortality induce highly variable transcriptional responses in the soybean aphid

Environmental stress affects basic organismal functioning and can cause physiological, developmental, and reproductive impairment. However, in many nonmodel organisms, the core molecular stress response remains poorly characterized and the extent to which stress-induced transcriptional changes differ across qualitatively different stress types is largely unexplored. The current study examines the molecular stress response of the soybean aphid (Aphis glycines) using RNA sequencing and compares transcriptional responses to multiple stressors (heat, starvation, and plant defenses) at a standardized stress level (27% adult mortality). Stress-induced transcriptional changes showed remarkable variation, with starvation, heat, and plant defensive stress altering the expression of 3985, 510, and 12 genes, respectively. Molecular responses showed little overlap across all three stressors. However, a common transcriptional stress response was identified under heat and starvation, involved with up-regulation of glycogen biosynthesis and molecular chaperones and down-regulation of bacterial endosymbiont cellular and insect cuticular components. Stressor-specific responses indicated heat affected expression of heat shock proteins and cuticular components, whereas starvation altered a diverse set of genes involved in primary metabolism, oxidative reductive processes, nucleosome and histone assembly, and the regulation of DNA repair and replication. Exposure to host plant defenses elicited the weakest response, of which half of the genes were of unknown function. This study highlights the need for standardizing stress levels when comparing across stress types and provides a basis for understanding the role of general vs. stressor specific molecular responses in aphids.

Circulative, "nonpropagative" virus transmission: an orchestra of virus-, insect-, and plant-derived instruments

Species of plant viruses within the Luteoviridae, Geminiviridae, and Nanoviridae are transmitted by phloem-feeding insects in a circulative, nonpropagative manner. The precise route of virus movement through the vector can differ across and within virus families, but these viruses all share many biological, biochemical, and ecological features. All share temporal and spatial constraints with respect to transmission efficiency. The viruses also induce physiological changes in their plant hosts resulting in behavioral changes in the insects that optimize the transmission of virus to new hosts. Virus proteins interact with insect, endosymbiont, and plant proteins to orchestrate, directly and indirectly, virus movement in insects and plants to facilitate transmission. Knowledge of these complex interactions allows for the development of new tools to reduce or prevent transmission, to quickly identify important vector populations, and to improve the management of these economically important viruses affecting agricultural and natural plant populations.

De novo transcriptome assembly for a non-model species, the blood-sucking bug Triatoma brasiliensis, a vector of Chagas disease

High throughput sequencing (HTS) provides new research opportunities for work on non-model organisms, such as differential expression studies between populations exposed to different environmental conditions. However, such transcriptomic studies first require the production of a reference assembly. The choice of sampling procedure, sequencing strategy and assembly workflow is crucial. To develop a reliable reference transcriptome for Triatoma brasiliensis, the major Chagas disease vector in Northeastern Brazil, different de novo assembly protocols were generated using various datasets and software. Both 454 and Illumina sequencing technologies were applied on RNA extracted from antennae and mouthparts from single or pooled individuals. The 454 library yielded 278 Mb. Fifteen Illumina libraries were constructed and yielded nearly 360 million RNA-seq single reads and 46 million RNA-seq paired-end reads for nearly 45 Gb. For the 454 reads, we used three assemblers, Newbler, CAP3 and/or MIRA and for the Illumina reads, the Trinity assembler. Ten assembly workflows were compared using these programs separately or in combination. To compare the assemblies obtained, quantitative and qualitative criteria were used, including contig length, N50, contig number and the percentage of chimeric contigs. Completeness of the assemblies was estimated using the CEGMA pipeline. The best assembly (57,657 contigs, completeness of 80 %, <1 % chimeric contigs) was a hybrid assembly leading to recommend the use of (1) a single individual with large representation of biological tissues, (2) merging both long reads and short paired-end Illumina reads, (3) several assemblers in order to combine the specific advantages of each.

Analysis of new aphid lethal paralysis virus (ALPV) isolates suggests evolution of two ALPV species

Aphid lethal paralysis virus (ALPV; family Dicistroviridae) was first isolated from the bird cherry-oat aphid, Rhopalosiphum padi. ALPV-like virus sequences have been reported from many insects and insect predators. We identified a new isolate of ALPV (ALPV-AP) from the pea aphid, Acyrthosiphon pisum, and a new isolate (ALPV-DvV) from western corn rootworm, Diabrotica virgifera virgifera. ALPV-AP has an ssRNA genome of 9940 nt. Based on phylogenetic analysis, ALPV-AP was closely related to ALPV-AM, an ALPV isolate from honeybees, Apis mellifera, in Spain and Brookings, SD, USA. The distinct evolutionary branches suggested the existence of two lineages of the ALPV virus. One consisted of ALPV-AP and ALPV-AM, whilst all other isolates of ALPV grouped into the other lineage. The similarity of ALPV-AP and ALPV-AM was up to 88 % at the RNA level, compared with 78-79 % between ALPV-AP and other ALPV isolates. The sequence identity of proteins between ALPV-AP and ALPV-AM was 98-99 % for both ORF1 and ORF2, whilst only 85-87 % for ORF1 and 91-92 % for ORF2 between ALPV-AP and other ALPV isolates. Sequencing of RACE (rapid amplification of cDNA ends) products and cDNA clones of the virus genome revealed sequence variation in the 5' UTRs and in ORF1, indicating that ALPV may be under strong selection pressure, which could have important biological implications for ALPV host range and infectivity. Our results indicated that ALPV-like viruses infect insects in the order Coleoptera, in addition to the orders Hemiptera and Hymenoptera, and we propose that ALPV isolates be classified as two separate viral species.

Variations in the identity and complexity of endosymbiont combinations in whitefly hosts

The target of natural selection is suggested to be the holobiont - the organism together with its associated symbiotic microorganisms. The well-defined endosymbiotic communities of insects make them a useful model for exploring the role of symbiotic interactions in shaping the functional repertoire of plants and animals. Here, we studied the variations in the symbiotic communities of the sweet potato whitefly Bemisia tabaci (Hemiptera: Aleyrodidae) by compiling a dataset of over 2000 individuals derived from several independent screenings. The secondary endosymbionts harbored by each individual were clustered into entities termed Facultative Endosymbiont Combinations (FECs), each representing a natural assemblage of co-occurring bacterial genera. The association of FECs with whitefly individuals stratified the otherwise homogeneous population into holobiont units. We both identified bacterial assemblages that are specific to whitefly groups sharing unique genetic backgrounds, and characterized the FEC variations within these groups. The analysis revealed that FEC complexity is positively correlated with both distance from the equator and specificity of the genetic clade of the host insect. These findings highlight the importance of symbiotic combinations in shaping the distribution patterns of B. tabaci and possibly other insect species.

Transcriptome analysis of "Candidatus Liberibacter solanacearum" in its psyllid vector, Bactericera cockerelli

"Candidatus Liberibacter solanacearum" (Lso) is an emergent pathogen of carrots in Europe and solanaceous plants in North and Central America and New Zealand. This bacterium is closely related to other pathogenic Candidatus Liberibacter spp., all vectored by psyllids. In order to understand the molecular interaction of this pathogen and its psyllid vector, Bactericera cockerelli, Illumina sequencing of psyllid harboring Lso was performed to determine if this approach could be used to assess the bacterial transcriptome in this association. Prior to sequencing, psyllid RNA was purified and insect and bacterial rRNA were removed. Mapping of reads to Lso genome revealed that over 92% of the bacterial genes were expressed in the vector, and that the COG categories Translation and Post-translational modification, protein turnover, chaperone functions were the most expressed functional categories. Expression levels of selected Lso genes were confirmed by RT-qPCR. The transcriptomic analysis also helped correct Lso genome annotation by identifying the expression of genes that were not predicted in the genome sequencing effort.

Reduction in fecundity and shifts in cellular processes by a native virus on an invasive insect

Pathogens and their vectors have coevolutionary histories that are intricately intertwined with their ecologies, environments, and genetic interactions. The soybean aphid, Aphis glycines, is native to East Asia but has quickly become one of the most important aphid pests in soybean-growing regions of North America. In this study, we used bioassays to examine the effects of feeding on soybean infected with a virus it vectors (Soybean mosaic virus [SMV]) and a virus it does not vector (Bean pod mottle virus [BPMV]) have on A. glycines survival and fecundity. The genetic underpinnings of the observed changes in fitness phenotype were explored using RNA-Seq. Aphids fed on SMV-infected soybean had transcriptome and fitness profiles that were similar to that of aphids fed on healthy control plants. Strikingly, a significant reduction in fecundity was seen in aphids fed on BPMV-infected soybean, concurrent with a large and persistent downregulation of A. glycines transcripts involved in regular cellular activities. Although molecular signatures suggested a small regulatory RNA pathway defense response was repressed in aphids feeding on infected plants, BPMV did not appear to be replicating in the vector. These results suggest that incompatibilities with BPMV or the effects of BPMV infection on soybean caused A. glycines to allot available energy resources to survival rather than reproduction and other core cellular processes. Ultimately, the detrimental impacts to A. glycines may reflect the short tritrophic evolutionary histories between the insect, plant, and virus.

Microbiome diversity of Aphis glycines with extensive superinfection in native and invasive populations

Associations among insects and microbes can lead to beneficial or parasitic interactions. Using 454 sequencing of 16S RNA genes, we compared microbiome diversity and abundance among field-collected (F) and laboratory-reared (L) populations of the invasive soybean aphid (Aphis glycines), a pest of soybean. Additionally, we screened A. glycines populations from native (Japan, South Korea and China) and invasive regions (North America) to broadly determine the microbiome diversity. Our results suggested that Arsenophonus (relative abundance of 54.6%), Buchnera (38.7%) and Wolbachia (3.7%) were the major bacteria associated with A. glycines. Arsenophonus was the most abundant in F populations but was significantly reduced in L populations; additional bacteria species also had lower relative abundances in L populations. Native and invasive populations were largely similar in bacteria communities and revealed substantial superinfection of Arsenophonus and Wolbachia. The lone exception was a lack of Arsenophonus in A. glycines from Japan. Divergent selection pressures among natural and laboratory populations were inferred as factors driving the differential bacterial communities observed. Our results will allow for improved comparative aphid-symbiont research and broaden our understanding of the interactions among insects, endosymbionts and their environments.

Sequencing and de novo assembly of the transcriptome of the glassy-winged sharpshooter (Homalodisca vitripennis)

BACKGROUND

The glassy-winged sharpshooter Homalodisca vitripennis (Hemiptera: Cicadellidae), is a xylem-feeding leafhopper and important vector of the bacterium Xylella fastidiosa; the causal agent of Pierce's disease of grapevines. The functional complexity of the transcriptome of H. vitripennis has not been elucidated thus far. It is a necessary blueprint for an understanding of the development of H. vitripennis and for designing efficient biorational control strategies including those based on RNA interference.

RESULTS

Here we elucidate and explore the transcriptome of adult H. vitripennis using high-throughput paired end deep sequencing and de novo assembly. A total of 32,803,656 paired-end reads were obtained with an average transcript length of 624 nucleotides. We assembled 32.9 Mb of the transcriptome of H. vitripennis that spanned across 47,265 loci and 52,708 transcripts. Comparison of our non-redundant database showed that 45% of the deduced proteins of H. vitripennis exhibit identity (e-value ≤1(-5)) with known proteins. We assigned Gene Ontology (GO) terms, Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations, and potential Pfam domains to each transcript isoform. In order to gain insight into the molecular basis of key regulatory genes of H. vitripennis, we characterized predicted proteins involved in the metabolism of juvenile hormone, and biogenesis of small RNAs (Dicer and Piwi sequences) from the transcriptomic sequences. Analysis of transposable element sequences of H. vitripennis indicated that the genome is less expanded in comparison to many other insects with approximately 1% of the transcriptome carrying transposable elements.

CONCLUSIONS

Our data significantly enhance the molecular resources available for future study and control of this economically important hemipteran. This transcriptional information not only provides a more nuanced understanding of the underlying biological and physiological mechanisms that govern H. vitripennis, but may also lead to the identification of novel targets for biorationally designed control strategies.

Plant-aphid interactions with a focus on legumes

Sap-sucking insects such as aphids cause substantial yield losses in agriculture by draining plant nutrients as well as vectoring viruses. The main method of control in agriculture is through the application of insecticides. However, aphids rapidly evolve mechanisms to detoxify these, so there is a need to develop durable plant resistance to these damaging insect pests. The focus of this review is on aphid interactions with legumes, but work on aphid interactions with other plants, particularly Arabidopsis and tomato is also discussed. This review covers advances on the plant side of the interaction, including the identification of major resistance genes and quantitative trait loci conferring aphid resistance in legumes, basal and resistance gene mediated defence signalling following aphid infestation and the role of specialised metabolites. On the aphid side of the interaction, this review covers what is known about aphid effector proteins and aphid detoxification enzymes. Recent advances in these areas have provided insight into mechanisms underlying resistance to aphids and the strategies used by aphids for successful infestations and have significant impacts for the delivery of durable resistance to aphids in legume crops.

Implementing an evolutionary framework for understanding genetic relationships of phenotypically defined insect biotypes in the invasive soybean aphid (Aphis glycines)

Adaptive evolution of pest insects in response to the introduction of resistant cultivars is well documented and commonly results in virulent (i.e., capable of feeding upon resistant cultivars) insect populations being labeled as distinct biotypes. Phenotypically defined, biotypes frequently remain evolutionarily indistinct, resulting in ineffective application of virulence control measures and shorter durability of resistant cultivars. Here, we utilize an evolutionary framework to discern the genetic relationship between biotypes of the soybean aphid (Aphis glycines, Matsumura). The soybean aphid is invasive in North America and is among the most destructive pests of commercial soybean on the continent. Attempts to breed host-plant-resistant soybean have been hampered by the emergence of virulent aphid biotypes that are unaffected by the plant's resistance mechanism(s). Comparative population genetic analysis of virulent and avirulent (i.e., unable to feed on resistant cultivars) biotypes found populations to be genetically indistinguishable across biotype and geographic distance, with high rates of interpopulation immigration and admixture. The lack of genetic distinction between biotypes coupled with elevated genotypic diversity within all populations suggested virulence has a nongenetic-based or includes a gene complex that is widely distributed throughout soybean aphid populations, which undergo regular dispersal and unimpeded sexual recombination.