Transient expression of SbDhr2 and MeHNL in Gossypium hirsutum for herbivore deterrence assay with Spodoptera litura

Abstract

Background

Spodoptera litura (Lepidoptera: Noctuidae), commonly known as tobacco cutworm or cotton leafworm, is a polyphagous pest which causes considerable damage to cotton (Gossypium hirsutum) and other crops. Herbivore-induced defence response is activated in plants against chewing pests, in which plant secondary metabolites play an important role. Dhurrinase2 (SbDhr2), a cyanogenic β-glucosidase from Sorghum bicolor, is the key enzyme responsible for the hydrolysis of dhurrin (cyanogenic β-glucosidic substrate) to p-hydroxymandelonitrile. Hydroxynitrile lyase (MeHNL) from Mannihot esculanta catalyses the dissociation of cyanohydrins to hydrogen cyanide and corresponding carbonyl compound, both enzymes play a pivotal role in plant defence mechanism.

Results

SbDhr2 and MeHNL genes were expressed individually and co-expressed transiently in cotton leaves. We examined the feeding response of S. litura to leaves in the choice assay. The S. litura population used in this study showed better feeding deterrence to leaves co-expressing both genes compared with the expression of an individual gene.

Conclusion

Our results suggest that co-expression of SbDhr2 and MeHNL genes in cotton leaves demonstrate feeding deterrence to S. litura. Engineering cyanogenic pathway in aerial parts of cotton would be an additional defence strategy against generalist pests and can be enhanced against specialist pests.

Background

Spodoptera litura (Lepidoptera: Noctuidae) commonly known as tobacco cutworm or cotton (Cheng et al. 2017) leafworm is a polyphagous pest, causes considerable damage to cotton (Gossypium hirsutum) and various other crops (Xue et al. 2010; Bragard et al. 2019). Failure to control S. litura and its resistance to various insecticides, lead to humongous economical loss (Ahmad and Gull 2017; Fand et al. 2015). Plant secondary metabolites, play a direct role in plant defence response (War et al. 2012; Schaller 2008) and in the adaptation of plants to abiotic/biotic stresses (Akula and Ravishankar 2011; Bartwal et al. 2013; Gleadow et al. 1998; Rosenthal and Berenbaum 1992). Engineering plant metabolic pathways would be a feasible alternative defence strategy against generalist insect pest. In two-component defence system, β-glucosidases and cyanogenic glucosides are separated by different subcellular compartments (Saunders and Conn 1978; Thayer and Conn 1981; Kesselmeier and Urban 1983; Poulton and Li 1994). In plant physiology, β-glucosidases play diverse roles (Morant et al. 2008), and more than 2 500 species of plants contain Cyanogenic glucoside (Panter 2018). Metabolic engineering of the whole cyanogenic pathway in different plants has been reported for insect herbivory deterrence (Franks et al. 2006; Blomstedt et al. 2016; Tattersall et al. 2001; Bak et al. 2000).

Transient gene expression in cotton using virus-induced gene silencing (VIGS) vector (Becker 2013) (TRV: Tobacco Rattle Virus) has already been performed (Li et al. 2018; Gao et al. 2011; Pang et al. 2013) for functional genomic studies. Expressing cyanogenic pathway enzymes in upland cotton can help to develop insect-pest resistant cotton varieties.

SbDhr2 (dhurrinase2/AF253508.1) from Sorghun bicolor and MeHNL (α-hydroxynitrile lyase/AY787210.1) from Mannihot esculanta are two genes of cyanogenic pathway. SbDhr2 is a monomeric unit (62 kDa) expressed in nodes and leaves. Due to broader substrate specificity SbDhr2 hydrolyses certain artificial substrates (4-Methylumbelliferyl β-D-galactopyranoside, 4-Nitrophenyl-β-D-glucopyranoside, triglochinin) in addition to its natural substrates, dhurrin and sambunigrin (Cicek and Esen 1998; Verdoucq et al. 2004). MeHNL monomer (29.3 kDa) is the key enzyme responsible for the release of hydrogen cyanide (HCN) from α -Hydroxynitrile glucosides (Dadashipour and Asano 2011), having broad substrate specificity (Chueskul and Chulavatnatol 1996; Yan et al. 2003; Cheng et al. 2001; Wajant and Pfizenmaier 1996; Lauble et al. 2002). This paper investigates feeding deterrence and settling preference choice assay with S. litura on cotton leaves, transiently expressing SbDhr2 and MeHNL genes, the former driven by phosphoenolpyruvate carboxylase (PEPC) and the latter by 2X CaMV 35S (duplicated CaMV 35S) promoter (Pauli et al. 2004).

