Deletion of the low-molecular-weight glutenin subunit allele Glu-A3a of wheat (Triticum aestivumL.) significantly reduces dough strength and breadmaking quality

Identification and characterization of seed proteins in CS and the Glu-A3deletion line CS-n

A Glu-3 deletion line of CS was screened and developed in our laboratory, and named CS-n. Compared to CS, the morphological characteristics of plants, spikes, and seeds, as well as the growth and development traits of CS-n showed no significant differences (Additional file 1: Figure S1, Additional file 2: Figure S2, and Additional file 3: Table S1). The grain protein compositions of CS and CS-n were identified by using various proteome approaches (Figure 1 and Additional file 4: Figure S3). The results indicated that CS-n showed the same albumin and globulin compositions as CS, while gliadins displayed minor differences between CS-n and CS; only one gliadin band obtained by acidic polyacrylamide gel electrophoresis (A-PAGE) was absent in CS-n (Additional file 4: Figure S3).

To obtain the accurate molecular mass of the Glu-A3a-encoded B-subunit, the expected protein band on the SDS-PAGE gel indicated in Figure 1a was collected and then analyzed by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS). As shown in Additional file 5: Figure S4, the Glu-A3a-encoded LMW-GS B-subunit was easily identified, and its molecular mass was determined to be 41,701.2 Da.

Confirmation of Glu-A3adeletion in CS-n with a sequence-tagged site polymerase chain reaction (STS-PCR) marker

To further confirm the deletion of the Glu-A3 locus in CS-n, a pair of STS primers developed from the single nucleotide polymorphisms (SNPs) in Glu-A3 allelic variants [24] were used to amplify the Glu-A3a gene. As shown in Figure 2, one specific PCR product of 529 bp was amplified in CS, the CS-1S^l/1B substitution line, the CS-1S^l addition line, and Aroona, which contain the Glu-A3a allele, whereas no such fragments were obtained in the other materials without Glu-A3a, such as CS-n. The specific amplified 529-bp fragment was sequenced, and the sequence was the same as those from the upstream 140–395 bp of the Glu-A3a-coding sequence shown in Additional file 6: Figure S5. Thus, these results confirmed that the Glu-A3 locus was deleted in CS-n.

Comparison of gluten quality properties between CS-n and CS

Farinograph analysis indicated that development time, stability time, tolerance index, and farinograph quality number in CS-n were significantly lower than those in CS (Table 2). These properties led to a decrease in loaf volume of CS-n from 760 to 735 cm³ (Table 2 and Figure 3). Bread texture analysis showed that the hardness and resilience of bread in CS were superior to those in CS-n. Further cell size analysis of the bread demonstrated that the quality in CS-n was significantly reduced (Table 3). For example, wrapper length, slice brightness, and wall thickness of CS-n bread slices were much lower than those of CS. The cell diameter and elongation in CS-n were also reduced as a result of Glu-A3a deletion.

Molecular characteristics of the LMW-GS allele Glu-A3a

To further understand the molecular mechanisms underlying the significant effects of Glu-A3a on gluten and breadmaking quality, the complete coding sequence of Glu-A3a was amplified and sequenced by allelic-specific (AS) PCR. Based on the previously characterized Glu-A3 genes, a pair of specific primers (A3-F and A3-R) for the Glu-A3 locus was designed and used to amplify the Glu-A3a allele from CS. As shown in Additional file 7: Figure S6, a single band of approximately 1100 bp was obtained from CS, whereas no product was amplified from CS-n. Since most of the complete coding sequences of LMW-GS genes vary in length between 909 and 1167 bp [6],[27]-[29], the size of the amplified band corresponded well to the known LMW-GS gene sizes. After sequencing of the amplified product, a complete open reading frame of 1134 bp was obtained. Sequence alignment showed that the cloned gene had no internal stop codons and contained typical structural features of LMW-GS, and therefore was named as Glu-A3a (Additional file 7: Figure S6). After searching the GenBank database, we found that the cloned Glu-A3a gene had the same sequence as GluA3-11 from cultivar Aroona-A3a (GenBank accession number FJ549928). The deduced amino acid sequence of Glu-A3a showed the presence of an isoleucine as the first amino acid residue in the N-terminal of the mature protein, indicating that it belongs to the LMW-i type subunit [6].

