De novo characterization of microRNAs in oriental fruit moth Grapholita molesta and selection of reference genes for normalization of microRNA expression

Molecular characterization and functional analysis of the ecdysone receptor isoform (EcR) from the oriental fruit moth Grapholita molesta (Lepidoptera: Tortricidae)

20-Hydroxyecdysone (20E) plays a vital role in a series of biological processes, via the nuclear receptors, EcR/USP by activating the ecdysone regulatory cascade. To clarify the role of EcR during the development of Grapholita molesta, the complementary DNA of ecdysone receptor isoform B1 (GmEcR-B1) was obtained from the transcriptome of G. molesta and verified by PCR. Alignment analysis revealed that the deduced protein sequence of GmEcR-B1 was highly homologous to EcR proteins identified in other lepidopteran species, especially the EcR-B1 isoform in Spodoptera litura. Quantitative real-time PCR showed that GmEcRs was expressed at all test developmental stages, and the expression level of GmEcRs was relatively higher during the period of the 3rd day of fifth instar larvae to 2nd of pupa than those in other stages. Moreover, the messenger RNA of GmEcRs was much more strongly expressed in the Malpighian tubule and epidermis than those in other tissues, which suggests that this gene may function in a tissue-specific manner during larval development. Silencing of GmEcRs could significantly downregulate the transcriptional level of ecdysone-inducible genes and result in increased mortality during metamorphosis and prolonged prepupal duration. Taken together, the present results indicate that GmEcRs may directly or indirectly affect the development of G. molesta.

Molecular characterization and functional analysis of a beta-1,3-glucan recognition protein from oriental fruit moth Grapholita molesta (Lepidoptera: Tortricidae)

The beta-1,3-glucan recognition protein (BGRP) is an important pattern recognition protein (PRP), which plays an important role in immune recognition and signaling pathway of insect innate immunity. Herein, a BGRP gene was obtained from the transcriptome of Grapholita molesta and its expression was verified by PCR. The full cDNA of the GmBGRP gene was 1691 bp encoding 486 amino acid residues. The calculated molecular mass of the mature protein was 54.83 kDa with an estimated pI of 6.14. The amino acid sequence of GmBGRP was highly homologous to BGRPs of other lepidopterans including Leguminivora glycinivorella BGRP-3. Expression profile of GmBGRP at different developmental stages and different tissues of 5th instar larvae showed that the expression level of this gene tends to slightly increase and then decrease at the adult stage, with the highest at the pupa stage; and mainly expressed in the epidermis, fat body and hemocytes compared with other tissues. In addition, we investigated the transcription levels of other immune-related genes, such as Serine-1, Serine-2, Serine-3, Serpin, SRCB (scavenger receptor gene), Toll, PPO (prophenoloxidase) upon GmBGRP gene silencing, indicating that GmBGRP expression is associated with immune responses of G. molesta. This was further supported by the upregulation of the mRNA level of GmBGRP following fungal infection. Taken together, these results provide experimental evidence for the role of GmBGRP gene in immune defense in G. molesta larvae.

Molecular characterization of insulin receptor (IR) in oriental fruit moth, Grapholita molesta (Lepidoptera: Tortricidae), and elucidation of its regulatory roles in glucolipid homeostasis and metamorphosis through interaction with miR-982490

As an intermediate molecule in the Insulin/Insulin-like growth factor signalling pathway (IIS), the insulin receptor (IR) plays vital roles linking nutritional signals to the downstream regulation of metabolic homeostasis, development, metamorphosis, reproduction and stress responses. In the present study, we describe the molecular characteristics of IR in the cosmopolitan fruit boring pest, Grapholita molesta, and its predicted posttranscription regulator miR-982490, and elucidate its regulatory roles in glucolipid homeostasis and metamorphosis. Phylogenetic and domain analyses indicate that lepidopteran IRs normally cluster within families, and that four main domains are conserved in GmIR and those of other Lepidoptera. Bio-informatic prediction, synchronic expression profile evaluation and dual luciferase reporter assays indicated negative regulation of GmIR by miR-982490. Injection of miR-982490 agomir into fifth instar larvae yielded effects similar to dsGmIR injection, resulting in enhanced levels of trehalose and triglyceride in haemolymph, and reduced pupation success and pupal weight, both of which could be rescued by co-injection of dsGmIR and miR-982490 antagomir. We infer that GmIR regulates glucolipid homeostasis and affects G. molesta metamorphosis via interactions with its posttranscriptional regulator miR-982490. This study expands our understanding of the regulatory network of IIS in insect nutritional homeostasis and development.

