Tandem Mass Tag-Based Quantitative Proteomics Reveals Implication of a Late Embryogenesis Abundant Protein (BnLEA57) in Seed Oil Accumulation in Brassica napus L.

Enhancing oil content is one of the major goals in Brassica napus breeding; however, genetic regulation of seed oil content in plants is complex and not fully elucidated. In this study, we report proteins that were differentially accumulated in immature seeds of 35 days after anthesis between two recombinant inbred lines with contrasting seed oil content, high oil content line (HOCL) and low oil content line (LOCL) using a multiplex isobaric tandem mass tags (TMT)-based quantitative proteomic approach. Over 4,600 proteins were quantified in seeds of the two lines, and 342 proteins showed differential accumulation between seeds of HOCL and LOCL. Gene Ontology enrichment analysis revealed that the differentially accumulated proteins were enriched in proteins involved in lipid biosynthesis and metabolism, photosynthesis, and nutrient reservoir activity. Western blot confirmed the increased abundance of a late embryogenesis abundant protein (BnLEA57) in HOCL seeds compared with LOCL seeds, and overexpression of either BnLEA57 gene or its homology BnLEA55 in transgenic Arabidopsis thaliana enhanced oil content in Arabidopsis seeds. Our work provides new insights into the molecular regulatory mechanism of seed oil content in B. napus.

Introduction

Rapeseed (Brassica napus L.) is an important oil-producing crop in the world because of large demand for both edible and industrial purposes (Dupont et al., 1989). Furthermore, B. napus oil has been pursued as an attractive alternative for renewable biofuels (Sharma et al., 2015). Improving seed oil content in B. napus is still one of the most important breeding objectives (Hua et al., 2016).

Brassica napus seeds accumulate oils during seed development, and the storage of seed oil is mainly in the embryo in the form of triacylglycerol (TAG) (Bai et al., 2020). The biosynthetic pathway of the fatty acid and TAG in higher plants has been well elucidated, including in Arabidopsis thaliana, Glycine max, and Brassica napus (Maeo et al., 2009; Woodfield et al., 2019; Chen et al., 2020); however, the regulatory mechanism of seed oil accumulation is only partly understood. Modification of key genes involved in fatty acid biosynthesis, lipid biosynthesis, or glycolysis is known to affect seed oil accumulation (Hua et al., 2016; Herman, 2017). Several transcriptional factors that include Leafy Cotyledon1 (LEC1) (Mu et al., 2008), Leafy Cotyledon2 (LEC2) (Stone et al., 2008), ABA-Insensitive3 (ABI3) (Giraudat et al., 1992), FUSCA3 (FUS3) (Wang and Perry, 2013), Transparent Testa2 (TT2) (Wang et al., 2014), and Wrinkled1 (WRI1) (Liu et al., 2010) were shown to regulate lipid biosynthesis and seed oil accumulation. In addition, a few studies have demonstrated the maternal effects on seed oil content, such as Apetala2 (AP2) (Jofuku et al., 2005) and Transparent Testa8 (TT8) (Chen et al., 2014) in Arabidopsis.

Molecular genetic studies have revealed a large number of loci involved in controlling B. napus oil content, but the identification of candidate genes from QTLs controlling seed oil content is rare because of the presence of large intervals and small effect of each QTL on oil accumulation (Rahman et al., 2013; Hua et al., 2016). The omics strategy such as transcriptome studies has identified a large number of genes associated with seed oil accumulation (Weselake et al., 2009; Tan et al., 2011; Huang et al., 2017; Tang et al., 2021). By combing comprehensive genome- and transcriptome-wide association studies using 505 inbred lines, Tang et al. (2021) identified two seed-oil-content QTLs as two homologous putative methyltransferase BnPMT6s and validated their functions using genetic studies on CRISPR/Cas9 mutants of each gene.

Seed oil accumulation in B. napus is accomplished by a complex array of physiological and biological processes, and members of late embryogenesis abundant (LEA) proteins are associated in this process. There are three main stages during seed oil accumulation in B. napus. The first stage is endosperm development. At this stage, considerable soluble sugars accumulate in the developing seeds mainly for osmotic regulation for seed expansion, substrate preparation for oil and protein biosynthesis, and signaling molecules for oil and other secondary metabolism (Lorenz et al., 2014). The second stage is embryo development. Approximately 20 days after anthesis, the embryo initiates, and then, the endosperm is quickly and thoroughly absorbed by the embryo. Then, the seed is full of the embryo and a layer of seed coat surrounding the embryo (Schwender and Ohlrogge, 2002). Most of the seed storages such as oil and proteins are mainly synthesized at this stage (Schwender and Hay, 2012). The third stage is seed dehydration and maturation. During seed maturation, considerable water in the seed will be dehydrated because excessive water will affect seed vitality. In B. napus seed, the oil does not exist as the formation of free fatty acids but in the form of oil body. The synthesized fatty acid normally is transported into endoplasmic reticulum membrane for further folding into oil body (He and Wu, 2009). The process is always accompanied by dehydration although the exact functions of dehydration are not very clear (Quettier and Eastmond, 2009). In previous investigations, a large family of small and highly hydrophilic proteins named late embryogenesis abundant (LEA) proteins was identified in plant tissues such as in maturating seeds (Liang et al., 2016, 2019). LEA proteins are mostly associated with the drought and freezing tolerance in plants (Liang et al., 2019) and are suggested to stabilize membranes and macromolecules by acting as water-binding molecules (Wang et al., 2003; Chakrabortee et al., 2012). But the exact functions of LEA proteins are largely unknown, especially during seed development and maturation in B. napus.

