Nutritional and tissue-specific regulation of cytochrome P450 CYP711A MAX1 homologues and strigolactone biosynthesis in wheat.

Strigolactones (SLs) are a class of phytohormones regulating branching/tillering, and their biosynthesis has been associated with nutritional signals and plant adaptation to nutrient-limiting conditions. The enzymes in the SL biosynthetic pathway downstream of carlactone are of interest as they are responsible for structural diversity in SLs, particularly cytochrome P450 CYP711A subfamily members, such as MORE AXILLARY GROWTH1 (MAX1) in Arabidopsis. We identified 13 MAX1 homologues in wheat, clustering in four clades and five homoeologous subgroups. The utilization of RNA-sequencing data revealed a distinct expression pattern of MAX1 homologues in above- and below-ground tissues, providing insights into the distinct roles of MAX1 homologues in wheat. In addition, a transcriptional analysis showed that SL biosynthetic genes were systematically regulated by nitrogen supply. Nitrogen limitation led to larger transcriptional changes in the basal nodes than phosphorus limitation, which was consistent with the observed tillering suppression, as wheat showed higher sensitivity to nitrogen. The opposite was observed in roots, with phosphorus limitation leading to stronger induction of most SL biosynthetic genes compared with nitrogen limitation. The observed tissue-specific regulation of SL biosynthetic genes in response to nutritional signals is likely to reflect the dual role of SLs as rhizosphere signals and branching inhibitors.

The Functional Diversity of the High-Affinity Nitrate Transporter Gene Family in Hexaploid Wheat: Insights from Distinct Expression Profiles.

High-affinity nitrate transporters (NRT) are key components for nitrogen (N) acquisition and distribution within plants. However, insights on these transporters in wheat are scarce. This study presents a comprehensive analysis of the NRT2 and NRT3 gene families, where the aim is to shed light on their functionality and to evaluate their responses to N availability. A total of 53 NRT2s and 11 NRT3s were identified in the bread wheat genome, and these were grouped into different clades and homoeologous subgroups. The transcriptional dynamics of the identified NRT2 and NRT3 genes, in response to N starvation and nitrate resupply, were examined by RT-qPCR in the roots and shoots of hydroponically grown wheat plants through a time course experiment. Additionally, the spatial expression patterns of these genes were explored within the plant. The NRT2s of clade 1, TaNRT2.1-2.6, showed a root-specific expression and significant upregulation in response to N starvation, thus emphasizing a role in N acquisition. However, most of the clade 2 NRT2s displayed reduced expression under N-starved conditions. Nitrate resupply after N starvation revealed rapid responsiveness in TaNRT2.1-2.6, while clade 2 genes exhibited gradual induction, primarily in the roots. TaNRT2.18 was highly expressed in above-ground tissues and exhibited distinct nitrate-related response patterns for roots and shoots. The TaNRT3 gene expression closely paralleled the profiles of TaNRT2.1-2.6 in response to nitrate induction. These findings enhance the understanding of NRT2 and NRT3 involvement in nitrogen uptake and utilization, and they could have practical implications for improving nitrogen use efficiency. The study also recommends a standardized nomenclature for wheat NRT2 genes, thereby addressing prior naming inconsistencies.

Identification of traits underpinning good breadmaking performance of wheat grown with reduced nitrogen fertilisation.

BACKGROUND: Nitrogen fertiliser is the major input and cost for wheat production, being required to support the development of the canopy to maximise yield and for the synthesis of the gluten proteins that are necessary for breadmaking. Consequently, current high-yielding cultivars require the use of nitrogen fertilisation levels above the yield optimum to achieve the grain protein content needed for breadmaking. This study aimed to reduce this requirement by identifying traits that allow the use of lower levels of nitrogen fertiliser to produce wheat for breadmaking. RESULTS: A range of commercial wheat genotypes (cultivars) were grown in multiple field trials (six sites over 3 years) in the UK with optimal (200 kg Ha-1 ) and suboptimal (150 kg Ha-1 ) application of nitrogen. Bulked grain samples from four sites per year were milled and white flours were baked using three types of breadmaking process. This identified five cultivars that consistently exhibited good breadmaking quality when grown with the lower nitrogen application. Chemical and biochemical analyses showed that the five cultivars were characterised by exhibiting grain protein deviation (GPD) and high dough elasticity. CONCLUSIONS: It is possible to develop novel types of wheat that exhibit good breadmaking quality by selecting for GPD and high dough strength. © 2023 The Authors. Journal of The Science of Food and Agriculture published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Time series canopy phenotyping enables the identification of genetic variants controlling dynamic phenotypes in soybean.