The deduced amino acid sequence of Glu-A3a had 376 amino acid residues with a predicted molecular mass of 41,346.1 Da, corresponding well to that determined by MALDI-TOF-MS (41,701.2 Da). Multiple alignment of the deduced amino acid sequences of Glu-A3a with the other 14 LMW-i type subunits (Figure 4) showed that all have conserved signal peptides and four domains in the mature protein sequences, including a repetitive domain, cysteine-rich region, glutamine-rich region, and C-terminal conservative region, as reported by Cassidy et al. [27]. Similar to other LMW-i type subunits, the Glu-A3a-encoded subunit contained eight cysteine residues at relatively conserved positions (Additional file 8: Table S2). It is speculated that the first and seventh cysteines of the LMW-GS form the inter-molecular disulfide bond, while the rest form three intra-molecular disulfide bonds [30],[31].

The number of repeats present in the repetitive domain is mainly responsible for the length variation and the general hydrophilic character of LMW-GS [30]. The Glu-A3a-encoded subunit contained the typical repeat motif of LMW-GS: P_1–2FP/SQ_2–6. Our results showed that Glu-A3a has a rather large and regular repeated sequence domain that includes a high proportion of glutamine residues (about 46%) in the repeats (consensus sequence PPFSQQQQ), and two polyglutamine stretches with 11 and 12 continuous glutamine residues in the repetitive and C-terminal domains, respectively. Repeat motif numbers in LMW-i subunits are much higher than those in the LMW-m and LMW-s subunits, ranking them the longest protein subunits among all Glu-3 loci.

Secondary structure and function prediction of the Glu-A3a-encoded protein

Phylogenetic analysis of Glu-A3aand other LMW-GS genes

A homology tree was constructed to reveal the phylogenetic relationships among 25 LMW-GS genes at Glu-3 loci from different species and genomes through nucleotide sequence alignment of their coding regions using MEGA5 software (Figure 5). These sequences comprised 21 LMW-GS genes from different genomes of Triticum diploid, tetraploid, and hexaploid species. The phylogenetic tree displayed two clear branches, which corresponded well to distinguishing the LMW-i type from the LMW-m and LMW-s type subunits. This demonstrated that LMW-i type genes have undergone greater divergence during evolution compared to LMW-s and LMW-m genes, as previously reported [38],[39]. Sine LMW-m and LMW-s type subunit genes generally show higher consistency, they showed close phylogenetic evolutionary relationships. Glu-A3a showed a closer relationship with other LMW-i type genes from common wheat. All of the LMW-i type subunit genes from common wheat and related species shared higher sequence identity, indicating their high evolutionary conservation.

Heterologous expression of Glu-A3a in Escherichia coli and determination of the corresponding native protein encoded by Glu-A3a

The Glu-A3a-coding region without signal peptides was expressed in E. coli. The expressed fusion protein was separated by both SDS-PAGE and 2-DE, and was further identified by LC-MS/MS. SDS-PAGE identification (Figure 6a) indicated that the relative mobility of the expressed protein was the same as that of the native Glu-A3a-encoded subunit of CS, confirming that Glu-A3a without the N-terminus can be expressed normally, similar to other LMW-i type genes [13]. Furthermore, 2-DE separation of the expressed protein (Figure 6b) demonstrated a similar pattern as that shown in Figure 1b. LC-MS/MS identification also confirmed that the expressed protein was the Glu-A3a-encoded subunit present in CS, as revealed by the previous tandem MS results (Table 1).

To verify the authenticity of the cloned sequence, LC-MS/MS was conducted by using the native Glu-A3a subunit digested by trypsin. We compared the results of LC-MS/MS of the SDS-PAGE band of CS, the heterologous protein, 2-DE spots, and the amino acid sequence of the Glu-A3a gene. This gives a coverage rate of 18.26% (65/356 amino acids of the mature polypeptide). These results revealed consistency in the peptide sequences among the samples, confirming the correspondence of the Glu-A3a gene and its native encoded subunit.

