Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Poplar protease inhibitor expression differs in an herbivore specific manner

BMC Plant Biology volume 21, Article number: 170 (2021) Cite this article

Abstract

Background

Protease inhibitors are defense proteins widely distributed in the plant kingdom. By reducing the activity of digestive enzymes in insect guts, they reduce the availability of nutrients and thus impair the growth and development of the attacking herbivore. One well-characterized class of protease inhibitors are Kunitz-type trypsin inhibitors (KTIs), which have been described in various plant species, including Populus spp. Long-lived woody perennials like poplar trees encounter a huge diversity of herbivores, but the specificity of tree defenses towards different herbivore species is hardly studied. We therefore aimed to investigate the induction of KTIs in black poplar (P. nigra) leaves upon herbivory by three different chewing herbivores, Lymantria dispar and Amata mogadorensis caterpillars, and Phratora vulgatissima beetles.

Results

We identified and generated full-length cDNA sequences of 17 KTIs that are upregulated upon herbivory in black poplar leaves, and analyzed the expression patterns of the eight most up-regulated KTIs via qRT-PCR. We found that beetles elicited higher transcriptional induction of KTIs than caterpillars, and that both caterpillar species induced similar KTI expression levels. Furthermore, KTI expression strongly correlated with the trypsin-inhibiting activity in the herbivore-damaged leaves, but was not dependent on damage severity, i.e. leaf area loss, for most of the genes.

Conclusions

We conclude that the induction of KTIs in black poplar is controlled at the transcriptional level in a threshold-based manner and is strongly influenced by the species identity of the herbivore. However, the underlying molecular mechanisms and ecological consequences of these patterns remain to be investigated.

Background

Over millions of years plants have developed numerous strategies to defend themselves against plant-feeding animals. Apart from indirect defenses, which involve the recruitment of an herbivore’s natural enemies, plants can harm their attackers directly by producing mechanical barriers, chemical toxins and deterrents, or by using biochemical defenses that interfere with the herbivore’s enzymatic machinery. Among chemical defenses, most emphasis has been placed on low molecular weight metabolites, but defensive proteins exist, such as protease inhibitors (PIs) that reduce the digestibility of plant tissue for the feeding herbivore. By inhibiting proteolytic enzymes in the midgut of the herbivore, PIs diminish protein digestion and hence lower the availability of free amino acids required for herbivore growth and development [15]. The PIs found in plants are numerous and diverse, with 99 different inhibitor families currently described [32]. Those PI families, as well as distinct members within a family, vary in their activity towards the four types of proteases found in herbivore guts, namely serine -, cysteine -, aspartic acid -, and metalloproteases. In herbivorous insects, the most abundant protein-degrading enzymes are the serine proteases [15]. It is therefore not surprising that serine PIs are widely distributed in the plant kingdom [20, 21]. One of the best characterized classes of serine PIs are the Kunitz-type trypsin inhibitors (KTIs; also Kunitz-type protease inhibitors, KPI), of which some are also able to inhibit cysteine proteases [2, 6]. KTIs are relatively small proteins with a mass of 20 to 25 kDa [39], with a β-trefoil structure, consisting of a β-barrel and several loops, of which one is binding to the active site of the target protease [42]. The biological activity of KTIs has been demonstrated by using gut extracts in in vitro assays [16, 29], as well as monitoring the fitness of herbivores feeding on KTI-enriched diets [2, 6, 22, 25, 26, 30]. Since the first description of a KTI in soybean [19, 24], most subsequent studies have also focused on KTIs from legume species [16, 17, 22, 30, 36, 43]. However, KTIs in trees have gained more attention in past years. In species of the genus Populus, several KTIs have been identified and characterized [7, 27, 29, 37, 39], and some shown to be inducible by mechanical wounding or insect herbivory [27,27,29, 39]. For example, feeding by the forest tent caterpillar, a generalist herbivore, increased KTI transcript abundance locally and systemically in hybrid poplar leaves [28]. In fact, genes encoding for KTIs belong to the most up-regulated ones in systemic poplar leaves upon mechanical wounding [9]. In a study by Philippe et al. [39] it was shown that the transcriptional induction triggered by wounding varies among the KTIs and in a time-dependent manner. So far, most studies used Malacosoma disstria, a generalist lepidopteran species, to investigate herbivore-triggered KTI responses in poplar [27, 28, 39]. To our knowledge, the specificity of poplar KTI induction towards other herbivore species has not yet been investigated.

Specificity of response to different herbivores may be especially important for large, long-lived woody perennials like trees, which encounter a vast diversity of herbivores in their lifetimes. For example, it is well known that plants react differently to leaf-chewing herbivores than herbivores feeding on phloem-sap [13, 23]. Specificity of anti-herbivore defenses can also be observed within the same feeding guild, and even within the same species depending on the insect’s developmental stage. For example, early instar generalist caterpillars induced a stronger defense reaction in black poplar leaves than late instar caterpillars of the same species [31]. The underlying mechanism might be explained by HAMPs or DAMPS (herbivore- or damage-associated molecular patterns, respectively) that plants perceive when being attacked [13]. These are influenced by the physical attributes of herbivory, such as leaf area removal or the timing of tissue damage, but also by chemical cues such as salivary compounds of the herbivores [33]. All of these traits can be herbivore species-specific and may allow plants to distinguish among attackers and mount adequate and effective defenses against specific herbivores. In black poplar trees, such herbivore-specific reactions could be shown for signaling molecules [14], as well as chemical defense traits such as volatile emission [14, 31, 47]. In a recent study by Fabisch et al. [14], total PI activity against trypsin was more strongly induced by beetle feeding than by caterpillar feeding on black poplar leaves. However, to date, we do not know which specific genes are responsible for the observed differences in PI activity and whether or not transcription of PI-encoding genes differs between beetle- and caterpillar-fed leaves.

In this study, we therefore tested the hypothesis that different herbivore species induce KTI genes in a species-specific manner. We identified 17 KTI genes from a transcriptome of black poplar and generated full-length cDNA sequences of the most up-regulated ones. Gene expression patterns of these KTI genes as determined by qRT-PCR upon herbivory by three different insect species (Fig. 1) show striking differences among the species.

