GhMYB18 confers Aphis gossypii Glover resistance through regulating the synthesis of salicylic acid and flavonoids in cotton plants

doi:10.21203/rs.3.rs-1963511/v1

Download PDF

Research Article

GhMYB18 confers Aphis gossypii Glover resistance through regulating the synthesis of salicylic acid and flavonoids in cotton plants

https://doi.org/10.21203/rs.3.rs-1963511/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Accumulating researches confirmed that plant MYB transcription factors are involved in the defense response to herbivores. However, little is known about the target genes and signaling pathways that MYB transcription factors regulate. In this research, we identified the R2R3 MYB transcription factor GhMYB18 from upland cotton (Gossypium hirsutum L.). GhMYB18 is localized throughout the cell and is a functional transcription activation factor. The expression level of GhMYB18 was significantly up-regulated under aphid infestation and exogenous salicylic acid (SA) treatment. Overexpression of GhMYB18 in cotton cotyledons conferred enhanced tolerance of cotton plants to aphid (Aphis gossypii Glover) feeding, expectedly silencing GhMYB18 increased susceptibility of cotton plants to aphid. Subsequent analysis showed that GhMYB18 positively regulates defense-related enzyme activities including catalase (CAT), peroxidase (POD), polyphenol oxidase (PPO) and phenylalanine ammonia-lyase (PAL), and increases the synthesis of salicylic acid and flavonoids to participate in the defense response against aphid. These results collectively suggest that the cotton MYB gene GhMYB18 is involved in the plant defense response to aphid attack, which may contribute to cotton aphid control strategies.

Gossypium hirsutum L.

Aphis gossypii Glover

R2R3-MYB

plant defense

R2R3 MYB transcription factor GhMYB18 is involved in the defense response to cotton aphid by participating in the synthesis of salicylic acid and flavonoids.

Plants defend themselves against herbivores via various resistance strategies, including constitutive and inducible defenses (Karban and Baldwin 2007). Specific plant defense mechanisms rely on signaling networks, incorporating phytohormones, secondary metabolites, transcription factors, etc (Erb and Reymond 2017). These transcription factors such as WRKY, NAC, bHLH, bZIP, ERF/AP2, and MYB families, play a crucial role in plant resistance to external biotic stresses (Singh et al. 2002). Transcription factors enhance plant adaptation to external stresses by regulating the expression of defense-related genes (Tsuda and Somssich 2015). MYB transcription factor is one of the most numerous and functional gene families in plants, with 51–52 amino acids at its N-terminal called MYB domain. The MYB domain is a highly conserved DNA-binding domain consisting of 1–4 imperfect amino acid sequence repeats (R), each forming three α-helices. The second and third helices of each repeat form a helix–turn–helix (HTH) structure required for binding to the promoters of the target genes (Dubos et al. 2010). Depending on the number of imperfect repeats, MYB proteins are split into four classes: R1-MYB, R2R3-MYB, 3R-MYB, and 4R-MYB, and most MYB transcription factors belong to R2R3-MYB (Stracke et al. 2001). R2R3-MYB transcription factors have been extensively studied in various plants, especially in Arabidopsis (Arabidopsis thaliana), to explore the role of MYBs response to abiotic stresses and biotic stresses (Ambawat et al. 2013). Some MYB proteins have been proved to participate in plant defense against herbivores, including AtMYB102 (De Vos et al. 2006), GsMYB15 (Shen et al. 2018), AtMYB44 (Lu et al. 2013), NaMYB8 (Kaur et al. 2010), TaMYB19, TaMYB29, TaMYB44 (Zhai et al. 2017), yet the regulatory mechanisms and signaling pathways mediated by MYB proteins in defense responses remain largely unknown.

Substantial researches have confirmed that MYB transcription factors contribute to the variation of reactive oxygen species (ROS) (Zhang et al. 2021; Chen et al. 2016; Li et al. 2019). ROS, working as the one-electron reduction forms of atmospheric oxygen, including O₂⁻, H₂O₂, and OH⁻, are generally not harmful to organisms unless they are in adverse circumstances (Dautréaux and Toledano 2007). ROS scavenging enzymes take effect when excess ROS are accumulated in plants. ROS scavenging enzymes, including catalase (CAT), and peroxidase (POD), are essential to maintain the balance of ROS in organisms, which may relate to plant immunity (Zhang et al. 2021). Additionally, polyphenol oxidase (PPO) is connected with the response of plants to diseases, pests, and environmental stresses (Jia et al. 2016; Thipyapong et al. 2004). It's capable of catalyzing the formation of lignin and other phenol oxides to form a defending shield to play a direct resisting role in the course of the defense response (Demeke and Morris 2002).

To date, different lines of evidence have indicated that SA interplays with ROS in stressed plants (Herrera-Vásquez et al. 2015). Salicylic acid (SA) plays a vital role in plant defense and activates resistance responses against external stresses. Mutants that are influenced in the accumulation of SA or insensitive to SA are more susceptive to abiotic and biotic stresses (Bari and Jones 2009). Protein arginine deiminase 4 (PAD4) and enhanced disease susceptibility 1 (EDS1) act upstream of SA to promote SA accumulation, and non-expressor of pathogenesis-related genes 1 (NPR1) is required for SA-induced expression of pathogenesis‐related (PR) genes and systemic acquired resistance (SAR) (Berrocal-Lobo and Molina 2004). Isochorismate synthase (ICS) pathway and the phenylalanine ammonia-lyase (PAL) pathway have been confirmed to be the main ways of SA synthesis (Bernal-Vicente et al. 2020). PAL is also the first key enzyme gene of the phenylpropane pathway, which can regulate the biosynthesis of flavonoids. Flavonoids are the most bioactive secondary metabolites among plants, which can eliminate reactive oxygen species by locating and neutralizing free radicals (Løvdal et al. 2010). The functions of flavonoids are diverse, which includes defense against pathogens and pests, protection against UV light damage and oxidative stress, regulation of auxin transport, and allelopathy (Lloyd et al. 2017). The flavonoid biosynthesis pathway in plants has been extensively studied, the main catalytic enzymes are phenylalanine ammonia-lyase (PAL), chalcone synthase (CHS), chalcone isomerase (CHI), flavonoid 3′-hydroxylase (F3′H), dihydroflavonol 4-reductase (DFR), anthocyanidin synthase (ANS), and flavonol synthase (FLS) (Boudet 2007; Dong and Lin 2021). Those genes which relate to the structural flavonoid pathway are controlled by transcription factors (TFs) including MYBs (Zhai et al. 2016).

The aphid, Aphis gossypii Glover (Hemiptera: Aphididae) is the most common phloem-sucking pest in cotton. It can cause direct damage by feeding or indirect prejudice through transmission of plant viruses (Claude et al. 2018). Existing studies have only confirmed the MYB protein-mediated resistance mechanism to aphid in some plants. For example, TaMYB19, TaMYB29 and TaMYB44, which are regulated by the ethylene signaling pathway, cooperatively participate in wheat phloem defense and influence the feeding behavior of English grain aphid (Zhai et al. 2017). CmMYB15 and CmMYB19 provide chrysanthemum resistance to aphid by regulating the biosynthesis of lignin (An et al. 2019; Wang et al. 2017). The Arabidopsis MYB transcription factors AtMYB102 and AtMYB44 regulate resistance to Myzus persicae by activating ethylene defense in Arabidopsis (Lu et al. 2013; Zhu et al. 2018). In general, the roles of MYBs in cotton defense against aphid are not well known.

