Transcriptional response of bean bug (Riptortus pedestris) upon infection with entomopathogenic fungus, Beauveria bassiana JEF-007

Elucidating the ecophysiology of soybean pod-sucking stinkbug Riptortus pedestris (Hemiptera: Alydidae) based on de novo genome assembly and transcriptome analysis

Food security is important for the ever-growing global population. Soybean, Glycine max (L.) Merr., is cultivated worldwide providing a key source of food, protein and oil. Hence, it is imperative to maintain or to increase its yield under different conditions including challenges caused by abiotic and biotic stresses. In recent years, the soybean pod-sucking stinkbug Riptortus pedestris has emerged as an important agricultural insect pest in East, South and Southeast Asia. Here, we present a genomics resource for R. pedestris including its genome assembly, messenger RNA (mRNA) and microRNA (miRNA) transcriptomes at different developmental stages and from different organs. As insect hormone biosynthesis genes (genes involved in metamorphosis) and their regulators such as miRNAs are potential targets for pest control, we analyzed the sesquiterpenoid (juvenile) and ecdysteroid (molting) hormone biosynthesis pathway genes including their miRNAs and relevant neuropeptides. Temporal gene expression changes of these insect hormone biosynthesis pathways were observed at different developmental stages. Similarly, a diet-specific response in gene expression was also observed in both head and salivary glands. Furthermore, we observed that microRNAs (bantam, miR-14, miR-316, and miR-263) of R. pedestris fed with different types of soybeans were differentially expressed in the salivary glands indicating a diet-specific response. Interestingly, the opposite arms of miR-281 (-5p and -3p), a miRNA involved in regulating development, were predicted to target Hmgs genes of R. pedestris and soybean, respectively. These observations among others highlight stinkbug's responses as a function of its interaction with soybean. In brief, the results of this study not only present salient findings that could be of potential use in pest management and mitigation but also provide an invaluable resource for R. pedestris as an insect model to facilitate studies on plant-pest interactions.

Identification and expression patterns of somatic piRNAs and PIWI genes in Riptortus pedestris (Hemiptera: Alydidae)

RNA interference (RNAi)-based gene silencing is a feasible and sustainable technology for the management of hemipteran pests by double-stranded RNA involvement, including small-interfering RNA, microRNA, and Piwi-interacting RNA (piRNA) pathways, that may help to decrease the usage of chemical insecticides. However, only a few data are available on the somatic piRNAs and their biogenesis genes in Riptortus pedestris, which serves as a significant pest of soybean (Glycine max). In this study, two family members of the PIWI gene were identified and characterized in R. pedestris, containing Argonaute3 (RpAgo3) and Aubergine (RpAub) genes with conserved protein domains, and their clusters were validated by phylogenetic analysis. In addition, they were widely expressed in all developmental stages of the whole body of R. pedestris and had lower expression levels in R. pedestris guts under different rearing conditions based on previous transcriptome sequencing. Furthermore, abundant clean reads were filtered to a total number of 45,998 piRNAs with uridine bias at the first nucleotide (nt) position and 26-32 nt in length by mapping onto the reference genome of R. pedestris according to our previous whole-transcriptome sequencing. Finally, our data revealed that gut bacterial changes were significantly positively or negatively associated with differentially expressed piRNAs among the five comparison groups with Pearson correlation analysis. In conclusion, these findings paved new avenues for the application of RNAi-based biopesticides for broad-spectrum hemipteran pest control.

Comprehensive insights into host-pathogen interaction between brown planthopper and a fungal entomopathogen by dual RNA sequencing

BACKGROUND

The brown planthopper (BPH) is one of the most destructive pests of rice, causing tremendous yield and economic losses every year. The fungal entomopathogen Metarhizium anisopliae was previously proved to have great potential for BPH biocontrol. Genome-wide insight into the insect-fungus interaction is crucial for genetic improvement of M. anisopliae to enhance its virulence to BPH but still has been poorly explored.

RESULTS

Using dual RNA-seq approach, we present here a global view of host and fungal gene expressions in BPH adults during the fungal infection. The results revealed that BPH could initiate strong defense responses against the fungal attack by upregulating the expressions of a large number of genes, including genes involved in cuticle formation, immune response, cell detoxification and biomacromolecule metabolism. Correspondingly, the fungal entomopathogen could induce a series of genes to infect and modulate BPH, including genes involved in fungal penetration, invasive growth, stress resistance and virulence. Three host defense-related genes (NlPCE4, NlPOD1 and NlCYP4DE1) were chosen for further function analysis. RNAi-mediated knockdown of NlPCE4 caused a significant decrease in BPH survival, but no obvious effects on the survival rates were detected by the suppression of NlPOD1 and NlCYP4DE1. Combination of dsRNA injection and fungal infection could significantly enhance the BPH-killing speed, as synergistic mortalities were observed in co-treatments of RNAi and M. anisopliae infection.

CONCLUSION

Our study provides a comprehensive insight into molecular mechanisms of host-pathogen interaction between BPH and M. anisopliae and contributes to future development of new efficient biocontrol strategies for BPH biocontrol.

