- Research article
- Open Access
Growth inhibition of Spodoptera frugiperda larvae by camptothecin correlates with alteration of the structures and gene expression profiles of the midgut
BMC Genomics volume 22, Article number: 391 (2021)
Spodoptera frugiperda is a serious pest that causes devastating losses to many major crops, including corn, rice, sugarcane, and peanut. Camptothecin (CPT) is a bioactive secondary metabolite of the woody plant Camptotheca acuminata, which has shown high toxicity to various pests. However, the effect of CPT against S. frugiperda remains unknown.
In this study, bioassays have been conducted on the growth inhibition of CPT on S. frugiperda larvae. Histological and cytological changes were examined in the midgut of larvae fed on an artificial diet supplemented with 1.0 and 5.0 µg/g CPT. The potential molecular mechanism was explored by comparative transcriptomic analyses among midgut samples obtained from larvae under different treatments. A total of 915 and 3560 differentially expressed genes (DEGs) were identified from samples treated with 1.0 and 5.0 µg/g CPT, respectively. Among the identified genes were those encoding detoxification-related proteins and components of peritrophic membrane such as mucins and cuticle proteins. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses indicated that part of DEGs were involved in DNA replication, digestion, immunity, endocrine system, and metabolism.
Our results provide useful information on the molecular basis for the impact of CPT on S. frugiperda and for future studies on potential practical application.
The fall armyworm, Spodoptera frugiperda (Lepidoptera: Noctuidae), is an important insect pest worldwide. The insect can feed on at least 353 plant species including major crops such as corn, rice, soybeans, sugar cane, and cotton [1,2,3]. S. frugiperda is native to tropical and subtropical regions of the Americas, but has been spread to Africa and Asia in recent years . The voracious caterpillar was first found in 2019 in Yunnan Province, China. Since then it has spread rapidly to most parts of the country . Previous studies have been focused on finding control strategies such as various monitoring methods, correct identification of species and strains by genotyping, biological control and chemical application . Effective insecticides for controlling S. frugiperda include pyrethroids, diacyl hydrazides, diamides, and benzoylureas . Extensive application of insecticides caused problems including arise of populations with resistance to insecticides, toxicity to beneficial animals, and harmful effects on human health.
Plants are considered the most abundant natural resource in the world for the identification of chemicals with insecticidal activity [7, 8]. Plant secondary metabolites protect plants from herbivores and are potential candidates for novel insecticides . In fact, many existing insecticides are derivatives of plant metabolites. For example, pyrethrum and nicotine are used as botanical pesticides for pest control for decades . In recent years, the effect of plant secondary metabolites against S. frugiperda has been investigated. For example, the botanical insecticide azadirachtin is very toxic to S. frugiperda, with LC50 values 0.59 and 0.46 mg/L for 2nd and 3rd instar larvae, respectively, under 0.3 % azadirachtin emulsifiable concentrate (EC) . Cedrelone, a metabolite isolated from the Australian red cedar Toona ciliata, is toxic to S. frugiperda larvae as well . The flavonoid rutin extracted from soybean prolongs the development of S. frugiperda larvae, causing reduced larval and pupal viability . The toxicity of extracts from Actinostemon concolor, Piper aduncum, and Ruta graveolens has also been tested against S. frugiperda caterpillars [14,15,16].
Camptothecin (CPT), a pentacyclic quinoline alkaloid isolated from the plant Camptotheca acuminata Decne, is a potent pharmaceutical secondary metabolite with antitumor activities in mammalian cells by targeting intracellular DNA topoisomerase I, resulting in inhibition of nucleic acid synthesis and induction of DNA strand breakage [17, 18]. CPT has also displayed insecticidal activity against several insects, including Drosophila melanogaster, Musca domestica, Mythimna separata, and Spodoptera exigua [10, 19, 20]. Field tests with 0.2 % camptothecin emulsifiable concentrate (EC) have shown high mortality on three important agricultural pests Nilaparvata lugens (Ståhl) Brevicoryne brassicae (L.), and Chilo suppressalis (Walker) . Due to its low water solubility properties, a series of CPT derivatives have been developed through structural modification . Phytophagous mites including Tetranychus urticae, Acaphylla theae and Brevipalpus obovatus were sensitive to the aqueous CPT-Na+ solution under laboratory and field conditions . The toxicity mechanism indicated that CPT inhibits DNA topoisomerase I (topo I) . Besides, CPT up-regulated the expression of programmed cell death protein 11 in Spodoptera litura, which could be involved in apoptosis induction . While the effects of CPT against S. frugiperda and relevant molecular mechanisms remain to be revealed.
The objective of this study is to investigate the adverse effect of CPT against S. frugiperda. Changes in the weight of S. frugiperda larvae were examined after treatments with different CPT concentrations. Histopathological and ultrastructural changes in the midgut of larvae fed on diets containing 1.0 and 5.0 µg/g CPT, respectively, were examined. In addition, comparative transcriptomic analyses were carried out with different midgut samples from larvae under different treatments. Our results indicated that CPT is a growth inhibitor of S. frugiperda larvae and has the potential as an insecticide for controlling this important insect pest in the field.
CPT inhibits S. frugiperda larval growth
To examine any adverse effect of CPT against S. frugiperda, third-instar larvae were fed on artificial diets containing 0, 1.0, 2.5, 5.0, 10, 20, and 30 µg/g CPT, respectively. The weight of larvae for each sample was recorded on 1, 3, 5, and 7 days after treatments. The average weight of larvae fed on CPT-diets for one day showed no significant difference compared with that of controls. Weight loss was observed in larvae fed on CPT diets for 3, 5, and 7 days (Fig. 1). Our results indicated that CPT inhibited the growth of S. frugiperda larvae in a dose-dependent manner.
CPT causes structural damages in S. frugiperda larval midgut
After 7 days of feeding, the larvae from the control group developed to sixth instar larvae, while the larvae treated with CPT grew slowly, and developed only into fourth or fifth instars. Histopathological changes were observed in the larval midgut fed on diets containing 1.0 and 5.0 µg/g CPT for 7 days based on hematoxylin-eosin (HE) staining. As shown in Fig. 2 A, midgut cells were tightly arranged in multiple layers with a thick intestinal wall in control insects. In comparison, many cells were disappeared and only a thin intestinal barrier was observed in the larval midgut fed on a 1.0 µg/g CPT diet (Fig. 2B). The severity of damage to the gut was dose-dependent. In larvae fed on a diet containing 5.0 µg/g CPT, only the basement membrane was left in the intestinal wall of the midgut, and nearly all functional cells disappeared (Fig. 2 C). Similar phenomena were observed in the gut structure under TEM. In control larvae, chromatin was evenly distributed in the nucleus. Mitochondria and endoplasmic reticulum were abundant and distributed evenly in the cytoplasm. Microvilli were ordinally distributed in the gut (Fig. 2D). In contrast, the number of mitochondria and endoplasmic reticulum decreased in midgut cells in larvae treated with 1.0 µg/g CPT. Microvilli were disorganized (Fig. 2E). In larvae fed on the 5.0 µg/g CPT diet, chromatin condensation occurred and chromatins were located close to the nuclear envelope. Microvilli decreased and deformed with large cavities (Fig. 2 F). Our results indicated that CPT had negative effects on the midgut structure of S. frugiperda larvae.
Midguts dissected from larvae treated with CPT (1.0 and 5.0 µg/g) for 7 days and control larvae were used for transcriptomic analyses. The number of raw reads from nine libraries ranged from 43,911,148 to 51,395,872. High quality reads ranged from 43,587,740 to 51,021,846 (Supplement Table 1). Q20 and Q30 refer to the percentage of bases with sequencing quality above 99 and 99.9 % to the total bases in the transcriptome. Values of Q20 and Q30 in each transcriptome were more than 98 and 96 %, respectively (Supplement Table 1). A total of 58,122 unigenes was obtained from the de novo assembling of all combined reads. The length of the unigenes ranged from 201 to 28,483 bp, with an average of 765.28 bp. N50 and GC content of the unigenes were 1358 bp and 40.55 %, respectively. Original data were deposited to the SRA database with the accession number of SRP242660. Transcripts assembled by Trinity were submitted to the TSA database with the accession number SUB8976341.
Functional annotation of unigenes
Unigenes were annotated by blasting six common databases, including NCBI non-redundant protein sequences (NR), Swiss-Prot, Protein family (Pfam), Cluster of Orthologous Groups of proteins (COG), Gene Ontology (GO), and KEGG. A total of 26,461 (45.53 %) unigenes were functionally annotated. The number of unigenes matched to NCBI NR, COG, GO, Pfam, Swiss-Prot, and KEGG databases were 24,937 (42.90 %), 22,844 (39.30 %), 18,265 (31.43 %), 16,888 (29.06 %), 16,335 (28.10 %), and 14,110 (24.28 %), respectively (Supplement Figure 1 A).
Based on BLAST results against the NR database, most (16,697 or 66.9 %) of annotated S. frugiperda unigenes shared the highest similarity (first hit) to sequences from S. litura, followed by Helicoverpa armigera (1447 unigenes, 5.80 %), Trichoplusia ni (764 unigenes, 3.06 %), Heliothis virescens (720 unigenes, 2.89 %), Eumeta japonica (388 unigenes, 1.56 %), and C. suppressalis (347 unigenes, 1.39 %) (Supplement Figure 1B). Only 303 unigenes (1.22 %) showed the highest similarity to sequences from S. frugiperda, suggesting that this important pest has been understudied genomically.
The 16,888 annotated unigenes were divided into three categories: biological process, cellular component, and molecular function. The GO terms of binding, catalytic activity, and cellular process were with the most numbers of unigenes, with 9310, 8554, and 6308 in each category, respectively (Supplement Figure 1 C). The unigenes with KEGG annotations could be classified into five major categories, including Metabolism (3182 unigenes), Genetic Information Processing (2507 unigenes), Environmental Information Processing (1889 unigenes), Cellular Processes (1987 unigenes), and Organismal Systems (2663 unigenes) (Supplement Figure 1D). For the secondary categories, the pathways of signal transduction, translation, and carbohydrate metabolism were ranked as the top three subcategories, with 1681, 1200, and 1079 unigenes, respectively, in each subcategory.
Identification of DEGs based on transcriptomes
A total of 915 unigenes were expressed differentially between controls and samples treated with 1.0 µg/g CPT. Compared to the control group, 612 unigenes were up-regulated and 291 unigenes were down-regulated in the group treated with 1.0 µg/g CPT (Fig. 3 A). The number of DEGs between control and 5.0 µg/g CPT-treated samples increased to 3560. Among the DEGs, 2201 were up-regulated and 1359 down-regulated (Fig. 3 A). Comparative analyses revealed that 683 unigenes were differentially expressed in both 1.0 and 5.0 µg/g CPT-treated samples when compared to control. Among the common DEGs, 464 were up-regulated and 217 down-regulated (Fig. 3B and C).
A large number of DEGs were genes involved in detoxification, including genes encoding cytochrome P450 monooxygenases (P450s), glutathione S-transferases (GSTs), carboxylesterases (COEs), UDP glucosyltransferases (UGTs), and ATP-binding cassette transporters (ABCs). As shown in Table 1 and 39 detoxification genes were differentially expressed between controls and samples treated with 1.0 µg/g CPT. These differentially expressed genes encode 20 P450s, 3 GSTs, 6 COEs, 7 UGTs, and 3 ABCs. Most of these DEGs were up-regulated, including 15 coding for P450s, 5 for COEs, and 5 for UGTs (Table 1). The number of detoxification genes expressed differentially between controls and samples treated with 5.0 µg/g CPT increased to 108, including genes encoding 57 P450s, 13 GSTs, 8 COEs, 11 UGTs, and 19 ABCs. Among the up-regulated DEGs, 23 were genes coding for P450s, 1 for GST, 2 for COEs, 4 for UGTs, and 11 for ABCs (Table 1).
In addition to DEGs with functions in detoxification, several genes encoding mucins were also expressed differentially among control and treated samples. Mucins are high molecular weight glycoproteins covering the surface of epithelial cells that respond to external environmental stimuli such as infection, dehydration, and physical and chemical injury . Two and 18 genes encoding mucins were differentially expressed between controls and samples treated with 1.0 and 5.0 µg/g CPT, respectively. The unigenes encoding mucin-5AC (DN34507_c0_g1) and mucin-17 (DN3811_c0_g1) were up-regulated in samples treated with 1.0 µg/g CPT when compared to control. Most (16) DEGs encoding mucin proteins were up-regulated in samples treated with 5.0 µg/g CPT, and only two mucin genes were down-regulated (Fig. 4 A).
The third major group of DEGs included genes encoding cuticle proteins (CPs), which are indispensable structural components for insect tissues such as cuticle and midgut peritrophic membrane. Specifically, Four genes encoding larval cuticle protein LCP-17 (DN2495_c0_g1), cuticle protein 6.4-like (DN2113_c0_g1), cuticle protein CP14.6-like (DN38621_c0_g1) and cuticular protein RR-2 (DN1220_c0_g1) were down-regulated in samples treated with 1.0 µg/g CPT. Interestingly, 26 unigenes encoding cuticle proteins were up-regulated in samples treated with 5.0 µg/g CPT, whereas there were only two cuticle protein-encoding genes that were down-regulated (Fig. 4B).
GO and KEGG analyses
A total of 553 DEGs between controls and samples treated with 1.0 µg/g CPT were assigned to 175 GO terms. Among these GO terms, 30 were enriched significantly (corrected P-values < 0.05). The enriched GO terms for biological process included “carbohydrate derivative metabolic process”, “aminoglycan metabolic process”, “chitin metabolic process”, “glucosamine-containing compound metabolic process”, and “amino sugar metabolic process”. The enriched GO terms for cell component included “membrance part”, “integral component of membrane” and “intrinsic component of membrane”. The enriched GO terms for molecular function included “oxidoreductase activity”, “transporter activity”, and “transmembrane transporter activity” (Supplement Figure 2 A).
A total of 2000 DEGs between controls and samples treated with 5.0 µg/g CPT were assigned to 215 GO terms. Among these GO terms, 43 were enriched significantly (corrected P-values < 0.05). The most significantly enriched GO term for biological process was “chitin metabolic process” (corrected P-value = 7.27901E-08, 50 DEGs). The most significantly enriched GO term for cellular component was “extracellular region” (corrected P-value = 8.61215E-08, 107 DEGs). The most significantly enriched GO term for molecular function was “structural constituent of cuticle” (corrected P-value = 7.27901E-08, 35 DEGs) (Supplement Figure 2B).
KEGG analysis revealed that 369 DEGs between controls and samples treated with 1.0 µg/g CPT were assigned to 178 pathways. Among the 178 pathways, five were enriched significantly (corrected P-values < 0.05), including “DNA replication” (15 DEGs), “Purine metabolism” (15 DEGs), and “Ribosome biogenesis in eukaryotes” (14 DEGs). KEGG analysis assigned 1401 DEGs between controls and samples treated with 5.0 µg/g CPT to 232 pathways. The most significantly enriched pathways included “Purine metabolism” (42 DEGs), “Ribosome biogenesis in eukaryotes” (41 DEGs), and “Peroxisome” (36 DEGs) (Fig. 5 C and 5D). “DNA replication” was the most significantly enriched KEGG pathway in both samples treated with either 1.0 or 5.0 µg/g CPT.
To confirm the results of transcriptomic analyses, 20 unigenes including genes involved in detoxification and DNA replication, genes encoding mucins and cuticle proteins genes, were selected for qRT-PCR validation. As shown in Fig. 6, the expression patterns of the selected genes in S. frugiperda midguts changed significantly after CPT treatments based on qRT-PCR analysis. The changes in gene expression levels based on qRT-PCR were largely consistent with the transcriptomic data.
S. frugiperda has become a serious insect pest in China in the past couple of years . Various chemicals such as chlorantraniliprole, spinetoram, emamectin benzoate, spinetoram, acephate, and pyraquinil have been evaluated to control this pest in the field [27,28,29]. Some bioactive compounds including azadirachtin isolated from Azadirachta indica and celangulins extracted from the medicinal plant Celastrus angulatus have been studied for potential control of this destructive insect [11, 30]. CPT is a natural indole alkaloid used for cancer therapy . The insecticidal activity of CPT against other insect pests has also been investigated. In this study, development delay was induced in S. frugiperda larvae treated with CPT, but no mortality was observed. This result may be due to the high number of detoxifying enzyme genes that are often in polyphagous pests . In addition, synergism between CPT and Bacillus thuringiensis (Bt) or nucleopolyhedroviruses exists against Trichoplusia ni and S. exigua . Our results showed that CPT can inhibit the growth of S. frugiperda larvae. Therefore, CPT might be used as an independent insecticide for controlling S. frugiperda. Alternatively, CPT may be used with other insecticides for enhanced efficiency. One limitation of CPT as an insecticide is its insolubility in water. More efficient derivatives with improved solubility and hydrophobicity may be developed in the future for pest control.
The insect midgut is an important organ responsible for food digestion and nutrient absorption [34, 35]. CPT has been reported to induce alterations in the midguts of Trichoplusia ni and S. exigua larvae, including the loss of the single layer of epithelial cells and the disruption of the peritrophic membrane . In this study, we observed the loss of epithelial cells, abnormal cell structure, and intestinal wall degradation in the midgut of S. frugiperda after CPT treatments. These observations are consistent with previous findings in other insects, suggesting that CPT holds the potential as an insecticide for controlling S. frugiperda and other insect pests.
Recently, transcriptomic analysis has become a routine method to identify the differentially expressed genes in insects in response to toxic compounds . For example, transcriptomic analyses have been carried out to identify DEGs in the Chinese populations of S. frugiperda in response to 23 pesticides [37, 38]. DEGs in S. frugiperda larvae treated with azadirachtin were also initially analyzed . In this study, DEGs in S. frugiperda larvae treated with CPT were analyzed for the first time. Our transcriptomic analyses of midguts from S. frugiperda larvae revealed that the expression levels of a large number of genes were affected by CPT treatments.
Among the up-regulated genes by CPT are genes involved in detoxification. Metabolic detoxification through the overexpression of metabolic genes is considered one of the main ways to handle toxic insecticides by pests . The insect midgut is an important tissue responsible for pesticide detoxification where a variety of detoxification enzymes are produced . Insect midgut often increases the expression of metabolic genes in response to pesticide treatments. For example, the transcription levels of detoxifying-related genes including P450s and GSTs were up-regulated by low-dose of acetamiprid in the midgut of B. mori . Sublethal concentrations of Cry1Ca protein altered the expressions of P450s, CarEs, and GSTs in S. exigua larval midgut . Detoxification-related genes including those encoding P450s, GSTs, and COEs are up-regulated in S. litura larval midguts after being treated with tomatine . The roles of several detoxification genes in pesticide resistance in insects have been validated by RNA interference, including CYP321A8, CYP321A9, and CYP321B1 in S. frugiperda, a GST in Ostrinia furnacalis, Pxae18 and Pxae28 in Plutella xylostella [43,44,45]. Here we found that a large number of detoxification-related genes were coordinately up-regulated in the midgut of S. frugiperda larvae after CPT treatments. We speculated that up-regulation of these detoxification genes such as the genes encoding P450s, GSTs, and CarEs could accelerate the removal of CPT, thus reducing the toxicity of this chemical to S. frugiperda. Our results indicated that high doses of CPT induced more detoxification-related genes. The dose-dependent upregulation of these genes further indicates their roles in the detoxification of CPTs in S. frugiperda larvae. Assuming the up-regulated detoxification genes are indeed critical to detoxify CPT, application of CPT together with detoxification enzyme inhibitor(s) might increase the efficacy of CPT as a controlling agent for the S. frugiperda pest.
ABC transporters are membrane proteins that are divided into eight subfamilies (A to H) based on the conserved nucleotide-binding domain (NBD) and their functions . Three ABCG subfamily genes are up-regulated in H. armigera after being treated with abamectin, indoxacarb, and lambda-cyhalothrin . Five ABC transporters are induced in B. mori treated with NaF . In the present study, eleven ABC transporters were up-regulated in S. frugiperda larvae treated with 5.0 µg/g CPT. These up-regulated genes belong to either ABCC or ABCG subfamilies. The ABCC transporters in most insects are composed of full transporters (FTs) and half transporters (HTs), while ABCG transporters are typical HTs, which need to form homo- or heterodimers to perform transport functions . There is evidence to suggest that ABC transporters from the ABCC and ABCG subfamilies are involved in insecticide resistance [49, 50]. Therefore, our results may suggest that these up-regulated ABC transporters might be involved in CPT metabolism in S. frugiperda.
UGTs perform diverse functions including substrate detoxification, ecdysteroid metabolism, cuticle formation, pigmentation, and olfaction in insects. Lipophilic compounds could conjugate with UDP-glucose of UGTs and be glycosylated as water-soluble products that can be easily excreted . Multiple UGTs were found involved in xenobiotics detoxification and insecticide resistance. Several UGTs are up-regulated in nymphs of N. lugens exposed to β-asarone . Ten UGTs are differentially expressed in harmine-treated Sf9 cells . In this study, five and four UGTs were identified as DEGs in S. frugiperda larvae treated with 1.0 and 5.0 µg/g CPT, respectively. Two UGT genes were previously found involved in the detoxification of benzoxazinoids from maize in S. frugiperda larvae . These findings indicated that UGTs could play roles in CPT detoxification.
Mucins are highly O-glycosylated proteins that are abundant in salivary glands, the midgut, and malpighian tubules of insects . Mucins play an important role in lubricating the epithelial cells and protecting tissues from physical and chemical injuries . Mucins are also involved in host adaptation and oviposition regulation [56, 57]. In S. frugiperda, mucins are located in the peritrophic membrane and thus may protect epithelial cells and immobilize proteolytic enzymes . Multiple genes encoding different mucins were mostly up-regulated in samples treated with CPT, suggesting that mucins may have participated in repairing damages caused by CPT in the peritrophic membrane as reported previously [58, 59]. In addition to mucins, cuticle proteins are also primary components of the peritrophic membrane . Several genes encoding different cuticle proteins were also induced by CPT treatments. The up-regulation of genes encoding both mucins and cuticle proteins suggested that the peritrophic membrane in the midgut of S. frugiperda larvae may have gained increased resistance in this insect to CPT. If that is the case, application of CPT with other compounds that can prevent induction of mucin or cuticle protein genes or both may result in enhanced efficiency of CPT in controlling this insect pest. Furthermore, the expressions of DNA topoisomerases were not changed after CPT treatments. The likely reason is that CPT could affect the activities of topoisomerases, but does not affect the transcription level of topoisomerases.
Carbohydrates, lipids, proteins, nucleotides, chitin, and vitamins are important resources of energy for insect growth and development . Our GO analyses revealed that many metabolic processes were enriched with DEGs significantly, including chitin metabolic process, glucosamine-containing compound metabolic process, aminoglycan metabolic process, amino sugar metabolic process, and carbohydrate derivative metabolic process (corrected P-values < 0.05). Innate immune response and immune system processes were enriched (corrected P-values < 0.05) as well. KEGG analyses demonstrated that some pathways associated with Carbohydrate metabolism, Amino acid metabolism, Lipid metabolism, Nucleotide metabolism, Xenobiotics biodegradation, and Metabolism of cofactors and vitamins were enriched significantly, as well as the pathways of Digestive system, Endocrine system, Cell growth and death (corrected P-values < 0.05). These results suggest that CPT inhibited the growth of S. frugiperda larvae by disrupting many digestive and metabolic processes.
In conclusion, growth inhibition was observed in S. frugiperda larvae treated with different concentrations of CPT. Damage in the midgut was found after CPT treatments based on HE staining and TEM observation. Comparative transcriptomic analyses identified a large number of DEGs caused by CPT treatments. Genes involved in detoxification, DNA replication, and structural components of the peritrophic membrane were up-regulated significantly. Our results suggest that CPT exerted its impact on the growth of S. frugiperda larvae potentially by disrupting DNA replication, the digestive system, and other metabolic processes. Our studies provided a foundation for future research on the basic mechanism and potential practical application of CPT on S. frugiperda control.
A laboratory colony of S. frugiperda was generated from larvae collected from a cornfield in Conghua District, Guangzhou City, Guangdong Province, China. The colony has been maintained on an artificial diet since then. The artificial diet of 1 kg was prepared with 100 g cornflour, 80 g soy flour, 26 g yeast powder, 26 g agar, 8.0 g vitamin C, 2.0 g sorbic acid, 1.0 g choline chloride, 0.2 g inositol, 0.2 g cholesterol, and 900 mL distilled water. Adults were fed on 10 % honey water. Insect cultures were kept in an incubator set at 25 ± 1 °C, 60–70 % relative humidity, and a 16:8 h light: dark cycle.
CPT was purchased from Selleck Chemicals (USA) and was dissolved in dimethyl sulfoxide (DMSO). Third-instar larvae were selected and fed on artificial diets containing CPT at 1.0, 2.5, 5.0, 10, 20, and 30 µg/g, respectively, for 7 days. Twenty larvae were selected as one sample for the experiment and three biological replicates of each concentration were performed. Weights of insects for each sample were recorded after feeding for 1, 3, 5, and 7 days, respectively.
After seven days of feeding, the midgut of S. frugiperda larvae treated with 1.0 and 5.0 µg/g CPT was dissected and washed with cold phosphate-buffered saline (PBS). Midgut tissues were fixed with 4 % paraformaldehyde (#G1101, Servicebio, Wuhan, China) at 4 °C for more than 24 h. Midgut samples were then embedded into paraffin wax and sliced into sections of 4 μm thickness. The sections were fixed onto a glass slide and stained with hematoxylin and eosin solution. Histopathological changes in the midgut of larvae fed on CPT were visualized on a microscope (Nikon, Japan).
Transmission electron microscope observation
Larval midguts exposed to 1.0 and 5.0 µg/g CPT for 7 days were dissected, washed with cold PBS, and then fixed in glutaraldehyde solution at 4 °C for 4 h. After washing with PBS three times, the tissues were postfixed in 1 % osmic acid solution at room temperature for 2 h and then washed with PBS three times again. The tissues were dehydrated in a graded series of ethanol solutions (50–100 %) and then kept in 100 % acetone. The midguts were then embedded in epoxy resin and sliced into sections of 60–80 nm thickness. After staining with 2 % uranium acetate saturated alcohol solution and 2 % lead citrate solution for 15 min each time, the sections were dried at room temperature overnight. The ultrastructures of larval midgut were observed on an HT7800 TEM (HITACHI, Japan).
RNA isolation and Illumina sequencing
Ten to fifteen larval midgut samples exposed to 1.0 and 5.0 µg/g CPT for 7 days were collected, Total RNA was extracted using TRIzol® Reagent (Invitrogen, USA). Potential genomic DNA contaminants were removed using DNase I (TAKARA, Japan). The quality and concentration of RNA samples were analyzed on a 2100 Bioanalyser (Agilent, USA) and quantified on an ND-2000 (NanoDrop Technologies, USA). cDNA library construction and RNA-seq were conducted at Shanghai Majorbio Bio-pharm Biotechnology Co., Ltd. (Shanghai, China) via a commercial contract. Briefly, mRNA was purified from 5 µg total RNA of each sample using oligo-dT-attached magnetic beads. Purified mRNA was fragmented using a fragmentation buffer and the resulting products were used as templates for double-stranded cDNA synthesis with a SuperScript double-stranded cDNA synthesis kit (Invitrogen, CA). The double-stranded cDNA was subjected to end-repair, phosphorylation, and ‘A’ base addition. cDNA fragments of 200–300 bp were selected by 2 % Low Range Ultra Agarose and amplified in PCR reaction with Phusion DNA polymerase (New England Biolabs, Boston, MA) for 15 cycles. Then the products were quantified on TBS380 and sequenced on the Illumina HiSeq 4000 platform (Illumina, USA).
De novo assembly and annotation
After sequencing, the base quality, base error rate, and A/T/G/C base content distribution statistics of raw paired-end reads from each sample were evaluated with the software of fastx_toolkit_0.0.14 (http://hannonlab.cshl.edu/fastx_toolkit/). The reads that contained adaptor contamination, low-quality bases, and undetermined bases in raw reads were removed with SeqPrep (https://github.com/jstjohn/SeqPrep) and Sickle (https://github.com/najoshi/sickle) under the default parameters. The clean reads were assembled using Trinity (https://github.com/trinityrnaseq/trinityrnaseq) with the parameters of k-mer = 25 and min_kmer_cov = 5. Then the false assembled transcripts include chimera, structural errors, incomplete assembly, base errors, etc. were filtered with TransRate (http://hibberdlab.com/transrate/). Besides, the redundant sequences were identified and removed by CD-HIT (http://weizhongli-lab.org/cd-hit/). Subsequently, the assembly integrity of these transcriptomes was assessed by BUSCO (Benchmarking Universal Single-Copy Orthologs, http://busco.ezlab.org). The transcripts that shared sequence content were clustered and the longest transcript in the cluster was selected as the unigene. Assembled transcripts and unigenes were annotated by searching the NCBI protein non-redundant (NR), Swiss-Port, and Clusters of Orthologous Groups of proteins (COG) databases with DIAMOND v0.8.37.99 software using BLASTX with a threshold E-value of 1 × 10− 5. Besides, the Pfam database was also used for unigenes and transcripts annotation with the software of HMMER3 3.1b2. The function of each transcript was assigned according to the first hit. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed using BLAST2GO 2.5.0 (http://www.blast2go.com/b2ghome) and KOBAS 2.1.1 (KEGG, http://www.genome.jp/kegg/) with default parameters, respectively.
Differentially expressed genes (DEGs)
The method of fragments per kilobase of exon per million mapped reads (FRKM) was used to analyze the expression level of each transcript in different samples. Differentially expressed genes (DEGs) were identified using R statistical package software named EdgeR (Empirical Analysis of Digital Gene Expression in R, http://www.bioconductor.org/packages/2.12/bioc/html/edgeR.html). The criteria for DEGs were |log2Fold change| ≥1 and P-value ≤ 0.05. GO functional enrichment analysis of DEGs was performed via Goatools (https://github.com/tanghaibao/Goatools) with Fisher’s exact test. KEGG pathway enrichment analysis was conducted by KOBAS software (http://kobas.cbi.pku.edu.cn/home.do) with Fisher’s exact test and the pathway enriched with DEGs significantly was identified with the judging indicator of corrected P-value < 0.05.
Two µg total RNA of each sample was used for qRT-PCR cDNA synthesis using a FastKing RT Kit (With gDNase) (KR116, TIANGEN, Beijing, China) following the manufacturer’s instruction. Primers were designed using the Primer 5.0 software and synthesized by Thermo Scientific. Primer sequences were listed in Supplement Table 2. The TransStart® Tip Green qPCR SuperMix (AQ141, TransGen Biotech, Beijing, China) was used for qRT-PCR in a 10 µL reaction solution on a LighCycler480 II machine (Roche, Basel, Switzerland). qRT-PCR has proceeded as follows: one cycle of denaturation at 94 °C for 3 min, followed with 45 cycles of denaturation at 94 °C for 15 s, annealing at 60 °C for 20 s, and elongation at 72 °C for 20 s, followed by a melting curve analysis. Two reference genes RPL3 (AY072287.1) and RPL13 (AF400183.1) were selected for normalization of qRT-PCR results. mRNA levels were analyzed by the 2−△△CT method. Each assay was repeated three times.
Each experiment was conducted with three biological replicates and all data were expressed as the mean ± standard deviation (SD). One-way ANOVA was applied for statistical analysis of single larval weight. Homogeneity of variance test was determined by F test and different capital letters in figures indicate significant differences between different doses as determined using ANOVA followed by Duncan’s new multiple range method (DMRT) (p < 0.01) in SPSS 17.0 (IBM, USA).
Availability of data and materials
The raw reads of transcriptomes in this study have been deposited in the NCBI SRA database with the accession number of SRP242660.
Differentially expressed genes
Kyoto Encyclopedia of Genes and Genomes
- Topo I:
DNA topoisomerase I
Cluster of Orthologous Groups of proteins
NCBI protein nonredundant
Cytochrome P450 monooxygenases
ATP-binding cassette transporters
- Bt :
Fragments per kilobase of exon per million mapped reads
quantitative real time polymerase chain reaction
Ribosomal Protein L3
Duncan’s new multiple range method
Wu QL, He LM, Shen XJ, Jiang YY, Liu J, Hu G, et al. Estimation of the potential infestation area of newly-invaded fall armyworm Spodoptera frugiperda in the yangtze river valley of China. Insects. 2019;10(9):298.
Eghrari K, de Brito AH, Baldassi A, Balbuena TS, Fernandes OA, Môro GV. Homozygosis of Bt locus increases Bt protein expression and the control of Spodoptera frugiperda (Lepidoptera: Noctuidae) in maize hybrids. Crop Prot. 2019;124:104871.
Reall T, Kraus S, Goodman CL, Ringbauer J Jr, Geibel S, Stanley D. Next-generation cell lines established from the fall armyworm, Spodoptera frugiperda (Lepidoptera: Noctuidae). In Vitro Cell Dev Biol Anim. 2019;55(9):686–93.
Rivero-Borja M, Rodríguez-Maciel JC, Urzúa Gutiérrez JA, Silva-Aguayo G, Chandrasena DI, Felix-Bermudez NC, et al. Baseline of susceptibility to the Cry1F protein in mexican populations of fall armyworm. J Econ Entomol. 2020;113(1):390–8.
Cui Y, Ren YD, Lyu M, Zheng SC, Feng QL, Xiang H. Genomic divergences between the two polyphagous Spodoptera relatives provide cues for successful invasion of the fall armyworm. Insect Sci. 2020;27(6):1257–65.
Ordóñez-García M, Rios-Velasco C, Ornelas-Paz JD, Bustillos-Rodríguez JC, Acosta- Muñiz CH, Berlanga-Reyes DI, et al. Molecular and morphological characterization of multiple nucleopolyhedrovirus from Mexico and their insecticidal activity against Spodoptera frugiperda (Lepidoptera: Noctuidae). J Appl Entomol. 2019;00:1–10.
Morteza-Semnani K, Ghanbarimasir Z. A review on traditional uses, phytochemistry and pharmacological activities of the genus Ballota. J Ethnopharmacol. 2019;233:197–217.
Gao Q, Shi Y, Liao M, Xiao J, Li X, Zhou L, et al. Laboratory and field evaluation of the aphidicidal activity of moso bamboo (Phyllostachys pubescens) leaf extract and identification of the active components. Pest Manag Sci. 2019;75(12):3167–74.
Inocente EA, Nguyen B, Manwill PK, Benatrehina A, Kweka E, Wu S, et al. Insecticidal and antifeedant activities of malagasy medicinal plant (Cinnamosma sp.) extracts and drimane-type sesquiterpenes against Aedes aegypti mosquitoes. Insects. 2019;10(11):373.
Yang F, Wang L, Zhang L, Zhang Y, Mao L, Jiang H. Synthesis and biological activities of two camptothecin derivatives against Spodoptera exigua. Sci Rep. 2019;9(1):18067.
Lin S, Liu K, Wang R, Liu B, Zhang Y, Wu J, et al. Indoor toxicity of azadirachtin to Spodoptera frugiperda and its control effect in field. J South China Agri Univ. 2020;41(1):22–7.
Giongo AM, Vendramim JD, Freitas SD, Silva MF. Toxicity of secondary metabolites from meliaceae against Spodoptera frugiperda (J. E. Smith) (Lepidoptera: Noctuidae). Neotrop Entomol. 2016;45(6):725–33.
Silva TRFB, Almeida ACD, Moura TD, da Silva AR, Freitas SD, Jesus FG. Effect of the flavonoid rutin on the biology of Spodoptera frugiperda (Lepidoptera: Noctuidae). Acta Sci-Agron. 2016;38(2):165–70.
Dejane SA, Geraldo AC, Denilson FO, Angelita DC. Screening of Brazilian plant extracts as candidates for the control of Spodoptera frugiperda (Lepidoptera: Noctuidae). Rev Colomb Entomol. 2018;44(1):32–8.
Lucena DC, Bertholdo-Vargas LR, Silva WC, Machado AF, Lopes TS, Moura S, et al. Biological Activity of Piper aduncum extracts on Anticarsia gemmatalis (Hubner) (Lepidoptera: Erebidae) and Spodoptera frugiperda (JE Smith) (Lepidoptera: Noctuidae). An Acad Bras Cienc. 2017;89(3):1869–79.
Ayil-Gutiérrez BA, Villegas-Mendoza JM, Santes-Hernndez Z, Paz-González AD, Mireles-Martínez M, Rosas-García NM, et al. Ruta graveolens extracts and metabolites against Spodoptera frugiperda. Nat Prod Commun. 2015;10(11):1955–8.
Wang K, Tu Y, Wan JB, Chen M, He C. Synergistic anti-breast cancer effect of pulsatilla saponin D and camptothecin through interrupting autophagic–lysosomal function and promoting p62-mediated ubiquitinated protein aggregation. Carcinogenesis. 2020;41(6):804–16.
Sun Y, Zhang N, Wang C, Wei Y, Liu J. Distribution of camptothecin biosynthetic intermediates and identification the rate-limiting step of camptothecin biosynthesis. Nat Prod Res. 2019;19:1–8.
Yan L, Nan X, Zhang C, Wang H, Huang X, Hu J, et al. Development of an enzyme–linked immunosorbent assay for camptothecin. Mol Med Rep. 2019;20(2):959–66.
Liu YQ, Liu Y, Zhao YL, Li HY. Synthesis of novel derivatives of camptothecin as potential insecticides. Pestic Biochem Physiol. 2010;98:219–23.
Ma J, Tong S, Wang P, Liao W, Liu H, Zhang L. Insecticidal activity of camptothecin against Nilaparvata lugens, Brevicoryne brassicae, and Chilo suppressalis. J Econ Entomol. 2010;103(2):492–6.
Zhang ZJ, Shang XF, Yang L, Shi YB, Liu YQ, Li JC, et al. Engineering of peglayted camptothecin into nanomicelles and supramolecular hydrogels for pesticide combination control. Front Chem. 2020;7:922.
Tong SM, Feng MG. Laboratory and field evaluations of camptothecin sodium salt against phytophagous mites. Pest Manag Sci. 2016;72(3):629–36.
Gong L, Wang H, Huang J, Hu M, Hu Z, Zhong G. Camptothecin-induced expression of programmed cell death gene 11 in Spodoptera litura. Pest Manag Sci. 2014;70(4):603–9.
Dias RO, Cardoso C, Pimentel AC, Damasceno TF, Ferreira C, Terra WR. The roles of mucus-forming mucins, peritrophins and peritrophins with mucin domains in the insect midgut. Insect Mol Biol. 2018;27(1):46–60.
Wang R, Jiang C, Guo X, Chen D, You C, Zhang Y, et al. Potential distribution of Spodoptera frugiperda (J.E. Smith) in China and the major factors influencing distribution. Glob Ecol Conserv. 2020;21:e00865.
Wang YQ, Ma QL, Tan YT, Zheng Q, Yan WJ, Yang S, et al. The toxicity and field efficacy of chlorantraniliprole against Spodoptera frugiperda. J Environ Entomol. 2019;41(4):782–8.
Huang SQ, Zhang ZX, Jiang XY, Xu HH. The bioactivity of fused heterocyclic compounds of pyraquinil against Spodoptera frugiperda. Chinese J Appl Entomol. 2019;56(4):780–7.
Zhao S, Yang X, Yang X, Song Y, Wang W, Wu K. Field efficacy of eight insecticides on fall armyworm, Spodoptera frugiperda. Plant Prot. 2019;45(4):74–8.
Lu Y, Tian J, Zheng X, Xu H, Yang Y, Yang T, et al. Laboratory toxicity of 26 insecticides against different instar larvae of Spodoptera frugiperda. Acta Agri Zhejiangensis. 2019;31(12):2052–9.
Shu B, Zhang J, Zeng J, Cui G, Zhong G. Stability of selected reference genes in Sf9 cells treated with extrinsic apoptotic agents. Sci Rep. 2019;9:14147.
Rane RV, Ghodke AB, Hoffmann AA, Edwards OR, Walsh TK, Oakeshott JG. Detoxifying enzyme complements and host use phenotypes in 160 insect species. Curr Opin Insect Sci. 2019;31:131–8.
Sun S, Cheng Z, Fan J, Cheng X, Pang Y. The utility of camptothecin as a synergist of Bacillus thuringiensis var. kurstaki and nucleopolyhedroviruses against Trichoplusia ni and Spodoptera exigua. J Econ Entomol. 2012;105(4):1164–70.
Shu B, Zhang J, Cui G, Sun R, Yi X, Zhong G. Azadirachtin affects the growth of Spodoptera litura Fabricius by inducing apoptosis in larval midgut. Front Physiol. 2018;9:137.
Zhao T, Lai D, Zhou Y, Xu H, Zhang Z, Kuang S, et al. Azadirachtin A inhibits the growth and development of Bactrocera dorsalis larvae by releasing cathepsin in the midgut. Ecotoxicol Environ Saf. 2019;183:109512.
Shu B, Yu H, Li Y, Zhong H, Li X, Cao L, et al. Identification of azadirachtin responsive genes in Spodoptera frugiperda larvae based on RNA-sEq. Pestic Biochem Physiol. 2021;172:104745.
Gui F, Lan T, Zhao Y, Guo W, Dong Y, Fang D, et al. Genomic and transcriptomic analysis unveils population evolution and development of pesticide resistance in fall armyworm Spodoptera frugiperda. Protein Cell. 2020. doi:https://doi.org/10.1007/s13238-020-00795-7.
Hafeez M, Li X, Zhang Z, Huang J, Wang L, Zhang J, et al. De novo transcriptomic analyses revealed some detoxification genes and related pathways responsive to noposion Yihaogong® 5 % EC (lambda-cyhalothrin 5 %) exposure in Spodoptera frugiperda third-instar larvae. Insects. 2021;12(2):132.
Kliot A, Ghanim M. Fitness costs associated with insecticide resistance. Pest Manag Sci. 2012;68(11):1431–7.
Wang H, Lu Z, Li M, Fang Y, Qu J, Mao T, Chen J, Li F, Sun H, Li B. Responses of detoxification enzymes in the midgut of Bombyx mori after exposure to low-dose of acetamiprid. Chemosphere. 2020;251:126438.
Ren X, Wang Y, Ma Y, Jiang W, Ma X, Hu H, Wang D, Ma Y. Midgut de novo transcriptome analysis and gene expression profiling of Spodoptera exigua larvae exposed with sublethal concentrations of Cry1Ca protein. 3 Biotech. 2020;10(3):138.
Li Q, Sun Z, Shi Q, Wang R, Xu C, Wang H, et al. RNA-seq analyses of midgut and fat body tissues reveal the molecular mechanism underlying Spodoptera litura resistance to tomatine. Front Physiol. 2019;10:8.
Zhang BZ, Xu S, Zhen GA, Lu LY, Li YS, Ge X, et al. Silencing of cytochrome P450 in Spodoptera frugiperda (Lepidoptera: Noctuidae) by RNA interference enhances susceptibility to chlorantraniliprole. J Insect Sci. 2020;20(3):12.
Zhang Y, Zhang Y, Fu M, Yin G, Sayre RT, Pennerman KK, et al. RNA interference to control Asian corn borer using dsRNA from a novel glutathione-S-transferase gene of Ostrinia furnacalis (Lepidoptera: Crambidae). J Insect Sci. 2018;18(5):16.
Xie M, Ren NN, You YC, Chen WJ, Song QS, You MS. Molecular characterisation of two α-esterase genes involving chlorpyrifos detoxification in the diamondback moth, Plutella xylostella. Pest Manag Sci. 2017;73(6):1204–12.
Tang W, Xiao Y, Li G, Zheng X, Yin Y, Wang L, et al. Analysis of digital gene expression profiling in the gonad of male silkworms (Bombyx mori) under fluoride stress. Ecotoxicol Environ Saf. 2018;153:127–34.
Jin M, Liao C, Chakrabarty S, Zheng W, Wu K, Xiao Y. Transcriptional response of ATP-binding cassette (ABC) transporters to insecticides in the cotton bollworm, Helicoverpa armigera. Pestic Biochem Physiol. 2019;154:46–59.
Wu C, Chakrabarty S, Jin M, Liu K, Xiao Y. Insect ATP-binding cassette (ABC) transporters: roles in xenobiotic detoxification and Bt insecticidal activity. Int J Mol Sci. 2019;20(11):2829.
Meng X, Yang X, Wu Z, Shen Q, Miao L, Zheng Y, Qian K, Wang J. Identification and transcriptional response of ATP-binding cassette transporters to chlorantraniliprole in the rice striped stem borer, Chilo suppressalis. Pest Manag Sci. 2020;76(11):3626–35.
Dermauw W, Van Leeuwen T. The ABC gene family in arthropods: comparative genomics and role in insecticide transport and resistance. Insect Biochem Mol Biol. 2014;45:89–110.
Chen X, Xia J, Shang Q, Song D, Gao X. UDP-glucosyltransferases potentially contribute to imidacloprid resistance in Aphis gossypii glover based on transcriptomic and proteomic analyses. Pestic Biochem Physiol. 2019;159:98–106.
Xu X, Li X, Wang F, Han K, Liu Z, Fan L, et al. Candidate detoxification-related genes in brown planthopper, Nilaparvata lugens, in response to β-asarone based on transcriptomic analysis. Ecotoxicol Environ Saf. 2019;185:109735.
Cui G, Sun R, Veeran S, Shu B, Yuan H, Zhong G. Combined transcriptomic and proteomic analysis of harmine on Spodoptera frugiperda Sf9 cells to reveal the potential resistance mechanism. J Proteomics. 2020;211:103573.
Israni B, Wouters FC, Luck K, Seibel E, Ahn SJ, Paetz C, et al. The fall armyworm Spodoptera frugiperda utilizes specific UDP-glycosyltransferases to inactivate maize defensive benzoxazinoids. Front Physiol. 2020;11:604754.
Huang HJ, Liu CW, Xu HJ, Bao YY, Zhang CX. Mucin-like protein, a saliva component involved in brown planthopper virulence and host adaptation. J Insect Physiol. 2017;98:223–30.
Lou YH, Shen Y, Li DT, Huang HJ, Lu JB, Zhang CX. A mucin-like protein is essential for oviposition in Nilaparvata lugens. Front Physiol. 2019;10:551.
Terra WR, Dias RO, Oliveira PL, Ferreira C, Venancio TM. Transcriptomic analyses uncover emerging roles of mucins, lysosome/secretory addressing and detoxification pathways in insect midguts. Curr Opin Insect Sci. 2018;29:34–40.
Kelkenberg M, Odman-Naresh J, Muthukrishnan S, Merzendorfer H. Chitin is a necessary component to maintain the barrier function of the peritrophic matrix in the insect midgut. Insect Biochem Mol Biol. 2015;56:21–8.
Alvarenga ES, Mansur JF, Justi SA, Figueira-Mansur J, Dos Santos VM, Lopez SG, et al. Chitin is a component of the Rhodnius prolixus midgut. Insect Biochem Mol Biol. 2016;69:61–70.
Wang YW, Li YZ, Li GQ, Wan PJ, Li C. Identification of cuticular protein genes in the Colorado potato beetle Leptinotarsa decemlineata (Coleoptera: Chrysomelidae). J Econ Entomol. 2019;112(2):912–23.
This work was financially supported by the fund from Key-Area Research and Development Program of Guangdong Province (No. 2020B020223004), the Innovation Team Project in Guangdong Provincial Department of Education (2017KCXTD018), and Guangzhou Science and Technology Plan Projects (Grants 201805010008 and 201903010043), and National Training Program of Innovation and Entrepreneurship for Undergraduates (202011347001). The authors sincerely thank Prof. Ming-Shun Chen of Kansas State University for revising this manuscript.
This work was financially supported by the fund from Key-Area Research and Development Program of Guangdong Province (No. 2020B020223004), the Innovation Team Project in Guangdong Provincial Department of Education (2017KCXTD018), and Guangzhou Science and Technology Plan Projects (Grants 201805010008, 201903010043 and 201904010135) and National Training Program of Innovation and Entrepreneurship for Undergraduates (202011347001).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Shu, B., Zou, Y., Yu, H. et al. Growth inhibition of Spodoptera frugiperda larvae by camptothecin correlates with alteration of the structures and gene expression profiles of the midgut. BMC Genomics 22, 391 (2021). https://doi.org/10.1186/s12864-021-07726-8
- Spodoptera frugiperda
- Camptothecin Adverse effects
- Transcriptome analysis