- Research Article
- Open Access
A PSTOL-like gene, TaPSTOL, controls a number of agronomically important traits in wheat
BMC Plant Biology volume 18, Article number: 115 (2018)
Phosphorus (P) is an essential macronutrient for plant growth, and is required in large quantities by elite varieties of crops to maintain yields. Approximately 70% of global cultivated land suffers from P deficiency, and it has recently been estimated that worldwide P resources will be exhausted by the end of this century, increasing the demand for crops more efficient in their P usage. A greater understanding of how plants are able to maintain yield with lower P inputs is, therefore, highly desirable to both breeders and farmers. Here, we clone the wheat (Triticum aestivum L.) homologue of the rice PSTOL gene (OsPSTOL), and characterize its role in phosphate nutrition plus other agronomically important traits.
TaPSTOL is a single copy gene located on the short arm of chromosome 5A, encoding a putative kinase protein, and shares a high level of sequence similarity to OsPSTOL. We re-sequenced TaPSTOL from 24 different wheat accessions and (3) three T. durum varieties. No sequence differences were detected in 26 of the accessions, whereas two indels were identified in the promoter region of one of the durum wheats. We characterised the expression of TaPSTOL under different P concentrations and demonstrated that the promoter was induced in root tips and hairs under P limiting conditions. Overexpression and RNAi silencing of TaPSTOL in transgenic wheat lines showed that there was a significant effect upon root biomass, flowering time independent of P treatment, tiller number and seed yield, correlating with the expression of TaPSTOL. However this did not increase PUE as elevated P concentration in the grain did not correspond to increased yields.
Manipulation of TaPSTOL expression in wheat shows it is responsible for many of the previously described phenotypic advantages as OsPSTOL except yield. Furthermore, we show TaPSTOL contributes to additional agronomically important traits including flowering time and grain size. Analysis of TaPSTOL sequences from a broad selection of wheat varieties, encompassing 91% of the genetic diversity in UK bread wheat, showed that there is very little genetic variation in this gene, which would suggest that this locus may have been under high selection pressure.
Phosphorus (P) is an essential macronutrient that is required for all major developmental processes in plants and is considered to be one of the most limiting plant nutrients to global agricultural production. Approximately 70% of global cultivated land suffers from phosphate deficiency, making research into phosphate nutrition a priority for both scientists and farmers [1,2,3]. Phosphate is usually taken up in the orthophosphate forms (Pi) which are available at low concentrations in the soil solution . This low availability can be compounded by the soil chemistry in which soils high in clay can bind the Pi requiring large amounts of fertilizer to be applied to maintain high yields [4,5,6]. To compound the problem of P limiting growth and yield, recent estimates suggest that only 20–30% of the inorganic phosphate (Pi) added to fields by farmers is assimilated by plants . An increased understanding of how plants use Pi and the current inefficiency of Pi assimilation would be of economic benefit to farmers and also potentially alleviate environmental problems resulting from large scale agricultural production. Wheat (Triticum aestivum L.) is one of the most important crops worldwide and grown on more than 240 million hectares in 2016 . P is one of the three most important nutrients for growth and yield improvement in wheat [9, 10]. Therefore an improvement in phosphate use efficiency (PUE) in wheat could potentially have a large global impact over 200 million hectares.
Recent advances in our understanding of the molecular mechanisms by which different plant species adapt to low-phosphate stress, the regulation and expression of phosphorus metabolism genes and alleles to deal with phosphate limitation, have enabled the design of more effective breeding strategies to produce highly phosphate efficient crops [11, 12]. A number of the underlying genes involved in the response of plants to low P are highly conserved and play similar roles in a number of diverse plant species including both model system and crop species. Thus, a greater understanding of the pathways involved in phosphate acquisition and signalling will allow breeders and plant molecular biologists to develop more efficient crops. Identification of these conserved pathways and genes from model organisms, and the subsequent transfer of this knowledge to crop species, would therefore allow farmers to optimize fertilizer use, resulting in increased food production efficiency with lower environmental cost.
One such locus which is believed to be important is the recent identification of the PUP1 locus from rice which contains a putative kinase gene called Phosphate Starvation Tolerance 1 (PSTOL) [13,14,15,16]. The PUP1 locus was originally identified in an upland variety of rice, Kasalath, yet is absent from most rice cultivars [17,18,19]. Rice varieties which have this genomic introgression containing the PSTOL gene, show increased biomass, increased root growth, increased tiller number and yield increases of up to 30% when grown under low P conditions whereas no deleterious consequences were seen when grown under normal soil fertility conditions [19, 20]. The identification of OsPSTOL and its role in helping rice tolerate low P conditions has led to the belief that we can engineer PUE into many crop species using translational science from gains in other model species.
Others have identified homologous PSTOL-like genes in both maize and sorghum, based upon QTL analysis and sequence homology. Further evidence of the role of PSTOL-like genes in PUE has been supported through QTL mapping rather than direct molecular characterization of candidate genes [21, 22]. Criteria to identify other potential PSTOL like genes has included identification of protein domains such as ATP kinase domains and on DNA sequence conservation meeting certain bioinformatic cut-offs for genes underlying these QTL. However some of these homologues appear to have differences in their gene structure such as number and length of introns and UTRs. Despite the lack of a highly conserved gene in most plant species, it is critical to understand whether other PSTOL genes exist in other important crop species and if they can be exploited.
Our study identifies a putative wheat PSTOL gene (TaPSTOL) and characterizes its role in PUE and other phenotypes of agronomic importance.
The OsPSTOL coding sequence from EMBL (AB458444.1) was used as a query for BLAST searches of the wheat genome (IWGSC 2014) and expressed sequence tag (EST) databases in ensembl, GenBank, Komugi and URGI [23,24,25,26,27]. Gene prediction was carried out using FGENESH . Protein domain prediction was found using IntroProScan .
Nullisomic (N) / Tetrasomic (T) wheat lines
PCR conditions for the verification of TaPSTOL homeologues were first optimized with wheat cv. Fielder DNA using a gradient from 53 to 63 °C Tm with and without 2% DMSO (final concentration) using FastStart™ Taq DNA Polymerase (SIGMA). All primers used in this study are listed in Additional file 1: Table S1. All three primer pairs amplified the predicted product under all conditions tested (Additional file 2: Figure S1). To reduce the stringency, the annealing temperature was then reduced to 55 °C to allow for small sequence variations and amplification of other distant homeologues using the nullisomic / tetrasomic wheat DNA .
Creation of plasmid constructs
TaPSTOL was PCR amplified from genomic DNA extracted from wheat cultivar Chinese Spring using primers TaPUP-K46-F and TaPUP-K46-R with Phusion Hotstart II DNA Polymerase (Thermofisher). Amplicons were cloned into pJET2.1 (Thermofisher) and fully sequenced. The full length gene was then reamplified from pJET-TaPSTOL with TaPUP-K46-F-PvuII and TaPUP-K46-R-XbaI to add PvuII and XbaI restriction sites plus a monocot ribosome binding site, CCACC , and cloned into pENTR-1A digested with DraI and XbaI (Thermofisher). TaPSTOL was then recombined into the binary vector pSc4ActR1R2  using a Gateway LR Clonase II Kit (Thermofisher) to create pRMH007. TaPSTOL was expressed in planta from the rice Actin promoter  and transcripts terminated by the A. tumefaciens nopaline synthase terminator (tNOS). The first 350 bp of the TaPSTOL sequence selected to trigger silencing by RNAi, was amplified from the genomic sequence using the primers RNAi-F GW and RNAi-R GW, and recombined directly into pDONR221 using a Gateway BP Clonase II kit (Invtrogen). The pDONR TaPSTOL insert was subsequently recombined into the binary vector pACT-IR2 in a Gateway LR reaction to create pMM2, with the TaPSTOL RNAi hairpin cassette expressed from the rice Actin promoter in planta.
The 2.4 kb promoter region upstream of the TaPSTOL start codon was amplified from genomic DNA extracted from Chinese Spring using primers TaPUP-prom–F and TaPUP-prom-R with Phusion Hotstart II Polymerase (Thermofisher). The resultant amplicon was ligated into pCR-Blunt (Thermofisher) and sequenced. The promoter region was then reamplified with primers TaPUP-prom-13R and TaPUP-prom-SalI. The amplicon was digested with SalI and ligated into pRMH013 digested with BmgBI and SalI. This intermediate vector containing the TaPSTOL promoter driving GUS and flanked by attL sites, was then recombined into pRLF12-R1R2-SCV to produce vector pRMH107. Additional file 3: Figure S2, shows the T-DNA regions from these constructs, which were transformed into wheat in this study.
Completed constructs were verified by restriction digest and sequencing before being electro-transformed into Agrobacterium tumefaciens. Plasmids were re-isolated from Agrobacterium cultures and verified by restriction digest prior to use in wheat experiments .
Wheat cv. Fielder plants were grown in controlled environment chambers (Conviron) at 20 °C day/15 °C night with a 16 h day photoperiod (approximately 400 μE m− 2 s− 1). Immature seeds were harvested for transformation experiments at 14–20 days post-anthesis (dpa). Isolated immature wheat embryos were co-cultivated with Agrobacterium tumefaciens for 2 days in the dark . Subsequent removal of the embryonic axis and tissue culture was performed as previously described . Individual plantlets were hardened off following transfer to Jiffy-7 pellets (LBS Horticulture), potted up into 9 cm plant pots containing M2 compost plus 5 g/l slow release fertilizer (Osmocote Exact 15:9:9) and grown on to maturity and seed harvest in controlled environment chambers, as above.
DNA analysis of transformed wheat plants
Plantlets which regenerated under G418 selection in tissue culture were transferred to Jiffy-7 pellets and validated using an nptII copy number assay relative to a single copy wheat gene amplicon, GaMyb, normalised to a known single copy wheat line. Primers and Taqman probes were used at a concentration of 10 μM in a 10 μl multiplex reaction using ABsolute Blue qPCR ROX mix (Thermofisher) with the standard run conditions for the ABI 7900 HT. The relative quantification, ΔΔCt, values were calculated to determine nptII copy number in the T0 and subsequent generations . Homozygous and null transgenic lines were identified on the basis of nptII copy number and segregation analysis. WT Fielder plants were null segregates.
Plant growth conditions
WT and transgenic lines were grown under low P conditions in sand and fertilized with nutrient solution  or under fertilized M2 compost conditions, as above. Total dry shoot weight, seed weight (yield per plant), seed number, seed size, tiller number and P concentration via ICP-MS were measured. Biological replicates each contained 14 plants per line and were grown until seed maturation. Tissues were allowed to dry for a further two weeks before harvesting. For low P conditions plants were germinated in sand and watered with 25 mL of a Magnavaca solution containing 3 μM KH2PO4, 1.3 mM NH4NO3, 3.52 mM Ca(NO3)2, 0.58 mM, KCl, 0.58 mM K2SO4, 0.56 mM KNO3, 0.86 mM Mg(NO3)2 0.13 mM H3BO3, 5 μM MnCl2, 0.4 μM Na2MoO4, 10 μM ZnSO4, 0.3 μM CuSO4, Fe(NO3)3 and 2 mM MES (pH 5.5) twice a week until maturity. Plants were grown in a controlled growth chamber under 16 h light and 20 °C/15 °C day night temperatures.
RNA expression analysis
Fielder seedlings were grown for seven days in 2.2 L pots containing Magnavaca solution as listed above and supplemented with either 1, 2.5, 5, 20 or 45 μM KH2PO4. Plants were grown for 7 days before harvesting tissue and separating the samples into root and shoot tissues for analysis. Total RNA was isolated from both roots and shoots for each P treatment using an RNeasy Kit (Qiagen) and treated with DNaseI (Thermofisher) prior to cDNA synthesis from 500 ng of total RNA using Omniscript RT Kit (Qiagen). The cDNA was diluted 1:2 with water and 0.5 μL was used as template in each RT-PCR reaction. Expression levels were quantified by quantitative PCR in triplicate reactions from three biological replications using SYBR Green JumpStartTaq ReadyMix (SIGMA) with the standard run conditions for the ABI 7900 HT. TaPSTOL expression was compared to two reference genes TaUbiquitin and TaEF1α. Primers used for amplification of transcripts were TaPSTOL-Q-F and TaPSTOL-Q-R, Ubi-F and Ubi-R for Ubiquitin  or EF1α-F and EF1α-R for EF1α . Data shown is in comparison to Ubiquitin for ease as both reference genes showed similar differences in expression.
Primers were designed to amplify a 3.3 kb fragment of TaPSTOL which included the 5′ upstream region and ORF from DNA extracted from Chinese Spring using the primers TaPUP-prom–F and TaPUP-prom-R. Genomic DNA was amplified from bread wheat varieties Alchemy, Banco, Bersee, Bridgadier Brompton, Claire, Copain, Cordiale, Fielder, Flamingo, Gladiator, Hereward, Holdfast, Kloka, Maris Fundin, Paragon, Rialto, Robigus, Slejpner, Soissons, Spark, Steadfast, Stetson, Xi19 and Chinese Spring, plus three wild T. turgidum ssp. dicoccoides accessions (PI 503314, PI 414722, and PI 428097) sourced from the USDA-ARS National Small Grains Collection. Genomic DNA was extracted using the Tanksley method .
Gus staining of wheat tissues
Whole 10-day seedlings grown on sand watered once with hydroponic solution containing 3 μM P, were stained in X-Gluc solution (0.1 M NaPO4, 10 mM EDTA pH 7.0, 0.5 mM K Ferricyanide, 0.5 mM K Ferrocyanide, 1.0 mM X-Glucuronide, 0.1% Triton X100, pH 7.0) overnight at 37 °C, then destained in 70% ethanol .
Digestion and elemental analysis of plant material
Dried samples of leaf, root and grain (~ 0.2 g) were digested and analysed for elemental content as described by Thomas et al. .
Phosphate use efficiency calculations
Definitions of PPUE and PER were taken from . PPUE was calculated as yield/ P concentration for an individual treatment. PER was calculated as yield / (P concentration * yield) for an individual treatment. ANOVAs were run using R and the aov and Tukey functions with the null hypothesis of no difference between lines. Tukey’s post hoc test was added to identify each significant interaction between the lines tested.
Identification of TaPSTOL
Our TaPSTOL candidate was identified following BLAST searches  using the OsPSTOL coding sequence to the first wheat genome release (IWGSC 2014 and EST collections; this has subsequently been annotated as gene model Traes_5AS_AA3DC6A5F on ENSEMBL (IWGSC, 2014 release, Ensembl). The identified homologue is located on the short arm of chromosome 5A and shares 90% homology at the DNA level with the rice gene sequence. Further analysis showed that the rice and wheat predicted proteins share 74.2% identity and 92.7% similarity at the amino acid level (Additional file 4: Figure S3). OsPSTOL and TaPSTOL both possess single exons, and their predicted proteins are similar in size: 324 amino acids and 289 amino acids, respectively. InterProScan identified two protein domains which are predicted in both the wheat and rice PSTOL genes which include a protein kinase ATP-binding region signature (PS00107) located from amino acids 9–31 and a Serine/Threonine protein kinase active-site signature (PS00108) located from amino acids 126–138. The second best BLASTp hit to OsPSTOL in the wheat genome is encoded by a gene located on chromosome 3B, which shared only 54% amino acid identity suggesting that three sub-genome specific homologues were not present in wheat. A putative phosphate starvation activation domain thought to be bound by PHR1, also known as a P1BS binding site, was identified at − 1156 to − 1149 bp in the predicted promoter region of TaPSTOL . To provide additional evidence that there is only a single gene present in wheat, the TaPSTOL sequence was also used to conduct a BLASTn search against Triticum urartu (A genome), Aegilops speltoides (B genome) and Aegilops tauschii (D genome) sequence databases. A single near perfect match was found for T.urartu. However no strong match was found in either the A.speltoides or the A. tauschii genome sequences. The best match for A. speltoides was 929 bp but with only 73% identity and for A. tauschii the best match was 266 bp with 85% identity [24, 27].
Due to the incomplete nature of the wheat genome assembly, we confirmed the absence of TaPSTOL homoeologues on chromosomes 5B and 5D by screening the Chinese Spring nulli-tetrasomic deletion lines . PCRs using three primer combinations were designed using the 5A TaPSTOL promoter and coding region. This revealed only one copy of the PSTOL gene on chromosome 5A and no homoeologues on 5B or 5D (Additional file 2: Figure S1). To allow for potential small sequence variations in primer hybridisation sites, the stringency of primer annealing was reduced to 55 °C, but no additional amplicons were amplified. This suggests that only a single copy of TaPSTOL exists in Chinese Spring and that it is located on the short arm of chromosome 5A.
To search for sequence variation which might exist in the TaPSTOL gene or regulatory regions, we re-sequenced TaPSTOL from different wheat accessions: (1) the 22 founders of two bread wheat multiparent advanced generation inter-cross (MAGIC) populations (NIAB Elite MAGIC, ; NIAB Diverse MAGIC, ) that collectively capture 91% of genetic diversity in UK wheat; (2) spring varieties Fielder (USA) and Chinese Spring (China); and (3) three T. durum accessions (Additional file 5: Table S2). No sequence differences in any of the bread wheat lines tested were detected. A small variant was found in T. durum accession PI503314 which contained two indels in the promoter region. The first was a 13 bp insertion at − 355 and the second was a 34 bp deletion at − 704 bp. The A-genome sequence from other two dicoccoides accessions was identical to the hexaploid wheat sequence. No differences were found in the coding region in any of the lines tested.
TaPSTOL gene expression
To further understand if Traes_5AS_AA3DC6A5F is the putative wheat PSTOL homologue, Fielder wheat plants were grown with varying levels of P (1, 5,10, 20 and 45 μM P) in hydroponics for seven days and expression of TaPSTOL was measured in both roots and shoots (Fig. 1). TaPSTOL is expressed relatively evenly in both root and shoot tissues, however TaPSTOL transcript abundance decreases with increasing P concentration in both tissue types. To further examine the expression of TaPSTOL, a 2.4 kb region directly upstream of the start codon was cloned and used to drive the expression of a GUS reporter gene. Characterization of 10 transformed wheat lines showed very low to no expression in roots, shoots or flowers of plants when grown under standard conditions in compost. When plants were grown under nutrient limiting conditions, TaPSTOL expression was found in the roots and shoots, with the highest expression in the root tips of both primary and lateral roots and coleoptiles (Fig. 2). Expression could also be seen in both the leaf trichomes and root hairs. Similar patterns of expression where observed in all ten lines tested.
Characterization of TaPSTOL overexpression and RNAi knockdown lines in wheat
In order to evaluate the function of TaPSTOL, and compare to that of known OsPSTOL phenotypes, both over expression (OE) and RNA interference (RNAi) wheat transgenic lines were created. The expression level of TaPSTOL in thirty three OE and forty RNAi independent T0 lines was measured. Lines were selected for further study with a high over-expression (30× higher than wild type), low level overexpression (3X higher than wild type), a highly knocked down line (90% knock down) and a lower level of knock down (80% knock down), line here to referred as OE-1, OE-2, RNAi-1 and RNAi-2 respectively. Homozygous T2 plants were selected on the basis of nptII qPCR copy number analysis and resulting lines were phenotypically assessed under both low P conditions (3 μM P) in sand to better control P levels and standard P growth conditions, (M2 compost with 109 mg/l available P).
For OE and RNAi lines grown on sand (fed twice a week with Magnavaca solution containing low P at 3 μM) a significant difference in the dry weight of roots from both the highest modified expression levels of OE and RNAi plants could be seen when grown to seed (ie OE-1 and RNAi-1, Fig. 3a). No significant differences in dry weight were found for the shoots of the OE or RNAi transgenic plants (Fig. 3b and e).
No difference in tiller number under low P conditions was seen in either the OE or RNAi lines (Fig. 3c). However, when grown on M2 compost, a significant difference in the tiller number was observed, with RNAi and OE lines showing a significantly lower and significantly higher tiller number, respectively (p val < 0.05) (Fig. 3f). Despite the difference in tiller number, modification of TaPSTOL expression showed no significant effect on total yield between the lines tested with growth in compost.
During these experiments a novel phenotype was observed under low P conditions, namely a difference in the time to flowering, defined by time until growth stage 49 in days . The RNAi lines were significantly earlier in flowering relative to the WT plants by almost two days (p val < 0.05)(Fig. 4a). As anticipated, the converse phenotype was seen for the OE lines as the time taken to flowering was significantly longer than the WT plants (p val. < 0.05). This difference in flowering time correlated with a significant negative effect on the total yield under low P conditions (3 μM) but not under M2 grown conditions for lines overexpressing TaPSTOL (Fig. 3d and g). The effect on flowering time was less pronounced under M2 compost growth conditions as only the highest overexpression line (OE-1) and highest knockdown line (RNAi-1) showed significant variation relative to WT Fielder (Fig. 4b).
A change in the overall root biomass production with a significant difference in root DW under low P for the both the highest (OE-1) and lowest (RNAi-1) TaPSTOL transgenic lines was observed (Fig. 3). While the overall total DW (root DW + shoot DW) under low P conditions was unaffected, there was an observable shift in DW from the root to shoot.
TaPSTOL transgenic lines also showed differences in seed size. OE lines grown under P stress produced a decreased number of seeds compared to both RNAi and WT plants. In contrast, while OE lines produced far fewer seeds, each seed was 40–50% larger than the RNAi seeds and 30–40% larger than the WT seeds (Fig. 5b and d). This effect was mainly observed by an increase in the width of the seeds (Fig. 6) while the length remained similar.
Total P concentration in the roots, shoot and seeds was also measured in five plants chosen randomly for each transgenic line grown under low P conditions and for shoots and seeds grown under compost conditions. No significant differences in P concentration were seen between root or shoot tissues grown under either condition in any of the sampled lines (Fig. 7a and b). There was a significant difference in the amount of P in the seeds of the RNAi-1 and OE-1 and OE-2 relative to WT (Fig. 7c). The differences in P concentration showed a lower P concentration in the RNAi line and higher Pi concentration in the OE lines (p val < 0.05). Other essential elements including K, Ca, Mg, Mn, Zn, Cu, Fe, Mo, S, and Co were also tested for significant differences between lines and no significant differences were found.
The role of TaPSTOL in phosphorus use efficiency:
To understand how modification of TaPSTOL expression changes the PUE of wheat. The Physiological P use efficiency (PPUE) and the P efficiency ratio (PER) were measured (Fig. 8). There was a significant difference in the PPUE and PER when plants were grown on limiting P (p val < 0.001). The OE lines showed a significantly lower PPUE and PER, probably due to a significant increase in P in the grain accompanied by a decrease in yield. These differences in PPUE and PER were not observed under compost growth conditions, as no difference in the P concentration in the grain or differences in yield were seen.
The need to drive efficiency in global agricultural production has led to the elucidation of a number of key genes in breeding for phosphate efficient crops. Here we set out to identify and characterize the wheat homologue of OsPSTOL, which we denote as TaPSTOL, and compare the effect of manipulating its expression on a number of agronomic traits both under compost and low P conditions. The rice gene OsPSTOL has been shown to confer phenotypes beneficial to growth on P limiting soils, including increased grain yield, root and shoot biomass and increased tillering [14, 15, 18, 47]. Our data indicates that many of these PSTOL associated phenotypes are conserved between rice and wheat, including altered root growth under low P conditions (Fig. 3a), increased P content in the grain (Fig. 7c), and increased tiller number (Fig. 3f). There are, however, subtle differences as to when these phenotypes are exhibited in wheat as compared to rice; for instance, the increased tillering phenotype in wheat was only observed in plants grown in compost, whereas in rice, it was also found under P limiting conditions . These differences in the conditions under which the increase in tillering occurs, might be a result of the stringent P limiting conditions in which the plants were grown in sand or the compost system used here to represent a replete P substrate. They could also be due to inherent differences between the two crop species, or the genetic background in which the PSTOL gene was studied in both rice and wheat. These differences are worthy of further study to identify the underlying mechanism of the PSTOL gene in both species and how TaPSTOL contributes to the regulation of tillering.
While many of the observed rice OsPSTOL mediated phenotypes were mirrored in wheat, further supporting the conserved nature of the PSTOL gene, changes in biomass and yield were not observed. This might be due the choice of germplasm used to test the effect of P efficiency in wheat or the P levels chosen might be too limiting or too generous to see the effects described in rice. There is currently no published data on P efficiency in the variety Fielder, but other work has suggested that natural variation in P efficiency does exist in wheat and this may alter the effect of TaPSTOL on P efficiency  within a specific variety. Recent findings in rice suggest that plants which contain the PSTOL gene are not the most efficient under low P conditions, which suggests that other genetic loci may play a greater role in increasing P efficiency . Further work in both rice and wheat is required to further understand the processes by which P efficiency is increased, and the role played by PSTOL in these differences. While we found PSTOL mediated phenotypic effects in common between rice and wheat, our findings also highlight that the need for more understanding of how PSTOL genes, and their interactions with other genetic loci, may make crops more efficient under P limiting conditions.
We also identified two phenotypes associated with the TaPSTOL locus which were not observed in rice, showing that TaPSTOL modulates wheat flowering time and seed size. The timing of flowering is of great agronomic importance, as it entrains fertilization and seed set with favourable environmental conditions, thus helping farmers reach full yield potential in target agricultural environments . Seed size is also an important component to driving yield gains, and we can now add another locus to those previously described as factors influencing grain size in wheat [51, 52]. It is interesting that the lack of genetic diversity at TaPSTOL in most of the wheat lines tested suggests that this gene may not have been found by traditional QTL mapping in most UK varieties, as no variation exists to be exploited in a mapping population. The lack of variation seen in the wheat cultivars tested, point to this region as having been under strong selection pressure. The majority of these bread wheat varieties originated from the UK, however additional varieties from Europe, USA and China were also included, yet did not show any variation at this locus. Only when the search is extended to the distantly related dicoccoides accessions, which originated from Israel, is some limited variation observed within TaPSTOL, limited to the promoter region. This finding is in contrast to OsPSTOL in rice where natural variation is found across species from a number of different continents in both the promoter and coding regions [17, 19]. The lack of DNA sequence variation might also be a consequence of a single PSTOL homoeologue present in the hexaploid wheat genome, located on chromosome 5A. Estimates for the number of genes present as just a single homoeologue in the hexaploid wheat genome are estimated to be in the range 3 to11% [53, 54].
Finally alteration of the expression of TaPSTOL led to changes in PUE as a decrease in yield and increase in P concentration when grown under P limiting conditions led to lower PUE scores for the over expression of TaPSTOL but not for the RNAi lines. As a positive increase in yield did not occur with the over expression of TaPSTOL to compensate for the increased P found in the grain, the PUE score was significantly lowered. This was not the case for the OsPSTOL where both an increase in yield along with an increase in P concentration of the tissues were found . This difference in rice versus wheat might be due to some other part of the signalling cascade needing to be identified and combined with increased PSTOL expression to help drive the yield increases under limiting P growth conditions.
We have identified and characterized the PSTOL gene in wheat and demonstrated that it controls a number of agronomic characteristics important to breeders and farmers, including altering root growth under P deficiency conditions and increased P content in the grain. While some of these phenotypes are conserved with PSTOL genes in other plant species, other phenotypes manifest themselves under slightly different P conditions. In addition, we show that while PSTOL may not influence yield as it does in rice, TaPSTOL modulates flowering time and seed size in wheat. These results demonstrate PSTOL genes in multiple crop species may be important targets for improving agronomic performance, both under P limiting and non P limiting conditions.
Days post anthesis
Expressed sequence tag
Inductively coupled plasma-mass spectrometry
Multiparent advanced generation inter-cross
P efficiency ratio
- Pi :
Physiological P use efficiency
Phosphate starvation tolerance
Phosphate use efficiency
Nopaline synthase terminator
Kirkby EA, Johnston A. The ecophysiology of plant-phosphorus interactions. Netherlands: Springer; 2008.
Hinsinger P. Bioavailability of soil inorganic P in the rhizosphere as affected by root-induced chemical changes: a review. Plant Soil. 2001;237:173–95.
Cakmak I. Plant nutrition research: Priorities to meet human needs for food in sustainable ways. Prog Plant Nutr Plenary Lect XIV Int Plant Nutr Colloq. Dordrecht: Springer Netherlands; 2002. p. 3–24.
Vance CP, Uhde-Stone C, Allan DL. Phosphorus acquisition and use: critical adaptations by plants for securing a nonrenewable resource. New Phytol. 2003;157:423–47.
Lynch JP, Brown KM. Topsoil foraging – an architectural adaptation of plants to low phosphorus availability. Plant Soil. 2001;237:225–37.
Shen J, Yuan L, Zhang J, Li H, Bai Z, Chen X, Zhang W, Zhang F. Phosphorus dynamics: from soil to plant. Plant Physiol. 2011;156(3):997–1005.
Herrera-Estrella L, López-Arredondo D. Phosphorus: the underrated element for feeding the world. Trends Plant Sci. 2016;21:461–3.
FAOSTAT [Internet]. [cited 2017 Dec 15]. Available from: http://www.fao.org/faostat/en/#home
López-Arredondo DL, Leyva-González MA, González-Morales SI, López-Bucio J, Herrera-Estrella L. Phosphate nutrition: improving low-phosphate tolerance in crops. Annu Rev Plant Biol. 2014;65:95–123.
Ozturk L, Eker S, Torun B, Cakmak I. Variation in phosphorus efficiency among 73 bread and durum wheat genotypes grown in a phosphorus-deficient calcareous soil. Plant Soil. 2005;269:69–80.
Malhi SS, Vera CL, Brandt SA. Seed yield potential of five wheat species/cultivars without and with phosphorus fertilizer application on a P-deficient soil in northeastern Saskatchewan. Agric Sci. 2015;6(02):224.
Zhang Z, Liao H, Lucas WJ. Molecular mechanisms underlying phosphate sensing, signaling, and adaptation in plants. J Integr Plant Biol. 2014;56:192–220.
Ni JJ, Wu P, Senadhira D, Huang N. Mapping QTLs for phosphorus deficiency tolerance in rice ( Oryza sativa L.). TAG Theor Appl Genet. 1998;97:1361–9.
Wissuwa M, Yano M, Ae N. Mapping of QTLs for phosphorus-deficiency tolerance in rice ( Oryza sativa L.). TAG Theor Appl Genet. 1998;97:777–83.
Wissuwa M, Wegner J, Ae N, Yano M. Substitution mapping of Pup1 : a major QTL increasing phosphorus uptake of rice from a phosphorus-deficient soil. TAG Theor Appl Genet. 2002;105:890–7.
Heuer S, Lu X, Chin JH, Tanaka JP, Kanamori H, Matsumoto T, et al. Comparative sequence analyses of the major quantitative trait locus p hosphorus up take 1 ( Pup1 ) reveal a complex genetic structure. Plant Biotechnol J. 2009;7:456–71.
Pariasca-Tanaka J, Chin JH, Dramé KN, Dalid C, Heuer S, Wissuwa M. A novel allele of the P-starvation tolerance gene OsPSTOL1 from African rice (Oryza glaberrima Steud) and its distribution in the genus Oryza. Theor Appl Genet. 2014;127:1387–98.
Gamuyao R, Chin JH, Pariasca-Tanaka J, Pesaresi P, Catausan S, Dalid C, et al. The protein kinase Pstol1 from traditional rice confers tolerance of phosphorus deficiency. Nature. 2012;488:535–9.
Vigueira CC, Small LL, Olsen KM. Long-term balancing selection at the phosphorus starvation tolerance 1 (PSTOL1) locus in wild, domesticated and weedy rice (Oryza). BMC Plant Biol. 2016;16:101.
Wissuwa M. Combining a modelling with a genetic approach in establishing associations between genetic and physiological effects in relation to phosphorus uptake. Plant Soil. 2005;269:57–68.
Hufnagel B, de Sousa SM, Assis L, Guimaraes CT, Leiser W, Azevedo GC, et al. Duplicate and conquer: multiple homologs of PHOSPHORUS-STARVATION TOLERANCE1 enhance phosphorus acquisition and Sorghum performance on low-phosphorus soils. Plant Physiol. 2014;166:659–77.
Azevedo GC, Cheavegatti-Gianotto A, Negri BF, Hufnagel B, e Silva L da C, Magalhaes J V, et al. Multiple interval QTL mapping and searching for PSTOL1 homologs associated with root morphology, biomass accumulation and phosphorus content in maize seedlings under low-P. BMC Plant Biol 2015;15:172.
The European Bioinformatics Institute < EMBL-EBI [Internet]. [cited 2012 Sep 1]. Available from: https://www.ebi.ac.uk/
Ensembl Plants [Internet]. [cited 2012 Sep 1]. Available from: http://plants.ensembl.org/index.html
BLAST: Basic Local Alignment Search Tool [Internet]. [cited 2012 Sep 1]. Available from: https://blast.ncbi.nlm.nih.gov/Blast.cgi
National BioResource Project (NBRP):KOMUGI -Integrated Wheat Science Database. - [Internet]. [cited 2012 Sep 1]. Available from: https://shigen.nig.ac.jp/wheat/komugi/
URGI BLAST [Internet]. [cited 2015 May 1]. Available from: https://urgi.versailles.inra.fr/blast/blast.php
Softberry Home Page [Internet]. [cited 2012 Sep 1]. Available from: http://www.softberry.com/
InterPro protein sequence analysis & classification < InterPro < EMBL-EBI [Internet]. [cited 2012 Sep 1]. Available from: https://www.ebi.ac.uk/interpro/
Rubio V, Linhares F, Solano R, Martín AC, Iglesias J, Leyva A, et al. A conserved MYB transcription factor involved in phosphate starvation signaling both in vascular plants and in unicellular algae. Genes Dev. 2001;15:2122–33.
Nakagawa S, Niimura Y, Gojobori T, Tanaka H, Miura K. Diversity of preferred nucleotide sequences around the translation initiation codon in eukaryote genomes. Nucleic Acids Res. 2008;36:861–71.
Perochon A, Jianguang J, Kahla A, Arunachalam C, Scofield SR, Bowden S, et al. TaFROG encodes a Pooideae orphan protein that interacts with SnRK1 and enhances resistance to the mycotoxigenic fungus fusarium graminearum. Plant Physiol. 2015;169:01056.2015.
McElroy D, Zhang W, Cao J, Wu R. Isolation of an efficient actin promoter for use in rice transformation. Plant Cell. 1990;2:163–71.
Bates R, Craze M, Wallington EJ. Agrobacterium -mediated transformation of oilseed rape (Brassica napus). Curr Protoc Plant Biol. 2017;2:287–298. Hoboken; Wiley. https://doi.org/10.1002/cppb.20060.
Ishida Y, Tsunashima M, Hiei Y, Komari T. Wheat (Triticum aestivum L.) transformation using immature embryos. In Agrobacterium Protocols. New York: Springer; 2015. p. 189–198.
Risacher T, Craze M, Bowden S, Paul W, Barsby T. Highly efficient agrobacterium-mediated transformation of wheat via in planta inoculation. In Transgenic Wheat, Barley and Oats. Humana Press; 2009. p. 115–124.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods. 2001;25:402–8.
Magnavaca R, Gardner CO, Clark RB. Evaluation of inbred maize lines for aluminum tolerance in nutrient solution. Genet Asp Plant Miner Nutr. Dordrecht: Springer Netherlands; 1987. p. 255–65.
Van Riet L, Nagaraj V, Van den Ende W, Clerens S, Wiemken A, Van Laere A. Purification, cloning and functional characterization of a fructan 6-exohydrolase from wheat (Triticum aestivum L.). J Exp Bot. 2006;57:213–23.
Xia N, Zhang G, Liu X-Y, Deng L, Cai G-L, Zhang Y, et al. Characterization of a novel wheat NAC transcription factor gene involved in defense response against stripe rust pathogen infection and abiotic stresses. Mol Biol Rep. 2010;37:3703–12.
Fulton TM, Chunwongse J, Tanksley SD. Microprep protocol for extraction of DNA from tomato and other herbaceous plants. Plant Mol Biol Report. 1995;13:207–9.
Stomp AM. Histochemical localization of b-glucuronidase. GUS protocols: using the GUS gene as a reporter of gene expression. 1992. p. 103–113.
Thomas CL, Alcock TD, Graham NS, Hayden R, Matterson S, Wilson L, Young SD, Dupuy LX, White PJ, Hammond JP, Danku JMC, Salt DE, Sweeney A, Bancroft I, Broadley MR. BMC Plant Biol. 2016;16:214.
Hammond JP, Broadley MR, White PJ, King GJ, Bowen HC, Hayden R, et al. Shoot yield drives phosphorus use efficiency in Brassica oleracea and correlates with root architecture traits. J Exp Bot. 2009;60:1953–68.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.
Sears ER. The aneuploids of common wheat. Univ Missouri, Coll Agric Agric Exp Stn Bull. 1954;572:1–59.
Mackay IJ, Bansept-Basler P, Barber T, Bentley AR, Cockram J, Gosman N, et al. An eight-parent multiparent advanced generation inter-cross population for winter-sown wheat: creation, properties, and validation. G3. 2014;4:1603–10.
White J, Law JR, MacKay I, Chalmers KJ, Smith JSC, Kilian A, et al. The genetic diversity of UK, US and Australian cultivars of Triticum aestivum measured by DArT markers and considered by genome. Theor Appl Genet. 2008;116:439–53.
Zadoks JC, Chang TT, Konzak CF. A decimal code for the growth stages of cereals. Weed Res. 1974;14:415–21.
Bentley AR, Horsnell R, Werner CP, Turner AS, Rose GA, Bedard C, et al. Short, natural, and extended photoperiod response in BC2F4 lines of bread wheat with different Photoperiod-1 (Ppd-1) alleles. J Exp Bot. 2013;64:1783–93.
Simmonds J, Scott P, Brinton J, Mestre TC, Bush M, del Blanco A, et al. A splice acceptor site mutation in TaGW2-A1 increases thousand grain weight in tetraploid and hexaploid wheat through wider and longer grains. Theor Appl Genet. 2016;129:1099–112.
Ma M, Zhao H, Li Z, Hu S, Song W, Liu X. TaCYP78A5 regulates seed size in wheat ( Triticum aestivum ). J Exp Bot. 2016;67:1397–410.
Akhunov ED, Akhunova AR, Linkiewicz AM, Dubcovsky J, Hummel D, Lazo G, et al. Synteny perturbations between wheat homoeologous chromosomes caused by locus duplications and deletions correlate with recombination rates. Proc Natl Acad Sci U S A. 2003;100:10836–41.
Qi LL, Echalier B, Chao S, Lazo GR, Butler GE, Anderson OD, et al. A chromosome bin map of 16,000 expressed sequence tag loci and distribution of genes among the three genomes of polyploid wheat. Genetics. 2004;168:701–12.
The authors would like to thank Olga Murshudova for additional tissue culture media support, Richard Horsnell for gDNA, Ian MacKay for information on the NIAB MAGIC populations, and Andy Greenland for his support and advice.
This work was supported by the NIAB Trust. The funders had no role in experimental design, data collection and analysis or preparation of the manuscript.
Availability of data and materials
Sequence data from this article has been deposited to NCBI, with the accession numbers MH043198 and MH043199. Constructs and plant materials are available under MTA from the corresponding author.
Ethics approval and consent to participate
The authors declare they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Primers used in this study. (DOCX 33 kb)
Figure S1. Location of gene-specific primers and PCR amplification of TaPSTOL genomic regions on DNA extracted from nullisomic (N) / tetrasomic (T) wheat lines. The PCR amplicons correspond to a promoter region (PCR1, 1458 bp), a promoter and coding region (PCR2, 553 bp), and a coding region (PCR3, 870 bp), to demonstrate that TaPSTOL is only present on chromosome 5A. (DOCX 63 kb)
Figure S2. T-DNA structure for the three TaPSTOL constructs used in this study. (DOCX 30 kb) (DOCX 80 kb)
Figure S3. ClustalW alignment of the rice and wheat PSTOL predicted proteins. Conserved amino acids are coloured in red. (DOCX 162 kb)
Table S2. Geographical origin of wheat accessions used in this study. (DOCX 30 kb)
About this article
Cite this article
Milner, M.J., Howells, R.M., Craze, M. et al. A PSTOL-like gene, TaPSTOL, controls a number of agronomically important traits in wheat. BMC Plant Biol 18, 115 (2018). https://doi.org/10.1186/s12870-018-1331-4
- Seed size
- Seed number
- Flowering time