Skip to main content

Genome-wide identification and functional characterization of LEA genes during seed development process in linseed flax (Linum usitatissimum L.)

Abstract

Background

LEA proteins are widely distributed in the plant and animal kingdoms, as well as in micro-organisms. LEA genes make up a large family and function in plant protection against a variety of adverse conditions.

Results

Bioinformatics approaches were adopted to identify LEA genes in the flax genome. In total, we found 50 LEA genes in the genome. We also conducted analyses of the physicochemical parameters and subcellular location of the genes and generated a phylogenetic tree. LuLEA genes were unevenly mapped among 15 flax chromosomes and 90% of the genes had less than two introns. Expression profiles of LuLEA showed that most LuLEA genes were expressed at a late stage of seed development. Functionally, the LuLEA1 gene reduced seed size and fatty acid contents in LuLEA1-overexpressed transgenic Arabidopsis lines.

Conclusion

Our study adds valuable knowledge about LEA genes in flax which can be used to improve related genes of seed development.

Peer Review reports

Background

Late embryogenesis abundant (LEA) proteins are widespread in multiple types of tissues of living organisms [1, 2]. These proteins have been observed in bacteria, cyanobacteria [3], fungi and animals [1, 3] but were first discovered in mature cotton seed by researchers in 1981 [4]. As the name implies, this protein accumulates during the late stage of seed maturation. Subsequent discoveries identified the protein in other plants, such as rice, Arabidopsis thaliana, maize [1, 5, 6], etc. [7,8,9]. In plants, LEA genes express in many different tissues, such as seeds, roots, stems, and buds [10], so their potential functions are not limited to the process of seed development. Scientists have identified that LEA proteins can be induced to express and function as protectants of proteins and membranes in unique ways when cells are under stress, in particular drought and desiccation. Most LEA proteins are low-weight molecules ranging in size from 10 to 30 kD.

Several classifications of LEA proteins have been identified according to different standards. A widely adopted classification sorts the LEA proteins into eight subgroups: LEA_1, LEA_2, LEA_3, LEA_4, LEA_5, LEA_6, dehydrin and seed maturation protein (SMP). This classification is based on the sequence homology and conserved motifs available in the Pfam database [2, 5]. Among the eight LEA subgroups, with the exception of a few atypical hydrophobic proteins in the LEA_2, LEA_3 and SMP subgroups [11], the proteins possess high contents of Arg/Lys, Glu, Ala, Thr and Gly [12]. All dehydrin proteins have K-segments that are rich with lysine, and some even have Y-segments or S-segments. These segments can exist in the form of tandem repeats [13]. Unlike other proteins, most LEA proteins that possess intrinsically disordered proteins (IDPs) have no three-dimensional structures [14, 15], which accords with their high hydrophilicity.

Seed development, a crucial part of the angiosperm life cycle, is regulated by a large intricate network involving multiple factors, including transcription, epigenes, hormones, peptides and sugar signaling regulators [16]. In general, seed development can be roughly divided into two phases, morphogenesis and maturation [17]. Of the latter phase, strong expression of LEA proteins is regarded as a clear indication of seed maturation [18, 19]. Previous studies indicate that LEA proteins might be related to seed longevity, desiccation tolerance, and viability [20,21,22,23]. A subset of LEA proteins are regulated by a network of transcription factors containing ABI3, ABI4, ABI5, EEL and DOG1, as evidenced by the down-regulation of LEA transcripts in abi3, abi5, leafy cotyledon1 and fusca3 mutants [18, 24]. The transcription factors LEC1, FUSCA3, and ABI3 are involved in fatty acid biosynthesis and lipid storage in seeds [25]. However, little evidence demonstrates that LEA proteins control seed traits directly or indirectly. This may be because most research has been focused on the contributions of LEA proteins to the tolerance of drought, heat, cold and other abiotic stresses [19, 26]. To our knowledge, only Liang et al. (2019) demonstrated that overexpression of LEA3 in Arabidopsis and Brassica napus enhanced seed, seed weight, and oil content [27]. Overall, our knowledge on how LEA proteins are involved in seed development and the lipid-regulated network still have many gaps to fill. Moreover, LEA proteins in every subfamily exhibit different functions, thus these potential functions are additional gaps of knowledge that need to be filled.

With the development of rapid sequencing technology, more and more plant genomic information has become available. In the last 20 years, many LEA proteins have been identified in different plant species, including rice [6], A. thaliana [5], maize [28], B. napus [29], sorghum [30], watermelon [7], and wheat (Triticum aestivum) [31]. Additionally, studies report most LEA proteins in plant species have many members, for example, the numbers of members are 51 in A, thaliana [5], 108 in B. napus [29], and 281 in wheat [32]; the relatively high numbers reflects their significant role in plants. However, still unknown are the precise functions of most LEA genes.

Flax (Linum usitatissimum L.), a self-pollinating annual herb, has a long history of domestication of 8000 years, originated in the Middle East, and now is widely distributed around the world [33]. Flax is classified into two types, fiber flax and linseed flax, based on how each are utilized. Current linseed flax varieties are able to accumulate up to 50% oil content in seeds, and the majority of the fatty acids are composed of palmitic acid (PAL; C16:0, ~ 6%), stearic acid (STE; C18:0, ~ 2.5%), oleic acid (OLE; C18:1,~ 19%), linoleic acid (LIO; C18:2, ~ 13%) and linolenic acid (LIN; C18:3, ~ 55%) [34]. Distinct from most oil-bearing crops, linseed contains a diversity of amino acids and vitamins and a much higher level of unsaturated fatty acids, in particular alpha-linolenic acid (ALA), which accounts for up to 64% of unsaturated fatty acids in flax seed oil [35]. The fatty acid ALA and its transformations such as DHA (docosahexaenioc acid), EPA (eicosapntemacnioc acid) are greatly benefited for people health care.

