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

Shared phylogeographic patterns and environmental responses of co-distributed soybean pests: Insights from comparative phylogeographic studies of Riptortus pedestris and Riptortus linearis in the subtropics of East Asia

Comparative phylogeographic studies of closely related species sharing co-distribution areas can elucidate the role of shared historical factors and environmental changes in shaping their phylogeographic pattern. The bean bugs, Riptortus pedestris and Riptortus linearis, which both inhabit subtropical regions in East Asia, are recognized as highly destructive soybean pests. Many previous studies have investigated the biological characteristics, pheromones, chemicals and control mechanisms of these two pests, but few studies have explored their phylogeographic patterns and underlying factors. In this study, we generated a double-digest restriction site-associated DNA sequencing (ddRAD-seq) dataset to investigate phylogeographic patterns and construct ecological niche models (ENM) for both Riptortus species. Our findings revealed similar niche occupancies and population genetic structures between the two species, with each comprising two phylogeographic lineages (i.e., the mainland China and the Indochina Peninsula clades) that diverged approximately 0.1 and 0.3 million years ago, respectively. This divergence likely resulted from the combined effects of temperatures variation and geographical barriers in the mountainous regions of Southwest China. Further demographic history and ENM analyses suggested that both pests underwent rapid expansion prior to the Last Glacial Maximum (LGM). Furthermore, ENM predicts a northward shift of both pests into new soybean-producing regions due to global warming. Our study indicated that co-distribution soybean pests with overlapping ecological niches and similar life histories in subtropical regions of East Asia exhibit congruent phylogeographic and demographic patterns in response to shared historical biogeographic drivers.

The ncRNA-mediated regulatory networks of defensins and lysozymes in Riptortus pedestris: involvement in response to gut bacterial disturbances

Insects depend on humoral immunity against intruders through the secretion of antimicrobial peptides (AMPs) and immune effectors via NF-κB transcription factors, and their fitness is improved by gut bacterial microbiota. Although there are growing numbers of reports on noncoding RNAs (ncRNAs) involving in immune responses against pathogens, comprehensive studies of ncRNA-AMP regulatory networks in Riptortus pedestris, which is one of the widely distributed pests in East Asia, are still not well understood under feeding environmental changes. The objective of this study employed the whole-transcriptome sequencing (WTS) to systematically identify the lncRNAs (long noncoding RNA) and circRNAs (circular RNA) and to obtain their differential expression from the R. pedestris gut under different feeding conditions. Functional annotation indicated that they were mainly enriched in various biological processes with the GO and KEGG databases, especially in immune signaling pathways. Five defensin (four novel members) and eleven lysozyme (nine novel members) family genes were identified and characterized from WTS data, and meanwhile, phylogenetic analysis confirmed their classification. Subsequently, the miRNA-mRNA interaction network of above two AMPs and lncRNA-involved ceRNA (competing endogenous RNA) regulatory network of one lysozyme were predicted and built based on bioinformatic prediction and calculation, and the expression patterns of differentially expressed (DE) defensins, and DE lysozymes and related DE ncRNAs were estimated and selected among all the comparison groups. Finally, to integrate the analyses of WTS and previous 16S rRNA amplicon sequencing, we conducted the Pearson correlation analysis to reveal the significantly positive or negative correlation between above DE AMPs and ncRNAs, as well as most changes in the gut bacterial microbiota at the genus level of R. pedestris. Taken together, the present observations provide great insights into the ncRNA regulatory networks of AMPs in response to rearing environmental changes in insects and uncover new potential strategies for pest control in the future.

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.

Evaluation of Reference Genes for Quantitative Real-Time PCR Analysis in the Bean Bug, Riptortus pedestris (Hemiptera: Alydidae)

Quantitative real-time PCR (qRT-PCR) is widely accepted as a precise and convenient method for quantitatively analyzing the expression of functional genes. The data normalization strongly depends upon stable reference genes. The bean bug, Riptortus pedestris (Hemiptera: Alydidae), is a significant pest of leguminous crops and broadly distributed across Southeast Asia. In this study, a total of 16 candidate reference genes (RPL32, RPS23, SDHA, UBQ, UCCR, GST, TATA-box, HSP70, GAPDH, RPL7A, SOD, RPS3, Actin, α-tubulin, AK, and EF1) were carefully chosen in R. pedestris, and their expression levels were assessed across various conditions, including different developmental stages, diverse tissues, temperature treatments, adult age, molting time, and mating status. Following this, the stability of these reference genes was evaluated using four algorithms (ΔCt, GeNorm, NormFinder, and BestKeeper). Ultimately, the comprehensive rankings were determined using the online tool RefFinder. Our results demonstrate that the reference gene for qRT-PCR analysis in R. pedestris is contingent upon the specific experimental conditions. RPL7A and EF1 are optimal reference genes for developmental stages. Furthermore, α-tubulin and EF1 exhibit the most stable expression across various adult tissues. RPL32 and RPL7A exhibit the most stable expression for adult age. For nymph age, RPL32 and SOD display the most stable expression. For temperature conditions, RPS23 and RPL7A were identified as the most suitable for monitoring gene expression. Lastly, we verified the practicability of evaluating expression levels of odorant-binding protein 37 (RpedOBP37) and cytochrome P450 6a2 (RpedCYP6) throughout developmental stages, tissues, and temperature conditions. These findings are a significant addition to the qRT-PCR analysis studies on R. pedestris, serving as a fundamental groundwork for future investigations on stable reference genes in R. pedestris as well as other organisms.