Methods

Plant materials

Seeds of cotton (G. hirsutum) were sown in pots containing peat moss and kept at 23 °C, 200 μmol^. m^-2. S^− 1 light, 65% relative humidity with 16 h/8 h day-night photoperiod in a growth room. After the emergence of four to five true leaves, cotyledons were used for infiltration.

Plasmid construction

SbDhr2 (Verdoucq et al. 2003) gene driven by PEPC promoter (S. vulgare Accession. No X63756.1) was digested with EcoRI and KpnI and ligated in TRV2 pYL156 (pTRV-RNA2) to get the recombinant plasmid pTD2 (Fig. 1a). MeHNL gene digested with EcoRI and SacI was ligated in TRV2 to get the second recombinant plasmid pTH2 (Fig. 1b). Verification of clones was done by restriction enzyme digestion.

Fig. 1

Schematic representation of plasmid vectors maps employed in this study, a pTD2 (pYL156 (pTRV-RNA2) + SbDhr2) here dhurrinase2 gene is under control of PEPC promoter, b pTH2 (pYL156 (pTRV-RNA2) + MeHNL), α-Hydroxynitrile lyase gene is under control of 2X CaMV 35S promoter

Agrobacterium tumefaciens LBA4404

TRV1, TRV2, pTH2, pTD2, independent clones were transformed in A. tumefaciens strain LBA4404 on LB agar plate containing kanamycin (50 mg·L^− 1), rifampicin (125 mg·L^− 1) and streptomycin (50 mg·L^− 1) using the freeze-thaw transformation method (Weigel and Glazebrook 2006). After growing at 28 °C for 2 days, polymerase chain reaction (PCR) was used to select positive transformants that would be used further in this study.

Agroinfiltration in cotton cotyledons

Five hundred microliter inoculum of freshly activated single colonies of A. tumefaciens carrying each binary vector of TRV1,TRV2, pTD2, pTH2 were transferred to flasks containing 50 mL of LB medium supplemented with kanamycin (50 mg·L^− 1), rifampicin (125 mg·L^− 1) and streptomycin (50 mg·L^− 1), 10 mmol·L^− 1 2-(4 morpholino)-ethane sulfonic acid (MES), 20 μmol·L^− 1 acetosyringone (Gao et al. 2011; Pang et al. 2013), then were grown overnight at 28 °C, 160 r·min^− 1. Cells were pelleted at 4 000 r·min^− 1 for 5 min and resuspended in an infiltration buffer containing 10 mmol·L^− 1 MgCl₂, 10 mmol·L^− 1 MES and 200 μmol·L^− 1 acetosyringone (Pang et al. 2013). At 600 nm, O.D. value of the culture was adjusted to 0.9 and suspension was kept for shaking 3 ~ 4 h at 25 °C under 10 r·min^− 1. Agrobacterial culture suspensions were prepared in the ratio of 1:1(Table 1). Then infiltrations were performed in 20 plants for each suspension, at the abaxial side of cotyledons with a needleless syringe (Gao et al. 2011; Senthil-Kumar and Mysore 2014) (Fig. 2).

Table 1 Sets of infiltrations performed

Fig. 2

Sites of Agrobacterium suspension infiltration in cotyledons underside

Gene detection in non-infilterated leaves using PCR

DNA was isolated from non-infiltrated leaves by the cetyl trimethyl ammonium bromide (CTAB) method (Healey et al. 2014). Then PCR analysis was performed using primers listed in Table 2 for SbDhr2 and MeHNL genes. PCR positive leaf samples were further analysed for protein expression.