Because the genome sequence of flax is available for study [36], researchers can more easily identify LEA genes in flax. In this study, several LuLEA genes were identified in the flax genome. Gene structure and phylogenic analyses showed that the genes could be classified into eight subgroups. Additionally, we determined gene expression levels during the seed development process. Lastly, from among the LuLEA genes that expressed abundantly at the late maturation stages, we selected one LEA gene, LuLEA1, to transform into Arabidopsis. The LuLEA1-over-expression lines produced seeds reduced in size and fatty acid contents compared to those in the WT (wild type). Our results will not only help improve understanding of the LEA family in the flax genome, but also provide insights into LEA functions correlating with oil metabolism in flax.

Results

Identification of LuLEA gene families in the flax genome

Combining the methods of local BLAST with HMM, 50 LuLEA gene members of the LEA family were identified the flax genome (Table 1). These genes were named in order from LuLEA1 to LuLEA50. Based on the sequence homology and conserved motifs in the Pfam database, these LuLEA genes were divided into eight subfamilies, the LuLEA_1, LuLEA_2, LuLEA_3, LuLEA_4, LuLEA_5, LuLEA_6, dehydrin, and SMP subfamilies. Among the subfamilies, the dehydrin subfamily had the highest number of genes, 10. Following the dehydrin group were the LuLEA_1, LuLEA_2, LuLEA_3 subfamilies with 9, 8, and 8 genes respectively. The smallest subfamilies were LuLEA_4 and LuLEA_6 in which each had two gene members.

Table 1 LEA genes in the linseed flax genome and their sequence characteristics and physicochemical parameters

According to the chromosomal locations of LuLEA genes noted in the NCBI database, we generated distribution profiles of 49 LuLEA genes for analysis (Supplemental Fig. 1). It was clear that chromosome1 had the largest number of LuLEA genes up to 14. Other chromosomes had fewer than 6 LuLEA genes. Except for the LuLEA_6 and LuSMP subfamilies, other subfamilies had 1 to 3 members located on chromosome1. We further mapped the other 14 chromosomes of flax and found they had one to six LuLEA genes. For example, chr4, chr5 and chr7 had only one LEA gene on each chromosome, while chr11, chr12, chr13 and chr15 had two LEA genes on each chromosome (Supplementary Fig. 1).

The physicochemical parameters of these 50 LuLEA genes were attained using ExPASy. With the exception of one gene fragment (LuLEA17) being 14,123 bp in length, the LuLEA gene fragments ranged from 165 bp (LuLEA32) to 1708 bp (LuLEA46). A majority of the members encoded less than 300 amino acids. Members in the same subgroup displayed similar features. For example, as the members of the LuLEA_4 group, both LuLEA26 (398 aa) and LuLEA27 (497 aa) encoded remarkably large numbers of amino acids, while members of the LuLEA_3 group (LuLEA18LuLEA25) encoded relatively small numbers of amino acids ranging from 81 to 109 aa. Likewise, molecular masses had the same pattern as amino acid numbers. Approximately two-thirds of the LuLEA proteins had relatively low isoelectric points (pI < 7), which consisted of all proteins in the LuLEA_2, LuLEA_4, LuLEA_5, LuLEA_6, and LuSMP subfamilies and some proteins in the dehydrin subfamily. The remaining proteins, in particular, both LuLEA_1 and LuLEA_3 subfamilies had pI > 7, meanwhile, LuLEA41 and LuLEA42 in dehydrin subfamily also had pI > 7. One-tenth of LuLEA proteins had relatively high values of grand average of hydropathicity (GRAVY > 0), and all of these proteins belonged to the LuLEA_2 subfamily. The data indicated that most LuLEA proteins were hydrophilic, especially those in LuLEA_5, while those in LuLEA_2 were determined as the most hydrophobic, which is consistent with the idea of atypical. Predictions of subcellular location showed that nearly 80% of LuLEA proteins were located in the nucleus. Only the LuLEA6 protein belonging to the LuLEA_1 subfamily was predicted to have a high possibility of being located in the plasma membrane. Interestingly, half of the LuLEA_3 proteins may be found in the chloroplast, and the other half of these members may be found in the mitochondrion. Moreover, LuLEA11 protein was also predicted to be in the chloroplast, and all of LuLEA_6 proteins with LuLEA2 were likely distributed in extracellular spaces (Supplementary Table 1).

Biological evolution and gene structure analysis of LuLEA genes

To investigate the homology and similarity for the identified LuLEA genes, an unrooted phylogenetic tree was constructed based on the alignment of all LuLEA protein sequences (Fig. 1). These genes divided into eight main clades, and the eight subfamilies of LuLEA protein sequences shared very low similarity. In contrast, high similarity was observed between a considerable number of proteins paired at the end of the branches, such as LuLEA1 and LuLEA4, LuLEA12 and LuLEA16, LuLEA35 and LuLEA37, which suggests there were still quite a few LuLEA members belonging to the same subfamily existing fair homology.

Fig. 1
figure 1

Phylogenetic analysis of the LuLEA genes in flax. LuLEA gene families are grouped by different colors. The unrooted tree was generated with the full-length amino acid sequences of the 50 LuLEA proteins using ClustalW in MEGA6 software

The distribution of exons and introns in the genetic sequences of the LuLEA genes are shown in Fig. 2. Approximately all genes longer than 400 bp contained both exons and introns. Those genes lacking introns were found in three subfamilies: LuLEA_2, LuLEA_3 and LuLEA_4. Most genes having introns had only one intron. Also worth noting is that LuLEA17, which grouped into the LuLEA_2 subfamily, had the longest length than any other gene, up to 14 kb, and it also had the largest number of introns (4) and exons (5). Furthermore, the longest intron in LuLEA17 was up to 10 kb in length.

Fig. 2
figure 2