In this study, we identified an aphid-responsive gene GhMYB18 from upland cotton (Gossypium hirsutum). Transient expression and Virus-induced gene silencing (VIGS) strategies were employed to clarify the role of GhMYB18 in plant response to Aphis gossypii infestation. Subsequent studies reveal that GhMYB18 can activate defense-related enzymes including CAT, POD, PPO and PAL, induce the synthesis of salicylic acid, and further accumulate flavonoids by regulating phenylpropanoid metabolic pathway in cotton for defense response against aphid.

2.1 Materials and culture condition

Cotton variety Zhongjixing 7 (Gossypium hirsutum L.), kindly provided by professor Dinguo Li from Yangtze University, was used in this study. Cotton and tobacco plants (Nicotiana benthamiana L.) were grown in pots (500 mL) filled with soil mix (vermiculite:humus, 1:1) at 25°C under 16 h light/8 h dark conditions. Cotton aphid was collected from field-grown Zhongjixing 7 plants near Jingzhou, China. The colony was inoculated on cotton seedlings and transferred to fresh plants every 2 weeks (25°C; 16 h/8 h light/dark).

2.2 Phylogenetic Analysis of GhMYB18

The homologous sequences of GhMYB18 were retrieved and analyzed using online BLAST tools from the NCBI (National Center for Biotechnology Information, USA). The sequences of MYB proteins were aligned using DNAMAN. A phylogenetic tree of GhMYB18 with its homologs was constructed using MEGA 7.0 with the maximum likelihood method (bootstraps = 1000).

2.3 RNA extraction and qPCR analysis

Total RNA was extracted using Spectrum™ plant total RNA kit (Sigma-Aldrich, USA) according to the manufacturer’s protocol. Total RNA was reverse transcribed into first-strand cDNA with reverse transcriptase (Promega, USA), which was used as a template for quantitative real-time PCR (qPCR). Cotton gene GhUBI1 (GenBank accession number: EU604080) was used as an internal control. The relative expression level of each gene was calculated using the 2^− MMCT method (Livak and Schmittgen 2001). The relative expression value was shown as mean values of three independent tests, and three replicates were performed for each test. All the primer sequences are shown in Table S1.

2.4 Subcellular localization of GhMYB18

The CDS of GhMYB18 was cloned into the pBI121-GFP vector for an overexpression under the control of the 35S promoter. The recombinant plasmid and empty vector (negative control) were separately transferred into Agrobacterium tumefaciens strain GV3101 competent cells, which were then injected into tobacco leaves. Fluorescence was detected with a confocal laser scanning microscope (Nikon, DS-Ri2, Japan) at 72 h after the injection (Ma et al. 2015).

2.5 Trans-activation activity assays

The full-length CDS of GhMYB18 together with the divided fragments nGhMYB18 (encoding 1-153 amino acids in the N-terminus of GhMYB18) and cGhMYB18 (encoding 154–327 amino acids in the C-terminus of GhMYB18) were ligated to the pGBKT7 vector (Clontech, Mountain View, CA, USA). Those resultant constructions which were designated as BD-GhMYB18, BD-nGhMYB18 and BD-cGhMYB18, were transferred into the Yeast strain AH109 (Clontech). The pGBKT7 empty vector was used as the negative control. The transformed strains were cultured on a medium at 28°C for 3 days before observation. Four types of media were used: SD/-Trp, SD/-Trp with aureobasidin A (AbA), SD/-Trp with X-α-Gal, and SD/-Trp with AbA and X-α-Gal (Jiang et al. 2021).

2.6 Expression profile of GhMYB18 in response to aphid infestation and phytohormone treatment

Uniformly grown cotton seedlings (14 days old) were used in this experiment and each cotton seedling cotyledon was inoculated with 20 adult cotton aphids of the same size. The expression levels of GhMYB18 in roots, stems and cotyledons were analyzed at 24 h, 48 h, and 72 h after the cotton aphid infestation. A solution of methyl jasmonate (MeJA) and salicylic acid (SA) with a concentration of 0.05 mmol/L was prepared and sprayed evenly on cotton cotyledons. The expression levels of GhMYB18 were analyzed in cotton cotyledons at 24 h, 48 h, and 72 h after MEJA and SA treatment. Untreated cotton seedlings were used as a negative control, three biological replicates were used for the above treatments.

2.7 Transient expression

The vector pBI121-GhMYB18 and empty vector were ligated into Agrobacterium tumefaciens GV3101 (Sun et al. 2015). The transformed strains were grown in a lysogeny broth (LB) liquid medium (pH = 7.0) containing 50 ug/mL kanamycin and rifampin at 28°C, 220 rpm. After centrifugation, the pelleted bacteria were resuspended (OD₆₀₀ = 0.8). The resuspension was injected into the cotton cotyledons by a needle-less syringe followed by the expression detection of GhMYB18 (Yue et al. 2012).

2.8 Virus-induced gene silencing

The pTRV1 and pTRV2 vectors were used for VIGS experiments. The specific fragment of GhMYB18 was inserted into pTRV2 (Liu et al. 2002). Agrobacterium tumefaciens cultures (OD₆₀₀ = 1.5) harboring pTRV1 and that containing pTRV2-GhMYB18 were mixed at a 1:1 ratio and inoculated into cotton plants by vacuum infiltration as described by Qu et al (Qu et al. 2012).

2.9 The resistance of GhMYB18-overexpressing and GhMYB18-silenced cotton against cotton aphid

No-choice and choice assays were performed using wild-type plants, control plants and treatment group plants. For no choice assay, 10 adult cotton aphids of the same individual size were inoculated on the cotton cotyledons, and the cotyledons were fixed with insect cages to prevent the aphids from escaping. The number of cotton aphids on the cotyledons was counted at 24 h, 48 h, and 72 h after inoculation. For the choice assay, the detached cotton cotyledons leaves were placed symmetrically on the inner edge of a petri dish, and the petioles are soaked with MS nutrient solution to prevent the leaves from withering. The leaves were inoculated with 30 adult cotton aphids of the same individual size and the number of aphids on the cotyledons was counted at 24 h after inoculation.

2.10 Quantification of aphid excreted honeydew

Sterilized agar medium was poured into a glass petri dish with a diameter of 9 cm, and the adaxial surface of the leaves were placed on the medium without leaving a gap. The cotyledons leaves were inoculated with 10 adult cotton aphids of the same individual size, the petri dish was inverted, and the aphid honeydew was collected with Whatman filter paper for three days. The determination method of aphid honeydew was as follows: the filter papers with collected honeydew were immersed in 0.1% (w/v) ninhydrin solution. After drying, the filter papers were torn into pieces, and stains were extracted in 1 mL of 90% (v/v) methanol for 1 h at 4°C with continuous oscillation. After centrifugation at 5000×g for 1 min, the absorbance was measured at 500 nm (Kim and Jander 2007; Nisbet et al. 1994).

2.11 Determination of enzyme activity

