Sucrose phosphate synthase (SPS) is an important component of the plant sucrose biosynthesis pathway. In the monocotyledonous Poaceae, five SPS genes have been identified. Here we present a detailed analysis of the wheat SPSII family in wheat. A set of homoeologue-specific primers was developed in order to permit both the detection of sequence variation, and the dissection of the individual contribution of each homoeologue to the global expression of SPSII.
The expression in bread wheat over the course of development of various sucrose biosynthesis genes monitored on an Affymetrix array showed that the SPS genes were regulated over time and space. SPSII homoeologue-specific assays were used to show that the three homoeologues contributed differentially to the global expression of SPSII. Genetic mapping placed the set of homoeoloci on the short arms of the homoeologous group 3 chromosomes. A resequencing of the A and B genome copies allowed the detection of four haplotypes at each locus. The 3B copy includes an unspliced intron. A comparison of the sequences of the wheat SPSII orthologues present in the diploid progenitors einkorn, goatgrass and Triticum speltoides, as well as in the more distantly related species barley, rice, sorghum and purple false brome demonstrated that intronic sequence was less well conserved than exonic. Comparative sequence and phylogenetic analysis of SPSII gene showed that false purple brome was more similar to Triticeae than to rice. Wheat - rice synteny was found to be perturbed at the SPS region.
The homoeologue-specific assays will be suitable to derive associations between SPS functionality and key phenotypic traits. The amplicon sequences derived from the homoeologue-specific primers are informative regarding the evolution of SPSII in a polyploid context.
Sucrose plays an important role in the plant life cycle. As the major photosynthetic product, it is essential for growth, the synthesis of biomass and as a carbon and energy source. In the cereals, it is converted into starch and storage proteins, while many dicotyledonous species use it to form lipids and/or storage proteins [1, 2]. Under stressful conditions (e.g, low temperature or drought ), the plant cell typically accumulates sucrose as a protective osomoticum [3, 4]. Sucrose phosphate synthase (SPS) is one of the main regulatory enzymes involved in sucrose biosynthesis pathway in wheat and many other crop species . It catalyzes the conversion of Fructose-6-Phosphate and UDP-glucose into Sucrose-6-phosphate, which is subsequently hydrolysed to sucrose by the action of sucrose phosphate phosphatase (SPP).
In important crop plants like maize, rice and sugarcane, plant growth and productivity have been correlated with SPS activity. SPS activity has been correlated with sucrose accumulation in sugarcane stems [6–8], while in maize, correlations have been demonstarted with vigour and biomass yield [9, 10]. Genetic studies in maize have also shown that grain yield QTL (quantitative trait loci) are linked to SPS activity QTL, as are ADP-glucose pyrophosphorylase (another starch biosynthesis enzyme) activity QTL [11–13]. In rice, the location of plant height QTL appears to coincide with that of OsSPS1, as the transgenic plants with increased SPS activity grew taller than the non-transformed control . In tobaccco lines engineered to over-express SPS, UDP- glucose pyrophosphorylase and sucrose synthase, plant height was increased and flowering time delayed , while the heterologous expression of AtSPS induced longer stems and greater biomass . Finally, the over-expression of ZmSPS in potato was shown to improve a number of yield-realted characters , while in tomato, it altered the pattern of carbohydrate partitioning in the leaf .
The regulation of SPS activity is rather complex, and involves fine tuning at both the transcriptional and the post-translational level . Among the dicotyledonous plant species, three SPS gene families have been recognized, while in the monocotyledonous Poaceae species wheat, maize and rice, there are five gene families . Here, we report the temporal patterns of expression of SPS and SPP, and the discovery of structural polymorphisms between homoeologous copies of SPSII. We have also identified the presence of unspliced introns in one of its homoeologues. Evolutionary relationships among the SPSII homoeologues have been illustrated by comparisons between the gene sequences present in hexaploid wheat and its progenitors. Finally, we extend the phylogeny to other grass species.
Gene expression analysis of sucrose biosynthesis genes during wheat plant ontogeny
As a first step, the wheat Affymetrix probe sets representing SPS and SPP members were identified, based on the criteria used to define the five SPS  and the three SPP  types. The expression patterns of these gene family members was extracted from the normalized expression data from imbibed embryo, coleoptile, seedling root, crown, seedling leaf, stem (at anthesis), developing anther, developing caryopsis (sampled at 3-5 dap, 8 dap, 18 dap and 22 dap), and the developing endosperm and embryo at 22 dap. Except for SPPIII, the SPP genes were expressed throughout plant development, although their level tended to be highest in the developing anther and developing caryopsis (Figure 1). Different SPS members show highly regulated spatio-temporal expression patterns during plant ontogeny, suggesting the existence of fine regulation at the first step of sucrose biosynthesis.
As SPSI was prominently expressed in the seedling leaf and the mature stem (Figure 1) its function may be associated with the re-synthesis of sucrose from fructan. SPSII (Ta.4256.1.S1_at), SPSIII (Ta.20568.1.S1_at), TaSPSIV and TaSPSVa were prominently expressed during microsporogenesis (at premeiosis, throught meiosis and in immature pollen) and in the mature anther. All the SPP families were highly expressed during anther development (Figure 1). Thus, these genes are clearly important for normal pollen/anther development. Interestingly, note that in the Arabidopsis thaliana mutant kns2-1, an altered SPS glycosyltransferase I domain disrupts the function of sucrose synthesis, and this leads to changes in the synthesis of the pollen primexine or callose . In the wild type, SPS5.2 (KNS2) is expressed in uni- and bicellular microspores as well in more mature anthers. The present data showed that probably all the SPS and SPP genes were constitutively expressed throughout caryopsis development, indicating that more than one SPS isoforms contributes to overall SPS activity (Figure 1).
Furthermore, it is important to know how these different SPS members are spatially expressed in different seed tissues and possible co-expression of these members in the same cell for making hetero-oligomeric complexes to be functional. The developing caryopsis imports sucrose to driving its growth and allow storage of metabolites. The expression of SPS and SPP members within the developing caryopsis suggests that sucrose cleavage, transport and re-synthesis during its internal transport can be differentially regulated within various distinct tissues. Of particular interest was the observation that both variants of SPSII are differentially expressed in the developing caryopsis. Thus, Ta.4256.1.S1_at was highly expressed both throughout anther development and during late embryo development, while TaAffx.8003.1.S1_at was expressed most strongly during early caryopsis development. The contribution of the individual SPSII homeologues and the occurance of alternative splicing are expanded upon in the following section. Besides alternative splicing, a contributory reason for the present exclusive focus on the SPSII family was that there were four SPSII cDNA sequences available in the public domain, which gave the opportunity to design homoeologue-specific primers for three wheat genomes. Insufficient sequence data are presently available for the other SPS families. Therefore, in order to study polyploidization and phylogenetics relationship between three wheat genomes, SPSII was selected as best suited candidate.
Structural analysis of SPSII genes in wheat and its relatives
Both introns and exons were targeted for the design of homoeologue-specific primers (Figure 2). The resulting amplicons allowed for the sequencing of portions of each of the SPSII homoeologues. In all, 7,364 bp of the cv. Chinese Spring A genome homoeologue sequence was acquired from a combination of homoeologue-specific amplicons and flow-sorted chromosome arm 3A sequences (kindly provided by Dr. E. Akhunov). We acquired 3,883 bp of the B genome homoeologue sequence, and 1,184 bp of the D genome homoeologue, both from cv. Chinese spring (Additional file 1). Sequence of the SPSII orthologues in Triticum urartu (einkorn wheat - genome AA; A genome progenitor), T. speltoides (genome SS ~ BB; putative B genome progenitor), Aegilops tauschii (goatgrass-D genome progenitor) and barley (genome HH) was obtained in a similar fashion. Their alignment is provided as Additional file 2. SPSII has 12 introns and 13 exons in wheat, barley, rice, sorghum and purple false brome (Brachypodiumdistachyon) (Additional file 1, Figure 3). In rice, sorghum and purple false brome, the gene length is 12,838 bp, 12,456 bp and 8,492 bp, respectively (Additional file 1). In einkorn and T. speltoides six exons and six introns were amplified, while in goatgrass, this was extended to eight exons and seven introns. In hexaploid wheat, certain regions could not be sequenced either due to the presence of large introns or due to unavailability of homoeologue-specific primers. The intronic sequences were more polymorphic than the exonic ones, which were highly conserved among the species compared (Additional files 3, 4, 5, 6). In the six exons for which sequence was available, the level of similarity level ranged from 83.1% (between barley and sorghum) to 99% (between the wheat A genome homoeologue and einkorn) (Additional file 6b). In contrast, for the six intronic regions compared, the similarity level ranged from 35.5% (between purple false brome and sorghum) to 99% (between the wheat A genome copy and that of einkorn) (Additional file 6c). Intron 12 of the wheat B genome homoeologue was unspliced, since it was present in the cDNA AF347066 . Sequence alignment and subsequent cluster analysis showed that the purple false brome sequence was more similar to those of wheat and barley in comparison to rice sequence, and the sorghum sequence was the one most highly diverged from that of wheat (Figure 4), as was the expectation on taxonomic grounds .
Genetic mapping of TaSPSII genes in wheat
The homoeologue-specific primers developed for each of the wheat SPSII copies were validated by an analysis of the nulli-tetrasomic lines of cv.'Chinese Spring'; this showed that each amplicon was derived from a member of the homoeologous group three chromosomes (Figure 2). A polymorphism for the A genome copy between the parents of the ITMI (International Triticeae Mapping Initiative) mapping population allowed it to be located 13.6 cM distal of Xtam61a in the sub-telomeric of the short arm of chromosome 3A (Figure 5). The B genome copy was placed in a matching location on chromosome 3B between Xksug53a and Xcdo460c (Figure 5). The D genome copy was not mappable in the ITMI mapping population, but was localised to D genome through anueploid stocks.
Detection of SNP diversity
The SPSII amplicons obtained from 27 diverse wheat accessions were sequenced in a search for single nucleotide polymorphisms (SNPs). For the A genome copy, the intron 6 to exon 12 region was targeted, while for the B genome, the target was exon 6 to the 3' UTR, and, for the D genome copy, intron 10 to exon 12. Three SNPs were uncovered in the 3A gene, all within an exonic region corresponding to the sucrose phosphate synthase/possible sucrose phosphate phosphatase domain. Of the three SNPs, two were non-synonymous SNPs; the first induced the residue change argeneine1554 to lysine1554, and the second from alanine1599 to valine1599 (Figure 3, Table 1). Of the four B genome copy SNPs detected, three were located within intronic sequences, and one within the 3' UTR (Figure 3, Table 1). No SNP was detected within the D genome copy. Overall, it was possible to define four distinct haplotypes for each of the A and the B genome copies (Table 1).
Haplotypes for SPSII gene obtained in wheat A and B genomes
Haplotypes -A Genome
SNP1 (Exonic, 1432 bp
SNP2* Exonic, 1554 bp)
SNP3* (Exonic, 1599 bp)
Cultivars per haplotype
SNP1 (Intronic, 451 bp)
SNP2 (Intronic, 523 bp)
SNP3, (Intronic 793 bp)
Cultivars per haplotype
Total = 27
Total = 27
Note: Position of the SNPs are based on sequence alignment (supplementary Figure 1). Asterisk showed the location location of non synonymous
Three genomes differentially contributes to the expression of TaSPSII
Quantitative real time PCR based on exonic sequences of the A and D genome copies, and within the unspliced intronic region of the B genome copy was applied to various tissues of, cv 'Prinz'. SPSII was expressed during germination as well as in the developing caryopsis at four and eight days after anthesis (DAF). A similar level of expression was measured in the germinating seed and in the developing caryopsis at four DAF, and this fell somewhat by eight DAF. The A genome copy dominated the global expression of SPSII; there was very little contribution from the B genome copy (Figure 6).
Genomic sequence analysis of SPSII genes
Sequence variation in sucrose biosynthesis pathway genes such as SPS is of particular interest because its consequence can be a change at the level of phenotype, sometimes involving vigour, biomass productivity and economic yield. Markers from the expressed portion of genome, that are responsible for some important physiological aspects, have been developed in wheat [22–26]. A particular complication in bread wheat is that most genes which are present in single copy in the diploid grasses such as rice and purple false brome are represented by three copies. In order to interpret wheat gene expression data, it is therefore necessary to derive specific assays for each homoeologous copy, as we have described here for SPSII (Table 2). These assays are also necessary for the application of pyrosequencing as a SNP discovery platform [27, 8].
Primer sequences used for the amplification of SPSII gene in wheat and other related grasses
*represents homoeologue-specific primers ; arepresents primer pairs used for quantitative real time PCR; # represents primer pairs used for genetic mapping. For details of PCR profile please see methods section
Frequency of SNPs
The low level of SNP within SPSII suggests that variation in its sequence is constrained by selection against changes in its physiological activity. Other starch metabolism enzyme encoding genes have similarly been associated with little sequence variation [22, 28]. The observed SNP frequency was 1 per 797 bp within the A genome copy, and 1 per 916 bp within the B genome copy, frequencies which are typical for bread wheat [29–31]. Most of the SPSII SNPs were located within either an intron or the 3' UTR, with only two exonic non-synonymous ones uncovered. These latter substitutions may be of physiological significance, given that an alanine to valine change compromised the processing of a 22 kDa protein in maize , and that an arginine to lysine substitution in the bZIP domain of the maize opaque gene abolished the specificity of its DNA binding .
Presence of an unspliced intron in TaSPSII cDNA is mapped to B genome
Intronic sequence was relatively variable at the inter-specific and inter-generic level (Additional files 4, 5c), while exonic sequence was well conserved (Additional files 3, 5a and 5b). This pattern is characteristic of gene sequences in the grasses [21, 34]. Analysis of the sequence of SPSII cDNA clone AF347066 has demonstrated alternative splicing for one of the SPSII genes . The development of homoeologue-specific assays has allowed the unspliced intron to be located to the B genome copy of SPSII. Both of the alternative transcripts (spliced and unspliced) occurred in the pre anthesis spike . cDNA of einkorn wheat (genome AA), goatgrass (genome DD), durum wheat (T. durum, genome constitution AABB) and bread wheat (genome AABBDD) were amplified using primer pairs w-spsII-iL & w-spsII-nR, which target the unspliced intron in the B genome copy of SPSII. Amplification was observed from bread and durum wheat template, but not from either einkorn or goatgrass (Figure 7b). An aneuploid analysis confirmed the B genome specificity of the unspliced intron (Figure 7c). Thus the unspliced intron is clearly restricted to the B genome copy of SPSII.
Contribution of A and D genome for the expression of TaSPSII
All five SPS members appeared to be both spatially and temporally regulated. Of special interest in SPSII with two different variants, Ta.4256.1.S1_at: highly expressed throughout anther development as well also during late seed development in embryo and TaAffx.8003.1.S1_at: expressed mainly during early seed development. The expression pattern of Ta.4256.1.S1_at was consistent with that previously exposed by Northern blot analysis , but the expression of the latter gene has not been documented before. Although microarray-based expression data reflect very clearly any differences between distinct genes, they are not informative at the level of homoeologues, because the array features are not homoeologue-specific. In this case, homoeologue specific primers could play a vital role in knowing the contribution of three different genomes in terms of gene expression pattern. Similar strategy was also used earlier for studying expression patterns of benzoxazinones biosynthesis in wheat contributed by different genomes . The present data have shown that the A genome copy of SPSII dominated the global expression of SPSII and that the B genome copy - despite its two alternative transcripts - the least.
Synteny to rice is perturbed at SPSII locus
At the cDNA sequence level, SPSII is highly similar to its rice orthologue, which is located on chromosome 8 . Gene order is largely conserved among the major cereal genomes [36–38], which has allowed for the identification of a syntenic relationship between rice chromosome 8 and the wheat chromosomes of homoeologous group seven. The prediction is therefore that wheat SPSII should map to wheat chromosomes 7A, 7B and 7D, and not, as is the case, to chromosomes 3A, 3B and 3D. The chromosomal location that we have deduced is secure, as it relies on a combination of anaeuploid analysis, genetic mapping, the placement of SPSII ESTs http://wheat.pw.usda.gov/cgi-bin/westsql/map_locus.cgi, and the presence of SPSII sequence on flow-sorted chromosome 3A template. Rice chromosome 1 has been reported to be highly syntenic to wheat chromosome three [38, 39]. In the present study, it appears that a segment at least of the group 3 chromosomes has lost synteny with rice chromosome 1. Failure of synteny and breakdown of colinearity is well documented [24, 38, 40, 41]. Prominent examples involve the regions around both GBSS1  and the Pbf orthologues in wheat and barley . Of 537 mapped ESTs mapping to the wheat group 3 chromosomes, ~15% shared homology with rice genes mapping outside of chromosome 1, and various complexities associated within co-linear regions suggests that the relationship between these chromosomes is not straightforward . Similarly, the analysis of 1 Gb of chromosome 3B has identified the presence of gene sequences whose rice orthologues do not map to chromosome 1 , while the annotation of 17 Mb of various regions of chromosome 3B has shown that only half of the predicted genes present had orthologues on rice chromosome 1 . The SPSII genes lie in the distal region of their respective wheat chromosome, which suggests that co-linearity with rice may be lost there. Complex duplications, rearrangements and polyploidization events have led to the erosion of synteny and colinearity between rice and wheat [45–47]. The imminent completion of the sequence of the entire chromosome 3B will allow the formulation of a better picture of orthologous relationships between wheat, rice, sorghum and false purple brome, as the genomes of later three species have already been fully sequenced.
Phylogenetic relationship of TaSPSII in the context of polyploidization among hexaploid wheat and its progeniotors
The sequence alignments and subsequent phylogenetic analyses was carried out considered either exons or introns or both (Figure 4, Additional files 2, 3, 4, 5, 6). The hexaploid wheat A genome copy consistently appeared closest to that of einkorn wheat, the A genome progenitor. At the level of exons only, the B genome copy appeared closely related to that of T. speltoides, a species considered to be taxonomically close to the B genome progenitor. Branching patterm in dendrogram was little disturbed when considered only introns depicting nonconserved nature of introns. With respect to the D genome copy, the comparison with its goatgrass orthologue was only possible on the basis of the sequence of one large exon. This showed that the D genome copy was more similar to the goatgrass orthologue than to any of the others. The level of homology in the exonic sequence between the A genome and the einkorn copy was 99.9%, decreasing to 99.0% when intronic sequences were compared (Additional file 6). The equivalent similarities between the B genome and the T. speltoides copies were 97.9% and 95.1% (Additional file 6). and that between the D genome and the goatgrass exonic sequence was 98.0%. Thus, polyploidization has left exonic sequence largely unchanged, while intronic sequence has been able to drift somewhat, as has been repeatedly observed in such comparisons [22, 34]. The origin of bread wheat has involved two separate allopolyploidization events, the first between einkorn and a species closely related to T. speltoides, and the second between the resulting wild tetraploid T. dicoccoides and goatgrass . The phylogeny predicted from the sequences of the SPSII homoeologues and wild species orthologues is fully consistent with this settled origin (Figure 4) and supports the previous published studies [49, 50].
Phylogenetic relationship of TaSPSII between different grasses
The phylogenetic analysis showed that the purple false brome SPSII sequence is more similar to those of wheat and barley than to that of rice, not only in terms of the sequence itself, but in terms of the gene's overall length. The purple false brome SPSII fragment extending from exon 8 to the 3'end of the gene (a region available for nine of the ten sequences) was 518 bp longer than the longest wheat or barley sequence, while this difference was 3,199 bp in rice and 3,663 bp in sorghum. A number of phylogenetic studies have shown that the taxonomic separation between the Triticeae and purple false brome is much less than between them and rice [51–53], which is the reason why purple false brome has been so readily adopted as a model genome for the Triticeae [53–55].
A study of the expression of sucrose biosynthesis gene family members has revealed that SPS genes are highly regulated both spatially and temporally. We have been able to locate the SPSII homoeoloci to the short arms of their respective group 3 chromosomes, and have shown that their divergence is consistent with the known evolutionary relationships between the component genomes of hexaploid wheat. A phylogenetic study confirmed that purple false brome is a more appropriate model for the Triticeae than rice. We have provided evidence that the syntenic relationship between rice chromosome 8 and wheat group 3 chromosomes is perturbed in the region where the SPSII genes map. We have shown that the transcript of the B genome copy of SPSII exists as two alternative forms, one with, and the other without an intron. The B genome copy contributed least to the global expression of SPSII expression. The development of homoeologue-specific markers for, and SNPs within the wheat SPSII genes will enable the definition of any associations between SPSII genotype and plant phenotype in terms of vigour, biomass productivity or source-sink relationships.
Plant materials and DNA isolation
Seed stocks of 27 bread wheat accessions and representative progenitors species were obtained from the gene bank of the Leibniz Institute of Plant Genetics and Crop Plant Research, Gatersleben (Table 3). A set of 21 nulli-tetrasomic lines in cv. 'chinese spring' originally obtained from late Dr. E.R. Sears, Coloumbia, MO, was used for chromosome arm assignment of PCR products. Total genomic DNA was extracted from 4-6 g young leaf material as described earlier  to be used as a template for PCR amplification.
Details of wheat accessions used for sequencing
Country of origin
Thatcher (TRI 1308)
Maris Huntsman (TRI 10287)
Chinese spring (TRI 23666)
W7984 (mapping parent)
- (synthetic line)
Triticum urartu (AA) (TRI 17140)
Aegilops speltoides (SS BB) (AE 413 )
Aegilops tauschii (DD) (TRI 145/96)
Triticum dicoccoides (AABB) (TRI 18495)
Triticum durum (AABB) (TRI 19197)
Expression analysis of SPS and SPP members
ESTs/cDNAs sequences for each of the five SPS gene families and three different SPP gene families in wheat were blasted against the probesets available on the wheat 61 K microarray platform on which corresponding wheat SPSI to SPSV and SPP1 to SPP3 cDNA gene sequences were identified with a significant e value of 10-10 or higher. To assess gene expression patterns of SPS and SPP gene family members during wheat plant development we downloaded CEL files of Affymetrix chip from the publicly available reference experiments (Series GSE12508; GSE9767; GSE6027) covering various tissues during germination (coleoptile, root, embryo), seedling stage (crown, leaf), stem at anthesis, reproductive tissues (immature inflorescence, floral bracts, pistil, mature anthers), various developmental stages of microsporogenesis and seed development. The transcriptome data of 61 K genes was subjected to RMA normalization with linear model using limma package. The log2 expression values were derived, differential expression is calculated with P value correction mode of Benzamini Hochberg and heat maps were generated for the differentially expressed SPS and SPP gene family members covering the complete wheat plant ontogeny using Genesis software .
The designing of SPSII family specific primers and subsequent PCR and cloning
The cDNA sequences of each SPS gene was retrieved from GenBank (accessions AF347064, AF347065, AF347066, AF347067, AF347068, AF347069, AF354298, AF534907 and AY4257109) and aligned using the Multalin  program. The SPSII gene family is represented by AF347064, AF347065, AF347066 and AF354298. Stretches of sequences limited to a particular gene family formed the basis for designing family specific primers, and in particular for SPSII (Table 2). To test the specificity of these family specific primers, the sequence of the resulting amplicons was used as a BLAST query against the GenBank sequence database http://blast.ncbi.nlm.nih.gov/Blast, applying high stringent conditions. Template DNA for the subsequent PCR was represented by cvs. 'Chinese Spring', 'M6' (a synthetic hexaploid wheat) and 'Opata'. Each 25 μl PCR comprised 50-100 ng template, 2.5 μl of buffer containing 1.5 mM MgCl2, 0.2 mM dNTP, 10 μM of each primer and 1 U Taq polymerase. The amplification conditions consisted of an initial denaturation step (94°C/4 min), followed by four cycles of 94°C/1 min, Ta/1 min, 72°C/2 min (where Ta began at 62°C and was reduced by 1°C per cycle), then 30 cycles of 94°C/1 min, 58°C/1 min, 72°C/2 min and completed by a final extension step of 72°C/5 min. This profile is referred to as "PCR profile 1". Amplicons were cloned into the pDrive cloning vector (QIAGEN, Hilden Germany). Plasmids were recovered from 20 clones from each of the three donor varieties using the QIAGEN plasmid purification kit (QIAGEN, Hilden Germany), and these were subjected to sequencing from both ends using a BigDye terminator v3.1 ready reaction cycle sequencing kit in combination with an ABI3730x1 sequencer (Applied Biosystems).
Phred-Phrap-Consed analysis of sequenced clones
Phred-Phrap-Consed analysis http://www.phrap.org/phredphrapconsed.html was applied to ensure the quality of the sequences taken forward [59, 60]. SNPs were detected by polyPhred  and visual inspection of Phrap assemblies. Two high quality contigs, assumed to represent two of the predicted three wheat homoeologues, were aligned using Sequencher software (Gene Codes Corporation, Ann Arbor, Michigan), and their consensus sequence formed a megablast [62, 63] query against the GenBank database. The outcome of this search was a confirmation that these two sequences were homologous to SPSII and not to any of the other SPS genes. Primer walking was employed to extend the sequence in both directions beyond the cloned region, and the sequence thereby acquired was subjected to megablast to ensure that it was SPSII specific.
Gene annotation, designing of homoeologue-specific primers, sequencing
A provisional gene structure was obtained by aligning the cDNA sequences with the rice genome sequence, and this was later confirmed by direct alignment with wheat genomic sequence. These alignments was performed using Genseqer http://www.plantgdb.org/PlantGDB-cgi/GeneSeqer/PlantGDBgs.cgi software, applying the the rice specific splice site model and a high stringency . Exon anchored primers were designed to span one or more intron using Primer3 software . Homoeologue-specific primers targeted indels and/or base substitutions mainly in the intronic sequence (Table 2), and their specificity was tested by the amplification of DNA from nulli-tetrasomic template. Four primer pairs specific for the 3A copy of SPSII were obtained, and six for the 3B copy (Table 2, Figure 2). However, no 3D copy primers were identified in this way. To obtain an assay for the 3D copy, primers designed to amplify an exonic fragment from bread wheat were applied to the goatgrass accession TRI 145/96, and a region of this amplicon which was specific to goatgrass (839-880 bp, see Additional file 3) was then targeted for successful primer design (Table 2, Figure 2). The A and B genome copies were amplified using PCR profile 1, while this was modified for the D genome copy to an initial denaturation step (94°C/4 min), followed by 35 cycles of 94°C/1 min, 65°C/1 min, 72°C/2 min and completed by a final extension step of 72°C/5 min (PCR profile 2). The primer sets were applied to obtain orthologous sequence from einkorn, T. speltoides, goatgrass and barley. For einkorn and T. speltoides it was possible to derive the sequence from intron 7 to the 3'end of the gene in this way, and from goatgrass and barley from exon 6 to the 3'end of the gene. For rice, sorghum and false purple brome, sequence information was retrieved from, respectively, TIGR http://rice.plantbiology.msu.edu/, Phytozome http://www.phytozome.net and BrachyBase http://www.brachypodium.org/. All sequences are shown in Additional file 1, and have been deposited in GenBank as accession numbers GU797178, GU797179, GU797180, GU797181, GU797182, GU797183 and GU929217.
Mapping of homeologue specific amplicons
Homeologue specific markers showing polymorphism between parents of International Triticeae Map ping Initiative (ITMI) population were used for genetic mapping. For genetic mapping 112 recombinant inbred (RI) lines from the ITMI population were used. This population was derived from the cross between W-7984, an amphi hexaploid wheat and Mexican wheat variety Opata 85 from CIMMYT (Centro Internacional de Me joramiento de Maizy Trigo) . The SPSII markers were integrated into the framework map composed of RFLP markers using software package MAPMAKER/Exp version 3.0b . Recombination fractions were converted to centimorgans with the Kosambi mapping function .
SNP identification, synteny study and cluster analysis
Primer pairs which produced a locus-specific amplicon were used amplify template DNA from 27 wheat cultivars, using the appropriate PCR profiles described above. The PCR products were purified with the MinElute™ UF PCR purification kit (QIAGEN, Hilden Germany) following the manufacturer's instructions, and subjected to sequencing from both ends. SNPs were identified by Sequencher v4.06 software (Gene Codes Cooperation, Ann Arbor, MI, USA). Position of the various SNPs shown in Additional file 3. Sequence alignment and phylogeny construction were performed using the MegAlign module of DNASTAR Lasergene® v8.0. software. The resulting sequence alignments and phylogenetic trees are shown in Additional files 2, 3, 4, 5 and Figure 4. Phylogenetic trees were subjected to bootstrap analysis, based on 10,000 replicates.
RNA isolation, cDNA synthesis and qRT-PCR
Total RNA was extracted from cv. 'Prinz' seedlings and developing caryopses at 4, 8 and 12 days after anthesis, and from seedlings of cv. 'Chinese Spring', einkorn, goatgrass and T.durum using the TRIzol reagent (Invitrogen GmbH, Karlsruhe, Germany) and RNAeasy columns (Qiagen, Hilden, Germany), as prescribed in the manufacturers' protocols, followed by in-column DNAse digestion. RNA concentration was measured using a NanoDrop photometer (Peqlab), following the manufacturer's instructions. A 2 μg sample of RNA was used as template for cDNA synthesis, using the SuperScript™III kit (Invitrogen GmbH), according to the manufacturer's instructions. cDNA yield was tested by quantitative PCR employing a serine/threonine protein phosphatase gene (TaPP2A) as reference. The cDNA was diluted to so that the concentration of all samples lay between ± 1 CT. A qRT-PCR assay of the A and D genome SPSII homoeologue targeted exonic regions, following recommended guidelines , while the B genome homoeologue targeted the unspliced intron. PCRs were performed in optical 384-well plates using the ABI PRISM® 7900 HT Sequence Detection System (Applied Biosystems, Foster City, CA, USA) and SYBR® Green to monitor dsDNA synthesis. Each 10 μl reaction contained 5 μl 2 × Power SYBR® Green Master mix reagent (Applied Biosystems), 1 μl cDNA and 200 nm of each homoeologue-specific primer. The PCR profile was 50°C/2 min and 95°C/10 min, followed by 45 cycles of 95°C/15 s, 60°C/60 s. Melting curves were recorded after cycle 45 by heating from 60°C to 95°C at 1.9°C min-1. Data were analysed using SDS2.2.1 software (Applied Biosystems). To generate a baseline-subtracted plot of the logarithmic increase in fluorescence signal (ΔRn) versus cycle number, the baseline data were collected between cycles 3 and 15. To obtain CT values, the amplification plots were analysed at an Rn threshold of 0.2. To compare data from different cDNA samples, the SPSII CT values were normalized against the expression of TaPP2A. PCR efficiency (E) was estimated from the exponential phase of each amplification, following . Since these lay between 1.8 and 2.0 for all samples, the expression level was calculated as 2-ΔCT, where ΔCT represented the difference between the CT of the SPSSII homoeologue and that of TaPP2A.
The financial assistance received from the Alexander von Humboldt foundation for carrying out this study is gratefully acknowledged. The research work was carried under Humboldt Research Fellowship awarded to Shailendra Sharma. Authors also thanks to Dr. Christof Pietsch for Phred-Phrap-Consed analysis.
Leibniz Institute of Plant Genetics and Crop Plant Research (IPK)
Sardar Vallabh Bhai Patel University of Agriculture and Technology
Iwate Biotechnology Research Center
Plant Breeding Institute, Christian-Albrechts University of Kiel
Department of Agronomy and Plant Breeding, College of Agriculture, Isfahan University of Technology
Department of Plant Pathology, Throckmorton Plant Sciences Center, Kansas State University
Lunn JE, Furbank RT: Sucrose biosynthesis in C4 plants.New Phytol 1999, 143:221–237.View Article
Lunn JE, MacRae E: New complexities in the synthesis of sucrose.Curr Opin Plant Biol 2003, 6:208–214.PubMedView Article
Strand A, Foyer CH, Gustafsson P, Gardestrom P, Hurry V: Altering flux through the sucrose biosynthesis pathway in transgenicArabidopsis thalianamodifies photosynthetic acclimation at low temperatures and the development of freezing tolerance.Plant Cell Environ 2003, 26:523–535.View Article
Yang J, Zhang J, Wang Z, Zhu Q: Activities of starch hydrolytic enzymes and sucrose-phosphate synthase in the stems of rice subjected to water stress during grain filling.J Exp Bot 2001, 52:2169–2179.PubMed
Castleden CK, Aoki N, Gillespie VJ, MacRae EA, Quick WP, Buchner P, Foyer CH, Furbank RT, Lunn JE: Evolution and function of the sucrose-phosphate synthase gene families in wheat and other grasses.Plant Physiol 2004, 135:1753–1764.PubMedView Article
Zhu YJ, Komor E, Moore PH: Sucrose accumulation in the sugarcane stem is regulated by the difference between the activities of soluble acid invertase and sucrose phosphate synthase.Plant Physiol 1997, 115:609–616.PubMedView Article
Grof CPL, Albertson PL, Bursle J, Perroux JM, Bonnett GD, Manners JM: Sucrose-phosphate synthase, a biochemical marker of high sucrose accumulation in sugarcane.Crop Sci 2007, 47:1530–1539.View Article
Mc Intyre CL, Jackson M, Cordeiro GM, Amouyal O, Hermann S, Aitken KS, Eliott F, Henry RJ, Casu RE, Bonnett GD: The identification and characterisation of alleles of sucrose phosphate synthase gene family III in sugarcane.Mol Breed 2006, 18:39–50.View Article
Causse M, Rocher JP, Henry AM, Charcosset A, Prioul JL, de Vienne D: Genetic dissection of the relationship between carbon metabolism and early growth in maize, with emphasis on key-enzyme loci.Mol Breed 1995, 1:259–272.View Article
Causse M, Rocher JP, Pelleschi S, Barriere Y, de Vienne D, Prioul JL: Sucrose-phosphate synthase: an enzyme with heterotic activity correlated with maize growth.Crop Sci 1995, 35:995–1001.View Article
Sarquis JI, Gonzalez H, Sánchez de Jimenez E, Dunlap JR: Physiological traits associated with mass selection for improved yield in a maize population.Field Crops Res 1998, 56:239–246.View Article
Prioul JL, Pelleschi S, Séne M, Thévenot C, Causse M, de Vienne D, Leonardi A: From QTLs for enzyme activity to candidate genes in maize.J Exp Bot 1999, 50:1281–1288.View Article
Bertin P, Gallais A: Genetic variation for nitrogen use efficiency in a set of recombinant inbred lines. II-QTL detection and coincidence.Maydica 2001, 46:53–68.
Ishimaru K, Ono K, Kashiwagi T: Identification of a new gene controlling plant height in rice using the candidate-gene approach.Planta 2003, 218:388–395.PubMedView Article
Coleman HD, Beamish L, Reid A, Park JY, Mansfield SD: Altered sucrose metabolism impacts plant biomass production and flower development.Transgenic Res 2009. DOI 10.1007/s11248009–9309–5.
Park JY, Canam T, Kang KY, Ellis DD, Mansfield SD: Over-expression of an arabidopsis family A sucrose phosphate synthase (SPS) gene alters plant growth and fibre development.Transgenic Res 2008, 17:181–192.PubMedView Article
Ishimaru K, Naoki H, Takayuki K, Yuka M, Kkiyoshi N, Kiyomi O, Ryu O: Overexpression of a MaizeSPSGene Improves Yield Characters of Potato under Field Conditions.Plant Prod Sci 2008, 11:104–107.View Article
Worrell AC, Bruneau JM, Summerfelt K, Boersig M, Voelker TA: Expression of a maize sucrose phosphate synthase in tomato alters leaf carbohydrate partitioning.Plant Cell 1991, 3:1121–1130.PubMedView Article
Lunn JE: Sucrose-phosphatase gene families in plants.Gene 2003, 303:187–196.PubMedView Article
Gale MD, Devos KM: Comparative genetics in the grasses.Proc Natl Acad Sci USA 1998, 95:1971–1974.PubMedView Article
Blake NK, Sherman JD, Dvorak J, Talbert LE: Genome-specific primer sets for starch biosynthesis genes in wheat.Theor Appl Genet 2004, 109:1295–1302.PubMedView Article
Caldwell KS, Dvorak J, Lagudah ES, Akhunov E, Luo MC, Wolters P, Powell W: Sequence polymorphism in polyploid wheat and their D-genome diploid ancestor.Genetics 2004, 167:941–947.PubMedView Article
Ravel C, Nagy IJ, Pierre M, Sourdille P, Mireille D, Balfourier F, Caroline P, Sandra G, Praud S, Charmet G: Single nucleotide polymorphism, genetic mapping, and expression of genes coding for the DOF wheat prolamin-box binding factor.Funct Integr Genomics2006:310–321.
Haudry A, Cenci A, Ravel C, Bataillon T, Brunel D, Poncet C, Hochu I, Poirier S, Santoni S, Glemin S, David J: Grinding up wheat: a massive loss of nucleotide diversity since domestication.Mol Biol Evol 2007, 24:1506–1517.PubMedView Article
Khlestkina EK, Röder MS, Salina EA: Relationship between homoeologous regulatory and structural genes in allopolyploid genome-A case study in bread wheat.BMC Plant Biol 2008, 8:88.PubMedView Article
Huang XQ, Röder MS: Development of SNP assays for genotyping thepuroindolineb gene for grain hardness in wheat using pyrosequencing.J Agric Food Chem 2005, 53:2070–2075.PubMedView Article
Whitt SR, Wilson LM, Tenaillon MI, Gaut BS, Buckler ES: Genetic diversity and selection in the maize starch pathway.Proc Natl Acad Sci USA 2002, 99:12959–12962.PubMedView Article
Ravel C, Praud S, Murigneux A, Canaguier A, Sapet F, Samson D, Balfourier F, Dufour P, Chalhoub B, Brunel D, Beckert M, Charmet G: Single-nucleotide polymorphism frequency in a set of selected lines of bread wheat (Triticum aestivumL.).Genome 2006, 49:1131–1139.PubMedView Article
Somers DJ, Kirkpatrick R, Moniwa M, Walsh A: Mining single-nucleotide polymorphisms from hexaploid wheat ESTs.Genome 2003, 46:431–7.PubMedView Article
Bryan GJ, Stephenson P, Collins A, Kirby J, Smith JB, Gale MD: Low levels of DNA sequence variation among adapted genotypes of hexaploid wheat.Theor Appl Genet 1999, 99:192–198.View Article
Coleman CE, Clore AM, Ranch JP, Higgins R, Lopes MA, Larkins BA: Expression of a mutant a-zein creates the floury2 phenotype in transgenic maize.Proc Natl Acad Sci USA 1997, 94:7094–7097.PubMedView Article
Aukerman MJ, Schmidt RJ, Burr B, Burr FA: An arginine to lysine substitution in the bZIP domain of an opaque-Z mutant in maize abolishes specific DNA binding.Genes Dev 1991, 5:310–320.PubMedView Article
Haga H, Yamada R, Ohnishi Y, Nakamura Y, Tanaka T: Gene based SNP discovery as part of the Japanese Millenium Genome Project: identification of 190,562 genetic variations in the human genome.J Hum Genet 2002, 47:605–610.PubMedView Article
Nomura T, Ishihara A, Yanagita RC, Endo TR, Iwamura H: Three genomes differentially contribute to the biosynthesis of benzoxazinones in hexaploid wheat.Proc Natl Acad Sci USA 2005, 102:16490–16495.PubMedView Article
Kurata N, Moore G, Nagamura Y, Foote T, Yano M, Minobe Y, Gale M: Conservation of genome structure between rice and wheat.Nat Biotech 1994, 12:276–278.View Article
Devos KM, Gale MD: Genome relationships: the grass model in current research.Plant Cell 2000, 12:637–646.PubMedView Article
Munkvold JD, Greene RA, Bermudez-Kandianis CE, La Rota CM, Edwards H, Sorrells SF, Dake T, Benscher D, Kantety R, Linkiewicz AM, Dubcovsky J, Akhunov ED, Dvorák J, Miftahudin , Gustafson JP, Pathan MS, Nguyen HT, Matthews DE, Chao S, Lazo GR, Hummel DD, Anderson OD, Anderson JA, Gonzalez-Hernandez JL, Peng JH, Lapitan N, Qi LL, Echalier B, Gill BS, Hossain KG, Kalavacharla V, Kianian SF, Sandhu D, Erayman M, Gill KS, McGuire PE, Qualset CO, Sorrells ME: Group 3 chromosome bin maps of wheat and their relationship to rice chromosome 1.Genetics 2004, 168:639–50.PubMedView Article
Sorrells ME, La Rota M, Bermudez-Kandianis CE, Greene RA, Kantety R, Munkvold JD, Miftahudin , Mahmoud A, Ma X, Gustafson PJ, Qi LL, Echalier B, Gill BS, Matthews DE, Lazo GR, Chao S, Anderson OD, Edwards H, Linkiewicz AM, Dubcovsky J, Akhunov ED, Dvorak J, Zhang D, Nguyen HT, Peng J, Lapitan NL, Gonzalez-Hernandez JL, Anderson JA, Hossain K, Kalavacharla V, Kianian SF, Choi DW, Close TJ, Dilbirligi M, Gill KS, Steber C, Walker-Simmons MK, McGuire PE, Qualset CO: Comparative DNA sequence analysis of wheat and rice genomes.Genome Res 2003, 13:1818–1827.PubMed
Li W, Gill BS: The colinearity of the Sh2/A1 orthologous region in rice,Sorghumand maize is interrupted and accompanied by genome expansion in the Triticeae.Genetics 2002, 160:1153–1162.PubMed
Isidore E, Scherrer B, Chalhoub B, Feuillet C, Keller B: Ancient haplotypes resulting from extensive molecular rearrangements in the wheat A genome have been maintained in species of three different ploidy levels.Genome Res 2005, 15:526–536.PubMedView Article
Li Z, Huang B, Rampling L, Wang J, Yu J, Morell M, Rahman S: Detailed comparison between the wheat chromosome group 7 short arms and the rice chromosome arms 6S and 8L with special reference to genes involved in starch biosynthesis.Funct Integr Genomics 2004, 4:231–240.PubMed
Paux E, Sourdille P, Salse J, Saintenac C, Choulet F, Leroy P, Korol A, Michalak M, Kianian S, Spielmeyer W, Lagudah E, Somers D, Kilian A, Alaux M, Vautrin S, Bergès H, Eversole K, Appels R, Safar J, Simkova H, Dolezel J, Bernard M, Feuillet C: A physical map of the 1-gigabase bread wheat chromosome 3B.Science 2008, 322:101–104.PubMedView Article
Choulet F, Paux E, Salse J, Leroy P, Ghislaine M, Sylvie S, Valerie B, Feullet C: Sequencing, annotation and charcterization of 17 MB of chromosome 3B contigs provide novel insights into the wheat genome organization and evolution [abstract]. In 19th International Triticae Mapping Initiative-3rd COST Tritigen 2009. Clermont-Ferrand, France; 2009:58.
Paterson AH, Bowers JE, Chapman BA: Ancient polyploidization predating divergence of the cereals, and its consequences for comparative genomics.Proc Natl Acad Sci USA 2004, 101:9903–9908.PubMedView Article
Lu HJ, Faris JD: Macro- and micro-colinearity between the genomic region of wheat chromosome 5B containing the Tsn1 gene and the rice genome.Funct Int Genomics 2006, 6:90–103.View Article
Singh NK, Dalal V, Batra K, Singh BK, Chitra G, Singh A, Ghazi IA, Yadav M, Pandit A, Dixit R, Singh PK, Singh H, Koundal KR, Gaikwad K, Mohapatra T, Sharma TR: Single-copy genes define a conserved order between rice and wheat for understanding differences caused by duplication, deletion, and transposition of genes.Funct Integr Genomics 2007, 7:17–35.PubMedView Article
Feldman M, Liu B, Segal G, Abbo S, Levy AA, Vega JM: Rapid elimination of low-copy DNA sequences in polyploidy wheat: A possible mechanism for differentiation of homoeologous chromosomes.Genetics 1997, 147:1381–1387.PubMed
Chalupska D, Lee HY, Faris JD, Evrard A, Chalhoub B, Haselkorn R, Gornicki P: Acc homoeoloci and the evolution of wheat genomes.Proc Natl Acad Sci USA 2008, 105:9691–9696.PubMedView Article
Dvorak J, Zhang HB: Variation in repeated nucleotide sequences sheds light on the origin of the wheat B and G genomes.Proc Natl Acad Sci USA 1990, 87:9640–9644.PubMedView Article
Vogel JP, Gu YQ, Twigg P, Lazo GR, Laudencia-Chingcuanco D, Hayden DM, Donze TJ, Vivian LA, Stamova B, Coleman-Derr D: EST sequencing and phylogenetic analysis of the model grassBrachypodiumdistachyon.Theor Appl Genet 2006, 113:186–195.PubMedView Article
Bossolini E, Wicker T, Knobel PA, Keller B: Comparison of orthologous loci from small grass genomesBrachypodiumand rice: implications for wheat genomics and grass genome annotation.Plant J 2007, 49:704–717.PubMedView Article
Huo N, Lazo GR, Vogel JP, You FM, Ma Y, Hayden DM, Coleman-Derr D, Hill TA, Dvorak J, Anderson OD, Luo MC, Gu YQ: The nuclear genome of Brachypodium distachyon: analysis of BAC end sequences.Funct Integr Genomics 2007, 8:135–47.PubMedView Article
Huo N, Vogel JP, Lazo GR, You FM, Ma Y, McMahon S, Dvorak J, Anderson OD, Luo MC, Gu YQ: Structural characterization ofBrachypodiumgenome and its syntenic relationship with rice and wheat.Plant Mol Biol 2009, 70:47–61.PubMedView Article
Kumar S, Mohan A, Balyan HS, Gupta PK: Orthology between genomes ofBrachypodium, wheat and rice.BMC Res Notes 2009, 2:93.PubMedView Article
Plaschke J, Ganal MW, Röder MS: Detection of genetic diversity in closely related bread wheat using microsatellite markers.Theor Appl Genet 1995, 91:1001–1007.View Article
Corpet F: Multiple sequence alignment with hierarchical clustering.Nucleic Acids Res 1988, 16:10881–10890.PubMedView Article
Ewing B, Green P: Base-calling of automated sequencer traces using Phred. II. Error probabilities.Genome Res 1998, 8:186–194.PubMed
Gordon D, Abajian C, Green P: Consed: a graphical tool for sequence finishing.Genome Res 1998, 8:195–202.PubMed
Nickerson D, Tobe V, Taylor S: PolyPhred: automating the detection and genotyping of single nucleotide substitutions using fluorescence-based resequencing.Nucl Acids Res 1997, 25:2745–2751.PubMedView Article
Altschul SF, Madden TL, Schäfer AA, Zhang JZ, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST; a new generation of protein database search programs.Nucleic Acids Res 1997, 25:3389–3402.PubMedView Article
Schlueter SD, Dong Q, Brendel V: GeneSeqer@PlantGDB: gene structure prediction in plant genomes.Nucleic Acids Res 2003, 31:3597–3600.PubMedView Article
Rozen S, Skaletsky H: Primer3 on the www for general users and for biologist programmers. In Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. Totowa, NJ: Humana Press; 2000:365–386.
Van Deynze AE, Dubcovsky J, Gill KS, Nelson JC, Sorrells ME, Dvorak J, Gill BS, Lagudah ES, McCouch SR, Appels R: Molecular genetic maps for group 1 chromosomes of Triticeae species and their relation to chromosomes in rice and oat.Genome 1995, 38:45–59.
Lander ES, Green P, Abrahamson J, Barlow A, Daly MJ, Lincoln SE, Newberry L: Mapmaker: an interactive computer package for constructing primary genetic linkage maps of experimental and natural populations.Genomics 1987, 1:174–181.PubMedView Article
Kosambi D: The calculation of map distances from recombination values.Ann Eugen 1944, 12:172–175.
Udvardi MK, Czechowski T, Scheible WR: Eleven Golden Rules of Quantitative RT-PCR.Plant Cell 2008, 20:1736–1737.PubMedView Article
Ramakers C, Ruijter JM, Deprez RH, Moorman AF: Assumption-free analysis of quantitative real-time PCR data.Neurosci Lett 2003, 339:62–66.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.