The distribution of exons and introns in LuLEA genes. Colored boxes represent exons, and grey lines indicate introns. The LEA subfamilies are distinguished by color. The lengths and positions of exons and introns in LuLEA genes are indicated by the scales at the bottom

In addition to a gene sequence structure analysis, the distribution of motifs of each protein sequence was analyzed (Fig. 3). A total of 50 LuLEA protein sequences were submitted to MEME tool to determine the characters of the motifs. In general, one to three motifs were found for each subfamily and the motifs differed greatly among subfamilies (Fig. 3). The LEA_4 and LEA_6 subfamilies had too few motifs in common with those of the other subfamilies and thus were not shown in the results. Nevertheless, much similarity was observed in the numbers and types of members with in the same subfamily, which reflects the credibility of the phylogenetic analysis. Remarkably, the dehydrin subfamily had plenty of conserved hydrophilic amino acids, such as G (Glycine) and K (Lysine), which implies a subfamily trait of hydropathy.

Fig. 3
figure 3

The distribution of motifs in LuLEA genes and the conserved amino acids in each subfamily exhibited by WebLogo plot. Different motifs are distinguished by different colored boxes. The maximum number of motifs in each sequence was set to 10. The representative motif of each subfamily is shown on the right

Gene expression pattern analysis of LuLEA genes during seed developing stages

The expression patterns showed that nearly all of the LuLEA genes expressed throughout all stages of seed development for both of our flax cultivars, Heiya No.14 and Macbeth. In Heiya No.14, a total of 42 LuLEA genes expressed during all stages, and 44 LuLEA genes expressed in Macbeth. In comparing the commonly expressed genes between these two cultivars, we found that 36 LuLEA genes expressed over 5 days, 10 days, 20 days 30 days after pollination (DAP); one gene express at the 30th day of seed development. Additionally, there was also only one gene that expressed at the 10th day in Macbeth but not in Heiya No.14, which signifies another difference between the two flax cultivars (Fig. 4a-b).

Fig. 4
figure 4

Expression profiles of LuLEA gene families in flax seed development. a Venn diagram of shared and non-shared numbers of genes of the cultivar Heiya No.14 expressed at 5, 10, 20, and 30 days after pollination; b similar to a but of the cultivar Macbeth; c comparison of LuLEA gene expression levels in every subfamily during seed maturation. Shades of blue color represent lower expression levels, and shades of red color represent higher levels

To confirm the observed variation in expression patterns among members in the LEA subfamilies, heat maps were produced for individual subfamilies. The trends of most LuLEA gene expression levels were consistent between Heiya No.14 and Macbeth. Some LuLEA genes, such as LuLEA1, LuLEA2 and LuLEA41, tended to highly express at late stages of seed development. On the contrary, expression of a few LuLEA genes, such as LuLEA15, LuLEA38, and LuLEA43, decreased from early to late developmental stages. Genes in the five main LuLEA subfamilies, LuLEA_1, LuLEA_4–LuLEA_6, and LuSMP, displayed similar expression patterns. An exception was observed in LuLEA32 where this gene highly expressed throughout our four sampling periods. The rest of the genes in the five subfamilies exhibited increased expression largely at days 20 and 30 (Fig. 4c).

LuLEA1 is responsible for seed development and fatty acid metabolism

Two independent overexpression lines, named LuLEA1–6 and LuLEA1–7, were generated and analyzed. Compared to those of the WT plants, both of the two overexpression lines had significantly lower values of the measured seed traits, seed weight, area and circumference. The results indicate that LuLEA1 may play a role in regulating seed size (Fig. 5a-c).

Fig. 5
figure 5

Effect of LuLEA1 on seed size and fatty acid content in transgenic Arabidopsis lines. a The 1000-seed weights of transgenic LuLEA1 lines and wild type grown under the “normal” condition; b the average area of transgenic and WT seeds; c the average circumference of transgenic and WT seeds. d Total fatty acid contents of overexpression lines and wild type grown under the “normal” condition; e 10 fatty acid components in overexpression lines and wild type grown under the “normal” condition; f proportions of each fatty acid component. The data represent means and standard deviations (STD) of at least three replicates. Statistically significant differences were determined by two-tailed paired Student’s t-tests comparing transgenic data with WT data of plants grown under the same condition. * p < 0.05; ** p < 0.01

To further explore whether LuLEA1 functions in controlling fatty acid metabolism during seed development, fatty acid content in the transgenic Arabidopsis lines were determined by GC-MS (Gas Chromatography-Mass Spectrometer). Total average fatty acid contents of the two overexpression lines were less than that of WT, and LuLEA1–6 was significantly reduced. Meanwhile, most contents of the different types of fatty acids of the transgenic lines were lower than those of WT, and the contents of C18:0, C18:3, C20:1 and C20:3 were significantly reduced. In addition, we found that the proportion of each fatty acid differed, too. The proportions of C18:1 and C18:2 in transgenic lines were markedly higher, while the proportions of C18:3 and C20:1 were lower than WT. These results suggest that LuLEA1 may block the process of transformations of C18:1 and C18:2 into C18:3 and C20:1 (Figs. 5d-f).

Discussion

The LEA gene family is a large and complicated family, having many members that belong to different subfamilies. Genes in the LEA family have been identified in many crops, such as rice, A. thaliana and wheat. Besides in plants, this family of genes has been reported in both animals and microorganisms. However, characterization and identification of the LEA protein family in flax has never been reported. In this research, 50 LEA genes were identified in the flax genome, nearly equal to the 51 LEA genes found in A. thaliana. Given that flax is diploid (2n = 30) and the number of LuLEA is close to that of A. thaliana, whole genome duplication events of LEA genes occurring in flax was doubtful and supported by many findings of the evolutionary conservation of LEA genes [22, 29]. Thus, it is easy to deduce that LEA must play a crucial role in the development of organisms.