Dynamic expression profiles of the Glu-A3agene and its encoded protein during grain development

The dynamic transcription expression profiles of the Glu-A3a gene at 5, 11, 14, 17, 20, 23, 26, and 29 days post anthesis (DPA) of grain development were detected by quantitative real-time (qRT)-PCR in both CS and CS-n. Real-time melting temperature curves for the gene showed a single peak. qRT-PCR efficiency was determined by five serial five-fold dilutions of cDNA, and the standard curve confirmed high RT-PCR efficiency rates (Additional file 9: Figure S7). As shown in Figure 7a, the Glu-A3a gene displayed an up-down expression pattern during grain development of CS, with peak expression occurring at 14 DPA. However, Glu-A3a mRNA could not be detected in CS-n, further confirming the deletion of the Glu-A3 locus. SDS-PAGE analysis showed that the Glu-A3a-encoded B-subunit exhibited a gradual up-regulated expression pattern, and it began to rapidly accumulate after 11 DPA (Figure 7b). At 5 DPA, no LWM-GS genes could be detected, and both LMW-GS and HMW-GS showed trace expression levels. After 14 DPA, the B-subunit as well as other LMW-GS and HMW-GS genes displayed significant up-regulation, and peak expression occurred at 17 DPA (Figure 7b).

Development and validation of an SNP-based molecular marker for Glu-A3a

An AS-PCR marker was developed based on the SNPs detected in Glu-A3 genes. A pair of specific primers for Glu-A3a (Glu-A3a F: GCAAAGAAGGAAAAGA GGTGG, R: GGTTGTTGTTGTTGCTGCA) was designed and tested in different genotypes and hybrid generations with different Glu-A3 alleles. The materials with different Glu-A3 alleles included 48 bread wheat cultivars, the CS-1S^l/1B substitution line, and the CS-1S^l addition line, as well as seven Aroona NILs and four recombinant inbred lines (RILs) derived from a cross between the CS substitution line CS-1S^l/1B with Glu-A3a and the bread wheat cultivar CB037A with Glu-A3c (Additional file 10: Table S3). The Glu-A3 allele compositions of all materials used were identified by SDS-PAGE (Figure 8a). The PCR results showed that one specific PCR product of 507 bp was amplified in all cultivars with Glu-A3a (Figure 8b). To validate the effectiveness of the STS marker, seven NILs and four RILs with different Glu-A3 allele compositions were used for PCR amplification. The results showed that the 507-bp fragment could be specifically amplified in the lines with the Glu-A3a allele, whereas no any amplification products were obtained from the lines with other Glu-A3 alleles, including CS-n without the Glu-A3 locus. These results confirmed that the developed AS-PCR marker could be used as an effective tool for rapidly screening the Glu-A3a allele in wheat quality improvement strategies through molecular marker-assisted selection.

Allelic variations at Glu-3loci and their effects on gluten quality

LMW-GS account for approximately 60% of glutenin proteins in mature seeds and play important roles in the formation of glutenin macropolymer and gluten quality [1],[45], particularly for dough extensibility and strength [3]-[6],[17]. LMW-GS genes belong to a multiple gene family and are found in multiple copies in Triticum aestivum; the copy number in hexaploid bread wheat was estimated to vary from 10–15 [46] to 35–40 [27],[47]. A recent study based on BAC library screening and proteomics analysis showed that Glu-A3, Glu-B3, and Glu-D3 in the Chinese bread wheat cultivar Xiaoyan 54 contain 4, 3, and 7 genes, respectively [16]. In addition, by using the LMW-GS gene marker system, at least 15 LMW-GS genes were identified in Aroona NILs [17].

