Genome-Wide Association Study and Identification of Candidate Genes for Nitrogen Use Efficiency in Barley (Hordeum vulgare L.)

Genome-Wide Comprehensive Identification and In Silico Characterization of Lectin Receptor-Like Kinase Gene Family in Barley (Hordeum vulgare L.)

Lectin receptor-like kinases (LecRLKs) are a significant subgroup of the receptor-like kinases (RLKs) protein family. They play crucial roles in plant growth, development, immune responses, signal transduction, and stress tolerance. However, the genome-wide identification and characterization of LecRLK genes and their regulatory elements have not been explored in a major cereal crop, barley (Hordeum vulgare L.). Therefore, in this study, integrated bioinformatics tools were used to identify and characterize the LecRLK gene family in barley. Based on the phylogenetic tree and domain organization, a total of 113 LecRLK genes were identified in the barley genome (referred to as HvlecRLK) corresponding to the LecRLK genes of Arabidopsis thaliana. These putative HvlecRLK genes were classified into three groups: 62 G-type LecRLKs, 1 C-type LecRLK, and 50 L-type LecRLKs. They were unevenly distributed across eight chromosomes, including one unknown chromosome, and were predominantly located in the plasma membrane (G-type HvlecRLK (96.8%), C-type HvlecRLK (100%), and L-type HvlecRLK (98%)). An analysis of motif composition and exon-intron configuration revealed remarkable homogeneity with the members of AtlecRLK. Notably, most of the HvlecRLKs (27 G-type, 43 L-type) have no intron, suggesting their rapid functionality. The Ka/Ks and syntenic analysis demonstrated that HvlecRLK gene pairs evolved through purifying selection and gene duplication was the major factor for the expansion of the HvlecRLK gene family. Exploration of gene ontology (GO) enrichment indicated that the identified HvlecRLK genes are associated with various cellular processes, metabolic pathways, defense mechanisms, kinase activity, catalytic activity, ion binding, and other essential pathways. The regulatory network analysis identified 29 transcription factor families (TFFs), with seven major TFFs including bZIP, C2H2, ERF, MIKC_MADS, MYB, NAC, and WRKY participating in the regulation of HvlecRLK gene functions. Most notably, eight TFFs were found to be linked to the promoter region of both L-type HvleckRLK64 and HvleckRLK86. The promoter cis-acting regulatory element (CARE) analysis of barley identified a total of 75 CARE motifs responsive to light responsiveness (LR), tissue-specific (TS), hormone responsiveness (HR), and stress responsiveness (SR). The maximum number of CAREs was identified in HvleckRLK11 (25 for LR), HvleckRLK69 (17 for TS), and HvleckRLK80 (12 for HR). Additionally, HvleckRLK14, HvleckRLK16, HvleckRLK33, HvleckRLK50, HvleckRLK52, HvleckRLK56, and HvleckRLK110 were predicted to exhibit higher responses in stress conditions. In addition, 46 putative miRNAs were predicted to target 81 HvlecRLK genes and HvlecRLK13 was the most targeted gene by 8 different miRNAs. Protein-protein interaction analysis demonstrated higher functional similarities of 63 HvlecRLKs with 7 Arabidopsis STRING proteins. Our overall findings provide valuable information on the LecRLK gene family which might pave the way to advanced research on the functional mechanism of the candidate genes as well as to develop new barley cultivars in breeding programs.

Ameliorating the effects of multiple stresses on agronomic traits in crops: modern biotechnological and omics approaches

While global climate change poses a significant environmental threat to agriculture, the increasing population is another big challenge to food security. To address this, developing crop varieties with increased productivity and tolerance to biotic and abiotic stresses is crucial. Breeders must identify traits to ensure higher and consistent yields under inconsistent environmental challenges, possess resilience against emerging biotic and abiotic stresses and satisfy customer demands for safer and more nutritious meals. With the advent of omics-based technologies, molecular tools are now integrated with breeding to understand the molecular genetics of genotype-based traits and develop better climate-smart crops. The rapid development of omics technologies offers an opportunity to generate novel datasets for crop species. Identifying genes and pathways responsible for significant agronomic traits has been made possible by integrating omics data with genetic and phenotypic information. This paper discusses the importance and use of omics-based strategies, including genomics, transcriptomics, proteomics and phenomics, for agricultural and horticultural crop improvement, which aligns with developing better adaptability in these crop species to the changing climate conditions.

