- 1State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing, China
- 2State Key Laboratory of Cotton Biology, Institute of Cotton Research, Chinese Academy of Agricultural Sciences (ICR, CAAS), Anyang, China
- 3Insect Ecology, Institute of Plant Protection, Jilin Academy of Agricultural Sciences, Changchun, China
- 4Institute of Plant Protection, Guangdong Academy of Agricultural Sciences, Guangzhou, China
- 5School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China
Background: Biocontrol strategies are the best possible and eco-friendly solution to develop resistance against O furnacalis and improve the maize yield. However, the knowledge about underlying molecular mechanisms, metabolic shifts, and hormonal signaling is limited.
Methods: Here, we used an axenic and a consortium of entomopathogenic Beauveria bassiana OFDH1-5 and a pathogen-antagonistic Trichoderma asperellum GDFS1009 in maize and observed that consortium applications resulted in higher chlorophyll contents and antioxidants activities [superoxide dismutase (SOD), peroxidase (POD), proline, protease, and polyphenol oxidase (PPO)] with a decrease in O. furnacalis survival. We performed a comprehensive transcriptome and an untargeted metabolome profiling for the first time at a vegetative stage in fungal inoculated maize leaves at 0-, 12-, 24-, 48-, and 72-h post insect infestation.
Results: The consortium of B. bassiana and T. asperellum leads to 80–95% of O. furnacalis mortality. A total of 13,156 differentially expressed genes were used for weighted gene coexpression network analysis. We identified the six significant modules containing thirteen candidate genes [protein kinase (GRMZM2G025459), acyl-CoA dehydrogenase (GRMZM5G864319), thioredoxin gene (GRMZM2G091481), glutathione S-transferase (GRMZM2G116273), patatin-like phospholipase gene (GRMZM2G154523), cytochrome P450 (GRMZM2G139874), protease inhibitor (GRMZM2G004466), (AC233926.1_FG002), chitinase (GRMZM2G453805), defensin (GRMZM2G392863), peroxidase (GRMZM2G144153), GDSL- like lipase (AC212068.4_FG005), and Beta-glucosidase (GRMZM2G031660)], which are not previously reported that are highly correlated with Jasmonic acid - Ethylene (JA-ET) signaling pathway and antioxidants. We detected a total of 130 negative and 491 positive metabolomic features using a ultrahigh-performance liquid chromatography ion trap time-of-flight mass spectrometry (UHPLC-QTOF-MS). Intramodular significance and real time-quantitative polymerase chain reaction (RT-qPCR) expressions showed that these genes are the true candidate genes. Consortium treated maize had higher jasmonic acid (JA), salicylic acid (SA), and ethylene (ET) levels.
Conclusion: Our results provide insights into the genetics, biochemicals, and metabolic diversity and are useful for future biocontrol strategies against ACB attacks.
Introduction
Maize plants undergo several challenges during their growth, including pathogenic infections, environmental stresses, and insect attacks (Batool et al., 2020). Asian corn borer, Ostrinia furnacalis (Guenée), is among the most destructive insect pests of maize in Southeast and East Asia, especially China, causing an estimated yield loss of 6–9 million tons annually (He et al., 2003). Based on the importance of sustainable and eco-friendly agricultural practices (Guo et al., 2019), O. furnacalis control strategies have been shifted toward the application of biological control agents rather than harmful chemical pesticides (Mitchell et al., 2016; Batool et al., 2020). The Trichoderma spp. is traditionally viewed as a kind of biocontrol fungi against a range of plant pathogens (Dou et al., 2020), however, Trichoderma spp. has also been found with some functions against plant insects in different ways. More interestingly, in our previous study, we proved that T. asperellum and Beauveria bassiana are great entomopathogenic against O. furnacalis in maize (Batool et al., 2020). The plant herbivory response is a dynamic process that temporally and spatially regulates proteins, metabolic profiles, and transcript patterns (Ehlting et al., 2008; Barah and Bones, 2015). This suggests that it is essential to study the plant herbivory defense mechanisms at transcriptional and metabolic levels in response to insect-feeding and pre-inoculated entomopathogenic fungal strains (Li et al., 2016).
To combat herbivory, plants have evolved a broad range of defense mechanisms (Abebe, 2021). The morphological structures such as trichomes, waxy leaf surface, and numerous secondary metabolites act as the first barrier or the constitutive defense against herbivore attack. Plants have also developed an induced defense response, which affects the herbivore’s survival, growth, and development directly through the production of defensive proteins and secondary metabolites, and indirectly through the secretion of volatile compounds to attract the pests (Guo et al., 2019). Plant-induced defense is initiated by the perception of insect-feeding-derived elicitor, and it is a sophisticated mechanism to balance the growth and defense. When plants are exposed to a herbivore attack, different receptors and transporters could be activated and pass the signals through different signal molecules to activate the genes involved in the biosynthesis of phytohormones such as jasmonic acid (JA), ethylene (ET), and salicylic acid (SA), secondary metabolites, and defensive enzymes, which are toxic, repellent, interfere in the digestive system, survival, and nutrient absorption of insects (Qi et al., 2018). JA and SA accumulate and exhibit a major role in regulating the plant defense mechanisms against herbivory (Tzin et al., 2015; Guo et al., 2019). JA-induced signaling is primarily regulated upon insect feeding, whereas salicylic acid regulates the response to phloem-feeding insects (Tzin et al., 2015; Heidel-Fischer et al., 2018). Herbivory-induced defensive enzymes, including enzymatic [superoxide dismutase (SOD), peroxidase (POD), polyphenol oxidase (PPO), and protease] and non-enzymatic (proline) antioxidants in maize plants, help to develop resistance against a corn borer and are responsible for the regulation of plant defense mechanism through the production of secondary metabolites and reactive oxygen species (ROS) scavenging (Batool et al., 2020). Trichoderma is also capable of stimulating an induced systemic resistance (ISR) plant defense response against aphids, caterpillars, and nematodes, modulated by crosstalk between SA and JA (Di Lelio et al., 2021).
Entomopathogenic fungi have traditionally been assumed to help regulate insect populations. So far, they have aroused a limited interest as plant growth promoters, alongside of protecting plant form herbivores (González-Guzmán et al., 2020, 2021). However, some hypocrealean ascomycetes, such as B. bassiana, play other poorly understood ecological roles that might be useful in developing novel strategies for both increased crop production and crop protection (Sánchez-Rodríguez et al., 2018). Sánchez-Rodríguez et al. (2018) concluded in their study that B. bassiana successfully colonized the plant and boosted the spike production in wheat inoculated with seed-dressing and soil drenching methods and caused a 30–75% larval mortality after consuming the endophytically-colonized leaves. In their work, González-Guzmán et al. (2020) also examined the performance of two wheat varieties inoculated with entomopathogenic fungi including B. bassiana and Metarhizium brunneum by soil drenching and seed coating method and found an increased yield with no significant difference in plant heights. Their results support the seed-dressing method of inoculating with entomopathogenic fungi as a promising green technology for crop protection and enhancement with no negative effect on plant performance. Landa et al. (2013) identified the endophytic colonization of the Opium poppy plant by B. bassiana and concluded that it can protect the plant from stem gall wasp. The B. bassiana, M, brunneum, and Isaria farinosa were also found to increase the iron (Fe) availability in soil differently depending on particle size, increased plant height, and inflorescences the production of sunflower grown (Raya-Díaz et al., 2017).
Advances in omics and quantitative biology offer several ways to identify gene networks and their regulatory mechanisms in living systems. One such promising approach is the RNA-Sequence-based weighted gene co-expression network analysis (WGCNA) (Langfelder and Horvath, 2008). The WGCNA is a useful method to identify modules/networks of coexpressed genes. Furthermore, the correlation of these modules with phenotypic traits was useful to detect the key genes within the networks (Umer et al., 2020). However, no such large-scale study regarding the synergistic interaction of B. bassiana OFDH1-5 and T. asperellum GDFS1009 in triggering the omics-based regulatory mechanisms of defense genes, biochemical enzymes, and phytohormones induced by the O. furnacalis herbivory, are available.
Here, we hypothesized that the consortium of entomopathogenic B. bassiana OFDH1-5 and pathogenic antagonistic T. asperellum GDFS1009 fungal inoculants can trigger the enhanced defense response in maize, then the normal herbivory-induced response by an enhanced expression of defense-related genes, antioxidants production, phytohormone synthesis, digestive enzyme inhibitors, etc. So, we performed the physiological, biochemical, transcriptome, metabolome, and phytohormonal level analysis in field-grown vegetative stage maize inoculated with B. bassiana OFDH1-5 and T. asperellum GDFS1009, singly and in the consortium at 0-, 12-, 24-, 48-, and 72-h under O. furnacalis feeding. By integrating differentially expressed genes (DEGs) and biochemical antioxidants through weighted gene coexpression network analysis (WGCNA), we revealed gene networks and key candidate genes linked to JA, SA, and ethylene biosynthesis pathway and production of antioxidants and insect cuticle and digestive enzyme inhibiters. Our results highlight the positive impact of B. bassiana OFDH1-5 and T. asperellum GDFS1009 on increasing the herbivory-induced defenses and restricting the survival and growth of O. furnacalis.
Materials and Methods
Plant Growth and Treatments
Maize variety “Jingke 968” was planted in an open field located at the Gongzhuling Experimental Station of the Institute of Plant Protection, CAAS, Gongzhuling, China (43°30′ N, 124°47′ E; 224.9 m above sea level) on May 11, 2018, and May 06, 2019. The treatments include: (CK) untreated control, (BB-1) seed treatment with B. bassiana OFDH1-5, (BB-2) soil drenching with B. bassiana OFDH1-5, (TH-1) seed treatment with T. asperellum GDFS1009, (TH-2) soil drenching with T. asperellum GDFS1009, (BT-1) seed treatment with consortium of B. bassiana OFDH1-5 and T. asperellum GDFS1009, (BT-2) soil drenching with a consortium of B. bassiana OFDH1-5 + T. asperellum GDFS1009, and (IC) no inoculation. All treatments except CK were infested with O. furnacalis larvae. Each treatment was replicated three times.
The B. bassiana OFDH1-5 (ACCC32726) sourced through Jilin Academy of Agricultural Sciences, Gongzhuling, China and the T. asperellum GDFS1009 (CGMCC NO. 9512) obtained from the School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China, were selected for this study based on the previous results (Batool et al., 2020). Fungal strains were grown and conidial suspensions of B. bassiana OFDH1-5 (1 × 109), T. asperellum GDFS1009 (1 × 109), and combined suspension of B. bassiana OFDH1-5 + T. asperellum GDFS1009 (5.9 × 105 + 8.4 × 108) in 1:1 were prepared as described by Batool et al. (2020). An interaction analysis to check the synergistic effect of two strains was performed based on laboratory pathogenicity analysis data (Supplementary Table 1). Surface-sterilized maize seeds were soaked in respective conidial suspension and placed in a shaker incubator (HZP-250) at 25 ± 2°C and 200 rpm for 12 h (Motholo, 2019), and then dried under a laminar air flow hood before planting (Cherry et al., 2004). For insect control-IC, control-CK, and soil drenching treatments BB-2, TH-2, and BT-2, seeds were soaked in a sterilized.1% (v/v) tween 80. For soil drenching treatments, 30 ml/plant fungal suspension was applied around the root zone after 1 week of germination.
Fifty plants per treatment were planted and replicated three times with 5 m spacing between rows. Plants of each treatment were planted 2.5 m apart from each other to hinder communication. The field conditions and management were consistent with the common local farming agricultural practices and covered with nylon netting to prevent insect immigration. All the plants selected for this experiment were healthy and developmentally similar.
Confirmation of Endophytic Colonization of Plants
To confirm the endophytic colonization of plants by fungal stains, a re-isolation method was used. Leaves were surface-sterilized in 1% sodium hypochlorite for 3 min, 70% ethanol for 1 min followed by washing with sterilized distilled water three times. Leaves were then cut into small segments by using sterilized scissors and placed on potato dextrose agar (PDA) plates. After incubation of 3–7 days at 25 ± 2°C, the presence or absence of fungal growth from leaves was recorded and the colonization % was calculated (Quesada-Moraga et al., 2006).
Insect Infestation and Plant Sampling
The O. furnacalis neonates were obtained from the Institute of Plant Protection, CAAS, and reared on an artificial diet. Insect infestation for larval survival and maize damage rating was carried out at the vegetative and silking stage, representing 1st and 2nd generation infestation, respectively. Healthy plants were infested artificially with 60 neonates (<24 h old) on June 30 and June 25 for 1st generation infestation, and on August 10 and 6 for the 2nd generation infestation in 2018 and 2019, respectively. The infestation was done during the evening to avoid exposure of neonate to direct sunlight and high temperature.
For antioxidants, chlorophyll, and molecular assays, 3rd instar larvae of O. furnacalis were placed on each leaf of healthy vegetative stage plants covered with small net cages to restrict their movement. Damaged leaves around the larval feeding site were collected in three replicates at 0-, 12-, 24-, 48-, and 72-h post infestation and immediately frozen in liquid nitrogen. Samples were stored at −80°C until further analysis (Guo et al., 2019).
Maize Damage Rating and Larval Mortality
When fifth instar O. furnacalis larvae were observed around the maize field of 1st generation on July 28 and 23 and of 2nd generation on September 9 and 4 (ready to harvest), maize plants were dissected to record the larval mortality rate (%), plant height, tunnel number, tunnel length, kernel number kernel weight, and fungal outgrowth on larval cadavers (Xie et al., 2015).
Measurement of Plant Antioxidants and Chlorophyll Content
Enzymatic antioxidants including POD, SOD, protease, polyphenol oxidase, proline content, and chlorophyll content were estimated in maize leaves at 0-, 12-, 24-, 48-, and 72-h post infestation. Detailed protocols are described in our previous article Batool et al. (2020).
Transcriptome Profiling and Bioinformatics Analysis
The RNA was isolated using the TIANGEN kit (Beijing, China), following the manufacturer’s instructions. Sequencing libraries were generated using NEBNext Ultra™ RNA Library Prep Kit for Illumina (New England Biolabs, Ipswich, MA, United States) according to the manufacturer’s instructions, and index codes were added to attribute sequences of each sample. Library quality was assessed on the Agilent Bioanalyzer 2100 system. Clean reads were then mapped with the reference genome sequence by Hisat2 tools soft (Kim et al., 2015). Gene function was annotated based on database information. DEGs were identified using the DESeq2 with P-value < 0.05 (Wang et al., 2010). We used KOBAS software to test the statistical enrichment of DEGs in kyoto encyclopedia of genes and genomes (KEGG) pathways.
Weighted Gene Co-expression Network Analysis for Identification of Key Candidate Genes Involved in Maize Defense
Weighted gene coexpression network analysis (WGCNA) was performed in R (version 4.1.2) using default parameters to simplify genes into expressed modules to identify the defense-related hub genes (Langfelder and Horvath, 2008). The fragments per kilobase of transcript per million mapped reads (FPKM) values were normalized, and an adjacency matrix was constructed. The plant antioxidant and the chlorophyll analysis data, used as phenotypic data, were imported into the WGCNA package, and correlation between antioxidants and chlorophyll data and gene modules were calculated using the default settings. The WGCNA package was used to convert the adjacency matrix into a topological overlap matrix (TOM). After constructing a network, the transcripts with identical expression patterns were grouped into one module, and eigengenes were also calculated for these modules. The genes from each module were exported using the default parameters for cytoscape export.
Validation of Intramodular Candidates Through Real Time-Quantitative Polymerase Chain Reaction
Validation of gene expression by RT-qPCR was carried out in three biological replicates on Applied Biosystems 7500 Fast Real-Time PCR System (Applied Biosystems, Foster City, CA, United States) using an SYBR Green (TAKARA Bio Inc., Japan) following the manufacturer’s guide. The RNA was isolated using a TIANGEN kit (Beijing, China) and the complementary DNAs (cDNAs) were synthesized using the One-Step gDNA Removal and cDNA Synthesis SuperMix (TransGen Biotech Co., Ltd., Beijing, China) by using the user’s manual. The amplification program followed was 95°C (15 s), followed by 40 cycles at 60°C (60 s), and 95°C (30 s). Actin (accession number-EU585777.1) was used as a reference gene. Fold changes of gene expression level were calculated using the 2–ΔΔCT method (Guo et al., 2019). Primers used are given in Supplementary Table 4.
Metabolome Profiling and Quantification of Phytohormones
Metabolite extraction and analysis were performed by using ultra-high-performance liquid chromatography-mass spectroscopy (UHPLC-MS, ExionLC, AB SCIEX; Waters, Manchester, United Kingdom) by following the protocol described by Yu et al. (2020). Quantification of SA, JA, and ET was done by following Guo et al. (2021). All the samples were collected as three biological replicates.
Data Analysis
Data collected from the maize physiological and biochemical experiments were analyzed using Statistix 8.1 software and the significance of p < 0.05 was applied. Alignment of RNA-seq data was done with HISAT2 software. Assembling of transcripts with mapped reads was done by using String Tie. To quantify the expression levels of transcripts for each sample, we used Asprofile software. The WGCNA was performed in R packages (version 4.1.2). Statistical comparisons of metabolite concentration were made using a SAS statistics package version 9.2 (SAS Institute Inc, 2015) and a significance level of P < 0.05 was applied.
Results
Confirmation of Endophytic Colonization of Plants
By a re-isolation method, the mycelia growth of B. bassiana OFDH1-5 and T. asperellum GDSF1009 were observed growing from leaf segments placed in PDA plates after 3 days, whereas, in the control treatment, no fungal growth was observed. In comparison with the soil drenching treatments, 19.6 and 15.9%, 42.2 and 39.44%, and 55 and 43.17% higher growth was recorded in seed coated treatment BB-1, TH-1, and BT-1 in 2018 and 2019, respectively (Table 1 and Figure 1).
Table 1. Percentage (%) endophytic colonization of plants by fungal isolates among different treatments in 2018 and 2019.
Figure 1. Confirmation of endophytic colonization of plants by fungal isolates through the re-isolation method.
Larval Survival Rate
The B. bassiana OFDH1-5 and T. asperellum GDFS1009 significantly affected the mortality rate of O. furnacalis larvae in both generations. In both generations, 89.38 and 92.28% of larval mortality was observed in BT-1 treatment (a consortium of both fungi) of 2018 and 2019, respectively, followed by 73.03 and 69.68% mortality was observed in BB-1, which was higher than IC. Plants inoculated with soil drenching method (BB-2, TH-2, and BT-2) showed relatively low mortality (60–65%, 40–50%, and 70–80%) in comparison with plants inoculated with seed treatment method (BB-1, TH-1, and BT-1) in 2018 and 2019, respectively (Table 2).
Table 2. Differential effect of single and consortium of B. bassiana OFDH1-5 and T. asperellum GDFS1009 inoculation on larval survival rate, number, and length of tunnels under O. furnacalis attack for two generations in 2018 and 2019.
Maize Growth and Damage Rating
The consortium of B. bassiana OFDH1-5 + T. asperellum GDFS1009 has a positive effect on plant growth and damage rating. In Generation 1, the highest number of tunnels and tunnel lengths were recorded in the insect control treatment (IC), whereas this damage was reduced up to 86.3 and 94.6% in BT-1 in 2018, and 92.2 and 96.7% in 2019, respectively. The second lowest damage was recorded in BB-1 with 77.7 and 81.8% reduction in the number of tunnels and 91.3 and 93.8% reduction in tunnel lengths in 2018 and 2019, respectively. Similar results were recorded in generation 2 with the highest number of tunnels and tunnel lengths in IC and lowest (82.9 and 93% in 2018, and 84.9 and 95.6% in 2019) in BT-1 (Table 2).
Over the 2-year average, plant heights were reduced up to 25.65% in the insect control treatment (IC) compared to the untreated control, while application of B. bassiana OFDH1-5 and T. asperellum GDFS1009 reduced the negative effect of O. furnacalis feeding and increased the plant heights up to 30.65, 29.1, and 34.5% in BB-1, TH-1, and BT-1, respectively. Less effect on plant height was seen in plants inoculated with the soil drenching method (Table 3). Additionally, the application of entomopathogenic fungi significantly increased the grain yield. The maximum increase in the average number of kernels per ear was recorded as 55.2 and 70.9% in BT-1, followed by 53.4 and 69.1% in BB-1, and 47.9 and 56.6% in TH-1 in 2018 and 2019, respectively. Similar patterns of increase in weight per 100 kernels were observed with 53.3 and 79.8% in BT-1, 46.2 and 65.7% in BB-1, and 36.1 and 43.4% in TH-1 in 2018 and 2019, respectively. In contrast, the maximum reduction in grain yield was recorded in insect control (IC) treatment with 20.6 and 29.1% in the number of kernels per ear, and 21.1 and 32.1% in weight per 100 kernels in the year 2018 and 2019, respectively, compared to untreated control (CK) (Table 3).
Table 3. Differential effect of single and consortium of B. bassiana OFDH1-5 and T. asperellum GDFS1009 inoculation on maize yield under O. furnacalis attack in 2018 and 2019.
Larval cadavers were collected from the field to determine the percentage (%) of cadavers showing a fungal outgrowth collected from the field. The highest number of cadavers (%) with fungal outgrowth was observed as 87.33 and 91.83% in seed treated consortium treatment (BT-1) in 2018 and 2019, respectively, at Generation 1 infestation. Similarly, in Generation 2 infestation, 89.68 and 84.07% cadavers were detected in BT-1 with fungal outgrowth in 2018 and 2019, respectively, followed by 71–74% in seed treated with B. bassiana treatment (BB-1) (Table 4 and Figure 2). Based on agronomic field data, samples from BB-1, TH-1, BT-1 (inoculated through seed treatment method), and IC were used for further analysis.
Figure 2. Field collected dead larvae showing the existence of fungal outgrowth from the larval cadavers.
Transcriptome Responses Triggered by B. bassiana OFDH1-5 and T. asperellum GDFS1009
For transcriptome profiling, four treatments (BB-1, TH-1, BT-1, and IC) in three replicates at 0-, 12-, 24-, 48-, and 72-h post O. furnacalis infestation were used to construct 60 cDNA libraries. A total of 458.04 Gb clean data was obtained and the clean data of each sample reached 5.71 Gb, and the Q30 base percentage was 93.24% and above. Clean reads were aligned with a designated reference genome (ZmB73_Ref-Gen_v4) and the comparison efficiency ranged from 79.13 to 89.36% (Supplementary Table 2). Based on the comparison results, the gene expression level analysis was performed and the differentially expressed genes were identified according to their expression levels | log2 (fold-change)| > 1 and an adjusted P-value < 0.05 in each pairwise comparison. The abundance of upregulated genes was higher than downregulated genes compared to all-time points (Supplementary Table 3). Identified genes in all samples were hierarchically clustered relative to control. The difference in color indicates high (red) and low (green) expressions (Figure 3A). The principal component analysis revealed the variability among transcriptome data of different samples. Clustering of samples away from each other and control group (IC) indicates a fungal inoculation and O. furnacalis-induced changes in gene expression specially in consortium treatment (BT) (Figure 3B). In total, 13,156 common DEG’s were identified in all the combinations at different time points.
Figure 3. Transcriptome response of maize inoculated with single and consortium of Beauveria bassiana OFDH1-5 and Trichoderma asperellum GDFS1009 to Ostrinia furnacalis. (A) Hierarchical Cluster dendrogram of differentially expressed genes, the difference in color indicate high (red) and low (green) expression. (B) PCA plot showing variability of transcriptome data; (C) top enriched resistance related pathways in IC vs. BB, IC vs. TH, IC vs. BT at 12-, 24-, 48-, and 72-h.
To explore the response of different processes in combined B. bassiana OFDH1-5 and T. asperellum GDFS1009 treatments, we analyzed the KEGG enrichment analysis of common DEG’s. The sets of DEG’s were assigned to significant KEGG Pathways (P > 0.05) (Supplementary Figure 1), and the top enriched resistance related pathways in IC vs. BB, IC vs. TH, and IC vs. BT at 12-, 24-, 48-, and 72-h were listed in Figure 3C, which show uniquely induced genes enriched in a higher number of metabolic processes in consortium treatment (BT).
B. bassiana OFDH1-5 and T. asperellum GDFS1009 Triggered Antioxidants and Chlorophyll Content in O. furnacalis-Induced Defense
The antioxidant enzyme activities started to increase at 12-h in all the infested treatments and reached a maximum point at 24-h of insect feeding and gradually decreased until 48-h, while polyphenol oxidase (PPO) reached the maximum at 48-h. The production of antioxidant enzymes was further triggered by entomopathogenic fungal treatments and the highest productions were observed in consortium treatment (BT-1). An increase of 132. 5-, 385. 3-, 100. 6-, and 277.7-folds were recorded in SOD, POD, proline, and protease activities at 24-h of O. furnacalis feeding. In the case of PPO, a 522.7-fold increase was recorded at 48-h and the lowest enzyme activities were observed in IC treatment (Figures 4A–E).
Figure 4. Beauveria bassiana OFDH1-5 and T. asperellum GDFS1009 induced response of maize antioxidants and chlorophyll content under O. furnacalis attack. (A) Superoxide dismutase (SOD); (B) peroxidase (POD); (C) polyphenol oxidase (PPO); (D) protease; (E) proline; (F) chlorophyll a; (G) chlorophyll b; (H) chlorophyll a + b; (I) carotenoids. Different lowercase letters above each bar indicate significant differences among treatments (p < 0.05). Treatments include (BB-1) seed treatment with B. bassiana OFDH1-5, (TH-1) seed treatment with T. asperellum GDFS1009, (BT-1) seed treatment with a consortium of B. bassiana OFDH1-5 and T. asperellum GDFS1009, (IC) no inoculation.
Chlorophyll content was observed to be gradually decreased from 0-h to 72-h due to O. furnacalis feeding, but fungal inoculation increased the chlorophyll content. The maximum content of chlorophyll-a and chlorophyll-b were recorded in BT-1, followed by TH-1, and BB-1 compared to insect control (IC). Chlorophyll-a, chlorophyll-b, and carotenoid were reduced up to 80-, 86-, and 86.3-folds, respectively, in the insect control (IC) (Figures 4F–I).
Co-expression Network Analysis for the Identification of Hub Genes Involved in Maize Defense
The FPKM values of 13,156 common DEGs and plant biochemical parameters (SOD, POD, proline, protease, PPO, Chl-a, chl-b, carotenoids, and chl-a + b) at different time points were used for the weighted gene coexpression network analysis and the module trait correlations were calculated. Thirty-three gene modules were identified based on coexpression patterns. Each module is represented by a different color and presented as a cluster dendrogram and network heatmap (Figures 5A–D).
Among the 33 identified modules, only six showed significant correlations with the phenotypic data. The blue module with 326 genes showed significant correlations with SOD, proline, protease, and PPO with correlation coefficients (r2) of 0.93, 0.82, 0.7, and 0.75, respectively. The darkgrey module possessing 149 genes showed positive correlations with POD (r2 = 0.84), PPO (r2 = 0.73), and protease (r2 = 0.95). The green-yellow module, having 223 genes, showed significant positive correlations with POD (r2 = 0.71), SOD (r2 = 0.81), proline (r2 = 0.73), protease (r2 = 0.71), and PPO (r2 = 0.9). Grey60 module consists of 57 genes, and it showed a positive correlation with carotenoids (r2 = 0.73). Magenta module with 98 genes showed positive correlations with POD (r2 = 0.82) and with PPO (r2 = 0.93). Salmon module with 145 genes showed a correlation with proline (r2 = 0.85) (Figure 5D). The hub genes from these modules were selected by the cytoscape built-in extension, namely “CytoHubba,” visualizing the gene networks (Figure 5E). Furthermore, to find out the genes involved in maize defense within these gene networks, gene annotation information was extracted from maize reference genome (ZmB73_Ref-Gen_v4).
Figure 5. Weighted gene coexpression network analysis (WGCNA) analysis to identify key candidate genes involved in maize defense. (A) Cluster dendrogram and network heatmap of genes subjected to coexpression module calculation. (B) Sample dendrogram and module trait heatmap of each time point. (C) Hierarchical clustering presenting eleven modules having coexpressed genes. Each leaflet in the tree corresponds to an individual gene. (D) Module-trait associations based on Pearson correlations. The color key from green to red represents r2 values ranging from –1 to 1. (E) Gene networks of six highly correlated modules with phenotypic traits. Genes in different red shape due to the highest weight within the module represents hub genes (key candidate genes) and are coded for gene descriptors based on annotation. Treatments include (BB-1) seed treatment with B. bassiana OFDH1-5, (TH-1) seed treatment with T. asperellum GDFS1009, (BT-1) seed treatment with a consortium of B. bassiana OFDH1-5 and T. asperellum GDFS1009, (IC) no inoculation.
Real Time-Quantitative Polymerase Chain Reaction Validation of Intramodular Hub Genes
Selected hub genes were further validated by performing RT-qPCR at different time points, thus, narrowing down the number of selected genes. Finally, from 28 hub genes, we selected 13 key candidate genes with stable and consistent expressions responsible for maize defense against O. furnacalis (Figure 6). The RT-qPCR results were consistent with the RNA-seq data (Supplementary Figure 2). All the key candidate genes, along with other hub genes, are shown in Table 5.
Figure 6. Schematic representation of transcriptional strategies employed by maize inoculated with a consortium of B. bassiana OFDH1-5 and T. asperellum GDFS1009 compared with single inoculation to induce enhanced expression of key genes to involve in initiating signaling pathway to activate phytohormone biosynthesis pathways which regulate the production and expression of defense genes to induce maize resistance against O. furnacalis. Treatments include (BB-1) seed treatment with B. bassiana OFDH1-5, (TH-1) seed treatment with T. asperellum GDFS1009, (BT-1) seed treatment with a consortium of B. bassiana OFDH1-5 and T. asperellum GDFS1009, (IC) no inoculation.
B. bassiana OFDH1-5 and T. asperellum GDFS1009-Enhanced Expression of Key Candidate Genes in Maize Defense Response Against O. furnacalis Herbivory
Three key genes in the blue module were identified as protein kinase gene (GRMZM2G025459), which participates in the phytohormone signaling pathway to initiate a defense response and to regulate the expression of genes involved in ion homeostasis and oxidative stress responses, acyl-CoA dehydrogenase (GRMZM5G864319), which is involved in an alpha-linolenic acid metabolism of jasmonic acid pathway and thioredoxin gene (GRMZM2G091481), which acted as an important component of the signaling pathway in plant antioxidant networks. Similarly, in the salmon module, the Glutathione S-transferase was identified as a key gene (GRMZM2G116273), which plays a crucial role in plant oxidative stress response against high ROS (reactive oxygen species) levels generated by the insect attack.
Four genes, (GRMZM2G144153, AC212068.4_FG005, GRMZM2G392863, and GRMZM2G004466) were identified as key genes in darkgrey module. These genes were identified as glutathione peroxidase, GDSL- like lipase, defensin, and protease inhibitor. Glutathione peroxidase gene (GRMZM2G144153) is known to be activated upon a cell wall damage by insect feeding. In addition, it is also involved in the ROS scavenging, lignification of the cell wall, and synthesis of phenolic compounds. The GDSL- like lipase (AC212068.4_FG005) acts as an elicitor of systemic resistance linked with ET pathways. Upon insect feeding, its expression depending upon the ethylene signaling generates the systemic signals which translocate to other healthy parts of plants and propagate a systemic resistance throughout the plant. The defensin gene (GRMZM2G392863) has a great insecticidal property of inhibiting insect digestive enzymes; α-amylase play a crucial role in breaking down plant starch and hydrolyze proteins. Protease inhibitor genes (GRMZM2G004466) were characterized as key genes. One more protease inhibitor gene (AC233926.1_FG002) was identified in the grey60 module as a key gene. They were found to participate in the reduction of insect growth, and, thereby, suppress the infection process by inhibiting extracellular proteases released by insects.
Patatin-like phospholipase gene (GRMZM2G154523) is responsible for the synthesis of linolenic acid, as a precursor of jasmonic acid biosynthesis pathway and chitinase gene (GRMZM2G453805), expressed by plants and acted as a self-defense mechanism against insect attack, were found as defense-related hub genes in the green-yellow module.
Beta-glucosidase 18 and cytochrome p450 (GRMZM2G031660 and GRMZM2G139874), which participate in various detoxification and biosynthetic pathways to protect the plant, were found as key genes in the magenta module. All genes showed differential expression patterns among all fungal-inoculated treatments, with the highest expression levels at 24- and 48-h post larval infestation specially in BT-1. The role and expression patterns of all key genes in maize defense response are shown in Table 5, Figure 6, and Supplementary Figure 2.
Metabolomic Adjustments Triggered by B. bassiana OFDH1-5 and T. asperellum GDFS1009
To obtain complete insights of B. bassiana OFDH1-5 and T. asperellum GDFS1009 inoculated maize under the O. furnacalis attack, we performed a comprehensive untargeted metabolomic analysis of all BB-1, TH-1, BT-1, and IC at 0-, 12-, 24-, 48-, and 72-h by UHPLC-QTOF-MS for the first time. The principal component analysis (PCA) analysis clearly showed the diversity among samples. Interestingly 0-h samples grouped indicating less diversity and other samples clustering away from control showed diversity in metabolomic data (Figure 7A). Through electrospray ionization (ESI) total of 7,367 NEG (negative) and 13,183 POS (positive) mass/retention time features were detected (Supplementary Data 4), and further 130 NEG and 491 POS differentially annotated and abundant metabolomic features were filtered based on different classes (Supplementary Data 5 and Figure 7B). Within the POS metabolic features, 9 fall in alkaloids and derivatives class, 26 in benzenoids, 3 in hydrocarbons, 194 in lipids and lipid-like molecule, 17 in nucleosides, nucleotides, and analogs, 47 in organic acids and derivatives, 11 in organic nitrogen compounds, 37 in organic oxygen compounds, 78 in organoheterocyclic compounds, 8 in organooxygen compounds, 57 in phenylpropanoids and polyketides and 1 each in hydrocarbon derivatives, organic compounds, organonitrogen compounds and lignans, neolignans, and related compounds. In the case of NEG metabolic features, 10 fall in benzenoids class, 15 each in nucleosides, nucleotide and analogs class and organic acid and derivative class, 1 in organic compounds, 16 in organic oxygen compounds, 13 in organoheterocyclic compounds, 6 in organooxygen compounds, 13 in phenylpropanoids and polyketides, and 41 in lipids and lipid-like molecules class (Figure 7B). Through hierarchical clustering, metabolites were divided into five clusters. Obvious changes were seen in the abundance of all metabolites belonging to different classes with a higher abundance of malonic acid, succinic acid, threonic acid, and gluconolactone in cluster 1, silicristin and mulberrin in cluster 2, etrogol and jasmonic acid in cluster 3, and betaine, L-threonine, L-carnitine, and thermophilin in cluster 5 were observed in a consortium treatment (BT-1). In cluster 4, all inoculated treatments have shown higher abundance as compared to control (IC) (Figure 7C).
Figure 7. Metabolome profiling of maize inoculated with single and consortium of B. bassiana OFDH1-5 and T. asperellum GDFS1009 and subjected to O. furnacalis attack. (A) PCA plot show divergence in metabolome data among treatments (MVBB, MVTH, MVBT, and MVI stands for BB-1, TH-1, BT-1, and IC, respectively). (B) Pie chart showing the percentage distribution of differentially annotated metabolites in different classes. (C) Hierarchical clustering of identified abundant metabolites. (D) Concentration of phytohormones in maize plants. Data is mean of three replicates and different lower-case letters above each bar indicate significant differences among treatments (p < 0.05). Treatments include (BB-1) seed treatment with B. bassiana OFDH1-5, (TH-1) seed treatment with T. asperellum GDFS1009, (BT-1) seed treatment with a consortium of B. bassiana OFDH1-5 and T. asperellum GDFS1009, (IC) no inoculation.
B.bassiana OFDH1-5 and T. asperellum GDFS1009-Enhanced Phytohormone Concentration Induced by O. furnacalis Herbivory
To determine changes in stress-related phytohormones, we quantified JA, ET, and SA. A significant increase in the contents of phytohormones was recorded in all fungal inoculated treatments, but highest in consortium treatment (BT-1) at 24-h post infestation with an increase of 32.2, 37.4, and 77.5% in JA, SA, and ET, respectively (Figure 7D).
Discussion
Ostrinia furnacalis is one of the most destructive insect pests of maize and it is responsible for major yield losses. It causes major damage at whorl stage plants by feeding in leaves and boring in maize stalk. Besides other pest control strategies, the biological control approach is of the highest interest because of its eco-friendly nature. The latest research proved that a consortium of biological control agents was effective against maize pests (Sarma et al., 2015; Batool et al., 2020). In our previous study (Batool et al., 2020), we have investigated and concluded that the entomopathogenic strain, B. bassiana, and plant pathogen antagonistic strain, T. asperellum, have both strong synergistic potential and can be used in the consortium to enhance their biological control activity. They inhibit the insect survival and damage on maize plants by inhibiting the immune response of O. furnacalis at the transcriptome level. Based on our previous findings, in the present investigation, we further investigated the effect of B. bassiana, T. asperellum, and their consortium on insect survival, maize damage, and growth, at field level under natural environmental conditions. Furthermore, WGCNA was used for the first time in an entomopathogenically inoculated maize to integrate transcriptome and biochemical data and an untargeted metabolite profiling at the vegetative stage of maize to identify the candidate genes and metabolites associated with the regulation of defensive biochemical enzymes, and induction of maize resistance and maize defense pathways against O. furnacalis.
Both entomopathogenic/antagonistic fungi can colonize the plant endophytically and can induce a defensive response upon insect feeding as a challenge, thereby killing or arresting the growth of insects (Batool et al., 2020). Plants can provide a suitable environment for entomopathogenic fungi and their insecticidal activity can suppress the infection and facilitate plants by increasing insect resistance and plant growth (Gupta et al., 2016). Plant colonization largely depends upon the method of application and the seed treatment method was proven to be more effective (Batool et al., 2020). The seed treatment has the advantage of colonizing radical and plumule at the same time, as they are close to seed and can colonize the whole plant during its growth (Muvea et al., 2014). Our results were also in accordance with the findings of pathogenicity against O. furnacalis and reducing its survival, tunnel number and length, and increasing grain yield and plant growth in the consortium-treated Plants (BT-1). Some previous research also showed that the growth of entomopathogenic fungi in plants can prevent plant herbivory, and the B. bassiana have insecticidal activity against several insects including corn borers (Wagner and Lewis, 2000). Trichoderma is generally known for the function to induce the ISR of the whole maize plant against plant diseases and pest insects through colonization of maize roots (Rivas-Franco et al., 2020). The activity of T. asperellum was relatively low as compared to others, but its pathogenic activity was enhanced by the co-cultivation with B. bassiana. Shakeri and Foster (2007) also stated that the Trichoderma pathogenicity against borer larvae was lower and it can be increased by modifying them at the genetic level. Trichoderma could generate synergistic action with other biocontrol organisms, including the entomopathogenic fungi (Netzker et al., 2015), describing the use of consortium technique to enhance the activity through intercommunication.
In response to an insect attack, plants produce several defensive antioxidants including, SOD, POD, protease, PPO, and proline to scavenge ROS. These enzymes were considered the most important plant defense-related enzymes against herbivores (Guo et al., 2017; Batool et al., 2020). We also observed that antioxidant activities such as SOD, POD, PPO, protease, and proline were increased in all insect-infested plants, with a maximum increase at 24-h post infestation, while these enzyme activities were relatively higher in entomopathogenic fungal inoculated plants. The results revealed that these biocontrol agents enhanced the defense mechanism in plants by enhancing the production of defensive enzymes. Consortium treatment showed a maximum increase among all treatments. Some previous studies also reported the induction of SOD, POD, proline, protease, and PPO by biological control agents (Bano and Muqarab, 2017). PPO and protease play a defensive role against insect attack and induce resistance in plants (Joe and Muthukumaran, 2008). The chlorophyll content is one of the important components for plant growth and we observed that the chlorophyll content was increased in plants treated with entomopathogenic fungi compared to the insect control (IC). The same trend was recorded in previous studies (Bano and Muqarab, 2017; Batool et al., 2020). SOD is considered as the first enzyme in the scavenging process of ROS and its augmented scavenging action and detoxified by POD and act as a signaling agent to activate defensive genes in plants (Bano and Muqarab, 2017). Proline is best known as hydroxyl radical scavenger also serves as an energy source in plants. PPO is an important enzyme, which induces insect resistance in plants by participating in cell lignification and oxidation of polyphenols. Protease also has a defensive role against insect attacks (Batool et al., 2020).
With rapid development in bioinformatics, the weighted gene co-expression analysis (WGCNA) approach has been developed to explore the functionality of transcriptome. WGCNA was used to detect the genes which co-express (Lin et al., 2019). The relationship between trait and module identifies the most vital hub genes, which supports an efficacious way to explore the in-depth mechanism of complex traits (Lv et al., 2020). Considering the importance of fungal inoculants against the control of O. furnacalis, plant yield, and maize defense through augmented production of enzymatic and non-enzymatic antioxidants. We also integrated these biochemical enzymes and compounds with plant transcriptome at 0-, 12-, 24-, 48-, and 72-h. We performed a transcriptome-based analysis to identify DEGs in different combinations. The DEGs were then analyzed by WGCNA along with biochemical parameters and 33 gene modules represented by different colors were identified. Among them, 6 were significantly correlated with SOD, POD, PPO, protease, proline, and chlorophyll content. All genes contained in modules were functionally annotated to identify hub genes and 28 hub genes were identified and validated through RT-qPCR. Finally, 13 genes were characterized as the key candidates that are involved in maize defense mechanism against O. furnacalis based on significant correlation of hub genes from the module with biochemical enzymes and chlorophyll content, and expression pattern of selected genes. Gamboa-Tuz et al. (2018) also identified stress related hub genes from papaya leaves and roots. Pan et al. (2020) correlated transcriptome through WGCNA to identify key genes in maize against herbivore.
Activity of Identified Candidate Genes in Phytohormone Signaling and Biosynthesis Pathways
Insect herbivores induce the expressions of genes associated with the phytohormone signaling pathway (Li et al., 2016). Plant recognizes herbivore damage as a primary signal and triggers the expression of defense related genes. Initially, when an insect starts feeding, the plant reacts to it by producing ROS and effector proteins to activate the expression of signaling molecules and initiate signaling pathways of systemic resistance (Almagro et al., 2009). We identified four such signaling related genes; protein kinase (GRMZM2G025459), GDSL-like Lipase (AC212068.4_FG005), Patatin-like phospholipase (GRMZM2G154523), and Cytochrome P450 (GRMZM2G139874), which showed a high response against O. furnacalis feeding in fungal inoculated treatments. The activity of GDSL- like lipase (GLIP) is highly dependent on the Ethylene ET- the pathway to generate systemic signals to produce GLIP in the uninfected part of the plants to protect it from the insect (Kwon et al., 2009). They also reported the reduction in the growth of Pseudomonas syringae by the accumulation of GLIP 1 in all parts of the plant. Patatin like phospholipase is a signaling molecule responsible for the release of linolenic acid from phospholipids to biosynthesize jasmonic acid (JA) (Ryu, 2004). Meijer and Munnik (2003) observed the presence of phospholipases in wounded tomatoes and many other plants. Dhondt et al. (2002) also found three lipase genes in response to external elicitor or infection before the accumulation of JA. Acyle co-A dehydrogenase gene (GRMZM5G864319) was also identified as a candidate gene that participates in the linolenic acid metabolism of the JA biosynthesis pathway (Tang et al., 2015). Cytochrome P450 (CYP) participates in plant growth promotion and protection from biotic and abiotic stresses through detoxification and biosynthesis pathways. To date, the involvement of cytochrome p450 in different metabolic pathways including allene oxide synthesis in the jasmonic acid pathway, hormone metabolism, and biosynthesis of phytoalexin has been confirmed (Jun et al., 2015). CYP74A (allene oxide synthase, AOS) catalyzes the dehydration of hydroperoxide to an instable allene oxide in the JA biosynthetic pathway and plays an important role in wound-induced defense against biotic attacks (Park et al., 2002).
Gene Expression Regulating Antioxidant Enzyme Production
Biosynthesis of JA- and ET- activate the production of PR-proteins, cell wall modifying enzymes, salicylic acid synthesis, and phytoalexins as a defense response to induce systemic resistance. In the present investigation, the genes encoding defensive enzymes and proteins such as Thioredoxin (GRMZM2G091481), Glutathione Peroxidase (GRMZM2G144153), Glutathione S-transferase (GRMZM2G116273), Defensin like protein (GRMZM2G392863), protease inhibitor (GRMZM2G004466), protease inhibitor (AC233926.1_FG002), Chitinase (GRMZM2G453805), Beta-glucosidase 18 (GRMZM2G031660) were upregulated with higher expressions in consortium treatment (BT-1). The expression patterns of these genes coincide with the levels of biochemical enzymes produced in all treatments, which indicates that the over expressions of these genes may regulate higher production of enzymatic (SOD, POD, PPO, and Protease) and non-enzymatic (proline) antioxidants to initiate maize defense responses, such as lignification, suberization, insect growth inhibition, and scavenging of ROS.
Genes Mediating Redox Signaling
Enzymatic and non-enzymatic antioxidants reduce the level of ROS and SOD is considered as the first and most important enzyme in the scavenging process as it catalyzes O2– into H2O2. The accumulation of the highest H2O2 is toxic for plants and POD detoxifies them (Bano and Muqarab, 2017; Batool et al., 2020). Peroxidases are the most important enzyme and they participate in different physiological processes of the plant (Almagro et al., 2009). Besides ROS scavenging, peroxidase encoding genes are also involved in lignin and suberin production to link the cell wall components. POD belongs to the PR-protein9 subfamily. The highest expression of POD increases the structure barrier and stops the infection. The expression of peroxidase genes in plants was induced by biotic and abiotic stresses (Sasaki et al., 2007). Thioredoxin genes were found to be associated with ROS scavenging enzymes such as ascorbate, POD, SOD, and catalase through disulfite bridge (Dos Santos and Rey, 2006), Moreover, Glutathione S-transferase detoxify electrophilic compounds in the scavenging process, also interact with plant thioredoxin (Yamazaki et al., 2004). Several stress conditions accumulate high ROS levels and resulted in the increased expression of the thioredoxin gene to participate in oxidative stress-linked signaling to induce high levels of antioxidant enzymes in potato plants (Dos Santos and Rey, 2006).
Genes Inhibiting Insect Digestive Enzyme
Polyphenol oxidase is a plant antioxidant enzyme important for regulating insect feeding, growth, development and defense response against herbivory and its production can be associated with the expression of different proteins and enzymes that can restrict the nutrient uptake of insects and attack their digestive track, e.g., chitinase, protease inhibitors and α-amylase inhibitors (Abebe). The chitinase enzyme is very important and it inhibits insect chewing and killing by degrading the insect skeleton (Zhong et al., 2021). The highest expression of chitinase reduces the insect induced damage and survival rate. Bordoloi et al. (2021) also identified two chitinase genes viz. TEA028279 and TEA019397 in tea plants with high expression under stress conditions with jasmonate and salicylic acid in their promoter upstream region to degrade chitin, of insect’s exoskeleton. Similarly, defensin gene also prevents the insect chewing by inhibiting insect’s digestive enzyme such as α-amylase and protease and makes it difficult to breakdown plant starch. Defensin express continuously in leaves stomatal cells, cell wall, peripheral layer of root, seed, and flowers, which indicate that they play an important role in protecting all entry points of plants from pest and pathogens. Khairutdinov et al. (2017) reported the inhibition of α-amylase activity by PsDef1 defensin gene in pine pest Panolis flammea. The inhibitory action of α-amylase by the defensin gene in Sorghum bicolor, Vigna unguiculata, and Tephrosia villosa have also been studied (Kovaleva et al., 2020). Similar to defensin, protease inhibitors and β-glucosidase are also important for inhibiting insect gut digestive enzymes like proteases to prevent its feeding. Guedes and Cutler (2014) reported the role of insecticidal protease inhibitors in tobacco, rice, oilseed rape, sugarcane, and brassica. Expressions of protease inhibitors have been detected in wounds and induce defense against plant herbivory and pathogen attack (Clemente et al., 2019). Maize, almond, and white mustard defensive β-glucosidases were shown to restrict the digestion of lepidopteran larvae (Vassão et al., 2018).
High accumulation of metabolites contributes to acclimation of maize pre-inoculated with entomopathogenic fungi in response to O. furnacalis stress. In our study, significant accumulation of organic acids, benzoxazinoids, phenylpropanoids and polyketides, and lipid-like molecules were observed with higher levels in consortium treatment. Benzoxazinoid is an important defensive secondary metabolite in maize defense against O. furnacalis (Guo et al., 2019). However, this untargeted metabolome profiling was performed for the first time in maize plants inoculated with a consortium of B. bassiana and T. asperellum in response to O. furnacalis attack. We found that the contents of benzoxazinoids were changed after insect attack but higher contents were seen in inoculated treatments as compared to control. Accumulation of organic compounds plays a key role in defense against biotic and abiotic stresses (Zhang et al., 2020). We found that the accumulation of metabolites belonging to benzoxazinoids, lipid-like molecules, organic compounds and phenylpropanoids can be a common mechanism in plant adaptation to different stresses. JA contents, expression levels of genes involved in JA biosynthetic pathway and JA-dependent metabolite adjustment in lipid-like molecules, which plays a defensive role, increased during O. furnacalis feeding in consortium inoculated plants, indicating its enhanced protective mechanism in JA-dependent maize resistance to O. furnacalis. Similar findings were reported by Guo et al. (2019) in O. furnacalis induced direct and in direct defenses. Jasmonates are responsible for regulating nearly all biosynthetic pathways that lead to metabolites (Martínez-Medina et al., 2021). According to the significant impact of B. bassiana and T. asperellum on JA, significant changes in secondary metabolites were triggered. Similarly, previous research also reported a strong impact on plant metabolites, primary metabolism, and changes in defensive compounds (Martínez-Medina et al., 2021).
Among consortium of biocontrol agents used against pest attack, one acts as a stress inducer and the other acts as a control agent. The combined application of BCAs by different mechanisms effectively controlled the incidence of Duponchelia fovealis in strawberry plants (Araujo et al., 2020).
Conclusion
By comprehensive physiological, biochemical, transcriptomic, metabolomic, and phytohormone analysis in maize inoculated with B. bassiana and T. asperellum, applied singly and in the consortium, we revealed enhanced maize defenses against O. furnacalis herbivory. Our study indicates that consortium of B. bassiana and T. asperellum can synergistically suppress O. furnacalis immune response and enhance maize immune response in the field by directly inducing enhanced expression of defense genes and indirectly activating antioxidants and phytohormone production as compared to normal herbivory-induced defenses. They can regulate the coordinated gene network expression related to defense signaling and antioxidants response and within these networks, certain genes make a key contribution to trait variation. Higher phytohormone content diverse metabolome abundance in consortium treatments was also highlighted in our results indicating the importance of JA, SA, and ethylene in plant defense. The correlation between the expression of these genes and enzymatic (SOD, POD, PPO, and protease) and non-enzymatic (proline) antioxidants, phytohormone synthesis, and metabolite adjustments were identified the first time in B. bassiana OFDH1-5, T. asperellum GDFS1009 consortium inoculated maize against O. furnacalis attack. The advantages of fungal BCAs are economic mass production, easy-to-use, sustainable control efficacy, and environmental safety suggest a bright future of mycoinsecticides in the world. Moreover, this study opens up new insight for future studies to focus on how certain changes in the gene networks and metabolites through biocontrol agents can affect above and below ground maize defense and other pathways. On the basis of current results, we believe that defense against ACB attack in maize is not a single gene-controlled process. Multiple genes are involved at the same time. Future research must focus on metabolic-gene clusters.
Data Availability Statement
The datasets presented in this study is submitted in (https://www.ncbi.nlm.nih.gov/bioproject) repository, BioProject: PRJNA797234.
Author Contributions
ZW and JC contributed to conceptualization and supervision. RB contributed to methodology, formal analysis, data curation, visualization, investigation, preparation, and writing of the original draft. MU contributed to software, formal analysis, review, and editing. YW provided resources. KH, MS, TZ, and SB were involved in the review and editing. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by Agricultural Science and Technology Innovation Program (CAAS-ZDRW202004) and China Agriculture Research System of MOF and MARA.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We acknowledge the lab members of State Key Laboratory for Biology of Plant Diseases and Insect Pest, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing, China and biomarker company for helping us during this study.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.790504/full#supplementary-material
References
Abebe, W. (2021). Review on plant defense mechanisms against insect pests. Int. J. Novel Res. Interdiscip. Stud. 8, 15–39.
Almagro, L., Gomez Ros, L. V., Belchi-Navarro, S., Bru, R., Ros Barcelo, A., and Pedreno, M. A. (2009). Class III peroxidases in plant defence reactions. J. Exp. Bot. 60, 377–390. doi: 10.1093/jxb/ern277
Araujo, E. S., Benatto, A., Rizzato, F. B., Poltronieri, A. S., Poitevin, C. G., Zawadneak, M. A., et al. (2020). Combining biocontrol agents with different mechanisms of action to control Duponchelia fovealis, an invasive pest in South America. J. Crop Prot. 134:105184. doi: 10.1016/j.cropro.2020.105184
Bano, A., and Muqarab, R. (2017). Plant defence induced by PGPR against Spodoptera litura in tomato (Solanum lycopersicum L.). Plant Biol. 19, 406–412. doi: 10.1111/plb.12535
Barah, P., and Bones, A. M. (2015). Multidimensional approaches for studying plant defence against insects: from ecology to omics and synthetic biology. J. Exp. Bot. 66, 479–493. doi: 10.1093/jxb/eru489
Batool, R., Umer, M. J., Wang, Y., He, K., Zhang, T., Bai, S., et al. (2020). Synergistic effect of Beauveria bassiana and Trichoderma asperellum to induce maize (Zea mays L.) defense against the Asian corn borer, Ostrinia furnacalis (Lepidoptera, Crambidae) and larval immune response. Int. J. Mol. Sci. 21:8215. doi: 10.3390/ijms21218215
Bordoloi, K., Dihingia, P., Krishnatreya, D., and Agarwala, N. (2021). Genome-wide identification, characterization and expression analysis of the expansin gene family under drought stress in tea (Camellia sinensis L.). J. Plant Sci. Today 8, 32–44.
Cherry, A. J., Banito, A., Djegui, D., and Lomer, C. (2004). Suppression of the stem-borer Sesamia calamistis (Lepidoptera; Noctuidae) in maize following seed dressing, topical application and stem injection with African isolates of Beauveria bassiana. Int. J. Pest Manag. 50, 67–73.
Clemente, M., Corigliano, M. G., Pariani, S. A., Sanchez-Lopez, E. F., Sander, V. A., and Ramos-Duarte, V. A. (2019). Plant serine protease inhibitors: biotechnology application in agriculture and molecular farming. Int. J. Mol. Sci. 20:1345. doi: 10.3390/ijms20061345
Dhondt, S., Gouzerh, G., Müller, A., Legrand, M., and Heitz, T. (2002). Spatio-temporal expression of patatin-like lipid acyl hydrolases and accumulation of jasmonates in elicitor-treated tobacco leaves are not affected by endogenous levels of salicylic acid. Plant J. 32, 749–762. doi: 10.1046/j.1365-313x.2002.01465.x
Di Lelio, I., Coppola, M., Comite, E., Molisso, D., Lorito, M., Woo, S. L., et al. (2021). Temperature differentially influences the capacity of Trichoderma species to induce plant defense responses in tomato against insect pests. Front. Plant Sci. 12:678830. doi: 10.3389/fpls.2021.678830
Dos Santos, C. V., and Rey, P. (2006). Plant thioredoxins are key actors in the oxidative stress response. Trends Plant Sci. 11, 329–334. doi: 10.1016/j.tplants.2006.05.005
Dou, K., Lu, Z., Wu, Q., Ni, M., Yu, C., Wang, M., et al. (2020). MIST: a multilocus identification system for Trichoderma. Appl. Environ. Microbiol. 86:e01532-20. doi: 10.1128/AEM.01532-20
Ehlting, J., Chowrira, S. G., Mattheus, N., Aeschliman, D. S., Arimura, G., and Bohlmann, J. (2008). Comparative transcriptome analysis of Arabidopsis thaliana infested by diamond back moth (Plutella xylostella) larvae reveals signatures of stress response, secondary metabolism, and signalling. BMC Genomics 9:154. doi: 10.1186/1471-2164-9-154
Gamboa-Tuz, S. D., Pereira-Santana, A., Zamora-Briseno, J. A., Castano, E., Espadas-Gil, F., Ayala-Sumuano, J. T., et al. (2018). Transcriptomics and co-expression networks reveal tissue-specific responses and regulatory hubs under mild and severe drought in papaya (Carica papaya L.). Sci Rep. 8:14539. doi: 10.1038/s41598-018-32904-2
González-Guzmán, A., Raya-Diaz, S., Sacristán, D., Yousef, M., Sánchez-Rodríguez, A. R., Barrón, V., et al. (2021). Effects of entomopathogenic fungi on durum wheat nutrition and growth in the field. Eur. J. Agron. 128:126282.
González-Guzmán, A., Sacristán, D., Quesada-Moraga, E., Torrent, J., del Campillo, M. C., and Sánchez-Rodríguez, A. R. (2020). Effects of entomopathogenic fungi on growth and nutrition in wheat grown on two calcareous soils: influence of the fungus application method. Ann. Appl. Biol. 177, 26–40. doi: 10.1111/aab.12596
Guedes, R. N., and Cutler, G. C. (2014). Insecticide-induced hormesis and arthropod pest management. Pest Manage. Sci. 70, 690–697. doi: 10.1002/ps.3669
Guo, J., Guo, J., He, K., Bai, S., Zhang, T., Zhao, J., et al. (2017). Physiological responses induced by Ostrinia furnacalis (Lepidoptera: Crambidae) feeding in maize and their effects on O. furnacalis performance. J. Econ. Entomol. 110, 739–747. doi: 10.1093/jee/tox060
Guo, J., Qi, J., He, K., Wu, J., Bai, S., Zhang, T., et al. (2019). The Asian corn borer Ostrinia furnacalis feeding increases the direct and indirect defence of mid-whorl stage commercial maize in the field. Plant Biotechnol. J. 17, 88–102. doi: 10.1111/pbi.12949
Guo, Q., Li, X., Niu, L., Jameson, P. E., and Zhou, W. (2021). Transcription-associated metabolomic adjustments in maize occur during combined drought and cold stress. Plant Physiol. 186, 677–695. doi: 10.1093/plphys/kiab050
Gupta, H., Saini, R., Pagadala, V., Kumar, N., Sharma, D., and Saini, A. J. (2016). Analysis of plant growth promoting potential of endophytes isolated from Echinacea purpurea and Lonicera japonica. J. Soil Sci. Plant Nutr. 16, 558–577.
He, K., Wang, Z., Zhou, D., Wen, L., Song, Y., and Yao, Z. (2003). Evaluation of transgenic Bt corn for resistance to the Asian corn borer (Lepidoptera: Pyralidae). J. Econ. Entomol. 96, 935–940. doi: 10.1093/jee/96.3.935
Heidel-Fischer, H. M., Musser, R. O., and Vogel, H. (2018). “Plant transcriptomic responses to herbivory,” in Annual Plant Reviews Online, Vol. 47, eds C. Voelckel and G. Jander, New York, NY: Wiley. 155–196. doi: 10.1002/9781118829783.ch5
Joe, M. M., and Muthukumaran, N. (2008). Role of certain elicitors on the chemical induction of resistance in tomato against the leaf caterpillar Spodoptera litura Fab. Not. Bot. Horti Agrobot. Cluj Napoca 36, 71–75.
Jun, X., Wang, X.-Y., and Guo, W.-Z. (2015). The cytochrome P450 superfamily: key players in plant development and defense. J. Integr. Agric. 14, 1673–1686. doi: 10.1016/s2095-3119(14)60980-1
Khairutdinov, B. I., Ermakova, E. A., Yusypovych, Y. M., Bessolicina, E. K., Tarasova, N. B., Toporkova, Y. Y., et al. (2017). NMR structure, conformational dynamics, and biological activity of PsDef1 defensin from Pinus sylvestris. Biochim. Biophys. Acta Proteins Proteom. 1865, 1085–1094. doi: 10.1016/j.bbapap.2017.05.012
Kim, D., Langmead, B., and Salzberg, S. L. (2015). HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12, 357–360. doi: 10.1038/nmeth.3317
Kovaleva, V., Bukhteeva, I., Kit, O. Y., and Nesmelova, I. V. (2020). Plant defensins from a structural perspective. Int. J. Mol. Sci. 21:5307. doi: 10.3390/ijms21155307
Kwon, S. J., Jin, H. C., Lee, S., Nam, M. H., Chung, J. H., Kwon, S. I., et al. (2009). GDSL lipase-like 1 regulates systemic resistance associated with ethylene signaling in Arabidopsis. Plant J. 58, 235–245. doi: 10.1111/j.1365-313X.2008.03772.x
Landa, B. B., López-Díaz, C., Jiménez-Fernández, D., Montes-Borrego, M., Muñoz-Ledesma, F. J., Ortiz-Urquiza, A., et al. (2013). In-planta detection and monitorization of endophytic colonization by a Beauveria bassiana strain using a new-developed nested and quantitative PCR-based assay and confocal laser scanning microscopy. J. Inverteb. Pathol. 114, 128–138. doi: 10.1016/j.jip.2013.06.007
Langfelder, P., and Horvath, S. (2008). WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics 9:559. doi: 10.1186/1471-2105-9-559
Li, J., Zhu, L., Hull, J. J., Liang, S., Daniell, H., Jin, S., et al. (2016). Transcriptome analysis reveals a comprehensive insect resistance response mechanism in cotton to infestation by the phloem feeding insect Bemisia tabaci (whitefly). Plant Biotechnol. J. 14, 1956–1975. doi: 10.1111/pbi.12554
Lin, C. T., Xu, T., Xing, S. L., Zhao, L., Sun, R. Z., Liu, Y., et al. (2019). Weighted gene co-expression network analysis (WGCNA) reveals the hub role of protein ubiquitination in the acquisition of desiccation tolerance in Boea hygrometrica. Plant Cell Physiol. 60, 2707–2719. doi: 10.1093/pcp/pcz160
Lv, L., Zhang, W., Sun, L., Zhao, A., Zhang, Y., Wang, L., et al. (2020). Gene co-expression network analysis to identify critical modules and candidate genes of drought-resistance in wheat. PLoS One 15:e0236186. doi: 10.1371/journal.pone.0236186
Martínez-Medina, A., Mbaluto, C. M., Maedicke, A., Weinhold, A., Vergara, F., and Van Dam, N. M. (2021). Leaf herbivory counteracts nematode-triggered repression of jasmonate-related defences in tomato roots. Plant Physiol. 187, 1762–1778. doi: 10.1093/plphys/kiab368
Meijer, H. J., and Munnik, T. (2003). Phospholipid-based signaling in plants. Annu. Rev. Plant Biol. 54, 265–306. doi: 10.1146/annurev.arplant.54.031902.134748
Mitchell, C., Brennan, R. M., Graham, J., and Karley, A. J. (2016). Plant defense against herbivorous pests: exploiting resistance and tolerance traits for sustainable crop protection. Front. Plant. Sci. 7:1132. doi: 10.3389/fpls.2016.01132
Motholo, L. F. (2019). Endophytic Establishment of Beauveria bassiana in Wheat (Triticum aestivum) and Its Impact on Diuraphis noxia. Potchefstroom: North-West University- Potchefstroom Campus.
Muvea, A. M., Meyhöfer, R., Subramanian, S., Poehling, H.-M., Ekesi, S., and Maniania, N. K. (2014). Colonization of onions by endophytic fungi and their impacts on the biology of Thrips tabaci. PLoS One 9:e108242. doi: 10.1371/journal.pone.0108242
Netzker, T., Fischer, J., Weber, J., Mattern, D. J., Konig, C. C., Valiante, V., et al. (2015). Microbial communication leading to the activation of silent fungal secondary metabolite gene clusters. Front. Microbiol. 6:299. doi: 10.3389/fmicb.2015.00299
Pan, Y., Zhao, S. W., Tang, X. L., Wang, S., Wang, X., Zhang, X. X., et al. (2020). Transcriptome analysis of maize reveals potential key genes involved in the response to belowground herbivore Holotrichia parallela larvae feeding. Genome 63, 1–12. doi: 10.1139/gen-2019-0043
Park, J. H., Halitschke, R., Kim, H. B., Baldwin, I. T., Feldmann, K. A., and Feyereisen, R. (2002). A knock-out mutation in allene oxide synthase results in male sterility and defective wound signal transduction in Arabidopsis due to a block in jasmonic acid biosynthesis. Plant J. 31, 1–12. doi: 10.1046/j.1365-313x.2002.01328.x
Qi, J., Malook, S. U., Shen, G., Gao, L., Zhang, C., Li, J., et al. (2018). Current understanding of maize and rice defense against insect herbivores. Plant Divers. 40, 189–195. doi: 10.1016/j.pld.2018.06.006
Quesada-Moraga, E., Landa, B. B., Muñoz-Ledesma, J., Jiménez-Diáz, R. M., and Santiago-Alvarez, C. (2006). Endophytic colonization of opium poppy, Papaver somniferum, by an entomopathogenic Beauveria bassiana strain. Mycopathologia 161, 323–329. doi: 10.1007/s11046-006-0014-0
Raya-Díaz, S., Quesada-Moraga, E., Barrón, V., Del Campillo, M. C., and Sánchez-Rodríguez, A. R. (2017). Redefining the dose of the entomopathogenic fungus Metarhizium brunneum (Ascomycota, Hypocreales) to increase Fe bioavailability and promote plant growth in calcareous and sandy soils. Plant Soil 418, 387–404. doi: 10.1007/s11104-017-3303-0
Rivas-Franco, F., Hampton, J. G., Narciso, J., Rostás, M., Wessman, P., Saville, D. J., et al. (2020). Effects of a maize root pest and fungal pathogen on entomopathogenic fungal rhizosphere colonization, endophytism and induction of plant hormones. Biol. Control 150:104347. doi: 10.1016/j.biocontrol.2020.104347
Ryu, S. B. (2004). Phospholipid-derived signaling mediated by phospholipase A in plants. Trends Plant Sci. 9, 229–235. doi: 10.1016/j.tplants.2004.03.004
Sánchez-Rodríguez, A. R., Raya-Díaz, S., Zamarreño, ÁM., García-Mina, J. M., del Campillo, M. C., and Quesada-Moraga, E. (2018). An endophytic Beauveria bassiana strain increases spike production in bread and durum wheat plants and effectively controls cotton leafworm (Spodoptera littoralis) larvae. Biol. Control 116, 90–102. doi: 10.1016/j.biocontrol.2017.01.012
Sarma, B. K., Yadav, S. K., Singh, S., and Singh, H. B. (2015). Microbial consortium-mediated plant defense against phytopathogens: readdressing for enhancing efficacy. Soil Biol. Biochem. 87, 25–33. doi: 10.1016/j.soilbio.2015.04.001
Sasaki, K., Yuichi, O., Hiraga, S., Gotoh, Y., Seo, S., Mitsuhara, I., et al. (2007). Characterization of two rice peroxidase promoters that respond to blast fungus-infection. Mol. Genet. Genomics 278, 709–722. doi: 10.1007/s00438-007-0286-1
Shakeri, J., and Foster, H. A. (2007). Proteolytic activity and antibiotic production by Trichoderma harzianum in relation to pathogenicity to insects. Enzyme Microb. Technol. 40, 961–968.
Tang, J. D., Perkins, A., Williams, W. P., and Warburton, M. L. (2015). Using genome-wide associations to identify metabolic pathways involved in maize aflatoxin accumulation resistance. BMC Genomics 16:673. doi: 10.1186/s12864-015-1874-9
Tzin, V., Fernandez-Pozo, N., Richter, A., Schmelz, E. A., Schoettner, M., Schafer, M., et al. (2015). Dynamic maize responses to aphid feeding are revealed by a time series of transcriptomic and metabolomic assays. Plant Physiol. 169, 1727–1743. doi: 10.1104/pp.15.01039
Umer, M. J., Bin Safdar, L., Gebremeskel, H., Zhao, S., Yuan, P., Zhu, H., et al. (2020). Identification of key gene networks controlling organic acid and sugar metabolism during watermelon fruit development by integrating metabolic phenotypes and gene expression profiles. Hortic. Res. 7:193. doi: 10.1038/s41438-020-00416-8
Vassão, D. G., Wielsch, N., Gomes, A. M. D. M. M., Gebauer-Jung, S., Hupfer, Y., Svatoš, A., et al. (2018). Plant defensive β-glucosidases resist digestion and sustain activity in the gut of a lepidopteran herbivore. Front. Plant Sci. 9:1389. doi: 10.3389/fpls.2018.01389
Wagner, B. L., and Lewis, L. C. (2000). Colonization of corn, Zea mays, by the entomopathogenic fungus Beauveria bassiana. Appl. Environ. Microbiol. 66, 3468–3473. doi: 10.1128/aem.66.8.3468-3473.2000
Wang, L., Feng, Z., Wang, X., Wang, X., and Zhang, X. (2010). DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics 26, 136–138. doi: 10.1093/bioinformatics/btp612
Xie, H., Liu, K., Sun, D., Wang, Z., Lu, X., and He, K. (2015). A field experiment with elevated atmospheric CO2-mediated changes to C4 crop-herbivore interactions. Sci. Rep. 5:13923. doi: 10.1038/srep13923
Yamazaki, D., Motohashi, K., Kasama, T., Hara, Y., and Hisabori, T. (2004). Target proteins of the cytosolic thioredoxins in Arabidopsis thaliana. Plant Cell Physiol. 45, 18–27. doi: 10.1093/pcp/pch019
Yu, X., Xiao, J., Chen, S., Yu, Y., Ma, J., Lin, Y., et al. (2020). Metabolite signatures of diverse Camellia sinensis tea populations. Nat. Commun. 11:5586.
Zhang, Y., Bouwmeester, H. J., and Kappers, I. F. (2020). Combined transcriptome and metabolome analysis identifies defence responses in spider mite-infested pepper (Capsicum annuum). J. Exp. Bot. 71, 330–343. doi: 10.1093/jxb/erz422
Keywords: myco-synergism, Ostrinia furnacalis, Beauveria bassiana, Trichoderma asperellum, antioxidants myco-synergism triggers herbivory-induced defense
Citation: Batool R, Umer MJ, Wang Y, He K, Shabbir MZ, Zhang T, Bai S, Chen J and Wang Z (2022) Myco-Synergism Boosts Herbivory-Induced Maize Defense by Triggering Antioxidants and Phytohormone Signaling. Front. Plant Sci. 13:790504. doi: 10.3389/fpls.2022.790504
Received: 06 October 2021; Accepted: 17 January 2022;
Published: 17 February 2022.
Edited by:
Maurizio Ruzzi, University of Tuscia, ItalyReviewed by:
Enrique Quesada Moraga, University of Córdoba, ColombiaYong Liu, Hunan Academy of Agricultural Sciences, Chinese Academy of Agricultural Sciences (CAAS), China
Copyright © 2022 Batool, Umer, Wang, He, Shabbir, Zhang, Bai, Chen and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jie Chen, amllY2hlbjU5QHNqdHUuZWR1LmNu; Zhenying Wang, d2FuZ3p5NjFAMTYzLmNvbQ==; enl3YW5nQGlwcGNhYXMuY24=