Screening Candidate Effectors of the Bean Bug Riptortus pedestris by Proteomic and Transcriptomic Analyses

The damage of Riptortus pedestris is exceptional by leading soybean plants to keep green in late autumn. Identification of the salivary proteins is essential to understand how the pest-plant interaction occurs. Here, we have tried to identify them by a combination of proteomic and transcriptomic analyses. The transcriptomes of salivary glands from R. pedestris males, females and nymphs showed about 28,000 unigenes, in which about 40% had open reading frames (ORFs). Therefore, the predicted proteins in the transcriptomes with secretion signals were obtained. Many of the top 1,000 expressed transcripts were involved in protein biosynthesis and transport, suggesting that the salivary glands produce a rich repertoire of proteins. In addition, saliva of R. pedestris males, females and nymphs was collected and proteins inside were identified. In total, 155, 20, and 11 proteins were, respectively, found in their saliva. We have tested the tissue-specific expression of 68 genes that are likely to be effectors, either because they are homologs of reported effectors of other sap-feeding arthropods, or because they are within the top 1,000 expressed genes or found in the salivary proteomes. Their potential functions in regulating plant defenses were discussed. The datasets reported here represent the first step in identifying effectors of R. pedestris.

Immune-Related Genes of Megalurothrips usitatus (Bagrall) Against Beauveria brongniartii and Akanthomyces attenuatus Identified Using RNA Sequencing

Megalurothrips usitatus (Bagrall) is an important pest of legumes worldwide, causing great economic loss every year. Beauveria brongniartii and Akanthomyces attenuatus have shown considerable pathogenicity against M. usitatus in our previous studies. The medial lethal concentration (LC₅₀) and the sublethal lethal concentration (LC₂₅) of B. brongniartii isolate SB010 against M. usitatus were 8.38 × 10⁵ and 1.73 × 10⁵ conidia mL⁻¹, respectively, whereas those of A. attenuatus isolate SCAUDCL-53 against M. usitatus were 4.37 × 10⁵ and 2.97 × 10⁴ conidia mL⁻¹, respectively. This study reports the transcriptome-based explanation of the stress responses of M. usitatus following the application of B. brongniartii and A. attenuatus. The analysis of the transcriptomic data revealed the expression of 254, 207, 195, and 234 immunity-related unigenes by M. usitatus in response to B. brongniartii LC₅₀ (SB1), B. brongniartii LC₂₅ (SB2), A. attenuatus LC₅₀ (V1), and A. attenuatus LC₂₅ (V2), respectively. The biological function and metabolic pathway analyses showed that these unigenes were mainly related to pattern recognition receptors, information transduction factors, and reaction factors, such as scavenger receptor, cytochrome b5, cuticle protein, lysozyme, and serine protease.

Comparative Transcriptomic Analysis of Riptortus pedestris (Hemiptera: Alydidae) to Characterize Wing Formation across All Developmental Stages

Simple Summary

Riptortus pedestris is a widely distributed pest insect in East Asia that causes considerable economic losses. In this study, we applied the Illumina HiSeq6000 platform to construct and sequence the transcriptome libraries of R. pedestris during all life stages. First, a total of 60,058 unigenes were assembled from raw data, and then annotated and classified with various databases. Furthermore, different numbers of differentially expressed genes were calculated by pairwise comparisons of all life stages, and some of these DEGs were associated with various functions by GO and KEGG analysis. Additionally, a total number of 35,158 SSRs and 715,604 SNPs were identified from all the transcriptome libraries. Finally, we analyzed ten wing formation-related signaling pathways, and detected the molecular and expression characterization of five wing development-related genes by qRT-PCR for all developmental stages of R. pedestris. Collectively, all these data may pave the avenue for exploring the developmental processes of hemimetabolous insects and pest management.

Abstract

Riptortus pedestris (Hemiptera: Alydidae) is a major agricultural pest in East Asia that causes considerable economic losses to the soybean crop each year. However, the molecular mechanisms governing the growth and development of R. pedestris have not been fully elucidated. In this study, the Illumina HiSeq6000 platform was employed to perform de novo transcriptome assembly and determine the gene expression profiles of this species across all developmental stages, including eggs, first-, second-, third-, fourth-, and fifth-instar nymphs, and adults. In this study, a total of 60,058 unigenes were assembled from numerous raw reads, exhibiting an N50 length of 2126 bp and an average length of 1199 bp, and the unigenes were annotated and classified with various databases, such as the Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Groups (COG), and Gene Ontology (GO). Furthermore, various numbers of differentially expressed genes (DEGs) were calculated through pairwise comparisons of all life stages, and some of these DEGs were associated with immunity, metabolism, and development by GO and KEGG enrichment. In addition, 35,158 simple sequence repeats (SSRs) and 715,604 potential single nucleotide polymorphisms (SNPs) were identified from the seven transcriptome libraries of R. pedestris. Finally, we identified and summarized ten wing formation-related signaling pathways, and the molecular properties and expression levels of five wing development-related genes were analyzed using quantitative real-time PCR for all developmental stages of R. pedestris. Taken together, the results of this study may establish a foundation for future research investigating developmental processes and wing formation in hemimetabolous insects and may provide valuable data for pest control efforts attempting to reduce the economic damage caused by this pest.