Advances in plant phenotyping technologies are dramatically reducing the marginal costs of collecting multiple phenotypic measurements across several time points. Yet, most current approaches and best statistical practices implemented to link genetic and phenotypic variation in plants have been developed in an era of single-time-point data. Here, we used time-series phenotypic data collected with an unmanned aircraft system for a large panel of soybean (Glycine max (L.) Merr.) varieties to identify previously uncharacterized loci. Specifically, we focused on the dissection of canopy coverage (CC) variation from this rich data set. We also inferred the speed of canopy closure, an additional dimension of CC, from the time-series data, as it may represent an important trait for weed control. Genome-wide association studies (GWASs) identified 35 loci exhibiting dynamic associations with CC across developmental stages. The time-series data enabled the identification of 10 known flowering time and plant height quantitative trait loci (QTLs) detected in previous studies of adult plants and the identification of novel QTLs influencing CC. These novel QTLs were disproportionately likely to act earlier in development, which may explain why they were missed in previous single-time-point studies. Moreover, this time-series data set contributed to the high accuracy of the GWASs, which we evaluated by permutation tests, as evidenced by the repeated identification of loci across multiple time points. Two novel loci showed evidence of adaptive selection during domestication, with different genotypes/haplotypes favored in different geographic regions. In summary, the time-series data, with soybean CC as an example, improved the accuracy and statistical power to dissect the genetic basis of traits and offered a promising opportunity for crop breeding with quantitative growth curves.

Hidden variation in polyploid wheat drives local adaptation.

Wheat has been domesticated into a large number of agricultural environments and has the ability to adapt to diverse environments. To understand this process, we survey genotype, repeat content, and DNA methylation across a bread wheat landrace collection representing global genetic diversity. We identify independent variation in methylation, genotype, and transposon copy number. We show that these, so far unexploited, sources of variation have had a significant impact on the wheat genome and that ancestral methylation states become preferentially "hard coded" as single nucleotide polymorphisms (SNPs) via 5-methylcytosine deamination. These mechanisms also drive local adaption, impacting important traits such as heading date and salt tolerance. Methylation and transposon diversity could therefore be used alongside SNP-based markers for breeding.

The interaction between wheat roots and soil pores in structured field soil.

Wheat (Triticum aestivum L.) root growth in the subsoil is usually constrained by soil strength, although roots can use macropores to elongate to deeper layers. The quantitative relationship between the elongation of wheat roots and the soil pore system, however, is still to be determined. We studied the depth distribution of roots of six wheat varieties and explored their relationship with soil macroporosity from samples with the field structure preserved. Undisturbed soil cores (to a depth of 100 cm) were collected from the field and then non-destructively imaged using X-ray computed tomography (at a spatial resolution of 90 µm) to quantify soil macropore structure and root number density (the number of roots cm-2 within a horizontal cross-section of a soil core). Soil macroporosity changed significantly with depth but not between the different wheat lines. There was no significant difference in root number density between wheat varieties. In the subsoil, wheat roots used macropores, especially biopores (i.e. former root or earthworm channels) to grow into deeper layers. Soil macroporosity explained 59% of the variance in root number density. Our data suggested that the development of the wheat root system in the field was more affected by the soil macropore system than by genotype. On this basis, management practices which enhance the porosity of the subsoil may therefore be an effective strategy to improve deep rooting of wheat.

Post-translational cleavage of HMW-GS Dy10 allele improves the cookie-making quality in common wheat (Triticum aestivum).

Wheat is a major staple food crop worldwide because of the unique properties of wheat flour. High molecular weight glutenin subunits (HMW-GSs), which are among the most critical determinants of wheat flour quality, are responsible for the formation of glutenin polymeric structures via interchain disulfide bonds. We herein describe the identification of a new HMW-GS Dy10 allele (Dy10-m619SN). The amino acid substitution (serine-to-asparagine) encoded in this allele resulted in a partial post-translational cleavage that produced two new peptides. These new peptides disrupted the interactions among gluten proteins because of the associated changes to the number of available cysteine residues for interchain disulfide bonds. Consequently, Dy10-m619SN expression decreased the size of glutenin polymers and weakened glutens, which resulted in wheat dough with improved cookie-making quality, without changes to the glutenin-to-gliadin ratio. In this study, we clarified the post-translational processing of HMW-GSs and revealed a new genetic resource useful for wheat breeding. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-021-01238-9.