Cotton leaves were frozen in liquid nitrogen and quickly ground into a fine powder. Samples of 100 mg were homogenized in extraction buffer and the mixture was collected in a centrifuge tube. After centrifugation at 4°C and 14000×g for 10 min, the supernatant was the crude enzyme solution. The activities of CAT, POD, PAL and PPO were determined according to the guidelines of the enzyme activity kit (Nanjing Jiancheng Bioengineering Institute, China).

2.12 Determination of salicylic acid content

The salicylic acid was extracted and purified according to the method in the previous study (Fu et al. 2012). 100 mg of fresh leaf tissue was fully ground in liquid nitrogen, suspended in 1 mL of methanol, and incubated at -20°C overnight. The mixture was centrifuged at 14000×g for 10 min at 4°C and purified with a 0.25 um needle nylon filter. 200 µL of the filtrate was taken and dried with nitrogen gas, then 200 µL of sodium acetate buffer (0.1 M, pH = 5.5) was added. The mixture was shaken at 4°C for 1 h and centrifuged at 14000 ×g for 10 min at 4°C. The suspension was treated with 10 µL of β-glucosidase, incubated in a water bath at 37°C for 2 h, and then boiled in water for 5 min to stop the enzymatic reaction. After that, the mixture was centrifuged at 14000×g for 10 min at 4°C, and then the supernatant was subjected to LC-MS to analyze the total SA content.

2.13 Determination of total flavonoids and free gossypol content

100 mg of leaf powder and 1 ml 80% (v/v) methanol were placed in a centrifuge tube and then shook overnight at 4 ℃. After centrifugation at 13000×g at 4 ℃ for 15 min, the supernatant was collected, and the residual pellet was reextracted with 1 mL of 80% (v/v) methanol. The supernatant of the two extractions was combined for subsequent research. Total flavonoids were determined according to the method described previously (Sultana et al. 2009). 0.1 mL of supernatant was mixed with 0.5 mL distilled water into a test tube, to which was added 30 µL of 5% NaNO₂ (w/v) solution. After 6 min, 60 µL of 10% (w/v) AlCl₃ solution was added. 5 min later, 200 µL of 1 M NaOH was added to terminate the reaction. The volume was brought to 1 mL with distilled water and 200 µL of them was used to measure the absorbance at 420 nm.

The gossypol equivalents were extracted and quantitated following a previously reported phloroglucinol/HCl method (Gao et al. 2013). A standard curve was made using gossypol (Sinopharm Chemical Reagent Co LTD, China).

2.14 Statistical Analysis

SPSS v26.0 was used for all statistical analysis, and all data was shown as mean ± SD of one representative experiment. The student’s t-test was used to compare the significant differences in GhMYB18 expression levels after aphid infestation and phytohormone treatment. The p-value was shown as p < 0.05 or p < 0.01 to indicate significant differences. In the case of other experiments, Tukey’s multiple range test (p < 0.05) and one-way analysis of variance was performed to determine the significance of the data. Different letters a, b and c indicated the significance levels (p < 0.05).

3.1 Sequence analysis of GhMYB18

GhMYB18 was isolated from upland cotton ‘Zhongjixing 7’, and sequence analysis showed that it was a new cotton MYB family gene. The GhMYB18 gene is composed of 981 bp and encodes a polypeptide of 327 amino acids (File S1). GhMYB18 belongs to the R2R3-MYB subfamily, as it has two conserved MYB domains including an R2MYB domain between aa 13–63 and an R3MYB domain between aa 66–114. Alignment analysis indicated that this gene has 67% similarity with GhMYB30 (Fig. 1A). Phylogenetic tree of GhMYB18 and MYB proteins from other plant species shows that GhMYB18 are closely related to GhMYB30 (Fig. 1B).

3.2 GhMYB18 was induced by cotton aphid infestation, SA and MeJA

GhMYB18 was expressed in roots, stems and cotyledons in cotton plants. At 24 h, 48 h and 72 h after the sustained infestation of cotton aphid, the expression levels of GhMYB18 were significantly higher than that of the control in roots, stems and cotyledons of cotton plants. Furthermore, GhMYB18 presented higher gene expression levels in roots, stems and cotyledons at 72 h under cotton aphid infestation. Their expression levels were respectively 4.03, 3.66, and 3.21 times more than that in the wild-type plants (WT) at 72 h (Fig. 2A). In addition, GhMYB18 transcripts were significantly up-regulated in leaves at 24 h, 48 h and 72 h after salicylic acid treatment (Fig. 2B). In contrast, after methyl jasmonic acid treatment, the expression level of GhMYB18 was inhibited and was significantly down-regulated at 48 h compared with control (Fig. 2C). These results show that aphid damage and exogenous SA induction can promote the expression level of GhMYB18, while exogenous MeJA induction may decrease the expression level of GhMYB18.

3.3 GhMYB18 localized throughout the cell and was a functional transcription activation factor

To determine the subcellular site where GhMYB18 functions, a subcellular localization assay was performed. Consistent with the control (35S-GFP) results, in the cells transferred into GhMYB18-GFP, the GFP signal was detected in the cell membrane, cytoplasm and nucleus (Fig. 3). For trans-activation activity assays, yeast harboring pGBKT7-GhMYB18 or pGBKT7 (negative control) grew normally on SD/-Trp medium (Fig. 4B). But on SD/-Trp with AbA medium, only the yeast strain of pGBKT7-GhMYB18 grew normally (Fig. 4C). At the same time, it was no more than yeast transferred into pGBKT7-GhMYB18 turned blue on the chromogenic medium (SD/-Trp with AbA and X-α-Gal medium) (Fig. 4D). To further explore the transcriptional activation domain (TAD), the CDs of GhMYB18 was incised into two fragments (N-terminal and C-terminal regions) based on its structural characteristics. As shown in Fig. S1, the BD-nGhMYB18 failed to grow on SD/-Trp with AbA and X-a-Gal medium, indicating that the transcriptional activation domains of GhMYB18 located in the C-terminal regions of GhMYB18. These results illustrate that GhMYB18 localizes throughout the cell and is a functional transcription activation factor.

3.4 Overexpression of GhMYB18 enhances cotton plants’ tolerance to Aphis gossypii

The GhMYB18 expression level in transiently overexpressed GhMYB18 cotton seedlings was determined using qPCR. Compared with wild-type cotton plants (WT) and cotton plants containing pBI121(pBI121 plants), the expression level of GhMYB18 in cotton plants containing pBI121-GhMYB18 (pBI121-GhMYB18 plants) significantly increased at 24 h, 48 h 72 h. The pBI121-GhMYB18 plants presented a higher expression level of GhMYB18 at 48 h, and its expression level was 6.26 times higher than that of pBI121 plants (Fig. 5A). The choice assay showed that the cotton aphid was more widely distributed on pBI121 plants and WT than pBI121-GhMYB18 plants at 24 h, 48 h and 72 h after sustained feeding of cotton aphid (Fig. 5B). Obviously, WT and pBI121 plants were more attractive to cotton aphid. Simultaneously, we conducted a no-choice assay to examine the aphid resistance of pBI121-GhMYB18 plants. Consistently, compared with WT and pBI121 plants, the number of cotton aphids on pBI121-GhMYB18 plants was lower at 24 h, 48 h, and 72 h after persistent feeding of cotton aphid (Fig. 5C). Assays of honeydew secretion from Aphis gossypii was consistent with the above conclusions (Fig. S2). Results showed that the amount of honeydew secreted on pBI121-GhMYB18 plants was significantly lower than that of WT and pBI121 plants at 48 h and 72 h (Fig. 5D). Above all, GhMYB18 overexpression may enhance cotton plant tolerance to Aphis gossypii.