Changes in the gut bacterial community affect miRNA profiles in Riptortus pedestris under different rearing conditions

Insects possess complex and dynamic gut microbial system, which contributes to host nutrient absorption, reproduction, energy metabolism, and protection against stress. However, there are limited data on interactions of host-gut bacterial microbiota through miRNA (microRNA) regulation in a significant pest, Riptortus pedestris. Here, we performed the 16S rRNA amplicon sequencing and small RNA sequencing from the R. pedestris gut under three environmental conditions and antibiotic treatment, suggesting that we obtained a large amount of reads by assembly, filtration and quality control. The 16S rRNA amplicon sequencing results showed that the abundance and diversity of gut bacterial microbiota were significantly changed between antibiotic treatment and other groups, and they are involved in metabolism and biosynthesis-related function based on functional prediction. Furthermore, we identified different numbers of differentially expressed unigenes (DEGs) and differentially expressed miRNAs (DEMs) based on high-quality mappable reads, which were enriched in various immune-related pathways, including Toll-like receptor, RIG-I-like receptor, NOD-like receptor, JAK/STAT, PI3K/Akt, NF-κB, MAPK signaling pathways, and so forth, using GO and KEGG enrichment analysis. Later on, the identified miRNAs and their target genes in the R. pedestris gut were predicted and randomly selected to construct an interaction network. Finally, our study indicated that alterations in the gut bacterial microbiota are significantly positively or negatively associated with DEMs of the Toll/Imd signaling pathway with Pearson correlation analysis. Taken together, the results of our study lay the foundation for further deeply understanding the interactions between the gut microbiota and immune responses in R. pedestris through miRNA regulation, and provide the new basis for pest management in hemipteran pests.

The lncRNA-mediated ceRNA network of Altica viridicyanea is involved in the regulation of the Toll/Imd signaling pathway under antibiotic treatment

Long noncoding RNAs (lncRNAs) play significant roles in the regulation of mRNA expression or in shaping the competing endogenous RNA (ceRNA) network by targeting miRNA. The insect gut is one of the most important tissues due to direct contact with external pathogens and functions in the immune defense against pathogen infection through the innate immune system and symbionts, but there are limited observations on the role of the lncRNA-involved ceRNA network of the Toll/Imd pathway and correlation analysis between this network and bacterial microbiota in the Altica viridicyanea gut. In this research, we constructed and sequenced six RNA sequencing libraries using normal and antibiotic-reared samples, generating a total of 17,193 lncRNAs and 26,361 mRNAs from massive clean data by quality control and bioinformatic analysis. Furthermore, a set of 8,539 differentially expressed lncRNAs (DELs) and 13,263 differentially expressed mRNAs (DEMs), of which related to various immune signaling pathways, such as the Toll/Imd, JAK/STAT, NF-κB, and PI3K-Akt signaling pathways, were obtained between the two experimental groups in A. viridicyanea. In addition, numerous GO and KEGG enrichment analyses were used to annotate the DELs and their target genes. Moreover, six Toll family members and nineteen signal genes from the Toll/Imd signaling pathway were identified and characterized using online tools, and phylogenetic analyses of the above genes proved their classification. Next, a lncRNA-miRNA-mRNA network of the Toll/Imd pathway was built, and it contained different numbers of DEMs in this pathway and related DELs based on prediction and annotation. In addition, qRT-PCR validation and sequencing data were conducted to show the expression patterns of the above DELs and DEMs related to the Toll/Imd signaling pathway. Finally, the correlated investigations between DELs or DEMs of the Toll/Imd signaling pathway and most changes in the gut bacterial microbiota revealed significantly positive or negative relationships between them. The present findings provide essential evidence for innate immune ceRNAs in the beetle gut and uncover new potential relationships between innate immune pathways and the gut bacterial microbiota in insects.