With the rapid improvement in the sensitivity, reproducibility and high throughput of liquid chromatography with tandem mass spectrometry (LC-MS/MS) -based quantitative proteomics, quantitative proteomics has become an important tool in identifying candidate proteins in improving plant productivity and stress responses (Li et al., 2011; Hu et al., 2015; Zhao et al., 2017; Wu et al., 2021; Zhu et al., 2022). In this study, TMT-based quantitative proteomics was used to uncover differentially accumulated proteins (DAPs) between two recombinant breeding lines differing in oil content, and an LEA protein BnLEA57 was found to be about two-fold in abundance in the higher oil content line (HOCL) compared to the low oil content line (LOCL). The function of the BnLEA57 gene and its close homolog BnLEA55 gene in seed oil accumulation was further analyzed through the ectopic expression in Arabidopsis thaliana. The result shed new insights on elucidating the molecular function of the LEA proteins on the oil biosynthesis in B. napus.

Results

Identification of Differentially Accumulated Proteins in Developing Seeds Between High Oil Content Line and Low Oil Content Line

To identify DAPs in the developing seeds of HOCL and LOCL in B. napus, seeds at 35 days after anthesis (DAA), which were at the critical stage of seed oil accumulation (Hua et al., 2014), were sampled for proteomic analysis. Tandem mass tag (TMT)-based quantitation was used to compare the seed proteome between HOCL and LOCL, following the general workflow in Figure 1. A total of 5,106 B. napus proteins were identified with FDR at 1% protein level, and 4,638 proteins were quantified across the three biological replicates for both HOCL and LOCL (no imputation was performed) (Supplementary Table 1). The principal component analysis showed that the three replicates of HOCL and LOCL were clustered together, respectively, but both clusters were well separated in component I, indicating the reproducible differences present in the proteome between both lines (Figure 2A). The data of DAPs were further analyzed according to the fold change and adjusted p-value. The result of volcano plot (Figure 2B) showed that 342 proteins were significantly different between HOCL and LOCL with a fold change >1.2 or <0.83 and adjusted p-value less than 0.05. Among those proteins, there were 199 proteins up-accumulated and 143 proteins down-accumulated in the HOCL compared to the LOCL (Supplementary Table 1).

FIGURE 1

Figure 1. Workflow of TMT labeling-based proteome profiling in B. napus seeds. Total proteins were extracted from 35-DAA seeds of two recombinant inbred lines of the B. napus with high oil content (HOCL) and low oil content (LOCL) (each line with three biological repeats), digested with trypsin after protein quality was checked on SDS-PAGE. The peptides were individually labeled with 6-plex TMT reagent (LOCL: 126, 127, and 128 labels; HOCL: 129, 130, and 131 labels) before combined, fractionated, and analyzed by LC-MS/MS.

FIGURE 2

Figure 2. Principal component analysis and volcano plot of the proteome in seeds of HOCL and LOC. (A) Principal component analysis of the proteome in seeds of HOCL and LOCL. Percentages (%) indicate the percentage of overall variance captured by each principal component. PC1: principal component 1; PC2: principal component 2. (B) Volcano plot of DAPs between the seeds of HOCL and LOCL. Threshold cutoff (>1.2 or <0.83) determined for log2 fold change ratios are represented by two vertical blue lines. Adjusted p-value cutoff of 0.05 is represented by a horizontal blue line. Red dots represent up-accumulated proteins, and green dots indicate down-accumulated proteins, whereas back dots represent proteins that do not meet the threshold or adjusted p-value, in HOCL compared to LOCL. The red cross indicates BnLEA57 protein.

Gene Ontology Enrichment Analysis of the Differentially Accumulated Proteins

The functions of the DAPs were further analyzed by Gene Ontology (GO) enrichment. The DAPs were classified by GO terms according to three categories, namely, biological process, molecular function, and cellular component. In biological process, the largest functional enriched groups were lipid biosynthesis and metabolism (Figure 3A). Among lipid metabolic proteins, two biotin carboxyl carrier proteins of acetyl-CoA carboxylases (BCCPs) were significantly up-accumulated in the seeds of HOCL (Table 1). There were two acyl carrier proteins downregulated whereas one upregulated in the seed cells of HOCL. A protein encoding beta-ketoacyl-(acyl-carrier-protein) synthase (KAS) was also markedly upregulated (Table 1). Another important type of the DAPs was related to photosynthesis. The proteins were involved in the whole process of photosynthesis in the seed cells including antenna proteins (LHCA1 and LHCB5), photosystem II proteins (PsbE, PsbQ, and PsbR), photosystem I proteins (PsaL and PsaN), ATPase (AtpH), and ribulose bisphosphate carboxylase (RbcL). Importantly, all the proteins except two ribulose bisphosphate carboxylase small subunits were significantly higher in the seed cells of HOCL (Table 1). For the DAPs of essential components of the photosynthetic apparatus, all of the proteins involving in isoprenoid biosynthesis were upregulated except one zeta-carotene desaturase (Supplementary Table 2). Amazingly, all the DAPs associated with cell wall modification, the polysaccharide metabolism, pectin catabolic–metabolic process, and galacturonan metabolic process were the members of pectinesterase family, suggesting a role of cell wall modification on oil accumulation (Supplementary Table 2).

FIGURE 3

Figure 3. Gene Ontology enrichment analysis for the differentially accumulated proteins between HOCL and LOCL. Gene Ontology enrichment analysis was categorized according to biology process [BP, (A)], cellular component [CC, (B)], and molecular function [MF, (C)].

TABLE 1

Table 1. Selected DAPs in immature seeds of 35DAA between HOCL and LOCL.

For cellular component, GO-enrichment-identified terms include extracellular region and chlorophyll thylakoid membrane. DAPs associated with chlorophyll thylakoid membrane term were mostly involved in photosynthesis, and all the 13 DAPs were upregulated in HOCL (Supplementary Table 2). A total of 18 DAPs were extracellular region proteins, 11 of which were up-accumulated in HOCL. Functions of these extracellular region proteins in seed oil accumulation are mostly unexplored, which include trypsin inhibitor, pectin acetylesterase, and peroxidase.