Glu-A3 and Glu-B3 alleles are known to play a major role in determining differences in processing qualities among the three Glu-3 loci, while Glu-D3 alleles play minor roles in determining quality variation in bread wheat [17]. In particular, the Glu-A3 locus was considered to have the biggest contribution to quality among all LMW-GS loci, in which Glu-A3f was found to have a strong positive effect on end-use quality [48]. In Australian wheat cultivars, LMW-GS provided better predictions of Rmax than HMW-GS [45]. The effects of different Glu-3 alleles on Rmax showed the following ranking: Glu-A3b > Glu-A3d > Glu-A3e > Glu-A3c, Glu-B3i > Glu-B3b = Glu-B3a > Glu-B3e = Glu-B3f = Glu-B3g = Glu-B3h > Glu-B3c, and Glu-D3e > Glu-D3b > Glu-D3a > Glu-D3c > Glu-D3d [17]. However, no studies of the effect of the Glu-A3a allele on gluten quality have been reported so far. In the present work, we found that the deletion of Glu-A3a significantly reduced dough strength and breadmaking quality, including most of the mixing and bread quality parameters (Tables 2 and 3). This indicates that Glu-A3a plays important roles in conferring high gluten quality to wheat.

Molecular basis of the relationship between Glu-A3aand gluten quality

The molecular structures of LMW-GS proteins play important roles in determining the dough strength and gluten quality; in particular, the distribution of cysteine residues could lead to functional protein differences [6]. The first and the seventh cysteines form the inter-molecular disulfide bond, while the remaining cysteines form three intra-molecular disulfide bonds [11],[30],[31]. Thus, the number and position of cysteines are important to the formation of the secondary protein structure and, consequently, dough quality. The presence of a long repetitive domain is also considered to have a positive influence on wheat flour quality [30],[49]. A repeated sequence domain could be helpful in increasing the viscosity and elasticity of the dough through increasing the inter-molecular interactions among the large number of glutamine side chains, which are both good hydrogen bond donors and acceptors [49]-[51]. According to Masci et al. [30], all the predicted α-helices in the 42 K LMW-GS seem to be located near the intra-molecular disulfide bonds. They also speculated that helix-helix interactions are involved in guiding the formation of the intra-molecular disulfide bonds. Therefore, a higher α-helix content may contribute to better quality of the dough [52]. The β-strands are generally considered to endow the protein with high elasticity and to improve the capability to resist distortion [38].

LMW-i type subunits contain eight highly conserved cysteine residues in the C-terminal domain (Table 5), which differ from the LMW-m and LMW-s subunits that have a cysteine residue in the N-terminus or in the repetitive domain [39]. Therefore, the secondary structures of LMW-i subunits are expected to be quite different from those of the LMW-m and LMW-s subunits. Previous work on LMW-GS AY542896 and AY263369 indicated that they have positive effects on quality properties [13],[37]. In the present study, comparative analysis with the secondary structures of AY542896 and AY263369 showed that the Glu-A3a-encoded subunit had more α-helices (Table 4). The secondary structure is the foundation for a highly complex spatial conformation and is composed of structural motifs, including α-helices, β-strands, and random coils. The higher α-helix number in the Glu-A3a-encoded subunit could guide the formation of the intra-molecular disulfide bonds and contribute to superior dough strength and breadmaking quality.

The sizes of most of the cloned LMW-GS genes range from 900 to 1000 bp, and the gene LMW-TD22 with 1167 bp is the longest complete gene among the cloned LMW-i genes analyzed to date [29]. The presence of a long repetitive domain is considered to have a positive influence on wheat flour quality because it can form more β-turns in the structure, thereby conferring elasticity to the protein molecule [30],[52]. The molecular structure of the deduced LMW-TD22 subunit indicated a long repetitive domain of 21 repeat motifs (consensus sequence P_1–2FP/SQ_2–6). In the present study, the Glu-A3a-encoded subunit was also found to have a long repeated sequence domain and a high proportion of glutamine residues (about 46%), which could improve the conformation of superior gluten structure and breadmaking quality.

Glu-A3aexpression and LMW-GS synthesis

It is well known that the expression levels as well as accumulation patterns of storage proteins are closely associated with gluten quality properties [53]-[55]. For example, the over-expressed HMW-GS 1Bx7^OE has positive effects on dough strength [53],[56]. In addition, the accumulation rates vary between different groups of proteins, suggesting differential regulation of protein biosynthesis and different quality performance. In particular, the wheat biotype with superior HMW-GS 5 + 10 subunits accumulated larger polymers more quickly than the biotype with poor allelic subunits 2 + 12 [54].