RNA Sequencing Reveals the Potential Adaptation Mechanism to Different Hosts of Grapholita molesta

Grapholita molesta is an important fruit tree worldwide pest which feeds on hosts extensively and does serious harm. In this paper, the growth and development parameters and protease activities of G. molesta fed on different hosts were compared. Using Illumina RNA sequencing technology, 18 midgut samples from five different hosts (apple, pear, plum, peach and peach shoots) and artificial diet were sequenced and compared with the reference genome, resulting in 15269 genes and 2785 predicted new genes. From 15 comparative combinations, DEGs were found from 286 to 4187 in each group, with up-regulated genes from 107 to 2395 and down-regulated genes from 83 to 2665. KEGG pathway analysis showed that DEGs were associated with amino acid metabolism, starch and sucrose metabolism, carbohydrate metabolism, and hydrolase activity. A total of 31 co-expression gene modules of different hosts were identified by WGCNA. qRT-PCR showed that the expression pattern of the trypsin gene was consistent with RNA sequencing. In this study, growth and development parameters, protease activity, DEGs, enrichment analysis and qRT-PCR were combined to reveal the adaptation process to different hosts of G. molesta in many aspects. The results of this study provide a basis for further exploration of the molecular mechanism of host adaptation of G. molesta.

Selection and Evaluation of Reference Genes for miRNA Expression Analysis in Bemisia tabaci Under Insecticide Tolerance

A growing number of studies have focused on the microRNA (miRNA) expression in Bemisia tabaci, one devastating agricultural insect pest of the tropical and subtropical areas for which the primary means of control are insecticides. In studying the genetic underpinnings of insecticide resistance, the choice of stable reference genes for normalizing data plays a key role to acquire unbiased expression profile results from quantitative real-time PCR (qPCR) analysis. Expression profiles of 11 selected reference genes were determined systematically in B. tabaci exposure to 11 insecticides. Furthermore, we assessed the stability of all the selected candidates in relation to other variables including sex, tissue type, and developmental stage. Candidate reference gene validation was conducted by analyzing the let-7-5p expression under various experimental treatments. Five programs BestKeeper, NormFinder, geNorm, △Ct, and RefFinder were applied to verify the accuracy of the selected candidates. Our results displayed that the best choices of the selected candidates for pymetrozine, sulfoxaflor, flonicamid, cyantraniliprole, afidopyropen, and deltamethrin treatment were miR-1-3p and miR-100-5p, U6 and miR-100-5p were best for chlorpyrifos and imidacloprid treatments, and U6 and miR-1-3p were best for flupyradifurone and β-cypermethrin treatments. The reference genes miR-624, miR-252, and miR-275 worked best in adult tissues, miR-100-5p and miR-1-3p worked best in either sex, and miR-624 and miR-11 were best to use across developmental stages. Not even one reference gene was found to be suitable for all experimental conditions. Our results contributed to the growing body of the literature on qPCR reference gene selection under various experimental conditions and facilitate further investigation on gene expression changes in B. tabaci, resulting from pesticide exposure.

The developmentally dynamic microRNA transcriptome of Glossina pallidipes tsetse flies, vectors of animal trypanosomiasis

Summary