Soil strength influences wheat root interactions with soil macropores.

Deep rooting is critical for access to water and nutrients found in subsoil. However, damage to soil structure and the natural increase in soil strength with depth, often impedes root penetration. Evidence suggests that roots use macropores (soil cavities greater than 75 µm) to bypass strong soil layers. If roots have to exploit structures, a key trait conferring deep rooting will be the ability to locate existing pore networks; a trait called trematotropism. In this study, artificial macropores were created in repacked soil columns at bulk densities of 1.6 g cm-3 and 1.2 g cm-3 , representing compact and loose soil. Near isogenic lines of wheat, Rht-B1a and Rht-B1c, were planted and root-macropore interactions were visualized and quantified using X-ray computed tomography. In compact soil, 68.8% of root-macropore interactions resulted in pore colonization, compared with 12.5% in loose soil. Changes in root growth trajectory following pore interaction were also quantified, with 21.0% of roots changing direction (±3°) in loose soil compared with 76.0% in compact soil. These results indicate that colonization of macropores is an important strategy of wheat roots in compacted subsoil. Management practices to reduce subsoil compaction and encourage macropore formation could offer significant advantage in helping wheat roots penetrate deeper into subsoil.

Phylogeny and gene expression of the complete NITRATE TRANSPORTER 1/PEPTIDE TRANSPORTER FAMILY in Triticum aestivum.

NPF genes encode membrane transporters involved in the transport of a large variety of substrates including nitrate and peptides. The NPF gene family has been described for many plants, but the whole NPF gene family for wheat has not been completely identified. The release of the wheat reference genome has enabled the identification of the entire wheat NPF gene family. A systematic analysis of the whole wheat NPF gene family was performed, including responses of specific gene expression to development and nitrogen supply. A total of 331 NPF genes (113 homoeologous groups) have been identified in wheat. The chromosomal location of the NPF genes is unevenly distributed, with predominant occurrence in the long arms of the chromosomes. The phylogenetic analysis indicated that wheat NPF genes are closely clustered with Arabidopsis, Brachypodium, and rice orthologues, and subdivided into eight subfamilies. The expression profiles of wheat NPF genes were examined using RNA-seq data, and a subset of 44 NPF genes (homoeologous groups) with contrasting expression responses to nitrogen and/or development in different tissues were identified. The systematic identification of gene composition, chromosomal locations, evolutionary relationships, and expression profiles contributes to a better understanding of the roles of the wheat NPF genes and lays the foundation for further functional analysis in wheat.

Functional QTL mapping and genomic prediction of canopy height in wheat measured using a robotic field phenotyping platform.

Genetic studies increasingly rely on high-throughput phenotyping, but the resulting longitudinal data pose analytical challenges. We used canopy height data from an automated field phenotyping platform to compare several approaches to scanning for quantitative trait loci (QTLs) and performing genomic prediction in a wheat recombinant inbred line mapping population based on up to 26 sampled time points (TPs). We detected four persistent QTLs (i.e. expressed for most of the growing season), with both empirical and simulation analyses demonstrating superior statistical power of detecting such QTLs through functional mapping approaches compared with conventional individual TP analyses. In contrast, even very simple individual TP approaches (e.g. interval mapping) had superior detection power for transient QTLs (i.e. expressed during very short periods). Using spline-smoothed phenotypic data resulted in improved genomic predictive abilities (5-8% higher than individual TP prediction), while the effect of including significant QTLs in prediction models was relatively minor (<1-4% improvement). Finally, although QTL detection power and predictive ability generally increased with the number of TPs analysed, gains beyond five or 10 TPs chosen based on phenological information had little practical significance. These results will inform the development of an integrated, semi-automated analytical pipeline, which will be more broadly applicable to similar data sets in wheat and other crops.

Genetic variation and heritability of grain protein deviation in European wheat genotypes.