The Genetic Dissection of Nitrogen Use-Related Traits in Flax (Linum usitatissimum L.) at the Seedling Stage through the Integration of Multi-Locus GWAS, RNA-seq and Genomic Selection

Nitrogen (N), the most important macro-nutrient for plant growth and development, is a key factor that determines crop yield. Yet its excessive applications pollute the environment and are expensive. Hence, studying nitrogen use efficiency (NUE) in crops is fundamental for sustainable agriculture. Here, an association panel consisting of 123 flax accessions was evaluated for 21 NUE-related traits at the seedling stage under optimum N (N+) and N deficiency (N-) treatments to dissect the genetic architecture of NUE-related traits using a multi-omics approach integrating genome-wide association studies (GWAS), transcriptome analysis and genomic selection (GS). Root traits exhibited significant and positive correlations with NUE under N- conditions (r = 0.33 to 0.43, p < 0.05). A total of 359 QTLs were identified, accounting for 0.11% to 23.1% of the phenotypic variation in NUE-related traits. Transcriptomic analysis identified 1034 differentially expressed genes (DEGs) under contrasting N conditions. DEGs involved in N metabolism, root development, amino acid transport and catabolism and others, were found near the QTLs. GS models to predict NUE stress tolerance index (NUE_STI) trait were tested using a random genome-wide SNP dataset and a GWAS-derived QTLs dataset. The latter produced superior prediction accuracy (r = 0.62 to 0.79) compared to the genome-wide SNP marker dataset (r = 0.11) for NUE_STI. Our results provide insights into the QTL architecture of NUE-related traits, identify candidate genes for further studies, and propose genomic breeding tools to achieve superior NUE in flax under low N input.

Comparative analysis of physiological traits and gene expression patterns in nitrogen deficiency among barley cultivars

BACKGROUND

Nitrogen is one of the most important mineral nutrients for plants and is absorbed by the root system mainly in the inorganic form (NH+4 and NO-3). Plants absorb nitrogen as a food source for growth, biomass production, and development. Nitrogen is mainly absorbed as nitrate, which is the most common source of nitrogen available to higher plants. One of the unique features of nitrate transport is that NO-3 is both a substrate for transport and an inducer of NO-3 transport systems in genes and at physiological levels.

METHODS

In the present study, morphological and physiological traits (chlorophyll a/b, total chlorophyll, and carotenoid, antioxidant enzymes, and protein content), correlation between traits and gene expression, and principle component analysis of traits among five barley cultivars were measured in response to nitrogen deficiency (ND). The starved plants were transferred to a nutrient solution containing 0.2 mM and 2 mM NO-3 up to 7 and 14 days after ND application and non-stressed conditions, respectively.

RESULTS

Gene expression analysis revealed that the 10 HvNRT2 genes were induced in the leaf and root tissues at 7 and 14 days after ND treatments in five barley cultivars. Expression of NRT2 genes by relative quantitative qRT-PCR analysis for 10 HvNRT2 genes were determined. Based on the gene expression, HvNRT2.1, HvNRT2.2, and HvNRT2.4 were strongly induced by NO-3, peaking at 7 and 14 days after ND treatment. In contrast, the HvNRT2.4 showed only moderate induction in both leaves and roots. From our results, the Reyhan cultivar showed a significant increase in root fresh weight (RFW), protein content, and antioxidant enzyme activity in roots at 7 and 14 days after ND treatment as compared to the non-stressed condition. A highly positive correlation was observed between root catalase (CATr) and HvNRT2.2/2.5/2.6 leaves.

CONCLUSION

The expression of HvNRT2.4 is increased during long-term nitrogen starvation, while the expression of HvNRT2.1 and HvNRT2.2 are transiently increased by ND. Based on physiological and morphological traits and molecular mechanisms, the Reyhan is considered a tolerant cultivar under ND condition.

Amand P, Diggle AJ, Rengel Z. Genetic characterization of root architectural traits in barley (Hordeum vulgare L.) using SNP markers