MicroRNAs (miRNAs) are single stranded gene regulators of 18-25 bp in length. They play a crucial role in regulating several biological processes in insects. However, the functions of miRNA in Glossina pallidipes, one of the biological vectors of African animal trypanosomosis in sub-Saharan Africa, remain poorly characterized. We used a combination of both molecular biology and bioinformatics techniques to identify miRNA genes at different developmental stages (larvae, pupae, teneral and reproductive unmated adults, gravid females) and sexes of G. pallidipes. We identified 157 mature miRNA genes, including 12 novel miRNAs unique to G. pallidipes. Moreover, we identified 93 miRNA genes that were differentially expressed by sex and/or in specific developmental stages. By combining both miRanda and RNAhybrid algorithms, we identified 5550 of their target genes. Further analyses with the Gene Ontology term and KEGG pathways for these predicted target genes suggested that the miRNAs may be involved in key developmental biological processes. Our results provide the first repository of G. pallidipes miRNAs across developmental stages, some of which appear to play crucial roles in tsetse fly development. Hence, our findings provide a better understanding of tsetse biology and a baseline for exploring miRNA genes in tsetse flies.

Availability and implementation

Raw sequence data are available from NCBI Sequence Read Archives (SRA) under Bioproject accession number PRJNA590626.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Laboratory and field studies supporting augmentation biological control of oriental fruit moth, Grapholita molesta (Lepidoptera: Tortricidae), using Trichogramma dendrolimi (Hymenoptera: Trichogrammatidae)

BACKGROUND

Oriental fruit moth, Grapholita molesta (Busck) (Lepidoptera: Tortricidae) is a cosmopolitan pests of fruit trees in temperate regions. Control of G. molesta is challenging because larvae feed in concealed locations and have evolved resistance to many pesticides. We compared three commercially available species of Trichogramma for efficacy against G. molesta in the laboratory assays and tested releases of the promising species in a pear orchard.

RESULTS

Laboratory assays indicated that Trichogramma dendrolimi parasitized G. molesta at the highest rate. Parasitoids took longer to oviposit in older host eggs, and fewer eggs were parasitized when they were more than 3 days old. Field tests produced ca 60% cumulative parasitism of sentinel G. molesta eggs with one release of Trichogramma dendrolimi, with most parasitism occurring within 24 h. Female wasps dispersed up to 12 m from release points with a bias toward upwind movement.

CONCLUSION

We calculated that, for each generation of G. molesta, 900 000 wasps per hectare, in three releases 3 days apart, each distributed among release points 10 m apart, would reduce fruit damage by half in an orchard where 50% of fruit would otherwise be damaged. Although augmentation of Trichogramma dendrolimi is a viable tactic for reducing G. molesta populations and fruit damage, it will require integration with other compatible control tactics in order to provide commercially acceptable levels of control in orchards experiencing significant pest pressure. © 2021 Society of Chemical Industry.

Reference gene selection for transcriptional profiling in Cryptocercus punctulatus, an evolutionary link between Isoptera and Blattodea

The subsocial life style and wood-feeding capability of Cryptocercus gives us an evolutionary key to unlock some outstanding questions in biology. With the advent of the Genomics Era, there is an unprecedented opportunity to address the evolution of eusociality and the acquisition of lignocellulases at the genetic level. However, to quantify gene expression, an appropriate normalization strategy is warranted to control for the non-specific variations among samples across different experimental conditions. To search for the internal references, 10 housekeeping genes from a gut transcriptome of a wood-feeding cockroach, Cryptocercus punctulatus, were selected as the candidates for the RT-qPCR analysis. The expression profiles of these candidates, including ACT, EF1α, GAPDH, HSP60, HSP70, αTUB, UBC, RPS18, ATPase and GST, were analyzed using a panel of analytical tools, including geNorm, NormFinder, BestKeeper, and comparative ΔC_T method. RefFinder, a comprehensive ranking system integrating all four above-mentioned algorithms, rated ACT as the most stable reference gene for different developmental stages and tissues. Expression analysis of the target genes, Hex-1 and Cell-1, using the most or the least appropriate reference genes and a single or multiple normalizers signified this research. Our finding is the first step toward establishing a standardized RT-qPCR analysis in Cryptocercus.