There is a well-established negative relationship between the yield and the concentration of protein in the mature wheat grain. However, some wheat genotypes consistently deviate from this relationship, a phenomenon known as Grain Protein Deviation (GPD). Positive GPD is therefore of considerable interest in relation to reducing the requirement for nitrogen fertilization for producing wheat for breadmaking. We have carried out two sets of field experiments on multiple sites in South East England. The first set comprised 11 field trials of 6 cultivars grown over three years (2008-2011) and the second comprised 9 field trials of 40 genotypes grown over two years (2015-2017) and 5 field trials of 30 genotypes grown in a single year (2017-2018). All trials comprised three replicate randomized plots of each genotype and nutrient regime. These studies showed strong genetic variation in GPD, which also differed in stability between genotypes, with cultivars bred in the UK generally having higher GPD and higher stability than those bred in other European countries. The heritability of GPD was estimated as 0.44, based on data from the field trials of 30 and 40 genotypes. The largest component contributing to the genetic variance was genotype (0.30), with a smaller contribution of the interaction between genotype and year/site (0.11) and a small (but statistically significant) contribution of nitrogen level. These studies suggest that selection for GPD is a viable target for breeders.

A breeding strategy targeting the secondary gene pool of bread wheat: introgression from a synthetic hexaploid wheat.

KEY MESSAGE: Introgressing one-eighth of synthetic hexaploid wheat genome through a double top-cross plus a two-phase selection is an effective strategy to develop high-yielding wheat varieties. The continued expansion of the world population and the likely onset of climate change combine to form a major crop breeding challenge. Genetic advances in most crop species to date have largely relied on recombination and reassortment within a relatively narrow gene pool. Here, we demonstrate an efficient wheat breeding strategy for improving yield potentials by introgression of multiple genomic regions of de novo synthesized wheat. The method relies on an initial double top-cross (DTC), in which one parent is synthetic hexaploid wheat (SHW), followed by a two-phase selection procedure. A genotypic analysis of three varieties (Shumai 580, Shumai 969 and Shumai 830) released from this program showed that each harbors a unique set of genomic regions inherited from the SHW parent. The first two varieties were generated from very small populations, whereas the third used a more conventional scale of selection since one of bread wheat parents was a pre-breeding material. The three varieties had remarkably enhanced yield potential compared to those developed by conventional breeding. A widely accepted consensus among crop breeders holds that introducing unadapted germplasm, such as landraces, as parents into a breeding program is a risky proposition, since the size of the breeding population required to overcome linkage drag becomes too daunting. However, the success of the proposed DTC strategy has demonstrated that novel variation harbored by SHWs can be accessed in a straightforward, effective manner. The strategy is in principle generalizable to any allopolyploid crop species where the identity of the progenitor species is known.

The role of ZIP transporters and group F bZIP transcription factors in the Zn-deficiency response of wheat (Triticum aestivum).

Understanding the molecular basis of zinc (Zn) uptake and transport in staple cereal crops is critical for improving both Zn content and tolerance to low-Zn soils. This study demonstrates the importance of group F bZIP transcription factors and ZIP transporters in responses to Zn deficiency in wheat (Triticum aestivum). Seven group F TabZIP genes and 14 ZIPs with homeologs were identified in hexaploid wheat. Promoter analysis revealed the presence of Zn-deficiency-response elements (ZDREs) in a number of the ZIPs. Functional complementation of the zrt1/zrt2 yeast mutant by TaZIP3, -6, -7, -9 and -13 supported an ability to transport Zn. Group F TabZIPs contain the group-defining cysteine-histidine-rich motifs, which are the predicted binding site of Zn2+ in the Zn-deficiency response. Conservation of these motifs varied between the TabZIPs suggesting that individual TabZIPs may have specific roles in the wheat Zn-homeostatic network. Increased expression in response to low Zn levels was observed for several of the wheat ZIPs and bZIPs; this varied temporally and spatially suggesting specific functions in the response mechanism. The ability of the group F TabZIPs to bind to specific ZDREs in the promoters of TaZIPs indicates a conserved mechanism in monocots and dicots in responding to Zn deficiency. In support of this, TabZIPF1-7DL and TabZIPF4-7AL afforded a strong level of rescue to the Arabidopsis hypersensitive bzip19 bzip23 double mutant under Zn deficiency. These results provide a greater understanding of Zn-homeostatic mechanisms in wheat, demonstrating an expanded repertoire of group F bZIP transcription factors, adding to the complexity of Zn homeostasis.

Temperature and nitrogen supply interact to determine protein distribution gradients in the wheat grain endosperm.