For the molecular function terms, the largest enriched group in DAPs was the nutrient reservoir activity, which included 12 DAPs (Table 1), nine of which were up-accumulated in the HOCL. These proteins included seed storage proteins on the endoplasmic reticulum and are the members of the cruciferin (11S globulin) or napin (1.7–2S albumin) family. The second largest group was enzyme inhibitor activity which included five pectinesterase, three pectinesterase inhibitors, one Kunitz-type soybean trypsin inhibitor (STI), and one cysteine proteinase inhibitor, four of which showed up-accumulation in HOCL (Supplementary Table 2).

Immunoblot Analysis of BnLEA57 in Developing Seeds of B. napus

Among the DAPs, one late embryogenesis abundant protein, BnLEA57 (BnaC05g37670D), was chosen for further analysis because of its much higher abundance in the HOCL (the expression ratio between HOCL and LOCL was 1.91) and the implication of LEA proteins in dehydration stress response (Liang et al., 2016). Its accumulation in the developing seeds in the two lines was further examined by immunoblot analysis. The result showed that BnLEA57 protein abundance rose from 35 to 45 DAA and sharply increased from 45 to 55 DAA in the developing seeds of both lines. In addition, BnLEA57 protein accumulation in 35- and 45-, DAA seeds was consistently higher in HOCL than in LOCL (Figure 4), supporting the proteomic analysis results.

FIGURE 4

Figure 4. Protein abundance of BnLEA57 in different seed development stages and in different lines (HOCL and LOCL). Seed protein samples were probed with an anti-BnLEA57 antibody (top lane). Coomassie blue-stained RuBisCO large subunit (RbcL) was used as a loading control (bottom lane).

Overexpression of the BnLEA57 and BnLEA55 Increased Seed Oil Content in Arabidopsis

To further characterize the BnLEA57 function on seed oil deposition, we overexpressed BnLEA57 (BnaC05g37670D) and its closest B. napus homolog BnLEA55 (BnaA05g23860D) in Arabidopsis plants, respectively. Less sequence similarity was observed among the 108 known LEA proteins in B. napus (Liang et al., 2016), but BnLEA57 and BnLEA55 that shared 96.1% of sequence identity (Figure 5) were clustered together in a subgroup in the LEA_4 family. BnLEA57 and BnLEA55 shared amino acid sequence identity of 84.3 and 86% with their ortholog LEA76 (AT3G15670) in Arabidopsis, respectively (Figure 5). Overexpression of either BnLEA57 or BnLEA57 in Arabidopsis increased seed oil content in Arabidopsis, which is consistent with the increased LEA57/LEA57 transcript levels compared with wild type (Figures 6A,B). The mean oil content was 38.9% in wild type seeds, 41.8% in 35S::BnLEA57 seeds, and 41.4% in 35S::BnLEA55 seeds. The results revealed the positive regulation of seed oil content by both BnLEA57 and BnLEA55 genes.

FIGURE 5

Figure 5. Alignment of the amino acid sequences of BnLEA57 (BnaC05g37670D), BnLEA55 (BnaA05g23860D) and their Arabidopsis ortholog AtLEA76 (AT3G15670). Residues identical in all the three LEAs are denoted with an asterisk (*), residues shared by two LEAs are indicated with a dot (.) or a colon (:).

FIGURE 6

Figure 6. Overexpression of BnLEA57 or BnLEA55 gene in Arabidopsis increased seed oil content. (A) Seed oil content in the Col-0 (WT) and different transgenic lines overexpressing BnLEA57 or BnLEA55 gene. (B) Relative gene expression level of BnLEA57 and BnLEA55 in different Arabidopsis transgenic lines. Error bars indicate standard deviation (n = 3). Significant difference between the transgenic plants and wild type (Col-0) is indicated: * for p < 0.05, ^** for p < 0.01, and ^*** for p < 0.001 for Student’s t-test significance.

Discussion

Fatty acid biosynthesis and oil body formation in B. napus seeds are the intricate and complicated networks of physiological and biochemical processes, and the underlying molecular mechanism remains to be further unveiled and characterized. In this study, TMT-based quantitative proteomics was applied to identify the proteins mediating seed oil accumulation via a pair of B. napus lines with contrasting seed oil content. GO enrichment analysis of DAPs suggests that proteins involved in lipid metabolism, photosynthesis, and nutrient reservoir activity are associated with contrasting seed oil accumulation between HOCL and LOCL (Figure 3). In addition, a role of LEA protein BnLEA57 and its homolog BnLEA55 in seed accumulation was further supported by their ectopic expression in Arabidopsis.

Our proteomic study suggests that higher photosynthesis capacity in the seeds contributes to higher lipid biosynthesis and oil accumulation in HOCL. Previous physiological studies comparing HOCL with LOCL, the same inbred lines used in this study, indicate higher chlorophyll content in the siliques of HOCL, supporting that higher photosynthetic capacity in siliques contribute to increased oil content in HOCL seeds. Results from this proteomic study support that photosynthesis from seeds could also contribute to oil accumulation in seeds. Among the DAPs, upregulation of protein abundance of two photosynthesis marker genes (chlorophyll A/B-binding proteins) (Martin et al., 2002) and two fatty acid biosynthesis marker genes (BCCP) (Ruuska et al., 2002) were observed (Table 1). In addition, among DAPs, GO-enriched terms in biological process include photosynthesis (9 of 11 proteins were up-accumulated in HOCL), and the enriched terms in the cellular process include thylakoid membranes (of which 13 proteins were all up-accumulated in HOCL, also involved in photosynthesis), supporting that higher photosynthesis activity contributes to higher seed oil content in HOCL. During fatty acid biosynthesis, considerable carbohydrates as substrate are required. Although it was estimated that about 70% of the photoassimilates were from silique wall (King et al., 1998), the remains were from transportation from other tissues such as leaf and stem and from developing seeds itself. A few lines of evidence showed that photosynthetic organelles were observed and the photosynthetic capability was found in developing seed cells (Eastmond et al., 1996; King et al., 1998; Ruuska et al., 2004; Houston et al., 2009). In addition, photosynthesis in seeds could enhance seed oil biosynthesis by increasing the carbon utilization efficiency (Schwender et al., 2004; Hay and Schwender, 2011) and by providing adenosine triphosphate (ATP), reductants NADPH, and NADH, which are the co-factors required for biosynthesis of fatty acids in chloroplasts (Goffman et al., 2005).