Comparison of gut bacterial communities and their associations with host diets in four fruit borers

BACKGROUND

Microbiota that live in the gut of insects have a wide range of effects on host nutrition, physiology, and behavior. They may shape the adaptation of their hosts to different habitats and lifestyles. To characterize the gut microbiota of fruit borers comprehensively, we compared bacterial communities among Grapholita molesta, Conogethes punctiferalis, Carposina sasakii, and Cydia pomonella, which are serious lepidopteran pests. We selected G. molesta as a representative pest to more explicitly test the influence of host dietary niche on the insect gut microbiome, and compared the bacterial microbial communities of G. molesta fed different diets (peach shoots and apple) using Illumina high-throughput sequencing technology.

RESULTS

The results show that Proteobacteria and Firmicutes are dominant in their gut microbiota. The C. sasakii had the highest richness values and G. molesta (shoot-feeding) had the highest diversity, whereas C. pomonella and G. molesta (fruit-feeding) had the lowest bacterial richness and diversity, respectively. The ANOSIM analysis revealed significant differences in the structure of gut microbiota among different insects. In addition, G. molesta with a different feeding diet had significant differences in gut microbiota composition. PICRUSt analysis indicated that most functional prediction categories were related to metabolism.

CONCLUSION

Our results show that gut microbiota composition is affected significantly not only by host species but also host diets. An enhanced understanding of these herbivore-associated microbial symbionts is essential for understanding the biology and ecology of the host insect, and may offer new possibilities to improve integrated pest-management strategies for efficient control of fruit borers. © 2019 Society of Chemical Industry.

Molecular Characterization of Primary Juvenile Hormone Responders Methoprene-Tolerant (Met) and Krüppel Homolog 1 (Kr-h1) in Grapholita molesta (Lepidoptera: Tortricidae) with Clarification of Their Roles in Metamorphosis and Reproduction

Methoprene-tolerant (Met) is a putative JH intracellular receptor that transduces JH signal by activation of the inducible Krüppel homolog 1 (Kr-h1). We analyzed the gene sequences of Met and Kr-h1 and their patterns of expression in Grapholita molesta (Busck) immature and adult stages in order to better understand the roles of these primary JH responders in regulating the metamorphosis and reproduction of this global pest of fruit crops. The deduced amino acid sequences of both GmMet and GmKr-h1 were highly homologous to those of other Lepidoptera, especially the cotton bollworm, Helicoverpa armigera (Hubner) (Lepidoptera: Noctuidae). Peak expression of GmMet occurred during the last 3 to 5 d of the final instar, followed by that of GmKr-h1, in the last 3 d of final instar. Similar patterns of GmMet and GmKr-h1 expression were detected across various tissue types in the fifth-instar larvae, with the highest expression observed in the head, followed by the epidermis, and the fat body. When expression of GmMet and GmKr-h1 was knocked down via dsRNA injection in the fifth instar, the results were increased larval mortality, abnormal pupation, delayed pupal duration, reduced adult emergence, extended preoviposition period, and reduced fecundity. We infer that both GmMet and GmKr-h1 participated in regulation of metamorphosis and reproduction in G. molesta, the former acting upstream of the latter, and could present biorational targets for novel pest control compounds.

High-throughput sequencing analysis of microRNAs in gills of red swamp crayfish, Procambarus clarkii infected with white spot syndrome virus

MicroRNAs (miRNAs) are important posttranscriptional regulators. They play an important role in the antiviral innate immunity of invertebrates. In the present study, high-throughput small RNAs Illumina sequencing systems were carried out to identify differentially expressed miRNAs (DEMs) in the gills of Procambarus clarkii, which was challenged with white spot syndrome virus (WSSV). Our results identified 11,617 known and 6 novel miRNAs in normal group (NG) and WSSV-challenged group (WG) small RNA libraries. Additionally, 27 DEMs were shown to participate in the antiviral innate immunity of P. clarkii and were significantly upregulated or downregulated. In addition, the results of the KEGG pathway prediction of the DEMs target genes showed that putative target genes of these 27 DEMs were related mainly to the RNA transport pathway, tight junction pathway, mRNA surveillance pathway, regulation actin cytoskeleton pathway, focal adhesion pathway, and MAPK signaling pathway. These results provide important information for future studies about the antiviral innate immunity of crustaceans.

