Identification of new copies of TaABI5 in wheat line 8425B
The full-length sequence of TaABI5 was isolated from Zhou 8425B using a universal primer set (Table S1). A total of 16 copies of TaABI5 from Zhou 8425B were detected based on the comparison of ABI5 with the sequences deposited in the IWGSC and the 2005 Supplement of the Wheat Gene Catalogue. These copies were named TaABI5-A1, TaABI5-A2, TaABI5-A3 and TaABI5-A4 on the short arm of chromosome 3A, TaABI5-B1, TaABI5-B2, TaABI5-B3, TaABI5-B4 and TaABI5-B5 on the short arm of chromosome 3B, TaABI5-D1, TaABI5-D2, TaABI5-D3, TaABI5-D4, TaABI5-D5, TaABI5-D6 and TaABI5-D7 on the short arm of chromosome 3D, respectively. The results showed that TaABI5 was present in Triticum aestivum as a multi-copy gene. Sequence analysis indicated that the 16 copies of TaABI5 had 84.43–94.06% similarity with sequence AB238934.1 at the nucleotide level; the length of these copies ranged from 876 to 1556 bp. All sequences were deposited in GenBank under accession numbers MK287847 to MK287860 and MK334234 to MK334236 (https://www.ncbi.nlm.nih.gov/gene/).
Selection of specific primers for identifying the allelic variation of TaABI5 and sequences analysis
Based on the 16 copies of TaABI5, gene specific primers were designed to amplify all copies of TaABI5 sequences. In total, three specific primer pairs, TaABI5D3F/R, TaABI5D6F/R, and TaABI5A4F/R, were identified to specifically amplify the sequences of TaABI5-D3, TaABI5-D6, and TaABI5-A4, respectively. However, no specific primers could be designed for the other 13 copies of TaABI5 because of very high sequence similarities among them.
Ten cultivars (Xiaobaiyuhua, Yumai 18, Yangxiaomai, Xiaoyuhua, Xiaoyan 6, Zhou 8425B, Jimai 19, Jing 411, Hengshui 7228, and Zhongyou 9507) with different levels of seed dormancy were selected to identify allelic variation of TaABI5-D3, TaABI5-D6, and TaABI5-A4. In these ten cultivars, four variants were detected in TaABI5-D3, designated TaABI5-D3a, TaABI5-D3b, TaABI5-D3c, and TaABI5-D3d, respectively. Compared with the sequence of TaABI5-D3a from Zhou 8425B, both TaABI5-D3c and TaABI5-D3d had two SNPs (A to G at position 449 bp changed His into Arg, and A to G at position 705 bp changed Asn into Asp), and TaABI5-D3b had two SNPs (A to G at position 449 bp changed His into Arg, and T to C at position 1270 bp in the intron). In TaABI5-D6, five allelic variants were detected, viz. TaABI5-D6a, TaABI5-D6b, TaABI5-D6c, TaABI5-D6d, and TaABI5-D6e. Compared with the sequence of TaABI5-D6a from Zhou 8425B, InDels and SNPs located in the first exon were identified in TaABI5-D6b, TaABI5-D6c, TaABI5-D6d, and TaABI5-D6e.
There were two allelic variants in TaABI5-A4, i.e. TaABI5-A4a and TaABI5-A4b. Compared with the sequence of TaABI5-A4a, TaABI5-A4b had three SNPs (A to G at positions 269 bp and 473 bp, and C to T at position 686 bp) (Fig. S1). TaABI5-A4a was detected in five PHS-resistant cultivars (Xiaobaiyuhua, Yumai 18, Yangxiaomai, Xiaoyuhua, and Xiaoyan 6), while TaABI5-A4b was identified in five PHS-sensitive cultivars (Zhou 8425B, Jimai 19, Jing 411, Hengshui 7228, and Zhongyou 9507) (Fig. 1).
In addition, the sequence analysis showed that TaABI5-A4a and TaABI5-A4b could not encode full-length proteins, compared with full length of TaABI5 protein, as there is a stop codon at the position of 266th amino acid.
Expression of TaABI5s and TaABI5-A4 from Wanxianbaimaizi and Jing 411 at different seed developmental stages
To evaluate the potential influence of TaABI5 and the specific copy of TaABI5-A4 in Jing 411 (GI = 64.2) and Wanxianbaimaizi (GI = 7.6), the expression patterns of Wanxianbaimaizi (TaABI5-A4a) and Jing 411 (TaABI5-A4b) were determined using semi-quantitative RT-PCR analysis. The transcript expression levels of TaABI5 showed an increasing trend in seeds from 10 to 30 days post-anthesis (DPA) with the highest level at 30 DPA and the lowest at 10 DPA (Fig. 2a). In addition, the transcript expression levels of TaABI5-A4 and TaABI5 had the same trend in Wanxianbaimaizi (TaABI5-A4a) and Jing 411 (TaABI5-A4b) at 10, 20 and 30 DPA (Figs. 2a and 2b), although Wanxianbaimaizi (TaABI5-A4a) had higher transcript levels at all stages than Jing 411 (Figs. 2a and 2b).
Three PHS-resistant (Xiaoyan 6, Wanxianbaimai and Xiaobaiyuhua, TaABI5-A4a) and three PHS-susceptible genotypes (Jing 411, Zhou 8425B and Hengshui 7228, TaABI5-A4b) were chosen to analyze the transcript level by qPCR at mature seeds. Significantly higher transcript levels of TaABI5 and TaABI5-A4 were observed in the genotype TaABI5-A4a than in TaABI5-A4b (Figs. 3a and 3b).
Development and validation of gene-specific markers TaABI5A4a and TaABI5A4b for PHS resistance
Based on the SNPs between the sequences of TaABI5-A4a and TaABI5-A4b, two complementary STS markers were developed, designated TaABI5A4a (primer set TaABI5A4aF/R) and TaABI5A4b (TaABI5A4bF/R), respectively, and used for association analysis with 103 wheat cultivars and advanced lines. It was first to amplify the cultivars with the primer set TaABI5A4F/R, and the PCR product was used as the template in the following PCR system. The primer set TaABI5A4aF/R produced a 494-bp fragment in the TaABI5-A4a genotype and no PCR product in TaABI5-A4b (Table 2 and Fig. 1a), while the primer set TaABI5A4bF/R generated a 703-bp fragment in the TaABI5-A4b genotype (Table 2 and Fig. 1b) and no PCR product in TaABI5-A4a. Among the 103 cultivars and lines tested, 51 had TaABI5-A4a presenting a 494-bp fragment and 52 contained TaABI5-A4b with a 703-bp fragment (Table 2 and Table S1); TaABI5-A4a and TaABI5-A4b accounted for 49.5% and 50.5% with average GI values of 22.2 and 49.6, respectively. The GI values of the 103 lines were consistent across two years (r = 0.966, P < 0.0001), with mean GI values and standard deviations of 36.1 ± 20 in 2006 and 35.9 ± 19 in 2007, respectively. The results indicated that the genotype TaABI5-A4a were more resistant to PHS than TaABI5-A4b (Table 2).
Table 2
Association analysis between TaABI5-A4 and GI values in 103 wheat cultivars and advanced lines with different PHS tolerance using the GIM model
Genotype | No. of lines | Average GI (%) values | F value | Phenotypic variance explained (%)R2 |
TaABI5-A4a | 51 | 22.2 | 114.054* | 52.6 |
TaABI5-A4b | 52 | 49.6 |
*Significant association between TaABI5-A4 and phenotypic values at P < 0.001 level |
To further confirm the association, 200 RILs from the cross of Yangxiaomai/Zhongyou 9507 were genotyped using TaABI5A4a and TaABI5A4b specific markers (Fig. 2). Statistical analysis confirmed the significant association (P < 0.001) of allelic variation of TaABI5-A4a/TaABI5-A4b with GI values and PHS tolerance. In this population, TaABI5-A4 explained 42.3%, 60.8%, and 63.7% of the phenotypic variations in Shijiazhuang, Beijing and the averaged GI values of two environments, based on the test of two complimentary STS markers TaABI5A4a and TaABI5A4b (Table 3).
Table 3
Association analysis between TaABI5-A4 and GI values in 200 RILs using the GIM model
Trait | TaABI5-A4 | No. of lines | Average GI% values | F value | Phenotypic variance explained (%) R2 |
Beijing GI value (%) | TaABI5-A4a | 130 | 20.28 | 310.001* | 60.8 |
TaABI5-A4b | 70 | 56.54 |
Shijiazhuang GI value (%) | TaABI5-A4a | 130 | 6.40 | 147.101* | 42.3 |
TaABI5-A4b | 70 | 29.08 |
Average GI value (%) | TaABI5-A4a | 130 | 13.34 | 349.687* | 63.7 |
TaABI5-A4b | 70 | 42.81 |
*Significant association between TaABI5-A4 and phenotypic values at P < 0.001 level |
Phenotypic evaluation of transgenic rice lines
Fifteen transgenic rice plants (T1) for each of TaABI5-A4a-GFP and TaABI5-A4b-GFP constructs were selected using hygromycin. The primers of hpt557F/R were used to identify transgenic events. The frequencies of positive TaABI5-A4a-GFP and TaABI5-A4b-GFP transgenic rice plants were 73% and 80%, respectively (Fig. S2).
Plant height, germination index, 100-grain weight, stem length and diameter of T2 lines were phenotyped and analyzed (Fig. 4 and Table 4). The plant height of TaABI5-A4b-GFP transgenic lines (84.76 cm) was significantly higher (P < 0.01) than that of the negative control with an empty vector (80.33 cm) and TaABI5-A4a-GFP transgenic lines (79.58 cm). The average GI values of TaABI5-A4a-GFP, TaABI5-A4b-GFP transgenic lines and control plants were 57.7, 67.5, and 62.1, respectively, with significant differences (P < 0.05). The lengths of the first internodes of TaABI5-A4a-GFP, TaABI5-A4b-GFP transgenic lines and control plants were 1.38 cm, 1.57 cm and 1.31 cm, and those of the second internodes were 3.09 cm, 3.65 cm and 3.32 cm, respectively. The data analysis showed that the lengths of the first and second internodes between TaABI5-A4a-GFP and TaABI5-A4b-GFP transgenic lines also were significantly different (P < 0.001). The length of the third internode of TaABI5-A4b-GFP transgenic line was 7.98 cm, with significant differences from the control (7.57 cm) and TaABI5-A4a-GFP (7.26 cm) (P < 0.05 and P < 0.01, respectively) (Fig. 4 and Table 4). The diameters of the second and third internodes of TaABI5-A4b-GFP and TaABI5-A4b-GFP transgenic lines were 2.91 mm and 2.76 mm, 2.83 mm and 2.76 mm, respectively, all with significant differences (P < 0.05 and P < 0.01, respectively). Analysis of variance indicated that the TaABI5-A4b-GFP genotype had significantly (P < 0.01) higher average GI values than the TaABI5-A4a-GFP (Fig. 4). These showed that the TaABI5-A4 gene affected not only GI values in transgenic rice lines, but also plant height, internode length and diameters (Fig. 4, Table 4).
Table 4
Phenotypic data of TaABI5-A4a-GFP and TaABI5-A4b-GFP transgenic rice lines
| Control | TaABI5-A4a-GFP | TaABI5-A4b-GFP |
Plant height (cm) | 80.33 | 79.58 | 84.76** |
GI% value | 62.13 | 57.74* | 67.54* |
100-grain weight (g) | 2.18 | 2.17 | 2.07 |
Length of first internode (cm) | 1.31 | 1.38** | 1.51**** |
Length of second internode (cm) | 3.32 | 3.09* | 3.65*** |
Length of third internode (cm) | 7.57 | 7.26 | 7.98* |
Diameter of first internode (mm) | 3.09 | 3.14 | 3.09 |
Diameter of second internode (mm) | 3.09 | 3.01 | 2.90* |
Diameter of third internode (mm) | 2.88 | 2.83 | 2.76** |
*, **, *** and **** represent significant differences between transgenic plants and controls at P < 0.05, P < 0.01, P < 0.001, P < 0.0001, respectively. Control stands for transgenic plants with an empty vector as a negative control |
Endogenous ABA and GA contents were examined in mature seeds of the control, TaABI5-A4a-GFP and TaABI5-A4b-GFP transgenic lines by HPLC. The endogenous ABA in seeds of TaABI5-A4a-GFP and TaABI5-A4b-GFP were 0.71 and 0.41 µg/g, respectively, and the endogenous GA contents were 9.82 and 12.14 µg/g, respectively, (Figs. 5a and 5b). These results indicated that the TaABI5-A4a-GFP transgenic rice line had significantly (P < 0.0001) higher endogenous ABA and significantly (P < 0.01) lower endogenous GA contents in mature seeds than TaABI5-A4b-GFP.
ABA sensitivity in transgenic rice lines
In order to determine the ABA responsiveness in the seeds of TaABI5-A4a-GFP and TaABI5-A4b-GFP lines, GI values of the seeds were evaluated in water and in 50 µM ABA solution, respectively. The results showed a much lower GI value in seeds of TaABI5-A4a-GFP than that in TaABI5-A4b-GFP lines, indicating that TaABI5-A4a-GFP transgenic rice lines had higher ABA sensitivity than TaABI5-A4b-GFP (Fig. 6).
To understand the expression level of TaABI5-A4 on abiotic stress response, we examined the relative gene expression in TaABI5-A4a-GFP and TaABI5-A4b-GFP transgenic lines under different stress conditions. The plants at the three-leaf stage were treated with ABA (100 µM), salt (150 mM NaCl) and drought (300 mM mannitol to mimic dehydration) at 0, 6, 12, 24, 48 h post treatments (Fig. 7). The expression levels of TaABI5-A4 in TaABI5-A4a-GFP transgenic rice were higher than those in TaABI5-A4b-GFP after 6 or 12 h treatment with NaCl, whereas the expression levels of TaABI5-A4 in the TaABI5-A4a-GFP lines were lower than those in the TaABI5-A4b-GFP lines after 12–48 h treatment with manntiol. After 24 h in the ABA treatment, the expression levels of TaABI5-A4 in TaABI5-A4a-GFP lines were higher than those in TaABI5-A4b-GFP, which is consistent with ABA sensitivity tests, indicating that the TaABI5-A4a-GFP transgenic plants were more sensitive to exogenous ABA than the TaABI5-A4b-GFP.
Effect of decreased accumulation of spliced mRNA caused by nonsense-mediated decay (NMD) on phenotypic changes
Although the sequences of TaABI5-A4a and TaABI5-A4b have premature nonsense codons at the position of 266th amino acid, TaABI5-A4a-GFP and TaABI5-A4b-GFP transgenic rice lines had significantly different GI values, ABA sensitivities, and contents of endogenous ABA and GA. TaABI5-D3 encodes a full-length protein and its DNA sequence has the most similarity (99.07%) with the sequences of TaABI5-A4a and TaABI5-A4b; therefore, TaABI5-D3 was chosen as the positive control to verify whether the premature nonsense codon affects the transcription level and even the translation level. The GFP transcription level and TaABI5-A4-GFP fusion protein level in transgenic rice lines were analyzed with qPCR and Western blotting (WB). In the two-week-old leaves of transgenic rice plants, the transcription expression level of GFP was very much higher in the positive control with the TaABI5-D3-GFP than those in TaABI5-A4a-GFP and TaABI5-A4b-GFP lines (P < 0.0001), and the expression in TaABI5-A4a-GFP lines were significantly (P < 0.01) higher than that in TaABI5-A4b-GFP lines (Fig. 8a). This result indicated that the much lower accumulation of fully spliced RNA in the TaABI5-A4a/b-GFP transgenic lines was caused by nonsense-mediated decay (NMD). Moreover, WB analysis was performed using an Anti-GFP monoclonal antibody to detect the expression level of GFP fusion protein in transgenic rice seeds. The result showed that GFP fusion protein was only detected in the positive control, rather than in transgenic TaABI5-A4a-GFP and TaABI5-A4b-GFP lines (Fig. 8b). These indicated that it was decreased accumulation of spliced mRNA, but fusion protein, that led to phenotypic changes between TaABI5-A4a-GFP and TaABI5-A4b-GFP transgenic rice lines.
Effect of mRNA accumulation of TaABI5-A4 in transgenic rice on the expression of level ERF13 and GSL-OH associated with seed dormancy
In order to identify the downstream genes regulated by transcription factor TaABI5, we performed RNA sequencing (RNA-seq) using mature seeds of TaABI5-A4a-GFP and TaABI5-A4b-GFP transgenic rice lines. A total of 68 DEGs were identified in TaABI5-A4a-GFP transgenic lines compared with TaABI5-A4b-GFP lines, including 26 up-regulated and 42 down-regulated genes (Fig. 9a). Compared with TaABI5-A4b-GFP transgenic lines, the up-regulated DEGs in TaABI5-A4a-GFP lines belong to the pathways of regulation for response to stress, abiotic stimulus and metabolic process based on GO and KEGG analysis, whereas the down-regulated DEGs in TaABI5-A4a-GFP transgenic lines are involved in the protein metabolic process (Figs. 9b, c and d). Among them, we found ERF113/At5g13330 (LOC_Os08g30100), GSTU16 (LOC_Os10g38360) and GSL-OH (LOC_Os03g08460) in the up-regulated DEGs, and ESD4 (LOC_Os01g16730) in the down-regulated DEGs (Fig. 9e), which may be involved in seed dormancy and germination.
We further examined the transcript expression levels of RAP2.6L, GSTU16 and GSL-OH with qPCR in mature seeds. The results confirmed significantly higher expression levels of these three genes in TaABI5-A4a-GFP transgenic rice lines than those in TaABI5-A4b-GFP (Fig. 9f). However, ESD4 could not be detected by qPCR due to the extremely low expression, which was a significantly down-regulated gene in RNA-seq. These results further verified the up/down-regulated DEGs in the RNA-seq from the mature seeds of TaABI5-A4a-GFP and TaABI5-A4b-GFP transgenic lines (Fig. 9).