The B-subunits of LMW-GS are the most abundant and have the greatest impact on wheat processing qualities [6]. In this work, RP-UPLC analysis revealed a higher expression level and greater proportion of Glu-A3a-encoded B-subunits, accounting for more than 22% of the total LMW-GS in CS (Figure 1c), indicating its major contribution to LMW-GS synthesis and its important roles in determining dough quality. A recent study also found that higher numbers of active LMW-GS genes at the Glu-A3 and Glu-D3 loci in Xiaoyan 54 tended to produce greater ZSVs, an important indicator of breadmaking quality [16]. Similarly, the decrease in the number of active LMW-GS genes in CS-n due to deletion of the Glu-A3 locus likely contributed to the significant reduction in dough strength and breadmaking quality (Table 2 and 3).

Wheat glutenin proteins generally display an up-regulated expression pattern during grain development (Figure 7a). Coordinated accumulation of transcripts from HMW-GS and LMW-GS genes, as well as α-, γ-, and ω-gliadin genes, occurs early in grain development [54],[57]. LMW-GS, HMW-GS, and ω-gliadins can be detected by gel electrophoresis as early as 7 DPA [54], and 10–18 DPA represents the key stage of storage protein synthesis [58]. In the present study, Glu-A3a transcripts demonstrated an up-down expression pattern during grain development, and the highest expression level occurred at 14 DPA (Figure 7a), similar to a previous report [59]. SDS-PAGE analysis revealed that the Glu-A3a-encoded B-subunit displayed an up-regulated expression pattern and showed rapid synthesis and accumulation at 11–17 DPA (Figure 7b), which is also generally in agreement with a previous report [55]. Thus, the Glu-A3a-encoded B-subunit has a higher accumulation rate during grain development similar to HMW-GS 5 + 10 [54], which could improve the conformation of the regular gluten structure. Some important genes related to storage protein folding and synthesis, such as protein disulfide isomerase (PDI) and binding protein (BiP) genes, generally have higher expression levels at the early grain developmental stages. For instance, the PDI genes PDIL1-1 and PDIL2-1, which are involved in disulfide bond formation, displayed a peak expression level in the early stages (about 10–15 DPA) of grain development [58]. The higher accumulation rate of the Glu-A3a-encoded B-subunit was accompanied by higher expression levels of the genes involved in storage protein synthesis and assembly during early grain development, suggesting that this subunit could improve the conformation of gluten macropolymers (GMP) and result in superior dough quality.

Potential application of Glu-A3ain wheat quality improvement through molecular marker-assisted selection

Characterization of the allelic variations of LMW-GS is important for improvement of wheat-processing quality. Some allelic variations of LMW-GS have greater positive effects on dough properties than others [3],[45],[60]. Marker-assisted selection is an effective supplement to conventional breeding practices. For LMW-GS, because of the low resolution of traditional SDS-PAGE and the tedious operation procedures involved in 2-DE, development of different markers is important for the study and application of target subunits.

Recently, with increasing numbers of LMW-GS alleles being cloned and sequenced from common wheat, different molecular markers have been developed to rapidly screen and select desirable Glu-3 alleles. Zhang et al. [41] developed a set of markers that can be used to discriminate the alleles Glu-A3a, b, c, d, e, f, and g. Long et al. [61] classified 69 LMW-GS genes registered in GenBank into nine groups and established nine group-specific primer sets to identify each group. Ikeda et al. [62] developed 12 specific PCR markers to distinguish 12 groups of LMW-GS genes in the wheat cultivar Norin 61. Ten allele-specific STS markers for Glu-D3 were developed by Zhao et al. [43],[63],[64]. Wang et al. [35] designed 10 allele-specific PCR markers for the Glu-B3 locus based on SNPs present in the sequences of different allelic variants. Wang et al. [24] reported an allele-specific marker for Glu-A3b, which was not reported by Zhang et al. [41].