Fig. 1
figure 1

Insects used in this study and their damage pattern after 2 d feeding on black poplar leaves. Amata mogadorensis and Lymantria dispar (gypsy moth) remove large areas from the leaves, whereas Phratora vulgatissima (blue willow beetle) causes small, but numerous lesions

Full size image

Methods

Plants and insects

Populus nigra L. (Salicaceae) trees were grown from cuttings obtained from trees in a common garden near Jena, Germany. These trees were originally derived from a single female genotype from a P. nigra population (species identified by Sybille Unsicker based on morphological features) located in Küstrin-Kietz, Germany (52°34′1“ N, 14°38’3” E). Since cuttings for this study were taken from trees in a common garden, no permission was necessary for collecting plant material; a voucher specimen will be deposited in spring 2021 in the Herbarium Haussknecht (JE) in Jena, Germany. The cuttings were potted in 2 L pots, grown in the greenhouse (18/20 °C, night/day, relative humidity 60%, natural light with 9–14 h photoperiod, supplemented light for 12 h) and transferred to a climate chamber (18/20 °C, night/day; relative humidity 60%; photoperiod 16 h) 2 days before the onset of the experiment. Trees were either grown for 4 months to approximately 0.5 m (Transcriptome samples) or grown to a height of 1.6 m (approximately 6 months) and pruned back to 0.8 m 4 weeks before treatment (Gene expression samples).

Lymantria dispar L. (Erebidae, Lepidoptera) caterpillars are generalist feeders with a broad host range, preferably deciduous trees. L. dispar caterpillars were hatched from eggs kindly provided by the US Department of agriculture (USDA, Buzzards Bay, MA, USA) and reared on artificial diet (MP Biomedicals LLC, Illkirch, France) in a climate chamber (14/10 h, light/dark, 20–23 °C, relative humidity 60%) until they reached the third instar, the stage used for the experiments. This species is reared continuously at the MPI-CE.

Amata mogadorensis Blachier (Erebidae, Lepidoptera) caterpillars are also generalists with a preference for woody plants and shrubs. A. mogadorensis caterpillars were hatched from eggs provided by a private breeder (www.entomologenportal.de) and reared on black poplar foliage until they reached the third instar, the stage used for the experiment. Individuals were reared until adult stage to confirm the species identity.

Phratora vulgatissima L. (Chrysomelidae, Coleoptera) beetles are specialists, feeding on a narrow range of hosts within the Salicaceae. Beetles (taxonomically determined by Lars Möckel; individuals in alcohol available at the MPI-CE) were reared in the laboratory on black poplar trees.

Experimental designs and sampling

Plant material from two different experiments was used to analyze the transcriptome (see Transcriptome samples) or the gene expression of Kunitz-type trypsin inhibitors (KTIs; see Gene expression samples).

Transcriptome samples