3.5 Silencing of GhMYB18 increased cotton plant susceptibility to Aphis gossypii.

The VIGS method was also employed to further clarify the role of GhMYB18 in cotton response to the infestation of Aphis gossypii. Compared with the wild-type plants (WT) and control group plants (TRV:00), the expression level of GhMYB18 was reduced in GhMYB18-silenced plants (TRV: GhMYB18) at 18 d and 24 d after Agrobacterium tumefaciens infiltration. The GhMYB18-silenced plants presented the lowest expression level of GhMYB18 at 18 d after the silencing treatment, which was only 49.5% of control group plants (Fig. 6A). The choice tests showed that cotton aphid was more inclined to feed on GhMYB18-silenced cottons than WT and control group plants at 18 d (Fig. 6B). In the no-choice assay, the number of cotton aphid in the GhMYB18-silenced plants was significantly higher than that of WT and control group plants after 18 d (Fig. 6C). The quantity of honeydew secreted on GhMYB18-silenced cotton was significantly higher than others at 18 d and 24 d after aphid infestation (Fig. 6D). These results suggest that silencing of GhMYB18 increased cotton plant susceptibility to Aphis gossypii.

3.6 Influence of altering GhMYB18 expression on enzyme activities

As a specific secondary metabolite in cotton, gossypol, conduce to the defense response against cotton aphid (Du et al. 2004). We first measured the content of free gossypol, there was no significant difference between treated and control plants (Fig. S3). Enzyme activities of PAL, CAT, PPO and POD go hand in hand with plant resistance to insect pests. Therefore, these enzyme activities and their gene expression levels were measured. The gene expression level and activities of PAL, CAT, PPO and POD in pBI121-GhMYB18 plants were always higher than those in pBI121 plants (Fig. 7A, B). In the same way, the gene expression levels and activities of these enzymes in GhMYB18-silenced plants were always lower than those in WT and control group plants (TRV:00) (Fig. 7C, D). These results prove that the expression level of GhMYB18 is closely related to PAL, CAT, PPO and POD.

3.7 GhMYB18 acts as a positive regulator in SA signaling pathway

Salicylic acid is irreplaceable in plant response to aphid infestation. Therefore, we analyzed the expression levels of the genes (EDS1, PAD4, ICS, PAL, NPR1) related to the salicylic acid signaling pathway and the content of salicylic acid in GhMYB18 transiently overexpressed cotton plants and GhMYB18-silenced cotton plants. Compared with WT and pBI121 plants, EDS1 and PAL transcripts were significantly higher in the GhMYB18-overexpressing plants, but the PAD4 transcripts had no significant difference. Moreover, the expressed levels of ICS and NPR1 decreased in the GhMYB18-overexpressing plants than in WT plants (Fig. 8A). A similar pattern was observed in GhMYB18-silenced plants, the expression levels of EDS1, ICS, PAL were down-regulated than that WT and control group plants (TRV:00), but the expression levels of PAD4 and NPR1 had no significant difference (Fig. 8B). The SA content was significantly increased in plants overexpressing GhMYB18, but decreased in GhMYB18-silenced plants (Fig. 8C, D). These results show that GhMYB18 may participate in the defense response to cotton aphid by regulating the SA signaling pathway.

3.8 GhMYB18 stimulates the production of total flavonoids by the PAL pathway

Previous studies have shown that PAL is involved in the biosynthesis of polyphenolic compounds, including flavonoids and lignin in plants (Rohde et al. 2004). To further explore the causal relationship between flavonoids or lignin content and expression level of GhMYB18, we detected related indicators in GhMYB18 transiently overexpressed cottons and GhMYB18-silenced cottons. Firstly, we examined the expression levels of genes related to flavonoids and lignin synthesis pathway. The expression level of CHS, CHI, FLS and ANS was up-regulated significantly in plants overexpressing GhMYB18, but compared with WT plants, the transcripts of F3’H and FLS significantly decreased in plants overexpressing GhMYB18. More interestingly all these genes were significantly down-regulated as expectedly in GhMYB18-silenced plants (Fig. 9A, B). Furthermore, it seems that GhMYB18 is not needed in lignin synthesis pathway (Fig. S4). At the same time, the total flavonoid content was significantly increased in plants overexpressing GhMYB18, but decreased in GhMYB18-silenced plants (Fig. 9C, D). These results indicate that GhMYB18 regulates the defense response of cotton to Aphis gossypii by influencing the synthesis of flavonoids.

As one of the largest transcription factors families in plants, MYBs play a vital role in regulating abundant plant-specific physiological processes (Jiang and Rao 2020). Compared with extensive studies on the MYB transcription factors in plant growth and abiotic stress (Loguercio et al. 1999; Lu et al. 2017), little is known about MYB proteins regarding their functions in resisting agriculturally pests such as aphid. In this work, a novel resistance-related R2R3-MYB gene, GhMYB18, was identified and cloned. Subsequent researches indicate that GhMYB18 confers Aphis gossypii resistance by regulating the defense-related enzyme activities, salicylic acid and phenylalanine ammonia-lyase pathways in cotton.

Unlike most MYB transcription factors, GhMYB18 is a functional transcription activation factor, not only localized in the nucleus, but also throughout the cell. Phylogenetic analysis showed that GhMYB18 is closely related to GhMYB30, AtMYB94, and GhMYB60, which are related to the plant immune system (Zhou et al. 2022). Moreso, the expression levels of GhMYB18 were significantly up-regulated after aphid infestation. To investigate the role of GhMYB18 in the defense response to the aphid, transient overexpression and VIGS methods were applied, and there were no phenotypic differences among wild-type, control and treatment plants (Fig. S5). Overexpression of GhMYB18 increased resistance in plants against aphid, reduced the reproduction and feeding rate of aphid. The opposite results were shown after silencing GhMYB18 in cotton plants. These results indicated that GhMYB18 contributes to the resistance of cotton aphid.

In response to aphid infestation, plants have evolved multiple strategies to resist aphid, including trichome formation, callose deposition, and the production of secondary metabolites such as phenols (Goggin 2007). Although there are many related studies, our understanding of plant defense mechanisms against aphid is far from enough (Nalam et al. 2019). Prior researches confirm the potential role of CAT, POD, PAL and PPO, in plant defense against aphid infestation, such as the activities of the POD and CAT enzyme in tobacco were enhanced to respond to aphid attacks (Ren et al. 2014); the resistant cultivars exhibit greater constitutive activities of POD, PPO and PAL than in susceptible cultivars following aphid infestation (Han et al. 2009); four defense-related enzymes including CAT, POD, PPO increased under high densities of aphid-infestation stress (Wang et al. 2014). In this study, we found that overexpression of GhMYB18 increased the activities and expression levels of PAL, CAT, PPO and POD in cotton seedlings. When the expression of GhMYB18 is suppressed, these results were opposite. Consistent with our findings, the cotton gene GhRac6 improves the plant defense response to aphid feeding by improving activities of PAL and CAT (Yang et al. 2018); overexpression of GhChi6 increases the activity of PPO, thus increasing plant immunity to aphid (Zhong et al. 2021).