Table 2 Primers used for PCR screening

Western blotting for confirmation of gene expression

Total protein was extracted from sorghum, cassava, control plants and PCR positive non-infiltrated cotton leaves (Fig. 3). For Western blotting (Trans Blot Turbo transfer system) 40 μg of total protein was transferred on Amersham Hybond- P 0.45 PVDF blotting membrane, as per manufacturer’s instructions. Amersham Hybond-P 0.45 PVDF, a 0.45 μm pore size polyvinylidene difluoride (PVDF) hydrophobic membrane, is used with standard colorimetric and chemiluminescent detection methods for proteins. Blots were probed with primary polyclonal antibodies raised in rabbit for SbDhr2 and MeHNL proteins, detected with horseradish peroxidase (HRP)-conjugated goat anti-rabbit IgG antibody, and chemiluminescence was performed using Pierce™ ECL Western blotting substrate as per manufacturer’s instruction. Membranes were exposed to X-ray film, then developed and fixed.

Fig. 3

Representative SDS-PAGE (mass fraction 10%) analysis of total protein, a lane 1- S. bicolor, M-marker, lane 2- cotton, lane 3-TRV1: TRV2 infiltrated cotton, lanes 4 ~ 7 cotton samples positive for both pTD2 & pTH2 plasmids. b Lane 1- M. esculenta as positive, lane 2- cotton, lane 3- TRV1: TRV2 infiltrated cotton, lane 4- cotton samples positive for both pTD2 & pTH2 plasmids

S. litura herbivory deterrence assay

Three settling preference choice tests (Krothapalli et al. 2013) were conducted in petri-plates containing a control leaf and a leaf transiently singly expressing and co-expressing SbDhr2 and MeHNL genes on wet germination paper (Table 3). Plates were kept at room temperature at a relative humidity of 68% and 16 h/8 h day-night photoperiod. Five of the third instar larvae of S. litura were released in the centre of each plate after being starved for 3 ~ 4 h, and the settling preference of insects was measured after every 24 h period till the 8^th day. Each choice assay was replicated four times. The number of insects on each leaf was used to measure the settling preference and t-test was performed with GraphPad prism-8 for insect count. Mean weight of the five larvae before feeding and post feeding for 2 days (48 h) and on the 8th day was recorded in all three set of tests.

Table 3 Three settling preference choice tests

Cyanogenic capacity (HCNc) in infiltrated leaves

Amount of hydrogen cyanide released per unit time is measured as Cyanogenic Capacity (HCNc) (Hay-Roe et al. 2011). In this experiment, we were interested in temporal comparisons of HCNc in three different experimental conditions. Leaves positive for Western blotting were crushed in 1.5 mL tubes for the qualitative test and Feigl-Anger cyanide test paper (Feigl and Anger 1966) was fixed inside the top portion of caps. Then changing in colour of test paper after every hour was monitored till 10 ~ 12 h ranked according to Hay-Roe et al. (2011).

Results

Validation of plasmid for transient expression

After verification of plasmids, pTH2 and pTD2 by restriction enzyme digestion, a 790 base pairs (bp) MeHNL gene fragment was released from pTH2 vector. A 2.3 kilobase (kb) SbDhr2 gene along with its PEPC promoter gene was released as a restriction digestion product from linearized pTD2 vector (Fig. 4). Colony PCR screening of A. tumefaciens (LBA4404) using gene specific primers confirmed the presence of pTH2 and pTD2 vectors (Fig. 5). Colony PCR gave amplification product of 770 bp for MeHNL gene (Fig. 5a) and a 2.3 kb amplicon of SbDhr2 gene along with its promoter (Fig. 5b).

Fig. 4

Agarose gel electrophoresis of restriction enzyme digestion of pTH2 by EcoRI and SacI, lane 1&3 are uncut pTH2 plasmids, lanes 2 & 4 pTH2 are digested, lanes 5 & 7 are uncut pTD2, lanes 6 & 8 pTD2 are digested with EcoRI and KpnI, M is a 2 000 bp DNA marker

Fig. 5