A leaf pool (8 leaves from the stem of a young black poplar tree (n = 4) was wrapped with gauze and then infested with L. dispar caterpillars (4 individuals per tree), adult P. vulgatissima beetles (6 individuals per tree), or left untreated (control). Due to time differences in the availability of the experimental insects, the beetle treatment was conducted two weeks earlier than the caterpillar treatment; both treatments had their own respective control group (n = 4), which was treated and sampled at the same time as the herbivore-treated plants, but was not exposed to herbivores. After 2 d, the treated leaves were flash-frozen in liquid nitrogen and stored at − 80 °C.

Gene expression samples

For gene expression analysis by qRT-PCR, leaf material from an experiment described in Fabisch et al. [14] was used, where further details on the methods are described. In short, a leaf pool (5 leaves) of black poplar trees (n = 10, but a random selection of 6 was used for gene expression analysis) was wrapped with PET bags (Bratschlauch, Toppits, Minden, Germany) and then infested with L. dispar caterpillars (10 per tree), A. mogadorensis (10 per tree), P. vulgatissima beetles (50 per tree), or left untreated (control). After 1 d, the number of caterpillars was reduced to prevent excessive leaf loss. After a total feeding period of 2 d, the leaves were photographed to assess the damage and subsequently flash-frozen in liquid nitrogen and stored at − 80 °C. The damage was quantified as leaf area loss from the photographs by reconstructing the original leaf area in the picture and counting the number of pixels representing the total and the removed leaf areas (Photoshop, Version 15.0.0, Adobe Systems Incorporated, San Francisco, USA). Pixels were converted to area (cm2) using a reference field in the photograph.

RNA isolation and cDNA synthesis

Frozen leaves were ground in liquid nitrogen and RNA was isolated using the InviTrap Spin Plant Mini Kit (Stratec Biomedical AG, Birkenfeld, Germany), including DNase digestion. RNA concentration was measured with a NanoDrop 2000c spectrophotometer (Peqlab Biotechnologie GmbH, Erlangen, Germany). For transcriptome samples, an additional quality check was conducted with the RNA 6000 Nano Kit on a Bioanalyzer (Agilent, Santa Clara, CA, USA). cDNA was synthesized from RNA using SuperScript-III reverse transcriptase and oligo-dT primers (Thermo Fisher Scientific, Waltham, MA, USA).

Transcriptome analysis

Sequencing was done at the Max Planck-Genome-Center (Köln, Germany) on a HiSeq 2500 (Illumina, San Diego, CA, USA) with 9 Mio reads per sample. Detailed information on quality control measures, the assembly of the de novo transcriptome and the annotation can be found in Eberl et al. [12], but the most relevant information will be summarized here. The annotation was done using, among others, BLAST, Gene Ontology (GO) and InterPro terms (InterProScan, EBI). Contigs encoding for potential KTI proteins were identified based on a positive BLAST hit against a known KTI in the NCBI nr database, GO terms associated with serine proteinase inhibitors and/or a hit against the Pfam domain PF00197 (Kunitz STI protease inhibitor), or InterPro domains IPR011065 (Kunitz inhibitor STI-like superfamily) and IPR002160 (Proteinase inhibitor I3, Kunitz legume). In order to identify further KTI candidates, the P. nigra transcriptome was uploaded in an internal database and used for BLAST analysis of poplar KTI sequences from NCBI (www.ncbi.nlm.nih.gov/) and Phytozome (https://phytozome.jgi.doe.gov/pz/portal.html). Digital gene expression analysis was carried out using CLC Genomics Workbench v9.1 to generate BAM (mapping) files, and expression levels were then estimated using QSeq Software (DNAStar Inc., Madison, WI, United States). The log2 (RPKM) values (normalized mapped read values; geometric means of the biological replicate samples) were used to calculate fold-change values. Differentially expressed genes were identified using the Student’s t-test (as implemented in Qseq) corrected for multiple testing using the Benjamini–Hochberg procedure to check the false discovery rate (FDR). With an FDR-corrected p-value less than 0.05 a gene was considered significantly differentially expressed.

In addition to the KTI gene sequences in the transcriptome of the herein described experiment, another KTI gene (PnKTI B1) was identified from an additional leaf transcriptome from the same P. nigra genotype and comparable L. dispar herbivory treatment (unpublished). Furthermore, another sequence encoding a KTI (PnKTI A4, or SQ33325–2), which was not present in the transcriptome, was identified during amplification from cDNA (see below) with primers originally designed for PnKTI A13 (SQ33325).

Cloning and sequencing of PI genes

Full-length open reading frames (ORF) were amplified from a mix of cDNA originating from herbivore-induced samples in a PCR using Phusion High Fidelity polymerase in HF-buffer according to the manufacturer’s manual (New England Biolabs GmbH, Frankfurt/Main, Germany). Primers were designed based on the putative ORF from the transcriptome whenever available, or with the ORF of the homologous genes retrieved from the NCBI data base (https://www.ncbi.nlm.nih.gov/). PCR products were cloned into a PCR4-blunt TOPO vector (Thermo Fisher) and fully sequenced using the Sanger protocol and capillary sequencing with an ABI Prism-Gene- Analyser 3130xl (Applied Biosystems).

Sequence alignments and phylogenetic analysis

Homologs of P. nigra KTI sequences were identified using the BLAST-search of the NCBI data base (https://blast.ncbi.nlm.nih.gov/Blast.cgi) and the P. trichocarpa genome v3.0 (https://phytozome.jgi.doe.gov/). Alignments and similarity calculations were done with Geneious software (Biomatters, Auckland, New Zealand).

An amino acid alignment of poplar KTI proteins was constructed using the MUSCLE algorithm implemented in MEGA6 [46]. Tree reconstruction was done with MEGA6 using the Neighbor-Joining method and the JTT matrix-based method. All positions with less than 80% site coverage were eliminated.

Gene expression analysis by qRT-PCR

cDNA (diluted 1:3 with water) from the Gene expression samples was used for quantitative real-time PCR (qRT-PCR), which was performed in a Brilliant III Ultra-Fast SYBR reaction mixture (Agilent) on a CFX Connect Real-Time PCR Detection System (Bio-Rad Laboratories, Hercules, CA, USA) with 40 2-step cycles (95 °C, 30s + 60 °C, 30s) and a melting curve from 53 to 95 °C. Primer sequences can be found in Table S2. The PCR products were verified by cloning and sequencing as described above. Gene expression was calculated using CFX Manager 3.1 (Bio-Rad) using the ΔΔcq method and taking primer efficiencies into account. Values were normalized to Actin as a reference gene [41] and expressed relative to a control sample.

Trypsin-inhibiting activity assay

In order to correlate gene expression with protease inhibitor activity, the trypsin-inhibiting activity assay was performed as described in Fabisch et al. [14]. In short, 10 mg freeze-dried leaf material was extracted with 400 μL buffer (25 mM Hepes-KOH, pH 7.2, 3% PVPP, 2% PVP, 1 mM EDTA) and the extract tested for trypsin-inhibiting activity in a colorimetric (cleavage of N-acetyl-DL-phenylalanine beta-naphthyl ester) in-gel diffusion assay.

Statistical analysis

All data were checked for statistical assumptions, i.e. homogeneity of variances and normal distribution. Gene expression data for all KTI genes had to be log10-transformed to meet the statistical assumptions for parametric testing. For gene expression data, a one-way MANOVA (multivariate analysis of variance) coupled to a Tukey’s post-hoc test was applied. All statistical analyzes were conducted using SPSS 17.0 (SPSS, Chicago, IL, USA).

Results

Identification of herbivore-induced Kunitz-type trypsin inhibitors

The transcriptome of black poplar leaves with and without herbivory by two different insect species, Lymantria dispar (Lepidoptera) and Phratora vulgatissima (Coleoptera), was used to identify genes encoding herbivore-induced Kunitz-type trypsin inhibitors (KTIs). Among all sequences in the transcriptome, 45 were identified as protease inhibitor genes (PIs), of which 30 were up-regulated upon both caterpillar and beetle herbivory, seven showed different regulation patterns depending on herbivore identity, and eight were down-regulated upon herbivory by either of the herbivores (Table S3). Among the 45 PI genes, 15 belong to the KTIs, and were all up-regulated upon herbivory (Fig. 2). These 15 KTI sequences, plus two additionally identified KTI genes, were compared to previously described poplar KTIs (Table S4) and named according to the nomenclature of Ma et al. [27].

Fig. 2
figure 2

Expression of Kunitz-type trypsin inhibitor (KTI) genes in black poplar leaves after feeding by gypsy moth caterpillars (L. dispar) or blue willow beetles (P. vulgatissima) compared to their respective controls (Control 1 and 2), and compared to actin (ACT) and elongation factor 1-α (EF1-α) as constitutively expressed ‘house-keeping genes’. Shown are the mean RPKM (reads per kilobase of transcript per million mapped reads; n = 4) as result of the transcriptome analysis

Full size image

A phylogenetic analysis based on the amino acid alignment revealed that the KTIs cluster into 4 subfamilies (Fig. 3). Most of the 17 KTIs belong to the subfamilies A and C, whereas only one protein belongs to subfamily D. Interestingly, all members of the C-subfamily showed a low expression and were only marginally up-regulated upon herbivory in comparison to members of the other three subfamilies (Fig. 2). Therefore, KTIs from the subfamily C were not considered in further analysis. Out of the remaining KTI genes, those with the highest expression levels in herbivore-induced samples were chosen for cDNA sequencing, yielding the full-length open reading frames of ten PnKTI genes (Fig. 3).

Fig. 3
figure 3

Phylogenetic tree of poplar KTI proteins. PnKTIs identified in this study are shown in blue and asterisks mark full-length cDNA sequences. The tree was inferred using the Neighbor-Joining method and the JTT matrix-based method. Bootstrap values (n = 500) are shown next to each node. The tree is drawn to scale, with branch lengths in the units of the number of amino acid substitutions per site

Full size image

Herbivore-specific induction of KTI gene expression

To study the specificity of KTI gene expression, we used three different herbivore species that exhibit either similar (L. dispar, Amata mogadorensis) or different (P. vulgatissima) damage patterns on black poplar leaves (Fig. 1), but all cause similar leaf area loss (Table S1). In a previous study, we showed that total trypsin inhibitor activity in black poplar leaves is induced upon herbivory by three different herbivores, especially by P. vulgatissima [14]. To study this phenomenon at the transcriptional level, the relative gene expression of nine candidate PnKTIs was analyzed by qRT-PCR, using randomly selected samples from this previous study.

While PnKTI A2 could not be amplified in the qPCR reaction and was therefore excluded from further analysis, all of the remaining eight PnKTI genes showed significant up-regulation upon herbivory by all of the tested insects (Fig. 4). A multivariate analysis including damage severity as covariate revealed that the herbivory treatment had the strongest effect on PnKTI gene expression (F(24) = 7.230; P < 0.001).

Fig. 4
figure 4

Transcript accumulation of Kunitz-type trypsin inhibitor genes (KTIs) of the a, b, and d subfamily in black poplar leaves after herbivory by two caterpillar species (L. dispar, A. mogadorensis) and one beetle species (P. vulgatissima). Shown are the gene expression normalized to Actin and relative to a control sample as boxplots (median with upper and lower quartile as bars; n = 6); results of the ANOVA are given in each graph. Different letters indicate significant differences among groups (P < 0.05; Tukey’s post-hoc test)

Full size image

Constitutive expression levels in undamaged leaves differed among the PnKTI genes, with members of the A subfamily generally displaying higher expression levels than those of the B and D subfamilies (Table S5). Upon herbivory, however, the genes showed even more apparent differences in their inducibility (Fig. 4). Caterpillar herbivory by L. dispar and A. mogadorensis resulted in an up-regulation of all KTI genes by approximately 10 (PnKTI A6) to 2000-fold (PnKTI D2) in comparison to the constitutive levels. All KTI genes were induced to similar levels by these two lepidopteran herbivores. Beetle herbivory by P. vulgatissima, however, caused a much stronger induction of KTI gene expression than caterpillar herbivory. Expression levels in beetle-damaged leaves increased up to 40,000-fold (PnKTI D2) compared to undamaged controls. Nevertheless, the induction levels differed substantially among the individual genes, ranging from approximately 40 (PnKTI A6) and several hundred (PnKTI A15, B1, B5) up to several thousand-fold (PnKTI A13, A14, D2). Interestingly, PnKTI D2, a gene with one of the lowest constitutive expression levels, showed by far the strongest relative induction upon both caterpillar and beetle herbivory (Fig. 4h).

When considering herbivore treatments only (excluding the undamaged control group), we found that the damage severity (% leaf area loss) did not have a substantial effect on expression levels of most of the PnKTIs. Only two genes, PnKTI A7 and PnKTI B1, were significantly influenced by this factor in their expression (ANCOVA; PnKTI A7: F(1) = 9.348, P = 0.009; PnKTI B1: F(1) = 5.012; P = 0.042). Accordingly, the expression of the PnKTIs also did not correlate with the damage severity, except for PnKTI A7, which showed a positive relationship with leaf area loss (Spearman’s rank correlation: ρ = 0.556, P = 0.017). The total trypsin-inhibiting activity (Table S1 [14];), on the other hand, strongly correlated with the expression of all PnKTIs with a positive relationship (Table 1).

Table 1 Correlations of individual PnKTI gene expression versus total foliar trypsin-inhibiting activity (μg g− 1 DW; data from [14]) in all herbivore-treated (L. dispar, A. mogadorensis, and P. vulgatissima feeding) samples of black poplar leaves. Spearman rank-correlation, significant values are highlighted in bold font
Full size table

Discussion

Here we describe sequence analyses and expression patterns of Kunitz-type trypsin inhibitors (KTIs) in black poplar (Populus nigra), including ten full-length cDNA sequences, of which six had not been described before in P. nigra. Eight of these PnKTIs were studied in the context of herbivore species-specific induction patterns in leaves and we could show that beetle herbivory elicits a much stronger transcriptional response than caterpillar herbivory of the same magnitude.

Expression levels and inducibility of individual black poplar KTI genes

The up-regulation of protease inhibitor (PI) transcription and activity has been described previously, also in black poplar [27, 38]. However, both constitutive expression levels and amplitude of induction vary between studies. In our study, members of the C-subfamily generally showed low expression levels and little or no up-regulation upon herbivory (Fig. 2). In contrast, Ma et al. [27] observed stronger herbivore induction for most of the genes in this subfamily. This suggests that the regulation of KTI transcription depends on more factors than herbivore feeding or wounding alone. Certain traits of the plants, such as age, genotype [44] or previously experienced damage may play a role, but also the experimental conditions such as abiotic conditions, timing [27, 39] or damage severity could potentially influence expression levels. However, there are also consistent patterns among the different studies. In our study, PnKTI D2, the only member of the D subfamily, showed the highest inducibility, i.e. relative change upon herbivory (Fig. 4). The same gene was amongst the most up-regulated KTIs upon herbivory and mechanical wounding in another black poplar study [27]. Similarly, the high herbivore-induced expression levels of PnKTI A14 in our experiments (Fig. 2; Table S5) match well with the results obtained for the corresponding ortholog in a hybrid poplar species (P. trichocarpa x deltoides) after herbivory and mechanical wounding [39]. However, this gene also showed relatively high transcript abundance in undamaged controls, assuming also a role in constitutive defense or primary metabolism. On the contrary, PnKTI D2, which displays minimal expression levels in undamaged tissue in our and a similar study [27], seems to act exclusively in induced anti-herbivore defense.

There was no correlation between gene expression for most of the PnKTI genes and damage severity, which suggests a threshold-based activation of PnKTI transcription rather than continuous control, in which more damage would lead to higher KTI transcript levels. Furthermore, we found a strong positive relationship between the trypsin-inhibiting activity in poplar leaves and the transcription levels for all PnKTI genes. This indicates that P. nigra KTI activity is predominantly controlled at the transcriptional level and hence by de novo biosynthesis. The importance of de novo biosynthesis of stress-induced PIs has already been demonstrated in rice [40].

Herbivore specificity in PnKTI induction

When we analyzed the transcription of KTIs in leaves damaged by different insect herbivores, it became evident that beetles elicited a much stronger induction of all tested KTIs than caterpillars (Fig. 4). Similar observations come from pine trees [35] and milkweed [1, 48], where beetle herbivory induced stronger defense responses (resins and terpenes, or latex, respectively) compared to caterpillar herbivory. Species-specificity has been reported for the induction of PIs in other systems, though not in poplar trees. In soybean, damage by fall armyworm caterpillars increased the activity of PIs, whereas thrips damage did not [43]. De Oliveira et al. [11] even observed varying response of tomato PIs to damage by herbivores of the same genus. They showed that PI activity was induced by the spider mite Tetranychus urticae, but was suppressed by T. evansi [11]. Interestingly, feeding damage by lepidopteran and coleopteran herbivores in tomato yielded opposite results to our study in black poplar. Here, gene expression and trypsin inhibiting activity was more strongly induced by the tobacco hornworm than by the Colorado potato beetle [10].

The difference in PnKTI expression between beetle (Phratora vulgatissima) and caterpillar (Lymantria dispar and Amata mogadorensis) herbivory might be based on the different damage pattern these insects cause, even though all three of them are leaf chewers and removed the same total leaf area. While caterpillars removed large chunks of the leaves, the beetles caused small but numerous lesions in the leaves (Fig. 1). The number of lesions was found to be a key factor determining the emission of volatiles, another important anti-herbivore defense trait in black poplar [31]. Other factors, such as the duration of damage or the chemical compounds deposited on the plant may also be important. When artificial damage was administered to lima bean with a mechanical caterpillar, changes in the amount of time that damage lasted as well as the area damaged affected the emission of volatiles [34]. Furthermore, species-specific compounds in the saliva could trigger distinct defense responses or the magnitude of response as reported here. The importance of insect-derived elicitors for PI induction has been demonstrated in another poplar species, where mechanical wounding and simultaneous application of oral secretions from forest tent caterpillars suppressed the induction of PIs [39]. It is likely that oral secretions of the insects used in this study also exhibit a suppressive effect, maybe with varying efficacy on PI induction. Whether herbivore host range, comparing generalists such as L. dispar and A. mogadorensis versus specialists such as P. vulgatissima, plays a role in the induction of PIs, is not clear. Specialists usually possess a higher tolerance towards specific chemical defenses of their hosts, such as salicinoids in black poplar trees [3]. An increased induction of a defense, such as the PIs, to specialist herbivores could therefore be a more effective way to defend against these insects. Future studies using more herbivore species, or generalists and specialists that are more closely related to each other and cause similar feeding patterns, are necessary to determine if herbivore host range influences PI induction.

Whether the herbivore specific induction patterns of PnKTIs have ecological relevance is another open question. One factor that plays an important role in this context is ‘effect specificity’ [20]. PIs possess varying effectiveness in defense against different herbivores, as could be observed in the performance of five different herbivores that had been reared on PI-supplemented diets [8]. Similarly, the cotton bollworm exhibited distinct preference and performance towards different classes of protease inhibitors [25]. This can be explained by the fact that PIs, on the one hand, vary in their ability to inhibit different proteases, i.e. trypsin, chymotrypsin and elastase [29], and that insects, on the other hand, vary in their gut protease activities [8, 20]. Additionally, the gut pH, which differs substantially between Lepidoptera and Coleoptera [20], also influences the inhibitory activity of PIs [49]. It would therefore be interesting to dissect the role of individual KTIs in black poplar towards different insect herbivores, for example by using transgenic trees or diet supplementation of recombinant KTIs. ‘Response specificity’ towards herbivore species is believed to be more cost-effective for a plant than a similar response to all herbivores [20]. Keeping in mind the fitness costs that are linked to the biosynthesis of PIs [18], a plant might aim to induce a subset of PIs to which a herbivore is most sensitive. In this context, PI activity should not be evaluated independently of other plant defense compounds. In tobacco, PIs function synergistically with the chemical defense compound nicotine, which becomes more toxic when herbivores have to compensate for nutritional deficits by increased feeding activity [45]. Black poplar contains toxic defense compounds called salicinoids, which have been shown to negatively influence herbivore performance and survival [4, 5]. Therefore, possible synergistic effects between salicinoids and PIs, the two main components of direct defense in this tree should be investigated in future studies.

Conclusion

Our major conclusion is that PI induction in black poplar leaves depends on the identity of the feeding herbivore, with beetles inducing a stronger response than caterpillars. Furthermore, PI activity is regulated at the level of transcription and most likely in a threshold-based fashion. However, most of the molecular mechanisms underlying the patterns observed and their ecological consequences remain to be elucidated.

Availability of data and materials

The short-read data have been deposited in the EBI short read archive (SRA) with the following sample accession numbers: ERS5844847- ERS5844862. The complete study can also be accessed directly using the following URL: http://www.ebi.ac.uk/ena/data/view/PRJEB43369.

Abbreviations

KTI:

Kunitz-type trypsin inhibitors

qRT-PCR:

Quantitative real-time PCR

cDNA:

Complementary DNA

PI:

Protease inhibitor

HAMPs:

Herbivore-associated molecular patterns

DAMPs:

Damage-associated molecular patterns

ORF:

Open reading frame

PCR:

Polymerase chain reaction

MANOVA:

Multivariate analysis of variance

References

  1. Agrawal AA, Hastings AP, Patrick ET, Knight AC. Specificity of herbivore-induced hormonal signaling and defensive traits in five closely related milkweeds (Asclepias spp.). J Chem Ecol. 2014;40(7):717–29. https://doi.org/10.1007/s10886-014-0449-6.

    Article  CAS  PubMed  Google Scholar 

  2. Arnaiz A, Talavera-Mateo L, Gonzalez-Melendi P, Martinez M, Diaz I, Santamaria ME. Arabidopsis Kunitz trypsin inhibitors in defense against spider mites. Front Plant Sci. 2018;9:986. https://doi.org/10.3389/fpls.2018.00986.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Boeckler GA, Gershenzon J, Unsicker SB. Phenolic glycosides of the Salicaceae and their role as anti-herbivore defenses. Phytochemistry. 2011;72(13):1497–509. https://doi.org/10.1016/j.phytochem.2011.01.038.

    Article  CAS  PubMed  Google Scholar 

  4. Boeckler GA, Paetz C, Feibicke P, Gershenzon J, Unsicker SB. Metabolism of poplar salicinoids by the generalist herbivore Lymantria dispar (Lepidoptera). Insect Biochem Mol Biol. 2016;78:39–49. https://doi.org/10.1016/j.ibmb.2016.08.001.

    Article  CAS  PubMed  Google Scholar 

  5. Boeckler GA, Towns M, Unsicker SB, Mellway RD, Yip L, Hilke I, et al. Transgenic upregulation of the condensed tannin pathway in poplar leads to a dramatic shift in leaf palatability for two tree-feeding Lepidoptera. J Chem Ecol. 2014;40(2):150–8. https://doi.org/10.1007/s10886-014-0383-7.

  6. Botelho-Júnior S, Machado OL, Fernandes KV, Lemos FJ, Perdizio VA, Oliveira AE, et al. Defense response in non-genomic model species: methyl jasmonate exposure reveals the passion fruit leaves’ ability to assemble a cocktail of functionally diversified Kunitz-type trypsin inhibitors and recruit two of them against papain. Planta. 2014;240(2):345–56. https://doi.org/10.1007/s00425-014-2085-3.

    Article  CAS  PubMed  Google Scholar 

  7. Bradshaw HD, Hollick JB, Parsons TJ, Clarke HR, Gordon MP. Systemically wound-responsive genes in poplar trees encode proteins similar to sweet potato sporamins and legume Kunitz trypsin inhibitors. Plant Mol Biol. 1990;14(1):51–9. https://doi.org/10.1007/BF00015654.

    Article  CAS  PubMed  Google Scholar 

  8. Broadway RM. Are insects resistant to plant proteinase inhibitors? J Insect Physiol. 1995;41(2):107–16. https://doi.org/10.1016/0022-1910(94)00101-L.

    Article  CAS  Google Scholar 

  9. Christopher ME, Miranda M, Major IT, Constabel CP. Gene expression profiling of systemically wound-induced defenses in hybrid poplar. Planta. 2004;219(6):936–47. https://doi.org/10.1007/s00425-004-1297-3.

    Article  CAS  PubMed  Google Scholar 

  10. Chung SH, Felton GW. "Specificity of induced resistance in tomato against specialist lepidopteran and coleopteran species." J Chem Ecol. 2011;37(4):378–86.

  11. De Oliveira EF, Pallini A, Janssen A. Herbivore performance and plant defense after sequential attacks by inducing and suppressing herbivores. Insect Sci. 2019;26(1):108–18. https://doi.org/10.1111/1744-7917.12499.

    Article  CAS  PubMed  Google Scholar 

  12. Eberl F, Perreca E, Vogel H, Wright LP, Hammerbacher A, Veit D, et al. Rust infection of black poplar trees reduces photosynthesis but does not affect isoprene biosynthesis or emission. Front Plant Sci. 2018;9:1733. https://doi.org/10.3389/fpls.2018.01733.

  13. Erb M, Meldau S, Howe GA. Role of phytohormones in insect-specific plant reactions. Trends Plant Sci. 2012;17(5):250–9. https://doi.org/10.1016/j.tplants.2012.01.003.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Fabisch T, Gershenzon J, Unsicker SB. Specificity of herbivore defense responses in a woody plant, black poplar (Populus nigra). J Chem Ecol. 2019;45(2):162–77. https://doi.org/10.1007/s10886-019-01050-y.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Fürstenberg-Hägg J, Zagrobelny M, Bak S. Plant defense against insect herbivores. Int J Mol Sci. 2013;14(5):10242–97. https://doi.org/10.3390/ijms140510242.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Garcia VA, Freire MDGM, Novello JC, Marangoni S, Macedo MLR. Trypsin inhibitor from Poecilanthe parviflora seeds: purification, characterization, and activity against pest proteases. Protein J. 2004;23(5):343–50. https://doi.org/10.1023/B:JOPC.0000032654.67733.d5.

    Article  CAS  PubMed  Google Scholar 

  17. Gatehouse AM, Boulter D. Assessment of the antimetabolic effects of trypsin inhibitors from cowpea (Vigna unguiculata) and other legumes on development of the bruchid beetle Callosobruchus maculatus. J Sci Food Agric. 1983;34(4):345–50. https://doi.org/10.1002/jsfa.2740340405.

    Article  CAS  Google Scholar 

  18. Glawe GA, Zavala JA, Kessler A, Van Dam NM, Baldwin IT. Ecological costs and benefits correlated with trypsin protease inhibitor production in Nicotiana attenuata. Ecology. 2003;84(1):79–90. https://doi.org/10.1890/0012-9658(2003)084[0079:ECABCW]2.0.CO;2.

    Article  Google Scholar 

  19. Ham WE, Sandstedt R. A proteolytic inhibiting substance in the extract from unheated soy bean meal. J Biol Chem. 1944;154(2):505–6. https://doi.org/10.1016/S0021-9258(18)71934-0.

    Article  CAS  Google Scholar 

  20. Haq SK, Atif SM, Khan RH. Protein proteinase inhibitor genes in combat against insects, pests, and pathogens: natural and engineered phytoprotection. Arch Biochem Biophys. 2004;431(1):145–59. https://doi.org/10.1016/j.abb.2004.07.022.

    Article  CAS  PubMed  Google Scholar 

  21. Heitz T, Geoffroy P, Fritig B, Legrand M. The PR-6 family: proteinase inhibitors in plant-microbe and plant-insects interactions. In: Datta SK, S. Muthukrishnan S (Eds.); Pathogenesis-Related Proteins in Plants. Boca Raton: CRC Press; 1999. p. 131–55.

  22. Jamal F, Pandey PK, Singh D, Ahmed W. A Kunitz-type serine protease inhibitor from Butea monosperma seed and its influence on developmental physiology of Helicoverpa armigera. Process Biochem. 2015;50(2):311–6. https://doi.org/10.1016/j.procbio.2014.12.003.

    Article  CAS  Google Scholar 

  23. Karban R, Baldwin IT. Induced responses to herbivory: University of Chicago Press; 2007.

  24. Kunitz M. Crystallization of a trypsin inhibitor from soybean. Science. 1945;101(2635):668–9. https://doi.org/10.1126/science.101.2635.668.

    Article  CAS  PubMed  Google Scholar 

  25. Kuwar SS, Pauchet Y, Heckel DG. Effects of class-specific, synthetic, and natural proteinase inhibitors on life-history traits of the cotton bollworm Helicoverpa armigera. Arch Insect Biochem Physiol. 2020;103:e21647.

    Article  CAS  PubMed  Google Scholar 

  26. Lee SI, Lee S-H, Koo JC, Chun HJ, Lim CO, Mun JH, et al. Soybean Kunitz trypsin inhibitor (SKTI) confers resistance to the brown planthopper (Nilaparvata lugens Stål) in transgenic rice. Mol Breed. 1999;5(1):1–9. https://doi.org/10.1023/A:1009660712382.

    Article  Google Scholar 

  27. Ma Y, Zhao Q, Lu M-Z, Wang J. Kunitz-type trypsin inhibitor gene family in Arabidopsis and Populus trichocarpa and its expression response to wounding and herbivore in Populus nigra. Tree Genet Genomes. 2011;7(2):431–41. https://doi.org/10.1007/s11295-010-0345-3.

    Article  Google Scholar 

  28. Major IT, Constabel CP. Molecular analysis of poplar defense against herbivory: comparison of wound-and insect elicitor-induced gene expression. New Phytol. 2006;172(4):617–35. https://doi.org/10.1111/j.1469-8137.2006.01877.x.

    Article  CAS  PubMed  Google Scholar 

  29. Major IT, Constabel CP. Functional analysis of the Kunitz trypsin inhibitor family in poplar reveals biochemical diversity and multiplicity in defense against herbivores. Plant Physiol. 2008;146(3):888–903. https://doi.org/10.1104/pp.107.106229.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Marchetti S, Delledonne M, Fogher C, Chiaba C, Chiesa F, Savazzini F, et al. Soybean Kunitz, C-II and PI-IV inhibitor genes confer different levels of insect resistance to tobacco and potato transgenic plants. Theor Appl Genet. 2000;101(4):519–26. https://doi.org/10.1007/s001220051511.

  31. McCormick AC, Boeckler GA, Köllner TG, Gershenzon J, Unsicker SB. The timing of herbivore-induced volatile emission in black poplar (Populus nigra) and the influence of herbivore age and identity affect the value of individual volatiles as cues for herbivore enemies. BMC Plant Biol. 2014;14(1):304. https://doi.org/10.1186/s12870-014-0304-5.

    Article  CAS  PubMed  Google Scholar 

  32. MEROPS database, Rawlings ND, Barrett AJ, Thomas PD, Huang X, Bateman A, et al. The MEROPS database of proteolytic enzymes, their substrates and inhibitors in 2017 and a comparison with peptidases in the PANTHER database. Nucleic Acids Res. 2017;46:D624–32 https://www.ebi.ac.uk/merops/.

  33. Mithöfer A, Boland W. Recognition of herbivory-associated molecular patterns. Plant Physiol. 2008;146(3):825–31. https://doi.org/10.1104/pp.107.113118.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Mithöfer A, Wanner G, Boland W. Effects of feeding Spodoptera littoralis on lima bean leaves. II. Continuous mechanical wounding resembling insect feeding is sufficient to elicit herbivory-related volatile emission. Plant Physiol. 2005;137(3):1160–8. https://doi.org/10.1104/pp.104.054460.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Moreira X, Lundborg L, Zas R, Carrillo-Gavilán A, Borg-Karlson A-K, Sampedro L. Inducibility of chemical defences by two chewing insect herbivores in pine trees is specific to targeted plant tissue, particular herbivore and defensive trait. Phytochemistry. 2013;94:113–22. https://doi.org/10.1016/j.phytochem.2013.05.008.

    Article  CAS  PubMed  Google Scholar 

  36. Negreiros ANM, Carvalho MM, Xavier Filho J, Blanco-Labra A, Shewry PR, Richardson M. The complete amino acid sequence of the major Kunitz trypsin inhibitor from the seeds of Prosopsis juliflora. Phytochemistry. 1991;30(9):2829–33. https://doi.org/10.1016/S0031-9422(00)98207-4.

    Article  CAS  PubMed  Google Scholar 

  37. Neiman M, Olson MS, Tiffin P. Selective histories of poplar protease inhibitors: elevated polymorphism, purifying selection, and positive selection driving divergence of recent duplicates. New Phytol. 2009;183(3):740–50. https://doi.org/10.1111/j.1469-8137.2009.02936.x.

    Article  CAS  PubMed  Google Scholar 

  38. Nishiguchi M, Yoshida K, Sumizono T, Tazaki K. A receptor-like protein kinase with a lectin-like domain from lombardy poplar: gene expression in response to wounding and characterization of phosphorylation activity. Mol Gen Genomics. 2002;267(4):506–14. https://doi.org/10.1007/s00438-002-0683-4.

    Article  CAS  Google Scholar 

  39. Philippe RN, Ralph SG, Külheim C, Jancsik SI, Bohlmann J. Poplar defense against insects: genome analysis, full-length cDNA cloning, and transcriptome and protein analysis of the poplar Kunitz-type protease inhibitor family. New Phytol. 2009;184(4):865–84. https://doi.org/10.1111/j.1469-8137.2009.03028.x.

    Article  CAS  PubMed  Google Scholar 

  40. Rakwal R, Agrawal GK, Jwa N-S. Characterization of a rice (Oryza sativa L.) Bowman–Birk proteinase inhibitor: tightly light regulated induction in response to cut, jasmonic acid, ethylene and protein phosphatase 2A inhibitors. Gene. 2001;263(1-2):189–98. https://doi.org/10.1016/S0378-1119(00)00573-4.

  41. Ramírez-Carvajal GA, Morse AM, Davis JM. Transcript profiles of the cytokinin response regulator gene family in Populus imply diverse roles in plant development. New Phytol. 2008;177(1):77–89. https://doi.org/10.1111/j.1469-8137.2007.02240.x.

    Article  CAS  PubMed  Google Scholar 

  42. Renko M, Sabotič J, Turk D. β-Trefoil inhibitors–from the work of Kunitz onward. Biol Chem. 2012;393(10):1043–54. https://doi.org/10.1515/hsz-2012-0159.

    Article  CAS  PubMed  Google Scholar 

  43. Romero B, Dillon FM, Zavala JA. Different soybean cultivars respond differentially to damage in a herbivore-specific manner and decreas herbivore performance. Arthropod Plant Interact. 2020;14(1):89–99. https://doi.org/10.1007/s11829-019-09730-y.

    Article  Google Scholar 

  44. Rubert-Nason KF, Couture JJ, Major IT, Constabel CP, Lindroth RL. Influence of genotype, environment, and gypsy moth herbivory on local and systemic chemical defenses in trembling aspen (Populus tremuloides). J Chem Ecol. 2015;41(7):651–61. https://doi.org/10.1007/s10886-015-0600-z.

    Article  CAS  PubMed  Google Scholar 

  45. Steppuhn A, Baldwin IT. Resistance management in a native plant: nicotine prevents herbivores from compensating for plant protease inhibitors. Ecol Lett. 2007;10(6):499–511. https://doi.org/10.1111/j.1461-0248.2007.01045.x.

    Article  PubMed  Google Scholar 

  46. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9. https://doi.org/10.1093/molbev/mst197.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Unsicker SB, Gershenzon J, Köllner TG. Beetle feeding induces a different volatile emission pattern from black poplar foliage than caterpillar herbivory. Plant Signal Behav. 2015;10(3):e987522. https://doi.org/10.4161/15592324.2014.987522.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Van Zandt PA, Agrawal AA. Specificity of induced plant responses to specialist herbivores of the common milkweed Asclepias syriaca. Oikos. 2004;104(2):401–9. https://doi.org/10.1111/j.0030-1299.2004.12964.x.

    Article  Google Scholar 

  49. Wolfson JL, Murdock LL. Diversity in digestive proteinase activity among insects. J Chem Ecol. 1990;16(4):1089–102. https://doi.org/10.1007/BF01021013.

Download references

Acknowledgements

We thank Almuth Hammerbacher and Chhana Ullah for assistance during the transcriptome experiment, Bettina Raguschke for the plasmid sequencing, and Beate Rothe for her help in sample analysis.

Funding

This research was funded by the Max Planck Society. Thomas Fabisch was funded by the Gerhard and Ellen Zeidler Stiftung. The funding bodies provided financial support to the research project, but were not involved in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript. Open Access funding enabled and organized by Projekt DEAL.

Author information

Authors and Affiliations

  1. Department of Biochemistry, Max Planck Institute for Chemical Ecology (MPI-CE), Hans-Knöll-Str. 8, 07745, Jena, Germany

    Franziska Eberl, Thomas Fabisch, Katrin Luck, Tobias G. Köllner, Jonathan Gershenzon & Sybille B. Unsicker

  2. Department of Entomology, MPI-CE, Hans-Knöll-Str. 8, 07745, Jena, Germany

    Heiko Vogel

Authors
  1. Franziska Eberl
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Thomas Fabisch
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Katrin Luck
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Tobias G. Köllner
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Heiko Vogel
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Jonathan Gershenzon
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Sybille B. Unsicker
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

FE, TF, JG, TGK, SBU conceived the project; TF conducted the biological experiments and leaf sample collection; FE analyzed the gene expression; KL and FE cloned and sequenced all candidate genes; TGK created the phylogenetic tree; HV assembled and analyzed the transcriptome; FE wrote the manuscript; all authors reviewed and commented on the manuscript. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Franziska Eberl.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1

. Feeding damage by the three herbivores and trypsin-inhibiting activity in poplar. Table S2. Primer sequences used for cloning and qRT-PCR. Table S3. Differential expression of contigs annotated as protease inhibitors in the transcriptome of black poplar leaves. Table S4. Nomenclature of KTI homologs in this and other studies. Table S5. Quantification cycles of the qRT-PCR analysis for individual KTI genes.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Eberl, F., Fabisch, T., Luck, K. et al. Poplar protease inhibitor expression differs in an herbivore specific manner. BMC Plant Biol 21, 170 (2021). https://doi.org/10.1186/s12870-021-02936-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-021-02936-4

Keywords

BMC Plant Biology

ISSN: 1471-2229

Contact us