NtAn1 genes were highly expressed in developing seed
Previous study reported that Arabidopsis TT8 gene has two homologs in tobacco genome, NtAn1a and NtAn1b [23]. Tobacco is a natural allotetraploid plant, sequence analysis revealed that NtAn1a originated from N. sylvestris, whereas NtAn1b derived from N. tomentosiformis. qPCR analysis showed that NtAn1a and NtAn1b expressed at developing flowers with highest expression level in corolla limb, which were consistent with the function of flower flavonoid biosynthesis regulation [23]. However, we noticed that the transcript level of both NtAn1a and NtAn1b were relatively high in developing ovary, which indicated that NtAn1 genes might play an important role in seed development.
To further confirm and characterize the expression pattern of NtAn1a and NtAn1b, their expression were assessed in different organs and various stages of seed development: 7, 14, 21, and 28 days after flowering (DAF). qRT-PCR results showed that NtAn1a and NtAn1b had a similar expression pattern, with the highest expression in developing seeds at 7 DAF and decreased at later stages (Fig. 1). The expression of both NtAn1 genes exhibited a high expression level in flower, which were consistent with the previous results. However, low transcript levels were detected in root, leaf, and stem (Fig. 1). These results suggested that NtAn1a and NtAn1b might regulate PAs and lipid accumulation during seed development in a way like in Arabidopsis and rapeseed.
Targeted mutagenesis of NtAn1 using the CRISPR-Cas9 system
NtAn1a and NtAn1b showed a high sequence identity of 92.95% and 90.36% at the nucleotide and protein levels, respectively (Additional file 1: Figure. S1 and Figure. S2). Furthermore, these two genes had similar expression patterns (Fig. 1). These results suggested that these two genes may have similar and redundant functions. Thus, two gRNAs recognizing both NtAn1 genes were designed to effectively knockout of them, and both of the gRNAs targeting the second extron of the coding sequence (Fig. 2a). The CRISPR-Cas9 construct containing these two gRNAs, which driven by the Arabidopsis U6-26 and U6-29 promoter, respectively (Fig. 2b), was produced based on the CRISPR-Cas9 multiplex genome-editing vector [24]. The resulting construct was transformed into wild type (WT) tobacco plant using Agrobacterium-mediated leaf disc transformation method. Through kanamycin selection, 12 kanamycin resistance T0 transgenic plants were generated. The targeting region of both NtAn1a and NtAn1b were amplified by a pair of primers at the same time. Two homozygous mutant lines (an1-1 and an1-2) were identified from the T0 transgenic plants by Sanger sequencing analysis of the gRNAs targeting region. The an1-1 mutant had one base insertion at both gRNA targeting sites of the NtAn1a and NtAn1b genes, while the an1-2 mutant line had a 105 base fragment deletion between the two gRNAs targeting sites (Fig. 2c and Fig. 2d). T-DNA free mutant plants were selected from the T1 progeny generated by the self-pollinated of the two independent T0 homozygous mutant lines. Twenty T-DNA free T1 generation plants from each mutant line were randomly selected for further analysis.
Mutation of NtAn1 genes resulted in yellow seed coat
The formation of seed color in most plant species is due to the deposition of PAs within the endothelial layer of the inner integument of the seed coat [13]. Previous studies demonstrated that TT8 played a key role in regulating PAs accumulation in various plants [18, 19]. In this paper, mutation of the NtAn1, the TT8 homolog genes in tobacco, generated yellow-seeded phenotype (Fig. 3a). This indicated that targeted mutation of the NtAn1 genes might disrupt the accumulation of PAs in tobacco seed coat. In order to check the PAs deposition visibly, the mutant tobacco seeds were dyed by DMACA reagent. The results further confirmed the defects of PAs accumulation in seed coat (Fig. 3b). The soluble and insoluble PAs contents were calculated quantitatively by spectrophotometric method. The results showed that the PAs were mainly stored in the insoluble form in tobacco seed coat, and targeted mutation of NtAn1 genes led to the significant decreases in both soluble and insoluble PAs content compared with those in WT tobacco seed coat (Fig. 3c).
ANR (Anthocyanidin reductase) and LAR (Leucoanthocyanidin reductase) are two key enzymes participated in PAs biosynthesis. ANR converts anthocyanidins to the epicatechin [25], whereas LAR could reduce leucocyanidin to catechin [26]. Catechin and epicatechin are considered to be the main building blocks for PAs biosynthesis. The encoding genes of ANR and LAR were directly regulated by TT8 gene in Arabidopsis [15]. The expression patterns of the homolog genes encoding these two enzymes during seed development were analyzed by qRT-PCR. Our results showed that both of these two genes had a similar expression profile as NtAn1, and the expression of them were significantly decreased in both an1-1 and an1-2 mutant seed (Fig. 3d). Taken together, these findings indicate that tobacco NtAn1 regulated the accumulation of PAs in a similar way like in Arabidopsis, and mutation of the tobacco NtAn1 genes could hinder the PAs deposition in the seed coat, which was consistent with the phenotypes observed in tt8 mutant seed in Arabidopsis and other Brassica species [12, 17, 21, 22].
Targeted mutagenesis of NtAn1 generated white flower
Previous study had demonstrated that the anthocyanin accumulation in transgenic tobacco flowers could be significantly elevated by overexpression of NtAn1a or NtAn1b gene [23]. The early biosynthesis genes (EBGs) and late biosynthesis genes (LBGs) in anthocyanin pathway, including CHS, CHI, F3H, DFR, and ANS, were dramatically induced by the overexpression of NtAn1 genes [23]. In this present paper, we found that targeted mutagenesis of NtAn1 genes generated white flower phenotype, which resulted from the defects in anthocyanin accumulation in the flower (Fig. 4a). TT8 genes regulate the biosynthesis of anthocyanin by manipulating the expression of the LBGs in anthocyanin pathway. The expression level of the downstream genes at different flower development stages were analyzed by qRT-PCR. Our results revealed that the examined anthocyanin biosynthesis genes expressed at all three developmental stages with expression level peaking at the late stage in WT plant flower (Fig. 4b). The expression patterns of these genes were consistent well with those of NtAn1 genes, which indicating the regulation relation between them [23]. By contrast, the expression of all examined anthocyanin biosynthesis genes at different developing stages were significantly repressed in an1-1 mutant line (Fig. 4b). Taken together, our results demonstrated again that NtAn1 genes played an essential role in the biosynthesis of anthocyanin in tobacco flower.
Targeted mutagenesis of NtAn1 increased seed lipid and protein contents
Both natural and targeted mutation of TT8 genes in Arabidopsis and Brassica species would result in a significant increases in seed lipid content. To characterize the effect of the NtAn1 genes targeted mutation on tobacco seed lipid accumulation, the seed lipid content was analyzed by GC-MS method. The results showed that WT tobacco seed lipid content was about 38.77 µg per seed, and the lipid content was approximately 45.91 µg per seed in an1-1 and 44.97 µg per seed in an1-2 mutant line, increased significantly by 18.42% and 15.99% relative to the WT seeds, respectively (Fig. 5a). These results indicated that the TT8 gene and its homologs regulated seed lipid accumulation in a conserved way among different plant species. Thus, the TT8 homologs from other oilseed plant could be used as a target to enhance the seed lipid content.
In most oilseed crops, the seed lipid content is negatively correlated with the protein content. Surprisingly, The BnTT8 mutant seeds showed simultaneous increases in both lipid and protein contents [12], which is different from the lower protein content in Arabidopsis tt8 mutant seed [21]. These results indicated that the seed protein accumulation was regulated by different mechanism in Arabidopsis and Brassica napus. To determine the effect of NtAn1 genes on tobacco seed protein accumulation, the seed protein contents in WT and an1 mutant lines were examined using the Pierce BCA Protein Assay Kit. The results showed that WT tobacco seed protein content was about 32.79 µg per seed, and the protein content was elevated to 36.56 µg per seed in an1-1 and to 35.97 µg per seed in an1-2 mutant line, increased significantly by 11.50% and 9.70% relative to the WT seeds, respectively (Fig. 5b). These results were consistent with those in Brassica napus. The increased protein content of the an1 mutant seeds could make the tobacco seed meal produced after lipid extraction to be more valuable for animal feed manufacture.
The property of biodiesel is partially determined by the fatty acid chain length, the position and number of the double bonds [27]. In Arabidopsis and Brassica napus, the mutation of TT8 gene resulted in an alteration in the seed fatty acid profile, including increases in palmitic acid, linoleic acid, and linolenic acid, while decreases in stearic acid and oleic acid compared with the WT seeds [12, 21]. Fatty acid composition of tobacco seed lipid shows a main presence of palmitic acid, stearic acid, oleic acid, and linoleic acid [3]. Possibly due the differences in the fatty acid composition, targeted mutation of NtAn1 genes just resulted in a significant decrease in stearic acid, and the other four main fatty acid components were not changed significantly compared with the WT tobacco seed (Table. 1). WT tobacco seed lipid shows a high ratio of linoleic acid (~ 72%), which would make biodiesel produced from tobacco seed lipid more susceptible to oxidation, and this will limit its use in traditional engine. High oleic acid content is a preferred trait for biodiesel feedstock production [27]. Our previous work had generated a high oleic acid tobacco variety through CRISPR-Cas9 mediated knockout of the NtFAD2- genes [28]. Next, a high lipid content tobacco seed with an ideal fatty acid profile could be expected by combining the high lipid content trait generated in this present study with the high oleic acid phenotype through hybridization or multiple gene editing.
Table 1
Fatty acid composition of the an1 mutant lines
| Fatty acid composition (mol%) |
16:0 | 18:0 | 18:1Δ9 | 18:2Δ9, 12 | 18:3Δ9, 12, 15 |
WT | 11.14 ± 0.37 | 4.36 ± 0.36 | 12.86 ± 0.35 | 71.14 ± 1.29 | 0.70 ± 0.02 |
an1-1 | 11.37 ± 0.28 | 2.77 ± 0.24** | 13.02 ± 0.38 | 72.69 ± 1.43 | 0.98 ± 0.03 |
an1-2 | 11.96 ± 0.39 | 2.35 ± 0.21** | 12.77 ± 0.29 | 73.08 ± 1.39 | 0.73 ± 0.01 |
16:0, palmitic acid; 18:0, stearic acid; 18:1Δ9, oleic acid; 18:2Δ9, 12, linoleic acid; 18:2Δ9, 12, 15, linolenic acid. The data are the mean ± SD (n = 3). |
Expression of genes involved in seed development and lipid biosynthesis were altered by NtAn1 targeted mutation
In Arabidopsis, TT8 protein could repress the lipid biosynthesis pathway through directly binding to the promoter region of the critical transcriptional factors important for seed development, such as LEAFY COTYLEDON1 (LEC1), LEC2, and FUSCA3 (FUS3). Thus, the Arabidopsis tt8 mutant showed increased seed lipid content [21]. In Brassica napus L., The expression levels of several genes involved in fatty acid biosynthesis during seed development were increased in BnTT8 mutant plants generated by CRISPR-Cas9 system [12]. In this present paper, targeted mutagenesis of NtAn1 genes led to the increases of both lipid and protein content, we suggested that NtAn1 genes might regulate seed development and storage component accumulation in a way similar to that in Brassica napus. To test this hypothesis, qRT-PCR was performed in the process of seed development to examine the expression of genes involved in seed development regulation and fatty acid biosynthesis, including LEC1, LEC2, FUS3, KASI, PI-PKβ1, and BCCP2. Our results showed that the expression of all examined genes were up-regulated in the an1-1 mutant at one or two developmental stages (Fig. 6). For example, the LEC1 and FUS3 genes were significantly increased at 14 DAF and 21 DAF seeds, while the KASI and PI-PKβ1 genes were up-regulated at 14 DAF stage (Fig. 6). The elevated expression levels of genes involved in seed development and lipid biosynthesis could explain the enhanced lipid content in mutant lines.
Seed yield related traits were not affected by the mutation of NtAn1
Besides the seed lipid content, seed yield is another important factor affecting the lipid yield. The seed yield related traits of the an1-1 and an1-2 mutant lines were also evaluated. The results showed that the yield related traits, including seed size (Fig. 7a and Fig. 7b), seed weight (Fig. 7c), fruit number per plant (Fig. 7d), and seed number per fruit (Fig. 7e) were not affected by targeted mutation of NtAn1 genes. Thus, the targeted mutation of NtAn1 genes could generate a useful tobacco variety with a high seed lipid yield and improved nutritional quality.