The 50 LEA genes in flax were divided into eight subfamilies. Among the subfamilies, the dehydrin subfamily has the greatest number of genes, 10, in the LuLEA family, while the LuLEA_6 subfamily has the least with 2. The distributions of the LEA_6 and dehydrin genes in flax are similar to those in A. thaliana. From multiple plant species comparisons, although some are largely occupied with LEA_4 subfamily or LEA_2 subfamilies, such as A. thaliana, B. napus, cotton (Gossypium hirsutum), tea (Camellia sinensis), dehydrin subfamily tends to share considerably part, which means dehydrin is relatively conserved and likely to provide more stable protection for cells during the evolution. Evidence shows that the LEA_6 subfamily is not found in algal and rice genomes [6, 22], which suggests LEA_6 was extended from other ancient LEA genes, and probably makes contribution to struggling with the water loss.

Based on our results, five LuLEA_2 genes likely encode hydrophobic proteins, while the others are hydrophilic proteins. This result is consistent with the results of past research on A. thaliana, Populus trichocarpa, and Solanum tuberosum. LEA_2 genes are thought to be heterologous to other subfamilies of LEA genes, which may explain the unique structures, atypical characters and even novel functions reported of members in the LEA_2 subfamily [3].

In flax, LEA_2 genes may only be present in chloroplasts and mitochondria, which indicates that LEA_2 may function in protecting proteins in these particular cellular organelles. There were also some LuLEA proteins in the nucleus and cytoplasm as well as cytoplasmic membranes. These results indicate that LEA proteins are widely distributed within cells, so these proteins having an important role such as protection of cellular compartments during stressful conditions is not without support. Moreover, most of our identified LuLEA proteins are hydrophilic according to their GRAVY values, which is quite similar to characterizations determined of LEA proteins in other higher plants [5, 9, 29]. Many studies have shown that the trait of high hydrophilicity is attributable to the presence of IDPs in LEA proteins, and high hydrophilicity facilitates their potential functions as protein and membrane protectants and molecular chaperones to ensure cellular survival in a variety of adverse environments.

The map of gene structures containing introns and exons clearly show a large number of LuLEA genes possessed less than 2 introns and relatively short gene lengths. One previous study showed that genes associated with stress response have few introns [37], which is supported by our results. Reports of many LEA genes with few introns in other plant species confirm this as well. In B. napus, 16/108 BnLEA genes have no introns, and the subfamily BnLEA_6 has five members that each have only one intron [29]. In wheat, 62% of its LEA genes have no introns [32]. In A. thaliana, 66.7% of its LEA genes contain only one intron [5]. In addition, similar conclusions of low intron numbers have been reported in other genes known to be involved in stress responses. For example, most StHsp20 genes (89.6%) with no or only one intron were demonstrated to respond to multiple abiotic stresses [38]. In another example, a high percentage (83.9%) of the zinc finger homeodomain genes that encode transcription factors involved in plant development and abiotic stress response in B. napus lack introns [39]. From the perspective of biomolecular activities, introns will be spliced out of the final sequence after transcription. Reduced introns of genes are benefit for the faster process from transcription to expression, which is convenient for cell to make a reaction to abiotic stresses and decrease the cost for transcription [40].

Different LuLEA subfamilies have various motif distributions. Proteins belonging to the same subfamily have similar numbers and types of motifs, which is illustrated by our phylogenetic tree. Maybe these characters imply the reasons for various functions of LEA proteins.

In most cases, gene expression analysis can help reveal important functions of target genes. According to the expression pattern of LuLEA genes at 5, 10, 20, and 30 DAP, only six LuLEA genes lacked expression during linseed maturation, while the other genes expressed throughout the entire process. These observations suggest that these genes play vital roles in the seed maturation process. Additionally, expression of members in several subfamilies, such as LuLEA_1, LuLEA_4 ~ LuLEA_6 and LuSMP, accumulated in abundance in late seed maturation, which is consistent with the reported data of previous studies [18, 19, 26]. These LuLEA genes were speculated to play an important part in seed maturation and desiccation. Meanwhile, some LuLEA genes, such as LuLEA_3 and Ludehydrin appeared no clearly regularity, which may explain the diversity of potential functions of LEA proteins and the correlations to the various structures.

Past studies have shown that LEA genes participate in the regulatory network of seed development [18], thus we investigated the phenotypes of seeds produced from LuLEA1-overexpressing transgenic Arabidopsis. The traits of seed weight, area and circumference were all reduced. Furthermore, fatty acid contents in seeds also declined. Based on those results, we conclude that the LEA_1 subfamily of genes negatively regulate seed size and fatty acid contents. Interestingly, Liang et al. [27] showed the opposite result: overexpression of a gene belonging to the LEA_4 subfamily, BnLEA3, could increase seed size and seed oil content in Arabidopsis. However, there is no evidence indicating the direct involvement of LEA genes in the regulatory mechanism of seed size and oil synthesis. Based on existing findings, LEA proteins are regulated by transcription factors ABI3, ABI4, ABI5 [18], and these factors have also been shown to affect seed size and lipid biosynthesis [23, 41, 42]. Thus, LEA proteins likely have a feedback relationship with these transcription factors, and different LEA families may have contrasting functions conferred by their different subfamilies to maintain a balance among functions in collectively protecting a plant.

Conclusions

In this research, a total of 50 LEA genes were identified in the flax genome, and they were divided into eight subfamilies based on their conserved domains. Genes from the same subfamily had similar structures, which is also supported by the results of phylogenetic analysis. All LuLEA genes were distributed on each chromosome. The overexpression of LuLEA1 in Arabidopsis decreased the traits of seed weight and size, as well as fatty acid contents.

Methods

Identification of LEA gene family members in the flax genome