Among those DAPs, we found a late embryogenesis abundant protein BnLEA57 that was up-accumulated in the HOCL and showed that BnLEA57 and its homolog BnLEA55 positively regulate seed oil accumulation. It is known that LEAs are involved in the abiotic stress such as water deficiency and salt resistance (Dalal et al., 2009; Cuevas-Velazquez et al., 2014), but the exact molecular mechanism is poorly understood. In B. napus, Liang et al. (2019) found that overexpression of members of LEA_3 family (BnLEA3.1, BnLEA3.2, BnLEA3.3, and BnLEA3.4) increased oil content in Arabidopsis and B. napus seeds and suggested that LEA proteins increase seed oil accumulation by enhancing photosynthetic efficiency and improving membrane stability. The enhancement of photosynthetic capacity may result in the increment of substrate for fatty acid biosynthesis, namely, carbohydrate. The result was in accordance with our proteomic analysis. The LEA_4 family proteins BnLEA57 and BnLEA55 in our experiment were also shown to improve seed oil content when overexpressed in Arabidopsis. There were 108 LEA genes identified in B. napus genome, and many of them showed increased expression at late-stage seeds (Liang et al., 2016). Further investigations are required to reveal the accurate molecular mechanism of BnLEA57 and BnLEA55 on seed oil accumulation and to characterize whether different members of LEA proteins work individually or together on seed maturity and oil accumulation.

Materials and Methods

Plant Materials and Sampling

A total of two recombinant inbred lines of B. napus with high oil content (50.4%, HOCL) and low oil content (41.4%, LOCL) lines, two recombinant inbred lines at the F8 generations derived from parental lines between Huyou 15 as male parents and Zheshuang 6 as female parents (Hua et al., 2012, 2014), were grown in the field on the campus of Zhejiang Academy of Agricultural Sciences as described (Hua et al., 2014). The seeds for western blot assay were collected at 35, 45, and 55 days after anthesis on the same position of main inflorescence, the samples were immediately frozen in liquid nitrogen and kept at –80°C until use, and the samples collected at 35 DAA were also used for proteomic analysis.

Arabidopsis thaliana ecotype Columbia-0 was used for transformation. Arabidopsis seeds were surface-sterilized with 20 mg g^–1 NaClO solution for 10 min and washed five times by sterilized distilled water, sown on agar medium (half-strength Murashige and Skoog basal salt mixture, 15 g kg^–1 sucrose, pH 5.7), and vernalized at 4°C for 3 days in the dark. After germination, seedlings were transferred into soil and then grown in a growth room at 20–22°C, 16-h/8-h (day/night) photoperiod, and 60–70% humidity. The seeds were harvested as matured for oil content analysis. A total of three biological replicates for each sample were taken for analysis.

Protein Extraction and Digestion for Liquid Chromatography With Tandem Mass Spectrometry Analysis

Total proteins were extracted from the seeds of the B. napus, as previously described (Deng et al., 2007). Washed and pelleted protein sample was resuspended in 8 M urea solution, and the protein concentration was determined by the Bradford method.

For trypsin digestion, the samples were reduced with 10 mM DTT for 1 h at 56°C and alkylated with 55 mM iodoacetamide for 45 min at room temperature in darkness. Proteins were precipitated with 5 volumes of cold acetone overnight and then dissolved in 50 mM TEAB. Proteins were digested with trypsin overnight at 37°C in a 1:50 trypsin-to-protein mass ratio.

Tandem Mass Tag Labeling and High Performance Liquid Chromatography Fractionation

After trypsin digestion, peptides were desalted by Strata X C18 SPE column (Phenomenex) and vacuum-dried. Peptides were reconstituted in 50 mM TEAB and then labeled with respective 6-plex TMT reagents (LOCL: 126, 127, and 128 tags; HOCL: 129, 130, and 131 tags). The samples were pooled and fractionated by high pH reversed-phase high performance liquid chromatography (HPLC) using Agilent 300Extend C18 column (5-μm particles, 4.6 mm ID, 250 mm length). Peptides were separated with a gradient of 2–60% acetonitrile in 10 mM ammonium bicarbonate (pH 8) over 80 min. The peptides were then combined into 14 fractions and dried by vacuum centrifuging.

Liquid Chromatography With Tandem Mass Spectrometry Analysis

Multiplexed TMT-labeled samples were dissolved in 0.1% formic acid (solvent A), directly loaded onto a reversed-phase column (360 μm OD × 75 μm ID, 25 cm length) packed in-house with 3-μm C18 beads (ReproSil-Pur C18-AQ, Dr. Maisch) and eluted with a linear gradient of 6–25% solvent B (0.1% formic acid in 98% acetonitrile) for 28 min and 25–35% solvent B for 6 min at a constant flow rate of 300 nl/min on an EASY-nLC 1000 UPLC system. The resulting peptides were analyzed by Q Exactive Plus hybrid Quadrupole-Orbitrap mass spectrometer (Thermo Fisher Scientific, Waltham, MA, United States). A data-dependent procedure that alternated between one MS scan at a resolution of 70,000 and AGC target 3E6 followed by 20 MS/MS scans was applied for the top 20 precursor ions above a threshold ion count of 5E4 in the MS survey scan with 30.0-s dynamic exclusion. For MS scans, the m/z scan range was 350–1,600 Da. MS/MS scans were performed at the resolution at 17,500 and NCE at 30 with 3% stepped energy. The fixed first mass was set at 100 m/z for TMT quantification.

Data Analysis