Salicylic acid (SA) plays a crucial role in minimizing infections and the effects of biotic and abiotic stresses (Zhang and Li 2019). The salicylic acid signaling pathway is activated when plants are exposed to piercing-sucking insects, such as Diuraphis noxia Mordvilko (Mohase and Westhuizen 2002), Nilaparvata lugens Stal (He et al. 2019), Aphis gossypii Glover (Zhong et al. 2021) and Myzus persicae Sulzer (Moran and Thompson 2001). Exogenous treatment with plant hormones on cotton seedlings resulted in significant up-regulation of GhMYB18 after SA induction, but MeJA induction down-regulated the expression level of GhMYB18. The reciprocal antagonism of SA and jasmonic acid (JA) signaling pathways has been demonstrated in numerous studies (Koornneef and Pieterse 2008). For example, when Bemisia tabaci Gennadius nymphs fed on Arabidopsis, JA-responsive defenses were suppressed downstream by inducing SA (Zarate et al. 2007); exogenous SA reduced the amount of JA that Egyptian cotton worm induced in both NPR1 mutant and wild-type Arabidopsis (Stotz et al. 2002). In this research, overexpressing of GhMYB18 in cotton plants increased the SA content and the expression levels of EDS1, ICS and PAL in the SA signaling pathway. By inhibiting the expression of GhMYB18, we then concluded that: GhMYB18 may be involved in the defense response to aphid through regulating the SA signal pathway and further research is still needed to know whether the expression of GhMYB18 affects the JA signaling pathway.

PAL is the first enzyme induced by the phenylpropanoid pathway in the biosynthesis of various phenolic compounds with structural and defense-related functions such as lignin, pigment, flavonoid, and phytoalexin (Solekha et al. 2020). These compounds have a significant effect on plant growth, development, reproduction, defense, and environmental responses (Qian et al. 2019). Different studies have reported the potential role of flavonoids and lignin in plant defense response to aphid. For example, when cotton plants are subjected to salt stress, the content of the flavonoids in the cotton will increase to promote the cotton’s defense response against aphid (Wang et al. 2015); more flavonoids and phenol will accumulate followed by aphid infestation in melon (Zahedi et al. 2019); overexpression of GhLac1 increased the lignin content in cotton associated with increased tolerance to aphid (Hu et al. 2018). Moreso, we also determined the content of lignin and flavonoids in GhMYB18-overexpressed and GhMYB18-silenced plants. In our findings, the expression of GhMYB18 did not affect the synthesis of lignin, but the content of total flavonoids and the expression levels of flavonoid biosynthetic pathway-related genes increased significantly in the plants overexpressing GhMYB18, but decreased in the plants silencing GhMYB18. Above all, the results suggested that GhMYB18 may regulate the defense response against aphid through the synthesis of flavonoids.

In summary, our research reveals the GhMYB18-mediated plant defense response against aphid. GhMYB18 expression was induced by aphid infestation and activates defense-related enzyme activities (CAT, POD, PPO and PAL). Moreover, GhMYB18 acts as a positive regulator in SA signaling pathway and flavonoid synthesis pathway to improve cotton resistance to aphid. Collectively, our research provides new insight into the mechanism of transcriptional regulation in plant defense response to aphid, offering valuable evidence of upland cotton breeding. However, on account of the complex interaction between plants and aphid, more studies on GhMYB18 response to cotton aphid will be performed in the future.

This study suggests that GhMYB18 functions as a positive regulator of aphid tolerance. GhMYB18 acted as a transcription activation factor and was only localized in the cell membrane, cytoplasm and nucleus. Exogenous salicylic acid treatment and aphid treatment significantly increased GhMYB18 expression levels, except for JA induction. The results of the choice assay, no-choice assay and assays of honeydew secretion from Aphis gossypii correlated to the expression level of GhMYB18 upon aphid infestation. Subsequent analysis showed that GhMYB18 not only activated defense-related enzymes such as CAT, POD, PPO and PAL, but also participated in salicylic acid and phenylpropane signaling pathways to improve upland cotton resistance to Aphis gossypii. Our results provide a significant advancement in figuring out the functions of R2R3-MYB transcription factors in aphid resistance.

Author contributions

J.M. Zhang, Y.Z. Yang, H.L. An, D. Xu, X.P. Wang and P. Wan conceived and designed the research. Z.W. Hu, X. Zhong and H.R. Zhang performed the experiments; Z.W. Hu, X.C. Luo and Y.X. Wang analyzed the data and wrote the manuscript. Y. Wang, T. Liu and Y. Zhang provided helpful comments and discussions. All authors read and approved the final manuscript.

Funding

Funding Open Access funding enabled and organized by JM Zhang. This work was supported by the National Natural Science Foundation of China (Grant No. 32172400) and open Fund from Key Laboratory of Integrated Pests Management on Crops in Central China, Hubei Key Laboratory of Crop Diseases, Insect Pests and Weeds Control.

Acknowledgements

This work was supported by the National Natural Science Foundation of China (Grant No. 32172400)

Conflict of interest

The authors do not have any conflict of interest.

Ethical approval

This article does not contain any studies with human participants or animals performed by any of the authors.