In the present study, we developed a new allele-specific PCR marker that can effectively discriminate Glu-A3a from other Glu-A3 alleles, which was validated using different cultivars, including RILs and NILs (Figure 8b). This Glu-A3a allele-specific marker can be used in marker-assisted breeding strategies aimed at the improvement of wheat quality. With help of this marker, it will be very convenient and effective for breeders to select this superior gene in early hybrid generations of a wheat quality program. However, this marker can only identify one specific gene in different generations and materials one time. The multiplex PCR systems showed to be more rapid and economic in identifying different desirable genes [24]. Thus, to improve the selection efficiency, it is needed to further develop multiplex PCR markers that can rapidly identify different desirable genes, including Glu-A3a and other quality-related genes.

Plant materials

Chinese Spring (Triticum aestivum L., 2n = 6x = 42, AABBDD) and its Glu-A3a deletion line (CS-n) developed in our laboratory were used in this study. Aroona-A3a (Glu-A3a) and its six NILs: Aroona-A3b (Glu-A3b), Aroona (Glu-A3c), Aroona-A3d (Glu-A3d), Aroona-A3e (Glu-A3e), Aroona-A3f (Glu-A3f), and Glenlea (Glu-A3g) were used for identifying Glu-A3 alleles, which were kindly provided by Dr. Xianchun Xia, Institute for Crop Science, Chinese Academy of Agricultural Science (CAAS). Bread wheat cultivars CB037A, Chinese Spring substitution line CS-1S^l/1B [58], Chinese Spring 1S^l genome addition line, and four recombination inbred lines (RILs) produced by crossing between CB037A and Chinese Spring substitution line CS-1S^l/1B were used for identifying Glu-A3a deletion in CS-n, and developing and validating AS-PCR marker for Glu-A3a. All materials used in this study were listed in Additional file 10: Table S3.

Identification of seed proteins

Gluten quality testing

Both CS and its CS-n were planted in Beijing experimental station of Chinese Academy of Agricultural Sciences during 2013–2014 growing season. The field design was three replications and each blot was 20 m². Gluten quality parameters were tested according to Sun et al. [73] at Academy of State Administration of Grain. Flour moisture and ash contents (% dry basis) were determined according to the American Association of Cereal Chemists Approved Methods (2000) 44-15A and 08–02, respectively. Protein content (%N 5.7, 14% moisture basis) was determined by nitrogen combustion analysis with a LECO (Model FP analyzer, St. Jopeph, MI) calibrated against EDTA.

Farinograph parameters were obtained by using 10 g Brabender Farinograh-E based on American Association of Cereal Chemists Approved Method (2000) 54–21.

Image analysis of crumb grain of bread was performed with a C-Cell, image analysing software and equipment (Calibre Control International Ltd.). Slice brightness and cell contrast were used to describe the brightness of slices. Number of cells, wall thickness, cell diameter, cell volume, coarse cell volume and average cell elongation were used to measure the cell properties. The sliced samples for textural analysis were prepared in the same way as those for the C-Cell.

DNA extraction and PCR amplification

Total genomic DNA from dry seeds was extracted according to McDonald et al. [74] and An et al. [39] with minor modifications. A pair of AS-PCR primers (A3-F and A3-R) was designed to amplify the coding regions of LMW glutenin gene based on the previously cloned sequences [13],[24]. The primer sequences were A3-F: 5’-GCCTTTCTTGTTTACGGCTG-3’, A3-R: 5’-TCAGATTG ACATCCACACAAT-3’ (synthesized by Sangong Inc., China). PCR amplifications were performed in 50 μl reaction volumes containing 2.5 U La Taq polymerase (TaKaRa), 100 ng of templet DNA, 25 ml of 2 × GC buffer I (MgCl₂⁺ plus), 0.4 mM dNTP, 0.5 μM of each primer, and double distilled H₂O added to 50 μl. The reactions were carried out in a PTC-100 (MJ Research, Watertown, MA, USA) thermocycler using the following protocol: 94°C for 2 min, followed by 35 cycles of 94°C for 45 s, 58°C for 70 s and 72°C for 2 min, finally extended at 72°C for 10 min.

Molecular cloning, DNA sequencing and sequences alignment