The MS/MS raw data were analyzed using Proteome Discoverer (Version 2.4.0.305, Thermo Fisher Scientific, Waltham, MA, United States). Peptide and protein identifications were performed using both SEQUEST and MSFragger search engines against Brassica napus proteome from UniprotKB database (total 62,904 entries, as of 8 August 2021). Both searches were performed using the following settings: enzyme as trypsin/P, carbamidomethyl (C), TMT6plex (N-term), and TMT6plex (K) were selected as fixed modifications and oxidation (M) as variable modifications, 2 missed cleavage allowed and precursor error tolerance at 10 ppm, fragment deviation at 0.02 Da. Principal component analysis (PCA) of the samples was performed with R package PCAtools¹. R package limma was used in screening the differentially accumulated proteins (DAPs) between HOCL and LOCL (Phipson et al., 2016). Adjusted p-value < 0.05 and fold change >1.20 or <0.83 were chosen as the cutoff threshold. The volcano plot was plotted using EnhancedVolcano package from R². Gene Ontology (GO) enrichment for proteins that met differential expression criteria was performed using clusterProfiler in R environment (Wu et al., 2021).

Protein Immunoblot Analysis

Immunoblot analysis was performed as described (Deng et al., 2014). The proteins were extracted and dissolved in 2-D DIGE buffer (6 M urea, 2 M thiourea, and 4% CHAPS) and quantified using the Bio-Rad protein assay. Then, proteins were electrophoresed on 12% SDS-PAGE gels in a Tris-glycine buffer system at a constant 120 V/gel. Proteins were transferred to a nitrocellulose membrane at room temperature in a Bio-Rad semidry transfer System (Bio-Rad Laboratories, Hercules, CA, United States) according to the manufacturer’s instructions. The transferred membrane was immunoblotted with anti-BnLEA57 antibody (Abmart, Shanghai, China) and then incubated with goat-anti-rabbit IgG-horseradish peroxidase (HRP) secondary antibody before development with ECL SuperSignal West Dura Extended Duration Substrate (Thermo Fisher Scientific, Waltham, MA, United States), and the intensities were captured by the ImageQuant LAS 4000 mini system (GE Healthcare Life Sciences, Piscataway, NJ, United States). Anti-BnLEA57 antibody was produced using the synthetic peptide C-TTQQKAQETAQATK (Abmart, Shanghai, China).

Vector Construction and Plant Transformation

cDNA sequences of BnLEA57 (BnaC05g37670D) and BnLEA55 (BnaA05g23860D) were downloaded from the Brassica napus genome database³. For both cloning constructs (35S::LEA57 and 35S::LEA55), target sequences were first cloned into the pENTR/D-TOPO vector (Invitrogen, Carlsbad, CA, United States) before being recombined into an destination vector pB7WGF2 using the gateway LR reaction (Invitrogen, Carlsbad, CA, United States). The following primers were used: LEA57-F, 5′-CACCATGGCGTCTAACCAACAGAGC-3′; LEA57-R, 5′-CC TCTGATAAGTATGATGAGTAGGA-3′; LEA55-F, 5′-CACCA TGGCGTCCAACCAACAA-3′; and LEA55-R, 5′-CTTCCTCTG ATAAGTATGATGAGTCG-3′. Genetic transformation was performed as described (Zhou et al., 2011). The plasmids were transformed into Agrobacterium tumefaciens GV3101 strain and then transformed into Columbia-0 genotype using the floral-dip method (Clough and Bent, 1998). Transgenic seedlings were first selected on agar growth media with appropriate antibiotics and then confirmed by PCR using the corresponding primers.

RNA Extraction and Real-Time PCR Quantification

Total RNAs were extracted from the 10-day-old Arabidopsis seedlings with TRIzol Reagent (Invitrogen, Carlsbad, CA, United States), and cDNA was synthesized from 2 μg of total RNA using Moloney murine leukemia virus reverse transcriptase (Promega) and random primers in a 15 μl reaction according to the manufacturer’s instructions. For real-time PCR, the cDNAs were diluted to 4 times, and 1 μl of the cDNAs was added to 5 μl of PowerUP SYBR Green PCR Master Mix (Thermo Fisher Scientific, Waltham, MA, United States) and 0.4 μl of each primer (100 nM final concentration) in 10 μl reactions. PCR amplification and detection were performed using a Bio-Rad CFX96 Tough Real-time PCR Detection System using the following cycling conditions: 50°C for 2 min, 95°C for 2 min, followed by 40 cycles of 95°C for 15 s, and 60°C for 1 min. Melting-curve analysis was used to confirm the absence of non-specific amplification products. UBQ5 transcripts were used as an endogenous control to normalize expression of other genes. Relative expression levels were calculated by subtracting the threshold cycle (Ct) values for UBQ5 from those of the target gene and then calculating 2^–ΔCt as described previously (Zhou et al., 2011). RT-PCR experiments were performed on at least three independent samples. The following primers were used: UBQ5-F, 5′-TCGACGCTTCATCTCGTCCT-3′; UBQ5-R, 5′-CGCTGAACCTTTCCAGATCC-3′; LEA57-QPCR-F, 5′-ACAAAGGAGACGGCTCAAGA-3′; LEA57-QPCR-R, 5′-TTCCAACAGTGTGCTTCACC-3′; LEA55-QPCR-F, 5′-ACAAAGGAGACGGCTCAAGA-3′; LEA55-QPCR-R, 5′-GGTCGTAGTTGTGCCTGGAT-3′.

Oil Content Measurement and Statistical Analysis

Oil content in B. napus seeds was determined as described (Hua et al., 2014). In this experiment, statistical analysis was carried out using two-tailed Student’s t-test with two-samples unequal variance (“*,” “^**,” and “^***” represents as p < 0.05, p < 0.01, and p < 0.001) for significance.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: ProteomeXchange, accession number PXD033185.

Author Contributions