Ambawat S, Sharma P, Yadav NR, Yadav RC (2013) MYB transcription factor genes as regulators for plant responses: an overview. Physiol Mol Biol Pla 19(3):307-321. http://dx.doi.org/10.1007/s12298-013-0179-1
An C, Sheng L, Du X, Wang Y, Zhang Y, Song A, Jiang J, Guan Z, Fang W, Chen F, Chen S (2019) Overexpression of CmMYB15 provides chrysanthemum resistance to aphids by regulating the biosynthesis of lignin. Hortic Res 6:84. https://doi.org/10.1038/s41438-019-0166-y
Bari R, Jones JD (2009) Role of plant hormones in plant defense responses. Plant Mol Biol 69(4):473-488. https://doi.org/10.1007/s11103-008-9435-0
Bernal-Vicente A, Petri C, Hernández JA, Diaz-Vivancos P (2020) Biochemical study of the effect of stress conditions on the mandelonitrile‐associated salicylic acid biosynthesis in peach. Plant Biol 22(2):277-286. https://doi.org/10.1111/plb.13066
Berrocal-Lobo M, Molina A (2004) Ethylene response factor 1 mediates Arabidopsis resistance to the soilborne fungus Fusarium oxysporum. Mol Plant Microbe Interact 17(7):763-770. https://doi.org/10.1094/MPMI.2004.17.7.763
Boudet AM (2007) Evolution and current status of research in phenolic compounds. Phytochemistry 68(22-24):2722-2735. https://doi.org/10.1016/j.phytochem.2007.06.012
Chen BS, Niu FF, Liu WZ, Yang B, Zhang JX, Ma JY, Cheng H, Han F, Jiang YQ (2016) Identification, cloning and characterization of R2R3-MYB gene family in canola (Brassica napus L.) identify a novel member modulating ROS accumulation and hypersensitive-like cell death. DNA Res 23(2):101-114. https://doi.org/10.1093/dnares/dsv040
Claudel P, Chesnais Q, Fouché Q, Krieger C, Halter D, Bogaert F, Meyer S, Boissinot S, Hugueney P, Ziegler-Graff V, Ameline A, Brault V (2018) The aphid-transmitted Turnip yellows virus differentially affects volatiles emission and subsequent vector behavior in two Brassicaceae plants. Int J Mol Sci 19(8):2316. https://doi.org/10.3390/ijms19082316
Dautréaux B, Toledano MB (2007) ROS as signalling molecules: mechanisms that generate specificity in ROS homeostasis. Nat Rev Mol Cell Biol 8(10):813-824. https://doi.org/10.1038/nrm2256.
Demeke T, Morris C (2002) Molecular characterization of wheat polyphenol oxidase (PPO). Theor Appl Genet 104(5):813-818. https://doi.org/10.1007/s00122-001-0847-3
De Vos M, Denekamp M, Dicke M, Vuylsteke M, Van Loon L, Smeekens SC, Pieterse CM (2006) The Arabidopsis thaliana transcription factor AtMYB102 functions in defense against the insect herbivore Pieris rapae. Plant Signal Behav 1(6):305-311. https://doi.org/10.4161/psb.1.6.3512
Dong NQ, Lin HX (2021) Contribution of phenylpropanoid metabolism to plant development and plant-environment interactions. J Integr Plant Biol 63(1):180-209. https://doi.org/10.1111/jipb.13054
Dubos C, Stracke R, Grotewold E, Weisshaar B, Martin C, Lepiniec L (2010) MYB transcription factors in Arabidopsis. Trends Plant Sci 15(10):573-581. https://doi.org/10.1016/j.tplants.2010.06.005
Du L, Ge F, Zhu S, Parajulee MN (2004) Effect of cotton cultivar on development and reproduction of Aphis gossypii (Homoptera: Aphididae) and its predator Propylaea japonica (Coleoptera: Coccinellidae). J Econ Entomo 97(4):1278-1283. https://doi.org/10.1093/jee/97.4.1278
Erb M, Reymond P (2019) Molecular interactions between plants and insect herbivores. Annu Rew Plant Biol 70(1):527-557. https://doi.org/10.1146/annurev-arplant-050718-095910
Fu J, Chu J, Sun X, Wang J, Yan C (2012) Simple, rapid, and simultaneous assay of multiple carboxyl containing phytohormones in wounded tomatoes by UPLC-MS/MS using single SPE purification and isotope dilution. Anal Sci 28(11):1081-1087. https://doi.org/10.2116/analsci.28.1081
Gao W, Long L, Zhu LF, Xu L, Gao WH, Sun LQ, Liu LL, Zhang XL (2013) Proteomic and virus-induced gene silencing (VIGS) analyses reveal that gossypol, brassinosteroids, and jasmonic acid contribute to the resistance of cotton to Verticillium dahliae. Mol Cell Proteomics 12(12):3690-3703. https://doi.org/10.1074/mcp.M113.031013
Goggin FL (2007) Plant–aphid interactions: molecular and ecological perspectives. Curr Opin Plant Biol 10(4):399-408. https://doi.org/10.1016/j.pbi.2007.06.004
Han Y, Wang Y, Bi JL, Yang XQ, Huang Y, Zhao X, Hu Y, Cai QN (2009) Constitutive and induced activities of defense-related enzymes in aphid-resistant and aphid-susceptible cultivars of wheat. J Chem Ecol 35(2):176-182. https://doi.org/10.1007/s10886-009-9589-5.
He J, Liu Y, Yuan D, Duan M, Liu Y, Shen Z, Yang C, Qiu Z, Liu D, Wen P, Huang J, Fan D, Xiao S, Xin Y, Chen X, Jiang L, Wang H, Yuan L, Wan J (2020) An R2R3 MYB transcription factor confers brown planthopper resistance by regulating the phenylalanine ammonia-lyase pathway in rice. Proc Natl Acad Sci 117(1):271-277. https://doi.org/10.1073/pnas.1902771116
Herrera-Vásquez A, Salinas P, Holuigue L (2015) Salicylic acid and reactive oxygen species interplay in the transcriptional control of defense genes expression. Front Plant Sci 19(6):171. https://doi.org/10.3389/fpls.2015.00171
Hu Q, Min L, Yang X, Jin S, Zhang L, Li Y, Li YY, Ma YZ, Qi XW, Li DQ, Liu HB, Lindsey K, Zhu LF, Zhang XL (2018) Laccase GhLac1 modulates broad-spectrum biotic stress tolerance via manipulating phenylpropanoid pathway and jasmonic acid synthesis. Plant physiol 176(2):1808-1823. https://doi.org/10.1104/pp.17.01628
Jia H, Zhao P, Wang B, Tariq P, Zhao F, Zhao M, Wang Q, Yang T, Fang J (2016) Overexpression of polyphenol oxidase gene in strawberry fruit delays the fungus infection process. Plant Mol Biol Rep 34(3):592-606. https://doi.org/10.1007/s11105-015-0946-y
Jiang CK, Rao GY (2020) Insights into the diversification and evolution of R2R3-MYB transcription factors in plants. Plant physiol 183(2):637-655. https://doi.org/10.1104/pp.19.01082
Jiang C, Wang D, Zhang J, Xu Y, Zhang C, Zhang J, Wang X, Wang Y (2021) VqMYB154 promotes polygene expression and enhances resistance to pathogens in Chinese wild grapevine. Hortic Res 8(1):151. https://doi.org/10.1038/s41438-021-00585-0
Karban R, Baldwin IT (2007) Induced responses to herbivory. University of Chicago Press, Chicago.
Kaur H, Heinzel N, Schöttner M, Baldwin IT, Gális I (2010) R2R3-NaMYB8 regulates the accumulation of phenylpropanoid-polyamine conjugates, which are essential for local and systemic defense against insect herbivores in Nicotiana attenuate. Plant Physiol 152(3):1731-1747. https://doi.org/10.1104/pp.109.151738
Kim JH, Jander G (2007) Myzus persicae (green peach aphid) feeding on Arabidopsis induces the formation of a deterrent indole glucosinolate. Plant J 49(6):1008-1019. https://doi.org/10.1111/j.1365-313X.2006.03019.x
Koornneef A, Pieterse CM (2008) Cross talk in defense signaling, Plant physiol 146(3):839-844. https://doi.org/10.1104/pp.107.112029
Li BZ, Fan RN, Guo SY, Wang PT, Zhu XH, Fan YT, Chen YX, He KY, Kumar A, Shi JP, Wang Y, Li LH, Hu ZB, Song CP (2019) The Arabidopsis MYB transcription factor, MYB111 modulates salt responses by regulating flavonoid biosynthesis. Environ Exp Bot 166(2):103807. https://doi.org/10.1016/j.envexpbot.2019.103807
Liu Y, Schiff M, Dinesh-Kumar SP (2002) Virus-induced gene silencing in tomato. Plant J 31(6):777-786. https://doi.org/10.1046/j.1365-313x.2002.01394.x
Livak KJ and Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods 25:402-408. https://doi.org/10.1006/meth.2001.1262
Lloyd A, Brockman A, Aguirre L, Campbell A, Bean A, Cantero A, Gonzalez A (2017) Advances in the MYB–bHLH–WD repeat (MBW) pigment regulatory model: addition of a WRKY factor and co-option of an anthocyanin MYB for betalain regulation. Plant Cell Physiol 58(9):1431-1441. https://doi.org/10.1093/pcp/pcx075
Loguercio LL, Zhang JQ, Wilkins TA (1999) Differential regulation of six novel MYB-domain genes defines two distinct expression patterns in allotetraploid cotton (Gossypium hirsutum L.). Mol Gen Genet 261(4):660-671. https://doi.org/10.1007/s004380050009.
Løvdal T, Olsen KM, Slimestad R, Verheul M, Lillo C (2010) Synergetic effects of nitrogen depletion, temperature, and light on the content of phenolic compounds and gene expression in leaves of tomato. Phytochemistry 71(5-6):605-613. https://doi.org/10.1016/j.phytochem.2009.12.014
Lu BB, Li XJ, Sun WW, Li L, Gao R, Zhu Q, Tian SM, Fu MQ, Yu HL, Tang XM, Zhang CL, Dong HS (2013) AtMYB44 regulates resistance to the green peach aphid and diamondback moth by activating EIN2‐affected defenses in Arabidopsis. Plant Biol 15(5):841–850. https://doi.org/10.1111/j.1438-8677.2012.00675.x
Lu N, Roldan M, Dixon RA (2017) Characterization of two TT2-type MYB transcription factors regulating proanthocyanidin biosynthesis in tetraploid cotton, Gossypium hirsutum. Planta 246(2):323-335. https://doi.org/10.1007/s00425-017-2682-z
Ma Y, Dai X, Xu Y, Luo W, Zheng X, Zeng D, Pan Y, Lin X, Liu H, Zhang D, Xiao J, Guo X, Xu S, Niu Y, Jin J, Zhang H, Xu X, Li L, Wang W, Qian Q, Ge S, Chong K (2015) COLD1 confers chilling tolerance in rice. Cell 160(6):1209-1221. https://doi.org/10.1016/j.cell.2015.01.046.
Mohase L, Westhuizen AJ (2002) Salicylic acid is involved in resistance responses in the Russian wheat aphid-wheat interaction. J Plant Physiol 159(6):585-590. https://doi.org/10.1078/0176-1617-0633
Moran PJ, Thompson GA (2001) Molecular responses to aphid feeding in Arabidopsis in relation to plant defense pathways. Plant Physiol 125(2):1074-1085. https://doi.org/10.1104/pp.125.2.1074
Nalam V, Louis J, Shah J (2019) Plant defense against aphids, the pest extraordinaire. Plant Sci 279:96-107. https://doi.org/10.1016/j.plantsci.2018.04.027
Nisbet AJ, Woodford JAT, Strang RHC (1994) Quantifying aphid feeding on non-radioactive food sources. Entomol Exp Appl 72:85-89.
Qian Y, Lynch JH, Guo L, Rhodes D, Morgan JA, Dudareva N (2019) Completion of the cytosolic post-chorismate phenylalanine biosynthetic pathway in plants. Nat commun 10(1):1-15. https://doi.org/10.1038/s41467-018-07969-2
Qu J, Ye J, Geng YF, Sun YW, Gao SQ, Zhang BP, Chen W, Chua NH (2012) Dissecting functions of KATANIN and WRINKLED1 in cotton fiber development by virus-induced gene silencing. Plant Physiol 160(2):738-748. https://doi.org/10.1104/pp.112.198564
Ren G, Wang X, Chen D, Wang, XR, Liu X (2014) Effects of aphids Myzus persicae on the changes of Ca²⁺ and H₂O₂ flux and enzyme activities in tobacco. J Plant Interact 9(1):883-888. https://doi.org/10.1080/17429145.2014.982221
Rohde A, Morreel K, Ralph J, Goeminne G, Hostyn V, De Rycke R, Kushnir S, Van Doorsselaere J, Joseleau JP, Vuylsteke M, Van Driessche G, Van Beeumen J, Messens E, Boerjan W (2004) Molecular phenotyping of the pal1 and pal2 mutants of Arabidopsis thaliana reveals far-reaching consequences on phenylpropanoid, amino acid, and carbohydrate metabolism. Plant Cell 16(10):2749-2771. https://doi.org/10.1105/tpc.104.023705.
Shen SJ, Wang YY, Zhang YX, Guo W, Jiao YQ, Zhou XA (2018) Overexpression of the wild soybean R2R3-MYB transcription factor GsMYB15 enhances resistance to salt stress and Helicoverpa Armigera in transgenic Arabidopsis. Int J Mol Sci 19(12):3958. https://doi.org/10.3390/ijms19123958
Singh KB, Foley RC, Oñate-Sánchez L (2002) Transcription factors in plant defense and stress responses. Curr Opin Plant Biol 5(5):430-436. https://doi.org/10.1016/S1369-5266(02)00289-3
Solekha R, Susanto FA, Joko T, Nuringtyas TR, Purwestri YA (2020) Phenylalanine ammonia lyase (PAL) contributes to the resistance of black rice against Xanthomonas oryzae pv. oryzae. J Plant Pathol 102(2):359-365. https://doi.org/10.1007/s42161-019-00426-z
Stotz HU, Koch T, Biedermann A, Weniger K, Boland W, Mitchell-Olds T (2002) Evidence for regulation of resistance in Arabidopsis to Egyptian cotton worm by salicylic and jasmonic acid signaling pathways. Planta 214(4):648-652. https://doi.org/10.1007/s004250100656
Stracke R, Werber M, Weisshaar B (2001) The R2R3-MYB gene family in Arabidopsis thaliana. Curr Opin Plant Biol 4(5):447-456. https://doi.org/10.1016/S1369-5266(00)00199-0
Sultana B, Anwar F, Ashraf M (2009) Effect of extraction solvent/technique on the antioxidant activity of selected medicinal plant extracts. Molecules 14(6):2167-2180. https://doi.org/10.3390/molecules14062167
Sun X, Gong SY, Nie XY, Li Y, Li W, Huang GQ, Li XB (2015) A R2R3‐MYB transcription factor that is specifically expressed in cotton (Gossypium hirsutum) fibers affects secondary cell wall biosynthesis and deposition in transgenic Arabidopsis. Physiol Plant 154(3):420-432. https://doi.org/10.1111/ppl.12317
Thipyapong P, Melkonian J, Wolfe D W, Steffens JC (2004) Suppression of polyphenol oxidases increases stress tolerance in tomato. Plant Sci 167(4):693-703. https://doi.org/10.1016/j.plantsci.2004.04.008
Tsuda K, Somssich IE (2015) Transcriptional networks in plant immunity. New Phytol 206(3):932-947. https://doi.org/10.1111/nph.13286
Wang Q, Eneji AE, Kong X, Wang K, Dong H (2015) Salt stress effects on secondary metabolites of cotton in relation to gene expression responsible for aphid development. PloS one 10(6):e0129541. https://doi.org/10.1371/journal.pone.0129541.
Wang XY, Zhou LH, Xu B, Xing X, Xu GQ (2014) Seasonal occurrence of Aphis glycines and physiological responses of soybean plants to its feeding. Insect Sci 21(3):342-351. https://doi.org/10.1111/1744-7917.12099
Wang Y, Sheng L, Zhang H, Du X, An C, Xia X, Chen F, Jiang J, Chen S (2017) CmMYB19 over-expression improves aphid tolerance in Chrysanthemum by promoting lignin synthesis. Int J Mol Sci 18(3):619 https://doi.org/10.3390/ijms18030619
Yang Z, Liu JJ, Luo L, Ye S, Yang YZ, Zhang GH, Wang XP, Zhang JM (2018) The cotton GhRac6 gene encoding a plant ROP/RAC protein improves the plant defense response to aphid feeding. Plant Mol Biol Rep 36(5):888-896. https://doi.org/10.1007/s11105-018-1127-6
Yue Y, Zhang M, Zhang J, Tian X, Duan L, Li Z (2012) Overexpression of the AtLOS5 gene increased abscisic acid level and drought tolerance in transgenic cotton. J Exp Bot 63(10):3741-3748. https://doi.org/10.1093/jxb/ers069
Zarate SI, Kempema LA, Walling LL (2007) Silverleaf whitefly induces salicylic acid defenses and suppresses effectual jasmonic acid defenses. Plant physiol 143(2): 866-875. https://doi.org/10.1104/pp.106.090035
Zahedi A, Razmjou J, Rafiee-Dastjerdi H, Leppla NC, Golizadeh A, Hassanpour M, Ebadollahi A (2019) Tritrophic interactions of cucumber cultivar, Aphis gossypii (Hemiptera: Aphididae), and its predator Hippodamia variegata (Coleoptera: Coccinellidae). J Econ Entomol 112(4):1774-1779. https://doi.org/10.1093/jee/toz072
Zhai R, Wang Z, Zhang S, Meng G, Song L, Wang Z, Li P, Ma F, Xu L (2016) Two MYB transcription factors regulate flavonoid biosynthesis in pear fruit (Pyrus bretschneideri Rehd.). J Exp Bot 67(5):1275-1284. https://doi.org/10.1093/jxb/erv524
Zhai Y, Li P, Mei Y, Chen M, Chen X, Xu H, Zhou X, Dong H, Zhang C, Jiang W (2017) Three MYB genes co-regulate the phloem-based defense against English grain aphid in wheat. J Exp Bot 68(15):4153-4169. https://doi.org/10.1093/jxb/erx204
Zhang Y, Li X (2019) Salicylic acid: biosynthesis, perception, and contributions to plant immunity. Curr Opin Plant Biol 50:29-36. https://doi.org/10.1016/j.pbi.2019.02.004
Zhang XY, He YQ, Li LY, Liu HR, Hong GJ (2021) Involvement of the R2R3-MYB transcription factor MYB21 and its homologs in regulating flavonol accumulation in Arabidopsis stamen. J Exp Bot 72(12): 4319-4332. https://doi.org/10.1093/jxb/erab156
Zhang S, Liu J, Xu B, Zhou J (2021) Differential responses of Cucurbita pepo to Podosphaera xanthii reveal the mechanism of powdery mildew disease resistance in pumpkin. Front Plant Sci 12:633221. https://doi.org/10.3389/fpls.2021.633221
Zhou YL, Lu L, Liu N, Cao H, Li H, Gui DP, Wang JH, Zhang CJ (2022) Analysis of MYB genes in four plant species and the detection of genes associated with drought resistance. Botany 99(99):1-14. https://doi.org/10.1139/cjb-2020-0227
Zhong X, Feng P, Ma QQ, Zhang Y, Yang YZ, Zhang JM (2021) Cotton Chitinase gene GhChi6 improves the Arabidopsis defense response to aphid attack. Plant Mol Biol Rep 39(1):251-261. https://doi.org/10.1007/s11105-020-01248-5
Zhu L, Guo J, Ma Z, Wang J, Zhou C (2018) Arabidopsis transcription factor MYB102 increases plant susceptibility to aphids by substantial activation of ethylene biosynthesis. Biomolecules 8(2):39. https://doi.org/10.3390/biom8020039