Increasing attention is paid to providing new tools to breeders for targeted breeding for specific root traits that are beneficial in low-fertility, drying soils; however, such information is not available for barley (Hordeum vulgare L.). A panel of 191 barley accessions (originating from Australia, Europe, and Africa) was phenotyped for 26 root and shoot traits using the semi-hydroponic system and genotyped using 21 062 high-quality single nucleotide polymorphism (SNP) markers generated by genotyping-by-sequencing (GBS). The population structure analysis of the barley panel identified six distinct groups. We detected 1199 significant (P<0.001) marker-trait associations (MTAs) with r2 values up to 0.41. The strongest MTAs were found for root diameter in the top 20 cm and the longest root length. Based on the physical locations of these MTAs in the barley reference genome, we identified 37 putative QTLs for the root traits, and three QTLs for shoot traits, with nine QTLs located in the same physical regions. The genomic region 640-653 Mb on chromosome 7H was significant for five root length-related traits, where 440 annotated genes were located. The putative QTLs for various root traits identified in this study may be useful for genetic improvement regarding the adaptation of new barley cultivars to suboptimal environments and abiotic stresses.

Multi-scale analysis provides insights into the roles of ureide permeases in wheat nitrogen use efficiency

The ureides allantoin and allantoate serve as nitrogen (N) transport compounds in plants, and more recently, allantoin has been shown to play a role in signaling. In planta, tissue ureide levels are controlled by the activity of enzymes of the purine degradation pathway and by ureide transporters called ureide permeases (UPS). Little is known about the physiological function of UPS proteins in crop plants, and especially in monocotyledon species. Here, we identified 13 TaUPS genes in the wheat (Triticum aestivum L.) genome. Phylogenetic and genome location analyses revealed a close relationship of wheat UPSs to orthologues in other grasses and a division into TaUPS1, TaUPS2.1, and TaUPS2.2 groups, each consisting of three homeologs, with a total of four tandem duplications. Expression, localization, and biochemical analyses resolved spatio-temporal expression patterns of TaUPS genes, transporter localization at the plasma membrane, and a role for TaUPS2.1 proteins in cellular import of ureides and phloem and seed loading. In addition, positive correlations between TaUPS1 and TaUPS2.1 transcripts and ureide levels were found. Together the data support that TaUPSs function in regulating ureide pools at source and sink, along with source-to-sink transport. Moreover, comparative studies between wheat cultivars grown at low and high N strengthened a role for TaUPS1 and TaUPS2.1 transporters in efficient N use and in controlling primary metabolism. Co-expression, protein-protein interaction, and haplotype analyses further support TaUPS involvement in N partitioning, N use efficiency, and domestication. Overall, this work provides a new understanding on UPS transporters in grasses as well as insights for breeding resilient wheat varieties with improved N use efficiency.

Use of GWAS analysis in deciphering the inability of barley seeds to germinate after hypoxia

This article comments on:

Gómez-Álvarez EM, Tondelli A, Nghi KN, Voloboeva V, Giordano G, Valè G, Perata P, Pucciariello C. 2023. The inability of barley to germinate after submergence depends on hypoxia-induced secondary dormancy. Journal of Experimental Botany 74, 4277–4289

Genetic resources and precise gene editing for targeted improvement of barley abiotic stress tolerance

Abiotic stresses, predominately drought, heat, salinity, cold, and waterlogging, adversely affect cereal crops. They limit barley production worldwide and cause huge economic losses. In barley, functional genes under various stresses have been identified over the years and genetic improvement to stress tolerance has taken a new turn with the introduction of modern gene-editing platforms. In particular, clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) is a robust and versatile tool for precise mutation creation and trait improvement. In this review, we highlight the stress-affected regions and the corresponding economic losses among the main barley producers. We collate about 150 key genes associated with stress tolerance and combine them into a single physical map for potential breeding practices. We also overview the applications of precise base editing, prime editing, and multiplexing technologies for targeted trait modification, and discuss current challenges including high-throughput mutant genotyping and genotype dependency in genetic transformation to promote commercial breeding. The listed genes counteract key stresses such as drought, salinity, and nutrient deficiency, and the potential application of the respective gene-editing technologies will provide insight into barley improvement for climate resilience.

Haplotype-Based Genome-Wide Association Analysis Using Exome Capture Assay and Digital Phenotyping Identifies Genetic Loci Underlying Salt Tolerance Mechanisms in Wheat