ZD, SH, and ZZ designed the experiment and wrote the manuscript. ZZ, ZD, and JT performed the proteomic analysis including data analysis. ZZ performed the western blotting and created the transgenic Arabidopsis lines. SH, BL, and PH prepared the B. napus samples and analyzed the seed oil content. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Key Research and Development Project (2018YFD1000900), Zhejiang Science and Technology Major Program on Agricultural New Variety Breeding (2021C02064), and Zhejiang Key Laboratory of Digital Dry Land Crops (2022E10012), and fundings from State Key Laboratory for Managing Biotic and Chemical Threats to the Quality and Safety of Agro-Products, and Zhejiang Academy of Agricultural Sciences.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We thank the technical support of the proteomic analysis provided by Jingjie PTM Biolabs, Inc. (Hangzhou, China). The mass spectrometry proteomic data have been deposited to the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the iProX partner repository (Ma et al., 2019) with the dataset identifier PXD033185.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.907244/full#supplementary-material

Supplementary Table 1 | Lists of identified, quantified, and differentially accumulated Brassica napus proteins between seeds of HOCL and LOCL.

Supplementary Table 2 | Functional protein groups identified by GO enrichment analysis of DAPs between HOCL and LOCL seed proteome.

Footnotes

1. ^ https://github.com/kevinblighe/PCAtools
2. ^ https://github.com/kevinblighe/EnhancedVolcano
3. ^ https://plants.ensembl.org/Brassica_napus/Info/Index

References

Bai, S., Wallis, J. G., Denolf, P., Engelen, S., Bengtsson, J. D., Van Thournout, M., et al. (2020). The biochemistry of headgroup exchange during triacylglycerol synthesis in canola. Plant J. 103, 83–94. doi: 10.1111/tpj.14709

PubMed Abstract | CrossRef Full Text | Google Scholar

Chakrabortee, S., Tripathi, R., Watson, M., Schierle, G. S., Kurniawan, D. P., Kaminski, C. F., et al. (2012). Intrinsically disordered proteins as molecular shields. Mol. biosyst. 8, 210–219. doi: 10.1039/c1mb05263b

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, B., Zhang, G., Li, P., Yang, J., Guo, L., Benning, C., et al. (2020). Multiple GmWRI1s are redundantly involved in seed filling and nodulation by regulating plastidic glycolysis, lipid biosynthesis and hormone signalling in soybean (Glycine max). Plant Biotechnol. J. 18, 155–171. doi: 10.1111/pbi.13183

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, M., Xuan, L., Wang, Z., Zhou, L., Li, Z., Du, X., et al. (2014). TRANSPARENT TESTA8 inhibits seed fatty acid accumulation by targeting several seed development regulators in Arabidopsis. Plant Physiol. 165, 905–916. doi: 10.1104/pp.114.235507

PubMed Abstract | CrossRef Full Text | Google Scholar

Clough, S. J., and Bent, A. F. (1998). Floral dip: a simplified method for agrobacterium-mediated transformation of Arabidopsis thaliana. Plant J. 16, 735–743. doi: 10.1046/j.1365-313x.1998.00343.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Cuevas-Velazquez, C. L., Rendon-Luna, D. F., and Covarrubias, A. A. (2014). Dissecting the cryoprotection mechanisms for dehydrins. Front. Plant Sci. 5:583. doi: 10.3389/fpls.2014.00583

PubMed Abstract | CrossRef Full Text | Google Scholar

Dalal, M., Tayal, D., Chinnusamy, V., and Bansal, K. C. (2009). Abiotic stress and ABA-inducible group 4 LEA from Brassica napus plays a key role in salt and drought tolerance. J. Biotechnol. 139, 137–145. doi: 10.1016/j.jbiotec.2008.09.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, Z., Oses-Prieto, J. A., Kutschera, U., Tseng, T. S., Hao, L., Burlingame, A. L., et al. (2014). Blue light-induced proteomic changes in etiolated Arabidopsis seedlings. J. Proteome Res. 13, 2524–2533. doi: 10.1021/pr500010z

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, Z., Zhang, X., Tang, W., Oses-Prieto, J. A., Suzuki, N., Gendron, J. M., et al. (2007). A proteomics study of brassinosteroid response in Arabidopsis. Mol. Cell Proteomics 6, 2058–2071. doi: 10.1074/mcp.M700123-MCP200

PubMed Abstract | CrossRef Full Text | Google Scholar

Dupont, J., White, P. J., Johnston, K. M., Heggtveit, H. A., Mcdonald, B. E., Grundy, S. M., et al. (1989). Food safety and health effects of canola oil. J. Am. Coll. Nutr. 8, 360–375. doi: 10.1080/07315724.1989.10720311

PubMed Abstract | CrossRef Full Text | Google Scholar

Eastmond, P., Koláčá, L., and Rawsthorne, S. (1996). Photosynthesis by developing embryos of oilseed rape (Brassica napus L.). J. Exp. Bot. 47, 1763–1769. doi: 10.1093/jxb/47.11.1763

CrossRef Full Text | Google Scholar

Giraudat, J., Hauge, B. M., Valon, C., Smalle, J., Parcy, F., and Goodman, H. M. (1992). Isolation of the Arabidopsis ABI3 gene by positional cloning. Plant Cell 4, 1251–1261. doi: 10.1105/tpc.4.10.1251

PubMed Abstract | CrossRef Full Text | Google Scholar

Goffman, F. D., Alonso, A. P., Schwender, J., Shachar-Hill, Y., and Ohlrogge, J. B. (2005). Light enables a very high efficiency of carbon storage in developing embryos of rapeseed. Plant Physiol. 138, 2269–2279. doi: 10.1104/pp.105.063628

PubMed Abstract | CrossRef Full Text | Google Scholar

Hay, J., and Schwender, J. (2011). Computational analysis of storage synthesis in developing Brassica napus L. (oilseed rape) embryos: flux variability analysis in relation to 13C metabolic flux analysis. Plant J. 67, 513–525. doi: 10.1111/j.1365-313X.2011.04611.x

PubMed Abstract | CrossRef Full Text | Google Scholar

He, Y. Q., and Wu, Y. (2009). Oil body biogenesis during Brassica napus embryogenesis. J. Integr. Plant Biol. 51, 792–799. doi: 10.1111/j.1744-7909.2009.00851.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Herman, E. M. (2017). “Cell and molecular biology of seed oil bodies,” in Seed Development and Germination, eds J. Kigel and G. Galili (Milton Park: Routledge), 195–214. doi: 10.1201/9780203740071