Fig.S1.tif
Fig. S1 Trans-activation assay of BD-nGhMYB18 and BD-cGhMYB18 in yeast. Yeast strains harboring BD control, BD-nGhMYB18 and BD-cGhMYB18 were grown on SD/-Trp (A), SD/-Trp with X-a-Gal (B) and SD/-Trp with AbA and X-a-Gal medium (C) and cultivated at 28 ℃ for 3 days.
FigS2.tif
Fig. S2 The aphid honeydew was collected with Whatman filter paper.
Fig.S3.tif
Fig. S3 The content of free gossypol in GhMYB18 transiently expressed cotton (A) and GhMYB18-silenced cotton (B). Those bars indicate the standard errors, the alphabets represent the level of significant difference(p＜0.05).
Fig.S4.tif
Fig. S4 The expression level of genes related to lignin synthesis pathway in GhMYB18 transiently expressed cotton (A) and GhMYB18-silenced cotton (B). Those bars indicate the standard errors, the alphabets represent the level of significant difference(p＜0.05)
Fig.S5.tif
Fig. S5 Phenotypic analysis of GhMYB18 transiently expressed cotton (A, B) and GhMYB18-silenced (C, D) cotton.
FileS1cottonMYBESTs.docx
TableS1Primersequences.xlsx