Soil salinity can impose substantial stress on plant growth and cause significant yield losses. Crop varieties tolerant to salinity stress are needed to sustain yields in saline soils. This requires effective genotyping and phenotyping of germplasm pools to identify novel genes and QTL conferring salt tolerance that can be utilised in crop breeding schemes. We investigated a globally diverse collection of 580 wheat accessions for their growth response to salinity using automated digital phenotyping performed under controlled environmental conditions. The results show that digitally collected plant traits, including digital shoot growth rate and digital senescence rate, can be used as proxy traits for selecting salinity-tolerant accessions. A haplotype-based genome-wide association study was conducted using 58,502 linkage disequilibrium-based haplotype blocks derived from 883,300 genome-wide SNPs and identified 95 QTL for salinity tolerance component traits, of which 54 were novel and 41 overlapped with previously reported QTL. Gene ontology analysis identified a suite of candidate genes for salinity tolerance, some of which are already known to play a role in stress tolerance in other plant species. This study identified wheat accessions that utilise different tolerance mechanisms and which can be used in future studies to investigate the genetic and genic basis of salinity tolerance. Our results suggest salinity tolerance has not arisen from or been bred into accessions from specific regions or groups. Rather, they suggest salinity tolerance is widespread, with small-effect genetic variants contributing to different levels of tolerance in diverse, locally adapted germplasm.

Mini review: Targeting below-ground plant performance to improve nitrogen use efficiency (NUE) in barley

Nitrogen (N) fertilizer is one of the major inputs for grain crops including barley and its usage is increasing globally. However, N use efficiency (NUE) is low in cereal crops, leading to higher production costs, unfulfilled grain yield potential and environmental hazards. N uptake is initiated from plant root tips but a very limited number of studies have been conducted on roots relevant to NUE specifically. In this review, we used barley, the fourth most important cereal crop, as the primary study plant to investigate this topic. We first highlighted the recent progress and study gaps in genetic analysis results, primarily, the genome-wide association study (GWAS) regarding both biological and statistical considerations. In addition, different factors contributing to NUE are discussed in terms of root morphological and anatomical traits, as well as physiological mechanisms such as N transporter activities and hormonal regulation.

Transcriptome and Metabolome Reveal the Molecular Mechanism of Barley Genotypes Underlying the Response to Low Nitrogen and Resupply

Nitrogen is one of the most important mineral elements for plant growth and development. Excessive nitrogen application not only pollutes the environment, but also reduces the quality of crops. However, are few studies on the mechanism of barley tolerance to low nitrogen at both the transcriptome and metabolomics levels. In this study, the nitrogen-efficient genotype (W26) and the nitrogen-sensitive genotype (W20) of barley were treated with low nitrogen (LN) for 3 days and 18 days, then treated with resupplied nitrogen (RN) from 18 to 21 days. Later, the biomass and the nitrogen content were measured, and RNA-seq and metabolites were analyzed. The nitrogen use efficiency (NUE) of W26 and W20 treated with LN for 21 days was estimated by nitrogen content and dry weight, and the values were 87.54% and 61.74%, respectively. It turned out to have a significant difference in the two genotypes under the LN condition. According to the transcriptome analysis, 7926 differentially expressed genes (DEGs) and 7537 DEGs were identified in the leaves of W26 and W20, respectively, and 6579 DEGs and 7128 DEGs were found in the roots of W26 and W20, respectively. After analysis of the metabolites, 458 differentially expressed metabolites (DAMs) and 425 DAMs were found in the leaves of W26 and W20, respectively, and 486 DAMs and 368 DAMs were found in the roots of W26 and W20, respectively. According to the KEGG joint analysis of DEGs and DAMs, it was discovered that glutathione (GSH) metabolism was the pathway of significant enrichment in the leaves of both W26 and W20. In this study, the metabolic pathways of nitrogen metabolism and GSH metabolism of barley under nitrogen were constructed based on the related DAMs and DEGs. In leaves, GSH, amino acids, and amides were the main identified DAMs, while in roots, GSH, amino acids, and phenylpropanes were mainly found DAMs. Finally, some nitrogen-efficient candidate genes and metabolites were selected based on the results of this study. The responses of W26 and W20 to low nitrogen stress were significantly different at the transcriptional and metabolic levels. The candidate genes that have been screened will be verified in future. These data not only provide new insights into how barley responds to LN, but also provide new directions for studying the molecular mechanisms of barley under abiotic stress.

The Genetic Architecture of Nitrogen Use Efficiency in Switchgrass (Panicum virgatum L.)