Fifty-one LEA gene sequences of A. thaliana were retrieved from the database TAIR (The Arabidopsis Information Resource, https://www.arabidopsis.org/), and then they were blasted using protein sequences of flax acquired from the genome database Phytozome (https://phytozome.jgi.doe.gov/pz/portal.html#!info?alias=Org_Lusitatissimum). We also used the Pfam database (https://pfam.xfam.org/search) and HMMER to search for the genes with the conserved LEA domain [43]. Combining BLAST with HMMER, the initial candidate LEA genes of flax were obtained after filtering the mismatched or redundant genes. Three website tools, CDD (Conserved Domain Database, https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/cdd/), Pfam and SMART (https://smart.embl-heidelberg.de/smart/set_mode.cgi?NORMAL=1) were used to confirm and ensure all candidate genes contained the LEA family domain. The final filtered genes were assigned new names in numbered order.

The number of amino acids and gene lengths were obtained through the Phytozome web portal (https://phytozome.jgi.doe.gov/pz/portal.html), and chromosome locations of the LuLEA genes were obtained from the NCBI database (National Center for Biotechnology Information, https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/). The physicochemical parameters, composed of molecular weight (kDa), GRAVY (grand average of hydropathy) and pI (isoelectric point), of each LuLEA protein were calculated by ExPASy (www.expasy.org/tools/). Subcellular location prediction was conducted using the BUSCA annotation system (https://busca.biocomp.unibo.it/).

Phylogenetic and sequence feature analysis of LuLEA family members

Multiple sequence alignment of 50 LuLEA protein sequences was performed using ClustalW [44], and these results were used to construct a phylogenetic tree with the MEGA7 software [45]. The method of maximum likelihood was adopted to construct the tree, and it had 1000 bootstrap replicates. To understand the structural features of LuLEA genes, the genetic sequences containing exons and introns were examined, and the distributions of motifs on each protein sequence were determined. Owing to the variation between each sequence, the maximum value of motif for each gene was set as 10. In the gene structure analysis of LuLEA genes, which was limited to the annotation of flax, UTRs (untranslated region) could not be displayed. The distribution of intron and exon fragments on each LuLEA gene were visualized by a diagram with the help of the Gene Structure Display Server (https://gsds.cbi.pku.edu.cn/). The relative locations of conserved amino acid motifs encoded by LuLEA family genes were determined using Multiple Expectation Maximization for Motif Elicitation tool (https://alternate.meme-suite.org/). The chromosomal locations of LuLEA genes were derived from the positional information available in the NCBI website. The distribution of LuLEA family members on the chromosomes were visualized using MG2C (https://mg2c.iask.in/mg2c_v2.0/).

RNA extraction and RNA-seq of developing seed samples

The flax cultivars Heiya No. 14 [46] and Macbeth were used as the plant materials for sample collection and RNA isolation. Heiya No.14 was bred for the purpose of better quality and high yields of fiber flax, and its seed oil content makes up about 25% of seed weight. Macbeth is an oilseed flax that produces about 40% seed oil content as well as large seed sizes. Plants were grown in a greenhouse under “normal” growth conditions of 24 °C and a 16 h daylight/8 h dark cycle. After plants reproduced, the siliques were collected at 5 days (DAP5), 10 days (DAP10), 20 days (DAP20), and 30 days after pollination (DAP30) and immediately frozen in liquid nitrogen before RNA isolation. Two replicates were prepared for the construction of a sequencing library per sample. Total RNA was isolated using TRIzol reagent (Invitrogen, 15,596–026), according to the manufacturer’s instructions. Then cDNA libraries were constructed and subsequently inspected. Based on sequencing by synthesis technology, the Illumina HiSeq 2500 platform was used to perform cDNA library sequencing and acquire a large amount of high-quality data.

Gene expression pattern analysis for LuLEA gene families with RNA-seq data

We used RNA-seq data to analyze the gene expression patterns of LuLEA genes. After filtering the sequenced raw data, the clean data were mapped to the flax reference genome (https://phytozome.jgi.doe.gov/pz/portal.html). Then, the FPKM (Fragments per Kilobase of Exon per Million Fragments Mapped) method [47] was applied to calculate gene expression levels based on the number of reads mapped to the reference sequence. A heatmap of gene expression profiles of all LuLEA genes was constructed using Mev 4.0 software [48] with Pearson’s correction and complete linkage clustering. The raw data have been submitted to the NCBI database with the GEO number GSE130378.

Vector construction, gene transformation, and phenotypic screening of transgenic plants

In order to test how LuLEA genes may affect plant development, we selected one LuLEA gene with high expression during late seed maturation for use in the genetic transformation of A. thaliana. The selected gene, LuLEA1, exhibited a level of expression at 30 DAP that was up to 10,000-fold that of the level at 5 DAP based on the RNA-Seq data. The RNA-Seq raw data is available in the NCBI database with the GEO number GSE130378 (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/geo/query/acc.cgi?acc=GSE130378). The full-length CDS of LuLEA1 was cloned into the CaMV 35S-Red vector. The plasmids were double digested with the restriction endonucleases XmaI and EcoRI and then ligated with the specific gene transcript fragment so that the gene expression of the target gene was under the control of the CaMV 35S promoter. The construct was transformed into Agrobacterium tumefaciens strain EHA105 using the freeze–thaw method. Arabidopsis Col-0 plants were then transformed using the floral dip method [49]. Untransformed Arabidopsis plants were used as WT controls. All plants were maintained in a greenhouse under standard conditions (24°Cday/18 °C night and 16 h light/8 h dark).

Transgenic plants were screened and cultivated to the T3 generation. Then the seeds were harvested, the size and weight of which were determined by a crop scanning test system (Wanshen SC-G, China) [27] and the Seed Count image analysis system [50]. And the fatty acid compositions in seed samples were quantified by gas chromatography mass spectrometry (GC-MS) [50].

Availability of data and materials

The raw RNA-seq data of cultivars Macbeth and Heiya No.14 of Linum usitatissimum L. obtained at different developmental stages of seeds are available in the NCBI database under the GEO number GSE130378 (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/geo/query/acc.cgi?acc=GSE130378). All data generated or analyzed during this study are included in this published article and its supplementary information files. The datasets used and analyzed for the current study are available from the corresponding author upon reasonable request.

Abbreviations

LEA:

Late embryogenesis abundant protein

SMP:

Seed maturation protein

ALA:

Alpha-linolenic acid

IDPs:

Intrinsically disordered proteins

PAL:

Palmitic acid

STE:

Stearic acid

OLE:

Oleic acid

LIO:

Linoleic acid

LIN:

Linolenic acid

DAP:

Day after pollination

GC-MS:

Gas Chromatography-Mass Spectrometer

References

  1. Hand SC, Menze MA, Toner M, Boswell LC, Moore DS. LEA proteins during water stress: not just for plants anymore. Annu Rev Physiol. 2011;73(1):115–34. https://0-doi-org.brum.beds.ac.uk/10.1146/annurev-physiol-012110-142203.

    Article  CAS  PubMed  Google Scholar 

  2. Hunault G, Jaspard E. LEAPdb: a database for the late embryogenesis abundant proteins. BMC Genomics. 2010;11(1):221. https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2164-11-221.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Tunnacliffe A, Wise MJ. The continuing conundrum of the LEA proteins. Naturwissenschaften. 2007;94(10):791–812. https://0-doi-org.brum.beds.ac.uk/10.1007/s00114-007-0254-y.

    Article  CAS  PubMed  Google Scholar 

  4. Dure L, Greenway S, Galau GA. Developmental biochemistry of cottonseed embryogenesis and germination: changing messenger ribonucleic acid populations as shown by in vitro and in vivo protein synthesis. Biochemistry. 1981;20(14):4162–8. https://0-doi-org.brum.beds.ac.uk/10.1021/bi00517a033.

    Article  CAS  PubMed  Google Scholar 

  5. Hundertmark M, Hincha DK. LEA (late embryogenesis abundant) proteins and their encoding genes in Arabidopsis thaliana. BMC Genomics. 2008;9(1):118. https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2164-9-118.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Wang X, Zhu H, Jin G, Liu H, Wu W, Zhu J. Genome-scale identification and analysis of LEA genes in rice (Oryza sativa L.). Plant Sci. 2007;172(2):414–20. https://0-doi-org.brum.beds.ac.uk/10.1016/j.plantsci.2006.10.004.

    Article  CAS  Google Scholar 

  7. Altunoglu YC, Baloglu MC, Baloglu P, Yer EN, Kara S. Genome-wide identification and comparative expression analysis of LEA genes in watermelon and melon genomes. Physiol Mol Biol Plants. 2017;23(1):5–21. https://0-doi-org.brum.beds.ac.uk/10.1007/s12298-016-0405-8.

    Article  CAS  Google Scholar 

  8. İbrahime M, Kibar U, Kazan K, Yüksel Özmen C, Mutaf F, Demirel Aşçı S, et al. Genome-wide identification of the LEA protein gene family in grapevine (Vitis vinifera L.). Tree Genet Genom. 2019;15:55.

    Article  Google Scholar 

  9. Wang W, Gao T, Chen J, Yang J, Huang H, Yu Y. The late embryogenesis abundant gene family in tea plant (Camellia sinensis): genome-wide characterization and expression analysis in response to cold and dehydration stress. Plant Physiol Biochem. 2019;135:277–86. https://0-doi-org.brum.beds.ac.uk/10.1016/j.plaphy.2018.12.009.

    Article  CAS  PubMed  Google Scholar 

  10. Hongbo S, Zongsuo L, Mingan S. LEA proteins in higher plants: structure, function, gene expression and regulation. Colloids Surf B Biointerfaces. 2005;45(3-4):131–5. https://0-doi-org.brum.beds.ac.uk/10.1016/j.colsurfb.2005.07.017.

    Article  CAS  Google Scholar 

  11. Wang M, Li P, Li C, Pan Y, Jiang X, Zhu D, et al. SiLEA14, a novel atypical LEA protein, confers abiotic stress resistance in foxtail millet. BMC Plant Biol. 2014;14(1):290. https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-014-0290-7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Battaglia M, Olvera-Carrillo Y, Garciarrubio A, Campos F, Covarrubias AA. The enigmatic LEA proteins and other hydrophilins. Plant Physiol. 2008;148(1):6–24. https://0-doi-org.brum.beds.ac.uk/10.1104/pp.108.120725.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Eriksson SK, Kutzer M, Procek J, Gröbner G, Harryson P. Tunable membrane binding of the intrinsically disordered dehydrin Lti30, a cold-induced plant stress protein. Plant Cell. 2011;23(6):2391–404. https://0-doi-org.brum.beds.ac.uk/10.1105/tpc.111.085183.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Wise MJ, Tunnacliffe A. POPP the question: what do LEA proteins do? Trends Plant Sci. 2004;9(1):13–7. https://0-doi-org.brum.beds.ac.uk/10.1016/j.tplants.2003.10.012.

    Article  CAS  PubMed  Google Scholar 

  15. Sun X, Rikkerink EHA, Jones WT, Uversky VN. Multifarious roles of intrinsic disorder in proteins illustrate its broad impact on plant biology. Plant Cell. 2013;25(1):38–55. https://0-doi-org.brum.beds.ac.uk/10.1105/tpc.112.106062.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Savadi S. Molecular regulation of seed development and strategies for engineering seed size in crop plants. Plant Growth Regul. 2018;84(3):401–22. https://0-doi-org.brum.beds.ac.uk/10.1007/s10725-017-0355-3.

    Article  CAS  Google Scholar 

  17. Locascio A, Roig-Villanova I, Bernardi J, Varotto S. Current perspectives on the hormonal control of seed development in Arabidopsis and maize: a focus on auxin. Front Plant Sci. 2014;5:412.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Olivier L, Anthoni P, Souha B, Julia B. Late seed maturation: drying without dying. J Exp Bot. 2017:827–41.

  19. Jin X, Cao D, Wang Z, Ma L, Li Y. Genome-wide identification and expression analyses of the LEA protein gene family in tea plant reveal their involvement in seed development and abiotic stress responses. Sci Rep. 2019;9(1):14123. https://0-doi-org.brum.beds.ac.uk/10.1038/s41598-019-50645-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Wu X, Liu H, Wang W, Chen S, Hu X, Li C. Proteomic analysis of seed viability in maize. Acta Physiol Plant. 2011;33(1):181–91. https://0-doi-org.brum.beds.ac.uk/10.1007/s11738-010-0536-4.

    Article  CAS  Google Scholar 

  21. Chatelain E, Hundertmark M, Leprince O, Gall SL, Satour P, Deligny-Penninck S, et al. Temporal profiling of the heat-stable proteome during late maturation of Medicago truncatula seeds identifies a restricted subset of late embryogenesis abundant proteins associated with longevity. Plant Cell Environ. 2012;35(8):1440–55. https://0-doi-org.brum.beds.ac.uk/10.1111/j.1365-3040.2012.02501.x.

    Article  CAS  PubMed  Google Scholar 

  22. Artur MAS, Zhao T, Ligterink W, Schranz E, Hilhorst HWM. Dissecting the genomic diversification of late embryogenesis abundant (LEA) protein gene families in plants. Genome Biol Evol. 2019;11(2):459–71. https://0-doi-org.brum.beds.ac.uk/10.1093/gbe/evy248.

    Article  CAS  PubMed  Google Scholar 

  23. Sano N, Rajjou L, North HM, Debeaujon I, Marion-Poll A, Seo M. Staying alive: molecular aspects of seed longevity. Plant Cell Physiol. 2016;57(4):660–74. https://0-doi-org.brum.beds.ac.uk/10.1093/pcp/pcv186.

    Article  CAS  PubMed  Google Scholar 

  24. Bies-Ethève N, Gaubier-Comella P, Debures A, Lasserre E, Jobet E, Raynal M, et al. Inventory, evolution and expression profiling diversity of the LEA (late embryogenesis abundant) protein gene family in Arabidopsis thaliana. Plant Mol Biol. 2008;67(1-2):107–24. https://0-doi-org.brum.beds.ac.uk/10.1007/s11103-008-9304-x.

    Article  CAS  PubMed  Google Scholar 

  25. Roscoe T. Devic, Martine. Seed maturation: Simplification of control networks in plants. Plant Sci. 2016;252:335–46.

    Article  PubMed  Google Scholar 

  26. Banerjee A, Roychoudhury A. Group II late embryogenesis abundant (LEA) proteins: structural and functional aspects in plant abiotic stress. Plant Growth Regul. 2016;79(1):1–17. https://0-doi-org.brum.beds.ac.uk/10.1007/s10725-015-0113-3.

    Article  CAS  Google Scholar 

  27. Liang Y, Kang K, Gan L, Ning S, Xiong J, Song S, et al. 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. 2019;17(11):2123–42. https://0-doi-org.brum.beds.ac.uk/10.1111/pbi.13127.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Li X, Cao J. Late embryogenesis abundant (LEA) gene family in maize: identification, evolution, and expression profiles. Plant Mol Biol Rep. 2016;34(1):15–28. https://0-doi-org.brum.beds.ac.uk/10.1007/s11105-015-0901-y.

    Article  CAS  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Nagaraju MSS, Kumar SA, Reddy PS, Kumar A, Rao DM, Kishor PBK. Genome-scale identification, classification, and tissue specific expression analysis of late embryogenesis abundant (LEA) genes under abiotic stress conditions in Sorghum bicolor L. PLoS One. 2019;14(1):e0209980. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pone.0209980.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Bhattacharya S, Dhar S, Banerjee A, Ray S. Structural, functional, and evolutionary analysis of late embryogenesis abundant proteins (LEA) in Triticum aestivum: a detailed molecular level biochemistry using in silico approach. Comput Biol Chem. 2019;82:9–24. https://0-doi-org.brum.beds.ac.uk/10.1016/j.compbiolchem.2019.06.005.

    Article  CAS  PubMed  Google Scholar 

  32. Zan T, Li L, Li J, Zhang L, Li X. Genome-wide identification and characterization of late embryogenesis abundant protein-encoding gene family in wheat: evolution and expression profiles during development and stress. Gene. 2020;736:144422. https://0-doi-org.brum.beds.ac.uk/10.1016/j.gene.2020.144422.

    Article  CAS  PubMed  Google Scholar 

  33. Kang QH, Jiang WD, Song XX, Sun ZY, Yuan HM, Yao YB, et al. Study Progress of Apomixis in flax (Linum usitatissimum L.). J Nat Fibers. 2019;18:1–11.

    Article  Google Scholar 

  34. Soto-Cerda BJ, Duguid S, Booker H, Rowland G, Cloutier S. Association mapping of seed quality traits using the Canadian flax (Linum usitatissimum L.) core collection. Theor Appl Genet. 2014;127:881–96.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Hall LM, Booker H, Siloto RMP, Jhala AJ, Weselake RJ. Chapter 6 - Flax (Linum usitatissimum L.). In: TA MK, Hayes DG, Hildebrand DF, Weselake RJ, editors. Industrial Oil Crops. USA: AOCS Press; 2016. p. 157–94.

    Chapter  Google Scholar 

  36. Wang Z, Hobson N, Galindo L, Zhu S, Shi D, McDill J, et al. The genome of flax (Linum usitatissimum) assembled de novo from short shotgun sequence reads. Plant J. 2012;72(3):461–73. https://0-doi-org.brum.beds.ac.uk/10.1111/j.1365-313X.2012.05093.x.

    Article  CAS  PubMed  Google Scholar 

  37. Lan T, Gao J, Zeng QY. Genome-wide analysis of the LEA (late embryogenesis abundant) protein gene family in Populus trichocarpa. Tree Genet Genom. 2013;9(1):253–64. https://0-doi-org.brum.beds.ac.uk/10.1007/s11295-012-0551-2.

    Article  Google Scholar 

  38. Zhao P, Wang D, Wang R, Kong N, Zhang C, Yang C, et al. Genome-wide analysis of the potato Hsp20 gene family: identification, genomic organization and expression profiles in response to heat stress. BMC Genomics. 2018;19(1):61. https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-018-4443-1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Song M, Zhang Y, Wang L, Peng X. Genome-wide identification and phylogenetic analysis of zinc finger Homeodomain family genes in Brassica napus. Chin Bull Botan. 2019;54:699–710.

    Google Scholar 

  40. Jeffares DC, Penkett CJ, Bahler J. Rapidly regulated genes are intron poor. Trends Genet. 2008;24(8):375–8. https://0-doi-org.brum.beds.ac.uk/10.1016/j.tig.2008.05.006.

    Article  CAS  PubMed  Google Scholar 

  41. Zafar S, Li YL, Li NN, Zhu KM, Tan XL. Recent advances in enhancement of oil content in oilseed crops. J Biotechnol. 2019;301:35–44. https://0-doi-org.brum.beds.ac.uk/10.1016/j.jbiotec.2019.05.307.

    Article  CAS  PubMed  Google Scholar 

  42. Li N, Li Y. Signaling pathways of seed size control in plants. Curr Opin Plant Biol. 2016;33:23–32. https://0-doi-org.brum.beds.ac.uk/10.1016/j.pbi.2016.05.008.

    Article  CAS  PubMed  Google Scholar 

  43. Potter SC, Aurélien L, Eddy SR, Youngmi P, Rodrigo L, Finn RD. HMMER web server: 2018 update. Nucleic Acids Res. 2018;46(W1):W200–4. https://0-doi-org.brum.beds.ac.uk/10.1093/nar/gky448.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Larkin MA, Blackshields G, Brown NP, Chenna RM, Mcgettigan PA, Mcwilliam H, et al. Clustal W. Clustal X version 2.0. Bioinformatics. 2007;23(21):2947–8. https://0-doi-org.brum.beds.ac.uk/10.1093/bioinformatics/btm404.

    Article  CAS  PubMed  Google Scholar 

  45. Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4. https://0-doi-org.brum.beds.ac.uk/10.1093/molbev/msw054.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Yu-Fu W, Yan L, Qing-Hua K, Ying LU, Xue Y, Feng-Zhi G, et al. The Breeding Report of Heiya No. 14 of New Fiber Flax Variety. China’s Fiber Crops. 2003;3(8–9):38.

    Google Scholar 

  47. Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106. https://0-doi-org.brum.beds.ac.uk/10.1186/gb-2010-11-10-r106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, et al. TM4: a free, open-source system for microarray data management and analysis. BioTechniques. 2003;34(2):374–8. https://0-doi-org.brum.beds.ac.uk/10.2144/03342mt01.

    Article  CAS  PubMed  Google Scholar 

  49. Clough SJ, Bent AF. Floral dip: a simplified method for agrobacterium-mediated transformation of Arabidopsis thaliana. Plant J. 1998;16(6):735–43. https://0-doi-org.brum.beds.ac.uk/10.1046/j.1365-313x.1998.00343.x.

    Article  CAS  PubMed  Google Scholar 

  50. Chen K, Yin Y, Liu S, Guo Z, Zhang K, Liang Y, et al. Genome-wide identification and functional analysis of oleosin genes in Brassica napus L. BMC Plant Biol. 2019;19(1):294. https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-019-1891-y.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgments

This work was conducted in the Central Laboratory of Biotechnology Research Institute, Chinese Academy of Agricultural Sciences.

Statements

All methods described above were carried out in accordance with relevant guidelines and regulations. The seeds of the flax cultivars of Heiya No. 14 and Macbeth were stored in Biotechnology Research Institute, and it is permitted to use these plant seeds in this research.

Funding

This research was financially supported in part by grants from the National Special Program of Transgenic Research (No. 2016ZX08011–001) and Science and Technology Department of Ningxia China (2021BBF02022).

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization, Y.L. and X.W.P.; methodology, Z.L.; software, Z.L. and C.Y.L.; validation, Z.L., H.C., T.B.Z. and C.Y.L.; data curation, L.D.H., H.C., C.Y.L., T.B.Z.; writing—original draft preparation, Z.L.; writing—review and editing, L.L. and Y.L.; funding acquisition, X.W.P.. All authors have read and agreed to the published version of the manuscript.

Corresponding authors

Correspondence to Xinwu Pei or Yan Long.

Ethics declarations

Ethics approval and consent to participate

Our research did not involve any human or animal subjects, material, or data. The plant materials used in this study were conserved by the Biotechnology Research Institute.

Consent for publication

Not applicable.

Competing interests

All of the authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Supplementary Figure 1 The positions of

LuLEA genes on chromosomes. Each box represents a chromosome, where the LuLEA genes are mapped with the slim bar. The genes in the same subfamily are marked by identical coloring. The scale to the left of the chromosome is in millions of bases (Mb).

Additional file 2: Supplementary Table 1

Subcellular localization prediction of all the 50 LuLEA genes.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, Z., Chi, H., Liu, C. et al. Genome-wide identification and functional characterization of LEA genes during seed development process in linseed flax (Linum usitatissimum L.). BMC Plant Biol 21, 193 (2021). https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-021-02972-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-021-02972-0

Keywords