Evaluation of Reference Genes for Real-Time Quantitative PCR Analysis in Larvae of Spodoptera litura Exposed to Azadirachtin Stress Conditions

Azadirachtin is an efficient and broad-spectrum botanical insecticide against more than 150 kinds of agricultural pests with the effects of mortality, antifeedant and growth regulation. Real-time quantitative polymerase chain reaction (RT-qPCR) could be one of the powerful tools to analyze the gene expression level and investigate the mechanism of azadirachtin at transcriptional level, however, the ideal reference genes are needed to normalize the expression profiling of target genes. In this present study, the fragments of eight candidate reference genes were cloned and identified from the pest Spodoptera litura. In addition, the expression stability of these genes in different samples from larvae of control and azadirachtin treatments were evaluated by the computational methods of NormFinder, BestKeeper, Delta CT, geNorm, and RefFinder. According to our results, two of the reference genes should be the optimal number for RT-qPCR analysis. Furthermore, the best reference genes for different samples were showed as followed: EF-1α and EF2 for cuticle, β-Tubulin and RPL7A for fat body, EF2 and Actin for midgut, EF2 and RPL13A for larva and RPL13A and RPL7A for all the samples. Our results established a reliable normalization for RT-qPCR experiments in S. litura and ensure the data more accurate for the mechanism analysis of azadirachtin.

Identification and characterization of intestine microRNAs and targets in red swamp crayfish, Procambarus clarkii infected with white spot syndrome virus

MicroRNAs (miRNAs) are small non-coding endogenous RNA molecules that play important roles in the innate immunity system of invertebrates, especially in the aspect of antivirus. In the present study, high-throughput small RNA Illumina sequencing systems were used to identify differentially expressed miRNAs (DEMs) from the intestines of Procambarus clarkii that were infected with white spot syndrome virus (WSSV). As a result, 39 known and 12 novel miRNAs were identified in both NG and WG small RNA libraries. Seven DEMs were determined to be involved in the antiviral innate immunity in the intestines of P. clarkii. The results of the target gene predictions of the DEMs showed that the putative target genes of these 7 DEMs are related to tight junctions, vascular smooth muscle contraction regulation of the actin cytoskeleton, focal adhesion, RNA transport, mRNA surveillance, viral carcinogenesis, and Salmonella infection. These results provide theoretical insights for future studies on the antiviral immunity of crustaceans.

Identification and characterization of lymph organ microRNAs in red swamp crayfish, Procambarus clarkii infected with white spot syndrome virus

MicroRNAs (miRNAs) were important post-transcriptional regulators and played vital roles in innate immunity system of invertebrates, especially in the aspect of antivirus. In this study, using high-throughput small RNAs Illumina sequencing system, differentially expressed miRNAs (DEMs) from lymph organs in red swamp crayfish, Procambarus clarkii, infected with white spot syndrome virus, were identified. As a result, 32 known miRNAs and 7 novel miRNAs were identified in crayfish lymph organ small RNAs library of NG and WG. Among them, 7 differentially expressed miRNAs (DEMs) were predicted to be involved in the lymph organ antiviral innate immunity of P. clarkii. Besides, the results showed that putative target genes of these DEMs were related with tight junction, RNA transport, regulation of actin cytoskeleton, focal adhesion, vascular smooth muscle contraction, mRNA surveillance pathway, NOD-like receptor signaling pathway, leukocyte transendothelial migration, and protein processing in endoplasmic reticulum. These results might provide the guiding theoretical foundation for future studies about crustaceans' antiviral innate immunity.