PCR products were separated on 1.2% agarose gels in Tris–acetic acid–EDTA buffer and expected fragments were purified from the gels using a Quick DNA extraction kit (Tiangen, Beijng, China). Subsequently, purified products were ligated into a PMD18–T Easy vector (TaKaRa, Dalian, China) and transformed into cells of E. coli strain DH5α according to Li et al. [75]. DNA sequencing was performed with three clones by Sino Geno Max, Beijing, China. Multiple sequence alignment of LMW glutenin nucleotide and protein sequences were completed by Bioedit 7.0.1.1.

SNPs and InDels identification and secondary structure prediction

Identification of SNPs and InDels present in LMW glutenin genes were based on multiple alignments and performed using Bioedit 7.0. Prediction of secondary structure of deduced amino acid sequences was carried out by PSIPRED server (http://bioinf.cs.ucl.ac.uk/psipred/) [34].

Phylogenetic analysis

MEGA 5 was used to construct a phylogenetic tree with the complete coding regions. Neighbor joining with Kimura two parameter correction methods and bootstrapping of 1,000 replicates were selected as working parameters [28],[34],[39],[76],[77].

Expression of the cloned LMW-GS gene in E. coli

The gene cloned was re-amplified to remove the signal peptides by designing a new pair of primers CS-F (5’-GGGCATATGATTTCACAGCAACAA-3’) and CS-R (5’-CTCGAGTCAGTAGACACCAACTCCGATG-3’), NdeI and XhoI sites (underlined) were incorporated into the 5’ ends of the CS-F and CS-R, respectively. After purification, the PCR products were ligated into the expression vector pET-30a (Novagen), and transformed into E.coli BL21 (DE3) plysS cells. And then we extract and separate the expressed protein from the E.coli, after that, we carried out them by SDS–PAGE according to Li et al. [75].

mRNA extraction, cDNA synthesis and qRT-PCR

Developmental seeds from three spikes were combined together to extract total RNA from endosperm of CS and CS-n, and cDNA synthesis, qRT-PCR were according to Wang et al. [56]. The primers were: LMW-i-F: TGAAGACCTTCCTCGTCTTTG, LMW-i-R: CTGTGAAATTTGCGCAACG. Gene-specific primers were designed using Primer 5.0 and their specificities were checked by the melting curves of the RT-PCR products. Each qRT-PCR reaction was performed in 20 μl volumes containing 10 μl 2 × SYBR® Premix Ex Taq™ (TaKaRa), 2 μl 50-fold diluted cDNA, 0.4 μl of each gene-specific primer and 7.2 μl ddH2O. PCR conditions were as follows: 95°C for 3 min, 45 cycles of 15 s at 95°C, 57°C for 20s and 72°C for 20s. Three replicates were used for each sample. Reactions were conducted in a CFX96 Real-Time PCR Detection System (Bio-Rad). All data were analyzed with CFX Manager Software (Bio-Rad).

Determination of Glu-A3adeletion in CS by STS-PCR marker

To identify the deletion of Glu-A3a in CS-n, we did the STS-PCR marker of the seven NILs, CS-1S^l/1B, S genome addition line, CB037A, CS, CS-n. We used the marker name of glu-A3a to do PCR as what Wang et al. [24] did before. The primer sets are LA1F: AAACAGAATTATTAAAGCCGG, and SA1R: GGTTGTTGTTGTTGCAG CA. Their PCR cycling conditions were 94°C for 4 min, followed by 35 cycles of 94°C for 35 s, 55°C for 45 s, 72°C for 40 s, and a final extension at 72°C for 8 min.

Development and validation of allele-specific PCR markers for Glu-A3a

To identify the Glu-A3a gene in different genotypes, based on the SNPs we detected in Glu-A3a, we designed the primer named Glu-A3a F: GCAAAGAAGGAAAAGAG GTGG, R: GGTTGTTGTTGTTGCTGCA as the primer to discriminate the gene Glu-A3a from others in CS and CS-n, this also validated in four RILs (CB037A and 1S^l/1B), CS-1S^l/1B, 1S^l genome addition line, 7 NILs of Aroona and 48 varieties. PCR cycling conditions were 94°C for 4 min, followed by 35 cycles of 94°C for 35 s, 60°C for 30 s, 72°C for 30 s, and a final extension at 72°C for 8 min.