Gradients exist in the distribution of storage proteins in the wheat (Triticum aestivum) endosperm and determine the milling properties and protein recovery rate of the grain. A novel image analysis technique was developed to quantify both the gradients in protein concentration, and the size distribution of protein bodies within the endosperm of wheat plants grown under two different (20 or 28 °C) post-anthesis temperatures, and supplied with a nutrient solution with either high or low nitrogen content. Under all treatment combinations, protein concentration was greater in the endosperm cells closest to the aleurone layer and decreased towards the centre of the two lobes of the grain, i.e. a negative gradient. This was accompanied by a decrease in size of protein bodies from the outer to the inner endosperm layers in all but one of the treatments. Elevated post-anthesis temperature had the effect of increasing the magnitude of the negative gradients in both protein concentration and protein body size, whilst limiting nitrogen supply decreased the gradients.

System analysis of metabolism and the transcriptome in Arabidopsis thaliana roots reveals differential co-regulation upon iron, sulfur and potassium deficiency.

Deprivation of mineral nutrients causes significant retardation of plant growth. This retardation is associated with nutrient-specific and general stress-induced transcriptional responses. In this study, we adjusted the external supply of iron, potassium and sulfur to cause the same retardation of shoot growth. Nevertheless, limitation by individual nutrients resulted in specific morphological adaptations and distinct shifts within the root metabolite fingerprint. The metabolic shifts affected key metabolites of primary metabolism and the stress-related phytohormones, jasmonic, salicylic and abscisic acid. These phytohormone signatures contributed to specific nutrient deficiency-induced transcriptional regulation. Limitation by the micronutrient iron caused the strongest regulation and affected 18% of the root transcriptome. Only 130 genes were regulated by all nutrients. Specific co-regulation between the iron and sulfur metabolic routes upon iron or sulfur deficiency was observed. Interestingly, iron deficiency caused regulation of a different set of genes of the sulfur assimilation pathway compared with sulfur deficiency itself, which demonstrates the presence of specific signal-transduction systems for the cross-regulation of the pathways. Combined iron and sulfur starvation experiments demonstrated that a requirement for a specific nutrient can overrule this cross-regulation. The comparative metabolomics and transcriptomics approach used dissected general stress from nutrient-specific regulation in roots of Arabidopsis.

Genetic variation in traits for nitrogen use efficiency in wheat.