Colony PCR for confirmation of genes in A. tumefaciens strain LBA4404, a lanes 1~10 pTH2 positive A. tumefaciens colonies, b lanes 1~10 pTD2 positive A. tumefaciens colonies, M: 2 000 bp marker, (+): positive control, B:blank

Screening of genes and expression in non-infiltrated cotton leaves

PCR analysis confirmed the presence of MeHNL (Fig. 6a) and SbDhr2 genes (Fig. 6b). Fifteen leaf samples were screened for presence of MeHNL and SbDhr2 genes, of which 13 samples were PCR positive for each gene. PCR results of negative control leaves confirmed the absence of either gene. Western blot analysis of total protein from five randomly selected PCR positive leaves confirmed the presence of MeHNL (29.3 kDa) (Fig. 7a) and SbDhr2 (62 kDa) (Fig. 7b) proteins bands. Un-infiltrated and empty vector infiltrated leaves of cotton served as negative control.

Fig. 6

PCR analysis of transgene in non-agroinfiltrated leaf samples, DNA isolated from non-infiltrated true leaves. a Lane 1: cotton plant as negative control, lane 2: negative control (TRV1: TRV2 infiltrated), lanes 3 ~ 9: plants infiltrated only with pTH2, lanes 10 ~ 17: plants infiltrated with both pTH2 + pTD2 binary vector, M: 2000 bp marker, H+: positive control for MeHNL gene. b Lane 1: cotton plant as negative control, lane 2: negative control (TRV1: TRV2 infiltrated), lanes 3 ~ 9: plants infiltrated only with pTD2, lanes 10 ~ 17: plants infiltrated with both pTH2 + pTD2 binary vector, D+: positive control for SbDhr2, M: marker (Hind III λ-DNA marker mixed with 2000 bp marker)

Fig. 7

Western Blot analysis of transiently expressed proteins. a Western blotting with polyclonal antibody against MeHNL protein, lane 1: total protein from M. esculanta as positive control, lane 2: cotton as negative control, lane 3: TRV1: TRV2 infiltrated cotton as negative control, lanes 4 ~ 7: cotton samples as positive (co-expressing MeHNL & SbDhr2 proteins) samples. b Western blotting with polyclonal antibody against SbDhr2 protein, lane 1: total protein from S. bicolor as positive control, lane 2: cotton as negative control, lane 3: TRV1 + TRV2 infiltrated cotton as negative control, lanes 4 ~ 7: cotton samples as positive (co-expressing MeHNL & SbDhr2 proteins) samples. Bands of our interest are highlighted; M: protein marker

Insect herbivory measurements

A two-tailed t-test (GraphPad prism-8) for the insect number on a leaf after 48 h indicated that a greater number of pests preferred settling on control leaf (Fig. 8a, c, e) as per observation. On the 8th day, settling preferences were recorded with no significant difference (P > 0.05) in plants singly expressing MeHNL and SbDhr2 genes with respect to the control leaves (Fig. 8b, d). Significant difference (P < 0.05) (Fig. 8f) was recorded in settling preference on the 8th day, where MeHNL and SbDhr2 genes were co-expressed compared with the control leaf. Settling preference results substantiate visual observation of damage caused by pests on control and independently expressed MeHNL and SbDhr2 proteins in leaves and was nearly homogeneous on the 8th day (Fig. 9a ~b) with no significant difference, whereas leaf samples co-expressing both SbDhr2 and MeHNL proteins demonstrated a better feeding deterrence till the 8th day (Fig. 9c ~d) as observed. The damage caused in co-expressed infiltrated leaf is far less compared with control and independently expressing proteins. Mean weight of five S. litura larvae before feeding was 0.020 mg, and the mean weight gained by larvae post feeding on all three sets of tests was recorded after 2 days (48 h) and on the 8th day (Table 4). There was no significant difference observed in weight gained by larvae (Fig. 10).

Fig. 8