CrossRef Full Text | Google Scholar

Houston, N. L., Hajduch, M., and Thelen, J. J. (2009). Quantitative proteomics of seed filling in castor: comparison with soybean and rapeseed reveals differences between photosynthetic and nonphotosynthetic seed metabolism. Plant Physiol. 151, 857–868. doi: 10.1104/pp.109.141622

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, J., Rampitsch, C., and Bykova, N. V. (2015). Advances in plant proteomics toward improvement of crop productivity and stress resistancex. Front. Plant Sci. 6:209. doi: 10.3389/fpls.2015.00209

PubMed Abstract | CrossRef Full Text | Google Scholar

Hua, S., Chen, Z. H., Zhang, Y., Yu, H., Lin, B., and Zhang, D. (2014). Chlorophyll and carbohydrate metabolism in developing silique and seed are prerequisite to seed oil content of Brassica napus L. Bot. Stud. 55:34. doi: 10.1186/1999-3110-55-34

PubMed Abstract | CrossRef Full Text | Google Scholar

Hua, S., Yu, H., Zhang, Y., Lin, B., and Chen, Z. (2012). Variation of carbohydrates and macronutrients during the flowering stage in canola (Brasscia napus L.) plants with contrasting seed oil content. Aust. J. Crop Sci. 6, 1275–1282. doi: 10.3316/informit.732097682474027

CrossRef Full Text | Google Scholar

Hua, W., Liu, J., and Wang, H. (2016). Molecular regulation and genetic improvement of seed oil content in Brassica napus L. Front. Agric. Sci. Eng. 3, 186–194. doi: 10.15302/J-FASE-2016107

CrossRef Full Text | Google Scholar

Huang, K. L., Zhang, M. L., Ma, G. J., Wu, H., Wu, X. M., Ren, F., et al. (2017). Transcriptome profiling analysis reveals the role of silique in controlling seed oil content in Brassica napus. PLoS One 12:e0179027. doi: 10.1371/journal.pone.0179027

PubMed Abstract | CrossRef Full Text | Google Scholar

Jofuku, K. D., Omidyar, P. K., Gee, Z., and Okamuro, J. K. (2005). Control of seed mass and seed yield by the floral homeotic gene APETALA2. Proc. Natl. Acad. Sci. U.S.A. 102, 3117–3122. doi: 10.1073/pnas.0409893102

PubMed Abstract | CrossRef Full Text | Google Scholar

King, S. P., Badger, M. R., and Furbank, R. T. (1998). CO₂ refixation characteristics of developing canola seeds and silique wall. Funct. Plant Biol. 25, 377–386. doi: 10.1071/PP97157

CrossRef Full Text | Google Scholar

Li, T., Xu, S. L., Oses-Prieto, J. A., Putil, S., Xu, P., Wang, R. J., et al. (2011). Proteomics analysis reveals post-translational mechanisms for cold-induced metabolic changes in Arabidopsis. Mol. Plant 4, 361–374. doi: 10.1093/mp/ssq078

PubMed Abstract | CrossRef Full Text | Google Scholar

Liang, Y., Kang, K., Gan, L., Ning, S., Xiong, J., Song, S., et al. (2019). Drought-responsive genes, late embryogenesis abundant group3 (LEA3) and vicinal oxygen chelate, function in lipid accumulation in Brassica napus and Arabidopsis mainly via enhancing photosynthetic efficiency and reducing ROS. Plant Biotechnol. J. 17, 2123–2142. doi: 10.1111/pbi.13127

PubMed Abstract | CrossRef Full Text | Google Scholar

Liang, Y., Xiong, Z., Zheng, J., Xu, D., Zhu, Z., Xiang, J., et al. (2016). Genome-wide identification, structural analysis and new insights into late embryogenesis abundant (LEA) gene family formation pattern in Brassica napus. Sci. Rep. 6:24265. doi: 10.1038/srep24265

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J., Hua, W., Zhan, G., Wei, F., Wang, X., Liu, G., et al. (2010). Increasing seed mass and oil content in transgenic Arabidopsis by the overexpression of wri1-like gene from Brassica napus. Plant Physiol. Biochem. 48, 9–15. doi: 10.1016/j.plaphy.2009.09.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Lorenz, C., Rolletschek, H., Sunderhaus, S., and Braun, H. P. (2014). Brassica napus seed endosperm - metabolism and signaling in a dead end tissue. J. Proteomics 108, 382–426. doi: 10.1016/j.jprot.2014.05.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, J., Chen, T., Wu, S., Yang, C., Bai, M., Shu, K., et al. (2019). iProX: an integrated proteome resource. Nucleic Acids Res. 47, D1211–D1217. doi: 10.1093/nar/gky869

PubMed Abstract | CrossRef Full Text | Google Scholar

Maeo, K., Tokuda, T., Ayame, A., Mitsui, N., Kawai, T., Tsukagoshi, H., et al. (2009). An AP2-type transcription factor, WRINKLED1, of Arabidopsis thaliana binds to the AW-box sequence conserved among proximal upstream regions of genes involved in fatty acid synthesis. Plant J. 60, 476–487. doi: 10.1111/j.1365-313X.2009.03967.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, T., Oswald, O., and Graham, I. A. (2002). Arabidopsis seedling growth, storage lipid mobilization, and photosynthetic gene expression are regulated by carbon:nitrogen availability. Plant Physiol. 128, 472–481. doi: 10.1104/pp.010475

PubMed Abstract | CrossRef Full Text | Google Scholar

Mu, J., Tan, H., Zheng, Q., Fu, F., Liang, Y., Zhang, J., et al. (2008). LEAFY COTYLEDON1 is a key regulator of fatty acid biosynthesis in Arabidopsis. Plant Physiol. 148, 1042–1054. doi: 10.1104/pp.108.126342

PubMed Abstract | CrossRef Full Text | Google Scholar