Switchgrass (Panicum virgatum L.) has immense potential as a bioenergy crop with the aim of producing biofuel as an end goal. Nitrogen (N)-related sustainability traits, such as nitrogen use efficiency (NUE) and nitrogen remobilization efficiency (NRE), are important factors affecting switchgrass quality and productivity. Hence, it is imperative to develop nitrogen use-efficient switchgrass accessions by exploring the genetic basis of NUE in switchgrass. For that, we used 331 diverse field-grown switchgrass accessions planted under low and moderate N fertility treatments. We performed a genome wide association study (GWAS) in a holistic manner where we not only considered NUE as a single trait but also used its related phenotypic traits, such as total dry biomass at low N and moderate N, and nitrogen use index, such as NRE. We have evaluated the phenotypic characterization of the NUE and the related traits, highlighted their relationship using correlation analysis, and identified the top ten nitrogen use-efficient switchgrass accessions. Our GWAS analysis identified 19 unique single nucleotide polymorphisms (SNPs) and 32 candidate genes. Two promising GWAS candidate genes, caffeoyl-CoA O-methyltransferase (CCoAOMT) and alfin-like 6 (AL6), were further supported by linkage disequilibrium (LD) analysis. Finally, we discussed the potential role of nitrogen in modulating the expression of these two genes. Our findings have opened avenues for the development of improved nitrogen use-efficient switchgrass lines.

CRISPR/Cas9 gene editing and natural variation analysis demonstrate the potential for HvARE1 in improvement of nitrogen use efficiency in barley

Nitrogen is a major determinant of grain yield and quality. As excessive use of nitrogen fertilizer leads to environmental pollution and high production costs, improving nitrogen use efficiency (NUE) is fundamental for a sustainable agriculture. Here, we dissected the role of the barley abnormal cytokinin response1 repressor 1 (HvARE1) gene, a candidate for involvement in NUE previously identified in a genome-wide association study, through natural variation analysis and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9)-mediated gene editing. HvARE1 was predominantly expressed in leaves and shoots, with very low expression in roots under low nitrogen conditions. Agrobacterium-mediated genetic transformation of immature embryos (cv. Golden Promise) with single guide RNAs targeting HvARE1 generated 22 T₀ plants, from which four T₁ lines harbored missense and/or frameshift mutations based on genotyping. Mutant are1 lines exhibited an increase in plant height, tiller number, grain protein content, and yield. Moreover, we observed a 1.5- to 2.8-fold increase in total chlorophyll content in the flag leaf at the grain filling stage. Delayed senescence by 10-14 d was also observed in mutant lines. Barley are1 mutants had high nitrogen content in shoots under low nitrogen conditions. These findings demonstrate the potential of ARE1 in NUE improvement in barley.

Genotypic Variation of Nitrogen Use Efficiency and Amino Acid Metabolism in Barley

Owing to the large genetic diversity of barley and its resilience under harsh environments, this crop is of great value for agroecological transition and the need for reduction of nitrogen (N) fertilizers inputs. In the present work, we investigated the diversity of a North African barley genotype collection in terms of growth under limiting N (LN) or ample N (HN) supply and in terms of physiological traits including amino acid content in young seedlings. We identified a Moroccan variety, Laanaceur, accumulating five times more lysine in its leaves than the others under both N nutritional regimes. Physiological characterization of the barley collection showed the genetic diversity of barley adaptation strategies to LN and highlighted a genotype x environment interaction. In all genotypes, N limitation resulted in global biomass reduction, an increase in C concentration, and a higher resource allocation to the roots, indicating that this organ undergoes important adaptive metabolic activity. The most important diversity concerned leaf nitrogen use efficiency (LNUE), root nitrogen use efficiency (RNUE), root nitrogen uptake efficiency (RNUpE), and leaf nitrogen uptake efficiency (LNUpE). Using LNUE as a target trait reflecting barley capacity to deal with N limitation, this trait was positively correlated with plant nitrogen uptake efficiency (PNUpE) and RNUpE. Based on the LNUE trait, we determined three classes showing high, moderate, or low tolerance to N limitation. The transcriptomic approach showed that signaling, ionic transport, immunity, and stress response were the major functions affected by N supply. A candidate gene encoding the HvNRT2.10 transporter was commonly up-regulated under LN in the three barley genotypes investigated. Genes encoding key enzymes required for lysine biosynthesis in plants, dihydrodipicolinate synthase (DHPS) and the catabolic enzyme, the bifunctional Lys-ketoglutarate reductase/saccharopine dehydrogenase are up-regulated in Laanaceur and likely account for a hyperaccumulation of lysine in this genotype. Our work provides key physiological markers of North African barley response to low N availability in the early developmental stages.