Insect (S. litura) settling preference on the control (TRV1: TRV2 infiltrated) leaf and infiltrated leaf, four replicates of each experiment with 5 insects in each plate were performed, to calculate two-way t-test for P –value. a there is a significant difference (P < 0.01) in 48 h, b no significant difference in settling preference (P > 0.14) of the 8th day, c there is a significant difference (P < 0.029) in 48 h, d no significant difference in settling preference (P > 0.099) of the 8th day, e there is a significant difference (P < 0.001) in 48 h, f there is a significant difference in settling preference (P < 0.001) on the 8th day

Fig. 9

Damage caused till the 8th day of feeding, Con: control (TRV1: TRV2), a pTH2 infiltrated, b pTD2 infiltrated, c-i pTH2 + pTD2 co-infiltrated (top view), c-ii pTH2 + pTD2 co-infiltrated (underside view)

Table 4 The mean larvae weight post feeding at different durations on all the three feeding preference tests

Fig. 10

Mean weight of five larvae of S. litura (third instar) fed on three different feeding preference tests and weights recorded post feeding after 2 days and on the 8th day

Cyanogenic capacity (HCNc) in leaves

Cyanide released from transiently expressed positive leaves was measured per unit time (Ballhorn et al. 2010; Alonso-Amelot and Oliveros-Bastidas 2005) to determine cyanogenic capacity (HCNc). No remarkable difference in MeHNL infiltrated or SbDhr2 infiltrated leaves was observed, whereas light colour change was observed in leaves co-expressing both enzymes after 9 ~ 10 h (Table 5).

Table 5 Cyanogenic capacity (HCNc) of crushed leaf tissue monitored over an hourly basis with Feigl-Anger paper

Discussion

Sources of HCN in cotton plants and its detoxification pathways

Cyanogenic glycoside (CNglcs), also known as specialized secondary metabolites, is derived from amino acids, L-tryosine, L-valine, L-leucine, L-isolucine, L-phenylalanine with oximes and cyanohydrins as important intermediates. Enzyme CYC79 family of cytochrome P450 (Andersen et al. 2000; Bak et al. 2006; Jørgensen et al. 2011; Morant et al. 2003) is responsible for the synthesis of majority of oximes with E or Z configuration. The derivatives of oximes play diverse roles in plant defence, growth regulation and communication. Oximes in plants exist in E or Z configurations, with E-oxime having broader biological activities (Sørensen et al. 2018). Cyanohydrins (α-hydroxynitrile) in plants are converted from E-oximes by the action of CYC71 or CYC736 (Jørgensen et al. 2011; Takos et al. 2011) family of enzymes. UDP-glucosyl transferase catalyses the last step in the conversion of a cyanohydrin to a cyanogenic glucoside, dhurrin (UGT85B1), linamarin (UGT85K), prunasin (UGT85A19), lotaustralin (UGT85K). The degradation of cyanogenic glucoside into HCN and aglycones is a two-step process (Jørgensen et al. 2011; Gleadow and Moller 2014) catalysed by cyanogenic β-glucosidases (EC: 3.2.1.21) (Morant et al. 2008; Cressey and Reeve 2019; Esen 1993; Vetter 2017) and α-hydroxynitrile lyases [dhurrin (EC: 4.1.2.11), lotaustralin (EC: 4.1.2.46), linamarin (EC: 4.1.2.46 & 4.1.2.47) and prunasin (EC: 4.1.2.10)] (Dadashipour and Asano 2011; Asano et al. 2005; Kassim and Rumbold 2014).