Crop nutrient and especially nitrogen use efficiency (NUE) is both an economically and an environmentally highly desirable trait. It has been estimated that only a third of nitrogen inputs to cereal crop worldwide are recovered in grain for consumption, resulting in a huge waste of resource with major negative impacts on the environment. Most measures of NUE in wheat and other cereals are based on field assessments of crop yields at given N inputs, performance responses to added N fertilizer, or by quantifying N fertilizer recovery rates. However, NUE is a complex trait comprising two key major components, N uptake and N utilization efficiency, both also complex traits in themselves, each involving many physiological processes and biochemical pathways. A deeper understanding of the processes involved in NUE has been a target of the UK Wheat Genetic Improvement Network project (http://www.wgin.org.uk/). This has enabled the breakdown of characteristics contributing to NUE and an assessment of the variation present in those characteristics, predominantly in modern cultivars; a total of 13 years of data have been obtained to date. Significant but limited variation suggests a requirement for broader germplasm screening such as older varieties, landraces, and wild relatives.

Genome-wide analysis of autophagy-associated genes in foxtail millet (Setaria italica L.) and characterization of the function of SiATG8a in conferring tolerance to nitrogen starvation in rice.

BACKGROUND: Autophagy is a cellular degradation process that is highly evolutionarily-conserved in yeast, plants, and animals. In plants, autophagy plays important roles in regulating intracellular degradation and recycling of amino acids in response to nutrient starvation, senescence, and other environmental stresses. Foxtail millet (Setaria italica) has strong resistance to stresses and has been proposed as an ideal material for use in the study of the physiological mechanisms of abiotic stress tolerance in plants. Although the genome sequence of foxtail millet (Setaria italica) is available, the characteristics and functions of abiotic stress-related genes remain largely unknown for this species. RESULTS: A total of 37 putative ATG (autophagy-associated genes) genes in the foxtail millet genome were identified. Gene duplication analysis revealed that both segmental and tandem duplication events have played significant roles in the expansion of the ATG gene family in foxtail millet. Comparative synteny mapping between the genomes of foxtail millet and rice suggested that the ATG genes in both species have common ancestors, as their ATG genes were primarily located in similar syntenic regions. Gene expression analysis revealed the induced expression of 31 SiATG genes by one or more phytohormone treatments, 26 SiATG genes by drought, salt and cold, 24 SiATG genes by darkness and 25 SiATG genes by nitrogen starvation. Results of qRT-PCR showing that among 37 SiATG genes, the expression level of SiATG8a was the highest after nitrogen starvation treatment 24 h, suggesting its potential role in tolerance to nutrient starvation. Moreover, the heterologous expression of SiATG8a in rice improved nitrogen starvation tolerance. Compared to wild type rice, the transgenic rice performed better and had higher aboveground total nitrogen content when the plants were grown under nitrogen starvation conditions. CONCLUSIONS: Our results deepen understanding about the characteristics and functions of ATG genes in foxtail millet and also identify promising new genetic resources that should be of use in future efforts to develop varieties of foxtail millet and other crop species that have resistance to nitrogen deficiency stress.

The dynamics of protein body formation in developing wheat grain.

Wheat is a major source of protein in the diets of humans and livestock but we know little about the mechanisms that determine the patterns of protein synthesis in the developing endosperm. We have used a combination of enrichment with (15) N glutamine and NanoSIMS imaging to establish that the substrate required for protein synthesis is transported radially from its point of entrance in the endosperm cavity across the starchy endosperm tissues, before becoming concentrated in the cells immediately below the aleurone layer. This transport occurs continuously during grain development but may be slower in the later stages. Although older starchy endosperm cells tend to contain larger protein deposits formed by the fusion of small protein bodies, small highly enriched protein bodies may also be present in the same cells. This shows a continuous process of protein body initiation, in both older and younger starchy endosperm cells and in all regions of the tissue. Immunolabeling with specific antibodies shows that the patterns of enrichment are not related to the contents of gluten proteins in the protein bodies. In addition to providing new information on the dynamics of protein deposition, the study demonstrates the wider utility of NanoSIMS and isotope labelling for studying complex developmental processes in plant tissues.

A novel approach to identify genes that determine grain protein deviation in cereals.

Grain yield and protein content were determined for six wheat cultivars grown over 3 years at multiple sites and at multiple nitrogen (N) fertilizer inputs. Although grain protein content was negatively correlated with yield, some grain samples had higher protein contents than expected based on their yields, a trait referred to as grain protein deviation (GPD). We used novel statistical approaches to identify gene transcripts significantly related to GPD across environments. The yield and protein content were initially adjusted for nitrogen fertilizer inputs and then adjusted for yield (to remove the negative correlation with protein content), resulting in a parameter termed corrected GPD. Significant genetic variation in corrected GPD was observed for six cultivars grown over a range of environmental conditions (a total of 584 samples). Gene transcript profiles were determined in a subset of 161 samples of developing grain to identify transcripts contributing to GPD. Principal component analysis (PCA), analysis of variance (ANOVA) and means of scores regression (MSR) were used to identify individual principal components (PCs) correlating with GPD alone. Scores of the selected PCs, which were significantly related to GPD and protein content but not to the yield and significantly affected by cultivar, were identified as reflecting a multivariate pattern of gene expression related to genetic variation in GPD. Transcripts with consistent variation along the selected PCs were identified by an approach hereby called one-block means of scores regression (one-block MSR).

Complex phylogeny and gene expression patterns of members of the NITRATE TRANSPORTER 1/PEPTIDE TRANSPORTER family (NPF) in wheat.

NPF (formerly referred to as low-affinity NRT1) and 'high-affinity' NRT2 nitrate transporter genes are involved in nitrate uptake by the root, and transport and distribution of nitrate within the plant. The NPF gene family consists of 53 members in Arabidopsis thaliana, however only 11 of these have been functionally characterized. Although homologous genes have been identified in genomes of different plant species including some cereals, there is little information available for wheat (Triticum aestivum). Sixteen genes were identified in wheat homologous to characterized Arabidopsis low-affinity nitrate transporter NPF genes, suggesting a complex wheat NPF gene family. The regulation of wheat NFP genes by plant N-status indicated involvement of these transporters in substrate transport in relation to N-metabolism. The complex expression pattern in relation to tissue specificity, nitrate availability and senescence may be associated with the complex growth patterns of wheat depending on sink/source demands, as well as remobilization during grain filling.