Skip to main content

Advertisement

  • Research article
  • Open Access
  • The biochemical composition and transcriptome of cotyledons from Brassica napus lines expressing the AtGL3 transcription factor and exhibiting reduced flea beetle feeding

    BMC Plant BiologyBMC series – open, inclusive and trusted201818:64

    https://doi.org/10.1186/s12870-018-1277-6

    • Received: 25 May 2017
    • Accepted: 29 March 2018
    • Published:

    Abstract

    Background

    Previously, transgenic trichome-bearing (hairy leaf) Brassica napus lines expressing either the Arabidopsis thaliana GL3 gene (line AtGL3+) [1] or the AtGL3 gene in combination with an RNAi construct to down-regulate TTG1 (line K-5-8) [2] were developed. The leaves of these lines exhibited altered insect feeding (flea beetle) and oviposition (diamondback moth) behaviour compared to the non-transgenic semi-glabrous leaves of B. napus cv. Westar. Interestingly, the cotyledons of these lines remained glabrous, but also showed reduced feeding by flea beetles. Here we examine the composition and global transcriptome of the glabrous cotyledons from these transgenic lines to ascertain the mechanism(s) underlying this unexpected phenomenon.

    Results

    Approximately, 7500 genes were up-regulated in cotyledons of each hairy line, compared with < 30 that were down-regulated. The up-regulated genes included those involved in cell wall synthesis, secondary metabolite production, redox, stress and hormone-related responses that have the potential to impact host plant cues required to elicit defense responses toward insect pests. In particular, the expression of glucosinolate biosynthetic and degradation genes were substantially altered in the glabrous cotyledons of the two hairy leaf lines. The transcriptomic data was supported by glucosinolate and cell wall composition profiles of the cotyledons. Changes in gene expression were much more extreme in the AtGL3+ line compared with the K-5-8 line in terms of diversity and intensity.

    Conclusions

    The study provides a roadmap for the isolation and identification of insect resistance compounds and proteins in the glabrous cotyledons of these hairy leaf lines. It also confirms the impact of mis-expression of GL3 and TTG1 on types of metabolism other than those associated with trichomes. Finally, the large number of up-regulated genes encoding heat shock proteins, PR proteins, protease inhibitors, glucosinolate synthesis/breakdown factors, abiotic stress factors, redox proteins, transcription factors, and proteins required for auxin metabolism also suggest that these cotyledons are now primed for resistance to other forms of biotic and abiotic stress.

    Keywords

    • Brassica napus
    • Trichomes
    • Cotyledons
    • Glucosinolates
    • Flea beetle
    • RNA sequencing

    Background

    Plants possess a variety of biochemical and morphological defences aimed at deterring insects from herbivory and oviposition. Defences can be either induced upon attack and localized to the site of attack, be systemic, or both (reviewed in [3]). Host plant resistance is categorized into three groups: i) antibiosis, resulting in increased mortality, reduced longevity or reduced reproduction of the insect; ii) antixenosis, affecting the behaviour of the insect and often expressed as non-preference for a resistant plant compared to a susceptible plant; and iii) tolerance, this being the ability of a plant to withstand or recover from insect damage and perform better than a susceptible plant grown under similar conditions [4]. Crops with insect resistance can reduce the accumulation of harmful chemical residues in the environment, as well as provide economic benefits to farmers and seed producers.

    Glabrous (smooth) and semi-glabrous (few trichomes) lines of Brassica napus L. and Brassica rapa L. canola (oilseed rape) and Brassica vegetable crops are susceptible to many specialist and generalist insect pests. Flea beetles (FBs) [Phyllotreta cruciferae (Goeze) and Phyllotreta striolata (Fab.)] are specialist Brassica pests in several parts of the world, including Canada [5], India, and Eastern Europe [6, 7]. FBs attack at the crucial seedling stage as well as at more mature stages, such as leaves and developing green pods, where they reduce seed yield and grade. FB have developed resistance to insecticides used for their control (reviewed in [8, 9]). This highlights the importance of developing insect-resistant cultivars to reduce the use of chemical protection methods.

    Trichomes (leaf hairs) have evolved as a physical defence against herbivore feeding and oviposition [10], and trichome density and length both negatively impact these processes in many insect species. Trichomes on Arabidopsis thaliana L. develop at the distal end of the developing leaf, thereby protecting the more supple, younger parts of the leaf that are preferred for diamondback moth oviposition and more vulnerable to FB and Psylloides sp. feeding damage [11, 12]. Mature Brassica villosa leaves have an extremely high density of trichomes (~ 4000 cm− 1) and are immune to FB damage in the field as the insects avoid the leaves [13].

    Several Arabidopsis thaliana glabrous (lacking trichomes) mutants are also deficient in the production of secondary metabolites, most notably anthocyanins. The link between trichome formation and secondary metabolite production occurs through sharing of components within the Myb/bHLH/TTG1 (MBW) transcriptional regulatory complex. In trichome formation, the basic helix-loop-helix transcription factors GL3, or the similar factor ENHANCER OF GL3, form a regulatory complex with the R2R3-MYB factor GL1 and the WD-repeat protein TTG1 which interacts with GL3 and EGL3. This regulatory complex activates the expression of genes encoding a secondary set of transcription factors comprising GL2, TTG2 and SIM to induce trichome formation. In addition, at least 6 other R3 MYB proteins (CPC, TRY, ETC1, ETC2, ETC3 and TCL1) can replace GL1 to disrupt and/or alter the specificity of the regulatory complex [1416]. Genes encoding enzymes involved in the later stages of anthocyanin biosynthesis are regulated by a similar complex consisting of TTG1 and GL3, but with different MYB factors [17].

    Earlier, we demonstrated that ectopic expression of A. thaliana GL3 dramatically increases trichome formation on B. napus leaves, but compromises plant development [1]. Trichome production is further enhanced, while the developmental abnormalities are alleviated, when the expression of the endogenous TTG1 gene is reduced through the introduction of an RNAi construct [2]. Like hairy B. villosa, FB feeding is reduced on seedlings of the hairy leaf transgenic B. napus lines [1, 2, 18] as the insects do not initiate feeding probes on the trichome-enhanced leaves [2, 19]. The leaves of these lines exhibit large changes to their transcriptomes and growth patterns [2]. Curiously, the glabrous cotyledons from the hairy-leaf transgenic B. napus lines also showed highly reduced FB feeding (30–50%) [18]. It was suggested that the antixenotic effect could be due to altered plant architecture as the cotyledons in both transgenic lines are initially vertically oriented before becoming horizontally oriented as the plant develops. However, they could not rule out the possibility that a change in biochemical composition was also involved. Attraction and stimulation of feeding in FBs is governed principally by glucosinolates and/or their breakdown products from the host plant [20]. This suggested that the composition and gene expression patterns of the AtGL3 transgenic glabrous cotyledons may also have been altered so as to impact FB behaviour. Here, we examine the composition of glabrous cotyledons of the two transgenic hairy B. napus lines: one hairy leaf line expressing the AtGL3 gene (AtGL3+) and one ultra-hairy leaf line (K-5-8) expressing AtGL3 and down-regulated in the expression of the BnTTG1 gene. We demonstrate that the glabrous cotyledons of these lines exhibited altered secondary metabolite (anthocyanins and glucosinolates) and lignin content, as well as altered expression of genes specifying secondary metabolite biosynthesis and degradation, cell wall biosynthesis, hormones, and redox proteins which may have contributed to changes in host plant cues for insect pests.

    Methods

    Plant material

    Untreated seeds of three plant entries were used in this study: the semi-glabrous leaf B. napus cv. Westar (parent line), a homozygous hairy-leaf AtGL3 transgenic B. napus developed from a cv. Westar parental background (line AtGL3+, 1), and a homozygous ultra-hairy transgenic leaf B. napus (line K-5-8) having highly reduced BnTTG1 expression within the AtGL3+ B. napus background [2]. Seeds were sterilized and placed onto solid MS media. Cotyledons for glucosinolates (GS), cell wall component analysis, and RNA sequencing were harvested from 7-day old seedlings grown in three replicates in magenta jars (10 seeds/replicate) under insecticide-free conditions in a controlled plant growth chamber (22/18o C; 16 h photoperiod; 60–80 μmoles.m− 2.s− 1). Cotyledons for qRT-PCR of trichome genes were grown at 22/24o C under a 16 h photoperiod with light at 400 μE. m− 2.s− 1. Cotyledons used to measure phenotype changes, anthocyanin content, and qRT-PCR of anthocyanin genes were grown for 10 days under continuous light (400 μE. m− 2.s− 1).

    Cotyledon composition analysis

    Seven-day-old cotyledons were extracted and glucosinolates (GS) converted to desulfoglucosinolates (DS-GS) based on the AOCS Official Method Ak 1–92. Specifically, freeze-dried cotyledons (~ 0.1 g) were agitated with steel rods (25 × 8 mm) on an Eberbach reciprocating shaker for 10 min at 280 rpm, and then 3 ml of methanol and 1 ml of 0.2 mM benzyl GS were added and shaking continued for 60 min. After centrifugation at 2300 g for 15 min, 3 ml of supernatant was loaded onto 0.3 ml pre-swollen DEAE-Sephadex resin (~ 30 mg) in Bio-Spin micro-columns (Bio-Rad, Mississauga Canada). The resin was rinsed with 1.5 ml of 2% acetic acid, 1.8 ml of water and 1.2 ml of 20 mM sodium acetate at pH 4.0, and then 100 μl of a purified sulfatase solution was added to the resin and the micro-columns sealed and incubated at 20 °C overnight. DS-GS were eluted with 1.2 ml of water, filtered, and separated using a Waters UPLC-PDA-TQD system and a BEH Shield RP18 column (2.1 × 50 mm; 1.7 μm) (100% water at 0.8 mL/min) for 0.3 min, followed by a linear gradient of 0% to 25% acetonitrile (v/v) over 6.7 min. The DS-GS were quantified at 229 nm and identified by monitoring the characteristic loss of 162.2 mass units using MS/MS constant neutral loss scans.

    Cell wall carbohydrates and lignin were measured on purified cell wall residue (CWR) as per methods in UpdeGraff [21], Brinkman et al. [22], and Aung et al. [23]. CWR was extracted in phosphate buffer with Triton X-100 using 10 mg of seven-day-old cotyledons. Lignin content was analysed on CWR based on the thioglycolate-alkaline hydrolysis assay [22] and quantified using UV/VIS spectrophotometry at 280 nm and a calibration curve developed with commercial lignin (Sigma-Aldrich, Oakville, ON, Canada). Total acid-releasable cellulosic glucose was determined from CWR using a commercial cellulose standard (Sigma-Aldrich) and anthrone reagent spectrophotometry at 280 nm after sugar and starch removal based on the method of Theander et al. [24].

    Cotyledon RNA sequencing

    Total RNA was extracted from cotyledons of seven-day-old seedlings using an RNAeasy Mini Kit with contaminating gDNA being removed using RNAse-free DNAse™ (Qiagene Inc., ON Canada). RNA samples were quantified and RNA integrity determined using an RNA6000 nano assay in an Agilent 2100 Bioanalyzer™ (Agilent Technologies, Palo Alto, CA USA). RNA library preparation and sequencing were carried out using the Illumina TrueSeq RNA sample preparation platform v.2 with multiplex labeling following the manufacturer’s protocols. Details on cDNA library development, RNA sequencing, and data analysis were identical to those found in Alahakoon et al. [2]. Cotyledon gene expression changes in the two transgenic hairy leaf lines relative to cv. Westar cotyledon expression were then introduced into MAPMAN for allocation into 36 functional barcode index numbered (BIN) categories [25].

    Quantitative-reverse transcription PCR (qRT-PCR) analysis

    qRT-PCR was conducted on 10-day-old cotyledons to test AtGL3 expression and the summed transcripts of all homeologues for each of five B. napus regulatory genes (BnGL1, BnGL2, BnGL3, BnTTG1 and BnTRY) encoding proteins known to be involved in the trichome MYB-βHLH − WD40 tri-protein initiation complex as outlined in Alahakoon et al. [2].

    Statistical analysis

    Cotyledon composition and qRT-PCR data were analyzed with either one-way or two-way ANOVA using a MIXED model in SAS 9.2 [26] or a t-test. Assumptions of ANOVA were tested using a Normality test (Shapiro-Wilk), an Equal variance test (Brown Forsythe) and a Levenes test. Means were compared using a Tukey test or pairwise using a Dunn’s method in SAS 9.2, and treatments were declared significant at P ≤ 0.05 and trends declared at P ≤ 0.1. Read counting and statistical analysis of the RNA-seq data were carried out using Cuffdiff in the Cufflinks software package [27].

    Results

    Cotyledon composition

    As with true leaves, the glabrous AtGL3+ cotyledons were much smaller than Westar and K-5-8 cotyledons (Fig. 1 insert), reflecting the smaller stature and lower vigor of the AtGL3+ line, while the K-5-8 line produced much larger cotyledons and subsequently plants with higher vigour [2, 18]. We compared the biochemical composition of the cotyledons expressing AtGL3 to that of the wild-type line to determine if particular specific secondary metabolites or polymers contributing to cell fortification could be correlated with altered insect behaviour. Anthocyanins are particularly useful in young plants as a means of insect defense as they alter the spectral properties of the plant, potentially making it less visible, and increase their phenolic content [28]. AtGL3+ cotyledons also exhibited more red coloration on the abaxial surface than Westar cotyledons, while the K-5-8 cotyledon abaxial surface was green when grown indoors under continuous light to stress the cotyledons (Fig. 1a insert). The red coloration was due to increased anthocyanin (Fig. 1a) and an increase in transcription of three anthocyanin genes (BnANS, BnDFR and BnGST) in AtGL3+ cotyledons was noted compared to the other two lines (Fig. 1b).
    Fig. 1
    Fig. 1

    Total seedling anthocyanins and qRT-PCR of anthocyanin gene expression in B. napus cv. Westar and transgenic lines (AtGL3+ and K-5-8) grown under 24 h continuous light. Panel A: Anthocyanins. Insert shows colour and morphology on the abaxial surface of the cotyledon. Panel B: qRT-PCR of anthocyanin genes. Expression of individual genes is relative to that of glabrous B. napus cv. Westar (set at 1), which has been normalized to the expression of the B. napus Act2 gene. A Tukey test was used to detect significant differences in total anthocyanins or expressed genes between the plant lines. Means (n = 3) + standard error with different letters differ significantly (p ≤ 0.05). FW = fresh weight

    Glucosinolates are sulfur-linked glucosides commonly found in species within the Brassicaceae/Cruciferae [29]. Feeding initiation in crucifer specialists, such as FB, are influenced by glucosinolates and their breakdown products derived from the host plant [20]. The level of GSs in glabrous cotyledons of the two transgenic lines and Westar was 1000-fold lower than GS levels in the seed (data not shown). Total GS levels in K-5-8 cotyledons were not significantly different from those of Westar, while total GS in the AtGL3+ line was lower (Fig. 2a). The GS profile of the lines expressing AtGL3 was also found to be different from that of Westar and different from one another. Progoitrin (2-OH-3-butenyl-GS), 4-hydroxyglucobrassicin (4-hydroxy-3-indolylmethyl-GS) and glucobrassicin (indol-3-ylmethyl-GS) were the most abundant cotyledon GSs. While progoitrin levels in the transgenic lines were similar to Westar, the level was significantly higher in the K-5-8 line when compared to AtGL3+. The level of glucoraphanin (4-methylsulfinylbutyl-GS) was not significantly different among the lines, whereas the level of gluconapin (3-butenyl-GS) was significantly higher in the transgenic lines and further elevated in the AtGL3+ line. 4-hydroxyglucobrassicin levels were reduced with the K-5-8 lines having only 60% of that in Westar. Glucobrassicin levels in the transgenic lines were statistically similar to Westar; however, the K-5-8 cotyledons had significantly less of this GS than those of AtGL3+. Both types of transgenic cotyledons had lower 4-methoxy-glucobrassicin (4-methoxy-indol-3-ylmethyl-GS) levels compared to Westar with the AtGL3+ having even lower levels than the K-5-8 line.
    Fig. 2
    Fig. 2

    Composition of 7-day-old glabrous cotyledons. a Glucosinolates. UPLC peaks represent: 1, progoitrin (2-OH-3-butenyl-GS); 2, glucoraphanin (4-methylsulfinylbutyl-GS); 3, 3-butenyl-GS; 4, 4-hydroxy-3-indolylmethyl-GS; 5, 4-methoxy-3-indolylmethyl-GS; 6, 4-methoxy-glucobrassicin (4-methoxy-indol-3-ylmethyl-GS). Different letters indicate pairwise significance difference of the means (± SD) for each line within each compound type (p < 0.05). b Cell wall residue carbohydrates. c Cell wall residue lignin. Asterisk (*) indicates significance difference of the means (± SD) relative to cv. Wester (p < 0.05)

    Changes in the composition of the plant cell wall (e.g. cellulose, hemicellulose, pectin and lignin) serves to fortify plant tissues and can make them more difficult to consume and more difficult from which to extract nutrition [30]. Cotyledon cell wall carbohydrate was similar between the two transgenic lines and cv. Westar (Fig. 2b). Lignin extracted from the cell wall residue was similar between Westar and the AtGL3+ line, but was significantly reduced in K-5-8 cotyledons (Fig. 2c).

    Cotyledon transcriptomes related to insect host cues, metabolism and regulation

    The glabrous cotyledons of the K-5-8 and AtGL3+ lines were as resistant to FB feeding as those from pesticide-treated Westar [18]. To further explore these phenomena, we conducted RNA-Seq analysis to identify transcription patterns in cotyledons which might affect insect behaviour. Cotyledons of the AtGL3+ line and K-5-8 showed up-regulation of 7924 and 7286 genes relative to Westar cotyledons, respectively, of which 4477 were common (Fig. 3). Curiously, fewer than 30 genes were down-regulated in either the hairy line or the ultra-hairy line relative to Westar (Fig. 3; Additional file  3: Figure S1 and Additional file 4: Figure S2). These genes could be organized within 36 MAPMAN functional categories (BINS), with the largest number of up-regulated genes falling within 17 functional categories (Fig. 3; Additional file 1: Table S1). Categories with moderate (ca. 20–200) numbers of up-regulated genes relative to Westar cotyledons, included those specifying photosynthesis (BIN #1), carbohydrates (BINs #2, 3), glycolysis (BIN #4), TCA/organic acid transformation (BIN #8), cell wall synthesis (BIN #10), lipid and amino acid metabolism (BIN #11, 13), metal handling (BIN #15), secondary metabolism (BIN #16), hormones (BIN #17), stress (BIN #20), redox (BIN #21), nucleotide-related (BIN #23), DNA-related (BIN #28), signalling (BIN #30), cell organization, etc. (BIN #31), development (BIN #31), and transport (BIN #34) (Fig. 3; Additional file 3: Figure S1 and Additional file  4: Figure S2). Categories with a large number of up-regulated genes (> 500) included those without known functions (BIN #26), RNA-related (BIN #27), protein-related (BIN #29) and miscellaneous functions (BIN #35; e.g. cytochrome P450 enzymes, carbohydrases, lipases) (Fig. 3). Moreover, ~ 25% of all genes that exhibited different expression patterns in the transgenic cotyledons could not be assigned a function (BIN #35).
    Fig. 3
    Fig. 3

    Overview of changes in glabrous cotyledon gene expression in 7-day-old hairy leaf (AtGL3+) and ultra-hairy leaf (K-5-8) B. napus lines relative to cv. Westar. Upper Panel: Venn diagram showing the number of up-regulated genes. Lower Panel: Mapman categories. The graphs show different data graphed on expanded or smaller Y axes. #/# above some pairs of bars indicate the number of down-regulated genes in that bin

    Within the functional categories above, a suite of genes potentially involved in insect resistance and specifying tissue toughness (wax, cell wall carbohydrates/proteins and lignin synthesis), metal handling, flavonoid glycosylation, as well as phenylpropanoid, alkaloid, and cyanogenic glycoside synthesis, were strongly up-regulated in glabrous cotyledons of both the AtGL3+ hairy leaf line and the K-5-8 ultra-hairy leaf line relative to Westar cotyledons (Table 1, Additional file 1: Table S1; Additional file 5: Figure S3). Genes involved in GS biosynthesis and degradation were also up-regulated in both types of transgenic cotyledons, but a much larger number of these genes tended to be represented in transcriptomes of AtGL3+ cotyledons than in K-5-8 cotyledons (Fig. 4; Additional file 2: Table S2). In particular, genes encoding proteins involved in core biosynthesis of methionine-based, indole and benzyl GSs, including cytochrome P450s and a sulfotransferase (SOT18), as well as a MYB factor, were > 2-fold up-regulated in the AtGL3+ line, although genes encoding SOT18 and an aliphatic acconitase were highly up-regulated in the K-5-8 line (Fig. 4; Additional file 1: Table S1 and Additional file 2: Table S2). Moreover, a large number of genes involved in GS degradation, including those encoding myrosinase, a number of myrosinase-associated proteins (MAP), a myrosinase binding protein (MBP), an AOP1 oxidoreductase, and nitrile-specifier and epithio-specifier proteins were also up-regulated more strongly in the AtGL3+ line (Fig. 4). Finally, a large number of stress response genes were up-regulated in both types of transgenic cotyledons, including those responsive to biotic, drought and salt stress, those encoding pathogenesis-related (PR) proteins, PR-related protease inhibitors, heat shock proteins, and redox proteins known to assist with protection against reactive oxygen species (ROS), as well as those involved in hormonal control of growth and development (Additional file 2: Table S2; Additional file 6: Figure S4 and Additional file 1: Table S1). Genes encoding a wide variety of metabolism and transcription factors, or enzymes involved in protein modification and degradation, were strongly up-regulated in both AtGL3+ and K-5-8 cotyledons (Additional file 6: Figure S4 and Additional file 7: Figure S5).
    Table 1

    Cotyledon genes (common to two hairy lines) with potential to impact host plant responses to flea beetles and diamondback moth

       

    Fold change relative to Westar

    Gene ID

    Sub-categories

    Description

    AtGL3 line

    K-5-8 linea

    Metal Handling

     bo7g039050

    general

    SBP3 (selenium-binding protein 3)

    1.07E+301

    1.07E+301

     bo9g123330

    binding, chelation, storage

    NAS2 (NICOTIANAMINE SYNTHASE 2)

    3.96E+08

    3.92E+09

     bo4g040020

    "

    ATSERAT2;1, SAT5, SAT1 (O-SERINE ACETYLTRANSFERASE 2;1)

    8.99

    352.70

     bra002851

    "

    NAS2 | NAS2 (NICOTIANAMINE SYNTHASE 2)

    47.80

    77.58

     bra015594

    "

    MT1C; copper ion binding

    4.42

    10.20

     bra009595

    "

    MT2A, ATMT-K, ATMT-1 (METALLOTHIONEIN 2A)

    3.23

    7.38

     bo5g008330

    "

    MT1C; copper ion binding

    4.41

    7.17

    Wax

     bo8g102710

    synthesis

    KCS4 (3-KETOACYL-COA SYNTHASE 4)

    3.33

    5.83

     bra033983

    "

    YBR159, KCR1 | YBR159; ketoreductase/oxidoreductase

    4.59

    2.81

     bra004034

    "

    CUT1, POP1, CER6, G2, KCS6 (3-KETOACYL-COA SYNTHASE 6)

    29.07

    2.58

     bo7g019710

    "

    KCS9 (3-KETOACYL-COA SYNTHASE 9)

    1.67

    2.42

     bra024749

    "

    CUT1, POP1, CER6, G2, KCS6 (3-KETOACYL-COA SYNTHASE 6)

    2.83

    1.51

    Cell Wall

     bra011899

    modification

    ATEXLA2, EXPL2, ATHEXP BETA 2.2 (EXPANSIN-LIKE A2)

    9.17

    12.64

     bra024848

    "

    EXGT-A3, XTH27, hydrolase, xyloglucan:xyloglucosyl transferase

    3.29

    4.03

     bo3g018260

    precursor synthesis

    UXS3 (UDP-GLUCURONIC ACID DECARBOXYLASE 3)

    4.72

    18.03

     bra006722

    "

    "

    3.87

    9.62

     bra021798

    "

    ATCSLB03, ATCSLB3 cellulose synthase/ glycosyl transferase

    6.43

    5.18

     bo1g037390

    "

    ATCSLG2 cellulose synthase/ glycosyl transferase

    3.73

    3.17

     bra016440

    cell wall proteins

    unknown protein

    8.59

    101.30

     bo4g194100

    "

    proline-rich extensin-like family protein

    5.50E+04

    11.24

     bra037731

    "

    "

    33328.67

    7.21

     bra026268

    "

    hydroxyproline-rich glycoprotein family protein

    20.79

    5.06

     bo1g051380

    "

    "

    9.65

    4.97

     bra003200

    degradation.mannan-xylose-arabinose-fucose

    PMR6 (powdery mildew resistant 6); pectate lyase

    5.41E+12

    190.90

     bra024089

    "

    MERI5B, MERI-5, SEN4 (meristem-5); hydrolyzing glycosyl bonds

    3.13

    66.59

     bo4g108130

    "

    pectinesterase family protein

    6.87

    43.88

     bra009234

    "

    ATPGIP1 (POLYGALACTURONASE INHIBITING PROTEIN 1); binding

    44.74

    37.33

     bo2g127040

    "

    pectinacetylesterase, putative

    3.30

    19.53

     bra010038

    "

    ATBXL1 (BETA-XYLOSIDASE 1); hydrolyzing O-glycosyl compounds

    1.70

    9.03

     bo2g013480

    degradation.pectate lyases/polygalacturonases

    glycoside hydrolase family 28 protein

    5.78

    7.76

     bo6g058470

    "

    PMR6 (powdery mildew resistant 6); pectate lyase

    1.96E+05

    4.92

    Lignin

     bo3g024650

    biosynthesis

    ATC4H, C4H, CYP73A5 (CINNAMATE-4-HYDROXYLASE)

    2717.78

    39.35

     bra036480

    "

    HCT (SHIKIMATE/QUINATE HYDROXYCINNAMOYL COA-TRANSFERASE)

    2.72

    2.68

    Flavonoids

     bra018364

    dihydroflavonols

    UGT71D1 (UDP-GLUCOSYL TRANSFERASE 71D1)

    7.75

    12.33

     bra037386

    "

    "

    7.15

    9.32

     bo9g003740

    "

    "

    65.69

    6.15

    Phenylpropanoids

     bra028893

    phenylpropanoids

    transferase family protein

    14.84

    17.57

     bra029364

    "

    NIC2 (NICOTINAMIDASE 2); catalytic nicotinamidase

    5.46

    8.84

     bo3g167180

    "

    O-methyltransferase family 2 protein

    10.07

    8.41

     bo7g064020

    "

    transferase family protein

    12.06

    4.93

    Alkaloids

     bra003263

    N misc. alkaloid-like

    strictosidine synthase family protein

    416.91

    3.94

     bo6g054230

    "

    "

    17.29

    2.92

    Cyanogenic glucosides

     bra014956

    cyanase

    CYN (CYANASE); DNA binding / cyanate hydratase/ hydro-lyase

    7.06

    6.21

    aCotyledon values are arranged from highest expression to lowest expression within each functional category using the K-5-8 line.

    Fig. 4
    Fig. 4

    Changes in expression of genes related to glucosinolate biosynthesis and degradation in 7-day-old cotyledons grown under a 16/8 h diurnal light cycle. Figure shows GS UPLC peaks and gene expression changes (up-regulation) for the AtGL3+ line (A) and the K-5-8 line (K), each relative to B. napus cv. Westar levels. Br, Brassica rapa A-genome homeologue. Bo, Brassica oleracea C-genome homeologue. Gene IDs and expression levels can be viewed in Table 2

    In general, more genes in each category were upregulated in the AtGL3+ line than in the K-5-8 line relative to Westar, for example, in hormones (BIN #17), stress (BIN #20), miscellaneous (BIN #26), protein (BIN #29), signalling (BIN #30), development (BIN #33), and transport (BIN #34) (Fig. 3; Additional file 1: Table S1). Up-regulated genes that were common to both hairy leaf lines often exhibited a greater level of expression in the AtGL3+ line than in the K-5-8 line (Tables 1, 2, and 3). This included genes involved in cell wall synthesis (cell wall proteins, pectate lyases, and lignin synthesis), and stress-response genes encoding chlorophyllase (COR1), PR protease inhibitors, and a SENESCENCE ASSOCIATED GENE 20, while the same genes were only moderately up-regulated within K-5-8 cotyledons. Up-regulated genes that were unique to AtGL3+ and K-5-8 cotyledons were linked to sulphur metabolism, metal-handling, secondary metabolism (anthocyanins, dihydroflavonols, carotenoids, non-mevalonate isoprenoids, alkaloids, phenylpropanoids, lignin, GS synthesis and GS degradation), hormones (auxin, brassinosteroids, cytokinin, ethylene, gibberellic acid, and jasmonate), and stress responses (general, biotic, PR proteins, protease inhibitors, defensins, heat shock, touch or wounding, drought and salt, and non-specified abiotic stress) (Additional file 1: Table S1). In addition, more redox-related genes potentially impacting the levels of reactive oxygen species (ROS) were up-regulated in AtGL3+ cotyledons than in K-5-8 cotyledons relative to Westar, although up-regulated redox-related genes unique to K-5-8 cotyledons had very strong transcriptional responses. In addition, genes involved in aspects of development, lipid metabolism, amino acid metabolism, photosynthesis, heat shock, and cold-response were more frequently and more strongly up-regulated in K-5-8 cotyledons (Additional file 1: Table S1).
    Table 2

    Strongly up-regulated stress-responsive genes changes common to transgenic hairy B. napus cotyledons

      

    Expression relative to Westar

    ID

    Description

    AtGL3 line

    K-5-8 line

    Biotic Stress

     bo5g027670

    CORI1, ATHCOR1, ATCLH1 (CORONATINE-INDUCED 1); chlorophyllase

    5.30E+07

    21.57

     bra008224

    pathogenesis-related thaumatin family protein

    2.83

    14.46

     bo2g069600

    MLP28 (MLP-LIKE PROTEIN 28)

    174.64

    13.00

     bra007947

    "

    119.90

    11.52

     bra016785

    RPS5 (RESISTANT TO P. SYRINGAE 5); nucleotide binding

    7.66

    8.11

     bra008667

    glycosyl hydrolase family 81 protein

    6.06

    7.46

     bo5g149860

    CHAT (acetyl CoA:(Z)-3-hexen-1-ol acetyltransferase)

    46.71

    5.28

     bra011734

    ATRCCR | ACD2 (ACCELERATED CELL DEATH 2); red chlorophyll catabolite reductase

    4.47

    4.47

     bo9g163710

    glycosyl hydrolase family 81 protein

    7.85

    3.96

     bo07027s010

    HRT, RCY1, RPP8 (RECOGNITION OF PERONOSPORA PARASITICA 8); binds nucleotides

    3.46

    3.84

     bra009184

    NHL3

    3.83

    3.57

     bo7g087120

    avirulence induced gene (AIG) protein, putative

    3.49

    3.26

    Biotic Stress Signalling

     bra031065

    TIFY10A | JAZ1 (JASMONATE-ZIM-DOMAIN PROTEIN 1); protein binding

    5.35E+14

    2.49E+08

     bo8g102890

    "

    369.12

    66.76

     bo5g027170

    "

    4.17

    12.21

    Biotic Stress Regulation of Transcription

     bra027377

    RSH2 (RELA-SPOT HOMOLOG 2); GTP diphosphokinase

    4.97

    5.09

     bo5g131760

    "

    3.24

    3.25

    PR proteins General

     bo8g091760

    disease resistance protein (TIR-NBS-LRR class), putative

    5.81

    6.67

     bra005378

    disease resistance family protein

    3.09

    6.40

     bo6g007620

    molecular_function unknown; LOCATED IN: endomembrane system

    3.94

    5.00

     bo1g048080

    disease resistance protein (NBS-LRR class), putative

    1.23

    3.75

    PR proteins Protease Inhibitors

     bo6g010170

    trypsin and protease inhibitor family protein / Kunitz family protein

    165.56

    3.67

     bo6g010250

    "

    11.99

    3.54

     bo6g010100

    "

    113.13

    3.33

     bra015999

    "

    78.81

    3.25

     bra037702

    trypsin inhibitor, putative

    18.46

    1.40

     bra016073

    trypsin and protease inhibitor family protein / Kunitz family protein

    49.47

    1.07

    Wounding

     bra034157

    WI12, SAG20 (SENESCENCE ASSOCIATED GENE 20)

    1.07E+301

    1.07E+301

     bra029887

    "

    2.93E+05

    4.53E+04

     bra010381

    wound-responsive protein-related

    6697.92

    76.02

     bo7g111980

    "

    5.51

    10.78

     bo8g054160

    "

    6.39

    9.56

    Abiotic Stress General

     bo01463s030

    benzodiazepine receptor-related

    1.90E+04

    4.85E+05

     bra021442

    "

    798.78

    32627.32

     bra007841

    unknown protein

    388.29

    266.92

     bo8g099690

    "

    96.36

    115.05

     bo4g154720

    "

    2.96

    6.57

    Heat responsive

     bo1g138440

    DNAJ heat shock N-terminal domain-containing protein

    1.25E+15

    7.51E+23

     bo9g026330

    DNAJ heat shock protein, putative

    6642.31

    1.68E+11

     bra037247

    "

    7768.93

    1.62E+10

     bo7g117750

    DNAJ heat shock N-terminal domain-containing protein (J11)

    1.46E+08

    7.57E+09

     bra039384

    DNAJ heat shock N-terminal domain-containing protein

    1238.27

    1.74E+06

     bra011656

    DNAJ heat shock N-terminal domain-containing protein (J11)

    216.64

    2.14E+04

     bo1g005990

    "

    415.52

    1.04E+04

     bra034691

    DNAJ heat shock N-terminal domain-containing protein

    12.28

    139.25

     bra017744

    DNAJ heat shock N-terminal domain-containing protein (J11)

    18.63

    77.88

     bra020505

    DNAJ heat shock protein, putative

    8.76E+04

    24.75

     bo5g132640

    DNAJ heat shock N-terminal domain-containing protein

    5.58

    23.31

     bra011735

    ATHSF4 | HSF4 (HEAT SHOCK FACTOR 4); DNA binding /transcription repressor

    7.14E+04

    18.15

     bra018216

    17.6 kDa class I small heat shock protein (HSP17.6C-CI) (AA 1-156)

    47.14

    12.92

     bo2g158600

    DNAJ heat shock protein, putative

    788.57

    7.39

     bo8g066630

    Hsp70b (heat shock protein 70B); ATP binding

    6.69

    7.32

     bo4g169420

    17.6 kDa class I small heat shock protein (HSP17.6B-CI)

    13.24

    7.18

     bo8g097710

    DNAJ heat shock N-terminal domain-containing protein

    3.07

    5.16

     bra016644

    Hsp70b (heat shock protein 70B); ATP binding

    4.04

    4.31

     bo1g134560

    DNAJ heat shock N-terminal domain-containing protein

    8.46

    3.95

     bo2g029130

    heat shock protein-related

    4.61

    3.81

     bra004457

    DNAJ heat shock N-terminal domain-containing protein

    2.72

    3.75

     bra003592

    J8; heat shock protein binding / unfolded protein binding

    4.72

    3.70

     bra020419

    heat shock protein-related

    4.29

    3.60

     bo9g176930

    DNAJ heat shock N-terminal domain-containing protein

    4.42

    3.46

     bo3g091720

    HSC70-1 (HEAT SHOCK COGNATE PROTEIN 70-1); ATP binding

    4.25

    3.29

     bo8g102330

    DNAJ heat shock family protein

    4.19

    3.10

    Cold responsive

     bra017742

    CSDP1 (cold shock domain protein 1); RNA/single/double-stranded DNA binding

    1160.12

    131.40

     bo7g117730

    "

    1.31E+14

    8.31

     bra013087

    WCOR413-LIKE, FL3-5A3 | COR413-PM1

    14.61

    2.84

    Drought and salt responsive

     bo3g052160

    early-responsive to dehydration protein-related / ERD protein-related

    177.05

    1.43E+11

     bra039623

    ATCOAD (4-phosphopantetheine adenylyltransferase)

    4.62

    5.42

     bo6g004950

    QUA2 | TSD2 (TUMOROUS SHOOT DEVELOPMENT 2); methyltransferase

    3.50

    3.64

     bo4g154160

    hydrophobic protein, putative / low temperature-salt responsive

    3.15

    3.24

    Abiotic Stress unspecified

     bo8g052730

    PHOS34 | universal stress protein (USP) family protein

    9.63

    35.62

     bo2g069600

    MLP28 (MLP-LIKE PROTEIN 28)

    174.64

    13.00

     bo7g111350

    PHOS34 | universal stress protein (USP) family protein

    6.84

    12.10

     bra007947

    MLP28 (MLP-LIKE PROTEIN 28)

    119.90

    11.52

     bra008745

    universal stress protein (USP) family protein

    4.45

    4.83

     bo5g002660

    ozone-responsive stress-related protein, putative

    4.21

    4.15

     bra022721

    PHOS34 | universal stress protein (USP) family protein

    4.15

    3.76

     bo9g165720

    "

    3.08

    3.60

     bo3g039740

    "

    3.31

    3.05

     bo6g020120

    "

    3.27

    2.94

    aCotyledon values are arranged from highest expression to lowest expression within each functional category using the K-5-8 line.

    Table 3

    Trichome-related glabrous cotyledon genes in the hairy leaf AtGL3+ and ultra-hairy leaf K-5-8 lines relative to glabrous leaf B. napus cv. Westar

    ID

    Bin Name

    Description

    Expression relative to Westar

    AtGL3 line

    K-5-8 line

    ESTs common to both transgenic lines

      

     bo5g002440

    Positive initiation

    AN (ANGUSTIFOLIA); protein binding

    5.39

    1.63

     bra024875

    "

    RGA1 (REPRESSOR OF GA1-3 1); protein binding / transcription factor

    3.29

    1.59

     bo9g070200

    "

    RGA1 (REPRESSOR OF GA1-3 1); protein binding / transcription factor

    1.49

    0.73

     bra007766

    "

    FDH, KCS10 (3-KETOACYL-COA SYNTHASE 10); acyltransferase

    1.51

    0.42

     bo1g116200

    Positive branching

    DER1, LSR2, ENL2 | ACT2 (ACTIN 2); structural constituent of cytoskeleton

    1.53

    1.46

     bra020572

    "

    TUA6; structural constituent of cytoskeleton

    1.54

    0.73

     bra039648

    "

    TUA6; structural constituent of cytoskeleton

    1.48

    0.55

     bo1g054880

    "

    SPK1 (SPIKE1); GTP binding / GTPase binding / guanyl-nucleotide exchange factor

    1.51

    0.47

     bra009451

    Multicellular trichomes

    SIM (SIAMESE); cyclin-dependent protein kinase inhibitor

    1978.24

    11.69

     bo9g178800

    "

    SIM (SIAMESE); cyclin-dependent protein kinase inhibitor

    935.76

    7.64

     bra027928

    Less developed

    unknown protein

    1.45

    0.42

     bo1g046580

    "

    MRH5, GPDL2 | SHV3 (SHAVEN 3); glycerophosphodiester phosphodiesterase/ kinase

    1.14

    0.37

     bra026409

    "

    MRH5, GPDL2 | SHV3 (SHAVEN 3); glycerophosphodiester phosphodiesterase/ kinase

    1.13

    0.35

     bo9g054590

    "

    unknown protein

    1.44

    0.24

    ESTs unique to the AtGL3+ line

      

     bra033258

    Positive initiation

    AN (ANGUSTIFOLIA); protein binding

    14.72

    NA

     bra017443

    "

    RGA1 (REPRESSOR OF GA1-3 1); protein binding / transcription factor

    1.47

    NA

     bra029388

    Distorted

    EMB3009 (embryo defective 3009); transferase/ transferase, transferring acyl groups

    1.07E+301

    NA

     bo7g097360

    "

    EMB3009 (embryo defective 3009); transferase/ transferase, transferring acyl groups

    1.05E+59

    NA

     bra004054

    Positive branching

    TPS6 | ATTPS6; alpha,alpha-trehalose-phosphate synthase (UDP-forming)/ transferase

    3.29

    NA

     bo3g071760

    "

    DER1, LSR2, ENL2 | ACT2 (ACTIN 2); structural constituent of cytoskeleton

    1.68

    NA

     bra022356

    "

    DER1, LSR2, ENL2 | ACT2 (ACTIN 2); structural constituent of cytoskeleton

    1.64

    NA

     bo5g117040

    "

    DER1, LSR2, ENL2 | ACT2 (ACTIN 2); structural constituent of cytoskeleton

    1.55

    NA

     bra037560

    "

    DER1, LSR2, ENL2 | ACT2 (ACTIN 2); structural constituent of cytoskeleton

    1.48

    NA

     bra008705

    Negative branching

    ATMIXTA | ATMYB16 (MYB DOMAIN PROTEIN 16); DNA binding / transcription factor

    5.94

    NA

     bo5g021100

    Endoreduplication

    CYCA2;3 (CYCLIN A2;3); cyclin-dependent protein kinase regulator

    1.39

    NA

     bo5g038710

    Trichome size

    ATSAC1 (suppressor of actin 1); phosphatidylinositol-4,5-bisphosphate 5-phosphatase

    7.73

    NA

     bo6g027740

    Less developed

    MRH5, GPDL2 | SHV3 (SHAVEN 3); glycerophosphodiester phosphodiesterase/ kinase

    38.59

    NA

    ESTs unique to the K-5-8 line

      

     bo8g115250

    Positive initiation

    HDG2 (HOMEODOMAIN GLABROUS 2); DNA binding / transcription factor

    NA

    0.43

     bo8g067530

    "

    HDG12 (HOMEODOMAIN GLABROUS 12); transcription factor

    NA

    0.27

     bra015401

    "

    HDG2 (HOMEODOMAIN GLABROUS 2); DNA binding / transcription factor

    NA

    0.13

     bo6g029470

    Positive branching

    TPS6 | ATTPS6; alpha,alpha-trehalose-phosphate synthase (UDP-forming)/ transferase

    NA

    1.65

     bo6g107020

    "

    TUA6; structural constituent of cytoskeleton

    NA

    0.63

     bo9g022270

    "

    PKCBP, KCBP | ZWI (ZWICHEL); calmodulin binding / microtubule motor

    NA

    0.57

     bra018825

    "

    TUA6; structural constituent of cytoskeleton

    NA

    0.53

     bra016164

    Negative initiation

    MYBL2 (ARABIDOPSIS MYB-LIKE 2); DNA binding / transcription factor

    NA

    1.56

     bo3g149420

    "

    UPL3 | KAK (KAKTUS); ubiquitin-protein ligase | chr4:18041031-18049292 REVERSE'

    NA

    0.68

     bo2g012520

    Negative branching

    ATMIXTA | ATMYB16 (MYB DOMAIN PROTEIN 16); DNA binding / transcription factor

    NA

    0.36

     bo3g010230

    "

    ATMIXTA | ATMYB16 (MYB DOMAIN PROTEIN 16); DNA binding / transcription factor

    NA

    0.29

     bo9g164230

    "

    ATMIXTA | ATMYB16 (MYB DOMAIN PROTEIN 16); DNA binding / transcription factor

    NA

    0.18

     bra024337

    Trichome size

    HYS1 | CPR5 (CONSTITUTIVE EXPRESSION OF PR GENES 5)

    NA

    1.94

     bo2g029050

    "

    FLP1, YRE, CER3, WAX2 | CER3 (ECERIFERUM 3); binding / catalytic/ iron ion binding /oxidoreductase

    NA

    0.56

     bra020412

    "

    FLP1, YRE, CER3, WAX2 | CER3 (ECERIFERUM 3); binding / catalytic/ iron ion binding /oxidoreductase

    NA

    0.41

     bo9g133540

    "

    FLP1, YRE, CER3, WAX2 | CER3 (ECERIFERUM 3); binding / catalytic/ iron ion binding /oxidoreductase

    NA

    0.20

    aCotyledon values are arranged from highest expression to lowest expression within each functional category using the K-5-8 line

    Expression of trichome genes in glabrous cotyledons of hairy-leaf and ultrahairy-leaf lines

    The transgenic cotyledons remained glabrous despite the expression of AtGL3 or in conjunction with the manipulation of BnTTG1 expression (2). Hence, transcript levels for all known trichome genes were compared in both types of transgenic cotyledons to those of cv. Westar. According to RNA sequencing, no trichome genes were down-regulated in the glabrous transgenic cotyledons relative to cv. Westar cotyledons, and up-regulation of most transcripts involved trichome genes that should impact trichome structure (Table 3). Overall, the greatest change in transcript level relative to cv. Westar cotyledons included two SIAMESE (SIM) genes (bra009451 and bo9g178890) specifying multi-cellular trichomes that were also more highly expressed in the AtGL3+ line than in the K-5-8 line (Table 3). As well, two embryo defective 3009 genes (EMB3009) encoding acyl transferases involved in trichome shape were highly up-regulated in the AtGL3+ line, while SAC1 (affecting trichome size), SHAVEN3, and the MYB16 MIXTA genes were also up-regulated in this line (Table 3). A paralogue of the BrRGA1 gene was uniquely up-regulated in the AtGL3+ cotyledons and two RGA1 genes were commonly up-regulated in both transgenic lines.

    Well-known trichome genes coding for the MBW tri-protein initiation complex were not differentially expressed in 7-day-old cotyledons. Only one well-known trichome regulatory gene, ANGUSTIFOLIA (AN, affecting trichome initiation), had a significantly different transcript profile in AtGL3+ or K-5-8 cotyledons compared with cv. Westar (Table 3). qRT-PCR with 10-day-old cotyledons showed that five well-known MBW trichome genes were weakly expressed. Although the BnGL2 trichome-initiation gene was expressed at a higher level in AtGL3+ cotyledons compared with K-5-8 and Westar cotyledons, BnGL3 and BnTRY were expressed at lower levels, while changes for BnTTG1 (which was manipulated) occurred only in the K-5-8 line (Fig. 5) as expected.
    Fig. 5
    Fig. 5

    qRT-PCR analysis of expression levels of five trichome regulatory genes encoding elements of the MBW tri-protein complex in glabrous cotyledons of 10-day-old hairy leaf (AtGL3+) and ultra-hairy leaf (K-5-8) B. napus lines relative to cv. Westar (W). Level of BnACT2 and AtGL3 are provided as controls. Plants were grown under a 16/8 h diurnal light cycle

    Discussion

    This study examined the composition of glabrous cotyledons from a hairy leaf line derived from the introduction of the AtGL3 gene (line AtGL3+) into B. napus [1] and an ultra-hairy leaf line derived by repression of BnTTG1 within the AtGL3+ B. napus background (line K-5-8) [2]. Expression of AtGL3 induces the formation of trichomes on young leaves of these lines [1, 2] which reduces the insect’s ability to physically interact with the host plant and reduces FB feeding damage; however, the glabrous cotyledons of these lines also deter FB feeding [21]. A comparison of the physiological properties and composition of the cotyledons (Table 4) provides a few clues as to whether a common mechanism(s) might be responsible for the antixenosis. In summary, the cotyledons of the AtGL3+ have an abnormal appearance, which is in keeping with vegetative tissue in this line, but develop normally in the K-5-8 line due to repression of BnTTG1 expression (2). The cotyledons in the wild-type line are horizontally-oriented; however, the cotyledons in both transgenic lines are vertically-oriented and it was suggested that this orientation may be less suitable for FBs to initiate feeding behaviour [19]. Anthocyanin content in the AtGL3+ line was significantly increased, but reduced in the K-5-8 line compared to wild-type. Both lines had similar cell wall carbohydrate contents; however, cell wall lignin was reduced in the K-5-8 line. While total GS content was reduced in the AtGL3+ line, the level of gluconapin was increased and 4-methoxy-glucobrassicin was decreased in both transgenic lines.
    Table 4

    Summary of physiological properties and composition of B. napus cotyledons expressing AtGL3

    A. Insect behaviour and cotyledon biochemistry relative to Westar

    Genotype

    FB Feeding

    Cotyledon Morphology

    Cotyledon Orientation

    Anthocyanin Accumulation

    Cell Wall Carbohydrates

    Cell Wall Lignin

    AtGL3+a

    Decreased

    Abnormal

    Vertical

    Increased

    No Difference

    No Difference

    K-5-8b

    Decreased

    Normal

    Vertical

    Decreased

    No Difference

    Decreased

    B. Glucosinolate composition of transgenic lines relative to Westar

    Genotype

    Total GS

    Progoitrinc

    Glucoraphaninc

    Gluconapinc

    4-hydroxy-glucobrassicinc

    Glucobrassicinc

    4-methoxy-glucobrassicinc

    AtGL3+

    Decreased

    No Diference

    No Difference

    Increased

    No Difference

    No Difference

    Decreased

    K-5-8

    No Difference

    No Difference

    No Difference

    Increased

    Decreased

    No Difference

    Decreased

    C. Glucosinolate composition of the transgenic lines relative to one another.

    Genotype

    Total GS

    Progoitrinc

    Glucoraphaninc

    Gluconapinc

    4-hydroxy-glucobrassicinc

    Glucobrassicinc

    4-methoxy-glucobrassicinc

    AtGL3+

    Decreased

    Decreased

    No Difference

    Increased

    No Difference

    Increased

    Decreased

    K-5-8

    -

    -

    -

    -

    -

    -

    -

    aAtGL3+, B. napus cv. Westar expressing AtGL3 under direction of CaMV 35S promoter

    bK-5-8, B. napus cv. Westar expressing AtGL3 under direction of CaMV 35S promoter as well as a TTG1 RNAi construct

    cProgoitrin, 2-hydroxy-3-butenyl-GS; Glucoraphanin, 4-methylsulfinylbutyl-GS; Gluconapin, 3-butenyl-GS; 4-hydroxyglucobrassicin, 4-hydroxy-3-indolylmethyl-GS; Glucobrassicin, 3-indolylmethyl-GS; 4-methoxyglucobrassicin, 4-methoxy-3-indolylmethyl-GS

    Role of secondary metabolites in cotyledon resistance

    Feeding is highly influenced by GS content for many crucifer specialists; however, there is little consensus as to which GSs stimulate and which deter feeding. Related to the GSs found in cotyledons in the current study, decreased quantities of glucoraphanin and increased levels of progoitrin correlated with increased FB feeding in Sinapis alba [31]. In a broad study of various Brassica species, progoitrin levels were also correlated with stimulation of FB feeding [32]; however, this does not explain the reduced FB feeding on the transgenic cotyledons in this study as progoitrin and glucoraphanin levels were similar to the wild-type line. The same study and studies with the stem FB (Psylliodes chrysocephala) [33, 34] reported that gluconapin is also a feeding stimulant; however, this GS was elevated in cotyledons of both transgenic lines in the current study. Interestingly, glucobrassicin was the most stimulatory of the GS tested in the study with P. chrysocephala [33], but levels of this GS were the same in the transgenic and wild-type lines. To our knowledge, no studies have implicated 4-methoxy-glucobrassicin in affecting FB feeding activity. Future experiments appear to warrant examining individual GS to determine, for example, if applying gluconapin further increases resistance or if adding 4-methoxy-glucobrassicin restores wild-type predation by FB in the transgenic cotyledons.

    It has been difficult to extrapolate the response to individual GSs in laboratory studies to damage occurring in the field [35, 36] as this is influenced by the presence of other GSs and other plant secondary metabolites, as well as the myriad GS-derived volatiles that result from the release of myrosinase and associated myrosinase-specifier proteins upon tissue disruption [37]. Communication and recruitment of con-specific insects must also be considered when assessing host susceptibility. In the field, allyl isothiocyanate derived from the GS sinagrin is a strong FB attractant [38]. This volatile also enhances the response of FBs to aggregation pheromone [39], which functions only in the presence of specific host plant volatiles [40]. A role for flavonoids in host recognition and acceptance by many adult insects has also been reported [41, 42]. Moreover, anthocyanins are known to impact insect feeding [28]. Since anthocyanin production was altered in the AtGL3+ cotyledons, these should also be tested for their impact on FB feeding.

    Impact of AtGL3 expression on the cotyledon transcriptome

    AtGL3+ cotyledon size and colour (small with a dark red abaxial side) differ from that of K-5-8 or Westar, and AtGL3+ plants are smaller and grow less vigorously [2, 18]. Not surprisingly, the alteration of the cotyledon transcriptome in the AtGL3+ plants was more extreme (both in transcript diversity and expression intensity) compared with K-5-8 cotyledons (both lines relative to the Westar cotyledon transcriptome). In contrast, changes to the K-5-8 cotyledon transcriptome were less extreme (relative to AtGL3+ cotyledons), which may reflect the healthier cotyledons that were similar in size to Westar cotyledons and exhibited more vigorous growth [2, 18]. The cell wall of K-5-8 cotyledons was also less lignified suggesting that K-5-8 cotyledons could be in a more vigorous growth phase than the AtGL3+ cotyledons. The red anthocyanin present in AtGL3+ cotyledons growing under continuous light appeared to be part of a stronger and more varied stress response as indicated by the transcriptome data. In addition to the up-regulation of genes directly involved in anthocyanin production, genes involved in hormone synthesis or signalling (jasmonate, auxins, gibberellins, brassinosteroids, and ethylene) were also more strongly up-regulated in the AtGL3+ cotyledons than in those from K-5-8, while genes involved in development, lipid synthesis, amino acid metabolism, and photosynthesis were more strongly up-regulated in K-5-8. These data suggest that AtGL3+ cotyledons are compromised and less able to fuel growth, while K-5-8 cotyledons are primed for better growth characteristics and possibly more able to tolerate FB damage.

    Expression of genes involved in trichome formation

    Cotyledons of lines expressing AtGL3 remained glabrous; however, RNA sequencing showed that both types of transgenic cotyledons had elevated levels of transcripts for many genes involved in trichome synthesis. SIM, which encodes a cyclin-dependent protein kinase inhibitor specifying multi-cellular trichomes, was up-regulated in these cotyledons. As well, the AtGL3+ line had an elevated level of transcripts for EMB3009 which specifies trichome shape, AUGUSTIFOLIA which is involved in trichome initiation, MIXTA which specifies negative trichome branching, SAC1 which specifies trichome size, and SHAVEN3 which specifies trichome number. These differences may underlie the smaller trichomes present on the AtGL3+ line [2]. Moreover, a paralogue of the BrRGA1 gene was uniquely up-regulated in the AtGL3+ cotyledons, and other RGA1 genes were up-regulated in cotyledons of both transgenic lines. RGA1 (repressor of gibberellic acid 1) is a GA-insensitive DELLA repressor protein that negatively impacts the activity of trichome transcriptional activators [43]. Hence, induction of this repressor could be one of the reasons why the transgenic cotyledons remained glabrous even in the presence of a high AtGL3 transcript levels. Induction of RGA1 also suggests that it could be a cotyledon-specific gene, since these repressors do not appear in the leaf transcriptome of either transgenic line, or in leaves of cv. Westar [2]. The data implies that the expression of the AtGL3 gene in the two transgenic lines may have a measure of positive control by stimulating over-expression of this trichome repressor, ensuring that cotyledons will be maintained as sources of nutrients for young seedlings rather than becoming trichome-bearing leaf-like organs. This is supported by the fact that five trichome regulatory genes related to the MBW initiation complex were only weakly expressed in cotyledons.

    Shared components within the MBW regulatory complex link trichome formation and anthocyanin biosynthesis. In trichome formation, GL3 or ENHANCER OF GL3 forms a complex with TTG1 and the MYB protein GL1. This complex activates the expression of genes encoding transcription factors that in turn induce trichome formation. Anthocyanin biosynthesis is regulated by a complex consisting of TTG1 and GL3, but with different MYB factors [17]. In the current study, the chemical profile and changes in the expression of genes in cotyledons expressing AtGL3 alone (AtGL3+) or accompanied by reduced BnTTG1 expression (K-5-8) suggest that the production of other secondary metabolites, such as GSs, may be subject to a similar type of regulation. At least 6 other MYB proteins (CPC, TRY, ETC1, ETC2, ETC3 and TCL1) can replace GL1 and alter the specificity of the regulatory complex [15]. In the future, it would interesting to establish if some type of MBW complex also regulates GS production and what MYB protein(s) might be involved.

    Conclusion

    An unexpected outcome of manipulating the expression of AtGL3 and BnTTG1 genes was cotyledon resistance to FB, in spite of the fact that this tissue remained glabrous [18, 19]. As such, these glabrous cotyledon transcriptomes represent roadmaps that can lead to the identification of insect-resistance compounds and properties. The altered host interactions could reflect the large up-regulation of genes affecting the synthesis of GSs, other secondary metabolites and tissue structural components, including cell wall carbohydrates, lignin, wax, metal handling systems, flavonoids, phenolics, and indole alkaloids. Quercetin (flavonoid) and chlorogenic acid (phenolic) derivatives accumulate strongly in cotyledons of the crucifer Camelina sativa, protecting them from FB feeding [41], while B. napus cotyledons normally accumulate only traces of these compounds (Gruber unpublished). Seeds of the crucifer Lunnaria annua L. (L. bioennis) accumulate toxic lunarine/lunaridin alkaloids [44], and alkaloids have been shown to deter insects. Plants that accumulate certain metals in the Brassicaceae can also be more resistant to specific insects [4551]. Finally, the induction of large numbers stress-responsive genes (specifying wounding, abiotic, and biotic stress responses), such as those encoding heat shock proteins, PR proteins, protease inhibitors, glucosinolate synthesis/breakdown factors, abiotic stress factors, redox proteins, transcription factors, and proteins required for auxin metabolism suggest that these cotyledons may now be primed for resistance to other forms of biotic and abiotic stress.

    Declarations

    Acknowledgements

    The authors acknowledge the assistance of Wesley Soroka in seed cleaning and L. Grenkow with statistical training. I. Parkin is gratefully acknowledged for sharing data from the B. oleracea genome prior to its publication.

    Funding

    Funding was available from the Agriculture Development Fund of the Government of Saskatchewan and SaskCanola. N. Nagubushana and A. Taheri received Visiting Fellowships from the Government of Canadian.

    Availability of data and materials

    Raw sequencing data files were deposited within NCBI as accession #SRP065063. Seeds of the three lines used in this study are available from D. Hegedus.

    Authors’ contributions

    MG and PB-S supervised the experiments. AU, AT, NN, MYG, and RZ conducted experiments and analysed data. PB-S, MG, and DD wrote and edited the manuscript. AS provided resources and instrumentation for RNA sequencing. AH provided resources for analysis of cotyledon structural composition. RZ provided resources and conducted glucosinolate analysis. All authors read and approved the final manuscript.

    Ethics approval and consent to participate

    Not applicable

    Competing interests

    The authors declare that they have no competing interests.

    Publisher’s Note

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

    Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

    Authors’ Affiliations

    (1)
    Saskatoon Research Centre, Agriculture and Agri-Food Canada, 107 Science Place, Saskatoon, SK, S7N0X2, Canada
    (2)
    Department of Biology, University of Saskatchewan, Saskatoon, SK, Canada
    (3)
    Present Address: DOW Agro-Sciences, Saskatoon, SK, Canada
    (4)
    Present Address: Department of Agricultural and Environmental Sciences, Tennessee State University, Nashville, TN, USA
    (5)
    Global Institute for Food Security, University of Saskatchewan, Saskatoon, SK, Canada
    (6)
    Agriculture and Agri-Food Canada, London, ON, Canada
    (7)
    Department of Biology, Western University, London, ON, Canada
    (8)
    Department of Food and Bio-Product Sciences, University of Saskatchewan, Saskatoon, SK, Canada

    References

    1. Gruber MY, Wang S, Ethier S, Holowachuk J, Bonham-Smith PC, Soroka J, Lloyd A. “Hairy canola”- Arabidopsis GL3 induces a dense covering of trichomes on Brassica napus seedlings. Plant Mol Biol. 2006;60:679–98.View ArticlePubMedGoogle Scholar
    2. Alahakoon U, Taheri A, Naghabushana N, Bonham-Smith PC, Gruber MY. Hairy canola re-visited: down-regulating TTG1 in an AtGL3-enhanced hairy leaf background improves growth, leaf trichome coverage, and metabolite gene expression diversity. BMC Plant Biol. 2016a;16:e12. https://doi.org/10.1186/s12870-015-0680-5.View ArticleGoogle Scholar
    3. Kessler A, Baldwin IT. Plant responses to insect herbivory: the emerging molecular analysis. Ann Rev Plant Biol. 2002;53:299–328.View ArticleGoogle Scholar
    4. Painter RH. Insect Resistance in Crop Plants 1951. New York: MacMillan.Google Scholar
    5. Canola Council of Canada. Flea beetles. Canola Encyclopedia. Canola Council of Canada. http://www.canolacouncil.org/canola-encyclopedia/insects [Feb. 10. 2015].
    6. Crop Protection Compendium. Phyllotreta cruciferae. CABI. http://www.cabi.org/cpc/datasheet/40780 [Feb. 10, 2015].
    7. Crop Protection Compendium. Phyllotreta striolata. CABI. http://www.cabi.org/cpc/datasheet/40784 [Feb. 2, 2015].
    8. Tansey JA, Dosdall LM, Keddie BA, Sarfraz RM. Differences in Phyllotreta cruciferae and Phyllotreta striolata (Coleoptera: Chrysomelidae) responses to neonicotinoid seed treatments. J Econ Entomol. 2008;101:159–67.View ArticlePubMedGoogle Scholar
    9. Heimbach U, Müller A. Incidence of pyrethroid-resistant oilseed rape pests in Germany. Pest Management Sci. 2012;69:209–16. https://doi.org/10.1002/ps.3351.View ArticleGoogle Scholar
    10. Woodman RL, Fernandes GW. Differential mechanical defense: herbivory, evapotranspiration, and leaf-hairs. Oikos. 1991;60:11–9.View ArticleGoogle Scholar
    11. Larkin JC, Young N, Prigge M, Marks D. The control of trichome spacing and number in Arabidopsis. Development. 1996;122:997–1005.PubMedGoogle Scholar
    12. Handley R, Ekbom B, Agren J. Variation in trichome density and resistance against a specialist insect herbivore in natural populations of Arabidopsis thaliana. Ecol Entomol. 2005;305:284–92.View ArticleGoogle Scholar
    13. Palaniswamy P, Bodnaryk RP. A wild Brassica from Sicily provides trichome-based resistance against flea beetles, Phyllotreta cruciferae (Goeze) (Coleoptera: Chrysomelidae). Can Entomol. 1994;126:1119–30.View ArticleGoogle Scholar
    14. Payne CT, Zhang F, Lloyd AM. GL3 encodes a bHLH protein that regulates trichome development in Arabidopsis through interaction with GL1 and TTG1. Genetics. 2000;156:1349–62.PubMedPubMed CentralGoogle Scholar
    15. Zhang F, Gonzalez A, Zhao M, Payne CT, Lloyd AMA. Network of redundant bHLH proteins in all TTG1-dependent pathways of Arabidopsis. Development. 2013;130:4859–69.View ArticleGoogle Scholar
    16. Pesch M, Schultheib I, Klopffleisch K, Uhrig JF, Koegl M, Clemen CS, Simon R, Weidtkamp-Peters S, Hulskamp M. TRANSPARENT TESTA GLABRA1 and GLABRA1 compete for binding to GLABRA3 in Arabidopsis. Plant Physiol. 2015;168:584–97.View ArticlePubMedPubMed CentralGoogle Scholar
    17. Gonzalez A, Zhao M, Leavitt JM, Lloyd AM. Regulation of the anthocyanin biosynthetic pathway by the TTG1/bHLH/Myb transcriptional complex in Arabidopsis seedlings. Plant J. 2008;53:814–27.View ArticlePubMedGoogle Scholar
    18. Alahakoon U, Adamson J, Grenkow L, Soroka J, Bonham-Smith P, Gruber M. Field growth traits and insect-host plant interactions of two transgenic canola (Brassicaceae) lines with elevated trichome numbers. Can Entomol. 2016b;1:1–13. https://doi.org/10.4039/tce.2016.9.Google Scholar
    19. Soroka JJ, Holowachuk JM, Gruber MY, Grenkow LF. Feeding by flea beetles (Coleoptera: Chrysomelidae; Phyllotreta spp.) is decreased on canola (Brassica napus) seedlings with increased trichome density. J Econ Entomol. 2011;104:125–36.View ArticlePubMedGoogle Scholar
    20. Soroka J, Grenkow L. Susceptibility of Brassicaceous plants to feeding by flea beetles, Phyllotreta spp. (Coleoptera:Chrysomelidae). J Econ Entomol. 2013;106:2557–67.View ArticlePubMedGoogle Scholar
    21. Updegraff DM. Semimicro determination of cellulose in biological materials. Anal Biochem. 1969;32:20–4.View ArticleGoogle Scholar
    22. Brinkmann K, Blaschke L, Polle AJ. Comparison of different methods for lignin determination as a basis for calibration of near-infrared reflectance spectroscopy and implications of lignoproteins. J Chem Ecol. 2002;28:2483–501.View ArticlePubMedGoogle Scholar
    23. Aung B, Gruber MY, Amyot L, Omari K, Bertrand A, Hannoufa A. MicroRNA156 as a promising tool for alfalfa improvement. Plant Biotechnol J. 2015;13:779–90.View ArticlePubMedGoogle Scholar
    24. Theander O, Aman P, Westerlund E, Andersson R, Pettersson D. Total dietary fiber determined as neutral sugar residues, uronic acid residues, and Klason lignin (the Uppsala method): collaborative study. AOAC Int. 1995;78:1030–44.Google Scholar
    25. Thimm O, Blasing O, Gibon Y, Nagel A, Meyer S, Kruger P, Selbig J, Muller LA, Rhee SY, Stitt M. MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004;37:914–39.View ArticlePubMedGoogle Scholar
    26. SAS Institute, SAS user's guide, version 9.2, SAS Institute. 2008; Cary. https://support.sas.com/documentation/cdl/en/statugintroduction/61750/PDF/default/statugintroduction.pdf
    27. Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and cufflinks. Nat Prot. 2012;7:562–78.View ArticleGoogle Scholar
    28. Karageorgou P, Manetas Y. The importance of being red when young: anthocyanins and the protection of young leaves of Quercus coccifer from insect herbivory and excess light. Tree Physiol. 2006;26:613–21.View ArticlePubMedGoogle Scholar
    29. Ishida M, Hara M, Fukino N, Kakizaki T, Morimitsu Y. Glucosinolate metabolism, functionality and breeding for the improvement of Brasicaceae vegetables. Breeding Sci. 2014;64:48–59.View ArticleGoogle Scholar
    30. Santiogo R, Barros-Rios J, Malvar RA. Impact of cell wall composition on maize resistance to pests and diseases. Int J Mol Sci. 2013;14:6960–80.View ArticleGoogle Scholar
    31. Soroka J, Grenkow L. Susceptibility of brassicaceous plants to feeding by flea beetles, Phyllotreta spp. (Coleoptera: Chrysomelidae). J Econ Entomol. 2013;106:2557–67.View ArticlePubMedGoogle Scholar
    32. Bohinc T, Kosir IJ, Trdan S. Glucosinolates as arsenal for defending Brassicas against cabbage flea beetle (Phyllotreta spp.) attack. Zemdirbyste-Agriculture. 2013;100:199–204.View ArticleGoogle Scholar
    33. Barlet E, Parsons D, Williams IH, Clark SJ. The influence of glucosinolates and sugars on feeding by the cabbage stem flea beetle, Psylloides chrysocephala. Entomol Exp Appl. 1994;73:7–83.Google Scholar
    34. Giamoustaris A, Mithen R. The effect of modifying the glucosinolate content of leaves of oilseed rape (Brassica napus spp. oleifera) on its interaction with specialist and generalist pests. Ann Appl Biol. 1995;126:347–63.View ArticleGoogle Scholar
    35. Bartlet E, Williams IH. Factors restricting the feeding of the cabbage stem flea beetle (Psylliodes chrysocephala). Entomol Exp Appl. 1991;60:233–8.View ArticleGoogle Scholar
    36. Barlet E, Mithen R, Clark SJ. Feeding of the cabbage stem flea beetle Psylloides chrysocephala on high and low glucosinolate cultivars of oilseed rape. Entomol Exp Appl. 1996;80:87–9.View ArticleGoogle Scholar
    37. Hopkins RJ, van Dam NM, van Loon JJA. Role of glucosinolates in insect-plant relationships and multitrophic interactions. Annu Rev Entomol. 2009;54:57–83.View ArticlePubMedGoogle Scholar
    38. Feeny P, Paauwe KL, Demong NJ. Flea beetles and mustard oils: host plant specificity of Phyllotreta cruciferae and P. striolata adults (Coleoptera: Chrysomelidae). Annal Entomol Soc Am. 1970;63:832–41.View ArticleGoogle Scholar
    39. Soroka JJ, Bartelt R, Zilkowski BW, Cosse AA. Responses of flea beetle Phyllotreta cruciferae to synthetic aggregation pheromone components and host plant volatiles in field trials. J Chem Ecol. 2005;31:1829–43.View ArticlePubMedGoogle Scholar
    40. Beran F, Mewis I, Srinivasan R, Svoboda J, Vial C, Mosimann H, Boland W, Buttner C, Ulrichs C, Hansson BS, Reinecke A. Male Phyllotreta striolata (F.) produce an aggregation pheromone: identification of male-specific compounds and interaction with host plant volatiles. J Chem Ecol. 2011;37:85–97.View ArticlePubMedGoogle Scholar
    41. Onyilagha JC, Gruber MY, Hallett RH, Holowachuk J, Soroka JJ. Constitutive flavonoids deter flea beetle insect feeding in Camelina sativa L. Biochem Syst Ecol. 2012;42:128–33.View ArticleGoogle Scholar
    42. Simmonds MSJ. Importance of flavonoids in insect-plant interactions: feeding and oviposition. Phytochemistry. 2001;56:245–52.View ArticlePubMedGoogle Scholar
    43. Gan Y. Yu H, Peng J Broun P..Genetic and molecular regulation by DELLA proteins of trichome development in Arabidopsis. Plant Physiol. 2007;145:1031–42.View ArticlePubMedPubMed CentralGoogle Scholar
    44. Henderson FG, Chen KK. The pharmacology of lunarine, the alkaloid of Lunaria biennis. J Amer Pharm Assoc. 1950;39:516–9.View ArticleGoogle Scholar
    45. Broadhurst CL, Chaney RL, Angle JS, Erbe EF, Maugel TK. Nickel localization and response to increasing Ni soil levels in leaves of the Ni hyperaccumulator Alyssum murale. Plant Soil. 2004;265:225–42.View ArticleGoogle Scholar
    46. Isaure MP, Fayard B, Sarret G, Pairis S, Bourguignon J. Localization and chemical forms of cadmium in plant samples by combining analytical electron microscopy and X-ray spectromicroscopy. Spectrochim Acta B At Spectrosc. 2006;61:1242–52.View ArticleGoogle Scholar
    47. Küpper H, Lombi E, Zhao FJ, Wieshammer G, McGrath SP. Cellular compartmentation of nickel in the hyperaccumulators Alyssum lesbiacum, Alyssum bertolonii and Thlaspi goesingense. J Exp Bot. 2001;52:2291–300.View ArticlePubMedGoogle Scholar
    48. Küpper H, Lombi E, Zhao FJ, McGrath SP. Cellular compartmentation of cadmium and zinc in relation to other elements in the hyperaccumulator Arabidopsis halleri. Planta. 2000;212:75–84.View ArticlePubMedGoogle Scholar
    49. McNear DH Jr, Peltier E, Everhart J, Chaney RL, Sutton S, Newville M, Rivers M, Sparks DL. Application of quantitative fluorescence and absorption-edge computed microtomography to image metal compartmentalization in Alyssum murale. Environ Sci Technol. 2005;39:2210–8.View ArticlePubMedGoogle Scholar
    50. Nagabushana N, Tan Y, Taheri A, Li X, Bjorndahl TC, Nowak J, Wishart DS, Hegedus D, Gruber MY. Brassica villosa, a system for studying non-glandular trichomes and genes in the brassicas. Plant Mol Biol. 2014;85:519–39.View ArticleGoogle Scholar
    51. Salt DE, Prince RC, Pickering IJ, Raskin I. Mechanisms of cadmium mobility and accumulation in Indian mustard. Plant Physiol. 1995;109:1427–33.View ArticlePubMedPubMed CentralGoogle Scholar

    Copyright

    © The Author(s). 2018

    Advertisement