Download PDF

Reviewers agreed at journal
25 Aug, 2022
Reviewers invited by journal
22 Aug, 2022
Editor assigned by journal
18 Aug, 2022
First submitted to journal
15 Aug, 2022

You are reading this latest preprint version

GhMYB18 confers Aphis gossypii Glover resistance through regulating the synthesis of salicylic acid and flavonoids in cotton plants

Status:

Version 1

Abstract

Figures

Key Message

1. Introduction

2. Materials And Methods

2.1 Materials and culture condition

2.2 Phylogenetic Analysis of GhMYB18

2.3 RNA extraction and qPCR analysis

2.4 Subcellular localization of GhMYB18

2.5 Trans-activation activity assays

2.6 Expression profile of GhMYB18 in response to aphid infestation and phytohormone treatment

2.7 Transient expression

2.8 Virus-induced gene silencing

2.9 The resistance of GhMYB18-overexpressing and GhMYB18-silenced cotton against cotton aphid

2.10 Quantification of aphid excreted honeydew

3. Results

3.1 Sequence analysis of GhMYB18

3.2 GhMYB18 was induced by cotton aphid infestation, SA and MeJA

3.3 GhMYB18 localized throughout the cell and was a functional transcription activation factor

3.4 Overexpression of GhMYB18 enhances cotton plants’ tolerance to Aphis gossypii

3.5 Silencing of GhMYB18 increased cotton plant susceptibility to Aphis gossypii.

3.6 Influence of altering GhMYB18 expression on enzyme activities

3.7 GhMYB18 acts as a positive regulator in SA signaling pathway

3.8 GhMYB18 stimulates the production of total flavonoids by the PAL pathway

Discussion

Conclusion

Declarations

References

Supplementary Files

Status:

Version 1