Phylogenomic and Microsynteny Analysis Provides Evidence of Genome Arrangements of High-Affinity Nitrate Transporter Gene Families of Plants

Nitrate transporter 2 (NRT2) and NRT3 or nitrate-assimilation-related 2 (NAR2) proteins families form a two-component, high-affinity nitrate transport system, which is essential for the acquisition of nitrate from soils with low N availability. An extensive phylogenomic analysis across land plants for these families has not been performed. In this study, we performed a microsynteny and orthology analysis on the NRT2 and NRT3 genes families across 132 plants (Sensu lato) to decipher their evolutionary history. We identified significant differences in the number of sequences per taxonomic group and different genomic contexts within the NRT2 family that might have contributed to N acquisition by the plants. We hypothesized that the greater losses of NRT2 sequences correlate with specialized ecological adaptations, such as aquatic, epiphytic, and carnivory lifestyles. We also detected expansion on the NRT2 family in specific lineages that could be a source of key innovations for colonizing contrasting niches in N availability. Microsyntenic analysis on NRT3 family showed a deep conservation on land plants, suggesting a high evolutionary constraint to preserve their function. Our study provides novel information that could be used as guide for functional characterization of these gene families across plant lineages.

Omics for the Improvement of Abiotic, Biotic, and Agronomic Traits in Major Cereal Crops: Applications, Challenges, and Prospects

Omics technologies, namely genomics, transcriptomics, proteomics, metabolomics, and phenomics, are becoming an integral part of virtually every commercial cereal crop breeding program, as they provide substantial dividends per unit time in both pre-breeding and breeding phases. Continuous advances in omics assure time efficiency and cost benefits to improve cereal crops. This review provides a comprehensive overview of the established omics methods in five major cereals, namely rice, sorghum, maize, barley, and bread wheat. We cover the evolution of technologies in each omics section independently and concentrate on their use to improve economically important agronomic as well as biotic and abiotic stress-related traits. Advancements in the (1) identification, mapping, and sequencing of molecular/structural variants; (2) high-density transcriptomics data to study gene expression patterns; (3) global and targeted proteome profiling to study protein structure and interaction; (4) metabolomic profiling to quantify organ-level, small-density metabolites, and their composition; and (5) high-resolution, high-throughput, image-based phenomics approaches are surveyed in this review.

Using chlorate as an analogue to nitrate to identify candidate genes for nitrogen use efficiency in barley

Nitrogen (N) is one of the most important macronutrients for crop growth and development. Large amounts of N fertilizers are applied exogenously to improve grain yield and quality, which has led to environmental pollution and high cost of production. Therefore, improvement of N use efficiency (NUE) is a very important aspect for sustainable agriculture. Here, a pilot experiment was firstly conducted with a set of barley genotypes with confirmed NUE to validate the fast NUE screening, using chlorate as an analogue to nitrate. High NUE genotypes were susceptible to chlorate-induced toxicity whereas the low NUE genotypes were tolerant. A total of 180 barley RILs derived from four parents (Compass, GrangeR, Lockyer and La Trobe) were further screened for NUE. Leaf chlorosis induced by chlorate toxicity was the key parameter observed which was later related to low-N tolerance of the RILs. There was a distinct variation in chlorate susceptibility of the RILs with leaf chlorosis in the oldest leaf ranging from 10 to 80%. A genome-wide association study (GWAS) identified 9 significant marker-trait associations (MTAs) conferring high chlorate sensitivity on chromosomes 2H (2), 3H (1), 4H (4), 5H (1) and Un (1). Genes flanking with these markers were retrieved as potential targets for genetic improvement of NUE. Genes encoding Ferredoxin 3, leucine-rich receptor-like protein kinase family protein and receptor kinase are responsive to N stress. MTA4H5468 which exhibits concordance with high NUE phenotype can further be explored under different genetic backgrounds and successfully applied in marker-assisted selection.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01239-8.