Increasing the amylose content of durum wheat through silencing of the SBEIIagenes
© Sestili et al; licensee BioMed Central Ltd. 2010
Received: 17 December 2009
Accepted: 14 July 2010
Published: 14 July 2010
High amylose starch has attracted particular interest because of its correlation with the amount of Resistant Starch (RS) in food. RS plays a role similar to fibre with beneficial effects for human health, providing protection from several diseases such as colon cancer, diabetes, obesity, osteoporosis and cardiovascular diseases. Amylose content can be modified by a targeted manipulation of the starch biosynthetic pathway. In particular, the inactivation of the enzymes involved in amylopectin synthesis can lead to the increase of amylose content. In this work, genes encoding starch branching enzymes of class II (SBEIIa) were silenced using the RNA interference (RNAi) technique in two cultivars of durum wheat, using two different methods of transformation (biolistic and Agrobacterium). Expression of RNAi transcripts was targeted to the seed endosperm using a tissue-specific promoter.
Amylose content was markedly increased in the durum wheat transgenic lines exhibiting SBEIIa gene silencing. Moreover the starch granules in these lines were deformed, possessing an irregular and deflated shape and being smaller than those present in the untransformed controls. Two novel granule bound proteins, identified by SDS-PAGE in SBEIIa RNAi lines, were investigated by mass spectrometry and shown to have strong homologies to the waxy proteins. RVA analysis showed new pasting properties associated with high amylose lines in comparison with untransformed controls. Finally, pleiotropic effects on other starch genes were found by semi-quantitative and Real-Time reverse transcription-polymerase chain reaction (RT-PCR).
We have found that the silencing of SBEIIa genes in durum wheat causes obvious alterations in granule morphology and starch composition, leading to high amylose wheat. Results obtained with two different methods of transformation and in two durum wheat cultivars were comparable.
Cereal grains contain a good balance of proteins, fats, carbohydrate, vitamins and minerals required for human growth and health. Unlike other cereals, wheat is rarely consumed in an unprocessed form but prepared into a wide range of end products. Common wheat (Triticum aestivum L.) is used in the preparation of bread, noodles, biscuits, and cakes. Durum wheat (T. turgidum L. var. durum) is used mainly for pasta production but also in an array of other regional foods in Italy, North Africa and West Asia (bread, cous cous, burghoul etc). The processing and end-use quality of wheat-based products depends on different factors such as protein content and composition, grain hardness and starch composition. Starch, the most important polysaccharide in human diet and is the major component of the wheat kernel, representing more than 70% of its dry weight. As well as its importance in the food industry, starch is also used as a raw material for the production of non-food products in the paper, plastic, adhesive, textile, medical and pharmaceutical industries .
Reserve starch is accumulated in the amyloplast organelles and is composed of two different glucosidic polymers, amylose and amylopectin. The main differences between these polymers are the degree of polymerization and the number of side branches. Amylose is a linear chain of D-glucose molecules with a low degree of polymerization (< 104 units), whereas amylopectin shows a higher degree of polymerization (105-106 units) and which has important implications for function. Amylopectin is the major constituent of starch in wheat endosperm and comprises about 70-80%; with amylose constituting the remaining 20-30%. Amylose and amylopectin are synthesized by two different pathways having a common substrate (ADP-glucose). A granule bound starch synthase (GBSSI) is involved in amylose synthesis, whereas amylopectin is produced by the concerted action of starch synthases (SSI, SSII, SSIII), starch branching enzymes (SBEI, SBEIIa and SBEIIb) and starch debranching enzymes of isoamylase- and limit dextrinase-type (ISA and LD) [2, 3].
SBEs are transglycosylase enzymes that catalyze the formation of α-1,6 linkages within the polymer by cleaving an internal alpha-1,4 linkage. In monocots, three starch branching isoforms are found: SBEI, SBEIIa and SBEIIb. In maize, rice and pea, suppression of SBEIIb leads to amylose-extender (ae) phenotype, with a very high amylose content (>50%) [4–6], in contrast suppression of SBEI or SBEIIa has no effect on the amount of amylose [7–9]. In wheat SBEIIa and SBEIIb genes have been characterized and found to be located on the long arm of the homoeologous group 2 chromosomes [10–12]. Regina et al.  demonstrated that wheat SBEIIa gene is syntenic to the corresponding gene in other cereals, in contrast the SBEIIb gene is not in a syntenic position. In wheat, SBEIIa is the predominant isoform present in the soluble phase of the endosperm , whereas in maize and rice endosperm SBEIIb is the predominant isoform involved in amylopectin biosynthesis [13, 14].
The role of SBEIIa and SBEIIb isoforms in bread wheat endosperm has been investigated by RNA interference technology . In contrast to other cereals, the silencing of SBEIIb genes has no effect on amylose content and starch granule shape; whereas silencing of SBEIIa genes results in a strong increase in amylose content (>70%) and granule deformation.
There is increasing interest in the manipulation of starch composition in wheat due to the recognition of its important role in food and non food applications and its uses in industry. In addition, the research is also focusing on the production of high amylose starch flours because derived foods have an increased amount of resistant starch which has been shown to have beneficial effects on human health. Resistant starch refers to the portion of starch that resists digestion in the stomach and small intestine and it can act as a substrate for microbial fermentation in the large intestine, the end-products being hydrogen, carbon dioxide, methane and short chain fatty acids (SCFA) [16, 17]. The nutritionists believe that the resistant starch has a role similar to dietary fibre inside the intestine, protecting against diseases as colon cancer, type II diabetes, obesity and osteoporosis [17–21].
A recent study demonstrated that high-amylose flours might be used to substitute for up to 50% normal wheat flour to make bread with both an acceptable bread quality and a significantly higher amount of resistant starch .
Furthermore, pasta produced with semolina containing a higher amylose content shows good cooking resistance and firmness, satisfying consumer preferences .
In the present study, a transgenic approach was used to increase the amylose content in durum wheat seeds. In particular, the effects of SBEIIa silencing were investigated in terms of amylose content, transcript accumulation and protein profile of the enzymes involved in starch biosynthesis.
Biolistic transformation and selection of transgenic plants
In total, 1954 immature embryos obtained from the Triticum durum cv. Svevo were co-transformed using pRDPT + SBEIIa(RNAi) and bar selectable marker plasmids. Forty eight T0 independent transgenic lines resistant to bialaphos and containing the RNAi cassette were obtained, with a co-transformation efficiency of 2.5%. The presence of the transgene was verified by polymerase chain reaction (PCR) analysis of genomic DNA from regenerated plants, by using two primer pairs, one specific for the promoter of the construct (PromDx5Fw/PromDx5R) and another for the bar gene (BarFw/BarR), to produce amplicons of 473 bp and 405 bp, respectively (data not shown). No significant differences in morphology and growth were observed between transgenic lines and either the corresponding "null" genotype that had lost the transgene by segregation (null segregant), or untransformed plants.
Agrobacterium-mediated transformation and selection of transgenic plants
In all, 1.759 immature embryos of the Triticum durum cv. Ofanto were transformed using pGUB-G + SBEIIa(RNAi) containing also the bar gene cassette in the T-DNA. Thirteen T0 transgenic lines were obtained with a transformation efficiency of 0.74%. The presence of the transgene was verified as described above. As with transgenic plants generated by biolistic transformation, no morphological differences were observed among the different lines or between transgenics and controls.
Molecular characterization of transgenic plants
The absence of SBEIIa transcript was investigated by semiquantitive RT-PCR and Real Time RT-PCR on total RNA, extracted from T2 immature kernels (18 DPA), using a primer pair specific for SBEIIa genes (SBEIIaFw/SBEIIaR). We focused on six transgenic lines, three derived from the biolistic method and three from Agrobacterium-mediated transformation. GADPH (glyceraldehyde 3-phosphate dehydrogenase) was used as housekeeping gene.
Peptides identified by mass-mass spectrometry on SGP-4 and SGP-5 proteins.
Pepti de N°
Actual mass (MH+)
Theor mass (MH+)
Determination of total starch, amylose and protein content
Seed weight, total starch, amylose and protein content of RNAi transgenic lines and control plants.
100 grain weight
Starch content %
Amylose content %
4.14 ± 0.14
59.80 ± 0.10
24.50 ± 0.70
18.36 ± 0.13
3.90 ± 0.12
50.98 ± 0.68
75.05 ± 0.36
18.95 ± 0.36
3.74 ± 0.10
44.91 ± 2.36
43.50 ± 2.52
16.16 ± 0.30
4.06 ± 0.22
45.17 ± 1.47
32.40 ± 0.39
17.38 ± 0.16
27.69 ± 1.49
56.44 ± 1.43
44.25 ± 0.46
30.87 ± 1.52
Starch granule analysis
RVA parameters of durum wheat cultivar Svevo and RNAi transgenic lines (MJ16-112 and MJ16-114).
Real Time RT-PCR on genes involved in starch biosynthesis
To understand if there were pleiotropic effects on the expression of other genes involved in starch biosynthesis, semiquantitative RT-PCR was performed on GBSSI, SSI, SSII, SSIII, SBEI, SBEIIb, ISA1 and LD genes. The expression levels of many starch biosynthetic genes were affected by the silencing of SBEIIa genes. Compared to the control lines, the transcripts of the genes encoding GBSSI, SSIII, LD and ISA1 were increased in transgenic plants, whereas significant changes were not found in SSI, SSII, SBEI and SBEIIb genes.
Transcript accumulation of GBSSI, SSIII, LD and ISA1 were investigated by quantitative Real-Time RT-PCR. As expected, all analyzed genes were up-regulated in comparison with untransformed plants but their level of up-regulation varied (Figure 2). The expression of GBSSI was strongly increased in transgenic seeds with a 12-28 fold increase compared to controls). The expression of SSIII, LD and ISA1 was up-regulated by 6-18, 1,5-3,5 and 4-6 fold respectively in Svevo transgenic plants compared to controls (Figure 2).
Different strategies can be used to modify amylose content in wheat starch. In the last twenty years the identification of null mutations in starch biosynthetic genes and the crossing of these materials to combine functional mutations have permitted the production of lines with low and high amylose content [24–26]. Unfortunately this strategy is not easy to apply for all starch biosynthetic proteins, due to the difficulty of obtaining the relevant mutants in polyploid species such as wheat.
Both transgenic and non-transgenic tools are available, and the latter, along with novel technologies, offers the possibility to manipulate starch composition and develop wheat varieties with new functionalities. Recently a reverse genetic strategy called TILLING (Target Induced Local Lesions in Genomes) has permitted the identification of induced mutations for genes involved in amylose  and amylopectin  biosynthesis. The TILLING strategy is an efficient, non-transgenic method to identify mutants in genes of interest, but the effect of silencing is often masked by the presence of other functional alleles in polyploid species. To obtain a null phenotype, it is usually necessary to generate double or triple mutants by crossing single null genotypes .
Recent studies have demonstrated that transgenic technology is a powerful research tool to modify amylose/amylopectin ratio. Starch composition can be manipulated through genetic engineering, either by overexpression or repression of starch biosynthetic genes. RNAi technology is a specific and powerful mean to silence a target gene. This technology is particularly useful in polyploid plants such as wheat because it increases the possibility to obtain transgenic plants with a novel phenotype.
In this study, an RNAi strategy has been followed to silence SBEIIa genes in order to increase the amylose content in two durum wheat cultivars, Svevo and Ofanto. Although two different methods have been used for genetic transformation (biolistic for cv Svevo and Agrobacterium for cv Ofanto), the silencing of SBEIIa genes produced similar effects. The Agrobacterium- or the biolistic-mediated genetic transformation showed various levels of silencing between independent transgenic wheat lines. The abundance of mRNAs of SBEIIa genes varied in the transgenic lines from nearly undetectable to wild type level confirming that the RNAi construct, designed to reduce the expression of SBEIIa genes, was effective in inducing gene silencing. The variation in the levels of gene silencing observed here is in agreement with previously reported results in wheat [30, 31] and other species, such for example tobacco , maize  and Arabidopsis .
The silencing of the SBEIIa gene caused marked modifications on amylose content, granule morphology and starch composition. An increase in amylose content was found in all transgenic lines but its percentage varied between lines, presumably due to differences in copy number and/or positional effects. The line MJ16-112 showed the exceptionally high value of 75% amylose content, while in the other lines it varied between 31-56%. A variable level of amylose content between RNAi lines was also found in sweetpotato  and potato .
No significant difference was observed in grain yield between transformed and control lines, although RNAi lines showed a lower total starch content. The reduction of total starch in RNAi lines could be due to the difficulty to estimate starch content in high amylose wheats. McClearly et al.  used different approaches to determine total starch content and found that in high amylose maize the starch content resulted always underestimated.
Morphological analysis of starch granules showed a marked difference between SBEIIa-silenced lines and control plants. Type-A granules appeared smaller and deflated, whereas type-B granules lost their normal spherical shape and became extended. Similar results were obtained by Regina et al.  who suppressed the activity of SBEIIa through RNAi technology in bread wheat.
RVA analysis, carried out on whole wheat flour, highlighted large alterations in chemical-physical properties of transgenic starch. RVA parameters, as mean trough, breakdown, setback and final viscosities of SBEIIa-silenced lines, were significantly lower than normal wheats. These results are consistent with those obtained on high amylose flour of wheat , corn  and barley . The decreased swelling and viscosity affects chemical-physical pasting properties of semolina, with subsequent changes in functional characteristics of derived food products. Although high amylose flour resulted in a smaller bread loaf volume , it has a positive impact in the pasta industry, reducing starch loss during the cooking and increasing pasta firmness. In addition, high amylose foods might have beneficial effects on human diet because of an increased amount of resistant starch, which has a role similar to dietary fibre inside the intestine, protecting against diseases as colon cancer, type II diabetes, obesity and osteoporosis [17, 18].
SDS-PAGE analysis of starch granule-bound proteins allowed the identification of two newly accumulated proteins localized between the SGP-2 and SGP-3 in the transgenic lines. These novel proteins are present in all transgenic plants, suggesting that they are a consistent result of SBEIIa silencing rather than an unintended effect of random transgene insertions. The novel proteins show a high similarity to the waxy isoforms. Since the molecular mass of the novel proteins is larger (85 and 78 KDa) than that of waxy proteins (58-60 KDa), one possibility is that the novel proteins are encoded by genes that are poorly or not expressed in wild type cultivar. An alternative explanation is that a mechanism of post-translation modification has somehow been triggered in these transgenic lines.
The staining intensity of the other starch granule proteins (SGP-140, SGP-145, SGP-1, SGP-3 and waxy) on SDS-PAGE is strongly increased in the silenced lines compared to the control. Real-Time RT-PCR analyses on genes involved in starch biosynthesis have confirmed possible interaction and feedback mechanisms between these genes. In particular, the absence of SBEIIa isoforms leads to an up-regulation of Wx, SSIII, LD and ISA1 genes.
Pleiotropic effects associated with up-regulation or suppression of specific genes present in the starch biosynthetic pathway have been reported previously. For example, Kosar-Hashemi et al.  observed a reduction of granule-bound branching enzymes and starch synthases in the starch granule of the SSII mutants and suggested that the lack of SSII would lead to the absence of the complex on starch granule. Similarly, the rice sugary-1 (null for isoamylase) and the maize dull-1 (not expressing SSIII) mutants had drastic changes in accumulation of pullulanase and SBEIIb isoforms, respectively [43, 44]. In a recent work, the down-regulation of a pullulanase-type DBE inhibitor produced pleiotropic effects on activities of soluble and granule-bound starch synthase, that were significantly reduced . Tetlow et al.  have discovered the formation of a heterocomplex comprising SBEI, SBEIIb, and starch phosphorylase in wheat. More recently, Tetlow et al.  have identified trimeric assemblages of SSI with SSIIa and either SBEIIa or SBEIIb. Hennen-Bierwagen et al. [48, 49] have demonstrated that SSs, SBEs and DBEs associate physically with each other in multisubunit complexes and the lack of one subunit can cause variation on the other isoforms. In maize the ae- mutation, in the gene coding for SBEIIb, was associated with the detection of four novel proteins not observed in wild-type granules. LC-MS/MS analysis showed that the four bands corresponded to SBEI, SBEIIa, SSIII and starch phosphorylase . Grimaud et al.  suggested that structural changes in amylopectin, due to the loss of SBEIIb, could offer a more efficient glucan substrate for binding of SBEI, SBEIIa, starch phosphorylase and SSIII. The pleiotropic effects observed in our SBEIIa-silenced lines represent additional evidence of a strict correlation between all starch enzymes and that the lack of one subunit can cause variation in the other isoforms.
The RNAi-mediated down-regulation of genes coding the SBEIIa isoforms produced remarkable changes in starch chemical-physical properties of durum wheat. No significant differences were identified in grain weight between transgenic plants and wild type controls. As well as a considerable increase in the amylose content, the silenced SBEIIa wheat lines also possessed altered starch granule structure and up-regulated transcriptional profiles of GBSSI, SSIII, LD and ISA1. The use of two different transformation methods, biolistic and Agrobacterium-mediated, and two different durum wheat cultivars, highlights that the effect of SBEIIa silencing on starch metabolism is genotype and protocol independent because similar effects have been detected in all the RNAi plants obtained.
The variation in amylose content and the different starch physical-chemical properties make these plant materials particularly valuable in a vast range of applications in both food and non food industries.
Durum wheat plants (cv Svevo and cv Ofanto) were vernalized at 4°C for one month and grown in a growth chamber at 20/18°C day/night temperature under a 16/8 h photoperiod.
Vectors for biolistic and Agrobacterium-mediated transformation
The first three exons, of the wheat SBEIIa gene were amplified by RT-PCR on total RNA extracted from immature seeds (18 dpa), while the third intron of the wheat SBEIIa gene was amplified from genomic DNA of Svevo. The primers used were: SalIXbaISBEII(cDNA)FW 5'-gtcgactctagaggcgggttgagtgagatctg; KpnIXhoISBEII(cDNA)R 5'-ggtaccctcgagtcggtagtcaagatggctcc; SalIKpnISBEII(INT)FW 5'-gtcgacggtaccgcagaaaatatacgagattg; XbaIXhoISBEII(INT)R 5'-acctctagactcgagccaccttcatgttggtcaatag.
cDNA synthesis was performed by using the QuantiTect Reverse Transcrition Kit (Qiagen) following the manufacturer instructions.
PCR reactions were carried out in 50 μl final volume using 50-100 ng of genomic DNA or 1 μl of cDNA, 2.5 units of FastStart High Fidelity PCR system (Roche Diagnostics), 1× Taq PCR buffer, 50 ng of each of the two primers and 100 μM of each deoxyribonucleotide. PCR conditions were: 1 cycle at 95°C for 2 min, 35 cycles at 95°C for 30 s, 60°C for 1 min, 72°C for 1 min, and a final step at 72°C for 5 min.
The intron was ligated in the vector pRDPT containing the promoter and terminator of the Dx5 high molecular weight glutenin subunit gene . Then taking advantage of restriction sites introduced on the primers, the sense and antisense sequences of the target region were inserted in the plasmid pRDPT separated each one by the intron. This construct was used for biolistic transformation.
The entire cassette (promoter-sense-intron-antisense-terminator) was extracted from pRDPT vector by digestion with SmaI and HpaI restriction enzymes and cloned into the SmaI site of the vector pG-UB carrying the bar gene that confers resistance to the bialaphos herbicide under the ubiquitin promoter. This construct was used for Agrobacterium-mediated transformation.
Biolistic and Agrobacterium transformation of durum wheat embryos
For the transformation experiments, the plasmid pAHC20 , carrying the bar gene, was co-bombarded with pRDPT + SBEIIa(RNAi) in a 1:3 molar ratio. Constructs were introduced into immature embryos of T. durum cv. Svevo excised from seeds in Zadoks stage 72, using a Model PDS-1000/He Biolistic particle delivery system (Bio-Rad, Hercules, CA, U.S.A.) and the protocols described in Okubara and associates . The protocol of Wu et al.  was used for transformation of durum wheat Ofanto by Agrobacterium. The presence of the construct in bialaphos-resistant plants and their progeny was verified by PCR on total DNA obtained from leaf sections of mature plants  with primer pairs specific for the promoter Dx5 and for the bar gene (PromDx5Fw catgcaggctaccttccac, PromDx5R cggtggactatcagtgaattg, BarFw catcgagacaagcacggtca and BarR gaaacccacgtcatgccagt). T2 non segregating progeny were used for molecular and biochemical analyses (additional file 1).
PCR reactions were carried out in 50 μl final volume using 50-100 ng of genomic DNA, 1× Red Taq ReadyMix PCR REACTION MIX (1.5 U Taq DNA Polymerase, 10 mM Tris-HCl, 50 mM KCl, 1.5 mM MgCl2, 0.001% gelatine, 0.2 mM dNTPs) and 0.5 μM of each of the two primers.
Amplification conditions included an initial denaturation step at 94°C for 3 min, followed by 35 cycles at 94°C for 1 min, 60°C for 1 min and 72°C for 1 min, followed by a final incubation at 72°C for 5 min.
Semiquantitative Reverse transcriptase-polymerase chain reaction
Total RNA was extracted from immature seeds (18 DPA) as reported in Laudencia-Chingcuanco et al.  with some modifications. The starting material was 0.1 g and all volumes of buffers and solutions were diluted 1 to 10. For reverse transcriptase-mediated PCR studies, cDNA was synthesised from 1 μg of total RNA using an oligo(dT) primer and Superscript Reverse Transcriptase III (Invitrogen). One of twentieth volume of each cDNA was used as a template for PCR amplification. PCR reactions were carried out in 30 μl final volume using 1 units of Ex-Taq (Takara), 1× buffer, 0,2 mM of each dNTPs, 0,5 μM of each primer. Amplification conditions included an initial denaturation step at 98°C, followed by 35 cycles at 98°C for 10 sec., 58°C for 1 min. and 72°C for 1 min, followed by a final extension at 72°C for 5 min. The following gene-specific primers were designed for GAPDH (GAPDHFw 5'-caacgctagctgcaccactaact and GAPDHR 5'-gactcctccttgatagcagcctt) , GBSSI (Wx-B1Fw 5'-cgaagcaacaaagccggaaag and Wx-B1R 5'-tcaccctctcgtactcgtccg), SSI (SSIFw 5'-agggtacagggtgggcgttct and SSIR 5'-gtagggttggtccacgaagg), SSII (SSIIFw 5'-ttcgaccccttcaaccactc and SSIIR 5'-acgtcctcgtagagcttggc), SSIII (5'-ggcgttggatgtgtatatgg and 5'-ggtgatgattccgacaatagg) , SBEI (SBEIFw 5'-aaacaaacttcaacaaccgcc and SBEIR 5'-acctcctgtagacgcctttt), SBEIIa (SBEIIaFw 5'-tgacgaatcttggaaaatgg and SBEIIaR 5'-ggcggcatttatcataactattg), SBEIIb (SBEIIbFw 5'-gtagatgcggtcgtttacttga and SBEIIbR 5'-ccagccaccttctgtttgtt), LD (LDFw 5'-tatgagtggaacagggattggt and LDR 5'-atttggtcagcgtaagtagcg) and ISA1 (ISA1Fw 5'-aagatgaaagacagggcgagat and ISA1R 5'-aatactaggatgaccgatgagtagct).
Real-Time RT-PCR (qRT-PCR)
One microlitre of the cDNA above described was used for real-time PCR in a 20 μL volume. For each sample three technical replicates were used for PCR amplification. The PCR reaction consisted of 10 μL of iQ™ SYBR Green Supermix 2X (BIO-RAD), which contained buffer, dNTPs and SYBR Green I. Concentrations of the forward and reverse oligodeoxynucleotide primers in the reaction were 500 nM for all the genes of interest.
qRT-PCR experiments were performed using the iCycler iQ (Bio-Rad Laboratories, Hercules, CA1, USA). Amplification conditions were as follows: initial 95°C for 15 min and 40 cycles of 95°C for 30 s, 60°C for 1 min and 72°C for 1 min each.
Relative expression analysis was determined by using the 2-ΔΔCT method ; Applied Biosystems User Bulletin No. 2-P/N 4303859). Calculation and statistical analyses were performed by Gene Expression Macro™ Version 1.1 (Bio-Rad Laboratories, Hercules, CA, USA). The efficiencies of target and housekeeping genes were determined by qRT-PCR on serial dilutions of RNA template over a 100-fold range , with similar results (data not shown). Amplified products were checked by gel electrophoresis and sequencing to verify primer specificity. Relative expression of each gene is reported as the number of fold increase of the transcript level at each time point, compared to the lowest transcript level.
As for semiquantitative RT-PCR, GAPDH was used as housekeeping gene.
Estimation of Total starch, Amylose and protein content
Total starch content of grounded kernels was determined by Megazymes Total Starch Assay Kit (AA/AMG, Megazyme Pty Ltd., Wicklow, Ireland).
Protein content was determined by the Dumas combustion method, using a Leco FP 428 analyzer (Leco Co., St. Joseph, MI) calibrated with an EDTA standard.
Protein extraction and analysis
The preparation of starch granules from half seeds and the separation of starch granule bound proteins by SDS-PAGE followed the method reported by Zhao and Sharp  with some modifications, as in Mohammadkhani et al. . Protein bands were visualized by silver staining.
Starch paste viscosity was measured by the rapid visco analyzer (RVA series 4, Newport Scientific, Sidney, Australia). Whole durum wheat flour (3,87 g) was mixed with 25,13 mL of 0,2% silver nitrate in a RVA canister. The starch suspension was, constantly, stirred at 160 rpm and heated from 50 to 95°C, held at 95°C for 4 min and then cooled to 50°C for 4 min. The following variables were recorded in duplicate: peak viscosity (PV), breakdown (BK), trough (TG), final paste viscosity (FV), peak time (Pt), pasting temperature (PT), and setback (SB). Viscosity was expressed in Rapid Visco Unit (RVU).
Starch was isolated from 5 seeds as described in Zhao and Sharp  with same modification. Starch was resuspended in 1000 μl extraction buffer (Tris 0,5 M; SDS 10%, glycerol 50%, sucrose and DTT 1%) and boiled for 15 min. Gelled starch solution was then cooled at -20°C for 1 hour and centrifuged at 13000 rpm for 10 min. Four volumes of acetone were added to the supernatant that was stored at -20°C overnight. The samples was centrifuged at 4°C at 8500 rpm for 40 min; pellet was washed with 2,5 ml of aceton for three times, and finally air-dried. Precipitated proteins were resuspended in 35 μl extraction buffer and loaded on SDS-PAGE (T = 10; C = 1,28), running at 200 V for 1 h and 35 minute. Protein bands were visualized by colloidal coomassie CBBG-250 , rinsed out with mill-Q water filter, and captured with sterile scalpel. All steps were conducted in sterile atmosphere, making attention to keratin contaminations. The captured bands were placed in a eppendorf and submitted by Mass spectrometry analysis (University of California, Davis, CA, U.S.A). Scaffold (version Scaffold-2.02.01) was used to analyze Tandem MS/MS based peptides and to identify the proteins. Peptide identifications were accepted if they could be established at greater than 95.0% probability as specified by the Peptide Prophet algorithm . Protein identifications were accepted if they could be established at greater than 99.0% probability and contained at least two identified peptides. Protein probabilities were assigned by the Protein Prophet algorithm .
The proteins were prepared for MS analysis using standard reduction, alkylation, and tryptic digest procedures . Peptides were dried down, then resolubilized in 2% acetonitrile/0.1% trifluoroacetic acid for LC-MS/MS analysis.
Digested peptides were analyzed by LC-MS/MS on a LTQ-FT with Waters nano Acuity UPLC. Peptides were separated with a 90 min gradient using a Waters 100 um × 100 mm 1.7 um BEH130 C18 reversed phase column at 2 μl/min. The LTQ-FT was operated with a top 4 method, where the first scan was acquired in the FT at 100,000 resolution and the top 4 ions in each survey scan were selected for low energy CID in the LTQ. Solvent A composition was 0.1% formic acid, and solvent B 100% ACN. Peptides were separated at 2 μl/min on the reversed-phase column with this gradient: 0 to 7% Solvent B in 5 minutes, 7-35% in 65 minutes, 35-70% in 5.2 minutes, held at 70% for 0.8 minutes, 70-1% in 0.2 minutes, with 13.8 minutes to equilibrate between runs.
X!Tandem was used to search the MS/MS spectra against both the provided database. A parent mass error of up to 10 ppm was allowed for the searches, and a product mass error of 0.4 Daltons. Complete modification of carbamidomethyl was used, as well as the potential modification of methionine oxidation.
Microscopic Analyses of Starch Granules
Starch granules were extracted from seeds as reported in Zhao and Sharp  with some modifications, as in Mohammadkhani et al. . Samples were critical-point-dried in a Balzer's apparatus equipped with a liquid CO2 inlet and metal-shadowed in a gold sputtering unit equipped with an argon inlet. Specimens were examined in a Jeol JSM 5200 scanning electron microscope.
Rothamsted Research receives grant-aided support from the Biotechnology and Biological Sciences Research Council (BBSRC) of the UK.
- Ellis RP, Cochrane MP, Dale MFB, Duffus CM, Lynn A, Morrison IM, Prentice RDM, Swanston JS, Tiller SA: Starch production and industrial use. J Sci Food Agric. 1998, 77: 289-311. 10.1002/(SICI)1097-0010(199807)77:3<289::AID-JSFA38>3.0.CO;2-D.View ArticleGoogle Scholar
- Ball SG, Morell MK: From Bacterial Glycogen to Starch: Understanding the Biogenesis of the Plant Starch Granule. Annu Rev Plant Biol. 2003, 54: 207-233. 10.1146/annurev.arplant.54.031902.134927.PubMedView ArticleGoogle Scholar
- Stone B, Morell MK: Carbohydrate. Wheat: Chemistry and Technology. Edited by: Khan K, Shewry PR, AACC International St. Paul MN,2009. 299-362. 4View ArticleGoogle Scholar
- Stinard PS, Robertson DS, Schnable PS: Genetic isolation, cloning, and analysis of a mutator-induced, dominant antimorph of the maize amylose extender1 locus. Plant Cell. 1993, 5: 1555-1566. 10.1105/tpc.5.11.1555.PubMedPubMed CentralView ArticleGoogle Scholar
- Mizuno K, Kawasaki T, Shimada H, Satoh H, Kobayashi E, Okumura S, Arai Y, Baba T: Alteration of the structural properties of starch components by the lack of an isoform of starch branching enzyme in rice seeds. J Biol Chem. 1993, 268: 19084-19091.PubMedGoogle Scholar
- Bhattacharyya MK, Smith AM, Ellis THN, Hedley C, Martin C: The wrinkled-seed character of pea described by Mendel is caused by a transposon-like insertion in a gene encoding starch-branching enzyme. Cell. 1990, 60: 115-122. 10.1016/0092-8674(90)90721-P.PubMedView ArticleGoogle Scholar
- Blauth SL, Yao Y, Klucinec JD, Shannon JC, Thompson DB, Guiltinan MJ: Identification of Mutator insertional mutants of starch-branching enzyme 2a in corn. Plant Physiol. 2001, 125: 1396-1405. 10.1104/pp.125.3.1396.PubMedPubMed CentralView ArticleGoogle Scholar
- Blauth SL, Kim KN, Klucinec J, Shannon JC, Thompson D, Guiltinan M: Identification of Mutator insertional mutants of starch-branching enzyme 1 (sbe1) in Zea mays L. Plant Mol Biol. 2002, 48: 287-297. 10.1023/A:1013335217744.PubMedView ArticleGoogle Scholar
- Satoh H, Nishi A, Yamashita K, Takemoto Y, Tanaka Y, Hosaka Y, Sakurai A, Fujita N, Nakamura Y: Starch-Branching Enzyme I-Deficient Mutation Specifically Affects the Structure and Properties of Starch in Rice Endosperm. Plant Physiol. 2003, 133: 1111-1121. 10.1104/pp.103.021527.PubMedPubMed CentralView ArticleGoogle Scholar
- Nair RB, Baga M, Scoles GJ, Kartha KK, Chibbar RN: Isolation, characterization and expression analysis of a starch branching enzyme II cDNA from wheat. Plant Sci. 1997, 122: 153-163. 10.1016/S0168-9452(96)04543-8.View ArticleGoogle Scholar
- Rahman S, Regina A, Li Z, Mukai Y, Yamamoto M, Kosar-Hashemi B, Abrahams S, Morell MK: Comparison of starch-branching enzyme genes reveals evolutionary relationships among isoforms. Characterization of a gene for starch-branching enzyme IIa from the wheat genome donor Aegilops tauschii. Plant Physiol. 2001, 125: 1314-1324. 10.1104/pp.125.3.1314.PubMedPubMed CentralView ArticleGoogle Scholar
- Regina A, Kosar-Hashemi B, Li Z, Pedler A, Mukai Y, Yamamoto M, Gale K, Sharp PJ, Morell MK, Rahman S: Starch branching enzyme IIb in wheat is expressed at low levels in the endosperm compared to other cereals and encoded at a non-syntenic locus. Planta. 2005, 222: 899-909. 10.1007/s00425-005-0032-z.PubMedView ArticleGoogle Scholar
- Gao M, Fisher DK, Kim KN, Shannon JC, Guiltinan MJ: Independent genetic control of maize starch-branching enzymes IIa and IIb - Isolation and characterization of a Sbe2a cDNA. Plant Physiol. 1997, 114: 69-78. 10.1104/pp.114.1.69.PubMedPubMed CentralView ArticleGoogle Scholar
- Yamanouchi H, Nakamura Y: Organ specificity of isoforms of starch branching enzyme (Q-enzyme) in rice. Plant Cell Physiol. 1992, 33: 985-991.Google Scholar
- Regina A, Bird A, Topping D, Bowden S, Freeman J, Barsby T, Kosar-Hashemi B, Li Z, Rahman S, Morell MK: High-amylose wheat generated by RNA interference improves indices of large-bowel health in rats. Proc Natl Acad Sci USA. 2006, 103: 3546-3551. 10.1073/pnas.0510737103.PubMedPubMed CentralView ArticleGoogle Scholar
- Englyst HN, Kingman SM, Hudson GJ, Cummings JH: Measurement of resistant starch in vitro and in vivo. Br J Nutr. 1996, 75: 749-755. 10.1079/BJN19960178.PubMedView ArticleGoogle Scholar
- Nugent AP: Health properties of resistant starch. Nutr Bull. 2005, 30: 27-54. 10.1111/j.1467-3010.2005.00481.x.View ArticleGoogle Scholar
- Topping DL, Morell MK, King RA, Li Z, Bird AR, Noakes M: Resistant Starch and Health - Himalaya 292, a Novel Barley Cultivar to Deliver Benefits to Consumers. Starch/Stärke. 2003, 55: 539-545. 10.1002/star.200300221.View ArticleGoogle Scholar
- Bingham SA, Day NE, Luben R, Ferrari P, Slimani N, Norat T, Clavel-Chapelon F, Kesse E, Nieters A, Boeing H, Tjønneland A, Overvad K, Martinez C, Dorronsoro M, Gonzalez CA, Key TJ, Trichopoulou A, Naska A, Vineis P, Tumino R, Krogh V, Bueno-de-Mesquita HB, Peeters PH, Berglund G, Hallmans G, Lund E, Skeie G, Kaaks R, Riboli E: Dietary fibre in food and protection against colorectal cancer in the European Prospective. Investigation into Cancer and Nutrition (EPIC): an observational study. Lancet. 2003, 361: 1496-1501. 10.1016/S0140-6736(03)13174-1.PubMedView ArticleGoogle Scholar
- Lintas C, Cappelloni M, Adorisio S, Clementi A, Del Toma E: Effect of ripening on resistant starch and total sugar in Musa paradisiaca sapientum: glycaemic and insulinaemic responses in normal subjects and NIDDM patients. Eur J Clin Nutr. 1995, 49: 303-306.Google Scholar
- Lintas C, Cappeloni M, Bonsmassar L, Clementi A, Del Toma E, Ceccarelli G: Dietary fibre, resistant starch and in vitro digestibility of cereal meals. Glycaemic and insulinaemic responses in NIDDM patients. Eur J Clin Nutr. 1995, 49: 264-267.Google Scholar
- Van Hung P, Yamamori M, Morita N: Formation of enzyme-resistant starch in bread as affected by high-amylose wheat flour substitutions. Cereal Chem. 2005, 82: 690-694. 10.1094/CC-82-0690.View ArticleGoogle Scholar
- Soh HN, Sisson MJ, Turner M: Effect of starch granule size distribution and elevated amylose content on durum dough rheology and spaghetti cooking quality. Cereal Chem. 2006, 83: 513-519. 10.1094/CC-83-0513.View ArticleGoogle Scholar
- Yamamori M, Nakamura T, Kuroda A: Variations in the content of starch-granule bound protein among several Japanese cultivars of common wheat (Triticum aestivum L.). Euphytica. 1992, 64: 215-219. 10.1007/BF00046051.View ArticleGoogle Scholar
- Urbano M, Margiotta B, Colaprico G, Lafiandra D: Waxy protein in diploid, tetraploid and hexaploid wheats. Plant Breeding. 2002, 121: 1-5. 10.1046/j.1439-0523.2002.00758.x.View ArticleGoogle Scholar
- Yamamori M, Fujita S, Hayakawa K, Matsuki J, Yasui T: Genetic elimination of a starch granule protein, SGP-1, of wheat generates an altered starch with apparent high amylose. Theor Appl Genet. 2000, 101: 21-29. 10.1007/s001220051444.View ArticleGoogle Scholar
- Slade AJ, Fuerstenberg SI, Loeffler D, Steine MN, Facciotti D: A reverse genetic, nontransgenic approach to wheat crop improvement by TILLING. Nat Biotechnol. 2005, 23: 75-81. 10.1038/nbt1043.PubMedView ArticleGoogle Scholar
- Sestili F, Botticella E, Bedo Z, Phillips A, Lafiandra D: Production of novel allelic variation for genes involved in starch biosynthesis through mutagenesis. Mol Breeding. 2010, 25: 145-154. 10.1007/s11032-009-9314-7.View ArticleGoogle Scholar
- Fu D, Uauy C, Blechl A, Dubcovsky J: RNA interference for wheat functional gene analysis. Transgenic Res. 2007, 16: 689-701. 10.1007/s11248-007-9150-7.PubMedView ArticleGoogle Scholar
- Yue SJ, Li H, Li YW, Zhu YF, Guo JK, Liu YJ, Chen Y, Jia X: Generation of transgenic wheat lines with altered expression levels of 1Dx5 high-molecular weight glutenin subunit by RNA interference. J Cereal Sci. 2008, 47: 153-161. 10.1016/j.jcs.2007.03.006.View ArticleGoogle Scholar
- Travella S, Klimm TE, Keller B: RNA interference-based gene silencing as an efficient tool for functional genomics in hexaploid bread wheat. Plant Physiol. 2006, 142: 6-20. 10.1104/pp.106.084517.PubMedPubMed CentralView ArticleGoogle Scholar
- Garcia-Perez RD, Houdt HV, Depicker A: Spreading of post-transcriptional gene silencing along the target gene promotes systemic silencing. Plant J. 2004, 38: 594-602. 10.1111/j.1365-313X.2004.02067.x.PubMedView ArticleGoogle Scholar
- McGinnis K, Murphy N, Carlson AR, Akula A, Akula C, Basinger H, Carlson M, Hermanson P, Kovacevic N, McGill MA, Seshadri V, Yoyokie J, Cone K, Kaeppler HF, Kaeppler SM, Springer NM: Assessing the efficiency of RNA interference for maize functional genomics. Plant Physiol. 2007, 143: 1441-1451. 10.1104/pp.106.094334.PubMedPubMed CentralView ArticleGoogle Scholar
- Kerschen A, Napoli CA, Jorgensen RA, Muller AE: Effectiveness of RNA interference in transgenic plants. FEBS Lett. 2004, 566: 223-228. 10.1016/j.febslet.2004.04.043.PubMedView ArticleGoogle Scholar
- Shimada T, Otani M, Hamada T, Kim SH: Increase of amylase content of sweetpotato starch by RNA interference of the starch branching enzyme II gene (IbSBEII). Plant Biotechnol. 2006, 23: 85-90.View ArticleGoogle Scholar
- Andersson M, Melander M, Pojmark P, Larsson H, Bulow L, Hofvander P: Targeted gene suppression by RNA interference: an efficient method for production of high-amylose potato lines. J Biotechnol. 2006, 123: 137-148. 10.1016/j.jbiotec.2005.11.001.PubMedView ArticleGoogle Scholar
- McClearly BV, Solah V, Gibson TS: Quantitative measurement of total starch in cereal flours and products. J Cereal Sci. 1994, 20: 51-58. 10.1006/jcrs.1994.1044.View ArticleGoogle Scholar
- Yamamori M, Kato M, Yui M, Kawasaki M: Resistant starch and starch pasting properties of a starch synthase IIa-deficient wheat with apparent high amylose. Aust J Agric Res. 2006, 57: 531-535. 10.1071/AR05176.View ArticleGoogle Scholar
- Juhász R, Salgó A: Pasting behavior of amylose, amylopectin and their mixtures as determined by RVA curves and first derivatives. Starch/Staerke. 2008, 60: 70-78. 10.1002/star.200700634.View ArticleGoogle Scholar
- You SG, Kim SM: Molecular Characteristics and Functional Properties of Barley Starches with Varying Amylose Content. J Food Sci Nutr. 2003, 10: 207-213.Google Scholar
- Morita N, Maeda T, Miyazaki M, Yamamori M, Miura H, Ohtsuka I: Dough and baking properties of high-amylose and waxy wheat flours. Cereal Chem. 2002, 79: 491-495. 10.1094/CCHEM.2002.79.4.491.View ArticleGoogle Scholar
- Kosar-Hashemi B, Li Z, Larroque O, Regina A, Yamamori M, Morell MK, Rahman S: Multiple effects of the starch synthase II mutation in developing wheat endosperm. Funct Plant Biol. 2007, 34: 431-438. 10.1071/FP06288.View ArticleGoogle Scholar
- Nakamura Y, Umemoto T, Takahata Y, Komae K, Amano E, Satoh H: Changes in structure of starch and enzyme activities affected by sugary mutations in developing rice endosperm. Physiol Plantarum. 1996, 97: 491-498. 10.1111/j.1399-3054.1996.tb00508.x.View ArticleGoogle Scholar
- Gao M, Wanat J, Stinard PS, James MG, Myers AM: Characterization of dull1, a maize gene coding for a novel starch synthase. Plant Cell. 1998, 10: 399-412. 10.1105/tpc.10.3.399.PubMedPubMed CentralView ArticleGoogle Scholar
- Stahl Y, Coates S, Bryce JH, Morris PC: Antisense down regulation of the barley limit dextrinase inhibitor modulates starch granule size distribution, starch composition and amylopectin structure. Plant J. 2004, 39: 599-611. 10.1111/j.1365-313X.2004.02159.x.PubMedView ArticleGoogle Scholar
- Tetlow IJ, Wait R, Lu Z, Akkasaeng R, Bowsher CG, Esposito S, Kosar-Hashemi B, Morell MK, Emes MJ: Protein phosphorylation in amyloplasts regulates starch branching enzyme activity and protein-protein interactions. Plant Cell. 2004, 16: 694-708. 10.1105/tpc.017400.PubMedPubMed CentralView ArticleGoogle Scholar
- Tetlow IJ, Beisel KG, Cameron S, Makhmoudova A, Liu F, Bresolin NS, Wait R, Morell MK, Emes MJ: Analysis of protein complexes in wheat amyloplasts reveals functional interactions among starch biosynthetic enzymes. Plant Physiol. 2008, 146: 1878-1891. 10.1104/pp.108.116244.PubMedPubMed CentralView ArticleGoogle Scholar
- Hennen-Bierwagen TA, Liu F, Marsh RS, Kim S, Gan Q, Tetlow IJ, Emes MJ, James MG, Myers AM: Starch biosynthetic enzymes from developing maize endosperm associate in multisubunit complexes. Plant Physiol. 2008, 146: 1892-1908. 10.1104/pp.108.116285.PubMedPubMed CentralView ArticleGoogle Scholar
- Hennen-Bierwagen TA, Lin Q, Grimaud F, Planchot V, Keeling PL, James MG, Myers AM: Proteins from Multiple Metabolic Pathways Associate with Starch Biosynthetic Enzymes in High Molecular Weight Complexes: A Model for Regulation of Carbon Allocation in Maize Amyloplasts. Plant Physiol. 2009, 149: 1541-1559. 10.1104/pp.109.135293.PubMedPubMed CentralView ArticleGoogle Scholar
- Grimaud F, Rogniaux H, James MG, Myers AM, Planchot V: Proteome and phosphoproteome analysis of starch granule-associated proteins from normal maize and mutants affected in starch biosynthesis. J Exp Bot. 2008, 59: 3395-3406. 10.1093/jxb/ern198.PubMedPubMed CentralView ArticleGoogle Scholar
- Tosi P, D'Ovidio R, Napier JA, Bekes F, Shewry PR: Expression of epitope-tagged LMW glutenin subunits in the starchy endosperm of transgenic wheat and their incorporation into glutenin polymers. Theor Appl Genet. 2004, 108: 468-476. 10.1007/s00122-003-1459-x.PubMedView ArticleGoogle Scholar
- Christensen AH, Quail PF: Ubiquitin promoter-based vectors for high-level expressin of selectable and/or screenable marker genes in monocotyledonous plants. Trans Res. 1996, 5: 213-218. 10.1007/BF01969712.View ArticleGoogle Scholar
- Okubara P, Blechl AE, McCormick SP, Alexander N, Dill-Macky R, Hohn T: Engineering deoxynivalenol metabolism in wheat through the expression of a fungal trichothecene acetyltransferase gene. Theor Appl Genet. 2002, 106: 74-83.PubMedGoogle Scholar
- Wu X, Doherty A, Jones HD: Efficient and rapid Agrobacterium-mediated genetic transformation of durum wheat (Triticum turgidum L. var. durum) using additional virulence genes. Transgenic Res. 2008, 17: 425-436. 10.1007/s11248-007-9116-9.PubMedView ArticleGoogle Scholar
- Tai TH, Tanksley SD: A rapid and inexpensive method for isolation of total DNA from dehydrated plant tissue. Plant Mol Biol Rep. 1990, 8: 297-303. 10.1007/BF02668766.View ArticleGoogle Scholar
- Laudencia-Chingcuanco DL, Stamova BS, You FM, Lazo GR, Beckles DM, Anderson OD: Transcriptional profiling of wheat caryopsis development using cDNA microarrays. Plant Mol Biol. 2006, 63: 651-668. 10.1007/s11103-006-9114-y.View ArticleGoogle Scholar
- Jin YF, Bian TF: Isolation and partial characterization of a novel pollen-specific cDNA with multiple polyadenylation sites from wheat. Acta Biochim Biophys Sin. 2004, 36: 467-476. 10.1093/abbs/36.7.467.PubMedView ArticleGoogle Scholar
- Hurkman WJ, McCue KF, Altenbach SB, Korn A, Tanaka CK, Kothari KM, Johnson EL, Bechtel DB, Wilson JD, Anderson OD, DuPont FM: Effect of temperature on expression of genes encoding enzymes for starch biosynthesis in developing wheat endosperm. Plant Sci. 2003, 164: 873-881. 10.1016/S0168-9452(03)00076-1.View ArticleGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.PubMedView ArticleGoogle Scholar
- Chrastil J: Improved colourimetric determination of amylose in starches or flours. Carbohydr Res. 1987, 159: 154-158. 10.1016/S0008-6215(00)90013-2.View ArticleGoogle Scholar
- Zhao XC, Sharp PJ: An improved 1-D SDS-PAGE method for the identification of three bread wheat waxy proteins. J Cereal Sci. 1996, 23: 191-193. 10.1006/jcrs.1996.0019.View ArticleGoogle Scholar
- Mohammadkhani A, Stoddard FL, Marshall DR, Uddin MN, Zhao X: Starch extraction and amylose analysis from half seeds. Starch/Stärke. 1999, 51: 62-66. 10.1002/(SICI)1521-379X(199903)51:2<62::AID-STAR62>3.0.CO;2-G.View ArticleGoogle Scholar
- Neuhoff V, Arold N, Taube D, Ehrhardt W: Improved staining of proteins in polyacrylamide gels including isoelectric focusing gels with clear background at nanogram sensitivity using Coomassie Brilliant Blue G-250 and R-250. Electrophoresis. 1998, 9: 255-262. 10.1002/elps.1150090603.View ArticleGoogle Scholar
- Keller A, Nesvizhskii AI, Kolker E, Aebersold R: Empirical statistical model to estimate the accuracy of peptide identifications made by MS/MS and database search. Anal Chem. 2002, 74: 5383-5392. 10.1021/ac025747h.PubMedView ArticleGoogle Scholar
- Nesvizhskii AI, Keller A, Kolker E, Aebersold RA: Statistical model for identifying proteins by tandem mass spectrometry. Anal Chem. 2003, 75: 4646-4658. 10.1021/ac0341261.PubMedView ArticleGoogle Scholar
- Rosenfeld J, Capdevielle J, Guillemot JC, Ferrara P: In-gel digestion of proteins for internal sequence analysis after one- or two-dimensional gel electrophoresis. Anal Biochem. 1992, 203: 173-179. 10.1016/0003-2697(92)90061-B.PubMedView ArticleGoogle Scholar
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.