Tiny amount of hydrogen cyanide is produced by all plants as a product or a co-product of a biosynthesis pathway. There are four reactions/ metabolic pathways that would liberate hydrogen cyanide in cotton on the basis of enzyme predictions (http://ptools.cottongen.org), i.e., ethylene biosynthesis I (plants) pathway (Xu and Zhang 2015), linustatin bioactivation (Schmidt et al. 2018; Jørgensen et al. 2005), neolinustatin bioactivation (Forslund et al. 2004; Lai et al. 2015) and vicianin bioactivation pathway (Mizutani et al. 2007).

HCN in plants is detoxified by two pathways. In the first pathway, HCN is converted to 3-cyano-L-alanine (Machingura et al. 2016), and is further metabolized to L-asparagine and L-aspartate (Asparagine pathway); in the second pathway, thiosulfate sulfurtransferase (rhodanese) (Nakajima 2015; Steiner et al. 2018) catalyses the conversion of thiosulfate and cyanide to thiocyanate and sulfite.

Heterologous expression of SbDhr2 and MeHNL in aerial parts of cotton

We have successfully demonstrated that transient co-expression of SbDhr2 and MeHNL could help to deter S. litura from feeding on cotton leaves. Higher expression of SbDhr2 compared with MeHNL gene was observed in Western blotting, which can be attributed to the choice of promoters, the use of PEPC (Matsuoka et al. 1994) and 2X CaMV 35S (Samac et al. 2004; Christensen et al. 1992; Weeks et al. 1993). According to hydrogen cyanide release detection by using Fiegl-Anger test paper, no colour change was observed in control and leaf tissue independently infiltrated with pTD2 or pTH2, whereas light colour change was observed after more than 9 h in leaf tissue samples co-infiltrated with pTD2 & pTH2 construct.

Bioinformatic analysis

Cyanogenesis in green tissue of cotton (Radin 1972), presence of (R)-mandelonitrile lyase-like enzyme (UniProtKB - A0A1U8PEZ9) E.C: 4.1.2.10 and β-glucosidae (UniProtKB - Q7XAS3) E.C: 3.2.1.21 in G. hirsutum have been reported, the former having 74% amino acid sequence similarity with M. esculanta (R)-mandelonitrile lyase and the latter having 76.1% identity with S. bicolor beta-glucosidase. Phenylalanine N-monooxygenase (E.C: 1.14.14.40) catalysing the conversion of L-phenylalanine to phenylacetaldoxime has been reported in G. hirsutum. Genome sequencing of G. hirsutum reported cyanohydrin beta-glucosyltransferase / uridine diphosphoglucose: aldehydecyanohydrin beta-glucosyltransferase (E.C: 2.4.1.85) gene sequence, the gene involved in the conversion of a cyanohydrin to cyanogenic glucoside. In our previous study (Mahajan et al. 2015), we have explored the possibilities of engineering SbDhr2 enzyme to broaden its substrate specificity, and extensive work on hydroxynitrile lyase regarding improving its substrate specificity and activity (Dadashipour and Asano 2011; Yan et al. 2003; Cheng et al. 2001; Dadashipour et al. 2011; Semba et al. 2008; Semba et al. 2010; Bühler et al. 2003; Lauble et al. 2001) have been published, which can aid in developing herbivory resistant variety of cotton (cyanogenic).

Lessons learned from transient expression

These findings along with the previous work (Pant et al. 2016) indicate that α-hydroxynitrile glucoside exists in G. hirsutum. There are possibilities that cyanide detoxification route/pathway (Gleadow and Moller 2014; Machingura et al. 2016; Ting and Zschoche 1970; Zagrobelny et al. 2004; Miller and Conn 1980; Sun et al. 2018; Nielsen et al. 2016; Pičmanová et al. 2015) is more active in cotton. The key to herbivore deterrence is the rapid release of HCN (Krothapalli et al. 2013) and HCNp (Miller and Conn 1980; Bokanga et al. 1994) against pests (Howe and Jander 2008; Gleadow et al. 2002). Engineering metabolic pathways in acyanogenic /cyanogenic plants are all about trails and errors (Morant et al. 2007).

Conclusion

This study was conducted to investigate whether transient expression of cyanogenic pathway enzymes in aerial parts of cotton protects plants against herbivory by S. litura. The results presented here clearly support the finding that transient co-expression of cyanoamino acid metabolism pathway enzymes can deter S. litura from feeding on cotton leaves. It has also demonstrated that strong green tissue-specific promoter of enzyme/transgene expression is a prerequisite for enhancing HCNp in cotton. These findings extrapolate novel opportunities for metabolic engineering of cyanogenesis in G. hirsutum, for which detailed knowledge of metabolic cross-talk, cyanogenic glucoside synthesis, transport, regulation and degradation is a prerequisite. Engineering cyanogenesis in cotton can be envisioned as an additional pest control strategy.