Phipson, B., Lee, S., Majewski, I. J., Alexander, W. S., and Smyth, G. K. (2016). Robust hyperparameter estimation protects against hypervariable genes and improves power to detect differential expression. Ann. Appl. Stat. 10, 946–963. doi: 10.1214/16-AOAS920

CrossRef Full Text | Google Scholar

Quettier, A. L., and Eastmond, P. J. (2009). Storage oil hydrolysis during early seedling growth. Plant Physiol. Biochem. 47, 485–490. doi: 10.1016/j.plaphy.2008.12.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Rahman, H., Harwood, J., and Weselake, R. (2013). Increasing seed oil content in Brassica species through breeding and biotechnology. Lipid Technol. 25, 182–185. doi: 10.1002/lite.201300291

CrossRef Full Text | Google Scholar

Ruuska, S. A., Girke, T., Benning, C., and Ohlrogge, J. B. (2002). Contrapuntal networks of gene expression during Arabidopsis seed filling. Plant Cell 14, 1191–1206. doi: 10.1105/tpc.000877

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruuska, S. A., Schwender, J., and Ohlrogge, J. B. (2004). The capacity of green oilseeds to utilize photosynthesis to drive biosynthetic processes. Plant Physiol. 136, 2700–2709. doi: 10.1104/pp.104.047977

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwender, J., Goffman, F., Ohlrogge, J. B., and Shachar-Hill, Y. (2004). Rubisco without the Calvin cycle improves the carbon efficiency of developing green seeds. Nature 432, 779–782. doi: 10.1038/nature03145

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwender, J., and Hay, J. O. (2012). Predictive modeling of biomass component tradeoffs in Brassica napus developing oilseeds based on in silico manipulation of storage metabolism. Plant Physiol. 160, 1218–1236. doi: 10.1104/pp.112.203927

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwender, J., and Ohlrogge, J. B. (2002). Probing in vivo metabolism by stable isotope labeling of storage lipids and proteins in developing Brassica napus embryos. Plant Physiol. 130, 347–361. doi: 10.1104/pp.004275

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, R. V., Somidi, A. K., and Dalai, A. K. (2015). Preparation and properties evaluation of biolubricants derived from canola oil and canola biodiesel. J. Agric. Food Chem. 63, 3235–3242. doi: 10.1021/jf505825k

PubMed Abstract | CrossRef Full Text | Google Scholar

Stone, S. L., Braybrook, S. A., Paula, S. L., Kwong, L. W., Meuser, J., Pelletier, J., et al. (2008). Arabidopsis LEAFY COTYLEDON2 induces maturation traits and auxin activity: implications for somatic embryogenesis. Proc. Natl. Acad. Sci. U.S.A. 105, 3151–3156. doi: 10.1073/pnas.0712364105

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, H., Yang, X., Zhang, F., Zheng, X., Qu, C., Mu, J., et al. (2011). Enhanced seed oil production in canola by conditional expression of Brassica napus LEAFY COTYLEDON1 and LEC1-LIKE in developing seeds. Plant Physiol. 156, 1577–1588. doi: 10.1104/pp.111.175000

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, S., Zhao, H., Lu, S., Yu, L., Zhang, G., Zhang, Y., et al. (2021). Genome- and transcriptome-wide association studies provide insights into the genetic basis of natural variation of seed oil content in Brassica napus. Mol. Plant 14, 470–487. doi: 10.1016/j.molp.2020.12.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, F., and Perry, S. E. (2013). Identification of direct targets of FUSCA3, a key regulator of Arabidopsis seed development. Plant Physiol. 161, 1251–1264. doi: 10.1104/pp.112.212282

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, W., Vinocur, B., and Altman, A. (2003). Plant responses to drought, salinity and extreme temperatures: towards genetic engineering for stress tolerance. Planta 218, 1–14. doi: 10.1007/s00425-003-1105-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Z., Chen, M., Chen, T., Xuan, L., Li, Z., Du, X., et al. (2014). TRANSPARENT TESTA2 regulates embryonic fatty acid biosynthesis by targeting FUSCA3 during the early developmental stage of Arabidopsis seeds. Plant J. 77, 757–769. doi: 10.1111/tpj.12426

PubMed Abstract | CrossRef Full Text | Google Scholar

Weselake, R. J., Taylor, D. C., Rahman, M. H., Shah, S., Laroche, A., Mcvetty, P. B. E., et al. (2009). Increasing the flow of carbon into seed oil. Biotechnol. Adv. 27, 866–878. doi: 10.1016/j.biotechadv.2009.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Woodfield, H. K., Fenyk, S., Wallington, E., Bates, R. E., Brown, A., Guschina, I. A., et al. (2019). Increase in lysophosphatidate acyltransferase activity in oilseed rape (Brassica napus) increases seed triacylglycerol content despite its low intrinsic flux control coefficient. New Phytol. 224, 700–711.

Google Scholar

Wu, T., Hu, E., Xu, S., Chen, M., Guo, P., Dai, Z., et al. (2021). clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. Innovation (N Y) 2:100141. doi: 10.1016/j.xinn.2021.100141

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, C., Wang, P., Si, T., Hsu, C. C., Wang, L., Zayed, O., et al. (2017). MAP kinase cascades regulate the cold response by modulating ICE1 protein stability. Dev. Cell 43, 618–629e5. doi: 10.1016/j.devcel.2017.09.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Z., An, L., Sun, L., Zhu, S., Xi, W., Broun, P., et al. (2011). Zinc finger protein5 is required for the control of trichome initiation by acting upstream of zinc finger protein8 in Arabidopsis. Plant Physiol. 157, 673–682. doi: 10.1104/pp.111.180281

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, Y., Qiu, W., Li, Y., Tan, J., Han, X., Wu, L., et al. (2022). Quantitative proteome analysis reveals changes of membrane transport proteins in Sedum plumbizincicola under cadmium stress. Chemosphere 287:132302. doi: 10.1016/j.chemosphere.2021.132302

PubMed Abstract | CrossRef Full Text | Google Scholar