Construction and Selection of an Entomopathogenic Fungal Library From Soil Samples for Controlling Spodoptera litura

Chemical pesticides have been used for pest control for many decades, but they cause serious problems, including insecticide resistance, secondary pest resurgence, and negative environmental impacts. Therefore, sustainable alternatives to chemical pesticides are necessary for pest control. Entomopathogenic fungi (EPF) are natural epizootic pathogens of insects, and some of them have been used as microbial biocontrol agents. Herein, we attempted to construct an entomopathogenic fungal library (EFLib) via the “Tenebrio molitor pathogenicity-based fungal collection method (TmPC)” to select EPF for control of the serious agricultural pest Spodoptera litura. A total of 172 soil samples were collected in northern and central Taiwan for the EFLib construction. The isolation efficiency of TmPC was 64.02%. The EFLib consisted of 101 isolates and was designated as the National Chung Hsing University (NCHU) EFLib. Among these isolates, 26 showed high virulence (mortality = 100%) to T. molitor larvae. Based on the results of molecular identification, the highly virulent isolates belonged to seven genera, including Beauveria, Clonostachys, Fusarium, Cordyceps, Penicillium, Purpureocillium, and Metarhizium. To evaluate the potential of these isolates for Spodoptera litura control, 12 isolates were selected for pathogenicity screening against S. litura larvae. A total of six EPF isolates belonging to the genera Beauveria and Metarhizium showed rapid eradication of the S. litura larvae. To rank the potential of these fungal strains for pest management, the six isolates were subjected to thermotolerance and conidial production assays, and a novel effective conidia number (ECN) formula was applied. The results indicated that the ECN index of Beauveria australis (NCHU-113) was much higher than that of Metarhizium isolates. However, among the Metarhizium isolates, NCHU-95 showed the highest ECN index. Altogether, NCHU-69 and NCHU-113 should be further tested in field trials. To our knowledge, this is the first attempt to integrate pathogenicity or virulence and ECN data into EPF screening and ranking, providing a baseline for mass selection of potential EPF strains for further applications.

Gene diversity explains variation in biological features of insect killing fungus, Beauveria bassiana

Beauveria bassiana is a species complex whose isolates show considerable natural genetic variability. However, little is known about how this genetic diversity affects the fungus performance. Herein, we characterized the diversity of genes involved in various mechanisms of the infective cycle of 42 isolates that have different growth rates, thermotolerance and virulence. The analysed genes showed general genetic diversity measured as non-synonymous changes (NSC) and copy number variation (CNV), with most of them being subjected to positive episodic diversifying selection. Correlation analyses between NSC or CNV and the isolate virulence, thermotolerance and growth rate revealed that various genes shaped the biological features of the fungus. Lectin-like, mucin signalling, Biotrophy associated and chitinase genes NSCs correlated with the three biological features of B. bassiana. In addition, other genes (i.e. DNA photolyase and cyclophilin B) that had relatively conserved sequences, had variable CNs across the isolates which were correlated with the variability of either virulence or thermotolerance of B. bassiana isolates. The data obtained is important for a better understanding of population structure, ecological and potential impact when isolates are used as mycoinsecticides and can justify industrialization of new isolates.

Beauveria bassiana ERL836 and JEF-007 with similar virulence show different gene expression when interacting with cuticles of western flower thrips, Frankniella occidentalis

Background

Insect-killing fungal species, Beauveria bassiana, is as an environment-friendly pest management tool, and many isolates are on the track of industrialization. However, some of B. bassiana isolates show similar morphology and virulence against insect pests, and so it is hard to differentiate them. Herein we used two patented isolates, ERL836 and JEF-007, and investigated their virulence against western flower thrips, Frankliniella occidentalis, and further analyzed genome structures and transcriptional responses when interacting with cuticles of thrips to see possible differences on the initial step of fungal infection.

Results

The two isolates showed no significant differences in fungal growth, conidial production, and virulence against thrips, and they were structurally similar in genome. But, in transcription level, ERL836 appeared to infect thrips easily, while JEF-007 appeared to have more difficulty. In the GO analysis of ERL836 DEGs (differentially expressed genes), the number of up-regulated genes was much larger than that of down-regulated genes, when compared to JEF-007 DEGs (more genes down-regulated). Interestingly, in the enrichment analysis using shared DEGs between two infecting isolates, plasma membrane-mediated transporter activity and fatty acid degradation pathway including cytochrome P450 were more active in infecting ERL836.

Conclusion

The two B. bassiana isolates had similar morphology and virulence as well as genome structure, but in transcription level they differently interacted with the cuticle of western flower thrips. This comparative approach using shared DEG analysis could be easily applied to characterize the difference of the two B. bassiana isolates, JEF-007 and ERL836.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07253-y.