Genome-wide identification of candidate chemosensory receptors in the bean bug Riptortus pedestris (Hemiptera: Alydidae) and the functional verification of its odorant receptor co-receptor (Orco) in recognizing aggregation pheromone

A sophisticated and sensitive olfactory system plays a vital role in the survival and reproduction of insects. Chemosensory receptors are indispensable for the molecular recognition and discrimination of semiochemicals. Riptortus pedestris is a notorious pest of legume plants, resulting in yield losses and quality decreases in soybeans. It is well accepted that R. pedestris highly relies on its olfactory system in detecting aggregation pheromones, host volatiles, and pesticides; however, little research focused on its chemosensory receptors. In the present study, we identified 237 odorant receptors (ORs), 42 gustatory receptors (GRs), and 31 ionotropic receptors (IRs) from the reported genome of R. pedestris, and analyzed their phylogenetic relationship with other hemipteran species. Through the results of RNA-seq and real-time quantitative PCR (qRT-PCR), we found that RpedORs displayed different expression levels in the antennae of R. pedestris at different development stages. To further verify the function of odorant receptor co-receptor (Orco), an obligate and unique insect OR, we silenced RpedOrco by RNA interference (RNAi) method. The results showed that silencing RpedOrco could significantly impair the response to aggregation pheromone in R. pedestris, indicating that RpedOrco plays an essential role in odorant detection. Our results can provide the theoretical foundations for revealing the olfactory recognition mechanism of R. pedestris and help explore and develop novel olfactory-based agents against this pest.

Modelling the current and future potential distribution of the bean bug Riptortus pedestris with increasingly serious damage to soybean

BACKGROUND

The bean bug, Riptortus pedestris, has received intense attention in recent years because of its involvement in increasing outbreaks of staygreen syndrome in soybean (Glycine max (L.)), often causing almost 100% loss of soybean yield in China. However, for this pest of great economic importance, potential current and future distribution patterns and their underlying driving factors remain unclear.

RESULTS

Maxent modelling under climate, elevation and land-use (including the distribution information of G. max) variables showed that the current potential distribution covered a vast geographic range, primarily including most parts of south, South East and east Asia. Under future environmental scenarios, suitable habitat expanded markedly. Areas that would become highly suitable for R. pedestris were primarily located in north-east China and west India. Five bioclimatic (BIO13, BIO08, BIO18, BIO02 and BIO07) and one land-use (C3 annual crops) predictors contributed approximately 95% to the modelling, and analyses of curve responses showed that to a certain extent, R. pedestris preferred relatively high temperature and precipitation. Our results indicate that a high risk of R. pedestris outbreaks is present in parts of Asia, especially in the soybean-growing regions of China, and this risk will continue in the future.

CONCLUSION

The predicted distribution pattern and key regulating factors identified herein could provide a vital reference for developing pest management policies and further alleviate the incidence of staygreen syndrome in soybean. © 2022 Society of Chemical Industry.

Selection and Validation of Reference Genes for RT-qPCR Normalization in Bradysia odoriphaga (Diptera: Sciaridae) Under Insecticides Stress

Bradysia odoriphaga (Diptera: Sciaridae) is the most serious root maggot pest which causes substantial damage to the Chinese chive. Organophosphate (OP) and neonicotinoid insecticides are widely used chemical pesticides and play important roles in controlling B. odoriphaga. However, a strong selection pressure following repeated pesticide applications has led to the development of resistant populations of this insect. To understand the insecticide resistance mechanism in B. odoriphaga, gene expression analysis might be required. Appropriate reference gene selection is a critical prerequisite for gene expression studies, as the expression stability of reference genes can be affected by experimental conditions, resulting in biased or erroneous results. The present study shows the expression profile of nine commonly used reference genes [elongation factor 1α (EF-1α), actin2 (ACT), elongation factor 2α (EF-2α), glucose-6-phosphate dehydrogenase (G6PDH), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal protein L10 (RPL10), ribosomal protein S3 (RPS3), ubiquitin-conjugating enzyme (UBC), and α-tubulin (TUB)] was systematically analyzed under insecticide stress. Moreover, we also evaluated their expression stability in other experimental conditions, including developmental stages, sexes, and tissues. Five programs (NormFinder, geNorm, BestKeeper, RefFinder, and ΔCt) were used to validate the suitability of candidate reference genes. The results revealed that the most appropriate sets of reference genes were RPL10 and ACT across phoxim; ACT and TUB across chlorpyrifos and chlorfluazuron; EF1α and TUB across imidacloprid; EF1α and EF2α across developmental stages; RPL10 and TUB across larvae; EF1α and ACT across tissues, and ACT and G6PDH across sex. These results will facilitate the standardization of RT-qPCR and contribute to further research on B. odoriphaga gene function under insecticides stress.