BnLACS9 is highly homologous with AtLACS9
In Arabidopsis thaliana, nine long-chain acyl-CoA synthetases (LACSs) belong to a large superfamily of acyl-activating enzymes [6, 8], and are involved in FA transport. In order to search the homologous proteins in B. napus L., the amino acid sequences of AtLACS1-9 (AT2G47240, AT1G49430, AT1G64400, AT4G23850, AT4G11030, AT3G05970, AT5G27600, AT2G04350, AT1G77590) were used as the query probe to search in Brassica genome sequences databases (http://www.genoscope.cns.fr/blat-server/cgi-bin/colza/webBlat). In total, 29 significant homologous proteins, BnaA01g13470D, BnaA01g23830D, BnaA02g21860D, BnaA03g29320D, BnaA03g57930D, BnaA05g00630D, BnaA05g00640D, BnaA05g07920D, BnaA05g16170D, BnaA05g31340D, BnaA06g28680D, BnaA07g20920D, BnaA09g22150D, BnaA10g24190D, BnaC01g15670D, BnaC02g28920D, BnaC03g25940D, BnaC03g34500D, BnaC03g44430D, BnaC04g00320D, BnaC04g26780D, BnaC04g51420D, BnaC05g45860D, BnaC05g51350D, BnaC06g20910D, BnaC07g28060D, BnaC09g11030D, BnaC09g26090D, were gained. The phylogenetic tree analysis was performed based on the similarities of the conserved domain sequences among these proteins. As shown in Fig. 1A, except AtLACS3, other AtLACSs have 2 or 4 homologous proteins in B napus, because the B. napus is allopolyploidy, and it has a duplicate genome [52]. All LACS proteins could be divided into three branches, LACS1 to LACS5 formed the first branch, LACS6 and LACS7 formed the second branch, and LACS8 and LACS9 formed the third branch (Fig. 1A), showed that these LACS proteins might have different function. The phylogenetic tree analysis showed that AtLACS9, BnaA07g20920D (BnLACS9-A07) and BnaC06g20910D (BnLACS9-C06) were in one clade (Fig. 1A), and they share 91.6% and 91.9% similarity in amino acid sequence, respectively (Fig. 1B). AtLACS9, BnLACS9-A07 and BnLACS9-C06 are closely related in terms of amino acid sequence and evolutionary relationship, which imply they likely have similar functions. The similarity between BnLACS9-A07 and BnLACS9-C06 was as high as 99.7%, so we selected BnLACS9-C06 as the target gene in this study.
Subcellular localization of BnLACS9
Localization studies were further performed on BnLACS9 to understand its mechanism of action. To ascertain the subcellular localization of BnLACS9, GFP was fused to the C terminus of BnLACS9 under the control of 35S promoter, and the fusion gene was transformed into N. benthamiana leaves and the cotyledons of B. napus, respectively. We found that the expression pattern of BnLACS9-GFP overlapped with the chlorophyll autofluorescence in N. benthamiana leaves and B. napus cotyledons (Fig. 3A-F), suggesting that BnLACS9 may be targeted to the chloroplasts. In order to further study the location of BnLACS9-GFP in chloroplasts, we expressed 35S:BnLACS9-GFP in N. benthamiana leaves for subsequent protoplast isolation and observed the N. benthamiana protoplasts by three-dimensional modeling. From the Fig. 3G-I, we observed circular structures around the chloroplast, and formed loop structures and thin tubules form the chloroplast. All these observations we could conclude that BnLACS9 may located on the chloroplast envelope.
BnLACS9 can complement LACS function in deficient yeast YB525
To determine whether BnLACS9 has LACS activity, a yeast vector pYES2-BnLACS9 was constructed and analyzed by yeast complementary expression system. Yeast strain YB525 is a strain with affected LACS activity [53], and this defective yeast can’t grow in the medium only using FA as the sole carbon source. Yeast complementarity test was used to determine the growth of yeast cells transferred into different vectors in drop-out medium with different FAs as the sole carbon source for 72 hours. The YB525 yeast cells transferred into pYES2 empty vector could not grow in drop-out medium with different FAs as the sole carbon source. While the YB525 yeast cells transferred into pYES2-BnLACS9 vector grow better and has different growth rates in drop-out medium with different FAs as the sole carbon source (Fig. 4). The OD value of YB525 yeast cells transferred into pYES2-BnLACS9 vector could reach to 0.38 in drop-out medium with only fatty acid C18:0 as the sole carbon source, while the OD value of yeast 0.05 in drop-out medium with only fatty acid C12:0 as the sole carbon source (Fig. 4). These results indicate that the heterologous expression of BnLACS9 gene can complement the deficient yeast and has the activity of LACS protease.
Transient expression of BnLACS9 in N. benthamiana leaves could enhance the content of FA, MGDG and total chlorophyll
In order to further study the function of BnLACS9 gene, BnLACS9 was transient expressed in N. benthamiana leaves. In plants, FAs are uniquely synthesized in plastids [54]. FA compositions from N. benthamiana leaves of both accessions were determined by gas chromatography-mass spectrometry (GC-MS). From Table 1, we could find that the transient expression of BnLACS9 gene in N. benthamiana leaves can increase the content of FAs, especially the C16:0 and C18:2, their contents were 221.3034 ± 14.0447 avfmol/mg and 59.2165 ± 9.703042 avfmol/mg, respectively. While the C16:0 and C18:2 are 134.9407 ± 7.264136 avfmol/mg and 23.56197 ± 3.008565 avfmol/mg in the control, respectively (Table 1). As a result, the total FA content in the transient expression of BnLACS9 gene in N. benthamiana leaves was 574.5768 ± 35.79813 avfmol/mg, while that of control was only 395.4662 ± 13.63986 avfmol/mg. These results indicated that BnLACS9 could increase the FA content in N. benthamiana leaves.
Table 1
The CoA content in N. benthamiana leaves
Acyl-CoA | P19 CoA content (avfmol/mg) | ±SD | pB2GW7.0-BnLACS9 CoA content (avfmol/mg) | ±SD |
C2:0 | 6.61083 | 1.39379 | 7.1831 | 1.358259 |
C4:0 | 7.387154 | 6.789183 | 14.32233 | 0.97188 |
C6:0 | 8.322458 | 2.74276 | 6.547055 | 1.406957 |
C8:0 | 36.06325 | 6.266329 | 42.23341 | 5.334781 |
C10:0 | 0.647754 | 0.567659 | 0.869375 | 0.580957 |
C12:0 | 1.564231 | 0.299744 | 0.770367 | 0.516654 |
C14:0 | 10.88591 | 0.777906 | 12.78282 | 0.763995 |
C18:3 | 47.59071 | 6.161552 | 41.7954 | 6.390999 |
C16:1 | 4.173502 | 1.654766 | 7.576914 | 0.711642 |
C18:2 | 23.56197 | 3.008565 | 59.2165 | 9.703042 |
C16:0 | 134.9407 | 7.264136 | 221.6034 | 14.0447 |
C18:1 | 6.365882 | 0.315986 | 7.941865 | 1.109978 |
C20:2 | 0.967406 | 0.186929 | 1.215209 | 0.579532 |
C18:0 | 30.6479 | 0.65538 | 36.49893 | 1.843828 |
C20:1 | 24.87717 | 1.749918 | 45.81828 | 3.616426 |
C20:0 | 12.76891 | 0.180914 | 16.41325 | 0.587022 |
C22:1 | 26.27505 | 0.531116 | 38.57146 | 1.448839 |
C22:0 | 11.81541 | 0.246871 | 13.21715 | 0.53581 |
TOTAL | 395.4662 | 13.63986 | 574.5768 | 35.79813 |
We also detected the chlorophyll content in N. benthamian leaves of transient expression of BnLACS9 and its control. We found that within 5 days after BnLACS9 gene was transferred into N. benthamiana leaves, there was no significant difference in chlorophyll content between the two groups. From the sixth day, the chlorophyll content of N. benthamiana leaves transfected with BnLACS9 gene was significantly higher than that of wild type, and reached the maximum on eighth day (Fig. 5A). MGDG is one of major component of chloroplast photosynthetic membrane [54]. We also detected the monogalactosyl-diaclyglycerol (MGDG) content in N. benthamian leaves between transient expression of BnLACS9 and its control. The content of MGDG in N. benthamiana leaves transfected with BnLACS9 gene was significantly higher than that of the control (Fig. 5B). All these results indicated that BnLACS9 could enhance the content of FA, MGDG and total chlorophyll.
Overexpression of BnLACS9 in B. napus could enhance the content of total chlorophyll
We found that the transient overexpression of BnLACS9 in N. benthamiana enhanced the content of FAs, MGDG and total chlorophyll (Fig. 5 and Table 1). To further address the function of BnLACS9 in rapeseed plant, transgenic rapeseed plants overexpressing BnLACS9, as well as those suppressing BnLACS9 (Bnlacs9), were generated (Fig. S1). Homozygous BnLACS9-overexpressed lines (#6, #12 and #18) were selected from the positive transgenic plants to evaluate the function of BnLACS9.
It was found that the cotyledons and leaves of BnLACS9-overexpressed plants were significantly greener than those of the wild type (Fig. 6A and B). Chlorophyll content of BnLACS9-overexpressed plants was also higher than that of wild type in 7-day-old cotyledons and 30-day-old leaves (Fig. 6C and D). We also found that the siliques of BnLACS9-overexpressed plants could delay senescence compared with control (Fig. 6E and F). So, we measured the chlorophyll content of the silique at 30 DAF and 50 DAF, and showed that the chlorophyll content in overexpressed plants was significantly higher than that in wild type plants (Fig. 6G and H). These results indicated that overexpression of BnLACS9 gene can increase chlorophyll content and delay senescence in B. napus.
Overexpression of BnLACS9 in B. napus could increase the number of thylakoid layer structure in chloroplast
The BnLACS9-overexpressed plants had higher chlorophyll content (Fig. 6), so we observed the chloroplast structure by transmission electron microscopy (TEM). Ultrastructural observation showed that the number of thylakoid layer structure in BnLACS9-overexpressed plants was more than that in wild type plants (Fig. 7). We also measured the photosynthetic efficiency of leaves of overexpressed plants and wild-type plants. We found that the photosynthetic efficiency of leaves of BnLACS9-overexpressed plants was significantly higher than that of wild-type plants (Fig. S2). These results indicated that overexpression of BnLACS9 gene could increase the number of chloroplast thylakoid layer structure, make leaves greener and improve photosynthetic efficiency.
In order to characterize the mechanism that BnLACS9 influenced the content of chlorophyll through regulating the genes in the pathway of the chlorophyll synthesis, we sequenced the transcriptome of the leaves of the overexpression of BnLACS9 plants and wide type. Transcriptome analysis showed 9 genes in chloroplast synthesis pathway (comp613754_c0, comp879850_c0, CHLD comp44567_c0, comp956830_c0, comp55669_c1, comp58218_c0, POR comp57696_c0, comp38830_c0, comp57073_c0) were up-regulated in BnLACS9-overexpressed plants (Fig. 8). We confirmed the transcript levels of genes HEMA, CHLD, PORB, CAO by RT-PCR (Fig. 8B), and showed that these genes were all up-regulated. The result suggested the BnLACS9 upgraded the key genes of chlorophyll synthesis pathway and further increased the content of the chlorophyll.
Overexpression of BnLACS9 in B. napus increased the content of the galactolipids
Our studies showed that the BnLACS9-overexpressed plants had more thylakoid layer structure than wild type plants (Fig. 7). It’s known that the main lipid component of the thylakoid layer structure is MGDG, DGDG, PG and SQDG. MGDG and DGDG account for approximately 84% of thylakoid membrane lipids [24]. PG is the main chloroplast phospholipid, represents about 8% of the chloroplast lipids. Another glycolipid, sulfoquinovosyldiacyl glycerol (SQDG), is present from 3–10% of the chloroplast membrane lipids [55]. Transient expression of BnLACS9 in N. benthamiana increased the synthesis of MGDG and chlorophyll. Overexpression of BnLACS9 also increased the content of the chlorophyll in B. napus (Fig. 9). In addition, the number of the thylakoid of the chloroplast was also increased.
Previous studies have shown that overexpression of BnLACS9 could improve the chlorophyll content and photosynthetic efficiency of leaves in rapeseed (Fig. 6 and S2). And overexpression of BnLACS9 in N. benthamiana leaves could increase TAG content (Table 1). We also found the transcript sequencing results of over-expressed BnLACS9 plants showed a significant increase in the expression of TAG synthesis pathway-related genes (Fig. 8). So, we measured the dry weight of the BnLACS9-overexpressed lines at 20 and 40 days after germination, and detected the seeds oil content of the transgenic plants by Nuclear Magnetic Resonance (NMR) and Near Infrared (NIR). At 20 days, we could see that the dry weights have been increased. At 40th day after germination, the dry weight of BnLACS9-overexpressed lines was particularly increased (Fig. 3S). As shown in Table 2, the oil content of BnLACS9-12 lines increased most obviously, reaching 45%, BnLACS9-6, BnLACS9-13 and BnLACS9-18 lines increased about 3%. In addition, the seed oil content of overexpressed BnLACS9 plants was determined by NIR, and the results also showed the same with NMR (Table 2). From the above, overexpression of BnLACS9 in plants could increase the biomass of the plant and increase the seed oil content in rapeseed.
Table 2
The content of seed oil of BnLACS9 overexpressed lines is measured by NMR and NIR
line | measured by NMR (%) | measured by NIR (%) |
NY12 | 39.70 | 38.99 |
BnLACS9-6 | 43.34 | 41.93 |
BnLACS9-12 | 45.64 | 44.97 |
BnLACS9-13 | 41.18 | 41.81 |
BnLACS9-17 | 41.42 | 39.25 |
BnLACS9-18 | 43.86 | 41.59 |
NMR: Nuclear magnetic resonance |
NIR: Near-infrared |
The BnLACS9 affected the biogenesis of the chlorophyll by regulating the chloroplast membrane lipid biosynthesis. In order to make a thorough inquiry about how the content of the glycolipid was increased in the overexpression of BnLACS9 plants and to elucidate the gene regulatory network of the pathway. Transcriptomes in the overexpression of BnLACS9 line and wild type were analyzed. The unigenes related to the glycolipid (Fig. 10A) synthesis were picked out from the transcriptomes data of the B. napus. The unigene was shown as ‘comp…_c0’. The expression level of the unigenes was represented by RPKM [56]. Nine unigenes for glycolipid biosynthesis were chosen for RT-PCR to confirm the expression level of the unigenes with RPKM (Fig. 10B). These unigenes were ATS1 (comp_775367c0), LPAAT (comp50451_c0), PP (comp30817_c0), PGPS (comp41948_c0), PGPP (comp56019_c0), MGDGS (comp55779_c0), DGDGS (comp56744_c0), SLS (comp52127_c0), CDP-DAGS (comp42015_c0), which the RPKM is the highest (Fig. 10). The diacylglycerol backbones LPA is assembled by the ATS1 [57], LPA in turn is then acylated by the lysophosphatidic acid acyltransferase (LPAAT) to phosphatidic acid (PA), PA is dephosphorylated to DAG by an inner-envelope localized enzyme PP [25]. CDP-DAGS, PGPS and PGPP are the genes involved in the PG synthesis [58, 59]. MGDGS and DGDGS are the enzymes responsible for the synthesis of MGDG and DGDG [29, 60]. The SLS takes charge of synthesis of SQDG [41, 61]. The result showed that ATS1, LPAAT, PP, PGPS, PGPP, PP, MGDGS, DGDGS, SLS, CDP-DAGS had higher expression in the overexpression of BnLACS9 plant than the wild type. In particular, the expression level of the ATS1 and SLS is four times higher more than the CK. Other genes in overexpression of BnLACS9 plant is exceeded double than the CK. Taken together, the results suggest that BnLACS9 regulated the expression of genes in the pathway of glycolipid biosynthesis. The higher expression level of these genes increased the content of the glycolipid.